脉冲响应不变法.docVIP

实验1用脉冲响应不变法设计IIR滤波器 实验目的 ?加深对脉冲响应不变法设计IIR滤波器数字滤波器基本方法的了解，熟 悉这 一设计的计算机编程。 ?观察用脉冲响应不变法设计的数字滤波器和响应模拟滤波器的时域特性 和频域特性，比较所涉及的数字滤波器和响应的模拟滤波器的频域特性， 观察脉 冲响应不变法设计中产生的频域混淆现象。 3?观察参数变化对滤波器性能的影响；了解脉冲响应不变法的优缺点，熟 悉该方法的应用范围。 实验原理与方法 脉冲响应不变法是实现模拟滤波器数字化的一种直观而常用的方法。它特别 适 合于对滤波器的时域特性有一定要求的场合。 一个模拟滤波器的传递函数可以用有理分式表达式表示为： TOC \o 1-5 \h \z M /八 (3)_ao as ????as ⑴ H a) rbobs ???bS 通过反拉普拉斯变换我们就可以得到它的冲激相应： ha(t)二L a(S)l ( 2) 脉冲响应不变法就是要保证脉冲响应不变，即： h(n)二 ha(nT) (3) 对上式的冲激相应序列h(n)作Z变换，就可以得到数字滤波器的传递函数： H ⑵二 ZT[ha(nT)] (4) 一般来说，“(s)的分母多项式阶次总是大于分子多项式的阶次。假定 H a(S)的没有多重极点，则式(1)就可分解为： xA (5) (5) 式中Si, A均为复数，Si是H a(s)的极点。其拉氏变换为 S「Tn ha(t)二 Ae u(n) i 土 般来说，Ha (S)的分母多项式阶次总是大于分子多项式的阶次。数字滤波 器的传递函数H(z)经过合并简化，成为一般形式的有理分式传递函数 Co _ C | ZJ 一??? 一 c：z丄do iz1 —dxz 在讨论采样序列z变换与模拟信号拉氏变换之间关系的有关章节中， 我们已 经知道 1 :: H(Z) |zj 二匚Ha(s j〒K) 一 T K 二二 按照z二M的关系，每一个s平面上宽度为2 /T的水平条带将重迭映射到z平面 上。因此脉冲响应不变法将s平面映射到z平面，不是一个简单的 对应 的关系。对于高采样频率(T小)的情况，数字滤波器在频域可能有极高的增益。为此 我们采用 TA H(z)；V ST (8) J e刃 N h(n)二 Tha(nT)二TAeSTnu (n) (9) i吕 在脉冲响应不变法设计中，模拟频率与数字频率之间的转换关系是线性的 (??■Vt)。同时，它可以保持脉冲响应不变，h(t)二ha(nT)。因此，这一方法往 往用 于低通时域数字滤波器设计及相应的模拟系统数字仿真设计。 3.实验内容及步骤： 1 ?复习脉冲响应不变法设计IIR滤波器和模拟滤波器设计的有关内容，认真 阅读 本实验原理与方法的内容，熟悉脉冲响应不变法的设计步骤。 2?编制一个用脉冲响应不变法设计IIR滤波器的程序： 采样频率f s二4Hz,周期T二0?25s,设计一个四阶Butterworth低通滤波器，其3dB截 止频率fc=lHz= 采样频率为1Hz,设计一个数字低通滤波器，要求其通带临界频率 仿二0. 2Hz, 通带内衰减小于ldB,阻带临界频率f s二0. 3Hz,阻带衰减大于25dBo a.用Butterworth逼近设计该滤波器 b.用Chebyshev逼近设计该滤波器。 在mat lab上建立M文件如下： T二0? 25;Fs二4;Q二2*pi;0megac=Q/T;N=4;t=[0:0. 01:3]; n=E0:l:3/T]; Lb, a]=butter (N, Omegac, s ) ; Ec, d]=impi nvar (b, a, T); [h, w]二freqs (b, a) ; Lhg, x, t]=impulse (b, a); mag 二 abs(h);mag 二 20*logl0(mag); figure (1);impz(ha); title( h(n) ) ;ylabel (，幅度，);grid on; figure (2) plot (t, hg) ; title( ，单位脉冲响应 h (t) * ) ;grid on ; figure (3) subplot (321) ;plot (w, mag) ; title ( ，幅度响应 Hg(j/0mega) * ) ;grid on; [h, w]=freqz(b, a);mag 二 abs(h);mag = 20*logl0(mag); subplot (322) ;plot(w, mag) ; title ( ，幅度卩冋应 H(eAjw) * ) ;grid on ; %%i%srworth 设计 %%%%%%%%%%%% Fs=l;0megap=0. 4:0megar=0. 6;Ap=l;Ar=25;T=l/