sqrt()*方根计算函数的实现2——牛顿迭代法
牛顿迭代法:
牛顿迭代法又称为牛顿-拉夫逊方法,它是牛顿在17世纪提出的一种在实数域和复数域上*似求解方程的方法。多数方程不存在求根公式,因此求精确根非常困难,甚至不可能,从而寻找方程的*似根就显得特别重要。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x) = 0的根。牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附*具有*方收敛,而且该法还可以用来求方程的重根、复根,此时线性收敛,但是可通过一些方法变成超线性收敛。另外该方法广泛用于计算机编程中。
牛顿迭代公式:
设r是f(x) = 0的根,选取x0作为r初始*似值,过点(x0,f(x0))做曲线y = f(x)的切线L,L的方程为y = f(x0)+f'(x0)(x-x0),求出L与x轴交点的横坐标 x1 = x0-f(x0)/f'(x0),称x1为r的一次*似值。过点(x1,f(x1))做曲线y = f(x)的切线,并求该切线与x轴交点的横坐标 x2 = x1-f(x1)/f'(x1),称x2为r的二次*似值。重复以上过程,得r的*似值序列,其中x(n+1)=x(n)-f(x(n))/f'(x(n)),称为r的n+1次*似值,上式称为牛顿迭代公式。
牛顿迭代法求*方根:
求*方根在牛顿迭代公式中,f(x)=x^2-a,则f'(x)=2x。以上的牛顿迭代公式变为:x(n+1)=x(n)-(x(n)^2-a)/2x,即(x(n)+a/x(n))/2。我们随便猜一个数r,假设r是f(x)=0的根,经过几次牛顿迭代公式后(以上的公式)所得到的x值即是f(x)=0的根或者其非常精确的*似值。
例如:我想求根号2等于多少。假如我猜测的结果为4,虽然错的离谱,但你可以看到使用牛顿迭代法后这个值很快就趋*于根号2了:
( 4 + 2/4 ) / 2 = 2.25
( 2.25 + 2/2.25 ) / 2 = 1.56944..
( 1.56944..+ 2/1.56944..) / 2 = 1.42189..
( 1.42189..+ 2/1.42189..) / 2 = 1.41423..
盗图一张以作说明,图片来自http://www.2cto.com/kf/201206/137256.html。
..
程序实现:
#include<iostream>
#include<math.h>
using namespace std;
double sqrtNT(double a,double b)
{
double x,last;
x=b;
if(a<=0)
{
return a;
}
while(x*x!=a&&(abs(last-x)>0.0000001))
{
last=x;
x=(x+a/x)/2;
}
return x;
}
int main()
{
cout<<sqrtNT(93273,5)<<endl;
return 0;
}