第三章 IIR滤波器设计方法02.ppt

第三章 无限长单位脉冲响应(IIR)滤波器的设计方法 3.3 从模拟滤波器低通原型到各种数字滤波器的频率变化 设计高通、带通、带阻等数字滤波器时,有两种方法： ① 先设计一个相应的高通、带通或带阻模拟滤波器，然后通过脉冲响应不变法或双线性变换法转换为数字滤波器。 模拟原型?模拟高通、带通、带阻?数字高通、带通、带阻，设计方法同上面讨论的低通滤波器的设计。 ② 直接利用模拟滤波器的低通原型，通过一定的频率变换关系，一步完成各种数字滤波器的设计。 频率变换 模拟原型 数字低通、高通、带通、带阻 由模拟低通到数字高通、带通、带阻 1、模拟低通—数字高通 将双线性变换中的s用1/s代替，即 2、模拟低通— 数字带通 由模拟低通到数字高通、带通、带阻 假定数字带通滤波器的中心频率为ω0 ，带通变换就是将?=0映射到数字域的ω0 上，而将?=∞映射到高低端 及 上。即应该将s的原点映射到 ，而将 点映射到 ， 满足这个要求的双线性变换为 2、模拟低通— 数字带通 由模拟低通到数字高通、带通、带阻 当 时 因此 （带通变换关系 ） 图中 点正好映射在 上，而 映射在 ， 两端，因此满足带通变换的要求。 带通变换的频率关系 设计带通时，一般只给出上、下边带的截止频率 作为设计要求。 为了应用以上变换，首先要将上下边带参数 换算成中心频率 及模拟低通截止频率 。 为此将 代入变换关系式： 设计： 由于 在模拟低通中是一对镜象频率， 代入上面两等式，求出 又 同时也就是模拟低通的截止频率 ， 有了这两个参数就可完成全部计算 例：采样 fs=400kHz，设计一巴特沃兹带通滤波器，其3dB边界频率分别为f2=90kHz,f1=110kHz,在阻带f3=120kHz处最小衰减大于10dB。 解：确定数字频域的上下边带的角频率 求中心频率： 求模拟低通的通带截止频率 与阻带边界频率 ： 从 频率增加了约1.05倍，衰减增加了(10-3）dB，故选用二阶巴特沃兹滤波器可满足指标（查表） 归一化的系统函数： 代入 ， 代入变换公式 带通滤波器设计 3dB 10dB f2 f1 f3 H(ejω) f 3、模拟低通—数字带阻 将带通的频率关系倒置 由模拟低通到数字高通、带通、带阻 给定 例 一数字滤波器采样频率fs=1kHz，要求滤除100Hz的干扰，其3dB的边界频率为95Hz和105Hz，原型归一化低通滤波器为 解 首先确定上下边界频率 求中心频率 则带阻数字滤波器为 需要求得的系统函数为 求模拟低通的通带边界频率 3.4 从低通数字滤波器到各种数字滤波器的频率变换 前面讨论由模拟低通滤波器原型来设计各种数字滤波器的方法，该设计方法同样也可直接在数字域上进行。 这种变换是由 Hp(z)所在的 z 平面到H(z)所在的 z 平面的一个映射变换。 为便于区分变换前后两个不同的z平面,我们把变换前的z 平面定义为u平面，并将这一映射关系用一个函数表示： 数字低通滤波器原型函数Hp(z) 变换 各种数字滤波器的H(z) 于是，其的原型变换可表为： 函数G(z -1) 的特性： 1) G(z -1) 是z -1 的有理函数。 2）希望变换以后的传递函数保持稳定性不变，因此要求u的单位圆内部必须对应于z 的单位圆内部。 3）G(z -1) 必须是全通函数。 为使两个函数的频率响应满足一定的变换要求，z 平面的单位圆应映射到 u 平面的单位圆上，若以 e jθ 和 e jω 分别表示u平面和z 平面的单位圆，则由上式为 且必有|G(e -jω)|=1，其中φ(ω)是e jω 的相位函数。 即函数在单位圆上的幅度必须恒为1， |G(e -jω)|称为全通函数。 全通函数的基本特性： 任何全通函数都可以表示为： 其中ai 为极点，可为实数，也可为共轭复数，但必须在单位圆以内，即 |ai| 1，以保证变换的稳定性不变，*