一种不使用导数的非线性方程求根算法.pdfVIP

山东农业大学学报 (自然科学版)，2012，43(2)：312—315 JournalofShandongAgriculturalUniversity (NaturalScience) 一 种不使用导数的非线性方程求根算法 李灵芝 ，王云诚 ，姚 磊2 (1．曲阜师范大学 运筹与管理学院，山东 日照 276826；2．山东农业大学信息科学与工程学院，山东 泰安 271018) 摘要：非线性方程求解是计算数学领域的一个基本问题。本文给出求解上述问题的一种新算法，该算法不使 用函数的导数，在一定条件下具有2阶收敛速度，并给出了初步的实验结果。 关键词：非线性方程；迭代法；收敛阶 中图分类号：TP391．41 文献标识码：A 文章编号：1000—2324 (2012)02—0312—04 ADERIVATIVE—FREEMETHODFORSOLVINGNONLINEAREQUATIONS LILing—zhi ，WANG Yun—cheng ，YAO Lei (1．CollegeofOperationsResearchandManagementScience，QufuNormalUniversity，Rizhao276826，China； 2．CollegeofInformationScienceandngineering，ShandongA culturalUniversiyt，Tman271018，China) Abstract：Findingthesolutionofanonlinearequationisoneofthefundamentalproblemsinhterealm ofnumer- icalmathematics．Inthispaper，wepresentanew approachtosolvesuchkindofproblems．Themehtodisde— rivative—freeand2orderconvergenceratecanbereachedundersomeproperassumptions．Preliminarynumed— calexamplesaregiven． Keywords：Nonlinearequation；iterativeformula；rateofconvergence 在科学与工程计算中，常需要求解如下形式的非线性方程 f() =0 (1．1) 这里 D∈R— 在开区间D上是一阶连续可微函数。除了少数具有特殊结构的方程，非线性方程 一 般不能使用解析方法求得其解析解，而只能用数值方法求得其近似解。在求解上述问题的算法中， Newton法简单易行，且具有局部二阶收敛速度，虽然近年来人们构造出许多收敛的阶更高的算法 (例 如 [1—5])，但Newton法仍然是最常用的算法之一。然而，由于Newton法用到厂 ()的导函数，这 对算法的应用带来不便。为了回避牛顿法中的导数计算而提出的弦截法，其收敛阶为P=1．618。 本文利用抛物插值多项式的切线的零点作为方程 (1．1)的近似根，得到一个新的求解非线性方 程的算法，在一定条件下，证明了其收敛速度与牛顿法相同，并通过数值试验验证了算法的有效性。 1 算法的构造与收敛性分析 1．1 算法的构造 设已有近似根 一。， ，取 = + ( 一 )。过这三个点作f()的抛物插值多项式： p()=厂(2xk-xk-1)+’。=-主 (一(2xk-x~_1)) + =_ 二生___l-)7_=—_(生=÷— (一(2一 ^， 一 )) x ) ． 、 、 、 ‘ … ， ， 收稿 日期：2011—09—28 作者简介：李灵芝 (1986一)，女，曲阜师范大学硕士研究生。 通讯作者：Authorforcorrespondence．E—