第20章--用脉冲响应不变法设计IIR数字滤波器.pptVIP

一、实验目的　　(1)加深对脉冲响应不变法设计IIR数字滤波器基本方法的了解。　　(2)掌握使用模拟滤波器原型进行脉冲响应变换的方法。　　(3)了解MATLAB有关脉冲响应变换的子函数。 二、实验涉及的MATLAB子函数　　impinvar　　功能：用脉冲响应不变法实现模拟到数字的滤波器变换。　　调用格式：　　［bd，ad］＝impinvar(b，a，Fs)；将模拟滤波器系数b、a变换成数字的滤波器系数bd、ad，两者的冲激响应不变。　　［bd，ad］＝impinvar(b，a)；采用Fs的缺省值1Hz。 三、实验原理　　1.脉冲响应不变法的基本知识　　脉冲响应不变法又称为冲激响应不变法，是将系统从s平面到z平面的一种映射方法，使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(n)。其变换关系式为z＝esT。 　　由于esT是一个周期函数，因而s平面虚轴上每一段2p/T的线段都映射到z平面单位圆上一周。由于重叠映射，因而冲激响应不变法是一种多值映射关系。数字滤波器的频率响应是原模拟滤波器的频率响应的周期延拓(如图20-1所示)。只有当模拟滤波器的频率响应是有限带宽的，且频带宽度 ，才能避免数字滤波器的频率响应发生混叠的现象。由于脉冲响应不变法只适用于限带的模拟滤波器，因此，在高频区幅频特性不等于零的高通和带阻滤波器不能采用脉冲响应不变法。 　　用MATLAB冲激响应不变法进行IIR数字滤波器设计的步骤如下(参见图19-1)：　　(1)输入给定的数字滤波器设计指标；　　(2)根据公式W＝w/T，将数字滤波器指标转换成模拟滤波器设计指标；　　(3)确定模拟滤波器的最小阶数和截止频率；　　(4)计算模拟低通原型滤波器的系统传递函数；　　(5)利用模拟域频率变换法，求解实际模拟滤波器的系统传递函数；　　(6)用脉冲响应不变法将模拟滤波器转换为数字滤波器。 　　 图20-1 数字滤波器的幅频响应 　　2.用脉冲响应不变法设计IIR数字低通滤波器　　例20-1 采用脉冲响应不变法设计一个巴特沃斯数字低通滤波器，要求：wp＝0.25p，Rp＝1 dB；ws＝0.4p，As＝15 dB，滤波器采样频率Fs＝2000 Hz。　　解 编写程序如下：　　%数字滤波器指标　　wp＝0.25*pi； %滤波器的通带截止频率　　ws＝0.4*pi；%滤波器的阻带截止频率　　Rp＝1；As＝15；%输入滤波器的通阻带衰减指标 　　ripple＝10^(－Rp/20)；%计算通带衰减对应的幅度值　　Attn＝10^(－As/20)；%计算阻带衰减对应的幅度值　　%转换为模拟滤波器指标　　Fs＝2000；T＝1/Fs；　　Omgp＝wp*Fs；Omgs＝ws*Fs；　　%模拟原型滤波器计算　　［n，Omgc］＝buttord(Omgp，Omgs，Rp，As，￠s￠)  %计算阶数n和截止频率　　［z0，p0，k0］＝buttap(n)； %设计归一化的巴特沃斯 模拟原型滤波器 　　ba1＝k0*real(poly(z0))； %求原型滤波器系数b　　aa1＝real(poly(p0))； %求原型滤波器系数a　　［ba，aa］＝lp2lp(ba1，aa1，Omgc)； %变换为模拟低 通滤波器　　%用脉冲响应不变法计算数字滤波器系数　　［bd，ad］＝impinvar(ba，aa，Fs)　　［C，B，A］＝dir2par(bd，ad) %转换成并联型　　%求数字系统的频率特性　　［H，w］＝freqz(bd，ad)； 　　dbH＝20*log10((abs(H)＋eps)/max(abs(H)))； %化为分贝值　　subplot(2，2，1)，plot(w/pi，abs(H))；　　subplot(2，2，2)，plot(w/pi，angle(H)/pi)；　　subplot(2，2，3)，plot(w/pi，dbH)；　　subplot(2，2，4)，zplane(bd，ad)； 　　程序结果如下：n＝ 6Omgc＝1.8897e＋003bd＝－0.0000 0.0031 0.0419 0.0569 0.0125 0.0003ad＝ 1.0000－2.5418 3.1813－2.3124 1.0072－0.2457 0.0260　　C＝ ［］　　B＝ 2.4935[KG*2]－0.5514　　－2.8792 [KG*2]1.1587　　 0.3857[KG*2]－0.5987 　　A＝ 1.0000[KG*2]－0.7790[KG*2]0.1612　　 1.0000[KG*2]－