毕业设计论文基于MATLAB的IIR滤波器的设计与仿真.doc

《毕业设计论文基于MATLAB的IIR滤波器的设计与仿真.doc》由会员分享，可在线阅读，更多相关《毕业设计论文基于MATLAB的IIR滤波器的设计与仿真.doc（41页珍藏版）》请在三一办公上搜索。

1、湖 北 省 高 等 教 育 自 学 考 试 本 科 毕 业 生 论 文评 审 意 见 表论文题目：基于MATLAB的IIR滤波器的设计与仿真 姓 名： 专 业： 电子信息工程 准考证号： 工作单位： 武汉纺织大学成教学院 填表日期：2011 年 2 月 15 日湖北省高等教育自学考试委员会办公室制指导老师单位职称指导教师评语： 指导老师： 年 月 日答辩小组评语： 组长签名： 年 月 日答辩委员会意见： 负责人签名： （签章） 年 月 日 基于MATLAB的IIR滤波器的设计与仿真I基于MATLAB的IIR滤波器的设计与仿真论文原创性声明我以诚信声明： 本人所呈交的论文是我个人在导师指导下进行

2、的研究工作及取得的研究成果。尽我所知，除了文中特别加以标注和致谢的地方外，论文中不包含其他人已经发表或撰写过的研究成果，也不包含为获得武汉纺织大学或其它教育机构的学位或证书而使用过的材料。与我一同工作的同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意。作者（签字）： 签字日期：201年 2 月15 日版权使用授权书武汉纺织大学有权保存学位论文的电子和纸质文档，可以借阅或上网公开本学位论文的全部或部分内容，允许论文被查阅，可以向有关部门或机构送交并授权其保存、借阅或上网公布本学位论文的全部或部分内容。对于保密论文，按保密的有关规定进行法律处理。作者（签字）： 签字日期： 2011

3、 年2 月15 日 指导老师（签字）： 签字日期： 年 月 日 摘 要在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号分析都是基于滤波器而进行的。而数字滤波器是通过数值运算实现滤波，具有处理精度高，稳定，灵活，不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应（IIR）数字滤波器和有限长冲激响应（FIR）数字滤波器。实现IIR数字滤波器的阶次较低，所用的储存单元较少，效率高，精度高，而且能够保留一些模拟滤波器的优良特性，因此应用很广。Matlab软件以矩阵运算为基础，把计算，可视化及程序设计有机融合到

4、交互式工作环境中，并且为数字滤波器的研究和应用提供了一个直观，高效，便捷的利器。尤其是Matlab中的信号处理工具箱使各个领域的研究人员可以直观方便的进行科学研究和工程应用。本文首先介绍了数字滤波器的概念，分类以及设计要求。接着利用Matlab函数语言编程，用信号处理图形界面FDATool来设计滤波器以及Sptool界面的设计方法，并用FDATool模拟IIR数字滤波器处理信号。重点设计Chebyshev I型和Chebyshev II型数字低通滤波器，并介绍最优化设计。关键词：IIR；滤波器；FDATool；Sptool；Simulink AbstractIn modern communic

5、ation systems, Because often mixed with various signal complex components, So many signal analysis is based on filters, and the digital filter is realized through numerical computation, digital filter with high precision, stability and flexibility, dont exist, can realize the impedance matching simu

6、lating the special filter cannot achieve filter function. Digital filter according to its impulse response function and characteristics of the time can be divided into two kinds, namely the infinite impulse response (IIR) digital filter and finite impulse response (FIR digital filters). The order of

7、 realizing IIR filter is used, low and high efficiency less storage unit, high precision, and can keep some simulation characteristics of filter, so it is widely used. Matlab software based on matrix computation, the calculation, visualization and program design of integration to intuitive, efficien

8、t and convenient tool. Especially in the Matlab signal processing to all areas of research toolbox personnel can easily for scientific and engineering application. This paper introduces the concept of digital filter, classification and design requirements. Then using MATLAB language programming, whi

9、t functions of signal processing FDATool graphical interface design of interface design and Sptool filter, and FDATool analog signal processing IIR digital filter. Key design Chebyshev type I and II digital chebyshev lowpass filter, and introduces optimization design. Key words: IIR; Filter; FDATool

