[信息与通信]数字信号处理第二版第五章.ppt

例5-5 用双线性变换法设计一个三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz（即采样周期为T=250 μs），其3dB截止频率为fc=1 kHz。 三阶模拟巴特沃思滤波器为 解 1、确定数字域截止频率ωc=2πfcT=0.5π。 2、确定预畸变的模拟滤波器的截止频率： 3、将Ωc代入三阶模拟巴特沃思滤波器Ha(s)，得： 4、将双线性变换关系代入就得到数字滤波器的系统函数： 图 5-16 用双线性变换法设计得到的三阶巴特沃思数字低通滤波器的频响 5.6 设计IIR滤波器的频率变换法 表5-3 截止频率为Ωc的模拟低通滤波器到其它频率选择性滤波器的转换公式 然后，通过脉冲响应不变法或双线性变换法，将此“样本”模拟低通滤波器数字化为所需的数字滤波器H(z)。 5.6.1 模拟低通滤波器变换成数字低通滤波器 首先，把数字滤波器的性能要求转换为与之相应的作为“样 本”的模拟滤波器的性能要求，根据此性能要求设计模拟滤波器 例 5-6 用脉冲响应不变法设计一个三阶巴特沃思数字低通滤波器，采样频率为fs=4 kHz（即采样周期为T=250μs），其3 dB截止频率为fc=1 kHz。 以s/Ωc代替其归一化频率，则可得三阶模拟巴特沃思低通滤波器的传递函数为： 解 查表可得归一化三阶模拟巴特沃思低通滤波器的传递函数： 式中，Ωc=2πfc。 将上式进行因式分解并表示成如下的部分分式形式： 将此部分分式系数代入(5-40)式就得到： 式中，ωc=ΩcT=2πfcT=0.5π。 将上式两项共轭复根合并，得： 将ωc=ΩcT=2πfcT=0.5π代入上式，得： 图 5-18 三阶巴特沃思数字低通滤波器的频响 由表5-3可知，由低通模拟原型到模拟高通的变换关系为： （5-52） 式中，Ωc为模拟低通滤波器的截止频率，Ωc′为实际高通滤波器的截止频率。 （5-53） 5.6.2 模拟低通滤波器变换成数字高通滤波器 模拟高通与数字高通之间的S平面与Z平面的关系仍为： 从模拟低通原型变换成数字高通滤波器公式： （5-54） 式中， 。由此得到数字高通系统函数为： 式中，Ha(s)为模拟低通滤波器传递函数。 模拟高通频率与数字高通频率之间的关系仍为： 则： 又因 故 （5-55） 模拟低通滤波器与数字高通滤波器频率之间的关系 或： （5-56） 令s=jΩ，z=ejω， 代入式（5-54），可得： 图 5-20 模拟低通变换到数字高通 例 5-7 设计一个巴特沃思高通数字滤波器，其通带截止频率（-3dB点处）为fc=3kHz，阻带上限截止频率fst=2kHz，通带衰减不大于3 dB， 阻带衰减不小于 14dB，采样频率fs=10kHz。 解 （1） 求对应的各数字域频率: （2）求常数C。采用归一化（Ωc=1）原型低通滤波器作为变换的低通原型，则低通到高通的变换中所需的C为（见表5-4） （3）求低通原型Ωst。设Ωst为满足数字高通滤波器的归一化原型模拟低通滤波器的阻带上限截止频率，可按Ω=C·cot(ω/2)的预畸变换关系来求，得： （4）求阶次N。由巴特沃思低通滤波器频率响应的公式|Ha(jΩst)|取对数， 即： 式中Ωc=1。解得： 取N=3。 （5） 求归一化巴特沃思低通原型的Ha(s)。取N=3，查表5-2可得Ha(s)为： （6） 求数字高通滤波器的系统函数H(z)， 有： 将C代入， 可求得： 5.6.3 模拟低通滤波器变换成数字带通滤波器 由表5-3可知， 由低通模拟原型到模拟高通的变换关系为 （5-57） 式中，Ωc为模拟低通滤波器的截止频率，Ωh、Ωl分别为实际带通滤波器的通带上、下截止频率。  根据双线性变换，模拟带通与数字带通之间的S平面与Z平面的关系仍为 （5-58） 把变换式（5-57）和变换式（5-58）结合起来，可得到直接从模拟低通原型变换成数字带通滤波器的表达式，也就是直接联系s与z之间的变换公式 经推导后得 式中 （5-60） （5-61） 根据双线性变换，模拟带通频率与数字带通频率之间的关系仍为 （5-62） 定义: 式中，Ω0为带通滤波器通带的中心频率，B为带通滤波器的通带宽度。 设数字带通的中心频率为ω0，数字带通滤波器的上、下边带的截止频率分别为ω2和ω1，则将式（5-62）代入式（5-63）、 式（5-64），