《时域窗口法设计与频域窗口法设计中设计指标与滤波器阶数的关系.docx》由会员分享,可在线阅读,更多相关《时域窗口法设计与频域窗口法设计中设计指标与滤波器阶数的关系.docx(9页珍藏版)》请在三一文库上搜索。
1、精品1、时域窗口法设计中设计指标与滤波器阶数的关系我主要用Matlab工具模拟了一下用Hamming窗口法设计一个指定设计指标的数字滤波器。并通过改变程序里的阶数N的值,观察并研究了时域窗口法设计中设计指标与滤波器阶数的关系。Matlab仿真程序如下:h=fiU(76,(4.43*2)/22.05 ,hamming(77);H,w=freqz(h,1,40000);plot(w*22.05/(2*pi),20*log10(abs(H),r);%axis(0,10,-100,10);%axis(3.5,4.5,-0.5,0.1);axis(4.5,5.5,-80,-20);可编辑N=76时系统的
2、频率响应当信号频率达到5K时,阻带衰减为-56.06dB当信号频率达到4K时,通带插损为0.09503dB进一步优化:h=fiU(64,(4.43*2)/22.05 ,hamming(65);H,w=freqz(h,1,40000);plot(w*22.05/(2*pi),20*log10(abs(H),r);%axis(0,10,-100,10);%axis(3.5,4.5,-0.5,0.1);axis(4.5,5.5,-80,-20);N=64时系统的频率响应当信号频率达到5K时,阻带衰减为-52.93dB当信号频率达到4K时,通带插损为0.03048同理N取为60,截止频率不变:取数字滤
3、波器的阶数N=60不能满足题目的设计指标。滤波器的通带插损随滤波器阶数增加而减小,减小而增加;阻带衰减随滤波器阶数增加而增加, 减小 而减小。2、频域窗口法设计中设计指标与滤波器阶数的关系Ws=1.425 ; Wp=1.14; wc=1.2823 d=0.14设定阶数N=40时编制仿真程序如下:N=40;wc=1.2823;d=0.14;n=1:N;h0=sin(wc*n)./(pi*n);x=2*sin(0.5*d*n)./(d*n);hp=h0.*x.*x.*x.*x;h=fliplr(hp) wc/pi hp;n1=-N:N;subplot(4,1,1),stem(n1,h),title
4、(冲激响应);H,w=freqz(h,1);subplot(4,1,2),plot(w,20*log10(abs(H/H(1),r),axis(0,3.5,-80,20),title(幅频特性); subplot(4,1,3),plot(w,20*log10(abs(H/H(1),r),axis(1,1.2,-0.5,0.1),title(通带特性);subplot(4,1,4),plot(w,20*log10(abs(H/H(1),r), axis(131.5,-80,-20),title( 阻带特性);当频率为阻带截止频率点时,频率响应为-27.9 ,不能达到阻带衰减为50dB的设计要求;
5、当频率为通带截止频率点的时候,频率响应为 -0.3592dB ,大于通 带插损,也不满足设计指标。取阶数为 30 阶,截止频率为 1.26 , d 为 0.118 时,程序如下:N=30;wc=1.26;d=0.118;n=1:N;h0=sin(wc*n)./(pi*n);x=2*sin(0.5*d*n)./(d*n);hp=h0.*x.*x.*x.*x;h=fliplr(hp) wc/pi hp;H,w=freqz(h,1,5000);n1=-N:N;subplot(4,1,1),stem(n1,h),title( 冲激响应 );H,w=freqz(h,1);subplot(4,1,2),p
6、lot(w,20*log10(abs(H/H(1),r),axis(0,3.5,-80,20),title( 幅频特性 );subplot(4,1,3),plot(w,20*log10(abs(H/H(1),r),axis(1,1.2,-0.5,0.1),title( 通带特性 );subplot(4,1,4),plot(w,20*log10(abs(H/H(1),r),axis(1.3,1.5,-80,-20),title( 阻带特性 );仿真图像如图所示:J %口;q*ia北:tajr;。1口 a4口二傅?中门下不1 .-:工土”公/a / 口守口二匕a %; /色 口。中 /门力::1凋
7、乾竹当频率为阻带截止频率点时,频率响应为-49.6 ,刚达到阻带衰减为 50dB的设计要求;当频率为通带截止频率点的时候,频率响应为 -0.3253dB取阶数为60阶,截止频率为1.26 , d为0.118时,程序如下:N=60;wc=1.26;d=0.118;n=1:N;h0=sin(wc*n)./(pi*n);x=2*sin(0.5*d*n)./(d*n);hp=h0.*x.*x.*x.*x;h=fliplr(hp) wc/pi hp;H,w=freqz(h,1,5000);n1=-N:N;subplot(4,1,1),stem(n1,h),title(冲激响应);H,w=freqz(h,
8、1);subplot(4,1,2),plot(w,20*log10(abs(H/H(1),r),axis(0,3.5,-80,20),title(幅频特性);subplot(4,1,3),plot(w,20*log10(abs(H/H(1),r),axis(1,1.2,-0.5,0.1),title(通带特性);subplot(4,1,4),plot(w,20*log10(abs(H/H(1),r),axis(131.5,-80,-20),title( 阻带特性);仿真图像如图所示:当频率为阻带截止频率点时,频率响应为-44.56 ,刚达到阻带衰减为50dB的设计要求;当频率为通带截止频率点的时候,频率响应为 -0.361dB。优化后能满足设计要求的指标。当截止频率和d保持该优化后的值不变,单 独增加滤波器的阶数时,设计指标不会明显的变化。结论时域窗口法设计的滤波器的阶数越高,结果与设计指标越接近。频域窗口法设计:当设计参数被设定在某一特定值之后,增加阶数对设计指 标的影响较小