dfs.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

function [Xk]=dfs(xn,N)
% calculate the dfs of xn
%[Xk]=dfs(xn,N)
%------------------------
%
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;