FFT c++语言实现

FFT公式如下：

用图表示整个过程：

/*
 * test.cpp
 *
 *  Created on: 2013-8-18
 *      Author: zhijian
 */




#include  <stdio.h>
#include  <math.h>


#define N 128            //N 点
#define RATE 128      //抽样频率
#define LEN 7            //1<<7=128
#define PI 3.141592653


float src[N];      //输入
float AFreal[N];//实部
float AFimg[N];      //虚部


int upset(int a,int len){
      int sum = 0;
      while(len--){
            sum <<=1;
            if(a&1)sum++;
            a >>=1;
      }
      return sum;
}
/*
 * 公式X(k) = F1(k) + WnkF2(k)  k =