10、; Sptool; Simulink 目 录摘 要.II前 言1第一章 数字滤波器21.1 数字滤波器的概念21.2 数字滤波器的作用31.3 数字滤波器的分类31.4 数字滤波器的设计要求51.4.1数字滤波器的设计过程51.4.2数字滤波器的技术要求5第二章 IIR数字滤波器的设计方法62.1 IIR数字滤波器的特点62.2 IIR数字滤波器的设计步骤72.3 用脉冲相应不变法设计IIR数字滤波器82.3.1设计原理82.3.2脉冲相应不变法优缺点122.3.3应用实例122.4 双线性变换法设计IIR数字滤波器152.4.1设计原理152.4.2双线性变换法优缺点172.4.3应用实例1

11、8第三章 IIR滤波器的MATLAB设计193.1 MATLAB简介193.2 IIR数字滤波器的典型设计法203.2.1 IIR数字滤波器的典型设计法的一般步骤203.2.2举例说明213.3 IIR数字滤波器的直接设计法233.3.1 IIR数字滤波器的直接设计法的特点233.3.2举例说明243.4 FDATOOL介绍和界面设计263.5 FDATOOL设计IIR数字滤波器27第四章 SIMULINK仿真IIR滤波器284.1 SIMULINK仿真284.1.1 SIMULINK仿真的功能284.1.2 SIMULINK仿真的特点284.2 SIMULINK仿真实例29第五章 设计总结3

12、1结语32参考文献33致 谢3434前 言随着信息时代和数字时代的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。目前数字信号处理在通信，语音，图像，自动控制，雷达，军事，航天，医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中起着重要作用并已获得广泛应用的是数字滤波器（DF, Digital Filter）, 根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIR (Infinite Impulse Response) 滤波器和有限冲激响应FIR （Finite Impulse Response） 滤波器。与FIR滤波器相比，IIR的实现采用的是递归结构，极点须在单

13、位圆内，在相同设计指标下，实现滤波器的阶次较低，即所用的储存单元少，从而经济效率高。MATLAB是英文MATrix LABoratory(矩阵实验室)的缩写。它是美国的MathWorks公司推出的一套用于科学计算和图像处理可视化，高性能语言与软件环境。MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常方便。 第一章 数字滤波器1.1 数字滤波器的概念凡是可以使信号

14、中特定的频率成分通过，而极大地衰减或抑制其他频率成分的装置或系统都称之为滤波器，相当于频率“筛子”。滤波器是指用来对输入信号进行滤波的硬件和软件。数字滤波器指输入输出均为数字信号，经过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器是一个离散时间系统（按预定的算法，将输入离散时间信号转换为所要求的输出离散时间信号的特定功能装置）。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即12抽样频率点呈镜像对称。为得到模拟信号，

15、数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。数字滤波器是一个离散时间系统，输入x(n)是一个时间序列，输出y(n)也是一个时间序列。如数字滤波器的系统函数为H(z)，其脉冲响应为h(n)，则在时间域内存在下列的关系： 在z域内，输入和输出存在下列关系： 式中，X(z)、Y(z)分别为输入x(n)和输出y(n)的z变换。同样在频率域内，输入和输出存在下列关系： 式中，为数字滤波器的频率特性; 和分别为x(n)和y(n)的频谱。为数字角

16、频率，单位rad。通常设计在某些频段的响应值为1，在某些频段的响应为0。和的乘积在频率响应为1的那些频段的值仍为，即在这些频段的振动可以无阻碍地通过滤波器，这些频带为通带。和的乘积在频率响应为0的那些频段的值不管大小如何均为零，即在这些频段里的振动不能通过滤波器，这些频带称为阻带。一个合适的数字滤波器系统函数H(z)可以根据需要改变输入x(n)的频率特性。经数字滤波器处理后的信号y(n)保留信号x(n)中的有用频率成分，去除无用频率成分1。1.2 数字滤波器的作用滤波器是一种对信号有处理作用的器件或电路。数字滤波器是数字信号处理技术的重要内容。和模拟滤波器一样，数字滤波器的主要功能是对数字信号

