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

 

原文链接: https://www.cnblogs.com/williamjie/p/11131314.html

欢迎关注

微信关注下方公众号,第一时间获取干货硬货;公众号内回复【pdf】免费获取数百本计算机经典书籍

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

原创文章受到原创版权保护。转载请注明出处:https://www.ccppcoding.com/archives/298262

非原创文章文中已经注明原地址,如有侵权,联系删除

关注公众号【高性能架构探索】,第一时间获取最新文章

转载文章受原作者版权保护。转载请注明原作者出处!

(0)
上一篇 2023年2月15日 下午7:24
下一篇 2023年2月15日 下午7:24

相关推荐