求平方根的几种方法

最新推荐文章于 2024-07-23 14:18:41 发布

zealerww

最新推荐文章于 2024-07-23 14:18:41 发布

阅读量7.5k

点赞数 2

分类专栏：算法文章标签： c语言 math 算法

本文链接：https://blog.csdn.net/mecephy/article/details/52916954

版权

算法专栏收录该内容

1 篇文章 0 订阅

订阅专栏

求解这里写图片描述是我们会常常遇到的，c++里一个sqrt(x)就解决了。这个值怎么得到的，倒是个有趣的问题。比较容易想到的是把这个问题转换为，接着就一个个试了……有人会说妈的智障，好吧，怪我没讲清楚，试要试得有水平，不能瞎试吧。其中，二分法就是一个试的有水平的方法，二分法简而言之就是找有可能范围中的中间的那个数来试。下面会谈到二分法、牛顿迭代法和一种神奇的方法。

二分法

以这里写图片描述为例，可能的范围是[0, 10]，取中间的值5，试下乘法口诀表五五二十五，大了。可能的范围的就缩小为了[0, 5]，接着取2.5，小了，取[2.5，5]……又有人说了，这得试到什么时候啊，这里就需要设定一个误差err，当这里写图片描述时就认为算对了。

c++代码如下：

float dichotomyMethod(float n, float err)
{

    float bigNum = n;
    float smallNum = 0;
    float midNum = 0;
    float errTemp = 0;

    do
    {
        if(errTemp > 0)
            bigNum = midNum;
        else
            smallNum = midNum;
        midNum = (bigNum + smallNum) / 2.0;
        errTemp = midNum * midNum - n;

    }while(fabs(errTemp) > err);

    return midNum;
}

牛顿迭代法

之前一直不知道这种方法，直到上了一门《数值分析》的课。当然是考试前复习的时候才直到的，O(╯□╰)O 。这种方法能够花比较少的次数接近结果，在单根附近有平方收敛。

对于这里写图片描述来说可以通过来迭代求解，直到达到精度的要求。

上面的式子很简单，几何上的意义是点Xn处的切线与x轴得到的交点就是 Xn+1。下面这张动图能比较好说明。

newton

代数上是对 f(x) 进行泰勒展开，并取前两项。这里写图片描述，由于 X0是零点，第一项则为0，即，和之前提到的迭代表达式一致。

c++代码如下：

float newtonIterationMethod(float n, float err)
{
    float x = n;
    float temp = x;
    do
    {
        temp = x;
        x = (x + n / x) / 2.0;

    }while(fabs(x * x - n) > err);

    return x;
}

对比下以上两种方法所花的时间。计算1000次总时间的对比：

二分法 625us
牛顿迭代法 193us

牛顿迭代法胜出，差距还是比较明显。

神奇方法

在博客上看到一种很神奇的方法，是来自游戏Quake III的源代码，求得是这里写图片描述。代码如下：

float InvSqrt (float x)
{
    float xhalf = 0.5f*x;
    int i = *(int*)&x;
    i = 0x5f3759df - (i>>1);
    float y = *(float*)&i;
    y = y*(1.5f - xhalf*y*y);
    return y;
}