Déconvolution via FFT, IFFT

Bonjour tout le monde est plus spécialement au habitué du forum :ccool:

Voila j'ai mon programme suivant qui a pour but de bruité un signal type sinus via une convolution, et d'utilisé une déconvolution pour retrouver mon signal de basse or je le retrouve bien mais avec un déphasage de pi/2

Code:

---

1 2 3 4 5

x= -pi:0.01:pi % Generate Sinusoidal signal s=sin(x); n=randn(1,629); % Generate Guassion Noise z= fft(b).*fft(n) y=ifft(z)

---

Y a t'il un moyen de retrouver notre signal non décalé ?

Merci de votre aide

Bonjour,

Tu parles de convolution/déconvolution, mais dans ton code, je ne vois qu'une convolution...
De plus que représente b dans z= fft(b).*fft(n)? Est-ce s?

Citation:

---

Envoyé par Winjerome

Bonjour,

Tu parles de convolution/déconvolution, mais dans ton code, je ne vois qu'une convolution...
De plus que représente b dans z= fft(b).*fft(n)? Est-ce s?

---

En effet c'est bien

Code:

---

z= fft(s).*fft(n)

---

La convolution de 2 signaux est la multiplication des 2 transformées de fourier
Donc la déconvolution est la Transformée inverse de fourier du produit des 2 signaux...

Code:

---

1 2 3

z= fft(s).fft(n) % convolution y=ifft(z) % déconvolution

---

D'ou mon problème...

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

x= -pi:0.01:pi % Generate Sinusoidal signal s=sin(x); n=randn(1,length(s)); z= fft(s).*fft(n) % convolution y=ifft(z) % déconvolution subplot(121) plot(s) title ('Signal réel S') subplot(122) plot(y) title ('signal déconvolué Y')

---

Citation:

---

Envoyé par jigsa

La convolution de 2 signaux est la multiplication des 2 transformées de fourier
Donc la déconvolution est la Transformée inverse de fourier du produit des 2 signaux...

---

Non ceci n'est que l'opération de convolution, avec l'ifft, tu ne fais que repasser dans le domaine temporel.

En effet je dois donc taper ceci

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

x= -pi:0.01:pi % Generate Sinusoidal signal s=sin(x); n=randn(1,length(s)); z= fft(s).*fft(n) % convolution Z=fft(z) y=fft(Z)./fft(n) % déconvolution y=ifft(z) % déconvolution subplot(121) plot(s) title ('Signal réel S') subplot(122) plot(y) title ('signal déconvolué Y')

---

Mais j'ai toujours un problème récurant, qui est un gain (environ 18) et un déphasage changeant...
Je ne vois vraiment pas comment pourquoi

Des petites erreurs:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

x= -pi:0.01:pi; % Generate Sinusoidal signal s=sin(x); n=randn(1,length(s)); z= fft(s).*fft(n); % convolution Z=ifft(z); y=fft(Z)./fft(n); % déconvolution y=ifft(y); % déconvolution subplot(121) plot(s) title ('Signal réel S') subplot(122) plot(y) title ('signal déconvolué Y')

---

Remarque: pense à mettre des ; à la fin de tes lignes, toutes les valeurs ne sont alors plus affichées dans ta ligne de commande ;)

Merci bcp,:ccool:
J'avais fait une fft de trop

ps: Il se peut que je vous sollicite de nouveau dans un avenir proche