main.cpp

分裂基算法

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
//#include "msp.h"
#include "complex.h"
#define pi 3.14159265359

void main()
{
	void msplfft(complex x[],int n,int isign);
	int count,power,i;
	count=256;
	power=4;
	complex *TD,*FD;
	TD=(complex *)malloc(sizeof(complex) * count);
	FD=(complex *)malloc(sizeof(complex) * count);
	for (i=0;i<count;i++)
	{
		TD[i].real=cos(2*pi*i/count)+cos(pi*i/count);
		TD[i].imag=sin(2*pi*i/count)+sin(pi*i/count);
	}
	FD=TD;
	msplfft(TD,count,-1);
	for (i=0;i<count;i++)
	{
		printf("%lf %lf\n",FD[i].real,FD[i].imag);
	}
	free(TD);
	free(FD);
}