Newton迭代法浅析

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

目录
第一章:绪论 (2)
第二章 Newton迭代原理 (3)
2.1 一般迭代思想的设计 (3)
2.2 Newton迭代法的原理 (3)
2.3小结: (5)
第三章 Newton迭代法的收敛性 (6)
3.1 Newton迭代法中不收敛的情况 (6)
3.2 定理证明 (7)
3.3 Newton迭代法的收敛性分析 (10)
3.4小结: (12)
第四章两种改进的Newton迭代法 (14)
4.1 改进初值
x的Newton下山法 (14)
4.2 一种新的Newton迭代法加速设计 (15)
4.3小结: (16)
第五章 Newton迭代法的应用 (17)
5.1 Newton迭代法的Matlab实现 (17)
5.2 数值举例 (17)
5.3小结: (20)
总论 (21)
参考文献 (22)
致谢 (23)
第一章绪论
在自然科学和工程技术中很多问题的解决常常归结为解非线性方程(组)或者线性方程(组)代数方程组,例如电学中的网络问题,船体数学放样中建立三次样条函数问题,用最小而乘求是实验数据的曲线拟合问题,用差分或者有限元方法解常微分方程等。

关于非线性方程(组)的求解,一般有两类解法:直接法和迭代法。

我们知道,只有一次、二次和三次方程有规范的求根公式,而高于三次的方程0
)
x
f是不存在求根公式的。

因此求根变得一异常的困难。

