不可微非線性方程的修正牛頓迭代法的收斂性分析*

金皓蘋， 徐秀斌

(浙江師范大學數理與信息工程學院，浙江金華 321004)

0 引言

令F是Banach空間X到Y的非線性算子，考慮如下一般的非線性方程：

求解非線性方程(1)的近似解是一個重要的問題，因為大量的不同類型的實際問題都可歸結為對非線性方程的求解．例如，微分方程、邊界值問題、積分方程等．目前，在F是Fréchet可導的條件下，牛頓法是求解非線性方程(1)的最有效方法之一，其迭代式為(初始點x0給定)

關于牛頓迭代法收斂性的研究目前已有許多，如文獻［1-3］等．然而，當F不可導時，牛頓迭代法就不能再用來解非線性方程．對于F不可導情形的修正牛頓迭代的研究，主要歸結為當F'(x)不存在時用什么來代替的問題．諸多文獻考慮將F分解為可導部分H和不可導部分G，即

如文獻［4-5］利用

對方程(1)進行求解．

另外，文獻［6-7］采用弦割法，用差商代替F'，利用

迭代對方程(1)進行近似求解，其中初始點x0，x-1給定．因為對G的限制，式(4)一般只能保持線性收

文獻［8］將以上2種迭代方法結合起來，構造了新的迭代式

最近，文獻［9］又利用差商［yn，xn;G］代替［xn-1，xn;G］，其中 yn=λxn+(1- λ)xn-1，λ∈［0，1］，并提出了修正的牛頓變形公式

式(7)中，x-1，x0∈D已知．此迭代推廣了迭代法(6)．文獻［9］利用ω條件證明了其半局部收斂性并給出收斂定理．本文主要目的是引入L-平均Lipschitz條件，并使用優序列的方法分析式(7)的收斂性．下面總設X，Y為Banach空間，且F=G+H，其中H和G如式(3)定義，H為一階Fréchet可導，G為連續但不可導函數．

1 引理

首先給出幾個重要的引理，然后在這些引理的基礎上證明修正牛頓迭代法式(7)的半局部收斂性、誤差估界及解的唯一性．

引理1 設優函數

式(8)中，L(u)和l(u)為非負非減可積的連續函數．記

而

則方程f(t)=0有2個正根r1，r2，顯然有r1＜R＜r2．引理1證畢．

引理2 設迭代序列{tn}滿足

式(10)中，h(t)，g(t)由式(8)定義．則{tn}單調遞增收斂到 r1．

證明 當 n=0 時，t1-t0=-(h'(t0)+［s0，t0;g］)-1(h(t0)+g(t0)) ＜ β．由 h(r1)+g(r1)=0，并根據式(8)，得 β＜ r1，則 t0＜t1＜r1．現假設 tk＜tk+1＜r1對 k≤n都成立．

當k=n+1時，

則由優函數表達式有：當 s，t∈(0，R)時，

當 t∈［0，r1)時，h(t)+g(t) ＞0．則由式(11)得 tn+1＜tn+2．因此，由歸納假設得{tn}單調遞增．記 I(t)=t-(h'(t)+［s，t;g］)-1(h(t)+g(t))，顯然 I(t)在［0，R)上是單調遞增的，則 tn+1＜ tn+2=I(tn+1)≤I(r1)=r1．所以{tn}有極限，不妨記為 t*，顯然 t*∈(0，r1］．現在證明 t*=r1．對迭代式

兩邊取極限，得h(t*)+g(t*)=0，即t*=r1．綜上所述，{tn}單調遞增收斂到r1．引理2證畢．

引理3 假設下列條件成立：

證明 當n=0時，以上5個結論顯然成立．現設它們從0到n都成立，則根據條件2)和3)有

所以，根據Banach引理知A－1n+1存在，且結論2)成立．由迭代式(7)得

所以，根據條件3)和4)得

因此，結論3)對n+1成立．

根據結論1)，2)和3)可得

則結論4)對n+1成立．

由結論4)得

所以，結論5)對n+1成立．引理3證畢．

2 半局部收斂性

下面給出迭代算法(7)的半局部收斂性定理，并加以證明．

顯然

所以A-1存在．故x*=y*．定理1證畢．

3 應用

當G=0時，迭代式(7)就是牛頓迭代式．下面選擇特殊的L和l對迭代式進行討論．

推論1 對于常數 γ ＞0，0≤c＜1，取

此時函數H和G滿足

其中：

結合迭代式(10)，易得

類似可得

因此

從而

引理5證畢．

結合定理1和以上引理，可得到定理2．

并滿足

［1］Kantorvich L V，Akilov G P．Functional analysis［M］．Oxford：Pergamon Press，1982．

［2］Smale S．Newton's method estimates from data at one point［C］//Ewing R，Gross K，Martin C．The merging of disciplines：New directions in pure applied and computational mathematics．New York：Spring-Verlag．1986：185-196．

［3］Wang Xinhua．Convergence of Newton's method and inverse functions theorem in Banach space［J］．Math Comput，1999，68(225)：169-186．

［4］Argyros I K．On the solution of equations with nondifferentiable and Ptak error estimates［J］．BIT Numer Math，1990，30(4)：752-754．

［5］Zabreko P P，Nguen D F．The majorant method in the theory of Newton-Kantorovich approximations and Ptak error estimates［J］．Numer Func Anal Optimiz，1987，9(5/6)：671-684．

［6］Ren Hongmin．New sufficient convergence conditions of the secant method for nondifferentiable operators［J］．J Comput Math Appl，2006，182(2)：1255-1259．

［7］Hernández M A，Rubio M J．A uniparametric family of iterative process for solving nondifferentiable equations［J］．J Math Anal Appl，2004，275(2)：821-834．

［8］Catinas E．On some iterative methods for solving nonlinear equations［J］．Rev D'Analy Numér Theorie L'Approx，1994，23(1)：47-53．

［9］Hernández M A，Rubio M J．A modification of Newton's method for nondifferentiable equations［J］．J Comput Appl Math，2004，164/165(1)：409-417．