非线性方程牛顿迭代法与斯特芬森迭代法的研究与比较

1、非线性方程牛顿迭代法与斯特芬森迭代法的研究与比较申林坚（南昌航空大学 测试与光电工程学院 江西 南昌 330063）摘要：本文针对一个具体的非线性方程进行研究，首先作出了了函数 的图像，大体判定其零点（即方程解）在（3,4）区间内，接着用牛顿迭代法和斯特芬森迭代法进行求解分析，牛顿法的迭代公式为 ，斯特芬森迭代法公式为记录两种方法求得指定精度解所需迭代次数及所需计算时间，并对其优缺点进行了分析。关键词：非线性方程；牛顿迭代法；斯特芬森迭代法引言非线性是实际问题中经常出现的，并且在科学与工程计算中的地位越来越重要，很多我们熟悉的线性模型都是在一定条件下由非线性问题简化得到的，为得到更符合实际的解

2、答，往往需要直接研究非线性模型，从而产生非线性科学，它是21世纪科学技术发展的重要支柱。本论文通过对特定非线性方程进行求解，介绍了两种常用的迭代法牛顿迭代法和斯特芬森迭代法，详尽阐述了其各自的数学几何原理及优缺点比较，从而更深入的理解非线性方程的迭代法求解。正文一.作出的图像，确定隔根区间在Matlab中输入以下指令并回车：x=(-10:0.001:10);y=3*x.2-exp(x);plot(x,y);grid on;图1得到图1所示的图像，易知，当及时，无零点将y轴方向放大，输入命令axis(-10 10 -2 2)，得到图2图2可知函数有三个零点，隔根区间为（-2,0），（0,2），（

3、2,4）将x轴方向放大，输入命令axis(-2 4 -2 2)，得到图3图3可将隔根区间进一步缩小为（-1,0），（0,1），（3,4）二.牛顿迭代法求区间（3,4）中的根对于方程，如果是线性函数，则它的求根是容易的。牛顿法实质上是一种线性化方法，其基本思想是将非线性方程逐步归结为某种非线性方程来求解。设已知方程有近似根（假定），将函数在点展开，有，于是方程可近似表示为.这是个线性方程，记其根为，则的计算公式为这就是牛顿法。牛顿法有明显的几何解释。方程的根可解释为曲线与x轴的交点的横坐标。设是根的某个近似值，过曲线上横坐标为的点引切线，并将该切线与x轴的交点的横坐标作为的新的近似值。注意到切线

