c++实现傅里叶变换的程序

void CShiyanView::FFT(double* dRe, double* dIm, int nDataNum, bool nFlag)
{
//为变量分配内存
double* ca = (double*)calloc(nDataNum, sizeof(double));
double* sa = (double*)calloc(nDataNum, sizeof(double));
if (ca == NULL || sa == NULL)
{
AfxMessageBox("分配存储空间失败！");
return;
}

//根据傅立叶变换标志进行计算
int i = 0;
double dSign;
if (nFlag)
{
dSign = 1.0;
for (i=0;i<nDataNum-1;i++)
{
dRe[i] /= nDataNum;
dIm[i] /= nDataNum;
}
}
else
dSign = -1.0;

int j = 0;
int k = 0;
for (i=0;i<nDataNum-2;i++)
{