Newton迭代法求解非线性方程

1、Newton迭代法求解非线性方程一、Newton迭代法概述构造迭代函数的一条重要途径是用近似方程来代替原方程去求根。因此，如 果能将非线性方程f(x)= 0用线性方程去代替，那么，求近似根问题就很容易解 决，而且十分方便。牛顿(Newton)法就是一种将非线性方程线化的一种方法。设兀是方程f(x) = 0的一个近似根，把如果/(x)在忑处作一阶Taylor展开,即:(1-1)(1-2)f(x)f(xk) + f(xk)(x-xk)于是我们得到如下近似方程：f(xk) + f(xk)(x-xk) = O设广(母)工0,则方程的解为:f(Xk)f(lk)(1-3)取刁作为原方程的新近似根X如，即令

2、:(1-4)上式称为牛顿迭代格式。用牛顿迭代格式求方程的根的方法就称为牛顿迭代法， 简称牛顿法。牛顿法具有明显的几何意义。方程：y = f(xk) + f(xk)(x-xk)(1-5)是曲线y = f(x)上点(Xk，f(xQ)处的切线方程。迭代格式(1-4)就是用切线式(1-5)的零点来代替曲线的零点。正因为如此，牛顿法也称为切线法。牛顿迭代法对单根至少是二阶局部收敛的，而对于重根是一阶局部收敛的。一般来说，牛顿法对初值X。的要求较髙，初值足够靠近x*时才能保证收敛。若 要保证初值在较大范围内收敛，则需对f(x)加一些条件。如果所加的条件不满足,而导致牛顿法不收敛时，则需对牛顿法作一些改时，

3、即可以釆用下面的迭代格式:x k+i = x k 九x J k = 0丄2, fg)9(1-6)上式中，0九1, 下山法。称为下山因子。因此，用这种方法求方程的根，也称为牛顿牛顿法对单根收敛速度快，但每迭代一次，除需计算f(xj之外，还要计算f(xk)的值。如果f(x)比较复杂，计算f(xk)的工作量就可能比较大。为了避免 计算导数值，我们可用差商来代替导数。通常用如下几种方法：1. 割线法如果用(乞)二门4 )代替f，(xQ,则得到割线法的迭代格式为：Xk-XkX k+1 = X kdf(X k )f(Xk)-f(Xk_J(1-7)2. 拟牛顿法如果用M 土代替仕小，则得到拟牛顿法的迭代格式

4、为: f(xQ=f2(xk)Xk+1Xk f(Xk)-f(Xk-f(Xk_,)(1-8)3. Steffenson 法如果用丄 代替仕,则得到拟牛顿法的迭代格式为： f(Xk)=f2(xk)Xk+1_Xk f(Xk+f(Xk)-f(Xk)(1-9)二、算法分析1. 割线法k-1割线法的迭代公式为:X“1 = Xkf(xQ ,k=0, 1,2,f(xQ-f(Xk_J 割线法是超线性收敛，其程序流程图为2. 拟牛顿法牛顿拟迭代法迭代公式为:fxk)k+lf(xk)-f(xk-f(xk_1)(1)对单根条件下，牛顿拟迭代法平方收敛，牛顿拟迭代法程序框图如下所 示：(2)对重根条件下，此时迭代公式修改

5、为:=Xk _mf(xk)-f(xk-f(xk_1)为根的重数此时，牛顿迭代法至少平方收敛。3. Steffenson 法Steffenson迭代法程序流程图与牛顿拟迭代法类似。三、牛顿法的程序给定初值久，用牛顿法格式Pk=P単丄，k = l,2,.,求解非线性方程 f (Pk-i)f(x) = o o *function Lpl,err,k,y = newton(f1041Tdf1041,pO,delta,maxi) % fl041是非线性函数。% dflO41 是 fl041 的微商。% p0是初始值。% delta是给定允许误差。% maxi是迭代的最大次数。% pl是牛顿法求得的方程的

6、近似解。% err是p0的误差估计。% k是迭代次数。% y = f(pl)p0, feval(ffl041p0)for k = 1:max1pl = p0 - feval(f1041, p0)/feval(dfl041, pO):err = abs (pl-pO)；pO = pl;pl, err, k, y 二 feval(r f1041F, pl)if (err prog1041计算结果如下：p0 =k =2ans =y =pl =pi =err =err =k =1k =3y =y =pi =pl =err =err =k =1k =3y =y =pl =pi =err =err =k

7、=2k =4y =y =pl =pi =err =err =k =4y =pi =err =k =5y =pi =err =k =5y =pi =err =k =6y =pi =err =k =6y =pi =err =k =7y =pi =err =k =7y =pi =err =y =pi = err =k =8y =Pl =err =k =9y =Pl = err =k =9y =pl = err =k =10y =Pl =err =k =10y =Pl =err =k =11y =pi =err =k =11y =k =8pl =y =err =pl =k =12err =y =k =

8、15pl =y =err =pi =k =12err =y =k =16pl =y =err =pi =k =13err =y =k =16pl =y =err =pl =k =13err =y =k =17pl =y =err =pi =k =14err =y =k =17pl =y =err =pi =k =14err =y =k =18pl =y =err =ans =k =15这说明,经过18次迭代得到满足精度要求的值。以下是程序运行截图:r uriii-xr i一 44 f “ O y mOM IT i XXZaA t Vm wW ftxx-dlow ItXlO兀,cF”V uvixx I t ys 1 O O” t IXX IT JO7 iU1 :rxhWX u. euxoo. rubryPCS mCImO *7Afllfilft a700 red- O GfOtSGOdOn _c.c TrcrzcnTrxO-O . moOUTQQO口 e vo .j/=dw 3 a z d3 *r 3Bt