实验四预作结果IIR滤波器设计.docVIP

实验四 IIR滤波器设计预做实验 实验目的 1．熟悉由模拟滤波器转换为IIR数字滤波器的原理与方法。 2．掌握数字滤波器的计算机仿真方法。 实验内容 1．用脉冲响应不变法设计巴特沃斯数字滤波器。 2．用双线性变换法设计切比雪夫数字滤波器。 3．用双线性变换法设计巴特沃斯数字滤波器。并将直接型结构转换成级联型结构。 4．比较分析两种不同方法设计的数字滤波器的频率特性的区别。 5．比较并分析变换前的模拟滤波器与变换后的数字滤波器的频率特性的区别，模拟滤波器的频率特性可以使用freqs()函数，注意模拟频率和数字频率之间的对应关系。 实验程序及结果 1．利用脉冲响应不变法设计巴特沃思滤波器 wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; OmegaP=wp/T;OmegaS=ws/T; [cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As); [b,a]=imp_invr(cs,ds,T) [db,mag,pha,w]=freqz_m(b,a); subplot(2,1,1);plot(w/pi,mag); title(digital filter Magnitude Response) axis([0,1,0,1.1]) subplot(2,1,2);plot(w/pi,db); title(digital filter Magnitude in DB) axis([0,1,-40,5]); 结果： N1 = 5.8858 N = 6 Butterworth Filter Order= 6 OmegaC = 0.7032 b = 0.0000 0.0006 0.0101 0.0161 0.0041 0.0001 a = 1.0000 -3.3635 5.0684 -4.2759 2.1066 -0.5706 0.0661 上图为利用脉冲响应不变法设计的巴特沃思低通滤波器频率响应的幅度特性，下图为相对幅度（衰减）特性，单位DB 2．利用双线性变换法设计切比雪夫数字滤波器 wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; OmegaP=(2/T)*tan(wp/2); OmegaS=(2/T)*tan(ws/2); ep=sqrt(10^(Rp/10)-1); Ripple=sqrt(1/(1+ep*ep)); Attn=1/(10^(As/20)); A1=1/Attn;a1=sqrt(A1*A1-1)/ep; a2=OmegaS/OmegaP; N=ceil(logm(a1+sqrt(a1*a1-1))/logm(a2+sqrt(a2*a2-1))); fprintf(\n Chebyshev Filter Order=%2.0f\n,N) [cs,ds]=u_chb1ap(N,Rp,OmegaP); [b,a]=bilinear(cs,ds,T) [db,mag,pha,w]=freqz_m(b,a); subplot(2,1,1);plot(w/pi,mag); title(digital filter Magnitude Response); axis([0,1,0,1.1]) subplot(2,1,2);plot(w/pi,db); title(digital filter Magnitude in DB); axis([0,1,-40,5]); 结果： Chebyshev Filter Order= 4 b = 0.0018 0.0073 0.0110 0.0073 0.0018 a = 1.0000 -3.0543 3.8290 -2.2925 0.5507 上图为利用双线性变换法设计的切比雪夫低通滤波器频率响应的幅度特性，下图为相对幅度（衰减）特性，单位DB 3．利用双线性变换法设计巴特沃思数字滤波器，并转换为级联型结构 wp=0.2*pi;ws=0.3*pi;Rp=1;As=15;T=1; OmegaP=(2/T)*tan(wp/2); OmegaS=(2/T)*tan(ws/2); [cs,ds]=afd_butt(OmegaP,OmegaS,Rp,As); [b,a]=bilinear(cs,ds,T) [db,mag,pha,w]=freqz_m(b,a); subplot(2,1,1);plot(w/pi,mag); title(digital filter Magnitude Response); axis([0,1,0,1.1]) subplot(2,1,2);plot(w/pi,db); tit