17、进行处理，最常见的处理是保留数字信号中的有用频率成分，去除信号中的无用频率成分。 主要作用是：让有用信号尽可能无衰减的通过，对无用信号尽可能大的衰减。 滤波器一般有两个端口，一个输入信号、一个输出信号 利用这个特性可以将通过滤波器的一个方波群或复合噪波，而得到一个特定频率的正弦波4。滤波器的作用：（1）将有用的信号与噪声分离，提高信号的抗干扰性及信噪比； （2）滤掉不感兴趣的频率成分，提高分析精度； （3）从复杂频率成分中分离出单一的频率分量。1.3 数字滤波器的分类数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。1从处理信号分：经

18、典滤波器、现代滤波器（1）经典滤波器：即一般滤波器（输入信号中有用的频率成分和希望滤除的频率成分各占有不同的频带，通过以合适的选频滤波器达到目的）。（2）现代滤波器：如维纳滤波器，卡尔曼滤波器，自适应滤波器等最佳滤波器（按随机信号内部的系统分部规律，从干扰中最佳提取信号）。2从数字滤波器的单元冲击响应来看，可分为俩大类：有限冲击响应（FIR）数字滤波器和无限冲击响应（IIR）数字滤波器。IIR滤波器系统函数的极点可以位于单位圆内的任何地方，因此可以用较低的阶数获得较高的选择性，所用的存储单元少，经济而效率高，但是系统函数的极点也可能位于单位圆外，可能会引起系统的不稳定。同时IIR滤波器的相位是

19、非线性的，且它的选择性越好，相位的非线性就越严重。相反FIR滤波器却可以得到严格的线性相位，然而由于FIR滤波器的系统函数的极点个定在原点，所以只能用较高的阶数来实现其高选择性，对于同样的滤波器设计指标，FIR滤波器所要求的阶数要比IIR高5至10倍，所以成本较高，信号延迟也较大。但是如果要求相同的线性，则IIR滤波器就必须加全通网络进行相位校正，同样也要增加滤波器网络的阶数和复杂性。FIR滤波器可以用非递归的方法实现，在有限精度下不会产生振荡，同时由于量化舍入以及系数的不确定性所引起的误差的影响要比IIR滤波器小的多，并且FIR滤波器可以采用FFT算法，运算速度快。但是不想IIR滤波器可以借

20、助模拟滤波器的成果，FIR滤波器没有现成的计算公式，必须要用计算机辅助设计软件（如MATLAB）来计算。由此可知，FIR滤波器应用比较广，而IIR滤波器应用比较广，而IIR滤波器应用比较广，FIR滤波器则用在相位要求不是很严格的场合。IIR滤波器： FIR滤波器：3从功能上分为如下4类：（1）低通滤波器（LPF）它允许信号中的低频或直流分量通过，抑制高频分量或干扰和噪声。常用的低通滤波器是用电感和电容组合而成的，电容并联在要滤波的信号线与信号地之间（滤除差模干扰电流）或信号线与机壳地或大地之间（滤除共模干扰电流）电感串联在要滤波的信号线上。按照电路结构分，有单电容型（C型），单电感型，L型和反

21、型，T型，型。（2）高通滤波器（HPF）它允许信号中的高频分量通过，抑制低频或直流分量。用于干扰频率比信号频率低的场合，如在一些靠近电源线的敏感信号线上滤除电源谐波造成的干扰。（3）带通滤波器（BPF）它允许一定频段的信号通过，抑制低于或高于该频段的信号、干扰和噪声。用于信号频率仅占较窄带宽的场合，如通信接收机的天线端口上要安装带通滤波器，仅允许通信信号通过。（4）带阻滤波器（BSF）它抑制一定频段内的信号，允许该频段以外的信号通过。用于干扰频率带宽较窄，而信号频率较宽的场合，如距离大功率电台很近的电缆端口处要安装带阻频率等于电台发射频率的带阻滤波器4。1.3 数字滤波器的设计要求1.4.1数

