中科因仑“3+1”工程特种兵精英论坛

标题: 信号数据的FFT变换 [打印本页]

作者: 张衍波 时间: 2015-4-23 20:35
标题: 信号数据的FFT变换
#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;
}

欢迎光临中科因仑“3+1”工程特种兵精英论坛 (http://bbs.enlern.com/)