dft.m

来自「matlab的一些基础程序」· M 代码 · 共 14 行

function [Xk]=dft(xn,N)
%Computes Discrete Fourier Transform
%------------------------------------------
%[Xk]=dft(xn,N)

%N=Length of DFT
%
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);%Wn factor
nk=n'*k;
WNnk=WN.^nk; 
Xk=xn*WNnk;