STM32

STM32 FFT之后计算模值比较慢的问题

2019-08-17 07:22发布

站内问答 / STM32/STM8
10653 15 1179
本帖最后由 xfdr0805 于 2017-3-8 22:29 编辑

STM32F103c8t6用DSP库进行 256点 FFT运算，仿真发现FFT运算倒挺快的，跟官方速度差不多，但是计算幅值的时间，时间特别长，有差不多40多ms，造成3264点阵刷新比较慢，速度差不多20帧吧，反正也过的去，这里想请教一下大家，这种情况求模运算比FFT运算慢那么多是怎么回事？有没有其它办法提高速度呢，请做过这方面的朋友给指点一下，谢谢了

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;
                if(i==0)
                {
                        lBUFMAG    = (u32)(Mag* 32768*31/256);
                }
                else
                {
                       lBUFMAG    = (u32)(Mag* 65536*31/150);
                }
    }
  }
}


友情提示: 此问题已得到解决，问题已经关闭,关闭后问题禁止继续编辑，回答。
15条回答
xfdr0805
2019-08-18 09:58
本帖最后由 xfdr0805 于 2017-3-9 21:09 编辑
xfdr0805 发表于 2017-3-9 19:22
根据大家回复，找到了一个快速开根号算法

  还有一个带神秘常数的
float   Q_rsqrt(   float   number   )   
  {   
      long   i;   
      float   x2,   y;   
      const   float   threehalfs   =   1.5F;   
   
      x2   =   number   *   0.5F;   
      y     =   number;   
      i     =   *   (   long   *   )   &y;     //   evil   floating   point   bit   level   hacking   
      i     =   0x5f3759df   -   (   i   >>   1   );   //   what   the   func?   
      y     =   *   (   float   *   )   &i;   
      y     =   y   *   (   threehalfs   -   (   x2   *   y   *   y   )   );   //   1st   iteration   
      //   y     =   y   *   (   threehalfs   -   (   x2   *   y   *   y   )   );   //   2nd   iteration,   this   can   be   removed   
   
      return   y;   
  }

一周热门 更多>