语音信号处理（二）：矩形窗和汉明窗时域波形和频域波形

一：通过代码实现

clear all;
h=zeros(10000,1);
h(2001:8000)=1;  % 在2001-8000上令值为1
subplot(2,2,1);
plot(h,'r-');
title('矩形窗时域波形');
xlabel('样点数');
ylabel('幅度');
hold on;
axis([0 10000 0 1.5]);

w=h(2001:2060);
W=fft(w,1024);  % 对w进行1024点fft计算
W2=W/W(1);
W3=20*log10(abs(W2));
W4=2*[0:511]/1024;
W5=[0:1023]/1024;
subplot(2,2,2);
plot(W4,W3(1:512));
% plot(W5,W3);
title('矩形窗幅度特性');
xlabel('归一化频率')
ylabel('幅度/dB');
hold on;

%haming 
ham=hamming(60);
subplot(2,2,3);
plot(ham,'b');
title('汉明窗时域波形');
xlabel('样点数')
ylabel('幅度');
hold on;

H=fft(ham,1024);
H1=H/H(1);
H2=20*log10(abs(H1));
H3=2*[0:511]/1024;
H5=[0:1023]/1024;
subplot(2,2,4);
plot(W4,W3(1:512));
% plot(H5,H2);
title('汉明窗幅度特性');
xlabel('归一化频率')
ylabel('幅度/dB');

二：通过wvtool来显示窗的形状和频域图形

再如：

wvtool(hamming(64),hann(64),gausswin(64));