用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w ....docVIP

用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w=rectwin(N); Wvtool(w); 仿真结果： 从仿真上来看：主瓣宽度为4π/N，第一副瓣比主瓣低13dB，阻带衰减21dB，能量主要集中在旁瓣上，而主瓣的能量较低。而且过渡带比较窄。 2、三角形窗程序： clear all; N=64; w=triang(N); wvtool(w); 仿真结果： 从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低26dB，阻带衰减25dB。 3、汉宁窗程序： clear all; N=64; w=hann(N); wvtool(w); 仿真结果为： 从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低31dB，阻带衰减44dB。汉宁窗的幅度函数由三部分相加，使能量更集中在主瓣中，但代价是主办的宽度加宽了8π/N 4、汉明窗程序： clear all; N=64; w=hamming(N); wvtool(w); 仿真结果为： 从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低40dB，阻带衰减53dB。 这种改进的升余弦窗，能量更加集中在主瓣中，主办的能量约占99.96%，第一旁瓣的峰值比主瓣校40 dB，但主瓣宽度和汉宁窗相同，仍为8π/N 5、布莱克曼窗程序： clear all; N=64; w=blackman(N); wvtool(w); 仿真结果为： 从上图上来看：主瓣宽度为12π/N，第一副瓣比主瓣低57dB，阻带衰减74dB。使旁瓣在进一步抵消，阻带衰减进一步增加，过渡带是矩形窗过渡带的3倍。 6、凯塞窗程序： clear all; N=64; w=kaiser(64,1.5); wvtool(w); 仿真结果为： 能量更加集中在主瓣上，过渡带增大。为依赖于N的参数，选择N可以产生各种过渡带宽和接近最优的阻带衰减。对于相同的N，此窗可以提供不同的过渡带宽，这是其他窗做不到的。 7、第一、二种情况理想低通滤波器的单位冲激响应hd(n)计算的程序： function hd=ideal_lp(wc,M); alpha=(M-1)/2; n=0:1:M-1; m=n-alpha+eps; hd=sin(wc*m)./(pi*m); 利用汉明窗设计第二种情况数字低通滤波器： Wp=0.2pi ws=0.4pi Ap=0.24dB As=50dB clear all; wp=0.2*pi; ws=0.4*pi; tr_width=ws-wp M=ceil(6.6*pi/tr_width)+1 n=0:1:M-1; wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_ham=(hamming(M)); h=hd.*w_ham; [db,mag,pha,w]=freqz_m2(h,[1]); delta_w=2*pi/1000; Ap=-(min(db(1:1:wp/delta_w+1))) As=-round(max(db(ws/delta_w+1:1:501))) subplot(221) stem(n,hd) title(理想单位脉冲响应hd(n)); grid subplot(222) stem(n,w_ham) title(汉明窗w(n)); grid subplot(223) stem(n,h) title(实际单位脉冲响应hd(n));grid subplot(224) plot(w/pi,db) title(幅度响应(dB)); grid axis([0,1,-100,10]) function[db,mag,pha,w]=freqz_m2(b,a) [H,w]=freqz(b,a,1000,whole); H=(H(1:1:501)); w=(w(1:1:501)); mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); 仿真结果为： 用布莱克曼窗设计上面的要求： 程序为： clear all; wp=0.2*pi; ws=0.4*pi; tr_width=ws-wp M=ceil(11*pi/tr_width)+1 n=0:1:M-1; wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_bman=(blackman(M)); h=hd.*w_bman [db,mag,pha,w]=freqz_m2(h,[1]); delta_w=2*pi/1000; Ap=-(min(db(1:1:wp/delta_w+1))) As=-round(max(db(ws/delta_w+1:1:501))) subplot(221) stem(n,hd) title(理想单位