%《數字信號處理教程——MATLAB釋義與實現》第二章例3.2.3程序hc323
% 信號序列的DTFT
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
M=4;N=2*M+1;T=0.5;n=-16:16; 
x=[zeros(1,12),ones(1,N),zeros(1,12)];  % 給出輸入序列
w=[-15:0.1:15]+1e-10;
X=sin(0.5*N*w*T)./sin(0.5*w*T);         % 給出頻譜序列
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/T,pi/T,v(3),v(4)])
xlabel('\Omega'),title('(c) 橫軸放大后幅特性')
set(gcf,'color','w')                    % 置圖形背景色為白