%《數字信號處理教程——MATLAB釋義與實現》第三章例3.2.7程序hc327
% 信號序列的DTFT
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
n = 0:10; x = (0.8*exp(j*pi/3)).^n;         % 給定輸入序列x
k = -200:200; w = (pi/100)*k;               % 很密的頻率下標和頻率向量
X = x * (exp(-j*pi/100)) .^ (n'*k);         % x的DTFT計算
x1 = (0.8).^n;                              % 給定輸入序列x的實部
X1 = x1 * (exp(-j*pi/100)) .^ (n'*k);       % x1的DTFT計算
subplot(4,2,1); stem(n,real(x));grid		% 以下繪制時間序列
title('復序列'),ylabel('實部')
subplot(4,2,3); stem(n,imag(x));grid
ylabel('虛部')
subplot(4,2,2); stem(n,real(x1));grid
title('實序列')
subplot(4,2,4); stem(n,imag(x1));grid
subplot(4,2,5); plot(w,abs(X));grid		    % 以下繪制頻譜圖
ylabel('幅度')
subplot(4,2,7); plot(w,angle(X));grid
ylabel('相角')
subplot(4,2,6); plot(w,abs(X1));grid
subplot(4,2,8); plot(w,angle(X1));grid
set(gcf,'color','w')                        % 置圖形背景色為白