matlab5基于MATLAB的IIRDF设计I.docVIP

实验五 基于MATLAB的IIR DF设计I 一、实验目的 1.练习并掌握脉冲响应不变法设计IIR数字滤波器的方法； 2.练习并掌握双线性变换法设计IIR数字滤波器的方法； 5. 参照实验内容1234,设计如下指标的Butteworth滤波器.画出模拟低通/高通/带通/带阻、数字低通/高通/带通/带阻的幅频和相频特性图。 1）低通指标：wp=0.2π,rp=1dB,ws=0.3π,rs=16dB。 2）高通指标：wp=0.3π,rp=1dB,ws=0.2π,rs=30dB。 3）带通指标： wp=[0.3π 0.5π],rp=1dB,ws=[0.2 0.6],rs=30dB。 4）带阻指标： wp=[0.2π 0.6π],rp=1dB,ws=[0.3 0.5],rs=30dB。 低通： Wp=0.2*pi; Ws=0.3*pi; Rp=1; Rs=16; Fs=1000; % 数字滤波器指标 %模拟滤波器设计 wp=Wp*Fs; ws=Ws*Fs; Rp=7; Rs=16; % 模拟滤波器指标:w=ΩT-Ω=w*Fs [n,Wn]=buttord(wp,ws,Rp,Rs,s); % 阶次估计 [z,p,k]=buttap(n); %模拟低通原型 [b,a]=zp2tf(z,p,k); %零极点-直接 [b,a]=lp2lp(b,a,Wn);%低通-低通(去归一化) [H,w]=freqs(b,a); figure(1) subplot(2,1,1),plot(w/2/pi,20*log10(abs(H)));grid xlabel(Frequence (Hz)),ylabel(Magnitude (dB)); subplot(2,1,2),plot(w/2/pi,angle(H)/pi*180);grid xlabel(Frequence (Hz)),ylabel(Phase (degress)); %冲激响应不变法设计数字滤波器 figure(2) [bz,az]=impinvar(b,a,Fs); freqz(bz,az,256); % 数字滤波器的直接设计 figure(3) [n,wn]=buttord(Wp/pi,Ws/pi,Rp,Rs) [b,a]=butter(n,wn); freqz(b,a,256); n = 3 wn = 0.1720 高通： wp=0.3*pi; ws=0.2*pi; rp=1; rs=30; Fs=2000; % Firstly to finish frequency prewarping; wap=2*Fs*tan(wp/2); was=2*Fs*tan(ws/2); [n,wn]=buttord(wap,was,rp,rs,s); [z,p,k]=buttap(n); [b,a]=zp2tf(z,p,k); [bt,at]=lp2hp(b,a,wap); %低通-高通(去归一化) w=0:pi:4000*pi; H=freqs(bt,at,w); subplot(3,1,1),plot(w/2/pi,20*log10(abs(H)));grid title( High-pass AF ); %双线性变换法设计数字滤波器 % Note: s=(2/T)(z-1)/(z+1);T=1,that is 2Fs=1,Fs=0.5; [bz,az]=bilinear(bt,at,Fs); [h,w]=freqz(bz,az,256,1); Subplot(3,1,2); plot(w,20*log10(abs(h)));grid; title( High-pass DF (1)); % Directly to design H(z) by butter.m [n,wn]=buttord(wp/pi,ws/pi,rp,rs); [b,a]=butter(n,wp/pi,high); [h1,w1]=freqz(b,a,256,1); Subplot(3,1,3); plot(w,20*log10(abs(h1)));grid;title( High-pass DF (2)); 带通： fp=[300 500]; fs=[200 600]; rp=1; rs=30;Fs=2000;%模拟技术指标 wp=fp*2*pi/Fs; ws=fs*2*pi/Fs;%数字技术指标 %模拟低通原型滤波器设计 wap=2*Fs*tan(wp./2);was=2*Fs*tan(ws./2);%预畸变 [n,wn]=buttord(wap,was,rp,rs,s); [z,p,k]=buttap(n); [bp,ap]=zp2tf(z,p,k);