数字信号处理窗函数法设计FIR数字滤波器实验.doc

N=15,33 N=15; wc=pi/4; a=(N-1)/2; n=0:(N-1); m=n-a+eps; hdn=sin(wc*m)./(pi*m); wn=hanning(N); hn=hdn.*(wn'); [H,w]=freqz(hn,[1],1024,'whole'); dbH=20*log10((abs(H)+eps)/max(abs(H))); figure(1);subplot(2,2,1); stem(n,hn,'.'); xlabel('n');ylabel('h(n)');title('N=15时设计汉宁窗h(n)'); subplot(2,2,2); plot(w,abs(H));