mynfft.m

来自「自编实现快速fft程序」· M 代码 · 共 17 行

function C=mynfft(A)
%i=1:256;j=1:256;
n=length(A);
B=zeros(n);
A=double(A);
% unitr=ones(1,n);
% unitl=ones(n,1);
% for i=1:n
%     unitr(1,i)=unitr(1,i)+i-2;
%     unitl(i,1)=unitl(i,1)+i-2;
% end
for u=1:n
    for v=1:n
        B(u,v)=exp(-j*2*pi*u*v/n);
    end
end
C=B*A*B;