資源描述：
《IIR數(shù)字濾波器的設(shè)計(jì)實(shí)驗(yàn).doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、附錄：實(shí)驗(yàn)2IIR數(shù)字濾波器的設(shè)計(jì)1、程序一：fc=300;fr=200;fs=1000;rp=0.8;rs=20;wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=cheb1ord(wc,wt,rp,rs,'s');[B,A]=cheby1(N,rp,wn,'high','s');[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az);f=w*fs/(2*pi);plot(f,20*log10(abs(h)));axis([
2、0,fs/2,-80,10]);grid;xlabel('頻率/Hz');ylabel('幅度/dB')圖2-1切比雪夫高通濾波器幅頻響應(yīng)2、程序二：fc=200;fr=300;fs=1000;rp=1;rs=25;wc=2*pi*fc;wt=2*pi*fr;[N,wn]=buttord(wc,wt,rp,rs,'s');[B,A]=butter(N,wn,'s');[bz1,az1]=impinvar(B,A,fs);[h1,w]=freqz(bz1,az1);wc2=2*fs*tan(2*pi*fc/(2*fs));wt2=2*
3、fs*tan(2*pi*fr/(2*fs));[N,wn]=buttord(wc2,wt2,rp,rs,'s');[B,A]=butter(N,wn,'s');[bz2,az2]=bilinear(B,A,fs);[h2,w]=freqz(bz2,az2);f=w*fs/(2*pi);plot(f,20*log10(abs(h1)),'-.',f,20*log10(abs(h2)),'-');axis([0,fs/2,-80,10]);grid;xlabel('頻率/Hz');ylabel('幅度/dB')圖2-2巴特沃思低通濾波器
4、幅頻響應(yīng)3、程序：（1）切比雪夫型數(shù)字低通濾波器fc=1200;fr=2000;fs=8000;rp=0.5;rs=40;wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=cheb1ord(wc,wt,rp,rs,'s');[B,A]=cheby1(N,rp,wn,'s');[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az);f=w*fs/(2*pi);plot(f,20*log10(abs(h)));axis([0,fs
5、/2,-80,10]);grid;xlabel('頻率/Hz');ylabel('幅度/dB')圖2-3（a）切比雪夫低通濾波器幅頻響應(yīng)（2）巴特沃思型數(shù)字低通濾波器fc=1200;fr=2000;fs=8000;rp=0.5;rs=40;wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=buttord(wc,wt,rp,rs,'s');[B,A]=butter(N,wn,'s');[bz,az]=bilinear(B,A,fs);[h,w]=freqz(b
6、z,az);f=w*fs/(2*pi);plot(f,20*log10(abs(h)));axis([0,fs/2,-80,10]);grid;xlabel('頻率/Hz');ylabel('幅度/dB')圖2-3（b）巴特沃思型數(shù)字低通濾波器幅頻響應(yīng)（3）橢圓型數(shù)字低通濾波器clear;fc=1200;fr=2000;fs=8000;rp=0.5;rs=40;wc=2*fs*tan(2*pi*fc/(2*fs));wt=2*fs*tan(2*pi*fr/(2*fs));[N,wn]=ellipord(wc,wt,rp,rs,'s'
7、);[B,A]=ellip(N,rp,wn,'low','s');[bz,az]=bilinear(B,A,fs);[h,w]=freqz(bz,az);f=w*fs/(2*pi);plot(f,20*log10(abs(h)));axis([0,fs/2,-80,10]);grid;xlabel('頻率/Hz');ylabel('幅度/dB')圖2-3（c）橢圓型數(shù)字低通濾波器幅頻響應(yīng)
