filtre Passe Bas

[invite17ae4981]

Bonjour tout le monde,

Je cherche a créer un filtre passe-bas sous matlab pour les données suivantes:
fréquence d'échantillonnage fe=5000Hz, Nombre échantillonnée N=50000,
j'ai trouvé des routines sur des forums mais comme je suis débutant sur matlab je comprend pas grande chose
comme par exemple ps :

%
% [num,den] = butter(n,Omega)
% n = nombre de poles
% Omega = cutoff frequency
%

%
% filter(num,den,X)
% where X is the data we want to filter
%

%test

n = 2;
Omega = 0.10; %Hz
[num,den] = butter(n,Omega);

X=[];W=[];
for i=-20:0.05:20
X=[X, atan(i)+rand(1)];
W=[W,atan(i)];
end

Y=-20:0.05:20;
Z=filter(num,den,X);

plot(Y,X,Y,Z,Y,W);

comment estimer le numérateur et le dénominateur ? c'est quoi la fréquence de coupure ?
-----

[lou_ibmix_xi]

comment estimer le numérateur et le dénominateur ? c'est quoi la fréquence de coupure ?

C'est bien la fonction "butter" qui détermine le numérateur et le dénominateur (les coefficients "auto-regressive" et "moving-average") pour un filtre de type "butterworth", pour un ordre et une fréquence de coupure donnés. "help butter". Si tout ça c'est du chinoix, il est alors grand temps de relire les cours de traitements du signal, la partie sur le filtrage...

[invite17ae4981]

Merci lou, je comprend ce que tu m'écris, bon je vais essayer de faire mon mieux pour m'améliorer
merci

[invite17ae4981]

Bonjour,
qlq'un peut m'expliquer pourquoi ça marche pas mon programme pour filtrer un signal de fechont=5000hz et N=50000

a=load('7.txt');
x=a(:,2);
N=500000;
dt=0.0002;
t=0:dtN-1)*dt;
[num,den] = butter(n,Omega)
% n = nombre de poles
% Omega = cutoff frequency

% filter(num,den,X)
% where X is the data we want to filter
%

%test

n = 2;
Omega = 0.10; %Hz
[num,den] = butter(n,Omega);

Z=filter(num,den,X);

plot(t,Z);

là je comprend pas pourquoi ça marche pas !

merci d'avance

[A voir en vidéo sur Futura]