而科学计算却
(
很好解决了这一问题,其中最基本的算迭代法了,它对于解决非线性方程(组)的根变得异常方面。

就迭代法而言,Newton迭代法可算是其经典之作。

Newton迭代法又称为Newton-Raphson迭代法,它是Newton在17世纪提出的一种在实数域和复数域上近似求解方程的方法。

牛顿迭代法是求非线性方程(组)根的重要方法之一,其迭代格式简单,且在单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根。

关于Newton迭代法,许多学者为之做了相当多的研究,并且留下了很多经典的文献([2-6])。

Newton迭代法在解决Banach空间中非线性方程或方程组的应用更为重要,如梯形Newton迭代法。

其中,Ksntorovich[7]关于Newton迭代法收敛的工作是解决方程算法现代研究的起点,并给出的Ksntorovich条件。

Newton迭代法是一个重要的计算方法和思想。

牛顿迭代法的主要功能:计算方程时可以比较快速方便的计算出来结果但并不影响计算出来结果的精确度,运用于多种工业设计和数学设计方面,其重要性可见一般。

在我们开设的由李庆扬等编《数值分析》(第五版)中,其第七章7.4节就专门对Newton迭代法进行了一定的讲解。

但是,其书上所讲内容甚微。

这对于一个初学者来说,算是一个障碍。

因此,本人想对Newton迭代法做一个系统的总结。

本文主要从Newton迭代法的基本思想、收敛性、及几种Newton迭代法的变形作一个全方面的介绍。

主要是想在总结前人的经典研究之上,系统的对Newton 迭代法作一下总结,以便更加快速的学习和掌握Newton迭代法。

第二章 Newton 迭代原理
2.1 一般迭代思想的设计
迭代法作为最常用的近似方法,对于求非线性方程(1.1.1)的近似根,我们考虑将方程(1.1.1)转化为等价方程
)(x x ϕ= (2.1.1)
其中,)(x ϕ称为迭代函数。

显然,若*x 是方程(2.1.1)的根,那么*x 必为方程(1.1.1)的;反之已然。

由方程(2.1.1)我们就可以建立迭代算法了。

我们去起根*x 的一个近似值0x 为初始猜测值,然后由
,2,1,0),(1==+k x x k k ϕ (2.1.2)
产生一个序列}{k x ,称之为迭代序列。

我们假定 )(x ϕ在],[b a 连续,且],[}{b a x k ⊂ ,若
}{k x 收敛与x ,则x 必满足方程(1.1.1)
,即为方程(1.1.1)的根。

我们也称方程(2.1.1)的根为函数的)(x ϕ的不动点,因此迭代算是(2.1.2)也称不动点迭代[9]。

2.2 Newton 迭代法的原理
我们从三个角度看Newton 迭代法:
(1)设*x 是方程0)(=x f 的根, 取*x x ≈*将)(x f 在k x 做一阶Taylor 展开:
2))((''2
1))((')()(k k k k x x f x x x f x f x f -+-+=ξ ])[())((')(2k k k k x x o x x x f x f -+-+=
其中ξ在x x k 和之间,于是
])[())((')()(02**k k k k x x o x x x f x f x f -+-+==
))((')(*k k k x x x f x f -+≈
解得)(')(*k k k x f x f x x -≈,用可得取代*1x x k +Newton 迭代法 ,2,1,0,)
(')(1=-=+k x f x f x x k k k k (2.2.1) 只要1C f ∈,每一步迭代都有0
)(lim ,lim ,0)('*==≠∞→∞→x f x x x x f k k k k k 就是方程则而且的根。

迭代函数(2.2.1)即为Newton 迭代格式。

(2) 从几何上看,方程(1.1.1)的根*x 是曲线)(x f y =与x 轴的交点。

设*x 为根x 的一个近似值,用过))(,(k k x f x 处的切线与x 轴的交点来近似(如图 1)
图2.1
因为过点))(,(k k x f x 处的切线方程为
))((')(k k k x x x f x f y -=-
它与x 轴的交点坐标即为 )
(')(1k k k k x f x f x x -=+, ,2,1,0=k (3) 对不动点方程(2.1.3),它导出的迭代过程有可能发散,也可能收敛得非常缓慢.这时,我们有没有办法改进不动点方程,让迭代过程收敛得快一些呢?
我们考虑到x x =和)(x g x =都是不动点方程,他们的加权平均
x x g x h )1()()(λλ-+= (2.3.1)
从而有λλ-+=1)(')('x g x h ,整理有
1)1)('()('+-=x g x h λ (2.3.2)
也是不动点方程,而且)(x h 与)(x g 有完全相同的不动点.适当选取λ的值,可以使发散的迭代过程变得收敛,使收敛慢的迭代过程变得收敛迅速.
由于)(x h 在不动点*x 附近的导数值在很大程度上决定了迭代过程的收敛性.)('x h 的绝对值越小,收敛性越好.因此,选择λ使得0)('=x h .由(2.3.2)式可得到理想的λ值为
)
('11x g -=
λ 将上式代入(2.3.1)得 )
('1)(')()(x g x xg x g x h --= (2.3.3) 在0)('≠x f 的情况下,为求解方程(1.1.1),可以使用不动点方程)(x f x x +=,相应的迭代函数为)()(x f x x g +=.对)(x g 进行迭代
)
(')()]('1[1)]('1[)]([)('1)(')()(x f x f x x f x f x x f x x g x xg x g x h -=+-+-+=--=……(2.3.4) 即为牛顿迭代格式。

从而,我们用Newton 迭代很好的解决了有关方程(1.1.1)的求根问题。

2.3小结:
不动点迭代是迭代法中的重点,它的基本思路与数列收敛差不多。

而Newton 迭代法则是从导数的角度来建立迭代格式,思路新异独到,关键在于起设计思路上。

第三章 Newton 迭代法的收敛性
3.1 Newton 迭代法中不收敛的情况
在Newton 迭代法中,不是所有情况Newton 迭代法都收敛的,若)(x f 不具有连续的有界的,那么Newton 迭代法有可能不收敛[9],如图2所示.
图中,P 为转向点,即在P 点处其导数为零。

并过点k x 的切线交点1+k x ,而过1+k x 的切线交点2+k x 与点k x 刚好重合,从而这个迭代就变成了一个死循环。

图 3.1
另外有一中情况是,其迭代初值的选取不恰当,从而导致迭代发散[12]。

如图3.1
图3.2
因而,在应用Newton 迭代法的时候,对其迭代条件的判断和迭代初值的选取十分重要,这关系到整个迭代成功与否。

++31k k +2k k
3.2 定理证明
定义 1[1] 设)(x ϕ有不动点*x ,如果存在*x 的某个领域R x x x R ∈≤-0,|*:|对任意的δ,迭代法(2.1.2)产生的序列R x k ∈}{,且收敛到*x ,则称迭代法(2.1.2)局部收敛; 定义 2 [9] 设序列.,}{0p p n p p n n n ≠∞=,且对所有的收敛与如
果存在常数01>≥c q 和,使得
c p p p p q
n n n =--=+∞←||||lim 1 (3.2.1) 成立,则称序列阶收敛的是q p n }{,或者说序列.}{p q p n 收敛与以阶称迭代法(2.1.2)是q 阶收敛的,如果相应的迭代序列阶收敛的。

是q x k }{特别地,当101<<=c q 且时称为线性收敛,1>q 称为超线性收敛,2=q 时称为平方收敛或二次收敛。

定理 1[9] 设迭代函数)(x ϕ在其不动点*x 附近有连续的二阶导数,且
1|*)('|<x ϕ
当0*)('',0*)('≠=x x ϕϕ时,不动点迭代(1.2.1)是平方收敛的。

证明 由于1|)('|<x ϕ,由连续性,存在闭领域1L )*,(<及常数δx U ,使得
)*,(,|)('|δϕx U x L x ∈∀<
故按照(2.2.1)产生的迭代序列均收敛于方程(1.1.1)的根*x 。

又由Lagrange 中值定理,
之间和位于**),)(('*)()(x x x x x x ξξϕϕϕ-=-
而**)(x x =ϕ,当)*,(δx U x ∈时,有)*,(δξx U ∈,故有
δϕ≤-<-≤-|*||*||*)(|x x x x L x x
即∈)(x ϕ)*,(δx U .
故定理中的迭代序列}{k x 是局部收敛的
当0*)('=x ϕ时,由Taylor 展开
*,*)(2)(''*)()(21x x x x x x x k k 和位于ηηϕϕϕ-+
==+ 即
21*)(2)
(''*x x x x k -=-+ηϕ (3.2.2)
所以 |2*)(''||*||*|lim 21x x x x x k
k n ϕ=---∞→ (3.2.3) 而又因为0*)(''≠x ϕ,不动点迭代(2.2.1)是平方收敛的。

推论 设迭代函数)(x ϕ在其不动点*x 附近有连续的n 阶导数,且
1|*)('|<x ϕ
当0,0*)(''*)(')()1(≠====-n n x x ϕϕϕϕ 时,不动点迭代(2.2.1)是n 阶收敛的。

证明: 将)(k x ϕ在*x 附近Taylor 展开
)*)((*)(!*)(*)('*)()()(n k n k n k k x x o x x n x x x x x -+-+
+-+=ϕϕϕϕ (3.2.4)
由(2.1.2)知 )(1k k x x ϕ=+,*)(*x x ϕ=

|*)(!*)
(*)*)(('||*)()(||*|)(1n k n k k k x x n x x x x x x x x -++-=-=-+ϕϕϕϕ (3.2.5)
因为0*)(''*)(')1(====-n x x ϕϕϕ ,所以
n k n k x x n x x x *)(!*)
(|*|)(1-=-+ϕ (3.2.6) 而0!*)
()(≠n x n ϕ且为常数的,而
|!*)(||*||*|lim )(1n x x x x x n n k
k k ϕ=--+∞→ (3.2.7) 由收敛阶的定义可知,给收敛是n 收敛的。

定理 2[10]设],[)(b a x f 是上连续可微的凸函数,满足条件
(1);0)()(<b f a f
(2)],[)('b a x f 在上不变号;
则有
(1)存在唯一的],[*b a x ∈,使得0*)(=x f
(2)对任意满足],[0)(00b a x x f ∈>的,记
,2,1,0,)
(')(1=-=+k x f x f x x k k k k (3.2.8) 则}{k x 是单调有界序列;
(3)对上述的序列}{k x ,当∞→k 时有*x x k →; 证明:首先证明(1)
因为0)()(<b f a f 且],[)('b a x f 在不变号,所以)(x f 是],[b a 上的单调函数,由函数的连续性可知,满足方程0)(=x f 的根*x 是唯一存在的。

(2) 由(3.2.8)式有
0))((')(1=-++k k k k x x x f x f (3.2.9)
不妨设0)('≥x f ,有因为
))((')()(k k k x x x f x f x f -+≥
即有0)(1≥+k x f ,故对所有的 ,3,2,1=n 均有
0)(≥n x f (3.2.10)
在由(3.2.9)式可知所有的 ,3,2,1=n 有
01≤-+n n x x (3.2.11)
又因为假定了],[0)(00b a x x f ∈>的,故有对任意的 ,2,1,0=k ,均有
01≤-+k k x x , ,2,1,0=k (3.2.12)
故}{k x 是单调有界序列;
(3) 由(2)可知}{k x 收敛,设其收敛到**x ,由Cuachy 收敛准则有
0||lim 1=-+∞
→k k k x x (3.2.13) 在联系(3.2.8)式有
0*)*(=x f
由(1)可知,若0*)(=x f ,则必有
***x x = (3.2.14)
故对于序列}{k x 有
*lim x x k k =∞
→ 3.3 Newton 迭代法的收敛性分析
(1) 单根情形
设方程(1.1.1)的函数)(x f 在零点 *x 的某领域)*,(δx U 内具有二阶连续导数,且
0*)(''≠x f
由Newton 迭代函数
)
(')()(x f x f x x -
=ϕ (3.3.1) 则 222|
)('|)('')(|)('|)('')(|)('|1)('x f x f x f x f x f x f x f x =--=ϕ (3.3.2) 因此 0*)('=x ϕ,有以上定理可知Newton 迭代法是至少二阶收敛。

(2) 重根情形[9]
以上两种情况,均为Newton 迭代法在方程(1.1.1)在*x 处有单根的情形,要是方程(1.1.1)在*x 处有重根呢?那是不一样的,这里我们考虑方程(1.1.1)在*x 处具有m )2(≥m 重根的情况。

由于方程(1.1.1)在*x 处具有m 重根,方程(1.1.1)的函数)(x f y =可表示成 )(*)()(x g x x x f m -=,且0)(≠x g ,在此,我们设函数)(x f 具有m 连续导数,则有0*)(,0*)(,,0*)(',0*)()1(≠===-x f x f x f x f m m 而 ,从而在*x 附近(*x 除
外)*)('*),(x f x f j 均不为零,Newton 迭代法仍然有效,但起收敛速度明显下降[4]。

在这种情况下,我们可以考虑方程(1.1.1)在*x 附近的m 阶Taylor 展式:
m n m m n m n n x x m f x x m x f x x x f x f x f *)(!)(*)()!1(*)(*)*)(('*)()()(1)1(-+--++-+=--ξ (3.3.3)
1)(2)1(*)()!
1()(*)()!2(*)(*)*)((''*)(')('-----+--++-+=m n m m n m n n x x m f x x m x f x x x f x f x f ξ (3.3.4)
因为)('),(n n x f x f 展开右端的前m 项都为零,而将其带入Newton 迭代格式有
*)()
()()!
1(*))((!*))(()(')()()(1)()(1x x mf f x m x x f m x x f x x f x f x x n m m n m n m m n m n n n n n --=----=-=-+ηξηξ (3.3.5) 其中ηξ与都介于*x x n 和之间。

与是
m mf f x x x x m m n n
n n 11])()(1[lim **lim )()(1-=-=--∞→+∞→ηξ (3.3.6) 当2≥m 时,111<-m
的,所以,Newton 迭代产生的迭代序列}{n x 是线性收敛与*x 的。

显然,在重根的时候,Newton 迭代的速度大大降低了,此时我们就考虑对其进行修正[9],提高其收敛速度。

(1)由以上可知,*x 为函数)(x f 的m )2(≥m 重根。

我们令)
(')()(x f x f x h =
,从而*x 是方程0)(=x h 的单根,故有0*)(≠x h 。

对)(x h 继续实施Newton 迭代
即 222])]
('[)('')(1)]('[)('')()('[)('x f x f x f x f x f x f x f x h -=-=………(2.3.7) 从而得到迭代公式:
)('')()]('[)(')()]('[)('')(1)(')()(')(2
1k k k k k k k k k k k k k k k k x f x f x f x f x f x x f x f x f x f x f x x h x h x x --=--=-=+ (3.3.8)
根据Newton 迭代可知,迭代函数(3.2.3)是二阶收敛的,从而修正后的重根时的Newton 迭代法是二阶收敛。

3.4小结:
Newton 迭代法不是什么情况下都收敛的,更不是什么情况下的收敛速度都相同。

本节着重就这两方面做了详细的分析,并知道改变Newton 迭代的格式还可以提高它的收敛速度。

第四章 两种改进的Newton 迭代法
4.1 改进初值0x 的Newton 下山法[11]
我们知道,Newton 迭代法对迭代初值0x 的选取是有很高要求的,若初值
的选取不好,很可能迭代就发散,为了改善这一苛刻的条件,我们设计了Newton 下山法。

Newton 下山法是在事先没有给出较好的初值0x 情况下, 求方程(1.1.1)
的根*x 的一种修正的Newton 迭代法。

为了防止迭代发散, 对迭代过程须再附加一项要求, 即要求|)(|k x f 单
调减, 也就是
|)(||)(|1k k x f x f <+ (4.1.1)
满足这项要求的算法称为下山法.将Newton 迭代法与下山法结合起来使
用的方法, 称为Newton 下山法, 即将Newton 迭代法的计算结果
)
(')(1k k k k x f x f x x -=+ (4.1.2) 与前一步的近似值k x 适当加权平均作为新的改进值 k k k x x x )1(11λλ-+=++ (4.1.3)
其中称为下)10(≤<λλ山因子, 下山因子的选取应保证)(k x f 单调减, 即
|)(||)(|1k k x f x f <+
(4.1.3)式即为
,3,2,1,)
(')(1=-=+k x f x f x x k k k k λ (4.1.4) 称为Newton 下山法。

下山因子的选择一般采用试算法。

即由迭代即由迭代得k x 后,取不同的λ值试算。

通常取 ,2
1,21,21,132=λ,用(4.1.3)式进行试算1+k x ,在计算)(1+k x f 。

每一次λ的试算,都要判断一次(4.1.1)式。

注:在进行算法设计时,若在整个的迭代过程中,存在一次试算不能满足(4.1.1),则迭代失败,在这种情况下,若对λ的取值不家限制,则整个算法进入死循环,为了防止算法的死循环,我们因对λ的选取作一个停机准则,对此我们可另选取初值0x 进行迭代。

4.2 一种新的Newton 迭代法加速设计[13]
对于非线性方程(1.1.1),若 *x 是该方程的单根,)(x f 在区间],[b a 满足二阶连续可导,由算术平均Newton 迭代法格式[8]
⎪⎪⎩⎪⎪⎨⎧+-=-=+++2))(')('()()(')(111n n n n n n n n n z f x f x f x x x f x f x z (4.2.1) 可知,我们对(4.2.1)中的第一式进行Newton 迭代法,则可的得到如下的迭代格式:
,2,1,0,)(')()(')(1=⎪⎪⎩
⎪⎪⎨⎧-=-=+k y f y f y x x f x f x y k k k k k k k k (4.2.2) 迭代格式要求初始点0x 要充分接近零点*x , )(x f 在区间],[b a 上具有二阶连续导数,且满足0)('≠x f ,则改迭代(4.2.2)三阶收敛。

证明:Newton 迭代格式为
)(')(1k k k k x f x f x x -
=+ 相应的迭代函数为
)
(')()(x f x f x x h -= (4.2.3) 由于Newton 迭代法在单根时至少是二阶收敛。

则有
0*)('',0*)('*,*)(≠==x h x h x x h
迭代格式(4.2.2)的迭代函数:
)
('))(()()(x f x h f x h x g -= (4.2.4) 故有**)(x x g =。

对(4.2.4)两边求导得:
2
))('()('')]([)(')(')]([')(')('x f x f x h f x f x h x h f x h x g --= (3.2.5) 由于0*)('*,*)(==x h x x h ,可知0*)('=x f
对(4.2.5)再求导有:
')')('()]([)]('[)('')]([('2))(')](['')(')(('))(')]([1)(('')(''12-----=x f x h f x f x f x h f x f x h f x h x h x f x h f x h x g 有因为 0*)('=x h ,0)]([=x h f ,故有
0)(''=x g (4.2.6)
由于*x 是方程(1.1.1)的单根,则Newton 迭代法至少二阶收敛。

由此,迭代格式(4.2.2)是三阶收敛
4.3小结:
在这里我们对Newton 迭代法作了一下推广和一种加速迭代的设计,同时也给出了一种在单根时有三阶收敛速度的Newton 迭代。

第五章 Newton 迭代法的应用
5.1 Newton 迭代法的Matlab 实现
下面我们给出Newton 迭代的算法设计步骤[1]:
步骤1 给定)(k x f 与x 的相允许误差1ε、2ε,给定最多迭代次数N ,给定初值0x 。

计算
)(00x f f =,)(''00x f f =。

步骤2 按Newton 迭代格式(1.1.3)计算0
001'f f x x -
=,再计算 )(11x f f = )(''11x f f = 步骤 3 若1x 满足2||εδ<或者11||ε<f ,则终止迭代,以1x 作为根输出。

否则转到步骤4。

这里有
⎪⎩
⎪⎨⎧≥-<-=,
,当,,当1||||||1||||1101101x x x x x x x δ
步骤4 若迭代次数达到N ,仍为求的满足允许误差的根*x ,或者0)('1=x f ,则停止程序:否则以)',,(111f f x 代替)',,(000f f x 转到步骤2。

5.2 数值举例
首先看Newton 迭代法应用的一个经典例子,开方计算。

即:
给定任何正数c ,有
02=-c x (5.2.1)
用Newton 迭代法求方程(5.2.1)近似根(计算机实现),其计算公式为:
)(211k
k k x c x x +=+ (5.2.2)
证明对给定的任意初值00>x 均二次收敛。

分析:由公式(5.2.2),设有迭代函数为
x
c x x 22)(+=ϕ (5.2.3) 从而 2221)('x c x -+=
ϕ 3)(''x
c x =ϕ 则有 0)('',0)('≠==c
c c c ϕϕ 从而对任意的初值0x ,由定理1可知,该迭代法均二阶收敛。

例1用Newton 迭代法解方程 03)(3=--=x x x f
分析:我们知道,该方程在区间]2,1[上有一实根(证明见附录(1))。

下面我们分别选取初值为0、1、2三种情况看看Newton 迭代法的效果
>> newton('x.^3-x-3','3*x.^2-1',0,5e-16,5e-16,1)
函数 x.^3-x-3 的Newton 迭代法
k f(x) dfdx x(k+1)
1 -3 -1 -3.00000000000000
2 -27 26 -1.96153846153846
3 -8.59 10.5 -1.14717596140355
4 -3.36 2.9
5 -0.00657937148071
5 -2.99 -1 -3.00038907407123
6 -2
7 26 -1.96181817566632
7 -8.59 10.5 -1.14743022848160
8 -3.36 2.95 -0.00725624755242
9 -2.99 -1 -3.00047318877322
10 -27 26 -1.96187864636024
11 -8.59 10.5 -1.14748519321677
12 -3.36 2.95 -0.00740250133287
13 -2.99 -1 -3.00049244291695
14 -27 26 -1.96189248824631
15 -8.59 10.5 -1.14749777454454
16 -3.36 2.95 -0.00743597525670
17 -2.99 -1 -3.00049690365354
18 -27 26 -1.96189569508551
19 -8.59 10.5 -1.14750068932977
20 -3.36 2.95 -0.00744373016891
Warning: 在达到指定的次数20 后仍未找到满足指定允许误差的根
>> newton('x.^3-x-3','3*x.^2-1',1,5e-16,5e-16,1)
函数x.^3-x-3 的Newton迭代法
k f(x) dfdx x(k+1)
1 -3
2 2.50000000000000
2 10.1 17.8 1.92957746478873
3 2.25 10.2 1.70786640021144
4 0.274 7.7
5 1.67255847335313
5 0.00634 7.39 1.67170038194364
6 3.69e-006 7.38 1.67169988165733
7 1.26e-012 7.38 1.67169988165716
8 -8.88e-016 7.38 1.67169988165716
ans =1.6717
>> newton('x.^3-x-3','3*x.^2-1',2,5e-16,5e-16,1)
函数x.^3-x-3 的Newton迭代法
k f(x) dfdx x(k+1)
1 3 11 1.72727272727273
2 0.426 7.95 1.67369117369117
3 0.0147 7.
4 1.67170256974750
4 1.98e-00
5 7.38 1.67169988166207
5 3.62e-011 7.38 1.67169988165716
6 -8.88e-016 7.38 1.67169988165716
ans =1.6717
由以上计算可知,当初值改变是,很可能会影响到Newton迭代法的结果和迭代过程。

对于初值去0的情况,我们可以从一下图示(图4)可以体会其发散的理由:
图5.1
并且我们可以看到,在区间]2,1[上函数)
f是连续单调递增的,故其Newton迭代
(x
法有很好的收敛效果。

一般地,如果)
f是二次连续可导且在根*
(x
x附近局部收敛,那么应用Newton 迭代法解决方程)
f的根效果提别的好,其中它有两大优点:一时单调收敛;二
(x
是具有二价收敛;但是,我们从科学计算的角度出发,Newton迭代法的算法并不简单。

因为,它的每一次迭代导函数和导函数的值,这对于科学计算来说是意见难事,就当今科学计算而言,还没有很好的解决计算机求导问题。

这使得Newton 迭代法在其算法的复杂性外,还存在这自身的不足。

但就其整体而言,它的缺点与不足没别要看的太严重。

不论是在科学计算上,还是在求非线性方程的理论上,它都有很重要的地位([2-6])。

应用Newton迭代法,我们必须注意它的收敛性和迭代初值的选取。

就如定理1与定理2所属条件。

其初值选取的重要性可有例1可看出。

在Newton迭代法收敛的情况下,其迭代初值直接影响到迭代的优越与否。

而非线性方程(1.1.1)我们有专们的方法实现确定它的有根区间[9],这也使得我们的Newton迭代法在使用上变得很方便。

5.3小结:
本章针对牛顿迭代存在的不足提出改进后的实际应用,选取适当的例题做一个实际应用上的检验。

对改进后的牛顿迭代在精确度上作实际的验证。

总结
Newton迭代法是一个重要的计算方法和思想。

牛顿迭代法的主要功能:计算方程时可以比较快速方便的计算出来结果但并不影响计算出来结果的精确度,运用于多种工业设计和数学设计方面,其重要性可见一般。

作为迭代法,其收敛性和收敛速度时期关键,收敛性保证了起迭代法的适用性,而收敛速度则是提高计算效率的根本,而适用和效率是我们解决问题的关键。

本文就主要在这两方面,在对前人所作贡献的基础上,集中的给出了Newton迭代法的收敛性和收敛速度,为我们的初学者提供了一个很好的参考平台。

在文中,我们给出了Newton迭代法的不足与优越性,总的来说,Newton迭代法是很重要的迭代法,很多的迭代法都是由它演变或启发而来,对它的全面掌握对我们学习迭代法是有很好的启发作用的。

参考文献
[1]李庆扬,王熊超,易大义.数值分析[M].北京:清华大学出版社,2009
[2] 吴新元,对牛顿迭代法的一个重要修改,应用数学和力学[M],1999,20(8):863-866
[3]J.M.Gutierrez and ,M.A.Hernandez , Newton's method under weak Kantorovich
conditions.IMA J.Numer.Anal.,20(2002)521-532
[4]王兴华,弱条件下Euler族迭代的收敛性[M],中国科学(A辑),30(2000),865-868
[5]王兴华,韩丹夫,Newton迭代的区域估计和点估计[M],计算数学,1(1990),47-53
[6]I.K.Argyros,Newton methods on Banach spaces with a convergence structure and
applications,Comput.Math.Appl.,40(2000)37-48
[7]Kantorovich L.On Newton's method(In Russian)[J],Trudy Mat Inst
Steklov,1949,28(1):104-144
[8]Weerakoon S and Fernando T G.I.A variant of Newton's method with acclerated
third-order convergence[J],Appl.Math.Lett,2000,13(8),87-93
[9]叶兴德,程晓良等著.数值分析基础[M].浙江:浙江大学出版社,2008
[10]甑西丰.实用数值计算方法[M].北京:清华大学出版社,2006
[11]石幸民,郝整清著.基于MATLAB的使用数值计算[M].北京:清华大学出版社·北京交通
大学出版社,2007
[12]姜启源,邢文训等著.大学数学实验[M].北京:清华大学出版社,2006
[13]倪健,马昌凤.解非线性方程牛顿迭代法的一种新的加速技巧[J].广西科学院学报,2010,
26(1):
致谢
衷心感谢我的指导老师令弧荣涛教授和班主任钱淑渠老师的指导。

你们的辛勤指导和谆谆教诲使我完成了本次论文。

你们精益求精的精神和言传身教的作风使我终身受益,我将谨记你们的教诲。

同时也深深的感谢在此过程中关心、支持和帮助过我的同学任明仲、蔡友、王盛兵以及本组其他成员,正是你们的关心和帮助,才使我的学习、工作和生活变得更加精彩,并促成本文的完成。

感谢院系领导及任课教师传授给我的专业知识,使我完成本论文的撰写,感谢他们教我做人的道路,对我终身受益。

感谢我的父母亲,我的兄弟姐妹,是他们的支持和鼓励,使我顺利完成大学四年的求学之路,感谢他们的养育之恩。

相关文档
最新文档