数字信号处理_lecture 25.pptVIP

数字信号处理第二十五讲 中国地质大学(北京) 地球物理与信息技术学院 电子信息工程教研室 制作 第十章 IIR数字滤波器的设计Chapter 10 Designs of IIR Digital Filter §10.4 双线性变换法10.4 Bilinear Transformation Method §10.5 IIR滤波器设计的频率转换法10.5 Frequency Transformation Method of IIR Filter Design §10.6 模拟频带转换法设计10.6 Analog Frequency Band Transform Method Design * * §10.1 常用模拟低通滤波器 §10.2 由模拟滤波器设计IIR数字滤波器 §10.3 冲激响应不变法 §10.4 双线性变换法 §10.5 IIR滤波器设计的频率转换法 §10.6 模拟频带转换法设计 §10.7 数字域频带转换法设计 §10.8 在数字域设计IIR数字滤波器 §10.9 FIR和IIR滤波器的比较 本讲的主要内容 IIR滤波器设计的频率转换法 模拟频带转换法设计 数字域频带转换法设计 映射原理 可以参考图 10-11，将式(10-38)改写成 (10-38) 将这种变换解析延拓到整个s平面内与s1平面内，并令j?＝s，j?1＝s1，那么 s平面和z平面的单值映射关系为 (10-39) (10-40) 图 10-11 双线性变换法的映射关系 引入待定常数c使式(10-38)和式(10-39)变成 将映射 代入式(10-42)，可得 (10-41) (10-42) (10-43) 常数c的选择 (1) 在低频处，当? ?0，?1?0时，有?≈?1，由式(10-41)得tan(?1T/2)≈?1T/2，所以 故 (10-44) (2) 采用数字滤波器的某一特定频率? 与模拟原型滤波器的一个特定频率?严格相对应，如截止频率?c＝?1cT，即 (10-45) 映射效果 图 10-12 双线性变换中模拟频带与数字频带之间非线性关系 (3) 双线性变换映射函 数是有理函数。 (4)避免了频率响应 混叠现象。 (1) (2) (10-46) 具体实现 (10-48) 假设模拟系统函数分解为级联的低阶子系统 通过双线性变换后，离散系统函数可表示为 (10-49) (10-50) (10-51) 假设模拟系统函数分解为并联子系统 通过双线性变换后，离散系统函数为 其中 (10-54) (10-53) (10-52) 模拟低通滤波器 归一化原型 Ωc=1 模拟低通、高通、 带通、带阻滤波器 数字低通 滤波器 数字低通、高通、 带通、带阻滤波器 2 1 3 数 字 化 数 字 化 数字频 带转换 模拟频 带转换 图 10-13 IIR数字滤波器设计的频率转换法 表 9?3 模拟低通滤波器到其它频率选择性滤波器的转换 ?1、?2分别为带阻的下、上通带截止频率 ?c 低通?带阻 ?1、?2分别为带通的下、上通带截止频率 ?c 低通?带通 ?c 低通?高通 ?c 低通? 低通 转换后的通 带截止频率 映射公式 原型3dB通带截止频率?c 转换 类型 ?c同上，?0、?1、?2分别为数字带阻的中心、下、上截止频率 带阻 ?c同上，?0、?1、?2分别为数字带通的中心、下、上截止频率 带通 ?c、?c分别为模拟低通原型和所设计数字高通的截止频率 高通 频带参数 设计参量的表达式 频率变换式 数字滤波器类型 表 10-4 根据模拟低通原型设计各类数字滤波器的频率转换式及有关参量的表达式 数字巴特沃斯低通滤波器 [例 10-5] 采用巴特沃斯低通滤波器作原型，用冲激响应不变法设计IIR数字低通滤波器，给定采样频率为fs=15kHz，通带截止频率为fp=1.5kHz，通带内的波纹小于1dB；阻带截止频率为fst＝2.25kHz，阻带衰减大于15dB。 解：(1) 首先将数字滤波器的技术指标转换为相应的模拟滤波器的技术指标。 阻带截止频率： 而Rp=1dB，As=15dB，所以 通带截止频率： 我们取 ?c =10548.1 rad/s。 (2) 用查表法 查表 10-2，当N＝6时，归一化原型模拟低通巴特沃斯滤波器的频率响应为 将s用s /?c = s/10548.1代入上式，可得截止频率为?c的模拟低通滤波器 (10-55) 另一种办法也可以将?c=10548.1 rad/s，N＝6代入到 得到z平面左半平面的三对极点为 (10-56) 利用式(10-56)，由 得到模拟滤波