• sqrt()平方根计算函数的实现1——二分法


    C语言标准库:

    http://www.cplusplus.com/reference/cmath/

    在一个区间中,每次拿中间数的平方来试验,如果大了,就再试左区间的中间数;如果小了,就再拿右区间的中间数来试。比如求sqrt(16)的结果,你先试(0+16)/2=8,8*8=64,64比16大,然后就向左移,试(0+8)/2=4,4*4=16刚好,得到了正确的结果sqrt(16)=4。

    实现时我第一次没有进行精度控制,导致一些不能完整开方的数运行不出来,也忘记了对非整数的判断。代码1为改正后的代码,后面代码2为第一次写的代码,贴在这里以作对比提醒。

    代码1:

    esp常在C++中用来控制迭代精度,是函数程序事先声明的常量,在微积分中是无限小值。

    float esp=0.000000;

    double sqrt(double x)

    {

           double down,up,n,last;

           down=0;

           up=x;

           n=(down+up)/2;

           if(x<=0)

           return x;

           while((n*n!=x)&&abs(last-n)>esp)

           {

                  if(n*n<x)

                  {

                         down=n;

                         last=n;

                         n=(down+up)/2;

                  }

                  else

                  {

                         up=n;

                         last=n;

                         n=(down+up)/2;

                  }

           }

           return n;

    }

    代码2:

    double sqrt(double x)

    {

           double down,up,n;

           down=0;

           up=x;

           n=(down+up)/2;

           while(n*n!=x)

           {

                  if(n*n<x)

                  {

                         down=n;

                         n=(down+up)/2;

                  }

                  else

                  {

                         up=n;

                         n=(down+up)/2;

                  }

           }

           return n;

    }

     参考:http://www.2cto.com/kf/201206/137256.html

  • 相关阅读:
    2021-06-22 总结
    【每日一题】13. 罗马数字转整数
    【每日一题】12. 整数转罗马数字
    【每日一题】1269. 停在原地的方案数
    【每日一题】1310. 子数组异或查询
    【每日一题】1734. 解码异或后的排列
    【每日一题】872. 叶子相似的树
    【每日一题】1482. 制作 m 束花所需的最少天数
    【每日一题】1723. 完成所有工作的最短时间
    【每日一题】1486. 数组异或操作
  • 原文地址:https://www.cnblogs.com/shrimp-can/p/4868141.html
Copyright © 2020-2023  润新知