（全新）数字信号处理--实验三.doc

实验三 IIR数字滤波器的设计 一、实验目的 (1) 掌握双线性变换法及脉冲响应不变法设计IIR数滤波器的具体设计方法及其原理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通IIR数字滤波器的计算机编程。 (2) 观察双线性交换及脉冲响应不变法设计的滤波器的频域特性，了解双线性交换法及脉冲响应不变法的特点。 (3) 熟悉巴特沃斯滤波器、切比雪夫滤波器和椭圆滤波器的频率特性。 二、实验内容 实验中有关变量的定义： 通带边界频率 阻带边界频率 通带波动 最小阻带衰减 采样频率 采样周期 上机实验内容： (1)；设计一切比雪夫高通滤波器，观察其通带损耗和阻带衰减是否满足要求。 程序: fc=300;fr=200;fs=1000;Rp=0.8;At=20; wp=2*pi*fc; ws=2*pi*fr; [N,wn]=cheb1ord(wp,ws,Rp,At,s); [b,a]=cheby1(N,Rp,wn,high,s); [bz,az]=bilinear(b,a,fs); [h,w]=freqz(bz,az); figure(1) f=2*w/pi*fs; plot(w,abs(h)) grid; xlabel(Hz); ylabel(dB); 程序图： 切比雪夫滤波器1型，在通带内有起伏波纹，在阻带内单调衰减。上图满足条件。 (2) ；分别用脉冲响应不变法及双线性变换法设计一巴特沃斯数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 程序： Rp=1;At=25; fc=200;fr=300; fs=1000; wp1=2*pi*fc; ws1=2*pi*fr; [N1,wn1]=buttord(wp1,ws1,Rp,At,s); [b1,a1]=butter(N1,wn1,s); [bz1,az1]=impinvar(b1,a1,fs); [h1,w1]=freqz(bz1,az1);%脉冲响应不变法 wp2=2*fs*tan(0.2*pi); ws2=2*fs*tan(0.3*pi); [N2,wn2]=buttord(wp2,ws2,Rp,At,s); [b2,a2]=butter(N2,wn2,s); [bz2,az2]=bilinear(b2,a2,fs); [h2,w2]=freqz(bz2,az2);%双线性变换法 figure f1=2*w1/pi*fs; f2=2*w2/pi*fs; subplot(2,1,1); plot(f1,abs(h1),r) grid; xlabel(频率/Hz); ylabel(幅度/dB); title(脉冲响应不变法); %plot(f,abs(h1),r,f,abs(h2),y); subplot(2,1,2); plot(f2,abs(h2), y) grid; xlabel(频率/Hz); ylabel(幅度/dB); title(双线性变换法); 程序图： 脉冲响应不变法是使数字滤波器的单位脉冲序列逼近模拟滤波器的冲激响应，但是有频谱周期延拓效应，存在频谱的交叠，即混淆，因此只能用于带限的频响特性 双线性变换法克服了脉冲响应不变法的缺点，使得s平面与z平面是一一对应关系，消除了多值变换性。 (3) 利用双线性变换法分别设计满足下列指标的巴特沃斯型、切比雪夫型和椭圆型数字低通滤波器，并作图验证设计结果： 。 程序： fc=1200;fr=2000;fs=8000; Rp=0.5;At=40;wp=2*pi*fc;ws=2*pi*fr; [N1,wn1]=buttord(wp,ws,Rp,At,s); [b1,a1]=butter(N1,wn1,s); [bz1,az1]=bilinear(b1,a1,fs); [h1,w]=freqz(bz1,az1);%巴特沃斯型，黑色 wp2=2*fc/fs; ws2=2*fr/fs; [N2,wn2]=cheb1ord(wp2,ws2,Rp,At); [b2,a2]=cheby1(N2,Rp,wn2); [h2,w]=freqz(b2,a2);%切比雪夫型,紫色 [N3,wn3]=ellipord(wp,ws,Rp,At,s); [b3,a3]=ellip(N3,Rp,At,wn3,s); [bz3,az3]=bilinear(b3,a3,fs); [h3,w]=freqz(bz3,az3);%椭圆型，红色 figure f=2*w/pi*fs; plot(f,abs(h1),f,abs(h2),f,abs(h3))