Matlab的窗函数,矩形窗,三角窗,汉明窗,汉宁窗,布莱克曼窗
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%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([0 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',[0 :0.2: 1])
set(gca,'XAxisLocation','top'); %设置X轴在上方set(gca,'Y AxisLocation','left');%设置Y轴在左方text(1,-108,'\pi');%gtext('\pi');
%求三角窗的频率响应图
wn1 = bartlett(51);
[h1,w1] = freqz(wn1,1);
figure(2);
plot(w/pi,20*log10(abs(h1/max(h1))));
axis([0 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',[0 :0.2: 1])
set(gca,'XAxisLocation','top');%设置X轴在上方set(gca,'Y AxisLocation','left'); %设置Y轴在左方%hanning
wn1 = hanning(51) ;
[h1,w1] = freqz(wn1,1);
figure(3);
plot(w/pi,20*log10(abs(h1/max(h1))));
axis([0 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',[0 :0.2: 1]);
set(gca,'XAxisLocation','top');%设置X轴在上方set(gca,'Y AxisLocation','left'); %设置Y轴在左方%hamming
wn1 = hamming(51) ;
[h1,w1] = freqz(wn1,1);
figure(4);
plot(w/pi,20*log10(abs(h1/max(h1))));
axis([0 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',[0 :0.2: 1])
set(gca,'XAxisLocation','top');%设置X轴在上方set(gca,'Y AxisLocation','left'); %设置Y轴在左方%Blackman
wn1 = blackman(51) ;
[h1,w1] = freqz(wn1,1);
figure(5);
plot(w/pi,20*log10(abs(h1/max(h1))));
axis([0 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',[0 :0.2: 1])
set(gca,'XAxisLocation','top');%设置X轴在上方set(gca,'Y AxisLocation','left'); %设置Y轴在左方
参考书：数字信号处理教程(第三版) 程佩青清华大学出版社2009年11月第9次印刷
P340的图7-11
矩形窗三角窗汉宁窗海明窗布拉克曼窗。