hamming01.m

来自「matla窗口制作」· M 代码 · 共 25 行

%___________________________常见窗函数的实现FIR滤波器_____________
%_______________________    条件：N=40，wc=0.2*pi----------------
N=10;
n=[0:1:N-1];
wc=0.5*pi;
hd=ideal_lp1(wc,N);
%____________________________海明窗函数___________________________
w_ham=(hamming(N))';
h=hd.*w_ham;                               %理想滤波器与窗函数相乘
[dp,mag,pha,grd,w]=freqz_m(h,[1]);         %求幅值响应
subplot(221);
stem(n,hd);
title('理想单位脉冲响应hd(n)');
subplot(222)
stem(n,w_ham);
title('海明窗w(n)');
subplot(223);
stem(n,h);
title('实际单位脉冲响应');
subplot(224);
%plot(ww/pi,Hr);
plot(w/pi,pha/pi);
title('幅度响应/db');
%axis([0,1,-200,10]);
grid on;