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