Discussion :

[Vendy12]

Bonjour,

Je suis sur un projet Scilab et je dois appliquer un filtre passe bas à un signal. Je bloque complètement malgré des heures de recherche.

Mes données sont simple, il s'agit de 2000 mesures effectué en 20 seconde. J'ai donc un vecteur temps ou le pas est de 0.01 et un vecteur contenant les 2000 mesure.

Je ne m'y connais pas en traitement du signal, et j'ai pour consigne d'appliquer un filtre passe bas sur ce signal (qui a priori ici est de 100Hz).

Je vous remercie pour votre aide.

[phryte]

Bonjour.
Un exemple de PB d'ordre un (1/(1+tau*p) associé au bloqueur d'ordre zéro. (Il y a plusieurs sortes de Passe-bas en numérique) :

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
//EQUATION RECURRENTE DU PASSE-BAS D'ORDRE UN
//-------------------------------------------
T=0.01;// Période d'échantillonnage
tau=0.9;// Constante de temps
beta=exp(-T/tau);
alpha=1-beta;
xmin=1;xmax=2000;
//
//INITIALISATION
//--------------
e=ones(1,xmax);// Signal d'entrée : Echelon
s=0;enm1=0;
//
//EQUATION RECURRENTE
//--------------------
//
// Filtre
for i=xmin:xmax;
s(i)=alpha*e(i)+beta*enm1;
enm1=s(i);// Vieillissement
end;
//
//TRACE
//-----
xset("colormap",[1 0 0]);
i=[1:xmax];
plot(i,s(i));
xset("colormap",[0 1 0]);
xgrid();
xtitle("PASSE BAS ORDRE UN","nT","Sn");

Tu peux aussi utiliser la fonction de transfert avec filter :

Code :

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

1
2
3
4
5
6
num=alpha;
den=[1 -beta];
[y,zf] = filter(num,den,e,[0])
figure(1)
plot(i,y,'r',i,e)
xgrid()

[Vendy12]

Merci pour ta réponse.

J'ai plus ou moins compris mais sa dépasse le niveau que j'ai dans ce domaine. Du coup je suis passé par la transformée de fourier et j'ai sélectionné les première valeur de mon vecteur pour reconstruire un signal. J'imagine que c'est une alternative valable ?

[phryte]

Bonjour,

Cela dépend du but recherché.
Voici une solution simple pour filtrer un signal temporel :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
clear
T=0.01;// Période d'échantillonnage
tau=0.9;// Constante de temps
t=0:T:20;
e=sin(2*%pi*0.1*t);// Signal d'entrée
s=poly(0,'s');
Sys=syslin('c',[1/(tau*s+1)]);// Définition du filtre en analogique F(p)
FTZ=ss2tf(dscr(tf2ss(Sys),T));// Fonction de transfert en Z
SF=flts(e,tf2ss(FTZ));   // Sortie du filtre
clf();
plot(t,e,'b');// Tracé du signal d'entrée
plot(t,SF,'r');// Tracé du signal de sortie
xgrid()
xtitle("PASSE BAS ORDRE UN")

[Vendy12]

Merci pour ce code.

Je pense l'utiliser car avec la transformé de Fourier la forme du signal est la même, filtré, mais les valeur sont totalement différente ce qui ne permet pas de superposer les deux signaux.