创龙官方例程单做ifft为啥算的不对

901 频率

[tr]/****************************************************************************/
/*                                                                         */
/*             基2 时间抽取快速傅里叶变换 / 快速傅里叶逆变换             */
/*                                                                         */
/*             2014年08月25日                                           */
/*                                                                         */
/****************************************************************************/
#include                // C 语言标准输入输出函数库
#include                // C 数学函数库

#include "mathlib.h"             // DSP 数学函数库
#include "dsplib.h"                // DSP 函数库

/****************************************************************************/
/*                                                                         */
/*             宏定义                                                    */
/*                                                                         */
/****************************************************************************/
// 软件断点
#define SW_BREAKPOINT    asm(" SWBP 0 ");

// 采用点数
#define Tn  16
// 采样频率
#define Fs  256.0

// π 及浮点数极小值
#define PI             3.14159
#define F_TOL          (1e-06)

/****************************************************************************/
/*                                                                         */
/*             全局变量                                                 */
/*                                                                         */
/****************************************************************************/
// 快速傅里叶变换测试
#pragma DATA_ALIGN(Input, 8);
float Input[Tn];

#pragma DATA_ALIGN(CFFT, 8);
float CFFT[2*Tn]={0,0,-9,9,9,-9,-9,-9,9,9,-9,9,9,-9,-9,-9,0,0,-9,9,9,9,-9,-9,9,-9,-9,9,9,9,-9,-9};

#pragma DATA_ALIGN(CFFT_In, 8);
float CFFT_In[2*Tn];

#pragma DATA_ALIGN(CFFT_Out, 8);
float CFFT_Out[2*Tn];

#pragma DATA_ALIGN(Cw, 8);
float Cw[2*Tn];

#pragma DATA_ALIGN(Cmo, 8);
float Cmo[Tn];

/****************************************************************************/
/*                                                                         */
/*             函数声明                                                 */
/*                                                                         */
/****************************************************************************/
void bit_rev(float* x, int n);
void gen_w_r2(float* w, int n);

void FFTTest(void);

/****************************************************************************/
/*                                                                         */
/*             主函数                                                    */
/*                                                                         */
/****************************************************************************/
int main(void)
{
      // FFT 测试
      FFTTest();

      // 断点
      SW_BREAKPOINT;
}

/****************************************************************************/
/*                                                                         */
/*             快速傅里叶变换测试                                        */
/*                                                                         */
/****************************************************************************/
// 比特位反转
void bit_rev(float* x, int n)
{
int i,j,k;
float rtemp,itemp;
j=0;

for(i=1;i<(n-1);i++)
{
      k=n>>1;

      while(k<=j)
      {
         j-=k;
         k>>=1;
      }

      j+=k;

      if(i       {
         rtemp=x[j*2];
         x[j*2]=x[i*2];
         x[i*2]=rtemp;
         itemp=x[j*2+1];
         x[j*2+1]=x[i*2+1];
         x[i*2+1]=itemp;
      }
}
}

// 产生旋转因子
void gen_w_r2(float* w, int n)
{
int i;

float e=PI*2.0/n;

for(i=0;i<(n>>1);i++)
{
      w[2*i]=cossp(i*e);
      w[2*i+1]=sinsp(i*e);
}
}

// 快速傅里叶变换
void FFTTest(void)
{


      // 保留一份 FFT 输入副本
memcpy(CFFT_In, CFFT, 2*Tn*sizeof(float));

gen_w_r2(Cw, Tn);

      bit_rev(CFFT_In, Tn);

      // IFFT 计算
      DSPF_sp_icfftr2_dif(CFFT_In, Cw, Tn);
      memcpy(CFFT_Out, CFFT_In, 2*Tn*sizeof(float));

      for(i=0;i       {
            CFFT_Out=CFFT_Out/Tn;
      }

}

[/tr]