Matlab fir iir滤波器的比较_设计程序及图片_原创.docxVIP

技术指标fpl=10000;fsl=11000;fsu=20000;fpu=21000;Fs=50000;Rp=1;Rs=50;时的滤波器 IIR 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); %求模拟带阻滤波器的阶数和通带截止频率 [b,a]=ellip(N,Rp,Rs,Wpo,stop); w=0:0.005*pi:pi; [h,w]=freqz(b,a,w); 阶数N =6 clc;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); %求模拟带阻滤波器的阶数和通带截止频率 [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));grid on;xlabel(\omega/\pi);ylabel(幅度/dB); %也可用plot(20*log10(abs(h)));grid on; thet=unwrap(angle(h)); figure; plot((0:199)/200,thet(1:200));grid on;xlabel(\omega/\pi);ylabel(相频响应); x=sin(2*pi*(1:256)*7000/Fs); %输入信号：0=n=255,f0=7kHz; yi=filter(b,a,x); figure; plot(yi(1:250));grid on;xlabel(n);ylabel(y(n)); 幅频响应 相频响应 系统输出 IIR f0=12kHz(在阻带) f0=7kHz（在通带） FIR clc;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));grid on;xlabel(\omega/\pi);ylabel(幅度/dB); %也可用plot(20*log10(abs(h)));grid on; thet2=unwrap(angle(hf)); figure; plot((0:199)/200,thet2(1:200));grid on;xlabel(\omega/\pi);ylabel(相频响应); x=sin(2*pi*(1:256)*12000/Fs); %输入信号：0=n=255,f0=7kHz; yi=filter(h1,1,x); figure; plot(yi(1:250));grid on;xlabel(n);ylabel(y(n)); 幅频响应 相频响应 阶数NF=166 系统输出 f0=12kHz f0=7kHz（在通带中） Fs和滤波器阶数的关系 Fs=50kHz IIR N=6 FIR 阶数NF=166 Fs=100kHz IIR N=7 FIR NF=331 （3）Fs=150kHz N=7 FIR NF=4