function [Xk] = dfs(xn,N)
% 計(jì)算離散Fourier序列的系數(shù)
% 調(diào)用方式：[Xk] = dfs(xn,N)
% Xk：離散Fourier變換的系數(shù)，序號(hào)：0 <= k <= N-1
% xn：取周期信號(hào)的一個(gè)周期，序號(hào)：0 <= n <= N-1
%  N： xn的數(shù)據(jù)個(gè)數(shù)
n = [0:1:N-1];                       % 數(shù)據(jù)數(shù)組序列
k = [0:1:N-1];                       % 頻率數(shù)組序列
WN = exp(-j*2*pi/N);                 % Wn 因子
nk = n'*k;                          % 產(chǎn)生N*N 的矩陣nk
WNnk = WN .^ nk;                     %離散Fourier變換矩陣
Xk = xn * WNnk;                      % Fourier變換系數(shù)