第六章IIR数字滤波器设计HD--8.pptVIP

Matlab的实现 调用的函数见教材P160 ： buttap, buttord, butter 设计程序见教材P161 习题 1 9 10 14 19 上式分母可以展开成为五阶多项式，或者将共轭极点放在一起，形成因式分解形式。不如直接查表简单，由N=5，直接查表得到： 极点：-0.3090±j0.9511、-0.8090±j0.5878、 -1.0000 其中： b0=1.0000,b1=3.2361,b2=5.2361,b3=5.2361,b4=3.2361 (3) 为将Ha(p)去归一化，先求3dB截止频率Ωc。 将p=s/Ωc代入Ha(p)中得到： 检验： 可以看出，满足?s=30dB 的真实fs在10.525kHz处，与12kHz比，还有富裕量。 第三节 用脉冲响应不变法 设计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(t)的采样信号用 表示， 二、S平面和Z平面之间的映射关系 拉氏变换： 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式表示： ? 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 多值映射 将s=jΩ代入上式，得 我们知道模拟信号ha(t)的傅里叶变换Ha(jΩ)和其采样信号 的傅里叶变换 之间的关系满足： 存在的问题：混叠失真 数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。 只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，才能使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不产生混叠失真。即： |ω|π 但是，任何一个实际的模拟滤波器频率响应都不是严格限带的（非理想）， 变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而是有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。 脉冲响应不变法的频率混叠现象 优点：