main.c

C语言DSPFFT算法

#include <math.h>
#define PI 3.1415926535897932384626433832795


float wav[2048];


void main(void){
	unsigned int i,j;
	for(i=0;i<2048;i++)
		wav[i]=0.0;
	for(j=1;j<100;j++)
		for(i=0;i<2048;i++)
			wav[i]+=sin(2*PI*(j+1)/1024*i);
	fft(wav,2048,0);
	fft(wav,2048,1);
	for(i=0;i<512;i++)
		wav[i]=wav[i]*wav[i]/512;
	j=i;
}