% 《數字信號處理及其MATLAB實現》第三章題3_1題解
function [X]=dtft(x,w)
% 計算離散時間傅立葉變換
% [X]=dtft(x,n,w)
% X = 在w頻率點上的DTFT數組
% x = 沿n的有限長度序列
% n = 樣本位置向量
% w = 頻率點位置向量

n=1:length(x);
ewn=exp(-n'*w*i);
X=x*ewn;