fir2.asv

FIR滤波器设计

wc=0.25*pi;N=33;
alpha=(N-1)/2;
n=0:1:N-1;
m=n-alpha+eps;
hd=sin(wc*m)./(pi*m);
wd1=hanning(N)';b1=hd.*wd1;
wd2=hamming(N)';b2=hd.*wd2;
wd3=hanning(N)';b3=hd.*wd3;
wd4=hanning(N)';b4=hd.*wd4;
[H1,w]=freqz(b1,1);
[H2,w]=freqz(b2,1);
[H3,w]=freqz(b3,1);
[H4,w]=freqz(b4,1);
subplot(3,1,1)
plot(w,abs(H1),':',w,abs(H2),':',w,abs(H3),':',w,abs(H4),':')
legend('汉宁窗低通滤波器','海明窗低通滤波器','布来克曼窗低通滤波器','凯泽窗低通滤波器')
subplot(3,1,2)
plot(w/pi,20*log10(abs(H1)),':',w/pi,20*log10(abs(H2)),':',w/pi,20*log10(abs(H3)),':',w/pi,20*log10(abs(H4)),':')
legend('汉宁窗幅度响应（dB）','海明窗幅度响应（dB）','布来克曼窗幅度响应（dB）','凯泽窗幅度响应（dB）')
subplot(3,1,3)
plot(n,b1,':',n,b2,':',n,b3,':',n,b4,':')
legend('汉宁窗h(n)','海明窗h(n)','布来克曼窗h(n)','凯泽窗h(n)')