C++实现离散傅里叶变换的快速算法(FFT)

前言：近期作为小白在学习FFT，故写此文，一来提高对FFT的领悟程度，而来求大神指导。



三、 FFT的C++代码设计的
1、 步骤
（1） 将原序列变化为序号为逆序的序列；
（2） 设总点数N为2的m次方，那么经过m级分解后就可以进行2点DFT运算。
（3） 从左到右找规律，对第i级分解（i=0、1、2…m），有规律：①要进行2^{(m-1-i)组蝶形运算②每组蝶形运算依次进行2}i次运算③2组蝶形运算之间的序号间隔为2^(i+1) ,④每组蝶形运算内部2个数之间的序号间隔为2^i。
（4） 由于蝶形运算是复数计算。因此可以选择将最终结果表达为模和辐角的形式，也可以表达为实部和虚部的形式，通常选择前者。本文采取的是输出2个数组，1数组个存储幅值（模），1个数组存储辐角。取出即将进行蝶形计算的2个数，分解成实部和虚部的形式，计算旋转因子。根据蝶形运算规则完成计算，再将实部和虚部形式转换成模和辐角的形式

/*
1、amp：幅值存储空间；phase：相位存储空间,
2、m：存储长度为2^m次（m >=1 ），在函数外保证
*/
void FFT(double *amp , double *phase , int m )
{
   
	double	pi = 3.14159265358979323846 ;  //定义π的值

	//将数据序列按逆序排列
	double	temp ;
	int		i_temp ;	//中间临时数据
	int		nixu ;		//数据的逆序下标

	// 1<<m 即 2^m
	for (int i = 0 ; i < (1<<m) ; i++)
	{
   
		i_temp	= i ; // i是顺序下标
		nixu	= 0 ; // 逆序下标
		for (int j = 0; j < m; j++)
		{
   
			nixu	= nixu << 1 ;				// 逆序下标向左移一位&#