dft.m

用于DSK5416的程序

function [Xk]=dft(xn,N)
% Discrete Fourier transform function
%	[Xk]=dft(xn,N)
% where xn is the time-domain signal x(n)
%   N is the length of sequence
%   Xk is the frequency-domain X(k)
n=[0:1:N-1];
k=[0:1:N-1];
WN=exp(-j*2*pi/N);
nk=n'*k;
WNnk=WN.^nk;
Xk=xn*WNnk;