%《數字信號處理教程——MATLAB釋義與實現》第三章例3.2.4程序hc324
% 信號序列的DTFT
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
x=[2,-1,1,1];nx=0:3;                % 給定x序列
K=64;dw=2*pi/K;                     % 確定頻率分辨率dw
k=floor((-K/2+0.5):(K/2-0.5));      % 設定對稱的頻率向量下標
X=x*exp(j*dw*nx'*k);                % 求序列x的DTFT
subplot(2,1,1),plot(k*dw,abs(X))    % 繪圖
xlabel('頻率'),ylabel('幅頻特性')
subplot(2,1,2),plot(k*dw,angle(X))
xlabel('頻率'),ylabel('相頻特性')
set(gcf,'color','w')                % 置圖形背景色為白