function?[Xk]=myfft（XnN）
%對(duì)Xn進(jìn)行L點(diǎn)的DFT
%用時(shí)間抽取的基2FFT-DFT進(jìn)行快速傅里葉變換
if?nargin<2
????N=length（Xn）;
end

%序列補(bǔ)零
Xn=[Xnzeros（1N-length（Xn））];???
%計(jì)算需分為蝶形運(yùn)算的級(jí)數(shù)M
tempbinary=dec2bin（N）;
M=length（tempbinary（2:end））;?
tempdata=zeros（NM+1）;?%定義中間操作存儲(chǔ)單元

%%對(duì)輸入進(jìn)行排序
for?n=0:N-1
????tempbinary=dec2bin（nM）;
????tempbinary=rot90（tempbinary2）;
????tempdata（n+11）=Xn（bin2dec（tempbinary）+1）;??%數(shù)組都從1開(kāi)始計(jì)數(shù)
end
%開(kāi)始計(jì)算
for?m=0:M-1
%產(chǎn)生W因子