IIR数字滤波器的原理及设计.ppt

第6章IIR数字滤波器的原理及设计6.1概述6.1.1IIR数字滤波器的差分方程和系统函数我们已经知道IIR数字滤波器是一类递归型的线性时不变因果系统，其差分方程可以写为：（6.1）,进行z变换，可得：于是得到IIR数字滤波器的系统函数：(6.2),6.1.2IIR数字滤波器的设计方法对(6.2)式的有理函数的分子、分母多项式进行因式分解，可以得到：(6.3)其中ci为零点而di为极点。H(z)的设计就是要确定系数、或者零极点、，以使滤波器满足给定的性能指标。一般有三种方法。,1.零极点位置累试法IIR系统函数在单位圆内的极点处出现峰值、在零点处出现谷值,因此可以根据此特点来设置H(z)的零极点以达到简单的性能要求。所谓累试，就是当特性尚未达到要求时，通过多次改变零极点的位置来达到要求。当然这种方法只适用于简单的、对性能要求不高的滤波器的设计。,2.借助于模拟滤波器的理论和设计方法来设计数字滤波器模拟滤波器的逼近和综合理论已经发展得相当成熟，产生了许多效率很高的设计方法，很多常用滤波器不仅有简单而严格的设计公式，而且设计参数已图表化，设计起来方便准确。,而数字滤波器就其滤波功能而言与模拟滤波器是相同的,因此，完全可以借助于模拟滤波器的理论和设计方法来设计数字滤波器。在IIR数字滤波器的设计中，较多地采用了这种方法。,3.用优化技术设计系统函数H(z)的系数、或者零极点、等参数，可以采用最优化设计方法来确定。最优化设计法的第一步是要选择一种误差判别准则，用来计算误差和误差梯度等。,第二步是最优化过程，这个过程的开始是赋予所设计的参数一组初值，以后就是一次次地改变这组参数，并一次次计算H(z)的特性与所要求的滤波器的特性之间的误差，当此误差达到最小值时，所得到的这组参数即为最优参数，设计过程也就到此完成。,这种方法能够精确地设计许多复杂的滤波器，但是往往计算很复杂，需要进行大量的迭代运算，故必须借助于计算机，因而优化设计又叫做IIR滤波器的计算机辅助设计(CAD)。第一种方法的算法简单、设计粗糙，在这里不具体讨论了；第三种方法所涉及的内容很多，并且需要最优化理论作为基础，因此在本章中只能作简要介绍；本章将着重讨论用得最多的第二种方法。,6.1.3借助于模拟滤波器的理论和方法的设计原理利用模拟滤波器来设计数字滤波器，要先根据滤波器的性能指标设计出相应的模拟滤波器的系统函数Ha(s)，然后由Ha(s)经变换而得到所需要的数字滤波器的系统函数H(z)。常用的变换方法有冲激响应不变法和双线性变换法。,6.2模拟低通滤波特性的逼近模拟滤波器的设计包括逼近和综合两大部分，其中逼近部分是与数字滤波器的设计有关的。本节要讨论的是，在已知模拟低通滤波器技术指标的情况下，如何设计其系统函数Ha(s)，使其逼近所要求的技术指标。,模拟系统的频率响应Ha(j)是冲激响应ha(t)的傅里叶变换，Ha(j)的模表征系统的幅频特性，下面要讨论如何根据幅频特性指标来设计系统函数。图6.1中用虚线画出的矩形表示一个理想的模拟低通滤波器的指标，是以平方幅度特性|Ha(j)|2来给出的。,c是截止频率，当0c时，|Ha(j)|2=0，是阻带。图6.1中的实的曲线表示一个实际的模拟低通滤波器的平方幅度特性，我们的设计工作就是要用近似特性来尽可能地逼近理想特性。通常采用的典型逼近有Butterworth逼近、Chebyshev逼近和Cauer逼近（也叫椭圆逼近。,6.2.1Butterworth低通滤波特性的逼近对于Butterworth滤波器有：（6.4）满足此平方幅度特性的滤波器又叫做B型滤波器。这里N为正整数，为B型滤波器的阶次，为截止频率。,6.2.1.1B型滤波特性1.最平坦函数B型滤波器的幅频特性是随增大而单调下降的。在=0附近以及很大时幅频特性都接近理想情况，而且在这两处曲线趋于平坦，因此B型特性又叫做最平坦特性。,2.3db带宽由(6.4)式可知，当=c时，=，而因此截止频率又叫做3db带宽或者半功率点。,图6.1Butterworth低通滤波器的平方幅度特性,3.N的影响在通带内，01，故N越大，随增大而下降越快。,因此，N越大，B型滤波器的幅频特性越接近理想的矩形形状；而不同的N所对应的特性曲线都经过c处的半功率点。离c越近，幅频特性与理想特性相差越大。,6.2.1.2由得到Ha(s),B型滤波器的极点由于Ha(s)是s的实系数有理函数，故有：，令s=j,则有：，而(6.5)由(6.4)式和(6.5)式有：用s代替上式中的j:(6.6),图6.2阶次N对B型特性的影响,(6.6)式的极点为：p=0,1,2N-1作为1的2N次方根，p均匀地分布在单位圆上，幅角间隔为/N；它们关于实轴对称，却没有一个在实轴上。显然，将的模乘上，再将其按逆时针方向旋转，就得到sp。因此，sp均匀地分布在半径为的圆周上，其位置关于虚轴对称，却没有一个在虚轴上，这就是说，2N个极点sp在s平面的左、右两半平面各有N个。,这2N个极点是Ha(s)Ha(-s)的极点，考虑到系统函数Ha(s)的极点必须在左半平面系统才是稳定的，因而将左半s平面的N个极点sk（k=0,1,N-1)分给Ha(s)，这样，右半平面的N个极点-sk就正好是Ha(s)的极点。因此有：(6.8),这个式子中的常数是为了使(6.5)式满足而加入的。这N个极点s0、s1、sN-1在s平面的左半平面而且以共轭形式成对出现，当N为奇数时,有一个在实轴上(为-)。,6.2.1.3一般情况下的B型低通滤波器,图6.3一般情况下低通滤波器的设计指标,此时，应该将角频率标称化，通常以1为基准频率，则标称化角频率为：=/1。于是通带边界的标称化角频率为1=1，并且在通带有01，在过渡带和阻带则有1。以下为了方便起见，仍用不带撇的表示标称化的角频率。频率标称化后，B型滤波器的平方幅度特性仍如(6.2)式所示，只是式中的参数和N都需要由图6.3给出的指标来确定。,（6.4）式可以写成：(6.10)当=1=1时，上式为：(6.11)令(6.12)则由(6.11)式可得：,当时有：(6.13)故(6.14)由(6.14)式可求出N，再将其代入(6.12)式，即可求得。,6.4冲激响应不变法本节和下一节所讨论的问题是，在已知模拟滤波器的系统函数Ha(s)的情况下，如何求相应的数字滤波器的系统函数H(z)。s是模拟复频率，Ha(s)也是模拟滤波器的冲激响应ha(t)的拉氏变换。,6.4.1冲激响应不变法的变换方法模拟滤波器的系统函数通常可以表示为：(6.62),而且一般都满足M0,r1；当0,r1。,(6.75)式既表示了数字角频率与模拟角频率之间的关系，也表示了z平面的幅角与s平面的虚部之间的关系。由(6.75)式还可以知道，s平面上由-/Ts到/Ts这一条状区域映射到z平面上由-到的区域，即整个z平面；s平面上的水平线=-/Ts映射到z平面上的射线=-，而当这条射线按逆时针方向旋转时，对应的s平面上的水平线就向上平移。,上面所阐述的不仅是模拟域s平面与数字域z平面之间的映射关系，而且也是模拟滤波器的频率与用冲激响应不变法所得到的数字滤波器的频率之间的关系。s平面与z平面的映射关系保证了将稳定的模拟滤波器变换为稳定的数字滤波器。,图6.14模拟复频率s与数字复频率z之间的映射关系,例6.6用冲激响应不变法设计一个三阶Butteworth数字低通滤波器，抽样频率为fs=1.2kHz,截止频率为=400Hz。解：此数字滤波器的截止频率：c=2fc=2400=800弧度/s这也是模拟滤波器的截止频率，于是可以写出模拟滤波器的系统函数:,其中，现在进行部分分式分解，令(*2),可以得到：根据(*1)式和(*2)式，再将A、B、C代入，便得到：,上式中Ts=1/fs=1/1200（秒）。,6.5双线性变换法6.5.1双线性变换关系的导出模拟滤波器的系统函数可以变换为：这里为了方便说明，已令M=N。,由此式可以看出，模拟滤波器的基本单元是积分器，因此，只要设法用某种数字网络来代替此基本单元，就能够将模拟滤波器转变成相应的数字滤波器。模拟滤波器基本单元的系统函数为：则其冲激响应为：,设有一信号(t0）输入到该积分器系统，则其输出也即对的响应为：设01,即单位圆外=0，即虚轴r=1,即单位圆0，即左半平面r1,即单位圆内,因此，用双线性变换法，稳定的模拟滤波器导出的数字滤波器也必定是稳定的。但是，与冲激响应不变法不同的是，在双线性变换下，模拟滤波器的复频率s与相应的数字滤波器的复频率z之间的映射是一一对应的关系。,图6.16双线性变换法s平面与z平面之间的映射关系,6.5.3频率预畸变下面讨论s平面的虚轴与z平面的单位圆的映射关系，也即模拟滤波器的角频率与相应的数字滤波器的角频率之间的关系。在(6.84)式中令=0，便可得到：或（6.85）,图6.17与之间的非线性关系,与的关系是非线性的，但是，s平面上的虚轴一一对应地映射到了z平面单位圆的一周之上，因此，采用双线性变换法，不存在频域混叠失真的问题。由双线性变换所引起的模拟滤波器频率与数字频率之间的非线性关系，使得所得到的数字滤波器的相位频率特性产生失真；,但对于幅度频率特性，可以通过频率预畸变来校正。实际上，只要首先根据所要求的数字滤波器的各关键频率，按照(6.85)式转变成相应的模拟频率，再根据这些频率指标来设计模拟滤波器，则最后转换成的数字滤波器的各关键频率就会正好映射到所要求的位置上。,6.5.4双线性变换法的特点1模拟滤波器经过双线性变换后，不存在频率特性的混叠失真，因而对模拟滤波器的频率响应函数Ha()无限带要求，而且能够直接用于设计低通、高通、带通、带阻等各种类型的数字滤波器。,2与冲激响应不变法中模拟频率与数字频率之间的线性关系=Ts不同的是，双线性变换法中模拟滤波器的频率与所转换成的数字滤波器的频率之间是非线性关系，但是，如果事先进行频率预畸变，这种非线性关系不会使所设计的数字滤波器的幅频特性受到影响。3双线性变换方法比较容易，不需要将模拟系统函数进行部分分式分解。,因此，双线性变换法是用得很普遍、并且很有效的一种方法；只是，由于频率的非线性关系会产生相频特性失真，所以若对数字滤波器的相位特性要求较严，则不宜采用这种变换方法。,最后必须强调说明一下用双线性变换法来设计数字滤波器时各种频率之间的关系。我们在考虑一个数字滤波器的频域特性时，所采用的频率变量可以是数字频率，也可以是模拟频率。模拟角频率=2f，f是以赫兹(Hz)为单位的真正具有物理意义的频率变量。与的关系为=Ts，Ts为抽样周期。,数字滤波器的频率响应=，是的周期函数，以2为周期；是的周期函数，周期为s=2/Ts。上述这些关系与数字滤波器的设计方法无关。如果数字滤波器是用冲激响应不变法设计的，则模拟滤波器的频率变量也就是数字滤波器的模拟频率变量；如果数字滤波器是用双线性变换法来设计的，那末模拟滤波器的频率变量并不是数字滤波器的模拟频率变量。,我们在6.5.2节中所述的双线性变换法s平面与z平面的映射关系实际上是被变换的模拟滤波器的复频率s与所得到的数字滤波器的复频率z之间的关系，(6.85)式中的也是此模拟滤波器的角频率，并不是数字滤波器的模拟角频率。为了便于区分，应该将(6.85)式中的模拟滤波器角频率用来表示，即为：；而数字滤波器的数字角频率与其本身的模拟角频率之间仍然是上面所述的那种线性关系，即有：=Ts。,