Tracé filtrage cheby1

**ben1015** · 11/01/2011, 14h40

Bonjour,

Je souhaite réalisé un filtre de Chebyshev avec la commande cheby1 et ensuite filtré un signal de la forme w=Asin(wt)+Bsin(nwt)+... J'arrive à réaliser le filtre et à tracé son diagramme de Bode mais je n'arrive pas à filtrer le signal par la fonction de transfert.
Voici mon code :

Code :

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

1
2
3
4
5
6
w=Asin(wt)+Bsin(nwt)+...
[b,a] = cheby1(n,R,Wp,'bandpass','s');
o = logspace(-1,1,5000);
figure(2);
freqs(b,a,o);
Hf=freqs(b,a,o);

J'ai essayé de faire la fft de mon signal w pour pouvoir le multiplier par Hf et ensuite faire la transformée inverse mais cela ne fonctionne pas.

Code :

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

1
2
3
4
wf=fft(w);
yf=Hf.*wf;
yt=ifft(yf);
plot(t,yt);

Quel opération je dois effectuer pour pouvoir visualiser le signal w filtré?

Merci d'avance
Ben

**phryte** · 12/01/2011, 08h35

Bonjour.

Quel opération je dois effectuer pour pouvoir visualiser le signal w filtré?

help filter

Un essai en signal échantillonné :

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
clear
Fe=1000;
t=0:1/Fe:0.1;
n=2;
w=2*pi*15;
10; R= 0.5;
A=1;B=0.5;
Wp = 2*pi*[10 20]/Fe/2;
w=A*sin(w*t)+B*sin(3*w*t);
[b,a] = cheby1(n,R,Wp,'bandpass');
o = logspace(-1,1,5000);
figure(1);
freqs(b,a,o);
Hf=freqs(b,a,o);
wf=fft(w);
yt=ifft(real(wf));
figure(2)
plot(t,real(yt));
sf=filter(b,a,w);
grid
figure(3)
plot(t,w,t,sf,'r')
grid

**ben1015** · 15/01/2011, 10h57

Bonjour,

J'ai regardé ton code et la fonction filter, j'arrive à avoir un résultat mais mon problème est que mon filtre est un passe-bas (voir diagramme de Bode) et non un passe-bande, pourtant je spécifie bien "passband", si quelqu'un voit le problème? Donc forcément mon signal de sortie est celui de 5 et 20 Hz alors que je devrais avoir un signal sinusoïdal de 20Hz.

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
clc 
clear all
pi = 3.14159;
f = 5;
Fe=1000;
t=0:1/Fe:1;
x = sin(2 * pi * 1*f*t);
y = sin(2* pi * 4*f*t);
z = sin(2 * pi * 7*f*t);
w = x+y+z;
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Wp = 2*pi*[10 30]/Fe/2;
n = 4;
R = 1;
[b,a] = cheby1(n,R,Wp,'bandpass');
o = logspace(-1,1,5000);
figure(1);
freqs(b,a,o);
 
s=filter(b,a,w);
figure(2);
plot(t,s);