第6部分无限冲激响应(IIR)滤波器设计课程.pptVIP

例6.7.1(例6.5.1) clear all; fp=100;fs=300;Fs=1000;rp=3;rs=20; wp=2*pi*fp/Fs;ws=2*pi*fs/Fs; Fs=Fs/Fs; % let Fs=1 % Firstly to finish frequency prewarping ; wap=tan(wp/2);was=tan(ws/2); % [n,wn]=buttord(wap,was,rp,rs,s) % Note: s! [z,p,k]=buttap(n); % [bp,ap]=zp2tf(z,p,k) % [bs,as]=lp2lp(bp,ap,wap) % % Note: s=(2/Ts)(z-1)/(z+1);Ts=1,that is 2Fs=1,Fs=0.5; [bz,az]=bilinear(bs,as,Fs/2) % [h,w]=freqz(bz,az,256,Fs*1000); plot(w,abs(h));grid on; 例6.7.1(例6.5.1) clear all; wp=.2*pi;ws=.6*pi;Fs=1000; rp=3;rs=20;% % Firstly to finish frequency prewarping; wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2); [n,wn]