Bonjour tous,

j'ai fais un peu le tour des discussions sur le forum et je pense avoir compris comment faire un filtre coupe bande de manière assez propre:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
 
clear all;close all;clc;
t=0:0.01:1;
f=sin(2*pi*5*t)+sin(2*pi*12*t);
nbPoint=length(t);
fe=nbPoint/(t(end)-t(1))
[b,a]=butter(2,[10 13]/(fe/2),'stop')
y=filtfilt(b,a,f);
figure (1)
plot(t,f,'b',t,y,'r')
ce bout de code fonctionne bien mais si j'ai un pas de temps non pas
de 0.01 mais de 0.1 alors je vais avoir un problème :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
The cutoff frequencies must be within the interval of (0,1).
si j'ai bien compris cela veut dire que je ne peux pas filtrer des frequences supérieurs à la frequence d'echantillonage divisé par 2 ?

Question 1
du coup quel recours j'ai ? si j'augmente l'ordre du filtre ça va changer quelque chose ? je vous demande ceci car dans un cas réel j'ai justement ce problème.. j'ai une frequence d'acquisition de 2Hz et je veux filtrer la frequence 50 Hz.

Question 2
ma deuxieme question porte que certains choix :
Code : Sélectionner tout - Visualiser dans une fenêtre à part
filtfilt
http://www.mathworks.fr/fr/help/sign.../filtfilt.html
j'ai utilisé 'filtfilt' car d'après ce que j'ai compris c'est la même chose que 'filter' sauf que ça compense le dephasage qui peut etre introduit par le filtrage ?

Code : Sélectionner tout - Visualiser dans une fenêtre à part
butter
http://www.mathworks.fr/fr/help/signal/ref/butter.html
butter est une technique qui permet de donner les paramètres du filtre mais si j'utilise "fir1" : http://www.mathworks.fr/fr/help/signal/ref/fir1.html
plutôt que "butter" qu'es ce que ça va changer ?

Code : Sélectionner tout - Visualiser dans une fenêtre à part
ordre
ma dernière question porte sur l'ordre du filtre: comment le choisit ton en pratique ? plus l'ordre est elevé plus on va filtrer nos pics de manière "etroite" ?

Pourriez vous m'aidez s'il vous plait ?