Matlab fir iir濾波器的比較-設(shè)計程序及圖片-原創(chuàng)

技術(shù)指標(biāo)fpl=10000;fsl=11000;fsu=20000;fpu=21000;Fs=50000;Rp=1;Rs=50;時的濾波器IIRfpl=10000;fsl=11000;fsu=20000;fpu=21000;Fs=50000;Rp=1;Rs=50;wpl=2*pi*fpl/Fs;wsl=2*pi*fsl/Fs;wsu=2*pi*fsu/Fs;wpu=2*pi*fpu/Fs;wp=[wpl/pi,wpu/pi];ws=[wsl/pi,wsu/pi];[N,Wpo]=ellipord(wp,ws,Rp,Rs);%求模擬帶阻濾波器的階數(shù)和通帶截止頻率[b,a]=ellip(N,Rp,Rs,Wpo,'stop');w=0:0.005*pi:pi;[h,w]=freqz(b,a,w);階數(shù)N=6clc;clear;fpl=10000;fsl=11000;fsu=20000;fpu=21000;Fs=50000;Rp=1;Rs=50;wpl=2*pi*fpl/Fs;wsl=2*pi*fsl/Fs;wsu=2*pi*fsu/Fs;wpu=2*pi*fpu/Fs;wp=[wpl/pi,wpu/pi];ws=[wsl/pi,wsu/pi];[N,Wpo]=ellipord(wp,ws,Rp,Rs);%求模擬帶阻濾波器的階數(shù)和通帶截止頻率[b,a]=ellip(N,Rp,Rs,Wpo,'stop');w=0:0.005*pi:pi;[h,w]=freqz(b,a,w);hw1=20*log10(abs(h));plot((0:199)/200,hw1(1:200));gridon;xlabel('\omega/\pi');ylabel('幅度/dB');%也可用plot(20*log10(abs(h)));gridon;thet=unwrap(angle(h));figure;plot((0:199)/200,thet(1:200));gridon;xlabel('\omega/\pi');ylabel('相頻響應(yīng)');x=sin(2*pi*(1:256)*7000/Fs);%輸入信號：0<=n<=255,f0=7kHz;yi=filter(b,a,x);figure;plot(yi(1:250));gridon;xlabel('n');ylabel('y(n)');幅頻響應(yīng)相頻響應(yīng)系統(tǒng)輸出IIRf0=12kHz(在阻帶)f0=7kHz（在通帶）FIRclc;clear;fpl=10000;fsl=11000;fsu=20000;fpu=21000;Fs=50000;Rp=1;Rs=50;wpl=2*pi*fpl/Fs;wsl=2*pi*fsl/Fs;wsu=2*pi*fsu/Fs;wpu=2*pi*fpu/Fs;wp=[wpl/pi,wpu/pi];ws=[wsl/pi,wsu/pi];w=0:0.005*pi:pi;B1=wpu-wsu;B2=wsl-wpl;B=min(B1,B2);NF=ceil(6.6*pi/B);wp1=[(fpl+fsl)/Fs,(fpu+fsu)/Fs];h1=fir1(NF-1,wp1,'stop');[hf,w]=freqz(h1,1,w);freqz(h1,1,w);hw2=20*log10(abs(hf));plot((0:199)/200,hw2(1:200));gridon;xlabel('\omega/\pi');ylabel('幅度/dB');%也可用plot(20*log10(abs(h)));gridon;thet2=unwrap(angle(hf));figure;plot((0:199)/200,thet2(1:200));gridon;xlabel('\omega/\pi');ylabel('相頻響應(yīng)');x=sin(2*pi*(1:256)*12000/Fs);%輸入信號：0<=n<=255,f0=7kHz;yi=filter(h1,1,x);figure;plot(yi(1:250));gridon;xlabel('n');ylabel('y(n)');幅頻響應(yīng)相頻響應(yīng)階數(shù)NF=166系統(tǒng)輸出f0=12kHzf0=7kHz（在通帶中）Fs與濾波器階數(shù)的關(guān)系Fs=50kHzIIRN=6FIR階數(shù)NF=166Fs=100kHzIIRN=7FIRNF=331（3）Fs=150kHzN=7FIRNF=496（4）Fs=2MHzN=7FIRNF=6601驗證了書上（1）P251中間幾行“無論是設(shè)計IIR數(shù)字濾波器還是FIR。。。。?！保?）P251第三行IIR數(shù)字濾波器隨著Fs的提高，階數(shù)提高不明顯。IIR濾波器頻率選擇性越好，其相位的非線性越嚴(yán)重。Fs=50000;Rp=1;Rs=50;fpl=10000;fsl=10250;fsu=20750;fpu=21000;fpl=100