基于matlab的数字滤波的设计与仿真论文.doc

《基于matlab的数字滤波的设计与仿真论文.doc》由会员分享，可在线阅读，更多相关《基于matlab的数字滤波的设计与仿真论文.doc（40页珍藏版）》请在三一文库上搜索。

1、目 录摘要第一章、绪 论1、1数字滤波器的研究背景与意义1、2数字滤波器的主要应用与发展趋势1、3数字滤波器的实现方法分析第二章、数字滤波器的概述2、1数字滤波器的基本结构IIR滤波器的基本结构FIR滤波器的基本结构2、2数字滤波器的基本原理滤波器的主要性能指标IIR数字滤波器的设计方法FIR数字滤波器的设计方法IIR滤波器与FIR滤波器的分析比较第三章、典型数字滤波器及仿真3、1由模拟滤波器设计IIR数字滤波器巴特奥兹滤波器切比雪夫滤波器椭圆滤波器3、2用matlab设计数字滤波器FDATool界面用Fdatool进行带通滤波器设计3、3将系统函数由直接型化成级联型二阶节系数的确定系数转换成

2、二进制码第四章、总结4、1滤波器功能和性能总结4、2研究心得和体会参考文献摘 要本文分析了国内外数字滤波技术的应用现状与发展趋势，介绍了数字滤波器的基本结构，在分别讨论了IIR与FIR数字滤波器的设计方法的基础上，指出了传统的数字滤波器设计方法过程的复杂、计算量大、滤波特性调整困难的不足，提出了一种基于matlab软件的数字滤波器设计方法，避免了这种的不足。文中深入分析了该滤波器系统设计的功能特点、实现原理以及主要参数对滤波性能影响，阐述了使用MATLAB进行带通滤波器设计及仿真的具体方法。最后用Matlab仿真来直观的说明该设计方法的可靠性好，效率高，极大的减轻了工作量，有利于滤波器设计的最

