function [Xk] = dft(xn,N)
% 計(jì)算離散付利葉變換
% -----------------------------------
% [Xk] = dft(xn,N)
% Xk = 在0 <= k <= N-1間的DFT 系數(shù)數(shù)組 
% xn = N點(diǎn)有限長度序列 
%  N = DFT的長度
%
n = [0:1:N-1];                       % n的行向量
k = [0:1:N-1];                       % k的行向量
WN = exp(-j*2*pi/N);                 % Wn 因子
nk = n'*k;                           % 產(chǎn)生一個(gè)含nk值的N 乘 N維矩陣
WNnk = WN .^ nk;                     % DFT 矩陣
Xk = xn * WNnk;                      % DFT 系數(shù)的行向量