matlab FFT源代码

function [Xk]=myfft（XnN）
%对Xn进行L点的DFT
%用时间抽取的基2FFT-DFT进行快速傅里叶变换
if nargin<2
    N=length（Xn）;
end

%序列补零
Xn=[Xnzeros（1N-length（Xn））];   
%计算需分为蝶形运算的级数M
tempbinary=dec2bin（N）;
M=length（tempbinary（2:end））; 
tempdata=zeros（NM+1）; %定义中间操作存储单元

%%对输入进行排序
for n=0:N-1
    tempbinary=dec2bin（nM）;
    tempbinary=rot90（tempbinary2）;
    tempdata（n+11）=Xn（bin2dec（tempbinary）+1）;  %数组都从1开始计数
end
%开始计算
for m=0:M-1
%产生W因子