第3章：无限长单位脉冲响应(IIR)滤波器设计-第一节.ppt

1、第三章 基本内容,设计数字滤波器的一般步骤：性能指标、逼近、实现（结构），理解 模拟低通滤波器特性：三种，理解 根据模拟滤波器设计数字滤波器的两种方法：脉冲响应不变法和双线性变换法，灵活应用 从模拟低通原型到各种数字滤波器的频率变换：4种，s-z，掌握 从低通数字滤波器到各种数字滤波器的频率变换：z平面变换法，掌握 IIR最优化设计，了解,第三章 无限长单位脉冲响应（IIR）滤波器设计,概述: 许多信息处理过程，如信号的过滤，检测、预测等都要用到滤波器，数字滤波器是数字信号处理中应用最广泛的一种线性系统，是数字信号处理的重要基础。 数字滤波器的功能（本质）是将一组输入的数字序列通过一定的运算后

2、转变为另一组输出的数字序列。实现方法主要有两种：数字信号处理硬件和计算机软件。,数字滤波器的数学描述： 1)差分方程,2)系统函数,（3）单位脉冲响应 递归系统 IIR 非递归系统 FIR （4）频率响应函数：注意数字频率以2为周期 高通 低通 带通 带阻,滤波器的幅频特性是分段常数的，以低通为例，在通带内 逼近于1，在阻带内逼近于0.,确定系数 、 或零极点 、 ，以使滤波器满足给定的性能要求第三章、四章讨论 3）用一个有限精度的运算去实现这个系统函数。包括 选择运算结构：如级联型、并联型、卷积型、频率采样型以及快速卷积（FFT）型等；选择合适的字长和有效数字的处理方法等,数字滤波器的设计步

3、骤： 1）按照实际需要确定滤波器的性能要求。 2）用一个因果稳定系统的 H(z) 或 h(n) 去逼近这个性能要求，即求 h(n) 的表达式。,IIR设计方法： 1）先设计一个合适的模拟滤波器，然后变换成满足预定指标的数字滤波器。 由于模拟的网络综合理论已经发展得很成熟，模拟滤波器有简单而严格的设计公式，设计起来方便、准确、可将这些理论推广到数字域，作为设计数字滤波器的工具。,此外还有其他多种误差最小准则， b)在此最佳准则下，求滤波的系数 和 通过不断地迭代运算，改变 、 ，直到 满足要求为止。,2）最优化设计方法 分两步： a)确定一种最优准则，如最小均方误差准则，即使设计出的实际频率响应

4、的幅度特性 与所要求的理想频率响应 的均方误差最小，,以上两种设计方法中，着重讲第一种，因为数字滤波器在很多场合所要完成的任务与模拟滤波器相同，如作低通、高通、带通及带阻网络等，这时数字滤波也可看作是“模仿”模拟滤波器。在IIR滤波器设计中，采用这种设计方法目前最普遍。由于计算机技术的发展，最优化设计方法的使用也逐渐增多。,3.1 根据模拟滤波器设计IIR滤波器,利用模拟滤波器设计数字滤波器，就是从已知的模拟滤波器传递函数Ha(s)设计数字滤波器传递函数H(z)，这归根到底是一个由s平面到z平面的变换，这种映射变换应遵循两个基本原则： 1）频响不变：H(z)的频响要能模仿Ha(s)的频响，即s

5、平面的虚轴应映射到z平面的单位圆 上。 2）因果稳定性不变：Ha(s) 的因果稳定性映射成 H(z)后保持不变，即S平面的左半平面 ReS0 应映射到Z平面的单位圆以内|Z|1。,利用模拟滤波器理论设计数字滤波器，也就是使数字滤波器能模仿模拟滤波的特性，这种模仿可从不同的角度出发。 学习：基本思路、转换步骤、s平面与z平面的映射关系、结论、Matlab实现 下面讨论两种常用的映射变换方法。 一、脉冲响应不变法 脉冲响应不变法基本理论 映射关系 混叠现象,一、脉冲响应不变法 （1）基本思路 脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)正好等于模拟滤波器的冲激响应

