matlab实现DFT(取了一些数据计算结果和fft相差很小，但未详细验证)

对标matlab中的fft函数，但不是快速算法，也没详细验证

function out=myDFT(in)
%in要求列向量
N=length(in);%时域信号元素个数
n=1:N;
X=[];
E=[];
for k=1:1:N
    M=@(n)exp(-1i*2*pi*(k-1)*(n-1)/N);
    E=arrayfun(M,n);%行向量
     X(k,1)=E*in;
end
out=X;

end