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

1、Matlab fir iir滤波器的比较设计程序及图片原创技术指标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); %求模拟带阻滤波器的阶数

2、和通带截止频率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(

3、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); %输入信号

4、：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

5、: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);figu

6、re;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系统输出（1） f0=12kHz（2） f0=7kHz（在通带中）Fs与滤波器阶数的关系（1） Fs=50kHzIIRN=6FIR阶数NF=166（2） Fs=100kHzIIRN=7

7、FIRNF=331（3）Fs=150kHzN=7FIRNF=496（4）Fs=2MHzN=7FIRNF=6601验证了书上（1）P251中间几行“无论是设计IIR数字滤波器还是FIR。”（2）P251第三行 IIR数字滤波器随着Fs的提高，阶数提高不明显。IIR滤波器频率选择性越好，其相位的非线性越严重。Fs=50000;Rp=1;Rs=50;fpl=10000;fsl=10250;fsu=20750;fpu=21000;fpl=10000;fsl=11000;fsu=20000;fpu=21000;fpl=10000;fsl=12000;fsu=19000;fpu=21000;fpl=10000;fsl=13000;fsu=18000;fpu=21000;