6、ha(t)的采样值，即 h(n)=ha(nT), T为采样周期。 （2）设计步骤 如以Ha(s)及H(z)分别表示ha(t)的拉氏变换及 h(n)的z变换，即 Ha(s)=Lha(t) , H(z)=Zh(n),设模拟滤波器的系统函数只有单阶极点，且分母的阶数高于分子阶数 NM，则可表达为部分分式形式； 其拉氏反变换为： u(t)为单位阶跃函数 对ha(t)采样得到数字滤波器的单位脉冲响应序列,对h(n)取Z变换，得到数字滤波器的传递函数： 第二个求和为等比级数之和，要收敛的话， 必有 即得到数字滤波器的系统函数：,1） 理想采样 的拉氏变换 与采样序列 的 z 变换 之间存在的 s 平面与

7、Z 平面的映射关系。,（3）S平面与z平面的映射关系,s平面与z平面的映射关系,根据理想采样的频谱特性，不难证明离散序列拉氏变换与模拟信号拉氏变换的关系：js,2） 理想采样 的拉氏变换 与模拟信号 的拉氏变换 之间的关系。,以上表明，采用脉冲响应不变法将模拟滤波器变换为数字滤波器时，它所完成的 s 平面到 z 平面的变换，正是拉氏变换到z变换的标准变换关系，即首先对Ha(s)作周期延拓，然后再经过 的映射关系映射到 z 平面上。,T采样周期,3）特性 稳定性 如果模拟滤波器是稳定的，则所有极点 Si 都在S左半平 面，即 Resi0 , 那么变换后H(Z)的极点 也都在单位圆以内，即 因此数

8、字滤波器保持稳定。,由此可见： 当=0时，r=1，这表明s平面虚轴映射为z平面的单位圆。 当0时，r1.这表明s左半平面映射为z平面的单位圆内部，而s平面右半平面映射为z平面单位圆外部。 由于=T，所以当自0至变化时，的对应值为0至/T。这表明，s平面上每一条宽为2 /T的横带部分，都将重叠地映射到z平面的整个平面上；而每一横带的左半部分映射到z平面单位圆以内，每一横带的右半部分映射到z平面单位圆以外。虽然j轴映射到单位圆上，但是j轴上的每一段2 /T都对应于绕单位圆一周。,映射关系 ：,S 平面,Z 平面,这时数字滤波器的频响才能不失真地重现模拟滤波器的频响（存在于折叠频率S/2以内），即

9、但任何一个实际的模拟滤波器，其频响都不可能是真正带限的，因此不可避免地存在频谱的交叠，即混淆，如图，这时，数字滤波器的频响将不同于原模拟滤波器的频响而带有一定的失真。模拟滤波器频响在折叠频率以上衰减越大，失真则越小，这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。,正如第一章的采样定律中所讨论的，如果模拟滤波器的频响带限于折叠频率S/2 以内，即,混叠 数字滤波器的频响并不是简单的重现模拟滤波器的频响，而是模拟滤波器频响的周期延拓：,脉冲响应不变法中的频响混淆,（4）几个 结论, 1)脉冲响应不变法的一个重要特点是频率坐标的变换是线性的，与是线性关系。 因此如果模拟滤波器的频响带限

10、于折叠频率以内的话，通过变换后数字滤波器的频响可不失真地反映原响应与频率的关系。 例如线性相位的贝塞尔低通滤波器，通过脉冲响应不变法得到的仍是线性相位的低通数字滤波器。 2)在某些场合，要求数字滤波器在时域上能模仿模拟滤波器的功能时，如要实现时域冲激响应的模仿，一般使用脉冲响应不变法。,3)如果Ha(s)是稳定的，即其极点在S左半平面，映射后得到的H(Z)也是稳定的。 4)脉冲响应不变法的最大缺点：存在周期延拓而导致的频谱混叠效应，因此只能用于带限的频响特性，例如衰减特性很好的低通或带通，而高频衰减越大，频响的混淆效应越小。至于高通和带阻滤波器,由于它们在高频部分不衰减，因此不宜采用脉冲响应不

