dfs.m

原名：数字信号处理基础及MATLAB实现,解压缩密码：[e-mail]

function [Xk]=dfs(xn,N)
%计算离散傅立叶级数（DFS）系数
%xn=周期信号在主值区间的单周期信号，N为xn的周期
n=[0:1:N-1];               %行向量数目
k=[0:1:N-1];               %列向量数目
WN=exp(-j*2*pi/N);         %WN的因子
nk=n'*k;                   %产生一个含nk值的N乘N维矩阵
WNnk=WN.^nk;               %DFS矩阵
Xk=xn*WNnk;                %DFS系数行向量