不用求導(dǎo)含參數(shù)的三階收斂迭代方法

裕靜靜，江平

(合肥工業(yè)大學(xué)數(shù)學(xué)學(xué)院，合肥230009)

[摘要]基于弦截法和Steffensen迭代法，本文提出了求解非線性方程f(x)=0的兩種帶有參數(shù)的迭代算法，這兩種算法不帶有導(dǎo)數(shù)，且經(jīng)過收斂性分析證明至少是三階收斂的，最后用數(shù)值試驗驗證了本文兩種迭代算法的有效性和優(yōu)越性.

[關(guān)鍵詞]牛頓迭代； 弦截法； Steffensen迭代法； 非線性方程

[收稿日期]2014-10-12

[基金項目]中央高校基本科研業(yè)務(wù)費專項(J2014HGXJ0077)

[中圖分類號]O241.7[文獻(xiàn)標(biāo)識碼]A

1引言

科學(xué)發(fā)展中經(jīng)常用到非線性方程，我們知道牛頓迭代法(CN)[1]是求解非線性方程f(x)=0的經(jīng)典算法.近年來，科學(xué)家們提出了很多種對牛頓迭代法進(jìn)行改進(jìn)的算法，如文獻(xiàn)[2-4]，但是這些算法需要調(diào)用導(dǎo)數(shù)值.對于復(fù)雜函數(shù)而言，求導(dǎo)比較困難，無疑增加了計算量，于是科學(xué)家們又研究了不用求導(dǎo)的改進(jìn)算法，例如對弦截法的改進(jìn)還有不動點的改進(jìn).由楊明波等人在文獻(xiàn)[5]中提出的求解非線性方程的二重弦截法，它的收斂階為2.618.還有吳新元[6]提出的兩族新的牛頓類方法.文獻(xiàn)[7]和[8]中，鄭權(quán)在文獻(xiàn)[6]的基礎(chǔ)上，提出了具有參數(shù)的不帶導(dǎo)數(shù)的平方收斂的歐拉型迭代法，它的收斂階是二階的.文獻(xiàn)[9]中Mehdi Dehghan和Masoud Hajarian所提出的幾種不用求導(dǎo)的平方收斂和三階收斂的迭代法.本文基于文獻(xiàn)[6]和[9]提出了兩種含有參數(shù)的收斂迭代算法，這兩種新算法不需要求導(dǎo)數(shù)，且收斂性分析證明它們至少三階收斂，也通過數(shù)值試驗證明了它們的有效性.

2不用求導(dǎo)含有參數(shù)的迭代格式

本文所提出的兩種算法，它們的迭代格式為

(1)

以及

(2)

其中μ∈R，下面給出新迭代格式的推導(dǎo)過程.

首先給出一個動力系統(tǒng)

(3)

其中U(x*)為x*的鄰域，可見方程f(x)=0的根x*是動力系統(tǒng)(3)平衡點，反之也是一樣的.利用不同的數(shù)值方法求解(3)可以得到不同的求根方法，現(xiàn)在用Euler方法對其進(jìn)行積分可得到

(4)

因為求導(dǎo)帶來不便，可以用差商

代替(4)中的導(dǎo)數(shù)f′(xn)，可得到

(5)

令(5)中h=1，則(5)變?yōu)?/p>

(6)

迭代格式(6)正是[4]中鄭權(quán)提出的帶參的不帶導(dǎo)數(shù)的平方收斂的迭代算法格式.根據(jù)上述的迭代格式，我們可以繼續(xù)改進(jìn)，與弦截法的迭代格式

結(jié)合，得到一個新的迭代格式

(1)

其中μ∈R.

類似地，如果采用差商

代替(4)中的導(dǎo)數(shù)f′(xn)，那么根據(jù)迭代格式(1)的推導(dǎo)過程，可以得到另一種不用求導(dǎo)含參數(shù)的迭代格式，如下所示

(2)

(1)和(2)就是本文所要討論的兩種不用求導(dǎo)含參數(shù)的三階收斂迭代格式，下面對它們進(jìn)行收斂性分析.

