DSP实验三【DOC精选】.doc

一、用双线性变换法设计IIR滤波器相应的程序和波形 （一）高通数字滤波器 clear t=1; fs=1/t fp=0.3 fr=0.2 ap=0.8 ar=20 wp=2*pi*fp/fs; wr=2*pi*fr/fs; wpp=2/t*tan(wp/2); wrr=2/t*tan(wr/2); [n,wn]=cheb1ord(wpp,wrr,ap,ar,s); [b,a]=cheby1(n, ap,wn, high,s); [bz,az]=bilinear(b,a,fs); freqz(bz,az,512,fs);  （二）、低通数字滤波器 （1）用buttworth设计 cleart=1; fs=1/t; fp=0.2; fr=0.3; ap=1; ar=25; wp=2*pi*fp/fs; wr=2*pi*fr/fs; wpp=2/t*tan(wp/2); wrr=2/t*tan(wr/2); [n,wn]=buttord(wpp,wrr,ap,ar,s); [b,a]=butter(n,wn,s); [bz,az]=bilinear(b,a,fs); freqz(bz,az,512,fs);  （2）用chebyshev设计 t=1;fs=1/t; fp=0.2; fr=0.3; ap=1; ar=25; wp=2*pi*fp; wr=2*pi*fr; wpp=2/t*tan(wp/2); wrr=2/t*tan(wr/2); [n,wn]=cheb1ord(wpp,wrr,ap,ar,s); [b,a]=cheby1(n,ap,wn,low,s); [bl,al]=bilinear(b,a,fs) freqz(bl,al,512,fs); （三）、带通数字滤波器 clear fp1=0.2; fp2=0.3; fr1=0.15; fr2=0.35; ap=1; ar=40; t=1; fs=1/t; wp1=2*pi*fp1/fs; wp2=2*pi*fp2/fs; wp=[wp1,wp2]; wr1=2*pi*fr1/fs; wr2=2*pi*fr2/fs; wr=[wr1,wr2]; wpp=2/t*tan(wp/2); wss=2/t*tan(wr/2); [n,wn]=buttord(wpp,wss,ap,ar,s); [b,a]=butter(n,wn,bandpass,s); [bz,az]=bilinear(b,a,fs); freqz(bz,az,512,fs);  二、思考题的回答： 1、双线性变换和脉冲响应不变法相比较，有哪些优点和缺点？为什么？ 答：双线性变换比起脉冲响应不变法来，其主要有点是S平面与Z平面式单值的一一对应关系。整个jΩ轴是单值的对应于单位圆一周，即。当时，ω终止于折叠频率ω=π处，因而双线性变换就不会现由于高频部分超过折叠频率而混乱到低频部分去的现象。单双线性不换法的优点是靠频谱的严重非线性关系而得到的，这种Ω与ω的非线性关系使数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变。例如，如果模拟滤波器是最平延迟滤波器，它的相位与频率是直线关系。但是，通过双线性变换后，所得的数字滤波器就不在保持相位与频率的直线关系了，因而这个数字滤波器也就不再是一个最平延迟特性的滤波器。 2、双线性变换是一种非线性变换，在实验中你观察到这种非线性关系了吗？应该怎样从哪种数字滤波器幅频特性曲线中可以观察到这种非线性关系？ 答：观察到了这种非线性关系。 我们知道双线性变换是靠频率的严重非线性关系得到的。这种Ω与ω 的非线性关系使得数字滤波器与模拟滤波器在响应与频率的对应关系上发生畸变。例如，如果模拟滤波器是最平延时滤波器，它的相位与频率是直线关系。但是通过双线性变换所得到的数字滤波器了。 三、实验总结 通过用双线性变换法设计IIR滤波器实验，我了解两种工程上最常用的变换方法：脉冲响应不变法和双线性变换法及了解多项式乘积和多项式乘方运算的计算机编程方法；掌握了双线性变换法设计IIR滤波器的原理及具体设计方法，熟悉用双线性设计法设计低通、带通和高通IIR数字滤波器的计算机程序以及熟悉用双线性变换法设计数字Butterworth和Chebyshev滤波器的全过程；知道了怎样观察用双线性变换法设计的滤波器的频域特性，并与脉冲响应不变法相比较，了解双线性变换法的特点。 通过此次实验，我加深了对用双线性变换法设计IIR滤波器这部分理论知识的理解及掌握，所以，说实践不仅是检验真理的标准，而且是唯一的标准。更有效的检验学生的学习效果。在此谢谢老师的指导！ 第 2 页 共 11 页