3、优化。关键词： 带通滤波器；IIR；FIR; Matlab仿真与设计;第一章、 绪 论1、1数字滤波器的研究背景与意义当今，数字信号处理1 (DSP：Digtal Signal Processing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号

4、、机械振动信号、遥感遥测信号，等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。因此，数字信号实际上是用数字序列表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号经采样和量化后，得到的数字信号是一个二维离散空间序列。数字信号处理，就是用数值计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字信号经行滤波以限制他的频带或滤除噪音和干扰，或

5、将他们与其他信号进行分离；对信号进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等等。数字滤波技术是数字信号分析、处理技术的重要分支2-3。无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传输是至关重要的。在所有的电子系统中，使用最多技术最复杂的要算数字滤波器了。数字滤波器的优劣直接决定产品的优劣。从事电子通信业而不能熟练操作使用matlab电子线路路设计软件，在工作和学习中将是寸步难行的。在数学、电子、金融等行业，使用matlab等计算机软件对产品

6、进行设计、仿真在很早以前就已经成为了一种趋势，这类软件的问世也极大地提高了设计人员在通信、电子等行业的产品设计质量与效率。众所周知，实际过程中信号传输都要经过调制与解调这一过程，由于消息传过来的原始信号即调制信号具有频谱较低的频谱分量，这种信号在许多信道中不宜传输。因而，在通信系统的发送端通常需要有调制过程，反之在接收端则需要有解调过程。1、2数字滤波器的主要应用与发展趋势在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统

7、中，数字滤波器应用极为广泛，这里只列举部分应用最成功的领域。(1) 语音处理 语音处理是最早应用数字滤波器的领域之一，也是最早推动数字信号处理理论发展的领域之一。该领域主要包括5个方面的内容：第一，语音信号分析。即对语音信号的波形特征、统计特性、模型参数等进行分析计算；第二，语音合成。即利用专用数字硬件或在通用计算机上运行软件来产生语音；第三，语音识别。即用专用硬件或计算机识别人讲的话，或者识别说话的人；第四，语音增强。即从噪音或干扰中提取被掩盖的语音信号。第五，语音编码。主要用于语音数据压缩，目前已经建立了一系列语音编码的国际标准，大量用于通信和音频处理。近年来，这5个方面都取得了不少研究成

8、果，并且，在市场上已出现了一些相关的软件和硬件产品，例如，盲人阅读机、哑人语音合成器、口授打印机、语音应答机，各种会说话的仪器和玩具，以及通信和视听产品大量使用的音频压缩编码技术。(2) 图像处理数字滤波技术以成功地应用于静止图像和活动图像的恢复和增强、数据压缩、去噪音和干扰、图像识别以及层析X射线摄影，还成功地应用于雷达、声纳、超声波和红外信号的可见图像成像。(3) 通信在现代通信技术领域内，几乎没有一个分支不受到数字滤波技术的影响。信源编码、信道编码、调制、多路复用、数据压缩以及自适应信道均衡等，都广泛地采用数字滤波器，特别是在数字通信、网络通信、图像通信、多媒体通信等应用中，离开了数字滤

9、波器,几乎是寸步难行。其中，被认为是通信技术未来发展方向的软件无线电技术，更是以数字滤波技术为基础。(4) 电视数字电视取代模拟电视已是必然趋势。高清晰度电视的普及指日可待，与之配套的视频光盘技术已形成具有巨大市场的产业；可视电话和会议电视产品不断更新换代。视频压缩和音频压缩技术所取得的成就和标准化工作，促成了电视领域产业的蓬勃发展，而数字滤波器及其相关技术是视频压缩和音频压缩技术的重要基础。(5) 雷达雷达信号占有的频带非常宽，数据传输速率也非常高，因而压缩数据量和降低数据传输速率是雷达信号数字处理面临的首要问题。告诉数字器件的出现促进了雷达信号处理技术的进步。在现代雷达系统中，数字信号处理

10、部分是不可缺少的，因为从信号的产生、滤波、加工到目标参数的估计和目标成像显示都离不开数字滤波技术。雷达信号的数字滤波器是当今十分活跃的研究领域之一。(6) 声纳声纳信号处理分为两大类，即有源声纳信号处理和无源声纳信号处理，有源声纳系统涉及的许多理论和技术与雷达系统相同。例如，他们都要产生和发射脉冲式探测信号，他们的信号处理任务都主要是对微弱的目标回波进行检测和分析，从而达到对目标进行探测、定位、跟踪、导航、成像显示等目的，他们要应用到的主要信号处理技术包括滤波、门限比较、谱估计等。 (7) 生物医学信号处理 数字滤波器在医学中的应用日益广泛，如对脑电图和心电图的分析、层析X射线摄影的计算机辅助

11、分析、胎儿心音的自适应检测等。(8) 音乐数字滤波器为音乐领域开辟了一个新局面，在对音乐信号进行编辑、合成、以及在音乐中加入交混回响、合声等特殊效果特殊方面，数字滤波技术都显示出了强大的威力。数字滤波器还可用于作曲、录音和播放，或对旧录音带的音质进行恢复等。(9) 其他领域5数字滤波器的应用领域如此广泛，以至于想完全列举他们是根本不可能的，除了以上几个领域外，还有很多其他的应用领域。例如，在军事上被大量应用于导航、制导、电子对抗、战场侦察；在电力系统中被应用于能源分布规划和自动检测；在环境保护中被应用于对空气污染和噪声干扰的自动监测，在经济领域中被应用于股票市场预测和经济效益分析，等等。1、3

12、数字滤波器的实现方法分析数字滤波器的实现6，大体上有如下几种方法：(1) 在通用的微型机上用软件来实现。软件可以由使用者自己编写或使用现成的。自IEEE DSP Comm.于1979年推出第一个信号处理软件包以来，国外的研究机构、公司也陆续推出不同语言不同用途的信号处理软件包。这种实现方法速度较慢，多用于教学与科研。(2) 用单片机来实现。目前单片机的发展速度很快，功能也很强依靠单片机的硬件环境和信号处理软件可用于工程实际，如数字控制、医疗仪器等。(3) 利用专门用于信号处理的DSP片来实现。DSP芯片较之单片机有着更为突出的优点，如内部带有乘法器、累加器，采用流水线工作方式及并行结构，多总线

13、，速度快，配有适于信号处理的指令等，DSP芯片的问世及飞速发展，为信号处理技术应用于工程实际提供了可能。第二章、数字滤波器的概述数字滤波器可以用差分方程、单位取样响应以及系统函数等表示。对于研究系统的实现方法，即它的运算结构来说，用框图表示最为直接。一个给定的输入输出关系，可以用多种不同的数字网络来实现。在不考虑量化影响时，这些不同的实现方法是等效的；但在考虑量化影响时，这些不同的实现方法性能上就有差异。因此，运算结构是很重要的，同一系统函数H(z)，运算结构的不同，将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。IIR(无限冲激响应)滤波器与FIR(有限冲激响应)滤波器在

14、结构上有自己不同的特点，在设计时需综合考虑。2、1数字滤波器的基本结构作为线形时不变系统的数字滤波器可以用系统函数来表示，而实现一个系统函数表达式所表示的系统可以用两种方法：一种方法是采用计算机软件实现；另一种方法是用加法器、乘法器、和延迟器等元件设计出专用的数字硬件系统，即硬件实现。不论软件实现还是硬件实现，在滤波器设计过程中，由同一系统函数可以构成很多不同的运算结构。对于无限精度的系数和变量，不同结构可能是等效的，与其输入和输出特性无关；但是在系数和变量精度是有限的情况下，不同运算结构的性能就有很大的差异。因此，有必要对离散时间系统的结构有一基本认识。2、1、1 IIR滤波器的基本结构一个

15、数字滤波器可以用系统函数表示为： (2-1) 由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为： (2-2) 可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的，其差分方程如(2-2)式所示，是递归式的，即结构上存在着输出信号到输入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0z)有极点存在。前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有

