我想给开发板一个三角波信号（3.3v，50Hz），然后信号经过ADC转换（开发板中自带的），转换后的数字信号我想再进行FFT运算，输出基波频率。我在网上找了段FFT程序，不知道如何加在ADC程序中？还有加进去后，程序中的一些参数该如何设置？请高人指点！谢谢
FFT程序如下：
/*****************fft programe*********************/
/*基二FFT算法*/
#include "math.h"
#include "stdio.h"
struct compx
{ double real;
double imag;
} compx ;
struct compx EE(struct compx b1,struct compx b2)
{
struct compx b3;
b3.real=b1.real*b2.real-b1.imag*b2.imag;
b3.imag=b1.real*b2.imag+b1.imag*b2.real;
return(b3);
}
void FFT(struct compx *xin,int N)
{
int f,m,LH,nm,i,k,j,L;
double p , ps ;
int le,B,ip;
float pi;
struct compx v,w,t;
LH=N/2;
f=N;
for(m=1;(f=f/2)!=1;m++){;}
nm=N-2;
j=N/2;
/*变址运算*/
for(i=1;i<=nm;i++)
{
if(i<j){t=xin[j];xin[j]=xin;xin=t;}
k=LH;
while(j>=k){j=j-k;k=k/2;}
j=j+k;
}
{
for(L=1;L<=m;L++)
{
le=pow(2,L);
B=le/2;
pi=3.14159;
for(j=0;j<=B-1;j++)
{
p=pow(2,m-L)*j;
ps=2*pi/N*p;
w.real=cos(ps);
w.imag=-sin(ps);
for(i=j;i<=N-1;i=i+le)
{
ip=i+B;
t=EE(xin[ip],w);
xin[ip].real=xin.real-t.real;
xin[ip].imag=xin.imag-t.imag;
xin.real=xin.real+t.real;
xin.imag=xin.imag+t.imag;
}
}
}
}
return ;
}
/*****************main programe********************/
float result[512];
struct compx s[512];
const int Num=256;/*计算点数*/
const float pi=3.14159;
long int f=12800,fs=50,fv;/*f为采样频率，fs为基频，fv为fft后计算得基频，fv结果应等于fs*/
int t;/*fft后模值最大点值序号*/
float max,module;/*最大模值及信号幅值*/
float door=0.01,module_test[512],mid_m;/*door为模值的阀值*/
int kav,fv_test[512],mid_f;
/*试验时，f,fs,Num的设定要保证下面s的取值是完整的周期点数*/
main()
{
int i,j;
for(i=0;i<Num;i++)
{
s.real=4*sin(2*pi*fs*i/f)+2.6*sin(6*pi*fs*i/f);
s.imag=0;
}
FFT(s,Num);
for(i=0;i<Num;i++)
{
//printf("%.4f",s.real);
//printf("+(%.4f)j\n",s.imag);
result=sqrt(pow(s.real,2)+pow(s.imag,2));
}
/**********************************下面为一个信号是的试验部分***************************/
/*
max=result[0];
for(i=1;i<Num;i++)
{
if(max<result)
{
max=result;
t=i;
}
}
module=max/(Num/2);//在此borland C 编译器中，要除以采样点数的一半（和matlab中一样）
printf("module=%f",module);
fv=f*t/Num;
printf("fv=%d",fv);
*/ | 
/2 
|手机版|OpenEdv-开源电子网
( 粤ICP备12000418号-1 )
发表于 2014-4-1 14:00:34
发表于 2014-4-1 22:08:45
楼主
