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

非线性方程牛顿迭代法与斯特芬森迭代法的研究与比较由刀豆文库小编整理，希望给你工作、学习、生活带来方便，猜你可能喜欢“牛顿迭代法相关研究”。

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

申林坚

（南昌航空大学 测试与光电工程学院 江西 南昌 330063）

摘要：本文针对一个具体的非线性方程3x2ex0进行研究，首先作出了了函数

f(x)3x2ex的图像，大体判定其零点（即方程解）在（3,4）区间内，接着用牛顿迭代法和斯特芬森迭代法进行求解分析，牛顿法的迭代公式为 xk1xkf(xk)，f(xk)斯特芬森迭代法公式为

yk(xk),zk(yk),(ykxk)2

xk1xk,zk2ykxk记录两种方法求得指定精度解所需迭代次数及所需计算时间，并对其优缺点 进行了分析。

关键词：非线性方程；牛顿迭代法；斯特芬森迭代法

引言

非线性是实际问题中经常出现的，并且在科学与工程计算中的地位越来越重要，很多我们熟悉的线性模型都是在一定条件下由非线性问题简化得到的，为得到 更符合实际的解答，往往需要直接研究非线性模型，从而产生非线性科学，它是21世纪科学技术发展的重要支柱。本论文通过对特定非线性方程3x2ex0进行求解，介绍了两种常用的迭代法牛顿迭代法和斯特芬森迭代法，详尽阐述了其各自的数学几何原理及优缺点比较，从而更深入的理解非线性方程的迭代法求解。

正文

一.作出f(x)的图像，确定隔根区间 在Matlab中输入以下指令并回车：

x=(-10:0.001:10);y=3*x.^2-exp(x);plot(x,y);

grid on;0.5x 1040-0.5-1-1.5-2-2.5-10-8-6-4-20246810

图1 得到图1所示f(x)的图像，易知，当x10及x10时，f(x)无零点 将y轴方向放大，输入命令axis([-10 10-2 2])，得到图2

21.510.50-0.5-1-1.5-2-10-8-6-4-20246810

图2 可知函数有三个零点，隔根区间为（-2,0），（0,2），（2,4）将x轴方向放大，输入命令axis([-2 4-2 2])，得到图3 21.510.50-0.5-1-1.5-2-2-101234

图3 可将隔根区间进一步缩小为（-1,0），（0,1），（3,4）

二.牛顿迭代法求区间（3,4）中的根

对于方程f(x)0，如果f(x)是线性函数，则它的求根是容易的。牛顿法实质上是一种线性化方法，其基本思想是将非线性方程f(x)0逐步归结为某种非线性方程来求解。

设已知方程f(x)0有近似根xk（假定f(xk)0），将函数在点xk展开，有

f(x)f(xk)f(xk)(xxk)，于是方程f(x)0可近似表示为

f(xk)f(xk)(xxk)0.这是个线性方程，记其根为xk1，则xk1的计算公式为

xk1xk这就是牛顿法。

f(xk),k0,1,..., f(xk)牛顿法有明显的几何解释。方程f(x)0的根x可解释为曲线yf(x)与x轴的交点的横坐标。设xk是根x的某个近似值，过曲线yf(x)上横坐标为xk的点Pk引切线，并将该切线与x轴的交点的横坐标xk1作为x的新的近似值。注意到切线方程为

***yf(xk)f(xk)(xxk)

这样求得的值xk1必满足f(xk)f(xk)(xxk)0。由于这种几何背景，牛顿法亦称切线法。

下面列出牛顿法的计算步骤：

步骤1 准备 选定初始近似值x0，计算f0f(x0)，f0f(x0)步骤2 迭代 按公式

x1x0f0 f0迭代一次，得到新的近似值x1，计算f1f(x1),f1f(x1)

步骤3 控制 如果x1满足哦1或f12，则终止迭代，以x1作为所求的根；否则转步骤4.此处1,2是允许误差，而

x1x0,当x1C时, x1x0,当x1C时，x1其中C是取绝对误差或相对误差的控制常数，一般可取C=1.步骤4 修改 如果迭代次数达到预先指定的次数N，或者f10，则方法失败；否则以(x1,f1,f1)代替(x0,f0,f0)转步骤2继续迭代

Matlab计算程序如下，取初始迭代值x0=3：

i=0;y=3;z=1;

while(i=10^(-8))x=y;

y=x-(3*x^2-exp(x))/(6*x-exp(x));

if abs(y)

else

z=abs((y-x)/y);

end i=i+1;end

format long;disp(y);disp(i);

输出结果为y= 3.***，i=9 可知，使用牛顿迭代法，初值为3时，需迭代9次可使eps

三.斯特芬森迭代法求区间(3,4)中的根

埃特金方法不管原序列{xk}是怎样产生的，对{xk}进行加速计算，得到序列{xk}。如果把埃特金加速技巧与不动点迭代结合，可得到如下的迭代法：

