第5章4节牛顿法.ppt_第1页
第5章4节牛顿法.ppt_第2页
第5章4节牛顿法.ppt_第3页
第5章4节牛顿法.ppt_第4页
第5章4节牛顿法.ppt_第5页
已阅读5页,还剩16页未读 继续免费阅读

下载本文档

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

文档简介

1、1,4 牛顿法,1 牛顿法及其收敛性,设已知方程 有近似根 (假定 ),,将函数 在点 展开,有,于是方程 可近似地表示为,2,这是个线性方程,记其根为 ,,则 的计算公式为,这就是牛顿(Newton)法.,牛顿法的几何解释.,方程 的根 可解释为,曲线 与 轴的交点的横坐标,图5-3,(图5-3).,3,设 是根 的某个近似值,过曲线 上横坐标为 的点 引切线,并将该切线与 轴的交点的横坐标 作为 的新的近似值.,由于这种几何背景,牛顿法也称切线法.,由定理4,可以直接得到牛顿法的收敛性.,4,由于,假定 是 的一个单根,,即 ,,则由上式知,牛顿法的迭代函数为,又因,5,可得,例7,解,取

2、迭代初值 ,迭代结果列于表5-5中.,用牛顿法解方程,这里牛顿公式为,所给方程实际上是方程 的等价形式.,在根 的邻近是平方收敛的.,牛顿法,6,牛顿法的计算步骤:,步骤1 准备,选定初始近似值 ,,步骤2 迭代,迭代一次,,得新的近似值 ,,若用不动点迭代到同一精度,可见牛顿法的收敛速度是很快的.,要迭代17次.,计算,按公式,计算,7,则终止迭代,以 作为所求的根;,此处 是允许误差,而,其中 是取绝对误差或相对误差的控制常数,,步骤4 修改,如果迭代次数达到预先指定的次数 ,,步骤3 控制,如果 满足 或 ,,否则转步骤4.,一般可取,8,或者 ,则方法失败;,否则以 代替 转步骤2继续

3、迭代.,9,2 牛顿法应用举例,对于给定的正数 ,应用牛顿法解二次方程,可导出求开方值 的计算程序,这种迭代公式对于任意初值 都是收敛的.,10,解,取初值 ,对,按牛顿迭代公式迭代3次,便得到精度为 的结果,例8,求 .,(见表5-6).,11,3 牛顿下山法,牛顿法的优点是收敛快,缺点一是每步迭代要计算 及 ,计算量较大且有时 计算较困难,,为克服这两个缺点,通常可用下述方法.,12,牛顿下山法.,牛顿法收敛性依赖初值 的选取.,如果 偏离所求根 较远,则牛顿法可能发散.,例如,用牛顿法求方程,在 附近的一个根 .,设取迭代初值 ,用牛顿法公式,计算得,13,迭代3次得到的结果 有6位有效

4、数字.,但如果改用 作为迭代初值,则依牛顿法公式 迭代一次得,这个结果反而比 更偏离了所求的根 .,为了防止迭代发散,对迭代过程再附加一项要求,即 具有单调性:,满足这项要求的算法称下山法.,14,将牛顿法与下山法结合起来使用,即在下山法保证函 数值稳定下降的前提下,用牛顿法加快收敛速度.,将牛顿法的计算结果,与前一步的近似值 适当加权平均作为新的改进值,其中 称为下山因子,,称为牛顿下山法.,即为,15,选择下山因子时从 开始,逐次将 减半进行试算,,若用此法解方程 ,当 时可求得,直到能使下降条件成立为止.,,它不满足下降条件.,通过 逐次取半进行试算,当 时可求得,此时有 ,,显然 .,

5、而,由 计算 时 , 均能使下降条件 成立. 计算结果如下 :,16,即为 的近似.,则可得到 ,从而使 收敛.,一般情况只要能使下降条件成立,,17,4 重根情形,设 ,整数 ,,只要 仍可用牛顿法计算,此时迭代函数,的导数为,此时有,且 ,所以牛顿法求重根只是线性收敛.,18,则 .,求 重根,则具有2阶收敛,但要知道 的重数 .,构造求重根的迭代法,还可令 ,,若 是 的 重根,则,若取,用迭代法,修正法1,修正法2,19,从而可构造迭代法,它是2阶收敛的.,故 是 的单根.,对它用牛顿法,其迭代函数为,20,例9,方程 的根 是二重根,,解,(1) 牛顿法,用上述三种方法求根.,先求出三种方法的迭代公式:,(2) 用修正法1),(3) 用修正法2),取初值

温馨提示

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

评论

0/150

提交评论