22、字滤波器的设计过程：（1）按设计任务，确定滤波器性能要求，制定技术指标（2）用一个因果稳定的离散LSI系统的系统函数H(z)逼近此性能指标（3）利用有限精度算法实现此系统函数：如运算结构、字长的选择等（4）实际技术实现：软件法、硬件法或DSP芯片法。1.4.2数字滤波器是技术要求：我们通常用的数字滤波器一般属于选频滤波器。假设数字滤波器的传输函数H(e j)用下式表示：幅频特性|H(ejw)|： 信号通过滤波器后的各频率成分衰减情况。相频特性j(w)： 各频率成分通过滤波器后在时间上的延时情况。理想滤波器不可实现，只能以实际滤波器逼近：通带内和阻带内允许的衰减一般用dB数表示，通带内允许的最大

23、衰减用ap表示，阻带内允许的最小衰减用as表示， ap和as分别定义为：滤波器设计的首要任务是要确定满足技术指标要求的系统函数，系统函数必须是物理可实现的，因而要满足因果性和稳定性。由于为复数则有由于集总系统的系统函数是s的有理函数，则是的有理函数，必为的有理函数1。 第二章 IIR数字滤波器的设计方法2.1 IIR数字滤波器的特点IIR数字滤波器的系统函数可以写成封闭函数的形式。IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不

24、断累积，有时会产生微弱的寄生振荡。 IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。 IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。 在MATLAB下设计IIR滤波器可使用Butterworth函数设计出巴特沃斯滤波器，使用Cheby1函数设计出契比雪夫I型滤波器，使用Cheby2设计出契比雪夫II型滤波器，使用ellipord函

25、数设计出椭圆滤波器。下面主要介绍前两个函数的使用。 与FIR滤波器的设计不同，IIR滤波器设计时的阶数不是由设计者指定，而是根据设计者输入的各个滤波器参数（截止频率、通带滤纹、阻带衰减等），由软件设计出满足这些参数的最低滤波器阶数。在MATLAB下设计不同类型IIR滤波器均有与之对应的函数用于阶数的选择。 IIR单位响应为无限脉冲序列FIR单位响应为有限的 IIR幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上； FIR幅频特性精度较之于IIR低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变。这是很好的性质。 另外有限的单位响应也有利于对数字信

26、号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要3。2.2 IIR数字滤波器的设计步骤IIR滤波器的设计就是根据滤波器的性能指标要求，设计滤波器的分子和分母多项式系数。它和FIR滤波器相比优点是在满足相同性能指标要求的前提下，IIR滤波器的阶数明显要低于FIR滤波器，但IIR滤波器的相位是非线性的。在数字滤波器设计中，首先让我们想到的是能否利用模拟滤波器的设计成果进行数字滤波器设计，答案是肯定的。现在我们所讲的IIR经典设计就是将已设计好的模拟滤波器按一定变换原理转换为数字滤波器。该方法先根据滤波器的技术指标设计出相应的模拟滤波器，然后再将设计好的模拟滤波器变换成数字滤波器。在

27、MATLAB中，经典法设计IIR数字滤波器主要采用以下步骤： 模拟离散化/bilinear,impinvar模拟低通滤波器原型/buttap,cheb1ap,cheb2ap,ellipap,besselap频率变换/lp2lp,lp2hp,lp2bp,lp2bsIIR数字滤波器由图可见，IIR滤波器设计利用了模拟滤波器的设计成果。第一步和第二步在前面已详细讨论过。第二步完成后，一个达到期望性能指标的模拟滤波器（低通、高通、带通、带阻等）已经设计出来。第三步模拟离散化的主要任务就是把模拟滤波器变换成数字滤波器，即把模拟滤波器的系统函数H(s)映射为数字滤波器的系统函数H(z)。这样，数字滤波器的

