Discussion :

[jigsa]

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

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 :

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

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

[Invité]

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?

[jigsa]

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 :

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

 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 :

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

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

D'ou mon problème...

Code :

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

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')

[Invité]

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.

[jigsa]

En effet je dois donc taper ceci

Code :

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

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

[Invité]

Des petites erreurs:

Code :

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

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

[jigsa]

Merci bcp,
J'avais fait une fft de trop

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