实验用的头文件 MYFFT.H
作用：为帮助小虎子做实验，这个头文件提供了完整的一维与二维 FFT 算法，
我想应改是够你折腾了吧！
#include <complex> // complex<oat>
using namespace std;
typedef complex<oat> Comp; // 复数类型定义
const oat _2PI_ = 2.0f * 3.14159265f; // 常数 2PI 定义
const int MAX_N = 256; // 最大 DFT 点数
/*----*----*----*----*----*----*----*----*----*----*----*----*
FFT 算法模块接口定义
*----*----*----*----*----*----*----*----*----*----*----*----*/
///////////////////////////////////////////
// Function name : BitReverse
// Description : 二进制倒序操作
// Return type : int
// Argument : int src 待倒读的数
// Argument : int size 二进制位数
int BitReverse(int src, int size)
{
int tmp = src;
int des = 0;
for (int i=size-1; i>=0; i--)
{
des = ((tmp & 0x1) << i) | des;
tmp = tmp >> 1;
}
return des;
}
//////////////////////////////////////////////////
// Function name : Reorder
// Description : 数据二进制整序
// Return type : void
// Argument : Comp x[MAX_N] 待整序数组
// Argument : int N FFT 点数
// Argument : int M 点数的 2 的幂次
void Reorder(Comp x[MAX_N], int N, int M)
{
Comp new_x[MAX_N];
for (int i=0; i<N; i++)