16、一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和延迟单元的结构。IIR滤波器实现的基本结构有：(1)IIR滤波器的直接型结构；优点：延迟线减少一半，变为N 个，可节省寄存器或存储单元；缺点：其它缺点同直接I型。 通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统(一、二阶)来实现。(2)IIR滤波器的级联型结构；特点： 系统实现简单，只需一个二阶节系统通过改变输入系数即可完成； 极点位置可单独调整； 运算速度快(可并行进行)； 各二阶网络

17、的误差互不影响，总的误差小，对字长要求低。 缺点：不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。(3)IIR滤波器的并联型结构。优点： 简化实现，用一个二阶节，通过变换系数就可实现整个系统； 极、零点可单独控制、调整，调整1i、2i只单独调整了第i对零点，调整1i、2i则单独调整了第i对极点； 各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差； 可流水线操作。 缺点： 二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。 a、直接型 b、并联型c、串联型图2-1、IIR滤波器的基本结构2、1、2 FIR滤波器的基本结构FIR滤

18、波器7的单位抽样响应为有限长度，一般采用非递归形式实现。通常的FIR数字滤波器有横截性和级联型两种。FIR滤波器实现的基本结构有：(1)FIR滤波器的横截型结构表示系统输入输出关系的差分方程可写作： (2-3) 直接由差分方程得出的实现结构如图2-2所示：图2-2、 横截型(直接型卷积型)若h(n)呈现对称特性，即此FIR滤波器具有线性相位，则可以简化加横截型结构，下面分情况讨论：图2-3、N为奇数时线形相位 图2-4、N为偶数时线形相 FIR滤波器实现结构 FIR滤波器实现结构(2)FIR滤波器的级联型结构将H(z)分解成实系数二阶因子的乘积形式： (2-4)这时FIR滤波器可用二阶节的级联

19、结构来实现，每个二阶节用横截型结构实现。如图所示：图2-5、 FIR滤波器的级联结构这种结构的每一节控制一对零点，因而在需要控制传输零点时可以采用这种结构。2、2数字滤波器的设计原理数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。IIR滤波器的特征是，具有无限持续时间冲激响应。这种滤波器一般需要用递归模型来实现，因而有时也称之为递归滤波器。FIR滤波器的冲激响应只能延续一定时间，在工程实际中可以采用递归的方式实现，也可以采用非递归的方式实现。数字滤波器的设计方法有多种，如双线性变换法、窗函数设计法、插值逼近法和Chebysh

20、ev逼近法等等。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善，不仅数字滤波器的计算机辅助设计有了可能，而且还可以使设计达到最优化。数字滤波器设计的基本步骤如下：(1)确定指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于FIR滤波器的设计。第二种指标是相对指标。它以分贝值的形式给出要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中具有线性

