本帖最后由 yijinxiaoyou 于 2017-3-9 16:33 编辑

ieeloleei 发表于 2017-3-9 13:56
第四点不赞同，按照你的说法，“频谱泄露”是ARM DSP库算法造成的，让人觉得这是ARM DSP库的缺陷，然而并 ...

不知道你测的是什么信号 我测发时候同样的数组经过dsp库与MATLAB比较高次有泄露后来就没深究 后来用自己的算法没有什么问题

我有一大堆矩阵运算，前阵子算了下计算速度，有些太慢了，就分别试用了有单精度浮点运算单元的片子和有有单精度浮点运算单元的片子F767。发现用单精度片子时，库里自带的函数确实比手敲的快太多了。
感觉原因是这样的：float型乘float型，理论上需要用double型进行存储，所以可能是用float做乘除法时它自动转成了double型进行计算，就没法用FPU了，所以速度会慢，比没有FPU的片子略快一些，但跟库函数这种直接用FPU进行运算的东西比还是慢了很多。

xfdr0805 发表于 2017-3-9 19:22
根据大家回复，找到了一个快速开根号算法

这里给出实现32位无符号整数开方得到16位无符号整数的C语言代码。


//------------------------------------------------------------------------------- -
/****************************************/
/*Function: 开根号处理 */
/*入口参数：被开方数，长整型 */ /*出口参数：开方结果，整型 */
/****************************************/
unsigned int sqrt_16(unsigned long M)
{
unsigned int N, i;
unsigned long tmp, ttp; // 结果、循环计数
if (M == 0) // 被开方数，开方结果也为0
     return 0;
N = 0;
tmp = (M >> 30); // 获取最高位：B[m-1]
M <<= 2;
if (tmp > 1) // 最高位为1
{
   N ++; // 结果当前位为1，否则为默认的0
   tmp -= N;
}
for (i=15; i>0; i--) // 求剩余的15位
{
   N <<= 1; // 左移一位
   tmp <<= 2;
   tmp += (M >> 30); // 假设
   ttp = N;
   ttp = (ttp<<1)+1;
   M <<= 2;
   if (tmp >= ttp) // 假设成立
   {
     tmp -= ttp; N ++;
   }
}
return N;
}

xfdr0805 发表于 2017-3-9 19:23
这里给出实现32位无符号整数开方得到16位无符号整数的C语言代码。

/*还有一个*/
/* 来至 Quake 3 的源码 */
float CarmSqrt(float x){
    union
    {
        int intPart;
        float floatPart;
    } convertor;
    union
    {
        int intPart;
        float floatPart;
    } convertor2;
    convertor.floatPart = x;
    convertor2.floatPart = x;
    convertor.intPart = 0x1FBCF800 + (convertor.intPart >> 1);
    convertor2.intPart = 0x5f3759df - (convertor2.intPart >> 1);
    return 0.5f*(convertor.floatPart + (x * convertor2.floatPart));
}

本帖最后由 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;   
  }