求高手帮忙解释一下这个程序

lhz815

　　2. 2 　请求高手给解释一下这个程序传统傅氏算法的仿真
　　在用MATLAB 进行仿真模拟时,要注意N 数目的选择,
　　这将对仿真效果产生影响. 仿真程序的流程图如图1.
　　仿真程序的部分代码如下:
　　%FFT
　　global step %运算步数
　　global bitbutterfly %每个蝴蝶结中所包含的点数
　　global frequency %每步运算中蝴蝶结的个数
　　for step = 1 :5
　　for bitbutterfly = 1 :2^(step - 1)
　　i = (2^(5 - step) ) 3 (bitbutterfly - 1) ;RW = cos(2 3 pi 3
　　i/ 32) ; % W的实部
　　IW= ( - 1) 3 sin(2 3 pi 3 i/ 32) ; % W的虚部
　　for frequency = 1 :2^(5 - step)
　　temp = (frequency - 1) 3 (2^step) + bitbutterfly ;
　　TR = dataR(temp) ;TI = dataI(temp) ;
　　dataR(temp) = dataR(temp) + RW3 dataR(temp + 2^( step -
　　1) ) - IW3 dataI(temp + 2^(step - 1) ) ;
　　dataI(temp) = dataI (temp) + RW 3 dataI (temp + 2^ ( step -
　　1) ) + IW3 dataR(temp + 2^(step - 1) ) ;
　　temp1 = dataR(temp + 2^(step - 1) ) ;
　　dataR(temp + 2^(step - 1) ) = TR - (RW3 dataR(temp + 2^
　　(step - 1) ) - IW3 dataI(temp + 2^(step - 1) ) ) ;
　　dataI(temp + 2^ ( step - 1) ) = TI - (RW 3 dataI (temp + 2^
　　(step - 1) ) + IW3 temp1) ;
　　end
　　end
　　end
　　[ 本帖最后由 zhangnan3509 于 2008-4-17 11:21 编辑 ]