% 程序3.3：haming.m
x=linspace(20,80,61);         %在20~80的横坐标间取61个值作为横坐标点
h=hamming(61);             %取61个点的哈明窗值为纵坐标值
figure(1);                   %画图
subplot(1,2,1);               %第一个子图
plot(x,h,'k');                 %横坐标为x，纵坐标为h，k表示黑色
title('Hamming窗时域波形');   %图标题
xlabel('样点数');              %横坐标名称
ylabel('幅度');             %纵坐标名称
w1=linspace(0,61,61);         %取窗长内的61个点
w1(1:61)=hamming(61);         %加哈明窗
w2=fft(w1,1024);               %对时域信号进行1024点傅立叶变换
w3=w2/w2(1);                 %幅度归一化
w4=20*log10(abs(w3))           %对归一化幅度取对数
w=2*[0:1023]/1024;            %频率归一化
subplot(1,2,2)                  %画第二个子图
plot(w,w4,'k')                  %画幅度特性图
axis([0,1,-100,0])                %限定横、纵坐标范围
title('Hamming窗幅度特性');      %图标题
xlabel('归一化频率 f/fs');      %横坐标名称
ylabel('幅度/dB');             %纵坐标名称
