%用二重循环实现DFT: function xk=dt_0(xn); %define a function N=length(xn); %caculate the length of the variable WN=exp(-j.*.*pi./N); xk=zeros(,N); %define a non-zero 一维矩阵 sum=zeros(,N); %define a non-zero 一维矩阵 :N %二重循环实现离散傅里叶变换DFT :N sum(n)=xn(n).*WN.^(k.*n