Matlab的窗函数,矩形窗,三角窗,汉明窗,汉宁窗,布莱克曼窗.pdf

%N=51

%求矩形窗的频率响应图

wn=rectwin(51);%矩形窗函数

%20*log10(abs(WN))

[h1,w]=freqz(wn,1);

figure(1);

plot(w/pi,20*log10(abs(h1/max(h1))));

axis([01-1000]);

xlabel(归一化频率

title(矩形窗的傅里叶变换);

set(gca,YTick,[-100-80-60-40-200])

set(gca,XTick,[0:0.2:1])

set(gca,XAxisLocation,top);%设置X轴在上方

set(gca,YAxisLocation,left);%设置Y轴在左方

%求三角窗的频率响应图

wn1=bartlett(51);

[h1,w1]=freqz(wn1,1);

figure(2);

plot(w/pi,20*log10(abs(h1/max(h1))));

axis([01-1000]);

1

xlabel(归一化频率

title(三角窗的傅里叶变换);

set(gca,YTick,[-100-80-60-40-200])

set(gca,XTick,[0:0.2:1])

set(gca,XAxisLocation,top);%设置X轴在上方

set(gca,YAxisLocation,left);%设置Y轴在左方

%hanning

wn1=hanning(51);

[h1,w1]=freqz(wn1,1);

figure(3);

plot(w/pi,20*log10(abs(h1/max(h1))));

axis([01-1000]);

xlabel(归一化频率

title(Hanning的傅里叶变换);

set(gca,YTick,[-100-80-60-40-200]);

set(gca,XTick,[0:0.2:1]);

set(gca,XAxisLocation,top);%设置X轴在上方

set(gca,YAxisLocation,left);%设置Y轴在左方

%hamming

wn1=hamming(51);

2

[h1,w1]=freqz(wn1,1);

figure(4);

plot(w/pi,20*log10(abs(h1/max(h1))));

axis([01-1000]);

xlabel(归一化频率

title(Hamming的傅里叶变换);

set(gca,YTick,[-100-80-60-40-200])

set(gca,XTick,[0:0.2:1])

set(gca,XAxisLocation,top);%设置X轴在上方

set(gca,YAxisLocation,left);%设置Y轴在左方

%Blackman

wn1=blackman(51);

[h1,w1]=freqz(wn1,1);

figure(5);

plot(w/pi,20*log10(abs(h1/max(h1))));

axis([01-1000]);

xlabel(归一化频率

title(Blackman的傅里叶变换);

set(gca,YTick,[-100-80-60-40-200])

set(gca,XTick,[0:0.2:1])

3

set(gca,XAxisLocation,top);%设置X轴在上方

set(gca,YAxisLocation,left);%设置Y轴在左方