21、相位。运用线性相位响应指标进行滤波器设计具有如下优点：只包含实数算法，不涉及复数运算；不存在延迟失真，只有固定数量的延迟；长度为N的滤波器(阶数为N-1)，计算量为N/2数量级。因此，本文中滤波器的设计就以线性相位FIR滤波器的设计为例。(2)逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。(3)性能分析和计算机仿真上两步的结果是得到以差分或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器

22、，再分析滤波结果来判断。2、2、1滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。以低通滤波器特性为例，频率响应有通带、过渡带及阻带三个范围。在通带内： 1- AP| H(ej)| 1 |c在阻带中， |H(ej)| Ast st | 图2-6、低通滤波器频率响应幅度特性的容限图其中c 为通带截止频率, st为阻带截止频率，Ap为通带误差, Ast为阻带误差。与模拟滤波器类似，数字滤波器按频率特性划分为低通、高通、带通、带阻、全通等类型，由于数字滤波器的频率响应是周期性的，周期为2。各种理想数字滤波器的幅度频率响应

23、如图所示：图2-7、 各种理想数字滤波器的幅度频率响应2、2、2 IIR数字滤波器的设计方法目前，IIR数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便，IIR数字滤波器的设计步骤是：(1)按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标；(2)根据转换后的技术指标设计模拟低通滤波器H(s)；(3)在按一定规则将H(s)转换为H(z)。若所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、带通或

24、者带阻滤波器，那么还有步骤：(4)将高通、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再由冲击响应不变法或双线性变换将H(s)转换为所需的H(z)。s - z 映射的方法有：冲激响应不变法、阶跃响应不变法、双线性变换法等。下面讨论双线性变换法。双线性变换法8是指首先把s 平面压缩变换到某一中介平面s1 的一条横带(宽度为2T,即从- T到T) ,然后再利用的关系把s1平面上的这条横带变换到整个z 平面。这样s 平面与z 平面是一一对应关系, 消除了多值变换性, 也就消除了频谱混叠现象。s 平面到z 平面的变换可采用 (2-5)

25、 (2-6)令 ，有： (2-7)从s1 平面到z 平面的变换,即 (2-8)代入上式，得到： (2-9) 一般来说，为使模拟滤波器的某一频率与数字滤波器的任一频率有一定的对应关系，可引入代定常数c， (2-10)则 (2-11) 这种s 平面与z 平面间的单值映射关系就是双线性变换。有了双线性变换，模拟滤波器的数字化只须用进行置换。2、2、3 FIR数字滤波器的设计方法IIR滤波器7的优点是可利用模拟滤波器设计的结果，缺点是相位是非线性的，若需要线性相位，则要用全通网络进行校正。FIR滤波器的优点是可方便地实现线性相位。 FIR滤波器单位冲激响应h(n)的特点：其单位冲激响应h(n)是有限长

26、(),系统函数为： (2-12)在有限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处，则滤波器具有准确的线性相位。窗函数设计法： 一般是先给定所要求的理想滤波器频率响应，由导出，我们知道理想滤波器的冲击响应是无限长的非因果序列，而我们要设计的是是有限长的FIR滤波器，所以要用有限长序列来逼近无限长序列，设：(2-13) 常用的方法是用有限长度的窗函数w(n)来截取 即：(2-14) 这里

27、窗函数就是矩形序列RN(n),加窗以后对理想低通滤波器的频率响应将产生什么样的影响呢?根据在时域是相乘关系,在频域则是卷积关系： (2-15) 其中,为矩形窗谱, 是FIR滤波器频率响应. 通过频域卷积过程看的幅度函数H()的起伏现象,可知,加窗处理后，对理想矩形的频率响应产生以下几点影响：(1)使理想频率特性不连续点处边沿加宽，形成一个过渡带，其宽度等于窗的频率响应的主瓣宽度。(2)在截止频率的两边的地方即过渡带的两边，出现最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，则取决于旁瓣的多少。(3)改变N，只能改变窗谱的主瓣宽度，改变的坐标比例以及改变的绝对

