【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');