第6章无限脉冲响应数字滤波器的设计讲述.ppt

图6.3.3 例6.3.1的幅度特性 6.4 用双线性变换法设计IIR数字低通滤波器 (6.4.1) 脉冲响应不变法产生频率混叠的原因：模拟低通的最高截止频率超过折叠频率π/T 克服方法：采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到 之间 设Ha(s), s=jΩ, 经过非线性频率压缩后得 Ha’(s1), s1=jΩ1, 用正切变换实现频率压缩： 再通过 转换到z平面上，得到： (6.4.3) (6.4.4) (6.4.2) 则 双线性变换 图6.4.1 双线性变换法的映射关系 下面分析模拟频率Ω和数字频率ω之间的关系： 令s=jΩ,z=e jω，并代入(6.4.3)式中，有 (6.4.5) 两次映射： 从s平面映射到s1平面，具有非线性频率压缩功能，消除频率混叠现象； 从s1平面映射到z平面，其 之间的水平带的左半部分映射到z平面的单位圆内，虚轴映射到单位圆 图6.4.2 双线性变换法的频率变换关系 ω与Ω之间呈非线性关系 优点：消除频率混叠 缺点：影响DF对AF频响模仿的逼真性 非线性的实质： （1）如果Ω的刻度是均匀的，映射到z平面ω的刻度不均匀，而是随Ω的增加越来越密； （2）若AF的频响具有片断常数特性，则DF的频响也具有片断常数特性，即适合设计具有片断常数特性的滤波器 —— 一般的滤波器都适合 设 系数Ak、Bk、ak、bk之间的关系见表6.4.1。 表6.4.1 系数关系表 例6.4.1 试分别用脉冲响应不变法和双线性不变法将图6.4.4所示的RC低通滤波器转换成数字滤波器。 解： 首先按照图6.4.4写出该滤波器的传输函数Ha(s)为 利用脉冲响应不变法转换，数字滤波器的系统函 数H1(z)为 利用双线性变换法转换，数字滤波器的系统函数H2(z)为 H1(z)和H2(z)的网络结构分别如图6.4.5(a),(b)所示。 图6.4.5 例6.4.1图——H1(z)和H2(z)的网络结构 (a)H1(z); (b)H2(z) 图6.4.6例6.4.1 图——数字滤波器H1(z)和H2(z)的幅频特性 利用模拟滤波器设计IIR数字低通滤波器的步骤： (1) 确定数字低通滤波器的技术指标：通带截止频率ωp、通带衰减αp、阻带截止频率ωs、阻带衰减αs (2) 将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。 脉冲响应不变法： 双线性变换法： (3) 按照模拟低通滤波器的技术指标设计模拟低通滤波器。 (4) 将模拟滤波器Ha(s)，从s平面转换到z平面，得到数字低通滤波器系统函数H(z)。 例6.4.2 设计低通数字滤波器，要求在通带内频率低于0.2πrad时，容许幅度误差在1dB以内；在频率0.3π到π之间的阻带衰减大于15dB。指定模拟滤波器采用巴特沃斯低通滤波器。试分别用脉冲响应不变法和双线性变换法设计滤波器。 采样间隔T的选择： 脉冲响应不变法：T=1； 双线性变化法：T任选 解 (1) 用脉冲响应不变法设计数字低通滤波器。 ① 数字低通的技术指标为 ωp=0.2πrad,αp=1dB; ωs=0.3πrad,αs=15dB ② 模拟低通的技术指标为 T=1s,Ωp=0.2πrad/s,αp=1dB; Ωs=0.3πrad/s,αs=15dB ③设计巴特沃斯低通滤波器。先计算阶数N及3dB截止频率Ωc。 取N=6。为求3dB截止频率Ωc,将Ωp和αp代入(6.2.17)式，得到Ωc=0.7032rad/s，显然此值满足通带技术要求，同时给阻带衰减留一定余量，这对防止频率混叠有一定好处。 根据阶数N=6，查表6.2.1，得到归一化传输函数为 为去归一化，将p=s/Ωc代入Ha(p)中，得到实际的 传输函数Ha(s), ④用脉冲响应不变法将Ha(s)转换成H(z)。首先将 Ha(s)进行部分分式，并按照(6.3.11)式、(6.3.12)式，或 者