用脉冲响应不变法设计IIR数字低通滤波器.ppt

* * 第三节 用脉冲响应不变法 设计IIR数字低通滤波器 利用模拟滤波器来设计数字滤波器，就是从已知的模拟滤波器传递函数Ha(s)设计数字滤波器的系统函数H(z)。 因此，它归根结底是一个由S平面映射到Z平面的变换，这个变换通常是复变函数的映射变换，为了保证转换后的H(z)稳定且满足技术要求，这个映射变换必须满足以下两条基本要求： (2)数字滤波器的频率响应模仿模拟滤波器的频 响，s平面的虚轴映射为z平面的单位圆，相应的频率之间成线性关系。 因果稳定的模拟滤波器转换成数字滤波器，仍 是因果稳定的。即 S平面的左半平面必须映射到Z平面单位圆的内部。 一、脉冲响应不变法的转换原理 核心原理： 通过对连续函数ha(t)等间隔采样得到离散序列ha(nT)。令h(n)= ha(nT) ，T为采样间隔。它是一种时域上的转换方法。 转换步骤： 设模拟滤波器的传输函数为Ha(s),相应的单位冲激响应是ha(t) 设模拟滤波器Ha(s)只有单阶极点，且分母多项式的阶次高于分子多项式的阶次，将Ha(s)用部分分式表示： 式中si为Ha(s)的单阶极点。将Ha(s)进行逆拉氏变换得到ha(t)： 式中u(t)是单位阶跃函数。对ha(t)进行等间隔采样，采样间隔为T，得到： 对上式进行Z变换，得到数字滤波器的系统函数H(z)： 一般Ha(s)的极点si是一个复数，且以共轭成对的形式出现，将一对复数共轭极点放在一起，可形成一个二阶基本节。如果模拟滤波器的二阶基本节的形式为： 可以推导出相应的数字滤波器二阶基本节(只有实数乘法)的形式为： 极点为 如果模拟滤波器二阶基本节的形式为： 极点为 相应的数字滤波器二阶基本节的形式为： 设ha(t)的采样信号用 表示， 二、存在的问题 拉氏变换： 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示： 将s=jΩ代入上式，得 我们知道模拟信号ha(t)的傅里叶变换Ha(jΩ)和其采样信号 的傅里叶变换 之间的关系满足： ? j ? S平面 Z平面 Re[z] jIm[z] ?=0 (s平面虚轴) r=1 (z平面单位圆)?0 (s左半平面) r1 (z平面单位圆内部)?0 (s右半平面) r1 (z平面单位圆外部) ① r与?的关系 ② ?与?的关系： ?= ?T ?=0 (s平面实轴) ?=0 (z平面正实轴) ?= ?0 ?=?0T  (s平面平行于实轴的直线) (z平面始于原点，辐角为?0T的辐射线) ?：从 -?/T ～ ?/T ?：从 -? ～ ?  (s平面为2?/T的一个水平带) (z平面辐角转了一周，覆盖整个z平面) Z平面 jIm[z] Re[z] ? j ? S平面 -?/T ?/T 3?/T -3?/T 多值映射 存在的问题：混叠失真 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。 只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不产生混叠失真。即： |ω|π 但是，任何一个实际的模拟滤波器频率响应都不是严格限带的（非理想）， 变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而是有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。 脉冲响应不变法的频率混叠现象 数字滤波器频率响应幅度还与采样间隔T成反比： |ω|π 如果采样频率很高，即T很小，数字滤波器可能具有太高的增益，这是不希望的。为了使数字滤波器增益不随采样频率而变化，可以作以下简单的修正，令： h(n)=Tha(nT) 则有: 很小的数在分子 比在分母安全！ 例： 已知模拟滤波器的传输函数Ha(s)为： 用脉冲响应不变法将Ha(s)转换成数字滤波器的系统函数H(z)。 解：首先将Ha(s)写成部分分式 极点为： 那么H(z