1、第4章 无限长单位脉冲响应数字滤波器设计第4章 无限长单位脉冲响应数字滤波器设计4.1 数字滤波器的基本概念4.2 模拟滤波器的设计4.3 利用模拟滤波器设计IIR数字滤波器4.4 从模拟低通原型到各种数字滤波器的频率变换第4章 无限长单位脉冲响应数字滤波器设计4.1 数字滤波器的基本概念数字滤波器的基本概念 许多信息处理过程，如通信系统中的信号检测、 预测、 变频、 调制等都要用到滤波器。数字滤波器是数字信号处理系统中最常用的一种线性系统，是数字信号处理的重要基础。 所谓数字滤波器是指输入、 输出均为数字信号，通过一定运算关系改变输入信号的频率或者相位特性的器件或者软件。数字滤波的概念和模拟
2、滤波基本相同，只是处理的信号形式不同，滤波实现方法不同。数字滤波器与模拟滤波器相比，具有精度高、 稳定、 体积小、 重量轻、 灵活、 不需要阻抗匹配、 能够实现模拟滤波器无法实现的特殊功能等优点。模拟信号可以通过A/DC、 D/AC转换匹配后，采用数字滤波器实现滤波。因此，数字滤波器得到广泛应用。第4章 无限长单位脉冲响应数字滤波器设计4.1.1 数字滤波器的分类数字滤波器的分类数字滤波器有不同的分类方法，主要有以下三种：(1) 根据滤波器的性能来分，可以分为经典数字滤波器和现代数字滤波器。经典数字滤波器即一般滤波器，特点是有用信号和干扰信号各占不同的频带，通过一个合适的滤波器选择出有用的频率
3、信号，滤除干扰信号，因此也称为选频滤波器。现代数字滤波器的特点是针对信号和干扰的频带相互重叠时，按照随机信号内部的一些统计分布规律，从干扰中提取有用信号，例如，维纳滤波器、 卡尔曼滤波器、 自适应滤波器等最佳滤波器。本课程仅介绍经典数字滤波器。第4章 无限长单位脉冲响应数字滤波器设计(2) 根据数字滤波器的功能来分，和模拟滤波器一样，可以分为低通、 高通、 带通、 带阻、 多带数字滤波器。它们的理想幅度特性如图4.1.1。从图4.1.1中可以看出，数字滤波器的传输函数H(ej)都是以2为周期的，滤波器的低通频带位于2的整数倍附近，而高通频带在的奇数倍附近。(3) 根据实现的网络结构或者从单位脉
4、冲响应来分，可以分成无限脉冲响应(IIR)滤波器和有限脉冲响应(FIR)滤波器。它们的系统函数分别为：(4.1.1)(4.1.2)第4章 无限长单位脉冲响应数字滤波器设计式(4.1.1)表示的滤波器为单位脉冲响应无限长，有反馈支路，即含有环路，也称为递归型；阶数为N，因此称为N阶IIR(Infinite Impulse Response)数字滤波器。式(4.1.2)表示的滤波器为单位脉冲响应有限长，无反馈支路，也称为非递归型；阶数为M，因此称为M阶FIR(Finite Impulse Response)数字滤波器。第4章 无限长单位脉冲响应数字滤波器设计图 4.1.1 理想低通、 高通、 带通
5、、 带阻、 多带数字滤波器幅度特性第4章 无限长单位脉冲响应数字滤波器设计4.1.2 数字滤波器的技术指标数字滤波器的技术指标常用的数字滤波器属于选频滤波器。假设数字滤波器的传输函数|H(ej)|用下式表示：(4.1.3)对于一般作为选频滤波器的IIR数字滤波器，其技术指标是由幅频特性给出的，而相频特性不作要求。本章主要研究的IIR数字滤波器的技术指标就是由幅频特性提出的。下面以低通IIR数字滤波器为例，介绍数字滤波器的技术指标。图4.1.2为低通数字滤波器的幅度特性，图中：0p为通带， p称为通带截止频率； p r为过渡带，r为阻带，r称为阻带截止频率。第4章 无限长单位脉冲响应数字滤波器设
6、计图 4.1.2 低通数字滤波器的幅度特性第4章 无限长单位脉冲响应数字滤波器设计通带内的技术指标为11|H(ej)|1 p (4.1.4)阻带内的技术指标为：|H(ej)|2r (4.1.5)通带内波动用dB数表示：(4.1.6)阻带内最小衰减用dB数表示：(4.1.7)第4章 无限长单位脉冲响应数字滤波器设计4.1.3 数字滤波器设计方法概述数字滤波器设计方法概述1) 直接设计法直接设计法是直接在频域或时域进行逼近设计，具体方法有零极点累试法、 频域逼近法、 时域逼近法，通常须借助计算机进行优化设计，因此也称为最优化设计方法。具体步骤分两步：(1) 确定一种最优准则，如最小均方误差准则，即
7、使设计出的实际频率响应的幅度特性与所要求的理想频率响应的均方误差最小，如:(4.1.8)第4章 无限长单位脉冲响应数字滤波器设计(2) 在此最佳准则下，通过不断地迭代运算求滤波器的系数，直到满足要求为止。此外还有其他多种误差最小准则。 随着计算机技术的发展，最优化设计方法的使用也在逐渐增多。但是，由于涉及到最优化设计算法，而且目前已经有商业设计程序可以使用，因此本章不介绍。第4章 无限长单位脉冲响应数字滤波器设计2) 间接设计法间接设计法是借助于模拟滤波器的设计方法进行的，其设计步骤是：先设计模拟滤波器得到传输函数Ha(s)，然后将Ha(s)按某种方法转换成数字滤波器的系统函数H(z)。由于模
8、拟滤波器的设计方法已经发展得很成熟，而且模拟滤波器有简单而严格的设计公式，设计起来方便、 准确，因此可将这些理论推广应用到数字域，作为设计数字滤波器的工具。本章后续内容将介绍这种设计方法。第4章 无限长单位脉冲响应数字滤波器设计4.2 模拟滤波器的设计模拟滤波器的设计 模拟滤波器的理论和设计方法已发展得相当成熟，且有若干典型的模拟滤波器供我们选择，如巴特沃斯滤波器、 切比雪夫滤波器、 椭圆滤波器、 贝塞尔滤波器等，这些滤波器都有严格的设计公式、 现成的曲线和图表供设计人员使用，因此设计方便。 模拟滤波器根据幅度特性可以分为低通、 高通、 带通和带阻滤波器，它们的理想特性如图4.2.1所示。由于
9、模拟滤波器的设计总是先设计低通滤波器，再通过频率变换将低通变换为希望类型的滤波器，因此，下面先介绍模拟低通滤波器的技术指标和逼近方法，然后介绍模拟滤波器的频率变换。 第4章 无限长单位脉冲响应数字滤波器设计图 4.2.1 各种理想模拟滤波器的幅频特性第4章 无限长单位脉冲响应数字滤波器设计4.2.1 模拟低通滤波器的技术指标及逼近方法模拟低通滤波器的技术指标及逼近方法 模拟低通滤波器的归一化特性曲线如图4.2.2所示，图中：p为通带截止频率；c为幅度下降到0.707(也就是下降3 dB)时的频率，因此通常被称为3 dB截止频率；r为阻带截止频率。与数字低通滤波器相似，模拟低通滤波器的设计指标也
10、用4个参数描述：、 p、 At和r。其中：是通带(= 0p)中的最大衰减系数；At是阻带r的最小衰减系数，和At一般用dB数表示，可表示成：(4.2.1)(4.2.2)第4章 无限长单位脉冲响应数字滤波器设计图 4.2.2 模拟低通滤波器的幅频特性曲线第4章 无限长单位脉冲响应数字滤波器设计如果模拟低通滤波器的特性是单调下降，则可表示为:=10 lg|Ha(jp)|2(4.2.3)At=10 lg|Ha(jr)|2 (4.2.4)滤波器的技术指标给定后，模拟滤波器的设计就是构造一个模拟系统函数Ha(s)，使其特性按照技术指标的要求去逼近某个理想滤波器特性。 对于一般滤波器的单位脉冲响应ha(t
11、)为实数的因果系统有(4.2.5)(4.2.6)将式(4.2.5)的右边展开得第4章 无限长单位脉冲响应数字滤波器设计由于ha(t)为实数，其频谱具有共轭对称性，因此有Ha(j)=H*a(j)(4.2.7)定义幅度平方函数：A(2)=|Ha(j)|2=Ha(j)H*a(j) (4.2.8)将式(4.2.7)代入式(4.2.8)得 A(2)=Ha(j)Ha(j)=Ha(s)Ha(s)|s=j (4.2.9)又由于s=j，2=s2，因此A(2)=A(s2)|s=j第4章 无限长单位脉冲响应数字滤波器设计对于给定的A(s2)，先在S复平面上标出A(s2)的极点和零点，由(4.2.9)式知，A(s2)
12、的极点和零点总是“成对出现”，且对称于S平面的实轴和虚轴，选用A(s2)的对称极、 零点的任一半作为Ha(s)的极、 零点，则可得到Ha(s)。 为了保证Ha(s)的稳定性，应选用A(S2)在S左半平面的极点作为Ha(s)的极点，零点可选用任一半。 可见只要给定幅度平方函数模型，模拟滤波器就可以设计出来，下面介绍三种模拟低通滤波器的设计。第4章 无限长单位脉冲响应数字滤波器设计4.2.2 巴特沃兹滤波器巴特沃兹滤波器巴特沃兹(Butterworth)模拟低通滤波器的幅度平方函数的表达式为(4.2.10)式中N为滤波器的阶数。巴特沃兹滤波器的幅度特性曲线如图4.2.3所示。由图可以看出，它的特点
13、是：通带内具有最大平坦的幅度特性，且随频率的增加，幅频特性单调下降。随着N的增加，通带和阻带的近似性越好，过渡带越陡。通带内，分母/c1，(/c)2N远远大于1，随着增加，A(2)快速减小。第4章 无限长单位脉冲响应数字滤波器设计当=c时，A(2)=1/2，因此，幅度衰减相当于3 dB衰减点。可以证明：随着N的增加，通带内频率响应变得更为平坦，阻带内的衰减更大；过渡带内频响更趋于斜率为6N dB/倍频程的渐近线。 由式(4.2.10)可以求得幅度平方函数的极点：式中：k=0，1，2，2N1 。可见，巴特沃兹滤波器的振幅平方函数有2N个极点，它们均匀对称地分布在|S|=c的圆周上。例如，N=3阶
14、时巴特沃兹幅度平方函数的极点分布如图4.2.4所示。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.4 三阶A(s2)的极点分布第4章 无限长单位脉冲响应数字滤波器设计考虑到系统的稳定性，可知模拟低通滤波器的系统函数应由S平面左半部分的极点(SP0, SP1, SP2)组成，它们分别为则式中：，令，称为归一化频率；再令一复变量p=j，得归一化的传输函数：第4章 无限长单位脉冲响应数字滤波器设计式中：，k=0，1，2，(2N1)，称为归一化的极点；它只与N有关，因此，当N确定时，归一化的极点可以事先算出并建立一个参数表，如表4.2.1所示。第4章 无限长单位脉冲响应数字滤波器设计第4章 无限
15、长单位脉冲响应数字滤波器设计将式(4.2.12)分母的因式展开为多项式，则Ha(p)可以写成以下形式：式中的系数bk(k=0，1，2，2N1)只与归一化极点有关，它可以事先算出并建立一个参数表，如表4.2.2所示。(4.2.13)第4章 无限长单位脉冲响应数字滤波器设计第4章 无限长单位脉冲响应数字滤波器设计这样，当N确定时，通过查表就可以得到滤波器的归一化原型，例如：当N=3时，查表得(4.2.14)令p=s/c，并去归一化后，得到所要设计的巴特沃兹模拟低通滤波器的传递函数为(4.2.15)因此，只要能够确定阶数N和3 dB截止频率c，就可以得到所要设计的滤波器。下面讨论如何确定巴特沃兹模拟
16、低通滤波器的阶数N和3 dB截止频率c。第4章 无限长单位脉冲响应数字滤波器设计巴特沃兹模拟低通滤波器的阶数N是根据滤波器的技术指标来确定的。将式(4.2.10)代入式(4.2.3)得整理得同样将式(4.2.10)代入式(4.2.4)，得将式(4.2.16)和式(4.2.17)联立求解，则有(4.2.16)(4.2.17)(4.2.18)第4章 无限长单位脉冲响应数字滤波器设计 计算出N后，由式(4.2.17)可以求得c的计算公式：(4.2.20)在这种情况下，设计结果通带指标满足要求，阻带有富余。 当然，也可以由式(4.2.18)求得c的计算公式：(4.2.21)在这种情况下，设计结果阻带指
17、标满足要求，通带有富余。在利用巴特沃斯模型设计模拟低通滤波器时，是选择式(4.2.20)还是选择式(4.2.21)，应根据实际情况而定。第4章 无限长单位脉冲响应数字滤波器设计【例【例4.2.1】已知模拟低通滤波器通带截止频率fp=0.5 Hz, 通带最大衰减=2 dB，阻带截止频率fr=1.2 Hz，阻带最小衰减At=30 dB。试采用巴特沃斯模型设计该模拟低通滤波器。解解 利用上述的方法设计，设计步骤为： (1) 确定模拟低通滤波器的阶数N。 (2) 由N=5，查表4.2.2得到5阶归一化低通原型滤波器：第4章 无限长单位脉冲响应数字滤波器设计(3) 求3 dB截止频率c。按照(4.2.2
18、0)式，得通带刚好满足要求，而阻带有富余。按照(4.2.21)式，得阻带刚好满足要求，而通带有富余。(4) 去归一化，得到所设计的低通滤波器：第4章 无限长单位脉冲响应数字滤波器设计上述设计直接计算相当复杂，可以采用MATLAB设计。MATLAB提供了三个函数用于ButterWorth模拟滤波器的设计。(1) 阶数及截止频率的求取函数buttord。buttord函数的语句格式为： n，Wn = buttord( Wp，Wr，Rp，Rr，s) 其中：Wp为通带截止频率，Wr为阻带截止频率，单位是rad/s；Rp为通带最大衰减；Rr为阻带最小衰减，单位是dB。 选项“s”表示设计模拟滤波器，省略
19、此参数为设计数字滤波器。当所设计的滤波器为带通、 带阻滤波器时，Wp、 Wr为2元数组。函数的返回值Wn为截止频率，n为满足技术指标的滤波器的最小阶数。第4章 无限长单位脉冲响应数字滤波器设计(2) 设计函数butter。butter函数的语句格式有两种： 当设计低通或者带通(带通的 Wn为2元数组)滤波器时，采用以下格式： b，a = butter(n，Wn，s) 当设计的是其他类型的滤波器时，可采用ftype说明滤波器的类型，语句格式为： b，a = butter(n，Wn，ftype，s)式中：ftype可以为“high”、 “low”、 “stop”，分别表示设计的滤波器为高通、 低通
20、、 带阻滤波器。函数的返回参数b，a为模拟滤波器的传递函数模型，可以根据下式写出传递函数：第4章 无限长单位脉冲响应数字滤波器设计式中：a(0)=1。 (3) 归一化模拟低通原型设计函数。语句格式为： z，p，k=buttap(n)n为阶数，z为零点，p为极点，k为放大系数。第4章 无限长单位脉冲响应数字滤波器设计【例【例4.2.2】 对于例4.2.1，这里给出采用MATLAB设计的程序代码。先设计模拟低通原型，然后通过频带转换设计模拟低通滤波器。程序设计如下：fp=0.5*2*pi；ap=2；fr=1.2*2*pi；as=30；n，Wn = buttord(fp，fr，ap，as，s)z，p
21、，k = buttap(n)；%低通原型设计b0=k*real(poly(z)；%poly函数用于把多项式根转换为多项式系数，real用于复数取实部运算a0=real(poly(p)；b，a=lp2lp(b0,a0,Wn) %模拟低通原型转换为模拟低通滤波器h，f=freqs(b，a)； %求模拟滤波器传输函数的频率响应mag = abs(h)；mag = 20*log10(mag)；phase = angle(h)；第4章 无限长单位脉冲响应数字滤波器设计subplot(2，1，1)；plot(f/(2*pi)，mag)；title(N=5 Butterworth Lowpass Filte
22、r)；axis(0 1.5 -35 0)；xlabel(f(Hz)；ylabel(幅度(dB)；grid；subplot(2，1，2)；plot(f/(2*pi)，phase)；axis(0 1.5 -4 4)；xlabel(f(Hz)；ylabel(相位(ard)；第4章 无限长单位脉冲响应数字滤波器设计也可跳过原型设计，直接设计模拟滤波器，程序设计如下：fp=0.5*2*pi；ap=2；fr=1.2*2*pi；as=30；n，Wn = buttord(fp，fr，ap，as，s)b，a = butter(n，Wn，s)h，f=freqs(b，a)；其他幅频特性、 相频特性的绘图程序与前面相
23、同，省略。程序运行结果：(1) 滤波器的阶数N = 5，与上述计算结果相同。(2) 截止频率Wn = 3.7792rad/sec= 0.6015 Hz，与按照式(4.2.21)计算的结果相同。第4章 无限长单位脉冲响应数字滤波器设计(3) 滤波器的传递函数，由程序运行结果可得：b = 770.9440a = 1.0000 12.2299 74.7850 282.6305 660.1398 770.9440则所设计的滤波器的传递函数为(4) 滤波器的特性分析结果如图4.2.5所示。图中上半部分的图形为幅频特性曲线，可以看出在通带截止频率0.5 Hz处的衰减小于2 dB，阻带截止频率1.2 Hz处
24、阻带最小衰减大于30 dB，满足指标要求。图中下半部分的图形为相位特性曲线，可以看出在通带内相位特性接近线性。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.5 巴特沃斯模型模拟低通滤波器特性第4章 无限长单位脉冲响应数字滤波器设计4.2.3 切比雪夫滤波器切比雪夫滤波器切比雪夫(chebyshev)滤波器具有波动性。它有两种类型：切比雪夫型，在通带中具有等波动响应；切比雪夫型，在阻带中具有等波动响应。切比雪夫型滤波器的振幅平方函数为(4.2.22)式中：表示|Ha(j)|波动范围的参数，VN(x)是一个N阶切比雪夫多项式，定义为(4.2.23)第4章 无限长单位脉冲响应数字滤波器设计由式
25、(4.2.23)可以看出，当|x|1时，|VN(x)|1，并按照余弦函数波动；当|x|1时，随着|x|增大，VN(x)将按双曲余弦函数快速增大。切比雪夫滤波器的振幅平方特性如图4.2.6所示。在通带内，由于c，满足|x|1，因此，切比雪夫滤波器的振幅平方函数|Ha(j)|2在1(1/1+2)范围内变化。当c时，随着/c增大，|x|1，因此|Ha(j)|2迅速趋于零。当=0时，由式(4.2.22)和式(4.2.23)可得第4章 无限长单位脉冲响应数字滤波器设计当N为奇数时，由于，可得|Ha(j)|2|=0=1，如图4.2.6(a)所示。当N为偶数时，由于 ,可得 ，如图4.2.6(b)所示。切比
26、雪夫型滤波器在通带范围内是等幅起伏的，所以同样的通带衰减，其阶数比巴特沃斯滤波器要小。可根据需要对通带内允许的衰减量(波动范围)提出的要求进行设计，然而由于切比雪夫多项式计算比较复杂，设计计算量比较大，因此，通常采用工具软件来设计。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.6 切比雪夫滤波器的振幅平方特性第4章 无限长单位脉冲响应数字滤波器设计切比雪夫型的平方幅度响应为切比雪夫型滤波器通带幅频特性是单调的，而阻带是等波动的。这种滤波器在S平面上既有极点，又有零点。它的计算比型更加复杂，因此，掌握MATLAB计算即可。MATLAB提供多个函数用于切比雪夫模拟滤波器的设计。下面以型为例进
27、行应用介绍(只要把函数中的“1”改为“2”，即为切比雪夫型滤波器设计函数)。(1) 阶数求取函数cheb1ord。cheb1ord函数的语句格式为： n，Wn = cheb1ord( Wp，Wr，Rp，Rr，s) 其中：Wp、 Wr为角频率，单位是rad/s。第4章 无限长单位脉冲响应数字滤波器设计当所设计的滤波器为带通、 带阻滤波器时，Wp、 Wr为2元数组。函数的返回值Wn为截止频率，n为满足技术指标的滤波器的最小阶数。(2) 设计函数cheby1。cheby1函数的语句格式有两种： 当设计低通或者带通(带通的 Wn为2元数组)滤波器时，采用以下格式： b，a = cheby1(n，Rp，
28、Wn，s) 当设计的是其他类型的滤波器时，可采用ftype说明滤波器的类型，语句格式为： b，a = cheby1(n，Rp，Wn，ftype，s)式中：ftype可以为high、 low、 stop，分别表示设计的滤波器为高通、 低通、 带阻滤波器。第4章 无限长单位脉冲响应数字滤波器设计函数的返回参数b，a为模拟滤波器的传递函数模型，可以根据下式写出传递函数:式中：a(0)=1。 (3) 归一化模拟低通原型设计函数。语句格式如下： z，p，k=cheb1ap(n，Rp)n为阶数，Rp为通带最大衰减，z为零点，p为极点，k为放大系数。【例【例4.2.3】 下面给出例4.2.1的MATLAB切
29、比雪夫型设计程序，供读者参考。先设计低通原型，然后通过频带转换设计模拟低通滤波器。程序设计如下：第4章 无限长单位脉冲响应数字滤波器设计fp=2*pi*0.5；Ap=2；fr=2*pi*1.2；Ar=30；N，Wn = cheb1ord(fp，fr，Ap，Ar，s)z，p，k = cheb1ap(N，Ap)；%低通原型A，B，C，D = zp2ss(z，p，k)；%零极表达式转换为状态空间表达式At，Bt，Ct，Dt=lp2lp(A，B，C，D，Wn)；%低通原型转换为低通滤波器b，a=ss2tf(At，Bt，Ct，Dt)；%状态空间表达式转换为传输函数表达式h，f=freqs(b，a)；ma
30、g = abs(h)；mag = 20*log10(mag)；phase = angle(h)；subplot(2，1，1)；第4章 无限长单位脉冲响应数字滤波器设计plot(f/(2*pi)，mag)；title(Cheby1 Lowpass Filter)；axis(0 1.5 -35 0)；xlabel(f(Hz)；ylabel(幅度(dB)；grid；subplot(2，1，2)；plot(f/(2*pi)，phase)；axis(0 1.5 -4 4)；xlabel(f(Hz)；ylabel(相位(ard)；第4章 无限长单位脉冲响应数字滤波器设计上述程序运行输出的幅频和相频特性曲线
31、如图4.2.7所示。根据程序运行结果可得：(1) 滤波器的阶数N=3，小于巴特沃斯滤波器的阶数。(2) 截止频率fn=3.1416 ard/s=0.5 Hz。(3) 滤波器的传递函数可由： b=0 0 0 10.1356 a=1 2.3179 10.0886 10.1356得到所设计的滤波器的传递函数为第4章 无限长单位脉冲响应数字滤波器设计(4) 由图4.2.7上半部分的幅频特性曲线可以看出在通带内等波纹波动并且在截止频率0.5 Hz处的衰减小于2 dB，阻带截止频率1.2 Hz处阻带最小衰减大于30 dB，满足指标要求；由图4.2.7中下半部分的相位特性曲线可以看出在通带内相位特性接近线性
32、。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.7 切比雪夫模型模拟低通滤波器特性第4章 无限长单位脉冲响应数字滤波器设计实际上，也可直接采用MATLAB函数设计模拟低通滤波器，将上述代码中的设计部分改为：fp=0.5*pi；ap=2；fr=1.2*2*pi；as=30n，Wn=cheblord(fp，fr，ap，as，s)b，a=chebyl(n，fp，Wn，s)h，f=freqs(b，a)请读者修改上述程序，加入绘制幅频特性、 相频特性曲线的程序，调试、 运行程序，并与前面结果比较，看是否相同。第4章 无限长单位脉冲响应数字滤波器设计4.2.4 椭圆滤波器椭圆滤波器(考尔滤波器考尔滤
33、波器)椭圆滤波器(Elliptic Filter)的振幅平方函数为(4.2.24)式中：RN(，L)为N阶雅可比椭圆函数，L是一个表示波纹性质的参量。图4.2.8为N=5时，R25(，L)的特性曲线，由图可见，在归一化通带内(11)，R25(，L)在(0，1)间振荡；而在超过L后，就在L2间振荡。这一特点使滤波器同时在通带和阻带内可以具有任意衰减量，并且在通带和阻带内都具有等波纹波动的特性。 第4章 无限长单位脉冲响应数字滤波器设计图 4.2.8 R25(，L)的特性曲线第4章 无限长单位脉冲响应数字滤波器设计图4.2.29为椭圆滤波器的振幅平方函数，图中和A的定义同切比雪夫滤波器。 其特点是
34、：幅值响应在通带和阻带内都是等波纹的，对于给定的阶数和给定的波纹要求，椭圆滤波器能获得较其他滤波器更窄的过渡带宽，就这点而言，椭圆滤波器是最优的。但是其设计计算同样很复杂，因此通常采用MATLAB设计。下面给出相应的设计函数： N，Wn = ellipord( Wp，Ws，Rp，Rs，s)b，a=ellip(N，Rp，Rs，Wn，s)z，p，k = ellipap( n，Rp，Rs)详细的设计程序请读者自己编写。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.9 椭圆滤波器的振幅平方函数第4章 无限长单位脉冲响应数字滤波器设计4.2.5 模拟高通、模拟高通、 带通、带通、 带阻滤波器设计带
35、阻滤波器设计1. 低通到高通的频率变换低通到高通的频率变换假设归一化的模拟低通滤波器G(j)和高通滤波器H(j)的幅度特性如图4.2.10所示。图中：p、 r分别为低通滤波器的归一化通带截止频率和归一化阻带截止频率，p、 r分别为高通滤波器的归一化通带截止频率和归一化阻带截止频率。由于|G(j)|和|H(j)|都是频率的偶函数，可以把|G(j)|左边的曲线和|H(j)|曲线对应起来，则低通滤波器的特性中从经过r、 p到0时，对应于高通滤波器的特性中从0经过r、p到。 因此和之间的关系为(4.2.25)第4章 无限长单位脉冲响应数字滤波器设计上式就是低通到高通的频率变换公式。如果已知低通滤波器G
36、(j)，则高通滤波器H(j)可用下式转换求解：(4.2.26)第4章 无限长单位脉冲响应数字滤波器设计图 4.2.10 低通与高通滤波器的幅度特性第4章 无限长单位脉冲响应数字滤波器设计模拟高通滤波器的设计步骤归纳如下：(1) 确定高通滤波器的技术指标：通带截止频率p，阻带截止频率r，通带内最大衰减，阻带最小衰减At。(2) 确定相应低通滤波器的设计指标。按照(4.2.25)式，将高通滤波器的边界频率转换成低通滤波器的边界频率，各项设计指标为：低通滤波器通带截止频率p=1/p；低通滤波器阻带截止频率r=1/r；通带最大衰减仍为，阻带最小衰减仍为At。(3) 设计归一化低通滤波器G(p)。第4章
37、 无限长单位脉冲响应数字滤波器设计(4) 求模拟高通的H(s)。将式(4.2.25)代入G(p)，得到归一化高通滤波器H(q)；为了去归一化，将q=s/c代入H(q)中，合并上述两步得由归一化低通原型到高通滤波器的变换关系：(4.2.27)第4章 无限长单位脉冲响应数字滤波器设计2. 低通到带通的频率变换低通到带通的频率变换假设归一化的模拟低通滤波器(G(j)和带通滤波器(H(j)的幅度特性如图4.2.11所示。图中：p、r分别为低通滤波器的归一化通带截止频率和归一化阻带截止频率；l、 u分别为带通滤波器的下通带截止频率和上通带截止频率；r1、r2分别为带通滤波器的下阻带截止频率和上阻带截止频
38、率；0为带通滤波器的中心频率。定义20=lu，称0为通带中心频率；定义B=ul，称B为通带带宽，一般B作为归一化参考频率。l、 u分别为带通滤波器的归一化下通带截止频率和归一化上通带截止频率，r1、r2分别为带通滤波器的归一化下阻带截止频率和归一化上阻带截止频率，0为带通滤波器的归一化中心频率。第4章 无限长单位脉冲响应数字滤波器设计带通滤波器归一化边界频率用以下式子计算：(4.2.28)将带通和低通滤波器的特性对应起来，得到与的对应关系如表4.2.3所示。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.11 带通与低通滤波器的幅度特性第4章 无限长单位脉冲响应数字滤波器设计第4章 无限长
39、单位脉冲响应数字滤波器设计由图4.2.11或表4.2.3，可得带通滤波器频率转换为归一化低通滤波器频率的关系式：(4.2.29)验证式(4.2.29)：例如，由表4.2.3知p对应u，代入式(4.2.29)中，得式(4.2.29)称为带通到低通的频率变换公式。利用该式可将带通的边界频率转换成低通的边界频率。下面推导由归一化低通传输函数到带通传输函数的转换公式。第4章 无限长单位脉冲响应数字滤波器设计p表示归一化低通的拉普拉斯变量，q表示归一化带通的拉普拉斯变量，s表示带通的拉普拉斯变量，因此p=j，q=j，s=j，代入式(4.2.29)，得将q=j代入，得到归一化低通到归一化带通的变换公式:将
40、带通去归一化，即将q=s/B代入，得因此(4.2.31)(4.2.30)第4章 无限长单位脉冲响应数字滤波器设计3. 低通到带阻的频率变换低通到带阻的频率变换低通与带阻滤波器的幅频特性如图4.2.12所示。图中，u、l分别为带阻滤波器通带上截止频率和通带下截止频率；令B=ul，称B为阻带宽度，一般B作为归一化参考频率。r1、r2分别为带阻滤波器的下阻带截止频率和上阻带截止频率。另外定义0为阻带中心频率，20=lu。相应的归一化边界频率为与式(4.2.28)完全相同。第4章 无限长单位脉冲响应数字滤波器设计图 4.2.12 低通与带阻滤波器的幅频特性第4章 无限长单位脉冲响应数字滤波器设计将带阻
41、滤波器和低通滤波器的幅度特性对应起来，便得到与的对应关系如表4.2.4所示。根据与的对应关系，可得且ul=1。p=1。式(4.2.32)称为低通到带阻的频率变换公式。将(4.2.32)式代入p=j，并去归一化，可得(4.2.32)(4.2.33)则将归一化的模拟低通滤波器变换为模拟带阻滤波器的变换公式为(4.2.34)第4章 无限长单位脉冲响应数字滤波器设计第4章 无限长单位脉冲响应数字滤波器设计下面总结设计模拟带阻滤波器的步骤：(1) 确定模拟带阻滤波器的技术要求：下通带截止频率l，上通带截止频率u，阻带下截止频率r1，阻带上截止频率r2，0为阻带中心频率，20=lu，阻带宽度B=ul。(2
42、) 确定归一化模拟低通技术指标:(3) 设计归一化模拟低通滤波器G(p)。(4) 由(4.2.34)式直接将G(p)转换成带通H(s)。 第4章 无限长单位脉冲响应数字滤波器设计4. 由模拟低通原型设计模拟高通、由模拟低通原型设计模拟高通、 带通、带通、 带阻滤波器的带阻滤波器的MATLAB函数函数 设计模拟低通、 高通、 带通、 带阻滤波器时，可先设计模拟低通滤波器原型，然后通过频率变换实现模拟低通、 高通、 带通、 带阻滤波器。MATLAB提供的函数如下：(1) 低通原型转换为低通：NUMT，DENT = LP2LP(NUM，DEN，Wo)，Wo为截止频率(rad/s)。(2) 低通原型转
43、换为高通：NUMT，DENT = LP2hP(NUM，DEN，Wo)，Wo为截止频率(rad/s)。(3) 低通原型转换为带通：NUMT，DENT = LP2bP(NUM，DEN，Wo，Bw)，Wo为中心频率(rad/s)，Bw为通带带宽(rad/s)。(4) 低通原型转换为带阻：NUMT，DENT = LP2bs(NUM，DEN，Wo，Bw) ，Wo为中心频率(rad/s)，Bw为阻带带宽(rad/s)。第4章 无限长单位脉冲响应数字滤波器设计5. 模拟高通、模拟高通、 带通、带通、 带阻滤波器直接设计的带阻滤波器直接设计的MATLAB函函数数为了简化设计过程，MATLAB提供了直接设计模拟
44、滤波器函数。同一类型的模拟低通、 高通、 带通、 带阻滤波器的MATLAB设计是同一函数，只是参数不同而已。设计时，只要选择不同参数，就可以设计出不同类型的模拟滤波器。以巴特沃兹滤波器为例，其设计函数为 b，a = butter(N，Wn，ftype，s)当ftype选项选择high时，设计的是高通滤波器；选择 low 时，设计的是低通滤波器；选择 stop 时，设计的是带阻滤波器；缺省时，设计的是带通滤波器。必须注意的是：设计带通、 带阻滤波器时，Wn为2元数组，设计结果为2*N阶。第4章 无限长单位脉冲响应数字滤波器设计切比雪夫型滤波器：b，a = cheby1(N，Rp，Wn，ftype
45、，s)，Rp为通带最大波动(dB)。切比雪夫型滤波器：b，a = cheby2(N，Rr，Wn，ftype，s)，Rr为阻带最小波动(dB)。椭圆滤波器：b，a = ellip(N，Rp，Rr，Wn，ftype，s)，Rp为通带最大波动(dB)，Rr为阻带最小波动(dB)。【例【例4.2.4】 3阶巴特沃斯型的归一化低通滤波器的传递函数为第4章 无限长单位脉冲响应数字滤波器设计试设计下列模拟滤波器：(1) 通带为10 Hz的低通滤波器；(2) 通带下边频为10 Hz的高通滤波器；(3) 中心频率为10 Hz，带宽为2 Hz的带通滤波器；(4) 中心频率为10 Hz，带宽为2 Hz的带阻滤波器；
46、又设采样周期为0.01秒，求相应的数字滤波器。解解 MATLAB程序如下：b=1；a=1 2 2 1；Wn=10*2*pi；B=2*2*pi；b1，a1=lp2lp(b，a，Wn)b2，a2=lp2hp(b，a，Wn)b3，a3=lp2bp(b，a，Wn，B)b4，a4=lp2bs(b，a，Wn，B)第4章 无限长单位脉冲响应数字滤波器设计运行结果：b1 = 2.4805e+005a1 = 1 125.66 7895.7 2.4805e+005b2 = 1 1.2095e-014 -4.2767e-013 -2.7202e-020a2 = 1 125.66 7895.7 2.4805e+005
47、b3 = 1984.4 -5.8081e-012 -7.7026e-009 -2.9869e-013a3 = 1 25.133 12159 2.0042e+005 4.8003e+007 3.9171e+008 6.1529e+010b4 = 1 3.5545e-013 11844 2.8096e-009 4.6756e+007 5.5449e-006 6.1529e+010a4 = 1 25.133 12159 2.0042e+005 4.8003e+007 3.9171e+008 6.1529e+010第4章 无限长单位脉冲响应数字滤波器设计所以，4个模拟滤波器的传输函数依次为：第4章 无
48、限长单位脉冲响应数字滤波器设计4.3 利用模拟滤波器设计利用模拟滤波器设计IIR数字滤波器数字滤波器利用模拟滤波器设计数字滤波器，实质是实现以下映射：将Ha(s)变换为H(z)，也即从S平面映射到Z平面。这种映射变换必须遵循两个基本原则：(1) H(z)的频响要能模仿Ha(s)的频响，即S平面的虚轴应映射到Z平面的单位圆上。(2) Ha(s)的因果稳定性映射成H(z)后保持不变，即S平面的左半平面 Res0 应映射到Z平面的单位圆以内，即|z|1。下面讨论两种常用的映射变换方法：脉冲响应不变法和双线性变换法。第4章 无限长单位脉冲响应数字滤波器设计4.3.1 脉冲响应不变法脉冲响应不变法1.
49、变换关系变换关系脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列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)。为了计算H(z)，进一步假设模拟滤波器的传递函数只有单阶极点，且分母的阶数高于分子阶数，NM，Ha(s)则可表达为部分分式形式：(4.3.1)第4章 无限长单位脉冲响应数字滤波器设计其拉氏反变换为(4.3.2)式中：u(t)为单位阶跃信号。对ha(t)采样得到数字滤波器的单位脉冲响应
50、序列为(4.3.3)对h(n)取Z变换，得到数字滤波器的传递函数为(4.3.4)(4.3.5)第4章 无限长单位脉冲响应数字滤波器设计(4.3.5)式就是脉冲响应不变法的变换结果。将它与(4.3.1)式比较可以看出：(1) S平面上的单阶极点s=si变换到Z平面上是极点：z=esiT(4.3.6) (2) Ha(s)与H(z)中部分分式所对应的系数不变。(3) 这种Ha(s)到H(z)的对应变换关系，只有将Ha(s)表达为部分分式形式才成立。第4章 无限长单位脉冲响应数字滤波器设计2. 稳定性分析稳定性分析 如果模拟滤波器是稳定的，则所有极点si都在S左半平面，即Resi0 ，那么变换后为|z