28、设计工作全部完成3。2.3 IIR数字滤波器的设计步骤脉冲响应不变法是一种时域上的转换方法，这种方法适用于系统函数可以用部分分式分解成单阶极点和滤波器是一个带限系统的情况，它使数字滤波器的冲击响应h(n）等于模拟滤波器的单位冲击响应ha(t)的采样值数字滤波器的脉冲响应与模拟滤波器的脉冲响应相似。以时问问隔T对其进行采样，得到h(n)， 经过z=esT 式的映射，可知：s平面的左半平面映射为z平面的单位圆内， 因此一个因果的和稳定的模拟滤波器映射成因果的和稳定的数字滤波器。脉冲响应不变法是使数字滤波器在时域上模仿模拟滤波器，由于s平面到z平面的多值映射关系，会导致数字滤波器频率响应出现混叠现象

29、3。2.3.1设计原理脉冲响应不变法的基本思路是直接设计一个数字滤波器并让它的时间特性逼近一个模拟滤波器。为了达到时间特性的最佳逼近，它把模拟滤波器冲激响应的均匀取值，作为数字滤波器的单位脉冲响应，即 然后将h(n)通过z变换即求得H（z）。脉冲响应不变法实即时域取样法，整个过程是先根据给定的指标设计一个模拟滤波器，进而按下列变换的顺序，最后求得数字滤波器的系统函数H（z），即： 例如已知一模拟滤波器其系统函数为： 则有 ， 故得 比较H1（z）与H1（s）可见，从模拟的系统函数变换为数字的系统函数，其极点从s=p1变换为z=ep1T，它们之间存在着： （8.3）的映射关系。对任一物理可实现的

30、模拟滤波器来说，由于系统函数是s的有理函数，所以从H(s)变换到H（z）存在有：即 式子所表示的映射关系，正是拉氏变换与z变换之间的映射关系。这说明脉冲响应不变法完全满足上节所述的设计IIR滤波器的两个基本条件。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n) 模仿模拟滤波器的冲击响应ha(t), 使h(n)正好等于ha(t)的采样值，即h(n)=ha(nT)T为采样周期。如以Ha(s)及H（z）分别表示ha(t)的拉氏变换及h(n)的z变换，即Ha(s)=Lha(t)H(z)=Zh(n) 则根据采样序列z变换与模拟信号拉氏变换的关系，得: 上式表明，采用脉冲响应

31、不变法将模拟滤波器变换为数字滤波器时，它所完成的S平面到Z平面的变换，正是以前讨论的拉氏变换到Z变换的标准变换关系，即首先对Ha(s)作周期延拓，然后再经过z=est的映射关系映射到Z平面上。 z=est的映射关系表明，S平面上每一条宽为2/T的横带部分，都将重叠地映射到Z平面的整个全部平面上。每一横带的左半部分映射到Z平面单位圆以内，每一横带的右半部分映射到Z平面单位圆以外，j轴映射在单位圆上，但j轴上的每一段2/T都对应于绕单位圆一周，如下图所示：图-1 脉冲响应不变法的映射关系应当指出，Z=est的映射关系反映的是Ha(s)的周期延拓与H（z）的关系，而不是Ha(s)本身与H（z）的关系

32、，因此，使用脉冲响应不变法时，从Ha(s)到H(z)并没有一个由S平面到Z平面的简单代数映射关系，即没有一个s=f(z)的代数关系式。另外，数字滤波器的频响也不是简单的重现模拟滤波器的频响，而是模拟滤波器频响的周期延拓，周期为S=2/T=2fs，即正如第一章采样定理中所讨论的，如果模拟滤波器的频响带限于折叠频率S/2以内，即 Ha(j)=0 |/T 这时数字滤波器的频响才能不失真地重现模拟滤波器的频响（在折叠频率以内） | 但任何一个实际的模拟滤波器，其频响都不可能是真正带限的，因此不可避免地存在频谱的交叠，即混淆，这时，数字滤波器的频响将不同于原模拟滤波器的频响而带有一定的失真。模拟滤波器频

33、响在折叠频率以上衰减越大，失真则越小，这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。 脉冲响应不变法特别适用于用部分分式表达的传递函数，模拟滤波器的传递函数若只有单阶极点，且分母的阶数高于分子阶数NM，则可表达为部分分式形式： 其拉氏反变换为： 其中u(t)为单位阶跃函数。对ha(t)采样就得到数字滤波器的单位脉冲响应序列再对h(n)取Z变换，得到数字滤波器的传递函数：第二个求和为等比级数之和， ，要收敛的话，必有 ，所以有 比较部分分式形式的Ha(s)和上式H(z)可以看到，S平面上的极点s=si,变换到Z平面上是极点 ,而Ha(s)与H(z)中部分分式所对应的系数不变。如果模