28、值大小，但不能改变主瓣与旁瓣的相对比例(此比例由窗函数的形状决定)。(4)对窗函数的要求a、窗谱主瓣尽可能窄，以获取较陡的过渡带；b、尽量减小窗谱的最大旁瓣的相对幅度；即能量集中于主瓣，使肩峰和波纹减小，增大阻带的衰减。频率采样法：窗函数设计法是从时域出发，把理想的hd(n)用一定形状的窗函数截取成有限长的h(n)，来近似理想的hd(n)，这样得到的频率响应逼近于所要求的理想的频率响应。频率抽样法则是从频域出发，把给定的理想频率响应 加以等间隔抽样得到 ，然后以此 作为实际FIR滤波器的频率特性的抽样值H(k)，即(2-16) 知道H(k)后，由DFT定义可唯一确定有限长序列 h(n)，利用这

29、N个频域抽样值H(k)同样利用频率内插公式可得FIR滤波器的系统函数H(z)，及频率响应 ，即： 频率抽样法内插公式： (2-17) 频率抽样法小结优点：可以在频域直接设计，并且适合于最优化设计。缺点：抽样频率只能等于 2/N 的整数倍，或等于2/N 的整数倍加上/N。因而不能确保截止频率 的自由取值，要想实现自由地选择截止频率，必须增加抽样点数N，但这又使计算量增大。 为了提高逼近质量，减少通带边缘由于抽样点的陡然变化而引起的起伏振荡。有目的地在理想频率响应的不连续点的边缘，加上一些过渡的抽样点，增加过渡带，减少起伏振荡。2、3 IIR滤波器与FIR滤波器的分析比较 (1) 选择数字滤波器是

30、必须考虑经济问题，通常将硬件的复杂性、芯片的面积或计算速度等作为衡量经济问题的因素。在相同的技术指标要求下，由于IIR数字滤波器存在输出对输入的反馈，因此可以用较少的阶数来满足要求，所用的存储单元少，运算次数少，较为经济。例如，用频率抽样法设计一个阻带衰减为20dB的FIR数字滤波器，要33阶才能达到要求，而用双线性变换法只需45阶的切比雪夫IIR滤波器就可达到同样的技术指标。这就是说FIR滤波器的阶数要高510倍左右。(2)很多情况下,FIR数字滤波器的线性相位与它的高阶数带来的额外成本相比是非常值得的。对于IIR滤波器，选择性越好，其相位的非线性越严重。如果要使IIR滤波器获得线性相位，又

31、满足幅度滤波器的技术要求，必须加全通网络进行相位校正，这同样将大大增加滤波器的阶数。 (3) FIR滤波器主要采用非递归结构，因而无论是理论上还是实际的有限精度运算中他都是稳定的，有限精度运算误差也较小。IIR滤波器必须采用递归结构，极点必须在z平面单位圆内才能稳定。 (4) 对于FIR滤波器，由于冲激响应是有限长的，因此可以用快速傅里叶变换算法，这样运算速度可以快得多。IIR滤波器不能进行这样的运算。(5) 从设计上看，IIR滤波器可以利用模拟滤波器设计的现成的闭合公式、数据和表格，可以用完整的设计公式来设计各种选频滤波器。一旦选定了已知的一种逼近方法(如巴特奥兹，切比雪夫等)，就可以直接把

32、技术指标带入一组设计方程计算出滤波器的阶次和系统函数的系数(或极点和零点)。FIR滤波器则一般没有现成的设计公式。窗函数法只给出了窗函数的计算公式，但计算通带和阻带衰减仍无显式表达式。(6) IIR滤波器主要是设计规格化、频率特性为分段常数的标准低通、高通、带通和带阻滤波器。FIR滤波器则灵活很多，适应性很广。表2-1、 两种滤波器特点比较分析FIR滤波器IIR滤波器运算误差一般无反馈，运算误差小，噪声功能小，运算速度快有反馈，由于运算中的四舍五入会产生极限环稳定性极点全部在原点(永远稳定)无稳定性问题有稳定性问题设计方法一般无解析的设计公式，要借助计算机程序完成利用AF的成果，可简单、有效地

