牛顿迭代法求数的平方根

最新推荐文章于 2022-09-30 18:21:19 发布

Penguinbupt

最新推荐文章于 2022-09-30 18:21:19 发布

阅读量1.6k

点赞数

分类专栏：数据结构与算法文章标签：迭代算法

本文链接：https://blog.csdn.net/u010566813/article/details/49247731

版权

数据结构与算法专栏收录该内容

26 篇文章 3 订阅

订阅专栏

牛顿法最初由艾萨克·牛顿在《流数法》（Method of Fluxions，1671年完成，在牛顿死后的1736年公开发表）。约瑟夫·拉弗森也曾于1690年在Analysis Aequationum中提出此方法。

方法说明

蓝线表示方程 f而红线表示切线. 可以看出 x _n+1比 x _n更靠近 f所要求的根 x.

首先，选择一个接近函数 $f(x)$ 零点的 $x_0$ ，计算相应的 $f(x_0)$ 和切线斜率 $f'(x_0)$ （这里 $f'$ 表示函数 $f$ 的导数）。然后我们计算穿过点 $(x_0, f(x_0))$ 并且斜率为 $f'(x_0)$ 的直线和 $x$ 轴的交点的 $x$ 坐标，也就是求如下方程的解：

f(x_0)= (x_0-x)\cdot f'(x_0)

我们将新求得的点的 $x$ 坐标命名为 $x_1$ ，通常 $x_1$ 会比 $x_0$ 更接近方程 $f(x)=0$ 的解。因此我们现在可以利用 $x_1$ 开始下一轮迭代。迭代公式可化简为如下所示：

x_{n+1} = x_n - \frac{f(x_n)}{f'(x_n)}

已经证明，如果 $f'$ 是连续的，并且待求的零点 $x$ 是孤立的，那么在零点 $x$ 周围存在一个区域，只要初始值 $x_0$ 位于这个邻近区域内，那么牛顿法必定收敛。并且，如果 $f'(x)$ 不为0, 那么牛顿法将具有平方收敛的性能. 粗略的说，这意味着每迭代一次，牛顿法结果的有效数字将增加一倍。

设： f(x)=x²-a=0,求 sqrt(a),使用牛顿迭代法。

#include "stdio.h"
#include "math.h"

int main(void)
{
            double n,y=1.0;
	    printf("请输入一个需要求其平方根的数：");
	    scanf("%lf",&n);
	 
	    // 反复代入 x(k+1) = 0.5[x(k)+n/x(k)]
	    while(fabs((1.0/2.0*(y+n/y))-y)>=1e-6)
	    {
	        y=1.0/2.0*(y+n/y);
	        printf( "y=%lf\n", y );
	    }
	    printf("平方根为%f\n",y);
	    return 0;
	}