yk(xk),zk(yk),kxk)xk1xkz(ky2ykxk,2k0,1,...*称为斯特芬森迭代法。它可以这样理解，我们要求x(x)的根x，令(x)(x)x，(x*)(x*)x*0，已知x*的近似值xk及yk，其误差分别为

(xk)(xk)xkykxk

(yk)(yk)ykzkyk

把误差“外推到零”，即过xk,(xk)及yk,(yk)两点做线性插值函数，它与x轴交点就是xk1，即方程

(xk)的解

(yk)(xk)ykxk(xxk)0

(xk)(ykxk)2xxk(ykxk)xkxk1

(yk)(xk)zk2ykxk斯特芬森迭代法的另一种表达方式如下：

xk1(xk),k0,1,...,其中

(x)xxx2((x))2(x)x.实验表明，即便用不动点迭代法不收敛，用斯特芬森迭代法仍可能收敛。

1.取(x)ln(3x)，迭代初值为3进行迭代，Matlab程序如下：

i=0;y=3;z=1;

while(i=10^(-8))x=y;

y1=log(3*x^2);

y=x-(y1-x)^2/(log(3*y1^2)-2*y1+x);

if abs(y)

else

z=abs((y-x)/y);

end i=i+1;end

format long;disp(y);disp(i);2输出结果为y= 3.***，i=4;可知，使用斯特芬森迭代法，取不动点函数为(x)ln(3x)，初值为3时，需迭代4次可

2使eps

2.取(x)3x2exx，迭代初值为3进行迭代，Matlab程序如下：

i=0;y=3;z=1;

while(i=10^(-8))x=y;

y1=3*x^2-exp(x)+x;

y=x-(y1-x)^2/(3*y1^2-exp(y1)+y1-2*y1+x);

if abs(y)

else

z=abs((y-x)/y);

end i=i+1;end

format long;disp(y);disp(i);

输出结果为y= 3.***，i=147;可知，使用斯特芬森迭代法，取不动点函数为(x)3xex，初值为3时，需迭代147次可使eps

2xex3.取(x)，迭代初值为3进行迭代，Matlab程序如下：

3i=0;y=3;z=1;

while(i=10^(-8))x=y;

y1=sqrt(exp(x)/3);

y=x-(y1-x)^2/(sqrt(exp(y1)/3)-2*y1+x);

if abs(y)

else

z=abs((y-x)/y);

end i=i+1;end

format long;disp(y);disp(i);输出结果为y= 3.***，i=10;

ex可知，使用斯特芬森迭代法，取不动点函数为(x)，初值为3时，需迭代10次可

3使eps

牛顿法的迭代函数固定，收敛速度较快，但求f(x)可能比较繁琐

斯特芬森迭代法收敛速度快，但要选择合适的迭代函数，因次如何构造一个合适的迭代函数成为关键。

参考文献

[1]蔡旭晖 刘卫国 蔡立燕 MATLAB基础与应用教程 北京:人民邮电出版社 2009 [2]李庆扬 王能超 易大义 数值分析(第5版)北京:清华大学出版社 2008 [3]高成 赖志国 Matlab图像处理与应用(第2版)北京:国防工业出版社 2007

The research and comparison of Newton's iterative method and Steffen Sen iteration method for nonlinear equation

Shen Linjian(Nanchang Institute of test and opto electronic engineering, Jiangxi University of Aeronautics 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 analysis, 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

个人心得体会

首先，我觉得课堂教学条件比较差，那么大的一个教室坐满了人，不能保证每个学生能够听清，看清每一个知识点。因此建议以后可开展小班教学。其次，坦率得说，关于这门课的知识点，我完全是通过自学获得的，从课堂上得到的少之又少，并不是老师讲的不好，只是数学本就枯索抽象，在课堂上并不能马上领悟，因而没有兴趣继续听下去。我相信，很大一部分同学是跟我有同感的。

另外，大作业这种考查形式还是不错的，并不是很难，但要花功夫。在完成大作业的过程中，其实是一个对所学知识深入在理解的过程，通过与matlab相结合，提升了自己运用所学知识解决实际问题的能力。

我是比较赞同老师对学生严格要求的，但并不是给一个高高在上的目标，让学生难以企及，而是要循循善诱，但凡有一点上进心都能有所收获，在这点上我感觉郑老师做得不错。

最后，我觉得在教学过程中可以引进给出问题，课堂小组讨论这种形式。让学生自己来学习。当然，这是建立在小班教学的基础上的。以上仅是个人的一点拙见，请老师参考。

非线性方程求根的方法简介与例题

非线性方程f(x)求根主要可以采用下面三种方法，下面简单介绍下，并附例题，让解法更一目了然。 1) 二分法简介：计算步骤如下：例题：2) 不动点迭代，也叫简单迭代。隐式化为显式，迭代法是......

比较优势与竞争优势研究

比较优势与竞争优势研究比较优势是经济学的一个重要概念，它最初产生在对于国际贸易模式的关注。竞争优势理论成熟于迈克尔·波特，其初始研究集中在企业战略与企业的竞争力问题......

德治与法治的比较研究

目录内容摘要„„„„„„„„„„„„„„„„„„„„„„„„„ 2 一．德治与法治的历史背景及其发展„„„„„„„„„„„„„ 2 (一).德治历史背景及其发展„„„„„......

可口可乐与王老吉的比较研究

第五章可口可乐与红罐王老吉营销策略实践比较研究可口可乐和红罐王老吉营销策略存在许多的共同点，都符合饮料营销的流程和特征，都主要包括产品、价格、分销、促销四个方面，都同......

李白与杜甫诗歌比较研究

李白与杜甫诗歌比较研究......