4、方程为这样求得的值必满足。由于这种几何背景，牛顿法亦称切线法。下面列出牛顿法的计算步骤：步骤1 准备 选定初始近似值，计算，步骤2 迭代 按公式迭代一次，得到新的近似值，计算步骤3 控制 如果满足哦或，则终止迭代，以作为所求的根；否则转步骤4.此处是允许误差，而其中C是取绝对误差或相对误差的控制常数，一般可取C=1.步骤4 修改 如果迭代次数达到预先指定的次数N，或者，则方法失败；否则以代替转步骤2继续迭代Matlab计算程序如下，取初始迭代值=3：i=0;y=3;z=1;while(i<=100&&z>=10(-8) x=y; y=x-(3*x2-exp(x)/(

5、6*x-exp(x); if abs(y)<1 z=abs(y-x); else z=abs(y-x)/y); end i=i+1;endformat long;disp(y);disp(i);输出结果为y= 3.733079028632816，i=9可知，使用牛顿迭代法，初值为3时，需迭代9次可使eps<10-8,近似解为3.73307903三.斯特芬森迭代法求区间(3,4)中的根埃特金方法不管原序列是怎样产生的，对进行加速计算，得到序列。如果把埃特金加速技巧与不动点迭代结合，可得到如下的迭代法：称为斯特芬森迭代法。它可以这样理解，我们要求的根，令，已知的近似值及，其误差分别为把

6、误差“外推到零”，即过及两点做线性插值函数，它与x轴交点就是，即方程的解斯特芬森迭代法的另一种表达方式如下：其中实验表明，即便用不动点迭代法不收敛，用斯特芬森迭代法仍可能收敛。1. 取，迭代初值为3进行迭代，Matlab程序如下：i=0;y=3;z=1;while(i<=1000&&z>=10(-8) x=y; y1=log(3*x2); y=x-(y1-x)2/(log(3*y12)-2*y1+x); if abs(y)<1 z=abs(y-x); else z=abs(y-x)/y); end i=i+1;endformat long;disp(y);di

7、sp(i);输出结果为y= 3.733079028632815，i=4;可知，使用斯特芬森迭代法，取不动点函数为，初值为3时，需迭代4次可使eps<10-8,近似解为3.73307903，与牛顿迭代法结果吻合。2.取，迭代初值为3进行迭代，Matlab程序如下：i=0;y=3;z=1;while(i<=1000&&z>=10(-8) x=y; y1=3*x2-exp(x)+x; y=x-(y1-x)2/(3*y12-exp(y1)+y1-2*y1+x); if abs(y)<1 z=abs(y-x); else z=abs(y-x)/y); end i=

8、i+1;endformat long;disp(y);disp(i);输出结果为y= 3.733079028632814，i=147;可知，使用斯特芬森迭代法，取不动点函数为，初值为3时，需迭代147次可使eps<10-8,近似解为3.733079033.取，迭代初值为3进行迭代，Matlab程序如下：i=0;y=3;z=1;while(i<=1000&&z>=10(-8) x=y; y1=sqrt(exp(x)/3); y=x-(y1-x)2/(sqrt(exp(y1)/3)-2*y1+x); if abs(y)<1 z=abs(y-x); else

9、z=abs(y-x)/y); end i=i+1;endformat long;disp(y);disp(i);输出结果为y= 3.733079028632815，i=10;可知，使用斯特芬森迭代法，取不动点函数为，初值为3时，需迭代10次可使eps<10-8,近似解为3.73307903由以上三种不同迭代函数可知，迭代函数不同，斯特芬森迭代法的迭代次数不同。综合比较牛顿法的迭代函数固定，收敛速度较快，但求可能比较繁琐斯特芬森迭代法收敛速度快，但要选择合适的迭代函数，因次如何构造一个合适的迭代函数成为关键。参考文献1蔡旭晖 刘卫国 蔡立燕 MATLAB基础与应用教程 北京:人民邮电出版社

10、 20092李庆扬 王能超 易大义 数值分析(第5版) 北京:清华大学出版社20083高成 赖志国 Matlab图像处理与应用(第2版) 北京:国防工业出版社 2007The research and comparison of Newton's iterative method and Steffen Sen iteration method for nonlinear equationShen Linjian(Nanchang Institute of test and opto electronic engineering, Jiangxi University of Aeron

11、autics and Astronautics, Nanchang 330063)Abstract: In this paper, a specific nonlinear equation is studied, firstly, the function of the image, generally determine its zero (the equation solution) in the (3,4) interval, followed by Newton iterative method and Steffen Sen iteration method for analysi

12、s, Newton iterative formula for, Steffen Sen iterative formula for the record of the two methods to obtain the specified accuracy of the required number of iterations and the required calculation time, and its advantages and disadvantages are analyzed.Key words: nonlinear equation; Newton iterative method; Steffen Sen iteration method个人心得体会首先，我觉得课堂教学条件比较差，那么大的一个教室坐满了人，不能保证每个学生能够听清，看清每一个知识点。因此建议以后可开展小班教学。其次，坦率得说，关于这门课的知识点，我完全是通过自学获得的，从课堂上得到的少之又少，并不是老师讲的不好，只是数学本就枯索抽象，在课堂上并不能马上领悟，因而没有兴趣继续听下去。我相信，很大一部分同学是跟我有同感的。另外，大作业这种考查形式还是不错的，并不是很难，但要花功夫。在完成大作业的过程中，其实是一个