fig12_13.m

来自「基于WEB与MATLAB编写的信号处理的源码,书名为"应用WEB和MATLAB的」· M 代码 · 共 11 行

% Figure 12.13
%
N = 10;
Omega = -pi:2*pi/300:pi;
Omega = Omega + (Omega==0)*eps;   % replaces Omega = 0 with eps
W = sin(Omega*(N)/2).*exp(-j*Omega*(N-1)/2)./sin(Omega/2);
plot(Omega,abs(W));
axis([-pi pi 0 10])
title('Figure 12.13')
xlabelg('W')
ylabel('|W|')