黃娜，馬昌鳳

（福建師范大學(xué) 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院，福建 福州 350007）

＊1求解非線性方程的一個(gè)新的三階迭代算法

黃娜，馬昌鳳＊

（福建師范大學(xué) 數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院，福建 福州 350007）

提出了求解非線性方程實(shí)根的一個(gè)新的迭代方法，并證明了這種方法是三次收斂的.特別地，當(dāng)函數(shù)在零點(diǎn)的三階導(dǎo)數(shù)值為零時(shí)，這種方法是超三次收斂的.此外，通過數(shù)值實(shí)驗(yàn)驗(yàn)證了所做的理論分析.給出了五個(gè)數(shù)值算例，從迭代次數(shù)，所用CPU時(shí)間，誤差以及收斂階這四個(gè)方面，將這個(gè)新的算法與經(jīng)典的牛頓法等三個(gè)算法進(jìn)行比較，數(shù)值結(jié)果表明文章提出的新算法是有效的.

非線性方程；迭代公式；收斂階；數(shù)值實(shí)驗(yàn)

在工程計(jì)算和科學(xué)研究中，如電路和電力系統(tǒng)計(jì)算、非線性力學(xué)、非線性積分和微分方程等許多領(lǐng)域都要遇到非線性方程的求根問題［1，5］.考慮一元非線性方程

本文提出求解非線性方程（1）的一個(gè)新的迭代公式，并證明了這個(gè)公式具有三次收斂速度.特別地，當(dāng)函數(shù)f（x）在零點(diǎn)的三階導(dǎo)數(shù)值為零時(shí)，這種方法是超三次收斂的.本文最后給出了五個(gè)數(shù)值算例，從迭代次數(shù)，所用CPU時(shí)間，誤差以及收斂階這四個(gè)方面，將這個(gè)新的算法及其離散形式與經(jīng)典的牛頓法等三個(gè)算法進(jìn)行比較，數(shù)值結(jié)果表明該算法是有效的.

1 算法及其收斂性分析

不妨設(shè)函數(shù)f∶［a，b］?R→R充分可微，x＊為非線性方程f（x）＝0的根.對(duì)任給的x k∈［a，b］，函數(shù)f（x）在點(diǎn)x k處的Taylor展開式為［2］：

由（5）式可以構(gòu)造出一種迭代方法來求解非線性方程（1），即經(jīng)典的牛頓法

已經(jīng)證明迭代公式（6）至少具有二階收斂速度.下面我們來推導(dǎo)具有更高收斂階的迭代公式.事實(shí)上，若用右矩形公式近似（3）式中的積分，則有

將（7）代入（3）并略去誤差項(xiàng)R［F（x）］（仍然使用等號(hào)“＝”），可得

令x k＋1是方程（8）的解，可得下面的迭代公式

由于迭代公式（9）是隱式格式，因此可以用牛頓迭代公式（6）進(jìn)行預(yù)報(bào)，于是可以得到求解非線性方程（1）的一個(gè)預(yù)報(bào)校正公式.

下面我們來考慮迭代格式（10）－（11）的收斂性和收斂速度.我們有如下的定理：

定理1.1 設(shè)函數(shù)f：［a，b］→R充分可微，x＊為非線性方程f（x）＝0的根，且x＊∈［a，b］，則迭代方法（10）－（11）是三次收斂的，且誤差滿足如下等式

等式兩邊同時(shí)乘以f′（y k），并整理得

在（2）式中令x＝x＊，并注意到f（x＊）＝0，有

證明記ek＝x k－x＊，由（11）可以得到

上式可以整理為

上式兩邊同時(shí)除以f′（x k），有

另一方面，在（2）式中令x＝y(tǒng) k，有

2 數(shù)值實(shí)驗(yàn)

