filtrage d'une fréquence particulière sous matlab

[invite9c7554e3]

Salut tous,

je vous écris car j'essai de faire du filtrage sous maltab.

dans le code ci dessous j'ai un signal qui est composé de deux sinus de fréquence 2 et 0.3Hz.

ce que j'aimerai faire c'est filtrer le sinus 2Hz mais le code ci dessous ne m'enleve pas cette composante du signal...

savez vous d'où cela vient ?

clear all;close all;clc;
%% signal origine
t = 0:0.01:10;
y = sin(2*pi*0.3*t)+sin(2*pi*2*t);
%% signal traité fir
fs = 1000;
fcoupure = 2*[1 3]/fs;
[b]=fir1(1,fcoupure,'stop');
yfilter=filter(b,1,y);
figure (1);hold on;plot(t,y);plot(t,yfilter,'r');
%% signal traité avec butter
fs = 1000;
fcoupure = 2*[1 3]/fs;
[b,a]=butter(1,fcoupure,'stop');
yfilter=filter(b,a,y);
figure (2);hold on;plot(t,y);plot(t,yfilter,'r');

-----

[phuphus]

Bonsoir,

ta fréquence d'échantillonnage n'est pas en accord avec ton pas de temps. Change plutôt la fréquence et garde le pas de temps actuel, il est toujours plus délicat de filtrer des fréquences très petites devant Fs.

Pour ta prochaine question : mets plus que "1".

[invite9c7554e3]

salut Phuphus et merci beaucoup d'avoir pris le temps de repondre !!!

ta fréquence d'échantillonnage n'est pas en accord avec ton pas de temps.

pourquoi ?
je ne comprends pas le lien entre mon pas de temps et la fréquence d'echantillonage qu'il faut choisir ...
=> le pas d'echantillonage doit etre soit supérieur soit egale à mon pas de temps ?

Change plutôt la fréquence et garde le pas de temps actuel, il est toujours plus délicat de filtrer des fréquences très petites devant Fs.

ça roule, merci.

Pour ta prochaine question : mets plus que "1".

mdr, tu parle de l'ordre du signal je pense
oui, en effet je vais essayer de changer ça

[phuphus]

Bonsoir,

pourquoi ?
je ne comprends pas le lien entre mon pas de temps et la fréquence d'echantillonage qu'il faut choisir ...
=> le pas d'echantillonage doit etre soit supérieur soit egale à mon pas de temps ?

En prenant un pas de temps de 10ms, tu construis un signal échantillonné à 100Hz. Il faut donc garder 100Hz.

mdr, tu parle de l'ordre du signal je pense
oui, en effet je vais essayer de changer ça

Tout juste !

[A voir en vidéo sur Futura]

[invite9c7554e3]

merci de ton aide !

[phuphus]

Bonsoir,

en référence au fil "inverse FFT (matlab)"

clear all;close all;clc;
%% signal origine
t = 0:0.001:10;
y = sin(2*pi*0.3*t)+sin(2*pi*2*t);
%% signal traité fir
fs = 1000;
fcoupure = 2*[1 10]/fs;
[b]=fir1(2000,fcoupure,'stop');
yfilter=filter(b,1,y);
figure (1);hold on;plot(t,y);plot(t,yfilter,'r');
%% signal traité avec butter
fs = 1000;
fcoupure = 2*[1 3]/fs;
[b,a]=butter(3,fcoupure,'stop');
yfilter=filter(b,a,y);
figure (2);hold on;plot(t,y);plot(t,yfilter,'r');

[phuphus]

Re,

légère amélioration (déphasage nul) :

clear all;close all;clc;
%% signal origine
t = 0:0.001:10;
y = sin(2*pi*0.3*t)+sin(2*pi*2*t);
%% signal traité fir
fs = 1000;
fcoupure = 2*[1 10]/fs;
[b]=fir1(2000,fcoupure,'stop');
yfilter=filtfilt(b,1,y);
figure (1);hold on;plot(t,y);plot(t,yfilter,'r');
%% signal traité avec butter
fs = 1000;
fcoupure = 2*[1 3]/fs;
[b,a]=butter(3,fcoupure,'stop');
yfilter=filtfilt(b,a,y);
figure (2);hold on;plot(t,y);plot(t,yfilter,'r');

[invite9c7554e3]

salut Phuphus !

merci beaucoup d'avoir pris le temps de repondre !!!!
j'ai bien compris ton exemple il est très bien !!!!!!!

je crois que je vais arriver à avoir des trucs pas mal avec ceci !

ps: je ne connaissais pas "filtfilt" c'est super (si j'ai bien compris c'est comme filter mais ça enlève les soucis de dephasage!?)

[invite9c7554e3]

Juste une dernier question :

=> comment as tu choisi l'ordre du filtre FIR1 et de BUTTER ?

par exemple pourquoi 2000 et pourquoi pas 200 ?
de meme pour butterworth ?

=> en fait tu choisi ceci à l'oeil ?

au fait lorsque je lance ceci sur mes courbes ça me donne :

??? Error using ==> fir1 at 77
Frequencies must fall in range between 0 and 1.

Error in ==> differentsFiltres at 49
b=fir1(2000,fcoupure,'stop');

as tu une diée d'où ça vient ?

ps: au fait je n'arrive pas à t'envoyer des MP ta boite de rececption est pleine
(je voulais t'envoyer mon fichier de données)