Discussion :

[sigma10]

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 :

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
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;

[Invité]

Bonjour,

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

Code :

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

1
2
3
z=abs(X);
figure(3) ; plot(abs(w));
k=fft(abs(w));

Qui devrait alors être:

Code :

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

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 :

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

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

[sigma10]

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

cette partie n'est claire ?

Code :

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

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 :

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

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

[Invité]

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

[sigma10]

voila ce que je fais mais ca marche pas ?

Code :

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

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 :

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

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

Il me donne pas le signal original

[Invité]

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

Code :

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

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

[sigma10]

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

[compaq610]

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

[Invité]

Tu as encore laissé un abs quelque part!
Copie simplement le code de mon premier message, c'est si compliqué que ç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
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;

[compaq610]

bonjour;
utilisé le code de "Winjerome" ça marche bien

[sigma10]

voila les resultats qui me donne :


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

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

[sigma10]

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 :

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