11、变法。,（5）Matlab实现,bz,az=impinvar(1,1,1,5,6,10) 10为采样频率 数字滤波器的分子分母系数,将一个具有如下系统函数 的模拟滤波器数字化。 解：,例1,模拟滤波器的频率响应为:,数字滤波器的频率响应为: 显然 与采样间隔T有关,如图b, T越小,衰减越大,混叠越小,当 fs=24Hz ,混叠可忽略不计 如何与T无关？ 看3.15式？ 改造h(n)=Tha(nT) 。,二、双线性变换法 （1）基本思路 脉冲响应不变法是使数字滤波器在时域上模仿模拟滤波器，主要缺点是频谱交叠产生的混淆，这是从S平面到Z平面的标准变换zesT的多值对应关系导致的,为了克服这一缺点

12、，设想变换分为两步： 第一步：将整个S平面压缩到S1平面的一条横带里； 第二步：通过标准变换关系将此横带变换到整个Z平面上去。 由此建立S平面与Z平面一一对应的单值关系，消除多值性，也就消除了混淆现象。,这里c是待定常数， c使模拟滤波器的频率特性与数字源波器的频率特性在不同频率点有对应关系。 经过这样的频率变换， 当由 时, 1由-/T经过变化到/T ，即S平面的整个j轴被压缩到S1平面的2/T 一段。,虚轴频率变换：为了将S平面的j轴压缩到S1平面j1轴上的-/T到/T一段上，可通过以下的正切变换实现：,（2）转换步骤,通常取c=2/T,再将 S1 平面通过标准变换关系映射到Z平面，即令,

13、将这一关系解析扩展至整个S平面 ，则得到S平面到S1平面的映射关系：,最后得S平面与Z平面的单值映射关系：,(1)式可以写为：,把变换式整理一下变成一个代数方程，可以看到方程以s为变量时是一个线性方程， 以z -1为变量时也是线性方程，因此由这个方程确定的变换称为双线性变换,现在我们看看，这一变换是否符合我们一开始提出的由模拟滤波器设计数字滤波器时，从 S平面到Z平面映射变换的二个基本要求： 单位圆映射：当 时，得：,时, 稳定性说明,即s左半平面映射在单位圆内，s右半平面映射在单位圆外，因此稳定的模拟滤波器通过双线性变换后，所得到的数字滤波器也是稳定的。,与脉冲响应不变法相比，双线性变换的主

14、要优点：S平 面与Z平面是单值的一一对应关系即整个j轴单值的对应于单位圆一周，避免了频谱混叠。 付出的代价是在频率轴上引进了失真，或者说靠频率的严重非线性关系得到频率的单值对应关系，因此，只有当容忍或能补偿这种失真时，使用双线性变换设计数字滤波器才是实用的。,考虑z = ej,频率映射关系,这是两个域关于频率的变换关系,图 双线性变换的频率非线性关系 由图中看到，在零频率附近，接近于线性关系，进一步增加时，增长变得缓慢， (终止于折叠频率处)，所以双线性变换不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象。,虽然双线性变换有这样的缺点，但它目前仍是使用得最普遍、最有成效的一种设计工具。这是因为大多数滤波器都具有分段常数的频响特性，如低通、高通、带通和带阻等，它们在通带内要求逼近一个衰减为零的常数特性，在阻带部分要求逼近一个衰减为的常数特性，这种特性的滤波器通过双线性变换后，虽然频率发生了非线性变化，但其幅频特性仍保持分段常数的特性。 ！各个分段边缘的临界频率点的畸变，可以校正。,预畸变：,即将模拟滤波器的临界频率事先加以畸变，然后通过双线性变换后正好映射到所需要的频率上。 利用关系式： 将所要设计的数字滤波器临界频率点 ，变换成对应的模拟域频率 ，利用此 设计模拟滤波器，再通过双线性变换，即可得到所需的数