cmpl.h

几个FFT算法

#if !defined(AFX_COMPLEX_H__667E486D_56ED_483D_B45D_68A191910917__INCLUDED_)
#define AFX_COMPLEX_H__667E486D_56ED_483D_B45D_68A191910917__INCLUDED_

#define PI_2 6.283185307179586476925286766559
#define PI   3.1415926535897932384626433832795

typedef unsigned long DWORD;

typedef int BOOL;

#define FALSE 0
#define TRUE  1

typedef struct 
{
	double Re;
	double Im;
}CMPL; 

//c_add(c,a,b),c=a+b

#define CMPL_ADD(c,a,b) c.Re=a.Re + b.Re; c.Im=a.Im+ b.Im;
#define CMPL_SUB(c,a,b) c.Re=a.Re-  b.Re; c.Im=a.Im- b.Im;
#define CMPL_MUL(c,a,b) c.Re=(a.Re*b.Re)-(a.Im*b.Im); c.Im=(a.Re*b.Im)+(a.Im*b.Re); 

#endif