汉宁窗设计带阻FIR滤波器.docVIP

wp1=0.3*pi;ws1=0.4*pi;ws2=0.6*pi;wp2=0.7*pi; N=ceil(8*pi/wdelta) N = 80 wdelta=ws1-wp1; wp1=0.3*pi;ws1=0.4*pi;ws2=0.6*pi;wp2=0.7*pi; subplot(3,1,2),impz(b,1) wp1=0.3*pi;ws1=0.4*pi;ws2=0.6*pi;wp2=0.7*pi; wdelta=ws1-wp1; N=ceil(8*pi/wdelta) N = 80 wc=[(ws1+wp1)/2,(ws2+wp2)/2]; n=0:N-1; b=fir1(N,wc/pi,stop,hanning(N+1)); [H,f]=freqz(b,1,512,100); figure(1) subplot(3,1,1),plot(f,20*log10(abs(H))) xlabel(频率/Hz);ylabel(振幅/dB);grid on; title(频率特性); subplot(3,1,2),impz(b,1) xlabel(n);ylabel(h(n));grid on; title(单位脉冲响应); subplot(3,1,3),zplane(b,1) xlabel(Re[z]);yla