求重根的一類三階迭代法*

潘云蘭

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

?

求重根的一類三階迭代法*

潘云蘭

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

給出了求非線性方程重根的一類迭代法,證明了這類方法的三階收斂性,獲得了迭代誤差,指出了這個類的廣泛性,即它包含了一些已知的方法.通過數值例子與一些已知方法進行比較,說明了新方法的有效性，即在某些情形下,新方法比一些已知方法收斂快,且在其他方法發散的情況下新方法還是以很快的速度收斂.

迭代法;收斂階;三階收斂性;重根

0 引 言

近似求解非線性方程f(x)=0在數學、物理和其他科學領域中具有非常重要的應用.除少數特例外,一般都通過迭代法求解這類問題,即從某個或某幾個初始點開始,產生一個逼近方程f(x)=0解的迭代序列.Newton法是最著名的方法之一，其公式為

這個方法具有2階收斂性，但它無法用于求重根.設非線性方程f(x)=0有m重根α,即f(j)(α)=0,j=0,1,…,m-1，但f(m)(α)≠0.當重數m事先不知時，由于f(x)的重根必為函數u(x)=f(x)/f′(x)的單根,故可通過對函數u(x)應用Newton法(1)來求f(x)的重根.但此時該方法只有1階收斂性,且需要求f(x)的2階導數，工作量增大，效率指數降低.當重數m事先知道時，可通過修改Newton法(1)來求α,即

這個方法也具有2階收斂性.

為提高求重根迭代法的收斂階，一些學者紛紛提出了新的方法.如文獻[1]提出3階Halley方法：

文獻[2]提出3階Euler-Chebyshev方法：

文獻[3]給出的3階方法：

文獻[4]給出的3階方法：

文獻[5]通過對式(4)和式(5)進行組合，給出了下面的一類3階方法：

式(7)中,θ∈R.這個方法類包含了下面2個新方法：

2)取θ=-1,得

文獻[6]也給出了一類3階方法：

式(10)中：C,D∈R;

受以上工作的啟發，本文首先引入一類更為廣泛的用于解非線性方程重根的3階迭代法,它包含了以上提到過的所有方法;然后，分析了它的收斂性，導出了它的誤差估計;最后,筆者給出了這個新類的幾種特殊情形,并利用數值例子,對引言中提到的方法和本文的新方法進行了比較，比較結果顯示:在大多數情形下,本文的新方法具有明顯的優勢.

1 一類新方法及其收斂性

為得到更多的解非線性方程重根的方法，引入如下更具一般性的方法類：

式(14)中：δ,β,γ,ρ,η,λ∈R;

下面定理給出由式(14)定義的方法類的收斂階和誤差估計:

定理1設D?R是開區間，f:D→R具有m+3階導數，α∈D是f(x)的一個m重根，x0充分靠近α.若

(ρ+λ+η)m2-(η+2ρ)m+ρ≠0,則由式(14)定義的方法類是3階收斂的，其誤差可表示為

式(18)中：

N是量m,γ,ρ,η,λ,c1,c2,s1,s2,d1,d2的多項式，而

證明 記en=xn-α.把f(x),f′(x)和f″(x)在α處Taylor展開,得

從而

(20)

(21)

式(20)和式(21)中：

所以

(23)

由式(20)、式(22)和式(23)，可把式(14)寫成

式(24)中：

(25)

5δλm2-βρm2-6δm3λ+βηm2+δm4ρ+δm4η+7δηm2-2βm3η-4δρm-3δηm+βρm+

2γm3ρ-5δηm3+βm4λ+βm4ρ+βm4η-3γm2ρ+γm3η+γm4λ+γm4ρ+γm4η]c1.

(26)

不難看出，要k1=0，只需

把式(27)代入式(26),再令k2=0,解出β,即得式(17).再把式(17)代入式(27)可得式(16).把式(16)和式 (17)代入式(24),經簡化可得M的表達式(19)和誤差公式(18).由條件知,M的分母不為零.定理1證畢.

2 特例與數值分析

通過調整δ,β,γ,ρ,η,λ的取值，式(14)可以給出各種各樣的迭代法.如由式(2)定義的Newton法(簡記NM);由式(3)定義的Halley法(簡記HM);由式(4)定義的Chebyshev法(簡記CM);由式(5)定義的Osada法(簡記OM);由式(6)定義的Chun-Neta法(簡記CNM);由式(7)定義的Chun-Bae-Neta法(簡記CHBM)及其2個特例——由式(8)定義的CHBM1法和由式(9)定義的CHBM2法;由式(10)定義的Biazar-Ghanbari法(簡記BGM)和它的特例——由式(13)定義的BGM1法；以及以下2個新的方法：

