求高手帮忙解释一下这个程序
本帖最后由 wdhd 于 2016-9-9 13:55 编辑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 编辑 ]
页:
[1]