kfft.asv

给出了数字信号处理课程设计的基本程序

function y=KFFT(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)=x(n);
n1=n1+1;
end
y(N)=x(N);
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