function [Xk]=dfs(xn,N)
%計(jì)算離散傅立葉級(jí)數(shù)（DFS）系數(shù)
%xn=周期信號(hào)在主值區(qū)間的單周期信號(hào)，N為xn的周期
n=[0:1:N-1];               %行向量數(shù)目
k=[0:1:N-1];               %列向量數(shù)目
WN=exp(-j*2*pi/N);         %WN的因子
nk=n'*k;                   %產(chǎn)生一個(gè)含nk值的N乘N維矩陣
WNnk=WN.^nk;               %DFS矩陣
Xk=xn*WNnk;                %DFS系數(shù)行向量