dft.m

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

function [Xk]=dft(xn,N)
%计算离散傅立叶变换（DFT）系数
%xn=N点有限长序列，N序列长度
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;               %DFT矩阵
Xk=xn*WNnk;                %DFT系数行向量