门爱东的数字信号处理教案4.ppt

第四章 IIR数字滤波器的设计 4.1 引言 4.2 模拟滤波特性的逼近 4.3 模拟系统的数字仿真 4.4 脉冲响应不变法 4.5 双线性变换法 4.6 数字滤波器的变换 4.7 IIR滤波器的优化设计 4.1 引言 1.数字滤波器的设计过程 一个数字滤波器的设计过程，大致可归纳为3个步骤。 (1) 按照实际需要确定滤波器的性能要求。 (2) 用一个因果稳定的系统函数去逼近这个性能要求，这种系统函数可以分为两类，即IIR系统函数与FIR系统函数。 (3) 用一个有限精度的运算去实现这个系统函数。这里包括选择运算结构，如级联型、并联型、直接型、正准型和频率采样型等，还包括合适的字长和有效数字处理的方法等。 2. IIR数字滤波器的设计方法 IIR数字滤波器的系统函数是z的有理函数，可表示为 这种设计一般有3种方法。 (1) 简单滤波器的零极点位置累试法。 (2) 用模拟滤波器的理论来设计数字滤波器。 (3) 用优化技术设计参数。 3.设计原则 利用模拟滤波器来设计数字滤波器，即从已知的模拟滤波器系统函数Ha(S)求数字滤波器的系统函数H(z)。其常用的方法有脉冲响应不变法和双线性变换法两种，但无论哪种方法都遵循数字仿真原则。 4.2 模拟滤波特性的逼近 1.问题的提出 其频带按作用不同分为以下3部分，如图4-1所示： 通带——滤波器中能使信号通过的频带； 阻带——抑制噪声通过的频带； 过渡带——从通带到阻带的过渡频率范围。 按划分的标准不同，滤波器分为很多种类。按处理信号的种类分为模拟滤波器(其输入输出皆为模拟量)和数字滤波器(其输入输出皆为数字量)；按其应用不同分为有源滤波器和无源滤波器；按其通频带划分为低通滤波器(只允许低频信号通过而抑制高频信号)、高通滤波器(只允许高频信号通过而抑制低频信号)、带通滤波器(只允许某一频带的信号通过)和带阻滤波器(只抑制某一频带的信号)。 2.巴特沃思(Butterworth)滤波特性的逼近 (1) Butterworth滤波器特性 Butterworth滤波器的特点是具有通带内最大平坦的振幅特性，而且在正频率范围内是随频率升高而单调下降的，它的振幅平方函数为 (2) Butterworth图表设计法 其设计步骤概括起来有以下几个方面： (3)巴特沃斯低通滤波器的 Matlab 实现 在 Matlab 中可以用 [n,Wn] = buttord(Wp,Ws,Rp,Rs,‘s’) 求出给定通带截止频率 wp、阻带起始频率 ws、通带最大衰减Rp、阻带最大衰减Rs 条件下，所需要的最小 Butterworth 模拟滤波器阶数和 3dB 频率。 Buttord(Wp,Ws,Rp,Rs,’s’) 也可用于高通、带通、阻带滤波器的设计。对于带通和带阻滤波器，Wp 和 Ws 说明了滤波器的两个边缘频率，其中第一个元素是低端的频率。 3.切比雪夫(Chebyshev) 滤波特性的逼近 (1) Chebyshev滤波器特性 Chebyshev滤波器频率特性在通带内是等波纹的，在阻带内是单调下降的，它比同阶的Butterworth滤波器下降快，如图4-5所示，而在通带内则是等波纹的，通带内波动的峰值个数与阶数N有关，N越大，波动越快 Chebyshev滤波器的幅度平方函数为 式中N是滤波器的阶数，Ωc是通带宽度，当ε=1时，它等于通带的3dB截止频率，如图4.5 所示。 (2) Chebyshev滤波器特性的极点 对于N阶Chebyshev滤波器，其幅度｜Ha(Ω)｜2平方函数在s平面上有2N个极点，这2N个极点都是成复共轭对出现的，且关于虚轴对称，没有点落在虚轴上。 极点可由作图法确定，即 (1) 以aΩc为半径画圆，将其圆周等分为2N个点； (2) 以bΩc为半径画圆，将其圆周等分为2N个点； (3) 每个极点的纵坐标由大圆的相应点的纵坐标确定，每个极点的横坐标由小圆的对应点的横坐标确定。 (3) Chebyshev滤波器的设计方法 如同Butterworth滤波器，Chebyshev滤波器也可以采用查表法进行设计，其步骤如下。 (1) 频率 归一化，求λp，λs。 (2) 根据αp，αs查归一化频率——幅度特性曲线，确定滤波器的阶数N。 (4) Chebyshev滤波器的 Matlab 实现 在 Matlab 中可以用 [n,Wn] = cheb1ord(Wp,Ws,Rp,Rs,s) 求出给定通带截止频率 wp、阻带起始频率 ws、通带最大衰减Rp、阻带最大衰减Rs 条件下，所需要的最小Chebyshev I 模拟滤波器阶数。在数字滤波器中，wp和ws 在0~1之间，而在模拟滤波器中，wp和ws可以大于1。 Che