2022年fir低通滤波器设计 .pdf
电子科技大学信息与软件工程学院学院标 准 实 验 报 告（实验）课程名称数字信号处理电子科技大学教务处制表名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 6 页 - - - - - - - - - 电 子 科 技 大 学实验报告学生姓名：学 号：指导教师：实验地点：实验时间： 14-18 一、实验室名称：计算机学院机房二、实验项目名称： fir 低通滤波器的设计三、实验学时：四、实验原理：1.FIR 滤波器FIR 滤波器是指在有限范围内系统的单位脉冲响应hk 仅有非零值的滤波器。M阶 FIR 滤波器的系统函数H(z) 为0()MkkHzh kz其中 H(z) 是kz的 M阶多项式，在有限的z 平面内 H(z) 有 M个零点，在z平面原点 z=0 有 M个极点 . FIR 滤波器的频率响应()jH e为0() MjjkkH eh k e它的另外一种表示方法为()()()jjjH eH ee名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 6 页 - - - - - - - - - 其中()jH e和()分别为系统的幅度响应和相位响应。若系统的相位响应()满足下面的条件()即系统的群延迟是一个与没有关系的常数，称为系统 H(z) 具有严格线性相位。由于严格线性相位条件在数学层面上处理起来较为困难，因此在FIR滤波器设计中一般使用广义线性相位。如果一个离散系统的频率响应()jH e可以表示为()()()jjH eAe其中和是与无关联的常数，()A是可正可负的实函数，则称系统是广义线性相位的。如果 M阶 FIR 滤波器的单位脉冲响应hk 是实数，则可以证明系统是线性相位的充要条件为 h kh Mk当 hk 满足 hk=hM-k,称 hk 偶对称。当hk 满足 hk=-hM-k,称 hk奇对称。按阶数hk 又可分为 M奇数和 M偶数，所以线性相位的FIR 滤波器可以有四种类型。2.窗函数法设计 FIR 滤波器窗函数设计法又称为傅里叶级数法。这种方法首先给出()jdHe，()jdHe表示要逼近的理想滤波器的频率响应，则由IDTFT 可得出滤波器的单位脉冲响应为1 ()2jjkddh kHeed由于是理想滤波器，故 dh k是无限长序列。但是我们所要设计的FIR 滤波器，其 hk 是有限长的。为了能用FIR 滤波器近似理想滤波器，需将理想滤波器的无线长单位脉冲响应 dh k分别从左右进行截断。当截断后的单位脉冲响应 dh k不是因果系统的时候，可将其右移从而获得因果的FIR 滤波器。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 6 页 - - - - - - - - - 另一种设计方案是将线性相位因子( 0.5)jMe加入到理想滤波器的频率响应中，然后利用IDTFT 计算出 dh k后，取 dh k在 0kM范围的值为FIR 滤波器单位脉冲响应。理想滤波器的频率响应()jdHe和设计出的滤波器的频率响应()jdHe的积分平方误差定义为221()()2jjdHeHed2也可以表示为22 dkh kh k122201 MdddkkkMh kh kh kh k上式中的第一项和第三项与所设计出的滤波器参数是没有关系的，为了使上式中的第二项达到最小，可选择 ,0dh khkkM所以用上面的方法得出的滤波器是在积分平方误差最小意义下的最佳滤波器。 Gibbs现象就是理想滤波器的单位脉冲响应 dh k截断获得的FIR 滤波器的幅度函数()A在通带和阻带都呈现出振荡现象。随着滤波器阶数的增加，幅度函数在通带和阻带振荡的波纹数量也随之增加，波纹的宽度随之减小，然而通带和阻带最大波纹的幅度与滤波器的阶数M 无关。窗函数的主瓣宽度决定了()jdHe过渡带的宽度，窗函数长度N增大，过渡带减小。五、实验目的：掌握FIR 低通滤波器设计以及在MATLAB 上的实现。六、实验程序及仿真：disp( 利用窗函数法设计FIR 滤波器 ); fp=100;%通带截止频率fs=200;%阻带截止频率Fs=1000;%采样频率名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 6 页 - - - - - - - - - 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(0 500 -200 20); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 6 页 - - - - - - - - - 七实验结论通过设计实现 FIR 低通滤波器，并不断的借鉴学习matlab的知识，编辑改善matlab程序，从而得到了滤波器较为理想的输出状态图。指导教师签字：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 6 页 - - - - - - - - -
