【matlab】矩形窗/三角窗/hanning窗/hamming窗/blackman窗的频率响应图

%{
--------------------------------------------------------------------------- 

File:Matlab的窗函数,矩形窗                               

功能：降低旁瓣水平
参数： 
--------------------------------------------------------------------------- 
%}
%N =51 
%==========================================================================
%求矩形窗的频率响应图  
%==========================================================================
W = linspace(-pi,pi,4096);
wn0 = rectwin(51)   %矩形窗函数 
%20*log10(abs(WN))  
[h1,w0] = freqz(wn0,1,W); 
%subplotfigure(5,1,1);  
subplot(511);
plot(w0/pi,20*log10(abs(h1/max(h1))));  
axis([-1 1 -100 0]); 
xlabel('归一化频率 /\pi');  
ylabel('20log_{10}|W(e^{j\omega})| /dB');  
title('矩形窗的傅里叶变换'); 
set(gca,'YTick',[-100 -80 -60 -40 -20 0])  
set(gca,'XTick',[-1 :0.2: 1])   
%set(gca,'XAxisLocation','top');%设置X轴在上方  
%set(gca,'XAxisLocation','buttom');%设置X轴在下方  
set(gca,'YAxisLocation','left'); %设置Y轴在左方  
text(1,-124,'\pi');%gtext('\pi');

%==========================================================================
%求三角窗的频率响应图  
%==========================================================================
wn1 = bartlett(51)
[h1,w1] = freqz(wn1,1,W);  
%figure(5,1,2); 
subplot(512);
plot(w1/pi,20*log10(abs(h1/max(h1))));  
%plot(w/pi,20*log10(h1/max(h1))); % 警告: 复数 X 和/或 Y 参数的虚部已忽略 
axis([-1 1 -100 0]); 
xlabel('归一化频率 /\pi');  
ylabel('20log_{10}|W(e^{j\omega})| /dB'); 
title('三角窗的傅里叶变换');   
set(gca,'YTick',[-100 -80 -60 -40 -20 0]) 
set(gca,'XTick',[-1 :0.2: 1])  
%set(gca,'XAxisLocation','top');%设置X轴在上方  
set(gca,'YAxisLocation','left'); %设置Y轴在左方  
text(1,-124,'\pi');%gtext('\pi');
%==========================================================================
%hanning 窗的频率响应图   
%==========================================================================
wn2 = hanning(51)  
[h1,w2] = freqz(wn2,1,W);  
%figure(5,1,3); 
subplot(513);
plot(w2/pi,20*log10(abs(h1/max(h1))));  
axis([-1 1 -100 0]); 
xlabel('归一化频率 /\pi');  
ylabel('20log_{10}|W(e^{j\omega})| /dB');  
title('Hanning的傅里叶变换');
set(gca,'YTick',[-100 -80 -60 -40 -20 0]);  
set(gca,'XTick',[-1 :0.2: 1]);  
%set(gca,'XAxisLocation','top');%设置X轴在上方  
set(gca,'YAxisLocation','left'); %设置Y轴在左方 
text(1,-124,'\pi');%gtext('\pi');
%==========================================================================
%hamming 窗的频率响应图    
%==========================================================================
wn3 = hamming(51) 
[h1,w3] = freqz(wn3,1,W); %离散系统频响特性的函数freqz()
%figure(5,1,4);  
subplot(514);
plot(w3/pi,20*log10(abs(h1/max(h1)))); 
axis([-1 1 -100 0]); 
xlabel('归一化频率 /\pi'); 
ylabel('20log_{10}|W(e^{j\omega})| /dB');  
title('Hamming的傅里叶变换');
set(gca,'YTick',[-100 -80 -60 -40 -20 0])  
set(gca,'XTick',[-1 :0.2: 1])  
%set(gca,'XAxisLocation','top');%设置X轴在上方  
set(gca,'YAxisLocation','left'); %设置Y轴在左方 
text(1,-124,'\pi');%gtext('\pi');
%==========================================================================
%Blackman   wn1 = blackman(51) 
%==========================================================================
wn4 = blackman(51)
[h1,w4] = freqz(wn4,1,W); 
%figure(5,1,5); 
subplot(515);
plot(w4/pi,20*log10(abs(h1/max(h1))));   
axis([-1 1 -100 0]);  
xlabel('归一化频率 /\pi'); 
ylabel('20log_{10}|W(e^{j\omega})| /dB');  
title('Blackman的傅里叶变换'); 
set(gca,'YTick',[-100 -80 -60 -40 -20 0]) 
set(gca,'XTick',[-1 :0.2: 1]) 
%set(gca,'XAxisLocation','top');%设置X轴在上方  
set(gca,'YAxisLocation','left'); %设置Y轴在左方 
text(1,-124,'\pi');%gtext('\pi');