33、完成设计，要借助成熟模拟滤波器设计成果来设计结构非递归系统递归系统阶数高低设计结果可得到幅频特性(可以多带)和线性相位(最大优点)只能得到幅频特性，相频特性未知，如需要线性相位，须用全通网络校准，但增加滤波器阶数和复杂性性能可得到严格的线性相位，可是成本高，信号延时较大可用较低阶数获得较高选择性，所用存储单元少，计算量小，经济高效第三章、典型数字滤波器及仿真3、1由模拟滤波器设计IIR数字滤波器在之前的部分中已经说明，理想的滤波器是非因果的，即物理上不可实现的系统。工程上常用的模拟滤波器都不是理想的滤波器。但按一定规则构成的实际滤波器的幅频特性可逼近理想滤波器的幅频特性，例如巴特奥兹(Butt

34、erworth)、切比雪夫(Chebyshev)滤波器和椭圆滤波器等。3、1、1巴特奥兹滤波器巴特沃兹滤波器9(Butterworth 滤波器)特点：具有通带内最大平坦的振幅特性，且随f，幅频特性单调。其幅度平方函数：(3-1)N为滤波器阶数，如图3-1图3-1、 巴特沃斯滤波器振幅平方特性通带： 使信号通过的频带阻带：抑制噪声通过的频带过渡带：通带到阻带间过渡的频率范围 c ：截止频率。 过渡带为零理想滤波器 阻带|H(j )|=0 通带内幅度|H(j)|=cons. H(j)的相位是线性的图3-1中，N增加，通带和阻带的近似性越好，过渡带越陡。通带内，分母/c1, ( /c)2N1， (

35、/c)2N1, 增加， A(2)快速减小。 =c, ,幅度衰减，相当于3db衰减点。 振幅平方函数的极点 (3-2) 可见，Butter worth滤波器 的振幅平方函数有2N个极点，它们均匀地分布在|S|=c的圆周上。考虑到系统的稳定性，知DF的系统函数是由S平面左半部分的极点(SP3，SP4，SP5)组成的，它们分别为：(3-3)系统函数为： (3-4)令 ，得归一化的三阶BF： (3-5)如果要还原的话，则有(3-6)3、1、2切比雪夫滤波器巴特奥兹低通滤波器的幅频特性随的增加而单调下降，当N较小时，阻带幅频特性下降较慢，要想使其幅频特性接近理想低通滤波器，就必须增加滤波器的阶数，这就将

36、导致模拟滤波器使用的原件增多，线路趋于复杂。切比雪夫滤波器10的阻带衰减特性则有所改善。特点：误差值在规定的频段上等幅变化。巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止频率 c处，幅度下降很多，或者说，为了使通常内的衰减足够小，需要的阶次(N)很高，为了克服这一缺点，采用切比雪夫多项式逼近所希望的 。 切比雪夫滤波器的 在通带范围内是等幅起伏的，所以同样的通带衰减，其阶数较巴特沃兹滤波器要小。可根据需要对通带内允许的衰减量(波动范围)提出要求，如要求波动范围小于1db。 振幅平方函数为 (3-7)式中 有效通带截止频率与通带波纹有关的参量,大,波纹大,0 1。Vn(x

37、)N阶切比雪夫多项式，定义为(3-8)(3-9)如图3-1,通带内 ,，变化范围1- c ，随/c ，0(迅速趋于零)当 =0时，(3-10)N为偶数，,(min) ， (3-11)N为奇数，,(max)， (3-12)图3-2、切比雪夫滤波器的振幅平方特性有关参数的确定：a. 通带截止频率 ，预先给定b. 由通带波纹表为 (3-13)(3-14)给定通带波纹值分贝数 后，可求。(3-15)c. 阶数N由阻带的边界条件确定。(,A事先给定)(3-16)(3-17)(3-18) (3-19)得 (3-20)3、1、3椭圆滤波器特点：幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和给定的波纹要

38、求，椭圆滤波器能获得较其它滤波器为窄的过渡带宽，就这点而言，椭圆滤波器11是最优的，其振幅平方函数为 (3-21) 式中，RN(，L)为雅可比椭圆函数，L是一个表示波纹性质的参量。 图3-3、 N=5时 的特性曲线 由图可见，在归一化通带内(-11)，在(0,1)间振荡，而超过L后，在L2， 间振荡。L越大，L也变大。这一特点使滤波器同时在通带和阻带具有任意衰减量。 下图为典型的椭园滤波器振幅平方函数：图3-4、 椭圆滤波器的振幅平方函数 图中和A的定义与切比雪夫滤波器相同。 当c、s、和A确定后，阶次N的确定方法为 ： 确定参数 确定参量 N=式中K(k)= 为第一类完全椭圆积分。3、2用m