34、拟滤波器是稳定的，则所有极点si都在S左半平面，即Resi0,那么变换后H(z)的极点 也都在单位圆以内，即| |= 1，因此数字滤波器保持稳定。值得注意的是，这种Ha(s)到H(z)的对应变换关系，只有将Ha(s)表达为部分分式形式才成立。 虽然脉冲响应不变法能保证S平面与Z平面的极点位置有一一对应的代数关系，但这并不是说整个S平面与Z平面就存在这种一一对应的关系，特别是数字滤波器的零点位置与S平面上的零点就没有一一对应关系，而是随着Ha(s)的极点si与系数Ai的不同而不同。H(ej) 是Ha(j)的周期延拓（周期为fs），因Ha(j)并不是带限，即在超过fs频率部分并不为0，所以就产生了

35、混迭。当为低通或带通滤波器时，fs越大，则Ha(j)的下一周期相隔越远，混迭也就越小。当为带阻或高通滤波器时，Ha(j)在超过fs/2频率部分全为通带，这样就不满足抽样定理，发生了完全的混迭，所以脉冲响应不变法不能设计带阻或高通滤波器4。2.3.2脉冲相应不变法优缺点在要求时域脉冲响应能模仿模拟滤波器的场合，一般使用脉冲响应不变法。 脉冲响应不变法的一个重要特点是频率坐标的变换是线性的，与是线性关系。因此如果模拟滤波的频响带限于折叠频率以内的话，通过变换后滤波器的频响可不失真地反映原响应与频率的关系。 例如线性相位的贝塞尔低通滤波器，通过脉冲响应不变法得到的仍是线性相位的低通数字滤波器。 如果

36、Ha(s)是稳定的，即其极点在S左半平面，映射到H(Z)也是稳定的。 脉冲响应不变法的最大缺点：有频谱周期延拓效应，因此只能用于带限的频响特性，如衰减特性很好的低通或带通。而高频衰减越大，频响的混淆效应越小，至于高通和带限滤波器,由于它们在高频部分不衰减，因此将完全混淆在低频响应中。所以用脉冲响应不变法实现高通和带限滤波器时，应增加一保护滤波器，滤掉高于折叠频率以上的频带，然后再用脉冲响应不变法转换为数字滤波器，这会增加设计的复杂性和滤波器的阶数，只有在一定需要频率线性关系或保持网络瞬态响应时才采用4。 2.3.3应用实例【例1】用脉冲响应不变法设计Butterworth低通数字滤波器，要求通

37、带频率为，通带波纹小于1dB，阻带在内，幅度衰减大于15dB，采样周期T=0.01s。假设一个信号，其中f1=5Hz,f2=30Hz。试将原信号与经过该滤波器的输出信号进行比较。%Samp6_6wp=0.2*pi;ws=0.3*pi;Rp=1;Rs=15; %数字滤波器截止频率、通带波纹和阻带衰减T=0.01;Nn=128; %采样间隔Wp=wp/T;Ws=ws/T; %得到模拟滤波器的频率采用脉冲响应不变法的频率转换形式N,Wn=buttord(Wp,Ws,Rp,Rs,s); %计算模拟滤波器的最小阶数z,p,k=buttap(N); %设计低通原型数字滤波器Bap,Aap=zp2tf(z,

38、p,k); %零点极点增益形式转换为传递函数形式b,a=lp2lp(Bap,Aap,Wn); %低通滤波器频率转换bz,az=impinvar(b,a,1/T); %脉冲响应不变法设计数字滤波器传递函数figure(1)H,f=freqz(bz,az,Nn,1/T); %输出幅频响应和相频响应subplot(2,1,1),plot(f,20*log10(abs(H);xlabel(频率/Hz);ylabel(振幅/dB);grid on;subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)xlabel(频率/Hz);ylabel(相位/o);grid on;

39、figure(2)f1=5;f2=30; %输入信号含有的频率N=100; %数据点数n=0:N-1;t=n*T; %时间序列x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t); %输入信号subplot(2,1,1),plot(t,x),title(输入信号)y=filtfilt(bz,az,x); %对信号进行滤波subplot(2,1,2),plot(t,y),title(输出信号),xlabel(时间/s) 程序的运行结果为图2和图3。该例所要求的通带频率为，而该滤波器的采样间隔为0.01s，采样频率为100Hz，即2对应于100Hz，则0.2对应于10Hz，即该滤

