fft.asv

来自「给出了数字信号处理课程设计的基本程序」· ASV 代码 · 共 32 行

function y=FFT(x,N)
M=log2(N);
long=length(x);
if(long<N)
     x=[x,zeros(1,N-long)];
end

n1=1;
for I=1:N-1 %100
 J=dec2bin(I,M);
 J=fliplr(J);
n=bin2dec(J);
y(n1)=
n1=n1+1;
end

for m=1:M
    a=2^m;b=(N/a)-1;
    u=exp(-j*2*pi/a);
       for lg=0:b
           c=lg*a;d=c+a/2;e=a/2;
           for k=1:e
               w=u^(k-1);
               T=w*y(c+k+a/2);
               D=y(c+k);
               y(c+k)=D+T;
               y(c+k+a/2)=D-T;
           end
       end
end