//在此采用的是时间抽选奇偶分解方式,所以在参加运算前首先要对时间序列进行倒序
ReverseOrder(N,pAData);
int L=1;
while(L<=M)
{
LE=(int)pow(2.0,L);
LE1=LE/2;
U.Re=1.0f;
U.Im=0.0f;
W.Re=(double)cos(PI/(1.0*LE1));//计算W算子的值
W.Im=(double)-1.0*sin(PI/(1.0*LE1));
if(fabs(W.Re)<1.0e-12)
W.Re=0.0f;
if(fabs(W.Im)<1.0e-12)
相关下载链接:
//download.csdn.net/download/cuitlc0710/4864230?utm_source=bbsseo