40、波器的通带范围为010Hz。观看图6-6上图，其通带范围最大衰减小于1dB，与该分析一致。题意要求阻带在内，幅度衰减大于15dB。其中0.3对应于15Hz，幅度衰减大于15dB。完全符合滤波器设计的要求。测试信号中含有5Hz和30Hz的信号，显然，5Hz可以通过该滤波器，而30Hz的信号不能通过该滤波器。从输出信号可以看出，滤波器滤去了30Hz的高频信号，达到了滤波要求。图-2例1设计的滤波器的频率特性上图：幅频特性；下图：相频特性图-3例1设计滤波器的输入和输出信号42.4 双线性变换法设计IIR数字滤波器2.4.1设计原理双线性变换法的基本思想是首先按给定的指标设计一个模拟滤波器，其次将这

41、个模拟滤波器的系统函数H (s)，通过适当的数学变换方法把无限宽的频带，变换成频带受限的系统函数 H (s)。最后再将H (s)进行常规z变换，求得数字滤波器的系统函数H (z)。这样由于在数字化以前已经对频带进行了压缩，所以数字化以后的频响可以做到无混叠效应。显然，这里寻找压缩频带而又能满足上述映射条件的变换式是个关键。设将s平面映射到s平面存在下列的关系式 （a） 式中s = + j，C为变换常数，在式的右边是以表示的周期函数，其周期为2 / T。如果考虑频率特性则分别以s = j， s = j代入式故得： （b） 现以为纵坐标，=T为横坐标，则模拟与数字频率变量间的关系将如图8.1所示。

42、结合式（a）不难看出；s的左半平面与s的左半平面相对应；s的右半平面与s的右半平面相对应；s平面的虚轴与s平面的虚轴相对应。它们之间主要的区别在于s平面-的无限频率范围，被映射到s平面的主值范围内,即：也就是说，通过变换式（a），把整个s平面映射到s平面以+-s/2为边界的水平窄区内。这避免数字化后可能出现的频谱混叠提供了必要条件。为了求出数字滤波器的系统函数，最后还得通过常规z变换将s平面变换到z平面上来，其关系式为： (c)显见，这时在s左半平面的窄区就被映射到z平面的单位圆内。现将式（c）代入式(a)最后求得： 或 (d)该式是两个线性函数之比，称为线性分式变换，若把它展开求z，则得：

43、（e）可见，其反变换也是线性分式函数，所以这种变换是双向的，因此叫做双线性变换。双线性变换仍然具有将s的左半平面映射到z平面单位圆内;j轴映射到单位圆上的基本性质。因为当=0时，z=1，说明s平面j轴映射到z平面单位圆上。当0时，上式中的分母大于分子z1，说明s左半平面映射到z平面单位圆内。因而它们一一对应，有着单值关系，是一种保角交换。因此一个稳定的模拟滤波器，通过双线性变换只能得到一个（唯一）稳定的数字滤波器。它不象脉冲响应不变法那样，由于在z平面与s平面之间的映射存在着多值关系，以至在s平面上许多不同的值，映射到z平面后都重叠在一个点子上，造成频谱混叠。由此可见，双线性变换法将；频带严格限制在s/2范围内，从根本上消除了频谱混叠。双线性变换法不仅解决了直接z变换中存在的混叠误差，而且还因为它是个代数变换，所以只要将s与z之间存在的简单代数关系，直接代入，模拟滤波器系统函数式中，就可以求得相应的数字滤波器的系统函数和频响。这样，模拟滤