compression d'un signal par FFT

Version imprimable

bonsoir
voila j'ai fait un programme qui fait la compression d'un signal par FFT mais le problème est que je ne retrouve pas le signal original lorsque je fais la reconstruction.

Code:

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
clear all
fe=300;
f1=10;
f2=60;
f3=90;
t=0:1/fe:1;
x=2*sin(2*pi*t*f1)+10*sin(2*pi*t*f2)+ 20*sin(2*pi*t*f3);
figure(1); plot(x)
X=fft(x);
z=abs(X);
figure(2) ;plot(z)
seuil=0.5;
Xseuil=seuillage(z,seuil);
w=ifft(Xseuil);
figure(3) ; plot(abs(w));
k=fft(abs(w));
figure(4); plot(abs(k))
 
------------------------------
%% function Xseuil=seuillage(z,seuil)
 
%en entree: 
%z:les coefficients a seuiller
%seuil: le seuil qu'on veut
%en sortie
%Xseuil: les coefficients seuillés
 
Xseuil=z;
u=find(abs(z)<seuil);
Xseuil(u)=0;

13/06/2011, 11h27
Invité
Bonjour,

Ton problème vient du calcul de la valeur absolue à certains endroits où tu ne devrais pas, en particulier là:
Code:

1 2 3 z=abs(X); figure(3) ; plot(abs(w)); k=fft(abs(w));
Qui devrait alors être:
Code:

1 2 3 4 5 6 7 8 X=fft(x); figure(2) ;plot(abs(X)) seuil=0.5; Xseuil=seuillage(X,seuil); w=ifft(Xseuil); figure(3) ; plot(real(w)); k=fft(w); figure(4); plot(abs(k))
Au passage, tu peux simplifier ta fonction de seuillage avec l'indexage logique et en sortant la même variable z d'entrée:
Code:

1 2 3 function z=seuillage(z,seuil) z(abs(z)<seuil)=0;

mais il me donne pas le signal origine lorsque je fais la reconstruction

cette partie n'est claire ?
Code:

1 2 3 4 5 6 7 X=fft(x); z=abs(X); figure(2) ;plot(z) seuil=0.5; Xseuil=seuillage(z,seuil); w=ifft(Xseuil); figure(3) ; plot(real(w));
%---------------------------
Code:

1 2 3 function Xseuil=seuillage(z,seuil) Xseuil=z; z(abs(z)<seuil)=0;

13/06/2011, 12h28
Invité

Il suffisait de recopier mon code, chose que tu n'as apparemment pas faite correctement:?

voila ce que je fais mais ca marche pas ?
Code:

1 2 3 4 5 6 7 8 9 10 X=fft(x); z=abs(X); figure(2) ;plot(z) seuil=5; z=seuillage(z,seuil); w=ifft(z); figure(3) ; plot(real(w)); k=fft(w); figure(4); plot(abs(k))
%-----------------------------
Code:

1 2 3 function z=seuillage(z,seuil) z(abs(z)<seuil)=0;
Il me donne pas le signal original

C'est pourtant pas bien compliqué de faire du copier/coller :cry:

Code:

1
2
3
4
5
6
7
8
9
X=fft(x);
z=abs(X); % z = X simplement
figure(2) ;plot(abs(z)) % seulement pour visualiser on prend la valeur absolue
seuil=5;
z=seuillage(z,seuil);
w=ifft(z);
figure(3) ; plot(real(w));
k=fft(w);
figure(4); plot(abs(k))

13/06/2011, 16h14
sigma10

1 pièce(s) jointe(s)

oui désolé j'ai fait le titre
le figure 1 c'est le signal original et 2 est le signal reconstruit
mon probleme : est ce que mon programme est fait la compression ?

j'attend votre reponse
13/06/2011, 16h19
compaq610

bonjour;
le signale reconstruite ne ressemble pas au signale originale je pense que vous avez un problème dans ton code.

Tu as encore laissé un abs quelque part!
Copie simplement le code de mon premier message, c'est si compliqué que ça? :mur:
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 function compression fe=300; f1=10; f2=60; f3=90; t=0:1/fe:1; x=2*sin(2*pi*t*f1)+10*sin(2*pi*t*f2)+ 20*sin(2*pi*t*f3); figure(1); plot(x) X=fft(x); figure(2) ;plot(abs(X)) seuil=5; Xseuil=seuillage(X,seuil); w=ifft(Xseuil); figure(3) ; plot(real(w)); function z=seuillage(z,seuil) z(abs(z)<seuil)=0;

13/06/2011, 16h36
compaq610

bonjour;
utilisé le code de "Winjerome" ça marche bien
13/06/2011, 16h43
sigma10

1 pièce(s) jointe(s)

voila les resultats qui me donne :
Pièce jointe 78895

Compap610 :Est ce que vous avez trouvé la mm chose

pour avoir une bonne compression , quelle valeur il faut prendre ?

Compression d'un signal par FFT

Bonjour
je veux savoir est ce que mon programme est ça marche pour la compression ?
et pour quel valeur de seuillage je vais trouvé une bonne compression ?
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 clear all close all fe=300; f1=10; f2=60; f3=90; t=0:1/fe:1; x=2*sin(2*pi*t*f1)+10*sin(2*pi*t*f2)+ 20*sin(2*pi*t*f3); figure(1); plot(x) title('sign original') X=fft(x); figure(2) ;plot(abs(X)) seuil=50; Xseuil=seuillage(X,seuil); w=ifft(Xseuil); figure(3) ; plot(real(w)); title('sign reconstruit') error=(x-w); figure(4) ; plot(error)