下面通過五個(gè)數(shù)值實(shí)驗(yàn)，從迭代次數(shù)、CPU時(shí)間、誤差以及收斂階這四個(gè)方面，將新的算法（10）－（11）（簡記為HM）與經(jīng)典的牛頓法（記為NM），Aslam Noor和 Waseem的算法［4］（記為NR1），Cordero和Torregrosa的算法［3］（記為CT）進(jìn)行比較.整個(gè)實(shí)驗(yàn)在Pentium（R）Dual－Core CPU T4400＠2.20 GHz的PC上執(zhí)行，軟件平臺(tái)為MATLAB 2009b.收斂階按照下列公式近似計(jì)算（參見文［3］）.

表1 取初值為x0＝0.5的數(shù)值結(jié)果Table 1 Numerical results of the example 1 for initial value x 0＝0.5

算例2函數(shù)f（x）＝sin（x）ex＋ln（x2＋1），其中x∈R，該函數(shù)有唯一零點(diǎn)x＊＝0.取初值為x0＝1，計(jì)算結(jié)果如表2.

表2 取初值為x 0＝1的數(shù)值結(jié)果Table 2 Numerical results of the example 2 for initial value x 0＝1

算例3函數(shù)f（x）＝x2－4，其中x∈R，該函數(shù)有兩個(gè)零點(diǎn)，分別為x＊＝－2.x＊＊＝2，分別取初值為x0＝－1.5（此時(shí)迭代序列收斂于x＊）和x0＝1.5（此時(shí)迭代序列收斂于x＊＊）計(jì)算結(jié)果分別如表3和表4.

表3 取初值為x0＝－1.5的數(shù)值結(jié)果Table 3 Numerical results of the example 3 for initial value x 0＝－1.5

算例4函數(shù)f（x）＝（x－1）6－1，其中x∈R，該函數(shù)有兩個(gè)零點(diǎn)，分別為x＊＝0，x＊＊＝2.分別取初值為x0＝0.5（此時(shí)迭代序列收斂于x＊）和x0＝4（此時(shí)迭代序列收斂于x＊＊）計(jì)算結(jié)果分別如表5和表6.

表5 取初值為x0＝0.5的數(shù)值結(jié)果Table 5 Numerical results of the example 4 for initial value x 0＝0.5

表6 取初值為x0＝4的數(shù)值結(jié)果Table 6 Numerical results of the example 4 for initial value x 0＝4.0

從以上四個(gè)算例可以看出，本文所提出的迭代格式（10）－（11）是有效的，無論是迭代次數(shù)或所用CPU時(shí)間均優(yōu)于其它算法.

［1］ Burden R L，F(xiàn)aires J D.Numerical Analysis（7th ed.）［M］.Boston：PWS Publishing Company，2001.

［2］ Ortega J M，Rheinboldt W C.Iterative Solution of Nonlinear Equations in Several variables［M］.Academic Press，1970.

［3］ Cordero A，Torregrosa J R.Variants of Newton’s Method Using Fifth－order Quadrature Formulas［J］.ApplMathComput，2007，190：686－698.

［4］ Aslam Noor M，Waseem M.Some Iterative Methods for Solving a System of Nonlinear Equations［J］.ApplMathComput，2009，57：101－106.

［5］ 馬昌鳳，林偉川.現(xiàn)代數(shù)值計(jì)算方法（MATLAB版）［M］.北京：科學(xué)出版社，2008：6.

A New Third－order Method for Solving Nonlinear Equation

HUANG Na，MA Chang－feng
（SchoolofMathematics＆ComputerScience，F(xiàn)ujianNormalUniversity，F(xiàn)uzhou350007，China）

A new iterative method for solving nonlinear equation is introduced.The method is proved to be cubic convergent.Especially，when the third order derivative of functionf（x）is equal to zero in the zero point of the function，this method will be supercubic convergent.In addition，five numerical examples are given to illustrate the efficiency and the performance of the newly developed method confirms the theoretical results.

nonlinear equation；iterative method；convergence order；numerical experiments

O241.7

A

0253－2395（2012）03－0460－05＊

2011－04－28；

2012－03－21

國家自然科學(xué)基金（11071041）

黃娜（1988－），女，福建閩清人，在讀博士研究生.＊通訊作者：macf＠fjnu.edu.cn