Design de filter avec la méthode de fenêtrage
Bonjour tout le monde :)
Je dois designer :
-un filtre passe haut telle que la fréquence de coupure est égale à 0.2 , on veut que N=15 en utilisant la fenêtre de Blackman.
-un filtre passe bande, telle que la bande est [20,30] Hz, S=100Hz, N=20 avec la fenêtre de Hanning.
Bon pour commencer.
1) filtre passe haut
Je change les spécifications du filtre passe haut en filtre passe bas (quelle méthode ? la transformation bilinéaire ? ) puis ayant les spécifications de mon passe-bas
Je cherche hw(n) = h(n).w(n)
w(n) étant la fenêtre de Blackman et h(n) la fonction représentant un passe-bas.
Il me semble que la meilleur représentation pour le passe bas est du type
2fc sinc (2pi fc n ) ?
Ensuite je calcul hw(f)= H(f)*W(f) ?
ensuite je sais pas quoi faire :'( ensuite
Code:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
|
clear all
N=15;
N1=1024
n=0:1:15;
fc=0.2
L=N+1;
fe3=1/(n(2)-n(1));
f=(0:N1-1)*fe3/N1;
w=blackman(L);
h=2*pi*fc*sinc(2*pi*fc*n)
figure(1)
plot(w);
H=fft(h,1024);
W=fft(w,1024);
s=conv(H,W);
figure(2)
plot(f,abs(H))
figure(3)
plot(f,abs(W))
figure(4)
plot(abs(s)) |
