dft.m

MATLAB仿真技术与实例应用教程

function [Xk]=dft(xn,N)
%计算离散傅立叶变换
%Xk为在0<=k<=N-1之间的dft系数数组
%xn为N个点的有限序列；N为dft的长度

n=0:N-1;
k=0:N-1;
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;