6-2 无限长单位冲激晌应 (IIR) 数字滤波器的设计方法.ppt

6.9 设计 IIR 滤波器的频率变换法 低通滤波器设计 直接2型的算法实现 同样, 实际应用中的数字滤波器也有低通、高通、带通、带阻等类型。这样设计各类数字滤波器可以有以下两种方法 :(1) 把一个归一化原型模拟低通滤波器经模拟频带变换成所需要类型 ( 包括高通,带通,带阻与另一截止频率的低通 ) 的模拟滤波器 , 然后再通过冲激响应不变法或双线性变换法数字化为所需类型的数字滤波器 , 如图 6-23(a) 所示。实际上是把这一方法中的两步合成一步来实现 , 即把模拟归一化低通原型变换到模拟低通、高通、带通、带阻等滤波器的公式与用双线性变换得到相应数字滤波器的公式合并, 6-10 先利用模拟域频带变换法, 再利用数字化法设计数字各型滤波器 一、模拟低通滤波器变换成数字低通滤波器 首先把数字滤波器的性能要求转换为与之相应的作为 “ 样本 ” 的模拟滤波器的性能要求,根据此性能要求设计模拟滤波器 , 这可以用查表的办法 , 也可以用解析的方法。然后 通过冲激响应不变法或双线性变换法 , 将此 “ 样本 ” 模拟低通滤波器 Ha(s) 数字化为所需的数字滤波器 H(z) 。下面用实例来加以说明。 1.冲激响应不变法设计数字巴特沃思低通滤波器 给定抽样频率 fs=10KHz, 要求在频率小于1kHz 的通带内 , 幅度特性下降小于 1dB;在频率大于fst =1.5kHz 的阻带内, 衰减大于15dB 。 解 (1) 讨论 f 与ω的关系及数字域性能的公式表示。模拟频率与数字频率之间为线性关系 , ω = ΩT=2πfT=2πf/ fs, T=1/fs=10-4, T为抽样周期, 故fp=1kHz 对应于ω p=2π×1×103×1/104 =0.2π 2. 双线性变换法设计数字切贝雪夫低通滤波器 这里仍利用上一实例的指标 , 但是直接由数字域给定指标 , 即在ω≤0.2π的通带范围内幅度特性下降小于 1dB, 在0.3π ≤ω≤π的阻带范围内衰减大于15dB 。 解 (1) 首先导出指标的公式表达式 , 仍为 (6-112a )式与(6-112b) 。 (2) 由于采用双线性变换, ω和Ω之间是非线性关系, 因而将数字域指标变为模拟域指标时, 首先要将数字截止频率按 (6-62) 式预畸为模拟 “ 样本 ” 滤波器的截止频率。 重写 (6-62) 式为 前面己说过 ,c 是调节模拟频带与数字频带间对应关系的一个常数 , 我们采用使模拟频率 特性与数字频率特性在低频频率处有较确切对应关系的常数 c=2/T, 则有 利用此式 , 求出对应的模拟 “ 样本 ” 滤波器的有关频率。 这样 , 模拟系统函数可写成 利用 (6-112a) 与 (6-112b) 两式 , 将数字域指标变成模拟 “ 样本 ” 的指标 , 即将 (6-125)式的 代替两式中的 ，且将ω =0.2π, ω =0.3 π分别代人 (6-125) 式 , 为了筒 化, 令T=1，则得到对模拟滤波器的指标要求为 由此看出 , 在把数字滤波器的性能要求转化为模拟滤波器的性能要求中 , 已经考虑了频率预畸。 （3）求ε。由于我们考虑的是波纹为 的切贝雪夫等波纹模拟滤波器 , 而 就是1dB 的截止频率 , 根据 (6-85) 式可求得 （4）根据（6-88)式计算滤波器阶次N。由于 ，故 ，故 代入（6－88）式，可得 选定N=4。 （5） 计算γ 及a，b。由（6-103）式得 由（6-104a）式得 由（6-104b）式得 从而求得 (6)求左半平面的两对极点。由 (6-101a ）式及 (6-101b ）式 , 当k =1,2,3,4 时 , 可得左半平面两对极点 (7) 求模拟滤波器的系统函数 Ha(s) 。按 (6-105) 式可得 分子中的常数 0.04381 是这样求出的 : 由图 6-19 可知 , N 为偶数时 , 在 s=0 处 ( 即 Ω =0 ), 它应与 s=0 处 ( 即 Ω =0) （6－130）式的值相等，从而求得分子中的常数。或者直接用 (6-105b) 式求出 (6-130) 式分子中的常数 , 结果是一样的。 实际上 , 在求出第 (4) 步以后 , 利用波纹分贝数 及所求得的 N=