分享FFT算法的应用程序及讨论

1万 · 2015-12-20 11:07

关于数据处理过程中用到的傅里叶变换的，大家一直都听说的，这里我是仅仅将我得到的数据用FFT变换处理成频率并输出的，，大家帮着看下还有啥需要改进的。。

1万 · 2015-12-20 11:08

我直接贴代码的了啊
vp=CalcVibraFreqAmp();    //频率求解
VibraParam CalcVibraFreqAmp(void)
{
float tmp=0,maxAcc=0; //
u16 i=0;  //
u8 max_i=0;

for(i=0;i<N;i++)
{
x.real=(float)accXYZ;  //*3.3/4096;;    //
x.img=0;
}
fft();

maxAcc =0;
   for(i=2;i<200;i++)// for(i=2;i<N/2;i++) //Fb:3.125~187.5Hz @ FS==400Hz
   {
tmp=x.real*x.real+x.img*x.img;
if(tmp>maxAcc)
{
   maxAcc=tmp;
// if(i ==50)       //ÆÁ±Î50Hz
// {;}
// else
   max_i =i;
}
   }
   i =(float)FS*max_i/N; //
   vp.VibraFrequency=i*100;
//ifft();
   vp.VibraAmplitude=0; //
   return vp;
}

1万 · 2015-12-20 11:09

这个程序中关键点还是fft()函数的，我继续给分享的了
void fft(void)
{
int i=0,j=0,k=0,l=0;
    complex up,down,product;
    change();
    for(i=0;i<log(size_x)/log(2);i++)/*Ò»¼¶µûÐÎÔËËã*/
    {
    l=1<<i;
    for(j=0;j<size_x;j+=2*l)/*Ò»×éµûÐÎÔËËã*/
    {
    for(k=0;k<l;k++)/*Ò»¸öµûÐÎÔËËã*/
    {
mul(x[j+k+l],W[size_x*k/2/l],&product);
add(x[j+k],product,&up);
sub(x[j+k],product,&down);
x[j+k]=up;
x[j+k+l]=down;
      }
   }
   }
}

1万 · 2015-12-20 11:10

这个是上边程序中的子函数的
void add(complex a,complex b,complex *c)
{
c->real=a.real+b.real;
c->img=a.img+b.img;
}

void mul(complex a,complex b,complex *c)
{
c->real=a.real*b.real - a.img*b.img;
c->img=a.real*b.img + a.img*b.real;
}

void sub(complex a,complex b,complex *c)
{
c->real=a.real-b.real;
c->img=a.img-b.img;
}

void divi(complex a,complex b,complex *c)
{
c->real=( a.real*b.real+a.img*b.img )/( b.real*b.real+b.img*b.img);
c->img=( a.img*b.real-a.real*b.img)/(b.real*b.real+b.img*b.img);
}