dft.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 12 行

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;