非线性方程求解的三类Newton法的比较.docVIP

【标题】非线性方程求解的三类Newton法的比较 【作者】唐 虹 【关键词】Newton法??非线性方程??迭代法??收敛阶 【指导老师】周 均 【专业】数学与应用数学 【正文】0引言经典Newton迭代法是求解非线性方程的常用算法?，求单实根时收敛阶为二阶.目前有诸多研究对求解非线性方程的方法进行改进?，使用的很多有效的迭代法都是以经典Newton法为基础或由它派生而来.有避免导数计算的超线性收敛的弦截法、抛物线法?；有避免二阶导数计算且具有至少三阶收敛速度的迭代法?和二阶敛速指数迭代法?等等.为了进一步提高算法的收敛速度，本文基于经典Newton法的思想研究一种促使牛顿迭代加速收敛的新方法.利用原函数与反函数的关系，将反函数利用Taylor公式展开，取前四项构造迭代公式；并为计算的简便性，利用一阶差商来近似代替?，利用二阶差商来近似代替?求解?，利用三阶差商来近似代替?求解?，构造新的迭代公式.通过计算实例，说明经典Newton法、改进的Newton法和修正的Newton法的计算量和收敛速度.1预备知识定义1??设?是方程?的根，当?充分接近?时，由迭代式?得到的序列?收敛于?.若常数?和正常数?，使得?，则称序列?是?阶收敛的，或称?的收敛阶是?阶，称?是?阶迭代函数.显然，收敛阶?的大小刻画了迭代序列?的收敛速度，?越大，收敛的速度越快.若?，称序列?是线性收敛的，若?，称序列?是超线性收敛的，若?，称序列?是平方收敛的.引理1?（收敛阶数的判定?）设迭代函数?满足下列条件:1） 在?的某一?邻域有?阶的连续导数；2） ?；3） ?，?；则当初始值?充分接近?时，由?产生的序列?是?阶收敛的，且??引理2??（利用插值多项式?作为?的近似函数建立的高阶导数数值微分公式关于其截断误差的结论?若?在?上有直到?阶导数，?为?在?个结点?上的?次插值多项式，则对任何?有?其中?，?且依赖于?.引理3??若函数?在点?存在直至?阶导数，则有?，即?2经典Newton法和改进的Newton法2.1经典Newton法解非线性方程?的经典Newton法的迭代公式??????????????????????????????式(2.1)且(2.1)式求单实根时，仅有二阶收敛速度.例1. 设??，求?的根，精确到?.解:取迭代初始值?，用经典Newton法：???????????求满足精度要求的近似根，计算结果见表2.1:表2.1???0 1.000000000 0.718281828 1 0.735758882 0.087065200 0.26424112 0.694042300 0.001791000 0.04171663 0.693147581 0.000000080 0.00089474 0.693147181 0.000000000 0.00000002.2改进的Newton法在日常生活中，我们会遇到一些有关数学方面的问题，并且很多问题的数学模型最后都转化为求解非线性方程，要想提高解决实际问题的速度，即提高求解非线性方程问题的速度，我们不妨提高迭代法的收敛速度.用?与?之和代替?，这样得到改进的Newton法???????????????????????????????????式(2.2)且至少是三阶收敛的.由迭代公式(2.2)可得迭代函数?，其中?.不妨设方程?的根为?，函数?在?的邻域内有四阶连续导数，且?.则有?，?，?.对?求导得，?因此有?，再对?求导得，???因此有，?.所以由引理1得，迭代式?在?邻近至少是三阶收敛的.现用迭代式?解例1，取迭代初始值?，求满足精度要求的近似根，计算结果见表2.2:表2.2???0 1.000000000 0.718281828 1 0.694042300 0.001791040 0.30595772 0.693147181 0.000000000 0.00089513 0.693147181 0.000000000 0.0000000若取迭代初始值?，迭代式?经过?次迭代获得要求精度的近似解，而迭代式?则要经过?次迭代才获得要求精度的近似解，迭代结果见表2.3:表2.3????????0 200 201 190.9999999 1 181.36787942 180.9999999 2 161.99999993 170.9999999 3 141.99999964 16.0000