%《數字信號處理教程——MATLAB釋義與實現》第二章例3.2.3程序hc323
% 信號序列的DTFT
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
M=4;N=2*M+1;n=-4*M:4*M; 
x=[zeros(1,3*M),ones(1,N),zeros(1,3*M)];  % 給出輸入序列
w=[-8:0.1:8]+1e-10;
X=sin(0.5*N*w)./sin(0.5*w);         % 給出頻譜序列
subplot(1,3,1),stem(n,x,'.')            % 畫出輸入序列
axis([-20,20,-0.1,1.1]),grid on
xlabel('n'),title('(a) 序列幅度')
subplot(1,3,2),plot(w,X),grid on        % 畫出頻譜序列
xlabel('\omega'),title('(b) 幅頻特性')
subplot(1,3,3),plot(w,X),grid on        % 改變橫軸比例，畫出頻譜序列
v=axis;axis([-pi,pi,v(3),v(4)])
xlabel('\omega'),title('(c) 橫軸放大后幅特性')
set(gcf,'color','w')                    % 置圖形背景色為白