%《數字信號處理教程——MATLAB釋義與實現》第七章演示程序fgp742
% 給定理想低通濾波器的13個頻點圖
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
y=[ones(1,3),zeros(1,8),ones(1,2)];     % 13個頻點的幅度
nx=0:12;                                % 分成的頻譜自變量序號
x=nx/13*2*pi;                           % 13點的頻率值
subplot(2,1,1)                          % 繪理想頻點頻譜圖
plot(x,y,'o')
line([0,pi/3,pi/3,5*pi/3,5*pi/3,2*pi],[1,1,0,0,1,1])  % 繪理想的連續頻譜圖
line([0,2*pi],[0,0])                    % 畫水平基線
line([pi,pi],[-0.1,0.1])                % 畫對稱中心線
axis([0,2*pi,-0.2,1.2])                 % 定坐標范圍
subplot(2,1,2)
theta=-12*pi*[0:12]/13;                 % 13個頻點的相位
plot(x,theta),hold on                   % 繪理想的連續相位圖
plot(x,theta,'o')                       % 繪理想頻點相位圖
axis([0,2*pi,-12*3.5,6])                % 定坐標范圍
grid
set(gcf,'color','w')                    % 置圖形背景色為白