金皓蘋, 徐秀斌
(浙江師范大學數(shù)理與信息工程學院,浙江金華 321004)
令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條件,并使用優(yōu)序列的方法分析式(7)的收斂性.下面總設X,Y為Banach空間,且F=G+H,其中H和G如式(3)定義,H為一階Fréchet可導,G為連續(xù)但不可導函數(shù).
首先給出幾個重要的引理,然后在這些引理的基礎上證明修正牛頓迭代法式(7)的半局部收斂性、誤差估界及解的唯一性.
引理1 設優(yōu)函數(shù)
式(8)中,L(u)和l(u)為非負非減可積的連續(xù)函數(shù).記
而
則方程f(t)=0有2個正根r1,r2,顯然有r1<R<r2.引理1證畢.
引理2 設迭代序列{tn}滿足
式(10)中,h(t),g(t)由式(8)定義.則{tn}單調(diào)遞增收斂到 r1.
證明 當 n=0 時,t1-t0=-(h'(t0)+[s0,t0;g])-1(h(t0)+g(t0)) < β.由 h(r1)+g(r1)=0,并根據(jù)式(8),得 β< r1,則 t0<t1<r1.現(xiàn)假設 tk<tk+1<r1對 k≤n都成立.
當k=n+1時,
則由優(yōu)函數(shù)表達式有:當 s,t∈(0,R)時,
當 t∈[0,r1)時,h(t)+g(t) >0.則由式(11)得 tn+1<tn+2.因此,由歸納假設得{tn}單調(diào)遞增.記 I(t)=t-(h'(t)+[s,t;g])-1(h(t)+g(t)),顯然 I(t)在[0,R)上是單調(diào)遞增的,則 tn+1< tn+2=I(tn+1)≤I(r1)=r1.所以{tn}有極限,不妨記為 t*,顯然 t*∈(0,r1].現(xiàn)在證明 t*=r1.對迭代式
兩邊取極限,得h(t*)+g(t*)=0,即t*=r1.綜上所述,{tn}單調(diào)遞增收斂到r1.引理2證畢.
引理3 假設下列條件成立:
證明 當n=0時,以上5個結論顯然成立.現(xiàn)設它們從0到n都成立,則根據(jù)條件2)和3)有
所以,根據(jù)Banach引理知A-1n+1存在,且結論2)成立.由迭代式(7)得
所以,根據(jù)條件3)和4)得
因此,結論3)對n+1成立.
根據(jù)結論1),2)和3)可得
則結論4)對n+1成立.
由結論4)得
所以,結論5)對n+1成立.引理3證畢.
下面給出迭代算法(7)的半局部收斂性定理,并加以證明.
顯然
所以A-1存在.故x*=y*.定理1證畢.
當G=0時,迭代式(7)就是牛頓迭代式.下面選擇特殊的L和l對迭代式進行討論.
推論1 對于常數(shù) γ >0,0≤c<1,取
此時函數(shù)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.