dft.m

来自「这是一个dsp的matlab包」· M 代码 · 共 15 行

function [Xk]=dft(xn,N)

%  Computes Discrete Fourier Transform
%  [Xk]=dft(xn,N)
%  0<=k<=N-1
%  xn=N-point finite-duration sequence
%  N=Length of DFT

m=[0:1:N-1];
K=[0:1:N-1];

WN=exp(-j*2*pi/N); % Wn factors
nk=m'*K;
WNnk=WN.^nk;       % DFT matrix
Xk=xn*WNnk;