1)在式(14)中,取δ=-m3,β=4m2,γ=m(m-1)2,ρ=-m2(3-m),η=-2m(m2-2m-1),λ=(m-1)2(m+1),則得以下新的3階方法(簡記YXJM1):

式(28)中：Lf(x)如式(15)所定義;φ1(x)=-m3x2+4m2x+m(m-1)2;φ2(x)=-m2(3-m)x2-2m(m2-2m-1)x+(m-1)2(m+1).

式(29)中:φ3(x)=2m2x2+m(4-7m)x+(5m2-5m+2);φ4(x)=2m2x2-4m2x+2m2.

為說明方法(14)的有效性,筆者用上面提到過的其中不含有任意參數的10種方法:NM法、HM法、CM法、OM法、CNM法、CHBM1法、CHBM2法、BGM1法、YXJM1法和YXJM2法分別求解如下2個方程的重根：

很明顯，f1(x)有4重根，f2(x)有3重根.

所有的數值計算都用數學軟件Maple 14.0在PC上進行,并且設置128位計算精度,即(Digits:=128).本文選擇如下的迭代終止條件:1)|f(xn+1)|≤ε;或2)迭代次數n≥1 000.

本文取ε=10-32.對f1(x)分別從2個初值x0=-0.5和x0=-2進行迭代;對f2(x)分別從2個初值x0=2和x0=1.5進行迭代.若1)滿足，則用x*xn+1作為精確解α的近似值;若2)滿足,則認為方法發散.把計算結果列在表2中.從表2結果可看出，本文方法比絕大多數其他方法收斂都快，且在困難情形下更有用.如當從初值x0=-0.5出發求f1(x)的重根時，有6種方法發散,而本文的2個方法只分別用了3或4步就求得了滿足精度要求的近似解.

3 結 語

本文提出了一族新的求解非線性方程重根的3階迭代方法，這個族具有很好的廣泛性，包含了一系列已知的方法,且具有很好的魯棒性，可很快求出其他方法無法求解的一些根.本文的這些方法是單點單步方法.還有一些學者從多點和多步2個方面構造求重根的高階迭代法.對于多點法，可參閱文獻[7]及其中的參考文獻；對于多步法，可參閱文獻[8]及其中的參考文獻.對事先不知重數的情形,高階迭代法的文獻并不多,筆者將在此方向做點工作.

[1]Neta B.New third order nonlinear solvers for multiple roots[J].Appl Math Comp,2008,202(1):162-170.

[2]Traub J F.Iterative methods for the solution of equations[M].New Jersey:Prentice Hall,1964.

[3]Osada N.An optimal multiple root-finding method of order three[J].J Comput Appl Math,1994,51(1):131-133.

[4]Chun C,Neta B.A third-order modification of Newton′s method for multiple roots[J].Appl Math Comput,2009,211(2):474-479.

[5]Chun C,Bae H,Neta B.New families of nonlinear third-order solvers for finding multiple roots[J].Comput Math Appl,2009,57(9):1574-1582.

[6]Biazar J,Ghanbari B.A new third-order family of nonlinear solvers for multiple roots[J].Comput Math Appl,2010,59(10):3315-3319.

[7]Kumar S,Kanwar V,Singh S.On some modified families of multipoint iterative methods for multiple roots of nonlinear equations[J].Appl Math Comp,2012,218(14):7382-7394.

[8]Chuna C,Neta B.Basins of attraction for several third order methods to find multiple roots of nonlinear equations[J].Appl Math Comp,2015,268(1):129-137.

(責任編輯 陶立方)

Afamilyofmultiplerootfindingmethodswithcubicalconvergence

PAN Yunlan

(CollegeofMathematics,PhysicsandInformationEngineering,ZhejiangNormalUniversity,JinhuaZhejiang321004,China)

A new family of iterative methods to find multiple roots of a nonlinear equation was obtained. Third order convergence was proved for these methods and iteration errors were given. The generality of the family was presented: the family includes, as particular cases, several well known families and methods. By comparing the proposed methods with some other methods through numerical experiments, the robustness and efficiency of the new methods were shown. It was showed that the presented methods converge faster than some other methods and even converge very fast at some cases while the other methods diverge.

iterative method;convergence order;cubical convergence;multiple root

10.16218/j.issn.1001-5051.2015.04.002

2015-06-21；

：2015-09-14

國家自然科學基金資助項目(61170109)

潘云蘭(1967-),女,浙江磐安人,副教授.研究方向:數值逼近;軟件工程.

O241

：A

：1001-5051(2015)04-0366-06