Discussion :

[kheira]

bonjour

Je veux appliquer une ondelette de Morlet sur un signal biomédical pour faire une extraction de caractéristiques (puissance de bande) et faire un seuillage après et je ne sais pas par quoi commencer.

Je voudrais de l'aide SVP

[Jerome Briot]

Quel est le problème et surtout quel est le rapport avec MATLAB ?

As-tu au moins compris la méthode ?

[kmaniche]

Pour la théorie tu as [ame="http://en.wikipedia.org/wiki/Morlet_wavelet"]WIKI >>[/ame]

• Continuous wavelet
• Mexican hat wavelet
• Hermitian wavelet
• Complex mexican hat wavelet
• Difference of Gaussians
• Hermitian hat wavelet
• Continuous wavelet transform
• Beta wavelet
• Modified Morlet wavelet
• Shannon wavelet

[kheira]

Je veux programmer ce qui est identique à la fonction CWT (continues wavelet transform 1D) de matlab. Je ne veux pas faire appel à cette fonction mais la faire moi même. Je n'y arrive pas !

ça ne donne rien avec ça :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
t = 0:1/250:7.5-1/250;
signal=load('C3.txt');
y=signal(76876:78750);
clf
ii=sqrt(-1);
%
%..signal a analyser........
%
n=length(y)-1;
p=fix(log(n)/log(2));
dx=1/n;
%t=0:dx:1;
x=-t;
%
% calcul de l'ondelette
%
a0=1/2;
a0=a0^(1/2);
p=2*p;
W=zeros(p,n+1);
for i=p-1:-1:0
	disp('echelle')
	a=a0^i;
	xx=x/a; 
	g=(1-2*xx.^2).*exp(-2*xx.^2)./sqrt(a);
	plot(x,g),title('ondelette chapeau mexicain')
	g=exp(-xx.^2/2).*exp(ii*5*xx)./sqrt(a);
   plot(x,real(g)),title('ondelette de Morlet-partie reelle ')
   disp('appuyer sur une touche pour continuer')
   pause
 
 calcul des coefficients d'ondelettes a l'echelle a
 
	wa=conv(y,g);
	W(i+1,1:n)=abs(wa(n+1:2*n));
end
%
% trace
%
colormap(jet);
subplot(3,1,1),plot(t(1:n),y(1:n)),title('fonction a analyser')
subplot(3,1,2),colormap(jet),title('coefficients d ondelettes')
subplot(3,1,3),contour(W),title('lignes isocontour')
 
%pcolor(W(1:p,1:n))

je veux bien une aide SVP

[Jerome Briot]

je veux programmer ce qui édentique a la fonction CWT (continues wavelet transform 1D) de matlab

Tu trouveras des informations utiles ici : How is the CWT function implemented in Wavelet Toolbox?

[kheira]

Pouvez-vous me dire où je peux trouver la fonction cwt.m utilisée par matlab pour que je m'inspire d'elle car je n'ai pas bien compris la documentation ?

svp

[kmaniche]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
>> wich cwt
C:\Program Files\MATLAB\R2009b\toolbox\wavelet\wavelet\cwt.m
 
>> edit(which('cwt'))

[kheira]

j'arrive toujour pas a calculer les coefissians de l'ondelette de morlet moi meme son passé par la fonction CWT

avec mon premiers code il donnepas le meme resultat comme la fonction cwt !!

ya pas quelqu'un qui peut maider

[kmaniche]

Est-ce que tu as réussi à comprendre le contenu de "CWT.m" ?

[kheira]

J’ai essayé de m’inspirer de la fonction de CWT pour faire ma propre CWT mai j’ai pas compris quelque routine, j’ai ou le même résultat

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
>> y=[ 0.123 1.45 0.67]
coefs = cwt(y,1:3,'morl')
 
coefs =
 
    0.1611   -0.1142   -0.2464
   -0.3223    0.0428    0.3312
   -0.6138    0.2950    0.6813

Lors de l’exécution de ma fonction MACWT j’ai ou les mêmes résultats

Mai je veux faire quelque modification par exemple dans la ces routine :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
% ----------- variable constante de morlet ---------  
precis = 10; 
[val_WAV,xWAV] = intwave('morl',precis);
stepWAV = xWAV(2)-xWAV(1);
wtype = wavemngr('type','morl');  % Wtype = 4 
if wtype==5 , val_WAV = conj(val_WAV); end
getINTEG = 0;
getWTYPE = 0;
xWAV = xWAV-xWAV(1);
xMaxWAV = xWAV(end);

je ponce que ce sont des variable fixé de mortel que on utilise Sour tout signal ?
es que je peut sauvegarder ces variable dans un fichiers text et travailler avec ?
je veux pas passer par :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
intwave
wavemngr ?!

Et je veux pas passer aussi les fonctions :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
floor
fliplr
wkeep1
..

Je veux les faire moi m’aime , alors je veux bien un aide svp … je fait ça car je veux reprogrammé la même chose après en C++ alors je veux bien un aide SVP , ma fonction et dans le piece jointe