数字信号处理第5章 无限脉冲响应数字滤波器的设计.ppt

00-7-10 5.1 数字滤波器的基本概念 5.1.1 数字滤波器的分类 1、按处理方式分 经典数字滤波器的主要功能就是频率选择功能，利用它可以实现对输入信号按频率范围进行过滤，最大限度地保留有用信号，滤除无用信号。 现代滤波器可以按照随机信号和干扰内部的一些统计特性，从干扰中最佳地提取出信号，这方面的研究属于随机信号处理的内容。 5.2 模拟滤波器的设计 模拟滤波器模型包括巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、贝塞尔(Bessel)滤波器、椭圆(Ellipse)滤波器等等。这里仅介绍巴特沃斯滤波器的设计过程。 在设计模拟滤波器时，总是首先设计低通滤波器，再通过频率变换将低通转换成所希望类型的滤波器。 5.3 IIR数字滤波器的设计 利用模拟滤波器成熟的理论，设计IIR数字滤波器是一种人们广泛采用的设计方法。转换后H(z)稳定且满足技术要求，转换关系提出了以下两个要求： (1)因果稳定的模拟滤波器转换成数字滤波器，仍应当是因果稳定的。 (2)数字滤波器的频率响应能够模仿模拟滤波器的频率响应，也就是使得S平面的虚轴映射到Z平面的单位圆上，相应的频率之间成线性关系。 将Ha(s)进行逆拉普拉斯变换就得到模拟系统的单位冲激响应 对ha(t)进行等间隔采样，采样间隔为T，得到 对h(n)进行Z变换，就可以得到数字滤波器的系统函数，可表示为 Ha(s)的极点si映射到Z平面，其极点变为 ，系数Ai不变。 以采样信号为桥梁，讨论S平面上任意点到Z平面上的映射关系。 进行拉普拉斯变换，可得 进一步表示为 拉普拉斯变换与相应的序列h(n)的Z变换之间的映射关系可用下式表示 通常要求所设计的系统是因果稳定 上式给出确定H(z)的具体过程：首先确定模拟滤波器系统函数Ha(s)，然后在S平面上沿虚轴（s=jΩ），按照周期Ωs=2π/T对Ha(s)进行延拓，将S平面映射到Z平面，就得到了H(z)。 根据S平面和Z平面的定义，有 关系 因此得到 需要注意 jΩ轴上(-π/T，π/T)映射到单位圆上，也就是说，当模拟频率Ω在 (-π/T，π/T)内变化时，对应于数字频率的变化范围是(-π，π)，ω＝ΩT，即ω和Ω之间成线性关系。 如果ha(t)的频带宽度比π/T大，则会在±π/T的奇数倍附近产生频率混叠，映射到Z平面，表现为在位置发生频率混叠。这就是脉冲响应不变法存在的问题，即频率混叠现象， 不发生频率混叠的情况时，需要满足条件 可以得到 用脉冲响应不变法所设计的数字滤波器虽然会引起频率混叠现象，但它却能很好地保持模拟和数字滤波器之间的线性关系，即ω=ΩT，因此，在不考虑频率混叠现象的情况下，用这种设计数字滤波器方法可以很好地重现原模拟滤波器的频率特性。 5.3.2 利用双线性变换法设计IIR数字低通滤波器 为了克服频率混叠现象，提出了双线性变换法，其处理过程分为两步：第一步，将整个模拟频率范围从(-∞，∞)压缩到(-π/T，π/T)之间，这一步也被称为非线性频率压缩；第二步，利用实现从S平面到Z平面的映射。 设模拟滤波器的系统函数为Ha(s)，其中s=jΩ，利用正切变换实现非线性频率压缩，得到Ha(s1)，其中s1=jΩ1，Ω和Ω1之间的关系可以表示为 当Ω1从-π/T经过0变化到π/T时，Ω则由-∞经过0变化到∞，实现了S平面上整个虚轴完全压缩到S1平面上虚轴的(-π/T，π/T)之间的转换。 利用脉冲响应不变法进行第二次映射，将S1平面映射成Z平面，则 采用这种变换方法不可能产生频率混叠现象，如果Ha(s)因果稳定，转换成的H(z)也是因果稳定。当考虑模拟频率Ω和数字频率ω之间的关系，可以得到 双线性变换法克服了脉冲响应不变法的多值对应而产生的混叠现象。然而为此付出的代价是频率映射发生了非线性畸变，这将直接影响数字滤波器的频率响应和模拟滤波器的频率响应的相似程度，造成数字滤波器设计的失真。 为了消除数字滤波器设计的失真，可以利用某些模拟滤波器频率响应所具有的片断常数特性，进行分段处理，因为，这样转换得到的Z平面数字滤波器也具有片断常数特性。 基于上述分析，利用模拟滤波器设计IIR数字低通滤波器的步骤可以总结如下。 (1)确定数字低通滤波器的技术指标：通带截止频率ωp和衰减αp，阻带截止频率ωs和衰减αs。 (2)将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。这里主要是对边界频率ωp和ωs进行转换，对通带衰减αp和阻带衰减αs不转换。 采用脉冲响应不变法，边界频率的转换关系为：