Discussion :

[sofiyya]

Bonjour,

ça fait deux mois que je galère avec la déconvolution de deux signaux et j'arrive pas à avoir un bon résultat. Est ce quelqu'un peut m'aider?
Voici le code matlab que j'utilise ainsi que le résulat:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
vect1=[1 2 3 4 5 6 7 8 9];
vect2=[1 2 3 4 0 0 0 0 0 0 0 0 0 0]; %Je cherche w telque conv(w,vect1)=vect2
Lx=length(vect2)-length(vect1)+1;
Lx2=pow2(nextpow2(Lx));
VECT1=fft(vect1,Lx2);
VECT2=fft(vect2,Lx2);
W=VECT2 ./ VECT1;
w=real(ifft(W,Lx2));
w=w(1:1:Lx);
w=w/max(abs(w));
conv(w,vec1)

donne

Code :

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

1
2
3
4
5
6
7
8
9
10
 
ans =
 
Columns 1 through 9
 
-0.1688 -0.5065 -1.0130 -1.6883 -1.3636 -0.9740 -0.5844 -0.1948 0.1948
 
Columns 10 through 14
 
2.2727 4.5195 6.9351 9.5195 0.5844

qui est totalement différent du résultat attendu qui est vect2


Merci beaucoup pour votre aide...

[phryte]

Bonjour.
Plus proche du résultat :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
clear
vect1=[1 2 3 4 5 6 7 8 9];
vect2=[1 2 3 4 0 0 0 0 0 0 0 0 0 0]; %Je cherche w telque conv(w,vect1)=vect2
Lx=length(vect2)-length(vect1)+1;
Lx2=pow2(nextpow2(Lx));
VECT1=abs(fft(vect1,Lx2));
VECT2=abs(fft(vect2,Lx2));
W=VECT2./VECT1;
w=(ifft(W,Lx2));
w=w(1:1:Lx);
w=w/max(abs(w));
conv(w,vect1)

[sofiyya]

Oui mais ça ne résout pas le problème, est ce que vous pensez que ça vient de la déconvolution? est ce que la déconvolution ne donne pas un bon résultat pour tout les coups?
ça fait 2 mois que je travaille dessous et je trouve pas de solution, et maintenant ça devient urgent!

Merci beaucoup pour votre aide

[phryte]

Bonjour.

ça vient de la déconvolution

Non, cette fonction est bonne :

Code :

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

1
2
3
4
vect1=[1 2 3 4 5 6 7 8 9];
vect2=[1 2 3 4 0 0 0 0 0 0 0 0 0 0];
c = conv(vect1,vect2)
[q,r] = deconv(c,vect1)

[sofiyya]

Merci pour ta réponse..

Oui ça explique bien que ça dépend des vecteurs à déconvoluer. Parce que pour mon exemple, même la fonction deconv ne donne pas un bon résultat! Je présume donc que pour mon cas, y'a pas une solution et c'est à cause des vecteurs utilisé et non pas de la manière de faire la déconvolution.