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 =