stm32上跑FFT 耗时50ms正常么

fqxwyf

做的是512点的 不是us级的么 怎么我测出来都快50ms了，频率和幅值算得还算正常，难道要用库，求使用库做FFT的例程

void test (void) 
{
 u16 i=0;

  for(i=0;i<512;i++)
  {
   x11.real=(ADCDataTab[i+fft_j])*1.0;
   x11.img=0;
   x22.real=0;
   x22.img=0;  
  }

  FFT1();
}
int FFT1(void)
{
 complex c;
 x1=x11;
 x2=x22; 
 //蝶形运算
 for(kk=0;kk<power;kk++)
 {
  for(jj=0;jj<(1<<kk);jj++)
  {
   bfsize=1<<(power-kk);
   for(ii=0;ii<bfsize/2;ii++)
   { 
    pp=jj*bfsize;
    x2[ii+pp].real=x1[ii+pp].real+x1[ii+pp+bfsize/2].real;
    x2[ii+pp].img=x1[ii+pp].img+x1[ii+pp+bfsize/2].img;
    c.real=x1[ii+pp].real-x1[ii+pp+bfsize/2].real;
    c.img =x1[ii+pp].img-x1[ii+pp+bfsize/2].img;
    x2[ii+pp+bfsize/2].real=c.real*W[ii*(1<<kk)].real-c.img*W[ii*(1<<kk)].img;
    x2[ii+pp+bfsize/2].img =c.real*W[ii*(1<<kk)].img+c.img*W[ii*(1<<kk)].real;
   }
  }
  x=x1;
  x1=x2;
  x2=x;
 }
 for(jj=0;jj<count;jj++)
 {
  pp=0;
  for(ii=0;ii<power;ii++)
  {
   if (jj&(1<<ii)) pp+=1<<(power-ii-1);
  }
  x2[jj]=x1[pp];
 }
    for (ii=0;ii<256;ii++)
 {
  mag[ii]=sqrt(x2[ii].real*x2[ii].real+x2[ii].img*x2[ii].img);
  //mag[ii+128]=atan2(x2[ii].real,x2[ii].img)*180/pi;
 }
 return ii;
}

XINSI

32的那个片子？

正点原子

代码要优化...

fqxwyf

回复【2楼】XINSI:
---------------------------------
103VET6

fqxwyf

回复【3楼】正点原子:
---------------------------------
我用了库，256点还是ms，5ms

void  dsp_asm_init()
{
  u16 i=0;
  for(i=0;i<NPT;i++)
  {
lBUFIN = ((s16)ADCDataTab[i+fft_j])<<16 ; //转存DMA的AD采集
  }
}
void dsp_asm_powerMag(void)
{
  s16 lX,lY;
  u32 i;
  for(i=0;i<NPT/2;i++)
  {
    lX  = (lBUFOUT ><< 16) >> 16;
    lY  = (lBUFOUT >> 16);
    {
    float X    = NPT * ((float)lX) /32768;
    float Y    = NPT * ((float)lY) /32768;
    float Mag = sqrt(X*X + Y*Y)/NPT;
    lBUFMAG    = (u32)(Mag * 65536);
    }
  }
}

/*
*********************************************************************************************************
*   Test FFT,calculate powermag
*   进行FFT变换，并计算各次谐波幅值
*********************************************************************************************************
*/
void  dsp_asm_test()
{
// 根据预定义选择合适的FFT函数
#ifdef NPT_64
  cr4_fft_64_stm32(lBUFOUT, lBUFIN, NPT);
#endif

#ifdef NPT_256
  cr4_fft_256_stm32(lBUFOUT, lBUFIN, NPT);
#endif

#ifdef NPT_1024
  cr4_fft_1024_stm32(lBUFOUT, lBUFIN, NPT);
#endif

  // 计算幅值
  dsp_asm_powerMag();

}
主函数
  while(1)
  {  
 
PD8_H;
delay_ms(5);
PD8_L;

dsp_asm_init();
dsp_asm_test();
PD8_H;
delay_ms(10); 
PD8_L;
delay_ms(10);
}

正点原子

回复【5楼】fqxwyf:
---------------------------------
快很多了...

Samouyang

我用STM32F407VG ，跑4096点的FFT，用的DSP的函数，怎么耗时要4.9秒呢，哪位大神可以解释一下吗？