#include <reg52.h>
#include <absacc.h>
#include <stdio.h>
#include <math.h>
#define uint unsigned int
#define uchar unsigned char
uint x0,x1,x2,x3,x4,x5,x6;
int L,i,j,k,b,p;
uchar w[3];
uint TR,TIs,temp;

void FFT(int dataR[],int dataI[])
{
// 以下是數據掉頭功能
for(i=0;i<128;i++)
{ 
//CHANGE();
}
for(i=0;i<128;i++)
{
 dataR[i]=dataI[i]; dataI[i]=0; 
}
// FFT算法
for(L=1;L<=7;L++) 						// 第一層循環
{
b=1; i=L-1;
while(i>0) 
{
b=b*2; i--;
} 
for(j=0;j<=b-1;j++)						 // 第二層循環
{ 
p=1; i=7-L;
while(i>0)					
{p=p*2; i--;}
p=p*j;
for(k=j;k<128;k=k+2*b) 				// 第三層循環
{ 
FFT(dataR,dataI);
}
}
} 
for(i=0;i<32;i++) 							// 32次以下的諧波分析
{
w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i]);
w[i]=w[i]/64;
}
w[0]=w[0]/2;
}