Matlabfiriir滤波器比较设计程序及图片原创.docxVIP

技术指标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); %求模拟带阻滤波器的阶数和通带截止频率[b,a]=ellip(N,Rp,Rs,Wpo,stop);w=0:0.005*pi:pi;[h,w]=freqz(b,a,w);阶数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); %求模拟带阻滤波器的阶数和通带截止频率[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));幅频响应相频响应系统输出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));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=12kHzf0=7kHz（在通带中）Fs与滤波器阶数的关系Fs=50kHzIIRN=6FIR阶数NF=166Fs=100kHzIIRN=7FIRNF=331（3）Fs=150kHzN=7FIRNF=496（4）Fs=2MHzN=7FIRNF=6601验证了书上（1）P251中间几行“无论是设计IIR数字滤波器还是FIR。。。。。”（2）P251第三