MATLAB的dffir,实验七基于MATLAB的FIRDF设计.ppt

实验七基于MATLAB的FIRDF设计.ppt

9.利用频率采样法设计FIR滤波器1 %1) 频率采样技术: 低通,最优法, T1 = 0.5 % wp = 0.2pi, ws=0.3pi, Rp=0.25dB, As=50dB M = 40; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1,5),0.5,zeros(1,29),0.5,ones(1,4)]; Hdr = [1,1,0,0]; wdl = [0,0.25,0.25,1]; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; H = Hrs.*exp(j*angH); h = real(ifft(H,M)); [db,mag,pha,grd,w] = freqz_m(h,1); [Hr,ww,a,L] = Hr_Type2(h); subplot(1,1,1) subplot(2,2,1);plot(wl(1:21)/pi,Hrs(1:21),'o',wdl,Hdr); axis([0,1,-0.1,1.1]); title('频率样本: M=40,T1=0.5') xlabel(''); ylabel('Hr(k)') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gca,'YTickMode','manual','YTick',[0,0.5,1]); grid subplot(2,2,2); stem(l,h); axis([-1,M,-0.1,0.3]) title('脉冲响应'); ylabel('h(n)');text(M+1,-0.1,'n') subplot(2,2,3); plot(ww/pi,Hr,wl(1:21)/pi,Hrs(1:21),'o'); axis([0,1,-0.1,1.1]); title('振幅响应') xlabel('频率(单位：pi)'); ylabel('Hr(w)') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gca,'YTickMode','manual','YTick',[0,0.5,1]); grid subplot(2,2,4);plot(w/pi,db); axis([0,1,-100,10]); grid title('幅度响应'); xlabel('频率(单位：pi)'); ylabel('分贝数'); set(gca,'XTickMode','Manual','XTick',[0;0.2;0.3;1]); set(gca,'YTickMode','Manual','YTick',[-30;0]); set(gca,'YTickLabelMode','manual','YTickLabels',['30';' 0']) % 2)频率采样技术: 低通, 最优法,T1 & T2 % wp = 0.2pi, ws=0.3pi, Rp=0.25dB, As=50dB % T1 = 0.5925, T2=0.1099 M = 60; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1,7),0.5925,0.11,zeros(1,43),0.11,0.5925,ones(1,6)]; Hdr = [1,1,0,0]; wdl = [0,0.2,0.3,1]; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; angH = [-alpha*(2*pi)/M*k1, alpha*(2*pi)/M*(M-k2)]; H = Hrs.*exp(j*angH); h = real(ifft(H,M)); [db,mag,pha,grd,w] = freqz_m(h,1); [Hr,ww,a,L] = Hr_Type2(h); subplot(1,1,1) subplot(2,2,1);plot(wl(1:31)/pi,Hrs(1:31),'o',wdl,Hdr); axis([0,1,-0.1,1.1]); title('低通: M=60,T1=0.59, T2=0.109') xlabel(''); ylabel('Hr(k)') set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1]) set(gc