%%%copyright:YT_wcj???///基2—按時間抽取
%date:2011/12/08
%corporation:ZYE
%?discription:FFT
%二進制數(shù)第一位是最高位
clear?all;
clc;
N?=?2048;%1024點FFT
m?=?log2（N）;%蝶形的級數(shù)?10
x?=?0:N-1;%FFT計算輸入的原始數(shù)據(jù)?注意是有規(guī)律的數(shù)據(jù)
for?add?=?0:N-1%
????new_add?=?bin2dec（fliplr（dec2bin（addm）））;%%%%%%%變址運算
????new_add1（add+1）=new_add;
????xx（new_add+1）?=?x（add+1）;%xx?相當于緩存器
end
%?for?i_2048=0:2047
%???addr_index（i_2048+1）=find（new_add1==i_2048）-1;
%?end
%?
%?i_1024=1:2:2048;
%?addr_up_1024=new_add1（i_1024）;

reverse_order?=?xx;%排序后的數(shù)據(jù)
%%%%?M=1第一級????旋轉因子的地址都是0
R?=?dec2bin（2^0m）;%-----用二進制表示級數(shù)--r9r8r7r6r5r4r3r2r1r0
butterfly_cnt?=?0:2^（m-1）-1;%-----蝶形計數(shù)器512個
c?=?0;
c?=?dec2bin（cm-1）;%0000?0000?0
wn?=?dec2bin（02*（m-1））;%0000?0000?0????0000?0000?0
for?i?=?1:（N/2）%第一級有512個蝶形運算，尋找512個蝶形運算的所需要的地址，存儲起來
????cnt?=?dec2bin（butterfly_cnt（i）m-1）;%b8b7b6b5b4b3b2b1b0---9bit??0：511，512個蝶形中計到第i個蝶形的計數(shù)值
????for?j?=?1:length（cnt）%erjinzhi?長度?9
????????wn（j）?=?cnt（j）;%wn的二進制長度也是設定的是18bit，這里賦值高9bit所以wn的大小和butte