数字信号处理第6章幻灯片.ppt

　　双线性变换法可由简单的代数公式（6.4.3）将Ha(s)直接转换成H(z)，这是该变换法的优点。但当阶数稍高时，将H(z)整理成需要的形式，也不是一件简单的工作。MATLAB 信号处理工具箱提供的几种典型的滤波器设计函数，用于设计数字滤波器时，就是采用双线性变换法。所以，只要掌握了基本设计原理，工程实际中设计就非常容易。 　　【例6.4.1】 试用脉冲响应不变法和双线性不变法将图6.4.4所示的RC低通滤波器转换成数字滤波器。 　　解 首先按照图6.4.4写出该滤波器的系统函数Ha(s)为   利用脉冲响应不变法转换，数字滤波器的系统函数H1(z)为   利用双线性变换法转换，数字滤波器的系统函数H2(z)为 图6.4.4 简单RC低通滤波器 　　设? =1000，T=0.001 s和0.002 s，H1(z)和H2(z)的归一化幅频特性分别如图6.4.5(b)和(c)所示。图6.4.5（a）是模拟滤波器幅频特性，是一个低通滤波器，但由于阶数低，选择性差，拖了很长的尾巴。图6.4.5（b）是采用脉冲响应不变法转换成的数字滤波幅频特性曲线，图中ω=π处对应的模拟频率与采样间隔T有关，当T=0.001 s时，对应的模拟频率为1/(2T)=500 Hz； 当T=0.002 s时，对应的模拟频率为250 Hz。对照图 6.4.5（a）， 图(b)所示的数字滤波器与原模拟滤波器的幅度特性差别很大，且频率愈高，差别愈大。 　　这是由频率混叠失真引起的。相对而言，T=0.001 s时混叠少一些。图6.4.5（c）是采用双线性变换法转换成的数字滤波器幅频特性曲线，由于该转换法的频率压缩作用，使 ω=π处的幅度降为零。但曲线的形状与原模拟滤波器幅度特性曲线的形状差别较大，这是由于该转换法的非线性造成的，T小一些，非线性的影响少一些。 图6.4.5 例6.4.1中Ha(s)、H1(z)和H2(z)的幅频特性 　　下面我们总结利用模拟滤波器设计IIR数字低通滤波器的步骤。 　　（1） 确定数字低通滤波器的技术指标： 通带边界频率ωp、通带最大衰减? p、阻带截止频率ωs、阻带最小衰减? s。 　　（2） 将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术指标。这里主要是边界频率ωp和ωs的转换， ? p和? s指标不变。如果采用脉冲响应不变法，边界频率的转换关系为 （6.4.6） 如果采用双线性变换法，边界频率的转换关系为   　　（3） 按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。设计方法及设计步骤参考本章6.2 节。 　　（4） 用所选的转换方法，将模拟滤波器Ha(s)转换成数字低通滤波器系统函数H(z)。 （6.4.7） 　　在设计过程中，要用到采样间隔T，下面介绍T的选择原则。如采用脉冲响应不变法，为避免产生频率混叠现象，要求所设计的模拟低通带限于－π/T~π/T区间。由于实际滤波器都是有限阶的，因此有一定宽度的过渡带，且频响特性不是带限于π/T。当给定模拟滤波器Ha(s)，要求单向转换成数字滤波器H(z), 且? s足够大时，选择T满足|Ωs|π/T，可使频谱混叠足够小，满足数字滤波器指标要求。 　　但如果先给定数字低通的技术指标时，情况则不一样。由于数字滤波器系统函数H(ejω)以2π为周期，最高频率在ω=π处，因此，ωsπ，按照线性关系Ωs=ωs/T，那么一定满足Ωsπ/T，这样T可以任选。一般选T=1 s。 这时，频谱混叠程度完全取决于? s， ? s越大，混叠越小。对双线性变换法不存在频谱混叠现象，尤其对于设计片断常数滤波器，T也可以任选。为了简化计算，一般取T=2 s。 　　【例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), 即