第6章-无限脉冲响应数字滤波器的设计.pptVIP

如果采用双线性变换法，边界频率的转换关系为： 　　（3） 按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。设计方法及设计步骤参考本章6.2 节。 　　（4） 用所选的转换方法，将模拟滤波器Ha(s)转换成数字低通滤波器系统函数H(z)。 　　【例6.4.2】 设计低通数字滤波器，要求频率低于0.2π rad时，容许幅度误差在1 dB以内； 在频率0.3 π到π之间的阻带衰减大于15 dB。指定模拟滤波器采用巴特沃斯低通滤波器。试用双线性变换法设计数字滤波器。 解 (1) 数字低通技术指标为 ωp=0.2π rad, ? p=1 dB ωs=0.3π rad, ? s=15 dB 　　(2) 为了计算简单，取T=1 s，预畸变校正计算相应模拟低通的技术指标为： 　　(3) 设计巴特沃斯低通模拟滤波器。根据（6.2.18）式阶数N计算如下： 取N=6。将Ωs和? s代入（6.2.20）式，求得Ωc=0.7663 rad/s。这样保证阻带技术指标满足要求，通带指标有富余。 　　根据N=6，查表6.2.1得到的归一化系统函数Ga(p)为 将p=s/Ωc代入Ga(p)，去归一化得到实际的Ha(s)为 　　(4) 用双线性变换法将Ha(s)转换成数字滤波器H(z), 即 本例设计的模拟和数字滤波器幅度特性分别如图6.4.6（a）和(b)所示。此图表明数字滤波器满足技术指标要求，且无频谱混叠。 图6.4.6 例6.4.2设计的模拟和数字滤波器幅度特性 　　为形成因果稳定的滤波器，2N个极点中只取s平面左半平面的N个极点构成Ha(s)，而右半平面的N个极点构成Ha(－s)。Ha(s)的表达式为　 设N=3，极点有6个，它们分别为： 取s平面左半平面的极点s0、s1、s2组成系统函数Ha(s)，即 巴特沃斯滤波器采用对3 dB截止频率Ωc归一化，归一化后的系统函数为 令p=s/Ωc，λ=Ω/Ωc，λ称为归一化频率， p称为归一化复变量，这样巴特沃斯滤波器的归一化低通原型系统函数为　 式中，pk=sk/Ωc, 为归一化极点，用下式表示： 这样，只要根据技术指标求出阶数N，再求出N个极点，然后得到归一化低通原型系统函数Ga(p)，如果给定Ωc，再去归一化，即将p=s/Ωc代入Ga(p)中(或sk=Ωcpk)，便得到期望设计的系统函数Ha(s)。 归一化原型系统函数Ga(p)的系数bk，k=0，1，…，N－1，以及极点pk，可以由表6.2.1得到。 表6.2.1 巴特沃斯归一化低通滤波器参数 巴特沃斯滤波器的设计实质上就是根据设计指标求阶数N和3 dB截止频率Ωc的过程。下面先介绍阶数N的确定方法。 阶数N的大小主要影响通带幅频特性的平坦程度和过渡带、阻带的幅度下降速度，它由技术指标Ωp、 ?p、Ωs和?s确定。 N可能有小数部分，应取大于或等于N的最小整数。 　　求3 dB截止频率Ωc， 请注意，如果采用(1)式确定Ωc，则通带指标刚好满足要求，阻带指标有富余； 如果采用(2)式确定Ωc，则阻带指标刚好满足要求，通带指标有富余。 　　总结，低通巴特沃斯滤波器的设计步骤如下： 　　(1) 根据技术指标Ωp、 ?p、Ωs和?s，用(6.2.18)式求出滤波器的阶数N。 　　(2) 按照(6.2.13)式，求出归一化极点pk，将pk代入(6.2.12)式，得到归一化低通原型系统函数Ga(p)。也可以根据阶数N直接查表6.2.1得到pk和Ga(p)。 (3) 将Ga(p)去归一化。将p=s/Ωc代入Ga(p)，得到实际的滤波器系统函数 　【例6.2.1】 已知通带截止频率fp=5 kHz，通带最大衰减? p=2 dB，阻带截止频率fs=12 kHz，阻带最小衰减? s=30 dB，按照以上技术指标设计巴特沃斯低通滤波器。 　　解 (1) 确定阶数N。 取N=5 (2) 按照(6.2.12)式，归一化低通原型系统函数为 由N=5直接查表得到： 极点： －0.3090±j0.9511， －0.8090±j0.5878, 　－1.0000 归一化低通原型系统函数为 式中, b0=1.0000，b1=3.2361，b2=5.2361，b3=5.2361，b4=3.2361 　　分母因式分解形式为 (3) 为将Ga(p)去归一化，先求3 dB截止频率Ωc。按照(6.2.19)式，得到： 将Ωc代入(6.2.20)式，得到: 　　此时算出的　　比题目中给的Ωs小，因此，过渡带小于指标要求。或者说，在Ωs=2π×12 krad/s时衰减大于30 dB，所以说阻带指标有富余量。 　　将p=s/Ωc代入Ga(p)中, 得到： 6.3 用脉冲