5-1IIR.ppt

信号分析与处理;第5章 数字滤波器;5.1 基本概念; 数字滤波器相对模拟滤波器而言，在体积、重量、精度、稳定性、可靠性、存储功能、灵活性以及性能价格比等方面都显示明显的优点，而且，数字滤波器除利用硬件电路实现之外还可借助计算机以软件编程方式实现。; 数字滤波器是一类重要的离散时间系统。可以从时域、复数域和频域的角度来描述。; 如果x(n)，y(n)和h(n)满足绝对可积条件，则相应地X(z)、Y(z)和H(z)均在单位圆上收敛，离散时间系统频率响应为;② 若1 ? k ? M，ak = 0，有;解：由差分方程可得;?H(e j? )? = sin(0.5?) ?(?) = ?/2? 0.5?; 【例5-2】已知一数字滤波器用差分方程表示为 y(n) ? y(n?1) = 0.5x(n) 试求其单位抽样响应h(n)。 解：由差分方程可得; 设计数字滤波器的第一步仍然是寻求满足性能要求的系统函数H(z)，因此这也是一逼近问题。;?; IIR和FIR数字滤波器的设计思路和方法有所不同。 IIR数字滤波器是对模拟滤波器的模仿，满足滤波幅频特性的要求。 而FIR数字滤波器是直接用理想数字滤波的单位抽样响应hd(n)的逼近来设计，能够同时保证信号传输中幅度和相位的要求。 ;5.2 IIR数字滤波器设计; 设计一般有以下两种方法： （1）计算机辅助设计法。 这是一种最优化设计法。先确定一种最优准则，例如设计出的实际频率响应幅度? H(e j?)?与所要求的理想频率响应幅度? Hd(e j? )?的均方误差最小准则等等，然后求出在此最优准则下滤波器系统函数的系数br和ak。这种设计一般得不到滤波器系数作为所要求的理想频率响应的闭合形式的函数表达式，而是需要进行大量的迭代运算，故离不开计算机。; （2）借助模拟原型滤波器导出所需数字滤波器。此方法由于模拟滤波器有简便的设计公式，有大量的设计图表可资利用，设计起来既方便又准确。这种设计方法的流程如图.; 此方法实际上是s域与z域之间的映射转换。为使数字滤波器保持模拟滤波器的特性，这种映射关系应满足下列条件： 1）H(z)的频率响应要能模仿Ha(s)的频率响应，即s平面的虚轴 j?必须映射到z平面的单位圆e j?上，也就是频率轴要对应。 2）因果稳定的Ha(s) 应能映射成因果稳定的H(z)。也就是s平面的左半平面（Re[s] 0）必须映射到z平面的单位圆内部（? z ? 1）。;5.2.2 冲激响应不变法;滤波器的系统函数 H(z) = z[h(n)];2．稳定性与逼近程度 如果模拟滤波器是稳定的，其Ha(s)的所有极点sk均在s左半平面，即Re[sk ]0。极点sk映射到z平面是位于;将s = j?，z = e sT= e j?T，? = ?T代入上式，得;Ha(j?) = 0 ; 结论：冲激响应不变法可以把稳定的模拟滤波器变换成稳定的数字滤波器，变换时频率间呈线性关系? = ?T，数字滤波器频率特性形状基本上与模拟滤波器相同（如果混叠不严重），在时域两者冲激响应形状一致。 由于混叠的结果，使频率特性高端严重失真，因而只适用于限带的模拟滤波器，所以高通和带阻滤波器不宜采用冲激响应不变法，否则要加保护滤波器，滤掉高于折叠频率以上的频率。对于低通和带通滤波器，需充分限带，若阻带衰减越大，则混叠效应越小。;数字滤波器的指标;【例5-3】设一模拟滤波器的系统函数为;26;27;28;29; 【例5-4】用冲激响应不变法设计一Butterworth数字滤波器，满足： 通带允许起伏：?3dB 0 ? ? ? 0.4π 阻带衰减： ? 15dB 0.8π ? ? ? π 设抽样频率为1000Hz。;3）用冲激响应不变法求H(z)。 ;将已知的?c = 400π，T=0.001代入，得 ;5.2.3 双线性变换法; 数值近似的方法很多，例如，用差分直接代替微分，但结果使时域、频域上产生的误差大； 另一种近似方法是采用将微分方程作积分，再对积分进行数值近似（梯形近似积分）。;用梯形法逼近积分项 ;2．稳定性与逼近程度; 说明s平面的虚轴映射到z平面单位圆上，而且? = ??为? = ??，?与?呈单值映射，可以避免混叠。;ω→0，tanω/2→ω/2 ，其频率转换关系近似成直线。; 双线性变换是非线性变换，频率频率特性会引起失真。双线性变换频率轴的这一非线性畸变问题，对于大量常见的一定频段内为常数幅频特性的滤波器来说，问题并不严重，如一般