数值分析7-4,5(牛顿法,弦截法)PPT课件_第1页
数值分析7-4,5(牛顿法,弦截法)PPT课件_第2页
数值分析7-4,5(牛顿法,弦截法)PPT课件_第3页
数值分析7-4,5(牛顿法,弦截法)PPT课件_第4页
数值分析7-4,5(牛顿法,弦截法)PPT课件_第5页
已阅读5页,还剩18页未读 继续免费阅读

下载本文档

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

文档简介

7.4-7.5牛顿法及其推广,/*NewtonMethod*/,一、牛顿迭代法的公式,二、牛顿迭代法的改进与推广,原理:将非线性方程线性化泰勒展开/*Taylorsexpansion*/,取x0x*,将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,迭代结果见下表,易见,故,k0123xk0.8800000.8846880.8846750.884675,例2计算的近似值,=10-6x0=0.88,解:令x=,问题转化为求f(x)=x2-0.78265=0的正根,由牛顿迭代公式,xk+1=xk-(xk)/(xk)=xk/2+0.78265/2xk,迭代结果,满足了精度要求,故,0.884675,设fC2a,b,若x*为f(x)在a,b上的根,且f(x*)0,则存在x*的邻域,NewtonsMethod产生的序列xk收敛到x*,且满足,使得任取初值,NewtonsMethod有,只要就有p2。重根是线性收敛的。,证明:NewtonsMethod事实上是一种特殊的不动点迭代,其中,收敛,则,由泰勒展开:,在单根/*simpleroot*/附近收敛快,注:NewtonsMethod收敛性依赖于x0的选取。,x*,x0,x0,x0,注(1)牛顿法要求初值充分接近根以保证局部收敛性。,(2)牛顿迭代法的主要优点是收敛较快,是平方收敛的缺点是公式中需要求f(x)的导数。若f(x)比较复杂,则使用牛顿公式就大为不便。,重根/*multipleroot*/加速收敛法:,问题1:若,NewtonsMethod是否仍收敛?,设x*是f的n重根,则:,因为NewtonsMethod事实上是一种特殊的不动点迭,其中,二、牛顿迭代法的改进与推广/*improvementandgeneralization*/,且,K1:有局部收敛性,但重数n越高,收敛越慢。,则,问题2:如何加速重根的收敛?,K2:将求f的重根转化为求另一函数的单根。,?,令,则f的重根=的单根。,下山法/*DescentMethod*/NewtonsMethod局部微调:,原理:若由xk得到的xk+1不能使|f|减小,则在xk和xk+1之间找一个更好的点,使得,注:=1时就是NewtonsMethod公式。当=1代入效果不好时,将减半计算。,弦截法/*SecantMethod*/,NewtonsMethod一步要计算f和f,相当于2个函数值,比较费时。现用f的值近似f,可少算一个函数值。,切线/*tangentline*/,割线/*secantline*/,切线斜率割线斜率,需要2个初值x0和x1。,收敛比NewtonsMethod慢,且对初值要求同样高。,弦截法与牛顿法相比较,相同之处:都是线性化方法,不同之处:牛顿法在计算xk+1时只用到前一步的值xk,故这种方法称为单点迭代法。,而弦截法在求xk+1时要用到前两步的值xk和xk-1,因此这种方法称为多点迭代法。,有关弦截法的收敛速度,与牛顿法相比,弦截法的收敛速度也是比较快的。可以证明,弦截法具有超线性收敛速度,收敛阶为,即,例,用弦截法求方程,在x=0.5附近的根。取,解,取x0=0.5,x1=0.6作为初始近似根,令,其弦截法迭代公式为,迭代结果见下表,易见,故,取初值x0=0.5,牛顿迭代结果见下表,抛物线法,本质

温馨提示

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

评论

0/150

提交评论