39、atlab设计数字滤波器传统的数字滤波器的设计过程复杂，计算工作量大，滤波特性调整困难。利用MATLAB信号处理工具箱(Signal Processing Toolbox)12-14可以快速有效的设计由软件组成的常规数字滤波器的设计方法。利用MATLAB设计滤波器，可以随时对比设计要求和滤波器特性调整参数，直观简便，极大的减轻了工作量，有利于滤波器设计的最优化。3、2、1 FDATool界面FDATool(Filter Design & Analysis Tool)是MATLAB信号处理工具箱里专用的滤波器设计分析工具，MATLAB6.0以上的版本还专门增加了滤波器设计工具箱(Filter D

40、esign Toolbox)。FDATool可以设计几乎所有的基本的常规滤波器，包括FIR和IIR的各种设计方法。FDATool界面总共分两大部分，一部分是Design Filter，在界面的下半部，用来设置滤波器的设计参数，另一部分则是特性区，在界面的上半部分，用来显示滤波器的各种特性。Design Filter部分主要分为：Filter Type(滤波器类型)选项，包括Lowpass(低通)、Highpass(高通)、Bandpass(带通)、Bandstop(带阻)和特殊的FIR滤波器。Design Method(设计方法)选项，包括IIR滤波器的Butterworth(巴特沃思)法、C

41、hebyshev Type I(切比雪夫I型)法、 Chebyshev Type II(切比雪夫II型) 法、Elliptic(椭圆滤波器)法和FIR滤波器的Equiripple法、Least-Squares(最小乘方)法、Window(窗函数)法。Filter Order(滤波器阶数)选项，定义滤波器的阶数，包括Specify Order(指定阶数)和Minimum Order(最小阶数)。在Specify Order中填入所要设计的滤波器的阶数(N阶滤波器，Specify OrderN-1)，如果选择Minimum Order则MATLAB根据所选择的滤波器类型自动使用最小阶数。Frenq

42、uency Specifications选项，可以详细定义频带的各参数，包括采样频率Fs和频带的截止频率。它的具体选项由Filter Type选项和Design Method选项决定，例如Bandpass(带通)滤波器需要定义Fstop1(下阻带截止频率)、Fpass1(通带下限截止频率)、Fpass2(通带上限截止频率)、Fstop2(上阻带截止频率)，而Lowpass(低通)滤波器只需要定义Fstop1、Fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定的，所以只需要定义通带截止频率，而不必定义阻带参数。Magnitude Specifications选项，可以定

43、义幅值衰减的情况。例如设计带通滤波器时，可以定义Wstop1(频率Fstop1处的幅值衰减)、Wpass(通带范围内的幅值衰减)、Wstop2(频率Fstop2处的幅值衰减)。当采用窗函数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。Window Specifications选项，当选取采用窗函数设计时，该选项可定义，它包含了各种窗函数，在通带内的衰减为6dB。 图3-5、 FDATool的操作页面通过菜单选项Analysis可以在特性区看到所设计滤波器的幅频响应、相频响应、零极点配置和滤波器系数等各种特性。设计完成后将结果保存为*.fda文件。在设计过程中，可以对比滤波器幅频相

44、频特性和设计要求，随时调整参数和滤波器类型，以便得到最佳效果。其它类型的FIR滤波器和IIR滤波器也都可以使用FDATool来设计。3、2、2 用Fdatool进行带通滤波器设计给定的数字带通滤波器的参数是：通带为45Hz 55Hz，低截止频率为40Hz，高截止频率为60Hz，通带内衰减不大于3db，阻带衰减大于80db。借助Matlab仿真，可以得到设计成不同类型的滤波器所需的阶数：表3-1、采用不同类型滤波器实现所需的阶数滤波器类型最低阶数稳定性FIR滤波器Kaiser窗函数201稳定Blackman窗函数180稳定Equirpple86稳定IIR滤波器Butterworth滤波器26稳定Chebishev滤波器16稳定椭圆滤波器12稳定 可见，FIR滤波器阶数过高，导致的直接结果是消耗的资源较多，成本增加。因此采用IIR滤波器实现比较合适。而在IIR滤波器