function [Xk] = dfs(xn,N)
% 計算離散付利葉級數(DFS)系數
% ---------------------------------------------
% [Xk] = dfs(xn,N)
% Xk = 在0 <= k <= N-1間的DFS 系數數組 
% xn = 周期信號在 0 <= n <= N-1之間的一個單周期信號
%  N = xn 的基本周期
%
n = [0:1:N-1];                       % n的行向量
k = [0:1:N-1];                       % k的行向量
WN = exp(-j*2*pi/N);                 % Wn 因子
nk = n'*k;                           % 產生一個含nk值的N 乘 N維矩陣
WNnk = WN .^ nk;                     % DFS 矩陣
Xk = xn * WNnk;                      % DFS 系數的行向量