FIR数字滤波器设计.doc

吉林建筑大学电气与电子信息工程学院数字信号处理课程设计报告设计题目： FIR数字滤波器设计 专业班级： 学生姓名： 学 号： 指导教师： 设计时间： 2015.1.122015.1.16 教师评语：成绩 评阅教师 日期 目 录 摘要3第1章绪论41.1设计的作用.51.2课程设计基本要求.51.3设计内容.5第2章FIR数字低通滤波器的原理.52.1设计原理.62.2设计方法62.3窗函数设计原理62.4各种窗函数比较72.5 FIR变换方法的原理7第3章FIR数字滤波器设计分析83.1用窗函数设计滤波器.83.2数字低通滤波器MATLAB编程.93.3哈明窗与汉宁窗优缺点.13总结.14参考文献.15 摘要摘要：本论文首先介绍了FIR数字滤波器的基础知识，然后讲述了本论文对于FIR数字滤波器的设计方法。设计思想为：以窗函数设计方法为基础，选定一个合理的选频滤波器，然后截取它的脉冲响应得到线性相位。仿真工具应用MATLAB软件设计FIR数字滤波器，并对所设计的滤波器进行了分析总结，可以得到较好的滤波效果。关键字：窗函数、数字滤波器、MATLAB.Abstract : this paper first introduced the basic knowledge of FIR digital filters, and then introduced for the design method of FIR digital filters, the design thought for to the window function design method as the foundation, the selected a ideal frequency selective filter and then intercept it get the pulse response of linear phase MATLAB simulation tool software design FIR digital filter ,and the design of filter is summarized can get good effect of filter.Key word :window function method MATLAB digital filter 第1章 绪 论1.1设计的作用课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。本次课程设计一方面通过MATLAB仿真设计内容，使我们加深对理论知识的理解，同时增强其逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充。1.2课程设计基本要求通过课程设计各环节的实践，应使学生达到如下要求：1.掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器以及窗函数法设计FIR数字滤波器的原理、具体方法及计算机编程；2.观察双线性变换法、脉冲响应不变法及窗函数法设计的滤波器的频域特性，了解各种方法的特点；3.用MATLAB画出三种方法设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查结果是否满足要求。1.3设计内容（1）设计一线性相位FIR数字低通滤波器，截止频率 , 过渡带宽度 阻带衰减.2）设计一线性相位FIR数字低通滤波器，截止频率 , 过渡带宽度 ，阻带衰减。第2章FIR数字低通波器的原理2.1 FIR数字低通滤波器的设计原理 FIR数字滤波器传统的设计方法有窗函数法、频率抽样法和等波纹逼近法，用窗函数设计FIR数字滤波器就是用有限长的脉冲相应逼近序列。 基本设计思想：首先选定一个理想的选项滤波器，然后截取它的脉冲响应得到线性相位。 滤波器：一种用于消除干扰杂讯的器件，将输入或输出经过过滤而得到纯净的直流电。对特定频率的频点或该频点以外的频率进行有效滤除的电路，滤波器其功能就是得到一个特定频率消除一个特定频率。就是允许某一部分频率的信号顺利的通过，而另一部分频率的信号受到较大的抑制，它实质上是一个选项电路。2.2 FIR数字低通滤波器设计方法 FIR数字滤波器的设计方法有窗函数法、频率采样法和基于 firls 函数和remez函数的最优化方法。MATLAB语言中的数字信号处理工具箱，提供的滤波函数有firl（），此函数以经典的方法实现加明窗线性相位FIR数字滤波器设计，可以设计低通、高通、带通和带阻滤波器；fir2函数设计的FIR滤波器，其滤波的频率特性由矢量f和m决定，f和m分别为滤波器的期望幅频响应的频率相量和幅值相量。2.3 用窗函数法设计原理 窗函数法比较简单，可应用现成的窗函数公式，在技术指标要求不高时世比灵活，窗函数法也称之为傅里叶级数变换法。FIR数字滤波器的设计首先给出的要求的理想滤波器的频率响应Hd()，设计一个FIR数字滤波器的频率响应H()，去逼近理想的滤波响应Hd()。然而，窗函数设计FIR数字滤波器是在时域进行的，因而必须由理想的频率响应Hd()推导出对应的hd（n），再设计一个FIR数字滤波器的单位取样响应h（n）去逼近hd（n）。设计过程如下：Hd( hd（n） h（n） H() （2-3-1） 加窗的作用是把理想滤波器的无限长脉冲响应hd（n）乘以窗函数w（n）来产生一个被截断的脉冲响应，即h（n）= hd（n）w（n）并且对频率响应进行平滑。FIR滤波器产生单位冲击响应h（n）；其单位冲激响应h（n）是无限长（1nN-1）系统函数为： H（z）= （2-3-2）在有限Z平面有（N-1）个零点，而它的（N-1）个极点均位于原点Z=0处。FIR滤波器线性相位的特点：如果FIR滤波器的单位抽样响应h（n）为实数，而满足以下任一条件：偶对称h（n）=h（N-1-n）奇对称h（n）=h-（N-1-n），其对称中心在n=（N-1）/2处，则滤波器具有准确的线性相位。窗函数应使设计的滤波器：（1）频率特性的主瓣宽度应该尽量窄，且尽可能将能量集中在主瓣内；（2）窗函数频率特性的旁瓣w趋于的过程中，其能量迅速减小为零。在实际工程中常用的窗函数有五种，即矩形窗、三角窗、汉宁窗、海明窗和凯泽窗。窗函数类型旁瓣峰值过渡带宽度阻带最小衰减近似值精确值矩形窗-13-21三角窗-25-25汉宁窗-31-44哈明窗-41-53 表2-32.4各种窗函数比较在各种窗函数中，三角窗阻带最小衰减最小，其次是汉宁窗，再次是哈明窗，布莱克曼窗阻带衰减最大。主瓣宽度方面，三角窗，汉宁窗以及哈明窗一样，布莱克曼窗比以上三种的主瓣要宽，表现在用布莱克曼窗的滤波器的过渡带比用其它三种窗函数要宽，主瓣高度与第一旁瓣高度的差值方面，三角窗高度差最小、其次是汉宁窗，再次是哈明窗，三角窗和哈明窗比较缓慢，汉宁窗相对前两种要快，布莱克曼窗很快旁瓣幅度汉宁窗和哈明窗比较缓慢，汉宁窗相对前两种要快，布莱克曼窗很快，旁瓣幅度汉宁窗和哈明窗较大，三角形窗和布莱克曼窗较小。2.5 FIR变换方法的原理FIR滤波器目前常用的设计方法有窗函数法和频率采样法，窗函数法是从时域进行设计，而频率采样法是从频域进行设计。窗函数法由于简单、物理意义清晰，因而得到了较为广泛的应用。窗函数法设计的基本思想是：首先根据技术指标要求，选取合适的阶数N窗函数的类型w(n)，使其幅频特性逼近理想滤波器幅频特性。其次，因为理想滤波器的是无限长的，所以需要对进行截断，数学上称这种方法为窗函数法。简而言之，用窗函数法设计FIR滤波器是在时域进行的，先用傅里叶变换求出理想滤波器单位抽样相应，然后加时间窗w(n)对其进行截断，以求得FIR滤波器的单位抽样响应h(n) 。第3章FIR数字滤波器仿真分析3.1用窗函数设计滤波器设计步骤（1）根据对阻带衰减及过渡带的指标要求，选择窗函数的类型，并估计窗口长度N。先按照阻带衰减选择窗函数类型。原则是在保证阻带衰减满足要求的情况下，尽量选择主瓣窄的窗函数。然后根据过渡带宽度估计窗口长度N。待求滤波器的过渡带宽度近似等于窗函数主瓣宽度，且近似与窗口长度N成反比，A取决于窗口类型，例如，矩形窗的,哈明窗的等，参数A的近似值和精确取值参考表3-1。（2）构造希望逼近的频率响应函数，即= 所谓的“标准窗函数法”就是选择为线性相位理想滤波器（理想低通、理想高通、理想带通、理想带阻）。以低通滤波器为例，应满足：= 理想滤波器的截止频率近似位于最终设计的FIRDF的过渡带的中心频率点，幅度函数衰减一半（约-6db）。所以如果设计指标给定通带边界频率的阻带边界频率和，一般取（3）计算。如果给出待求滤波器的频响函数为，那么单位脉冲响应应用下式求出：=如果较复杂，或者不能用封闭公式表示，则不能用上式求出。我们可以对从w=0和w=采样M点，采样值，k=0,1,2,3，M-1，进行M点IDFT(IFFT)得到：根据频域采样理论，与应满足如下关系：=因此，如果M选得较大，可以保证在窗口内有效逼近。=为保证线性相位特性，。（4）加窗得到设计结果：h（n）=。3.2数字低通滤波器MATLAB编程及幅频特性曲线 3.2.1窗函数法的MATLAB设计函数简介（1）firl 是用窗函数设计线性相位FIR数字滤波器的工具箱函数，实现线性相位FIR数字滤波器的标准窗函数法设计。这里所谓的“标准”，只是指在设计的低通、带通、高通、带通、和带阻FIR滤波器时，分别表示相应的线性相位理想低通、高通、带通、和带阻FIR滤波器的频率响应函数。因而将所设计的滤波器的频率响应称为标准频率响应。（2）fir2为任意形状幅度特性的窗函数法设计函数，用fir2设计时，可以指定任意形状的，它实质是一种频率采样法与窗函数法的综合设计函数。主要用于设计幅度特性形状特殊的滤波器。用help命令查阅其调用格式及调用参数的含义。（3）汉宁窗(Harming Window)升余弦窗窗函数：频率响应函数：当N1时，N-1N幅度函数：参数为：n=-31dB；Bg=8/N；s=-44dB。图(3)汉宁窗的滤波器频率响应 图 (3) 汉宁窗函数的脉冲响应（4）哈明窗(Hamming Window)函数窗函数：频谱函数：幅度函数：参数为：n=-41dB；Bg=8/N；s=-53dB。图(4)哈明窗的滤波器频率响应 图 (4) 哈明窗函数的脉冲响应3.2.2 MATLAB简介MATLAB是一种面向科学与工程计算的高级语言。它的语句比Basic、C等 更加接近我们所熟悉的数学公式的表达方式，所以它的编程效率高，易学易懂。MATLA提供了大量的图形功能，能很方便的将计算结果用图形显示出。MATLAB不但有很强的矩阵运算功能，还提供了大量的工具箱(toolbox)。 利用这些工具箱，给人们进一步开发和研究信号处理、图象处理、通讯等提供了极大的便利。在我们的课程中，主要利用信号处理工具箱(signalprocessingtoolbox)和控制系统工具箱。3.2.3窗函数法的MATLAB设计（1）汉宁窗的MATLAB语言编程delta=0.4*pi;%阻带截止频率Wc=0.2*pi;通带截止频率 AS=30;阻带最小衰减 N=ceil(8*pi/delta)+1;计算窗体长度 win=hanning(N);汉宁窗 h=fir1(N-1,Wc/pi,win);调用firl计算h（n），关于pi归一化 omega=linspace(0,pi,512);到pi分成512格mag=freqz(h,1,omega);数字滤波器频率响应函数 magdb=20*log10(abs(mag);%abs为取绝对值 plot(omega/pi,magdb);axis(0 pi-100 0);%axis（x、y的范围为最小值与最大值）grid，%网格图xlabel（频率）ylable（幅度/db）；（2）哈明窗的MATLAB语言编程delta=0.4*pi;wc=0.2*pi;as=30;n=ceil(8*pi/delta)+1;win=hamming(n);%哈明窗h=fir1(n-1,wc/pi,win);omega=linspace(0,pi,512);mag=freqz(h,1,omega);magdb=20*log10(abs(mag);plot(omega,magdb); axis(0 pi -100 0);grid;xlabel(归一化频率) ylabel(幅度/db);hammingdelta=0.4*pi; wc=0.2*pi; as=30; n=ceil(8*pi/delta)+1; win=hamming(n);%h=fir1(n-1,wc/pi,win); o mega=linspace(0,pi,512); mag=freqz(h,1,omega); magdb=20*log10(abs(mag); plot(omega,magdb); axis(0 pi -100 0); grid; xlabel(归一化频率) ylabel(幅度/db);3.2.4幅频特性曲线（1）汉宁窗特性曲线 图3.2.2（1）（2）哈明窗特性曲线 图3.2.2（2）3.3哈明窗与汉宁窗优缺点从图3.2.2中我们看到，对于同一种窗N的值越大主瓣宽度越小，相应地滤波器的过渡带就越窄，主瓣高度与第一旁瓣高度的差值方面基本不变，而旁瓣的个数越来越多，但旁瓣的宽度随着N的加大而减小。设计的FIR数字滤波器的阶数N不同，所得到的滤波器的过渡带会不一样，N越小，过渡带就越宽。哈明窗的FIR滤波器使能量主要集中在主瓣内，旁瓣峰值小于主瓣峰值，是各种窗函数中性能最好的一种，因此被广泛应用。采用不同的窗函数都会造成频谱泄露问题，采用哈明窗可以将能量集中在主瓣内，有效地消除频谱泄漏，并在一定程度上减少栅栏效应。汉宁（Hanning）窗又称升余弦窗，汉宁窗可以看作是3个矩形时间窗的频谱之和，或者说是3个sinc（t）型函数之和，而括号中的两项相对于第一个谱窗向左、右各移动了/T，从而使旁瓣互相抵消，消去高频干扰和漏能，是能量更加集中在主瓣中。适用于非周期性的连续信号。 哈明窗本质上和汉宁窗是一样的，只是系数不同。哈明窗比汉宁窗消除旁瓣的效果好一些而且主瓣稍窄，但是旁瓣衰减较慢是不利的方面。适当地改变系数，可得到不同特性的窗函数。其主瓣宽度与汉宁窗相同。可见哈明窗是一种高效窗函数。哈明窗和汉宁窗的主瓣具有最小的旁瓣和较大的衰减速度，是较为常用的窗函数。窗函数法的主要缺点是：一、不容易设计预先给定截止频率的滤波器；二、满足同样设计指标的情况下所设计出的滤波器的阶数通常偏大。 总结 总结课程设计是培养学生综合运用所学知识发现、提出、分析和解决实际问题，是锻炼实践能力的重要环节，是对学生实际工作能力的具体训练和考察过程。通过这次课程设计使我懂得了理论与实际相结合是很重要的，只有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理中得出结论，才能真正为社会服务，从而提高自己的实际动手能力和独立思考的能力。本次课程设计使我们巩固有关数字信号处理相关内容的基础上，提高了我们的实践能力。经过这次 的学习和对MATLAB在数字信号系统和其他方面的实训，我充分的了解了MATLAB的实用性和便捷性，对我以后在其他方面的学习大有好处。我觉得学习首先要把基础掌握，然后是绘图，各种画图函数，函数图形的加工，最后是编程，这是最重要的，能解决很多问题，也是后面学习的基础。这些方面掌握后，就可以做实践锻炼了，有很多技巧在实践中可以不断地积累。这次课程设计对我的帮助很多，同时加深了对课堂上学习的知识的理解。通过这次课程设计使我充分理解到数字信号处理的重要性和实用性，特别加深了对及其操作各方面的了解和设计，对实际操作设计中所涉及的各个方面要注意问题都有所了解。通过这次对FIR数字滤波器的设计，不仅让我将所学的知识应用到实际中，而且对知识也是一种巩固和提升充实。在老师和同学的帮助下，及时的按要求完成了设计任务，通过这次课程设计，使我获得了很多重要的知识，同时也提高了自己的实际动手和知识的灵活运用能力。通过这一周的学习，学到了不少的知识！不仅纠正了课程学习过程中出现的许多错误，还在运行MATLAB时验证了自己的一些猜想。在学习的过程中有失败，当然也有困惑，有成功，当然就有喜悦。虽然只是课程设计，但我拿出了自己的全部精力去对待，能学到知识固然值得骄傲，能认识到自己的过错和不足不也是一件幸事吗！做学问也是做人，在做学问的过程中体味做人的道理不也是一种收获吗？记得古语中说：“学，然后知不足”！希我必将更加力的学习它完善自己。想这就是我这次做课程设计最大感受吧。 总之通过本次课程设计，锻炼了自己的动手能力，使自己将课本上学到的理论知识联系到具体的实践过程中去，理论联系实践，巩固了自己所学的知识，同时也暴漏出自己的不少题。通过过本次课程设计对自己有了很大的进步，希望以后自己有更多的机会参与这样的课程践在实践中提升自己的能力与专业素质。参考文献1.高西全数字信号处理西安电子科技大学出版社20082.丁玉美数字信号处理西安电子科技大学出版社2008 3.张亚妮基于MATLAB的数字滤波器设计辽宁工程技术大学出版社20054数字信号处理实验指导书电气与电子实验中心吉林建筑工程学院20115.刘波MATLAB信号处理电子工业出版社20066.刘卫国，陈昭民，张颖Matlab程序设计与应用高等教育出版社2002 16