fir低通滤波器设计(完整版)

电子科技大学信息与软件工程学院学院标准实验报告（实验）课程名称数字信号处理电子科技大学教务处制表电孑科技大学实验报告学生姓名： 学号：指导教师：实验地点：实验时间：14-18一、 实验室名称：计算机学院机房二、 实验项目名称：fir低通滤波器的设计三、 实验学时：四、 实验原理：FIR滤波器FIR滤波器是指在有限范围内系统的单位脉冲响应h[k]仅有非零值的滤波器。M阶FIR滤波器的系统函数H(z)为H(z)=产h[k]z-kk=0其中H(z)是少的M阶多项式，在有限的z平面内H(z)有M个零点，在z平面原点z=0有M个极点.FIR滤波器的频率响应H(eJQ)为H(ejQ)=费h[k]e-jk。k=0它的另外一种表示方法为H(ejQ)=H(ejQ)1ej^(^)其中I""')和'(。)分别为系统的幅度响应和相位响应。若系统的相位响应'(。)满足下面的条件NQ)=—aQ即系统的群延迟是一个与。没有关系的常数a，称为系统H(z)具有严格线性相位。由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR滤波器设计中一般使用广义线性相位。如果一个离散系统的频率响应H(幻Q)可以表示为H(ejQ)=A(Q)ej(-aQ+&)其中a和3是与。无关联的常数，A(Q)是可正可负的实函数，则称系统是广义线性相位的。如果M阶FIR滤波器的单位脉冲响应h[k]是实数，则可以证明系统是线性相位的充要条件为h[k]=±h[M—k]当h[k]满足h[k]=h[M-k],称h[k]偶对称。当h[k]满足h[k]=-h[M-k],称h[k]奇对称。按阶数h[k]又可分为M奇数和M偶数，所以线性相位的FIR滤波器可以有四种类型。窗函数法设计FIR滤波器窗函数设计法又称为傅里叶级数法。这种方法首先给出"*0)，H*^)表示要逼近的理想滤波器的频率响应，则由IDTFT可得出滤波器的单位脉冲响应为1"h[k]=^—JH(ejQ)ejkqdQ—由于是理想滤波器，故"d伙]是无限长序列。但是我们所要设计的FIR滤波器，其h[k]是有限长的。为了能用FIR滤波器近似理想滤波器，需将理想滤波器的无线长单位脉冲响应"d伙]分别从左右进行截断。当截断后的单位脉冲响应"d伙]不是因果系统的时候，可将其右移从而获得因果的FIR滤波器。另一种设计方案是将线性相位因子幻(ewQ+P)加入到理想滤波器的频率响应中，然后利用IDTFT计算出hd伙]后，取hd伙]在OMkMM范围的值为FIR滤波器单位脉冲响应。理想滤波器的频率响应Hd(e；Q)和设计出的滤波器的频率响应Hd(°j)的积分平方误差定义为&2=2—J|h^(eg)—H(e0)|2dQ-ne2也可以表示为£2=2°\h^[k]—h[k]2k=—°=2\hd[k]|2+!^|hd[k]—h[k]|2+2\hd[k]|2k=-° k=0 k=M+1上式中的第一项和第三项与所设计出的滤波器参数是没有关系的，为了使上式中的第二项达到最小，可选择h[k]=h*],0<k<M所以用上面的方法得出的滤波器是在积分平方误差最小意义下的最佳滤波器。Gibbs现象就是理想滤波器的单位脉冲响应hd伙]截断获得的FIR滤波器的幅度函数A(Q)在通带和阻带都呈现出振荡现象。随着滤波器阶数的增加，幅度函数在通带和阻带振荡的波纹数量也随之增加，波纹的宽度随之减小，然而通带和阻带最大波纹的幅度与滤波器的阶数M无关。窗函数的主瓣宽度决定了Hd@j)过渡带的宽度，窗函数长度N增大，过渡带减小。五、 实验目的：掌握FIR低通滤波器设计以及在MATLAB上的实现。六、 实验程序及仿真：disp(-利用窗函数法设计FIR滤波器')；fp=100；%通带截止频率fs=200;%阻带截止频率Fs=1000；%采样频率Wc=(fs+fp)/Fs%确定截止频率m=(0.3*fp)/Fs;%定义过度带宽M=round(8/m);%定义窗函数的长度N=M-1;%定义滤波器的阶数box=fir1(N,Wc,boxcar(N+1));%使用fir1函数设计滤波器ham=fir1(N,Wc,hamming(N+1));han=fir1(N,Wc,hanning(N+1));black=fir1(N,Wc,blackman(N+1));[h1,f1]=freqz(box,1,1024);%滤波器的幅频特性图[h2,f2]=freqz(ham,1,1024);[h3,f3]=freqz(han,1,1024);[h4,f4]=freqz(black,1,1024);plot(f1*Fs/(2*pi),20*log10(abs(h1)),':',f2*Fs/(2*pi),20*log10(abs(h2)),'--',f3*Fs/(2*pi),20*log10(abs(h3)),'-.',f4*Fs/(2*pi),20*log10(abs(h4)));legend('矩阵窗'，'哈明窗'，'汉宁窗'，'布莱克曼窗');grid;xlabel('频率/赫兹')；ylabel('增益/分贝')；title('滤波器的增益响应,);axis([0500-20020]);disp。’利用鲨函故沽渔轩FIEi?直聂)：段品童止建天fs-叫。遇阻希琴止餐本FE=Laoo.sSftfi!l曾皓fMH,-Fs*有走IS止耳车PhtaundI日占l：I:羿走：乂定画熟的枝度旺L【:■表必旧裱在的用可boK=firI torci.ril+i)I.*irJ®feiRd-整波露h»-fLrlCW/fc?hwnin：〔H+lIJ:h«i"Juf1I.M；,fc；!hann.Lnf；I.H4-LJ.1；black=firl.ClT,Vc-.,bl-k-:-U>an ;D4.fL]=f戏眠IdLr]0241做St吾的国羲帮性13D12).f2]=ft:eq.zI1IU+L>]1)2^1.iharuL>]024):Di4,id=fEeqj:tbla£k7h1叩5.Pint-If1*7=/<2*pi)P2O*lnELi<abs：(hLIJ,'i: <2*pi)F2t>*lnELOIa^s-IlJ))/—,f3*1?3/l.2-=piI?^O+losiOCate<h3)),'\£4>FV 10ihi)]I.loEeadr拒傅冒・/咆用石：/醍宇苴/■rESMta").gridMlab«LlJKiS/帷)jyLfitosiirt5&/#JjqrltiEl*C占检嚣的号月响成’)理5QQ-2002d]):

』,回FileEditViewInsertToolsDesktopWindowHelp□蛉I nioi^n滤波器的增益响应50 100 150 200 260 300 350 400 450 500频率f赫兹 矩阵窗滤波器的增益响应50