数字带阻滤波器
这里仍然选择凯泽窗来完成带阻滤波器的设计。
带阻滤波器的性能指标为：
flp=200; fup=3500; flc=300; fuc=3300; As=46; Ap=3; Fs=22050;
具体程序如下：
>>
As=46; Ap=3; Fs=22050;
beta=0.1102*(As-8.7); %由阻带最小衰减指标确定形状参数beta
wp1=2*pi*200/Fs; ws1=2*pi*300/Fs;
wp2=2*pi*3500/Fs;ws2=2*pi*3300/Fs;
wc1=(wp1+ws1)/2; wc2=(wp2+ws2)/2;wc=[wc1,wc2];
BW=ws1-wp1; N2=12.8*pi/BW; N3=ceil(N2);
b=fir1(N3,wc/pi,'stop',kaiser((N3+1),beta));
figure;[h,w]=freqz(b,1,512);
plot(w/pi,20*log10(abs(h)));grid;axis([0,0.5,-180,20]);
xlabel('归一化频率');ylabel('幅度/dB');
title('凯泽窗设计FIR数字带阻滤波器频率响应曲线');
figure;freqz(b,1,512,Fs);axis([0,5000,-180,20]);
title('凯泽窗设计FIR数字带阻滤波器幅频与相频曲线');
运行结果如图20所示：
评论0