数字信号处理实验报告——用窗函数法设计-FIR数字滤波器.docx
实验四用窗函数法设计FIR数字滤波器一、实验目的(1)掌握用窗函数法设计FlR数字滤波器的原理和方法。(2)熟悉线性相位FIR数字滤波器特性。(3)了解各种窗函数对滤波特性的影响。二、实验原理、滤波器的理想频率响应函数为Hd(>),那么其对应的单位脉冲响应为hd(n)二白瓦32咽窗函数设计法的根本原理是用有限长单位脉冲响应序列h(n)逼hd(n)t>由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数。w(n)将hd(n)截断，并进行加权处理：h(n)=hcj(n)w(n)h(n)就作为实际设计的FIR数字滤波器的单位脉冲响应序列，其频率响应函数HG")为N-IH(ej)=Yjh(n)ejn=G用窗函数法设计的滤波器性能取决于窗函数w(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择适宜的窗函数类型和窗口长度N。三、实验内容1 .MATALAB程序任务一Nl=15;N2=33;bl=firl(14,l4,hanning(15);b2=fir1(32,1/4,hanning(33);Hl,W=freqz(bl,l);HLdb=20*1Og10(abs(Hl);magHl=abs(Hl);phaHl=angle(Hl);H2,W=freqz(b2,l);H2.db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);figure(l);subplot(2,1,1);stem(b1);title('N=15时,汉宁窗得到的FIR滤波器脉冲响应subplot(2,1,2);stem(b2);title('N=33时,汉宁窗得到的FIR滤波器脉冲响应)；figure(2);subplot(3,1,1);plot(W/pi,magH1);title('N=15时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,Hl_db);title('N=15时汉宁窗得到的FIR滤波器的对数幅频特性)SUbPIOt(3,3);PIot(Wpi,phaH1);title('N=15时汉宁窗得到的FIR滤波器的相频特性)figure(3);subplot(3,1,1);plot(W/pi,magH2);title('N=33时汉宁窗得到的FIR滤波器的幅频特性)subplot(3,1,2);plot(W/pi,H2_db);title('N=33时汉宁窗得到的FIR滤波器的对数幅频特性)subplot(3,1,3)；plot(Wpi,phaH2);title('N=33时汉宁窗得到的FIR滤波器的相频特性)任务二N=33;b1=firl(32,l4,boxcar(33);b2=firl(32,l4,hanning(33);b3=fir1(32,1/4,bartlett(33);b4=firl(32,l4,blackman(33);Hl,W=freqz(bl,l);HLdb=20*1Og10(abs(H1);magHl=abs(Hl);phaHl=angle(Hl);H2,W=freqz(b2,l);H2_db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);H3,W=freqz(b3,l);H3_db=20*log10(abs(H3);magH3=abs(H3);phaH3=angle(H3);H4,W=freqz(b4,l);H4_db=20*log10(abs(H4);magH4=abs(H4);phaH4=angle(H4);figure(l);subplot(2,2,1);stem(b1);titleC矩形窗得到的FIR滤波器脉冲响应)subplot(2,2,2);stem(b2);title，汉宁窗得到的FIR滤波器脉冲响应)SUbPIOt(2,2,3);Stem(b3);title。三角窗得到的FIR滤波器脉冲响应)SUbPlOt(2,2,4);Stem(b4);title，布拉克曼窗得到的FIR滤波器脉冲响应)figure(2);subplot(2,2,1);plot(W/pi,magH1);titleC矩形窗得到的FIR滤波器幅频特性)SUbPIOt(2,2,2);PIOt(Wpi,magH2);title，汉宁窗得到的FIR滤波器幅频特性)SUbPlOt(2,2,3);PlOt(Wpi,magH3);titleC三角窗得到的FIR滤波器幅频特性)SUbPIOt(2,2,4);PIOt(Wpi,magH4);title，布拉克曼窗得到的FIR滤波器幅频特性)figure(3);subplot(2,2,1);plot(W/pi,Hl_db);出世矩形窗得到的FIR滤菠器相频特性)subplot(2,2,2);plot(W/pi,H2_db);title,汉宁得到的FIR滤波器相频特性)subplot(2,2,3)；plot(W/pi,H3_db);出世三角窗得到的FIR滤菠器相频特性)SUbPIOt(2,2,4);Plot(W/pi,H4_db);title,布拉克曼得到的FIR禧波器相频特性)figure(4);subplot(2,2,1);plot(Wpi,phaH1);title,矩形窗得到的FIR滤波器相频特性)SUbPIOt(2,2,2);PIOt(Wpi,phaH2);titleC汉宁窗得到的FIR滤波器相频特性)SUbPIOt(2,2,3);PIOt(Wpi,phaH3);titleC三角窗得到的FIR滤波器相频特性)SUbPIOt(2,2,4);PlOt(Wpi,phaH4);title，布拉克曼窗得到的FIR滤波器相频特性)2 .实验波形图任务一0.30.250.20.150.10050N=15时,汉宇苜得到的Fl磁波器肽冲响应0250.2015N=33时双宁窗得到的Fl磁波器脉冲响应0.33035N=15时汉宁亩得到的FIR滤波器的幅频特性任务二布拽克曼由得到的Fl礴源器脉冲雨应0.30250.20.150.10.05四、实验结论1.N的大小决定了窗谱的主瓣宽度，N越大，窗谱的主瓣宽度越大2 .最小阻带衰减只有窗行决定，不受N的影响，过渡带宽度与N和窗形都有关，N越大，过渡带宽越小3 .由实验可知滤波特性：布拉克曼窗汉宁窗三角窗矩形窗五、思考题(1)如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。答：根据A°二°s一%计算出所要设计的滤波器的过渡带查表找到满足的阻带最小衰减AS所需的窗函数利用已选的窗函数的过渡带宽除以AG得到低通滤波器单位脉冲响应的长度N根据/)二吊("()求出h(n)Af-I再根据H(ei)=2%求出响应函数H(ei)n=0如果要求用窗函数法设计带通滤波器，且给定上、下边带截止频率为必和0)2,试求理想带通的单位脉冲响应hd(n)o答：由=竺言算出纭，求出滤波器过渡带宽Ag=q-牡和滤波器长度N,最后由hd(n)=7-_cSinl<.(一理声c或()=幺,=7求出儿()，其中丁=-n-)2