3收斂性分析

(7)

如果所選取的參數(shù)

證設(shè)xn在x*的去心鄰域內(nèi)，令en=xn-x*，由帶Peano余項的Taylor公式得

(8)

(9)

由(8)以及泰勒公式展開可得

(10)

由(9)，(10)可得

μf2(xn)+f(xn+f(xn))-f(xn)

(11)

因此利用輾轉(zhuǎn)相除法，由(9)和(11)可得

(12)

由(12)可知

所以可以得到f(yn)在x*處的泰勒展開式為

所以可得

(13)

以及

(14)

由(13)和(14)利用輾轉(zhuǎn)相除法可得

由此可得

(15)

根據(jù)誤差方程的定義以及(15)可知

(16)

(17)

定理2的證明過程參照定理1即可得出結(jié)論.

4數(shù)值試驗

(i)f1(x)=cosx+e-x,x*=1.74613953040801;

(ii)f2(x)=cosx-x,x*=0.73908513321516;

(iii)f3(x)=2sinx-1,x*=0.52359877559830;

(iv)f4(x)=x3-1,x*=1.00000000000000.

測試結(jié)果如表1所示：

表1　迭代次數(shù)的比較表

5結(jié)束語

通過表1中的迭代次數(shù)可知，本文的兩種迭代算法不管μ取何值，都比經(jīng)典牛頓法和Steffensen加速法的收斂效果要好，而且當(dāng)取得的μ值使得斂速估計r*=0時，收斂效果更好.本文所給出的兩種迭代算法不需要求導(dǎo)數(shù)值，而且可以調(diào)節(jié)參數(shù)μ使得收斂速度改變.通過比較可知，第二種迭代算法的效果比第一種略好.通過計算效率指數(shù)可知，格式(1)的效率指數(shù)是1.4422，格式(2)的效率指數(shù)是1.316，而牛頓迭代法和弦截法的效率指數(shù)是1.414，更加驗證了本文迭代算法的有效性.

[參考文獻(xiàn)]

[1]朱曉臨. 數(shù)值分析[M].合肥:中國科學(xué)技術(shù)大學(xué)出版社，2010.

[2]Ozban A Y. Some new variants of Newton’s method[J]. Appl Math Lett， 2004， 17:667-682.

[3]Weerakoon S， Fernando T G L. A variant of Newton’s method with accelerated third-order convergence [J]. Appl Math Lett， 2000， 13:87-93.

[4]Changbum Chun. Construction of third-order modifications of Newton’s method[J]. Appl Math and Comp， 2007， 189:662-668.

[5]楊明波. 求解非線性方程的二重弦截法[J]. 河南師范大學(xué)學(xué)報，2010， 38(3): 14-16.

[6]吳新元. 對牛頓迭代法的一個重要修改[J]. 應(yīng)用數(shù)學(xué)與力學(xué)，1999， 20(8): 863-866.

[7]鄭權(quán). 具有參數(shù)的不帶導(dǎo)數(shù)的平方收斂的迭代法[J]. 計算數(shù)學(xué)，2003， 25(1):107-112.

[8]鄭權(quán). 斯蒂芬森-牛頓類迭代法的二階收斂性[J]. 吉林大學(xué)學(xué)報，2003，41(2): 134-139.

[9]Mehdi Dehghan， Masoud Hajarian. Some derivative free quadratic and cubic convergence iterative formulas for solving nonlinear equations[J].Appl Math and Comp， 2010， 29:19-30.

Third-order Convergent Iterative Method with

Parameter without Derivation

YUJing-jing，JIANGPing

(School of Mathematics， Hefei University of Technology， Hefei 230009， China)

Abstract:According to the secant method and Steffensen iterative method， this paper puts forward two iterative methods with parameters for solving the approximate solution of nonlinear equations f(x)=0.The two algorithms are without derivative， and we know that the two methods are at least convergent to order three through the convergence analysis and proof. Their effectiveness is validated with numerical test.

Key words: Newton’s iterative method; secant method; Steffensen iterative method; nonlinear equation