求方程实根的一种迭代方法.docVIP

PAGE PAGE 29 关于求方程实根的具有大范围收敛的迭代方法，近年来受到广泛的重视，出现了许多有 意义的成果.但在这些方法中，有的需要计算函数的二阶导数，有的对函数的性质要求 高.本文给出了一个不需要计算函数二阶导数的迭代公式，它对函数性质要求低，并论证了 该公式的大范围收敛性及收敛速度?最后给出了数值例子? 　　1 迭代公式的建立 　　设 /(a) 6 cz[_ajb] 　　由台劳展开式得： 　　/(?r) = f{x~ ) + f (x~ ) (x - X7 ) +- x~ )2/2，a i b(1) 　　或/(工）=f{Xn ) + /，- x广）+ 尸(f+)(?r - X； )2/2,a c+ 6 　　其中，6 (?6),且 xt，令： 　　giCr) = ) + f (x^)(j: - x；f ) + D+ (x - xtY/2(3) 　　g*2(x) - f{x~ ) + /; {x； ) (x - x;) + D~ (x - x~ y/2(4) 　　\ 一 M± f (xt) 0 , 　　其中，Z)士I广(6士）|，比较（1)，（4)和（2)，（3)式有： 　　[M±/(xr) 0 　　当/?)0时， 当{fix) - gx(x) ^ 0(5) 　　\f(x) - g2(x) ^ 0/Cr?) 0 时， 　　(6) 　　{/(x) - ^i(^) 0,(7) 　　I/O) - g2{x) 0(8) 　　因此，用方程幻U) = 0大于xt的根作为/Or) = 0大于W的最小根的近似值，用方 　　收稿日期：1992-10-09程a(x = 0小于■r；的根?r+i作为/(x) = 0小于工:的最大根的近似值.由一元二次方 程求根方法可得： 　　x:+i = xt + 2f(.xt)/[- /?) + (sign/(x.+ ) VA.+ ](9) 　　x~+l = x~ + 2/(x7 )/[- f (.x~ ) - (sign/Cr;)) VaF](10) 　　其中，△?士 =[尸 U,士)]2 - 2/U?士 )D士 _ 　　2 M±的确定 　　定理l设/u) e ^0，6]，且满足条件： 　　(1)/(a)/(6)0； 　　(2)尸(x)在|，6]上不变号； 　　(3)尸（工）#0，:c 6 [a,6]. 　　则 I 尸(0丨 2 丨一尸(6) + 尸(a)|/|5 -a|，x* 　　其中，5=min{a- (b - a)/(a)/[/(A) -/(a)] ,a -/(a)//(a)} .x* 为/(x)在U，6)上 的零点. 　　证明由台劳展开式得： 　　f(x) = f(a) + f(a)U - a)+尸-a)2/2 = 0,x* 　　/(a) = /U*)+/,(7)(a-^), 　　则尸(f) =一 2[/(?) + 尸(a)]/(z* - a)(11) 　　由文献[4]可知，/(x)在[a,6]上有四种不同情况，但只需对其中的一种情况证明即可 .故设在定理的条件下，函数/U)在卩，6]上有下列性质： 　　(1)/(a) 0, /(? 0； 　　(2)尸(x)0，x€ 0，6]，即/U)单调下降; 　　(3)尸（6)尸（2) 0. 　　设 a = a ~ (b - a)/(a)/[/(6) - /(a)]，即 a = a + /(a) {b - a)/[/(a) - /(6)]. 　　由于0 /(a)/[/(a) - /(A)] 1，所以a ? 6.由于fix)在0，A]上是单调下降的， 　　所以只需证/(?) 0即可得《Z. 　　令 p、x) = /⑷ + [/(*) - /(a)](x - a)Kb 一 a)，显然《 是 p{x)的零点? 　　由插值法可得： 　　fix) = /(a) + [/(6) - /(a)](x - a)/(b - a) 　　+ /(f)(x - a)(x - b)/2, aii ib 则 /U) -/?Cr)=尸(f)U - a)(x - 6)/2，/(?) = fX?) (a - a) (a - b)/2 0. 　　故有 a ?:?:?. 　　令I? = a -/(?)//(a)，则由台劳展开式得： 　　P = a- [/U* +f(S)(a~ r-)]//(a) 　　=a + f (f)(a:* - a)//(a、_，a c ^ x* 　　由于 / (f) f (a) 0，所以 f (a) 1.故有 f .因此，取S = min{a，仍，必有 axx*.另一方面，由于尸(6)/(7)/G)0，故 I -/^幻+尸^糾彡丨 一/(7) + /(a) I，将上述结论代入(11)式可得： 　　I尸(f)| 2丨-fib) + f