《牛顿法,弦截法》PPT课件.ppt_第1页
《牛顿法,弦截法》PPT课件.ppt_第2页
《牛顿法,弦截法》PPT课件.ppt_第3页
《牛顿法,弦截法》PPT课件.ppt_第4页
《牛顿法,弦截法》PPT课件.ppt_第5页
已阅读5页,还剩17页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

7.4 -7.5 牛顿法及其推广,/* Newton Method */,一、牛顿迭代法的公式,二、牛顿迭代法的改进与推广,原理:将非线性方程线性化泰勒展开 /* Taylors expansion */,取 x0 x* ,将 f(x*) 在 x0 做一阶泰勒展开:, 在 x0 和 x*之间。,将 (x* x0)2看成高阶小量,则有:,一、牛顿迭代法的公式,线性 /* linear */,x0,只要 每一步迭代都有f ( xk ) 0, 而且 ,则 x*就是 f 的根。,牛顿迭代法的基本思想,将非线性方程 f(x)=0 的求根问题归结为计算一系列线性方程的求根问题。,牛顿迭代法的计算步骤,(1)给出初始近似根 x0 及精度;,(3)若 ,转向(4), 否则 ,转向(2);,(4)输出满足精度的根 x1 ,结束。,(2)计算,例,用牛顿迭代法求方程,在 x=0.5 附近的根。取,解,其牛顿迭代公式为,取初值 x0=0.5 ,迭代结果见下表,易见,故,k 0 1 2 3 xk 0.880000 0.884688 0.884675 0.884675,例 2 计算 的近似值, =10-6 x0=0.88,解:令 x=,问题转化为求f(x)=x2-0.78265=0 的正根,由牛顿迭代公式,xk+1= xk-(xk)/(xk)= xk/2+0.78265/2xk,迭代结果,满足了精度要求,故,0.884675,设 f C2a, b,若 x* 为 f (x) 在a, b上的根,且 f (x*) 0,则存在 x* 的邻域,Newtons Method产生的序列 xk 收敛到x*,且满足,使得任取初值, Newtons Method 有 ,只要 就有 p 2。重根是线性收敛的。,证明: Newtons Method 事实上是一种特殊的不动点迭代,其中,局部收敛,则,由泰勒展开:,在单根 /*simple root */ 附近收敛快,注: Newtons Method 收敛性依赖于x0 的选取。,x*,x0,x0,x0,注 (1) 牛顿法要求初值充分接近根以保证局部收敛性。,(2)牛顿迭代法的主要优点是收敛较快,是平方收敛的缺点是公式中需要求 f(x) 的导数。若 f(x)比较复杂,则使用牛顿公式就大为不便。, 重根 /* multiple root */ 加速收敛法:,问题1: 若 ,Newtons Method 是否仍收敛?,设 x*是 f 的 n 重根,则:,因为 Newtons Method 事实上是一种特殊的不动点迭, 其中,二、牛顿迭代法的改进与推广 /* improvement and generalization */,且,K1: 有局部收敛性,但重数 n 越高,收敛越慢。,则,问题2: 如何加速重根的收敛?,K2: 将求 f 的重根转化为求另一函数的单根。,?,令,则 f 的重根 = 的单根。, 下山法 /* Descent Method */ Newtons Method 局部微调:,原理:若由 xk 得到的 xk+1 不能使 | f | 减小,则在 xk 和 xk+1 之间找一个更好的点 ,使得,注: = 1 时就是Newtons Method 公式。 当 = 1 代入效果不好时,将 减半计算。 。, 弦截法 /* Secant Method */,Newtons Method 一步要计算 f 和 f ,相当于2个函数值,比较费时。现用 f 的值近似 f ,可少算一个函数值。,切线 /* tangent line */,割线 /* secant line */,切线斜率 割线斜率,需要 2 个初值 x0 和 x1 。,收敛比Newtons Method 慢,且对初值要求同样高。,弦截法与牛顿法相比较,相同之处:都是线性化方法,不同之处:牛顿法在计算xk+1时只用到前一步的值xk,故这种方法称为单点迭代法。,而弦截法在求xk+1时要用到前两步的值xk和xk-1,因此这种方法称为多点迭代法。,有关弦截法的收敛速度,与牛顿法相比,弦截法的收敛速度也是比较快的。可以证明,弦截法具有超线性收敛速度,收敛阶为,即,例,用弦截法求方程,在x=0.5附近的根。取,解,取x0=0.5,x1=0.6作为初始近似根,令,其弦截法

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论