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因子