原创 单片机开方算法的实现

long temp0;

unsigned long kaifang(unsigned long ss){

  uint8 dat_count,count;

  long temp1;

  uint8 i;

  i=17;

  dat_count=0;

  temp1=ss>>14;

  do{

if((temp1&0x3ffff)!=0x0000){//判断被开方数的最18位是否是1，是1就做右移操作，

//但是必须保证是2的偶数次，当判到有效位小于等于14位的时候进行查表，结束后

//在向右移n位乘回来，可以根据自己要求的精度来建表。

       ss>>=1;

       temp1>>=1;

       dat_count++;

    }

    else {

      if(dat_count%2==0){

       break;

      }

      else{

       ss>>=1;

       dat_count++;

      }

    }

  }while(i-->0);

  temp0=sqrt_32[ss];

  count=dat_count/2;

  for(uint8 j=0;j<count;j++){

    temp0<<=1;

  }

  return temp0;

}

该方法在单片机中使用可以大大加快程序的速度，还可以保证在比较高的精度上，缺点太耗rom了。

方法2：在程序rom比较小，但是对精度要求不高的场合可以使用

unsigned int sqrt_16(unsigned long dat){

     unsigned int NN, i;

     unsigned long tmp, ttp;    // 结果、循环计数

     if (dat == 0)                // 被开方数，开方结果也为0

         return 0;

     NN = 0;

     tmp = (dat >> 30);           // 获取最高位：B[m-1]

     dat <<= 2;

     if (tmp > 1) {              // 最高位为1

         NN ++;                  // 结果当前位为1，否则为默认的0

         tmp -= NN;

     }

     for (i=15; i>0; i--) {      // 求剩余的15位

         NN <<= 1;               // 左移一位

         tmp <<= 2;

         tmp += (dat >> 30);      // 假设

         ttp = NN;

         ttp = (ttp<<1)+1;

         dat <<= 2;

//////////////////////////////////版权所有，如需转载注明出处 鱼香茄子/////////////////////////////////////////////