DFT离散傅立叶变换C++实现

DFT的执行效率是O(n^2)，FFT为O(log2n)，但是它对点数没有限制。

/*
     离散傅立叶算法V1.0
	   含有：DFT,IDFT
         made by xyt
		  2015/7/5
*/
#ifndef _DFT_H
#define _DFT_H
#include
#include
using namespace std;
#define PI 3.14159265354
struct complex{
	double r,i;
};
complex multi(complex a,complex b){
	complex tmp;
	tmp.r=a.r*b.r-a.i*b.i;
	tmp.i=a.r*b.i+a.i*b.r;
	return tmp;
}
int fi(double in){
	if((in-(int)in)>0.5) return (int)in+1;
	else return (int)in;
}
/* 离散傅立叶正变换,输出[][2]数组实部在前,采样容量n可以任意正整数 */
void DFT(int *in,double **out,const int &n)
{
	int i,j;
	complex **W=new complex*[n];
	for(i=0;i