dft.m

Signal and Sytem Using Matlab 书的源代码 注释详尽 m文件

%discrete fourier transform
%page 133
%DFT
function Xk=dft(x)
[N,M]=size(x);  %make sure that x is a column vector
if M~=1
    x=x';
    N=M;
end
Xk=zeros(N,1);
n=0:N-1;             %有些疑问，求和过程体现在哪里？？？？？？？？？？？
for k=0:N-1
    Xk(k+1)=exp(-j*2*pi*k*n/N)*x;
end