用matlab自己编写代码实现基2—FFT

function [xk] = my_fft(x,n)
longx=length(x);
if n==0%如果输入序列只有一位，fft是他本身
   xk=x(1)*exp(-2*1j*pi*0*0);
else
   x1=[];
   for i=1:longx/2
       x1=[x1,x(2*i-1)];
   end% 提取离散序列下标为 0 2 4 6...的值，构成偶序列
   x2=[];
   for i=1:longx/2
       x2=[x2,x(2*i)];
   end% 提取离散序列下标为 1 3 5 7...的值，构成偶序列
   x1k=my_fft(x1,n-1);%嵌套自己
   x2k=my_fft(x2,n-1);
   wnk=[];
   for i=1:longx/2
       wnk=[wnk,exp(-2*1j*pi/(longx)*(i-1))];
   end
   x1kk=x1k+wnk.*x2k;%蝶形运算
   x2kk=x1k-wnk.*x2k;
   xk=[x1kk x2kk];
end
end

适用于x为长2^n长度的序列