西安电子科技大学(高西全丁美玉第三版)数字信号处理第6章.pptVIP

第5章　无限脉冲响应(IIR)数字滤波器的设计 　　3． IIR数字滤波器的设计方法　　关于滤波器的设计原理与具体的设计方法， 教材及其他《数字信号处理》书中都有详细叙述， 本章不再赘述。 下面仅对IIR数字滤波器的几种设计方法及设计步骤作简要归纳， 并指出学习要点， 并通过例题及习题与上机题解答说明各种设计方法的具体设计过程及相关设计公式， 以便读者能有条理地解答滤波器设计题目， 设计实际应用滤波器。 　　为了使初学者对IIR数字滤波器设计方法有一个整体概念， 先抛开繁杂的设计过程和设计公式， 用图5.1.3归纳IIR数字滤波器的一般设计方法。 　　下面对图5.1.3中给出的五种设计方法及其学习要点进行简要归纳。 对频域直接逼近法和时域波形逼近法都必须借助计算机设计（即CAD设计）， 且已有商业设计程序， 所以只简要介绍其设计思想及逼近准则。 5.1.2 模拟滤波器的设计　　为了叙述方便， 用AF表示模拟滤波器， 用DF表示数字滤波器。 从教材中的详细介绍知道， 从AF入手设计DF时， 首先要设计一个“相应的AF”， 所以下面以流程图形式给出AF的设计步骤， 如图5.1.4所示。 　　由于AF设计手册中给出了各种典型AF归一化低通原型的设计公式和图表及系统函数G(p)， 因此设计AF很方便。 所以， 对需要设计的实际AF（低通、 高通、 带通和带阻AF)，首先将其指标参数转换成相应的归一化(λp=1)低通指标参数， 将设计各种实际AF转化为设计归一化低通AF； 最后将设计好的归一化低通G(p)转换成实际滤波器H(s)。 为此，下面归纳四种实际AF系统函数H(s)与其相应的归一化低通AF系统函数G(p)的相互转换关系， 如图5.1.5所示。 图中总结出了图5.1.4中第（2）步和第（4）步所涉及的所有转换关系和转换公式。 图5.1.6 图5.2.1 图5.2.2 　　在实际工作中， 调用IIR滤波器阶数来计算程序函数， 很容易求出满足要求的最小阶数N值。 如前述， buttord函数用于计算butterworth滤波器阶数。 本例中， 求阶数的MATLAB程序如下： 　　　Wp=800； Rp=3； 　　　Ws=1931.37； Rs=20； 　　　［N、 Wc］=buttord(Wp， Ws， Rp， Rs， ′s′)； 　　运行结果： 　　　N=3， Wc=897.9654 　　如上结果中，Ωc的值未代入相乘， 这样使读者能清楚地看到去归一化后，3 dB截止频率对归一化系统函数的改变作用。 　　  2． 设计一个切比雪夫低通滤波器， 要求通带截止频率fp=3 kHz，通带最大衰减αp=0.2 dB，阻带截止频率fs=12 kHz, 阻带最小衰减αs=50 dB。 求出滤波器归一化系统函数G(p)和实际的Ha(s)。 　　解: （1） 确定滤波器技术指标。 　　　αp=0.2 dB, Ωp=2πfp=6π×103 rad/s　　　αs=50 dB,　Ωs=2πfs=24π×103 rad/s　　　 λp=1, 　　解: 该题所给Ha(s)正是模拟滤波器二阶基本节的两种典型形式。 所以， 求解该题具有代表性， 解该题的过程， 就是导出这两种典型形式的Ha(s)的脉冲响应不变法转换公式。 设采样周期为T。 题6解图 题8解图 　　11． 设计一个数字高通滤波器， 要求通带截止频率ωp=0.8π rad， 通带衰减不大于3 dB，阻带截止频率ωs=0.5π rad， 阻带衰减不小于18 dB。 希望采用巴特沃斯型滤波器。 　　解: （1） 确定数字高通滤波器技术指标： 　　　　　ωp=0.8π rad, αp=3 dB　　　　　ωs=0.5π rad, αs=18 dB （2） 确定相应模拟高通滤波器技术指标。 由于设计的是高通数字滤波器， 所以应选用双线性变换法， 因此进行预畸变校正求模拟高通边界频率（假定采样间隔T=2 s）： 　　12． 设计一个数字带通滤波器， 通带范围为0.25 π rad到0.45π rad， 通带内最大衰减为3 dB， 0.15 π rad以下和0.55 π rad以上为阻带， 阻带内最小衰减为15 dB。 要求采用巴特沃斯型模拟低通滤波器。 　　解: （1） 确定数字带通滤波器技术指标： 　　　　ωpl=0.25π rad, ωpu=0.45π rad　　　　ωsl=0.15π rad, ωsu=0.55π rad　　通带内最大衰减αp=3 dB， 阻带内最小衰减αs=15 dB。