实验11-用MATLAB设计FIR数字滤波器

;ylabel('h(n)');subplot(2,2,2);stem(n,windows);axis([0,N,0,1.1]);title('窗函数特性');xlabel('n');ylabel('wd(n)');subplot(2,2,3);plot(w/pi,db);axis([0,1,-100,2]);title('幅度频率响应');xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');set(gca,'XTickMode','manual','XTick',[0,ws/pi,wp/pi,1]);set(gca,'YTickMode','manual','YTick',[,-40-20,-3,0]);gridsubplot(2,2,4);plot(w/pi,pha);axis([0,1,-4,4]);title('相位频率响应');xlabel('频率(单位:\pi)');ylabel('\phi(\omega)');set(gca,'XTickMode','manual','XTick',[0,ws/pi,wp/pi,1]);set(gca,'YTickMode','manual','YTick',[-pi,0,pi]);gridN=67Rp=0.2321As=324、选择合适的窗函数设计一个FIR数字带通滤波器,要求:fp1=3.5kHz,fp2=6.5kHz,Rp=0.05dB;fs1=2.5kHz,fs2=7.5kHz,As=60dB。滤波器采样频率Fs=20kHz。描绘该滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相频响应曲线。分析:根据设计指标应选择布莱克曼窗。程序清单如下:function[db,mag,pha,grd,w]=freqz_m(b,a);[H,w]=freqz(b,a,1000,'whole');