现在的位置: 首页 > 综合 > 正文

牛顿迭代法

2018年12月20日 ⁄ 综合 ⁄ 共 537字 ⁄ 字号 评论关闭

牛顿迭代

详细介绍见维基百科 http://zh.wikipedia.org/wiki/牛顿法

 


举一个栗子


找平方根

给一个数a,求其平方根。

 

设其平方根为则有 x^2 - a = 0,设函数f(x) = x^2 - a

x0的初值尽量靠近a的平方根(因为初值的选择影响迭代的次数)

根据 f(x0)  = (x0 - x).f(x0)

--->  x = x0 -  f(x0)/f(x0)   

--->  x0 = x                        

重复①、②直到abs(x-x0)<=eps,达到精度要求即可。

代码:

a = input();

x0 = Init_Value;

While( fabs(x-x0)>eps){

       x0 = x;

       x = x0 - f(x0)/f(x0);

}

 

Eg.  a  = 2 , eps = 1e-8

X0 = 1.0

迭代过程中x的变化

1.5000000000000000

1.4166666666666667

1.4142156862745099

1.4142135623746899

1.4142135623730951

 

总的来说

F(x)函数连续,且在零点范围内收敛。

迭代的次数受x0和精度影响。

迭代次数很少相对(二分),(精度倍增)

抱歉!评论已关闭.