定理Newton迭代法xk+1=xkf(xk)f′(xk)在f(x)=0的单根x临近为平方收敛.
YongChengComputingMethodsCh04方程求根的迭代法20/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCTNewton迭代法分析Newton迭代法优缺点:Newton迭代法逻辑结构简单、收敛速度很快(平方收敛),但它通常依赖初值x0的选取,如果初值x0选择不当,将导致迭代发散或产生无限循环;此外,每一步迭代都需要计算导数值f′(x),有时计算f′(x)是不方便的.
基于这两点,产生了几种Newton迭代法的变形形式.
1牛顿下山法;2弦截法;3快速弦截法;YongChengComputingMethodsCh04方程求根的迭代法21/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT牛顿法例子一般地说,Newton法的收敛性依赖于初值x0的选取,如果x0偏离解x较远,则Newton法可能发散或产生无限循环.
例题:用Newton求方程x3x1=0在x=1.
5附近的一个根.
解:因f′(x)=3x21可得牛顿迭代公式:xk+1=xkf(xk)f′(xk)=xkx3kxk13x2k1YongChengComputingMethodsCh04方程求根的迭代法22/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT牛顿法例子(续)分别取x0=1.
5和x0=0.
6,计算结果如下表.
kxkxk01.
50.
611.
3478317.
9000021.
3252011.
9468031.
324727.
98551941.
32472由上表可知道,当x0=0.
6时结果偏离所求的根,不收敛(发散)或收敛较慢.
YongChengComputingMethodsCh04方程求根的迭代法23/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCTNewton下山法为了防止迭代发散,通常对迭代过程再附加一项要求,即保证函数值单调下降:|f(xk+1)|<|f(xk)|满足这项要求的算法称下山法.
将Newton法与下山法结合使用,即在下山法保证迭代函数值稳定下降的前提下,用Newton法加快速度,即可得到如下Newton下山法:xk+1=xkλf(xk)f′(xk)其中0<λ<1,称下山因子,在迭代过程中通过适当地选取λ以使下山条件|f(xk+1)|<|f(xk)|满足.
下山因子的选择是个逐步探索的过程,从λ=1开始反复将因子λ的值减半进行试算,一旦单调条件|f(xk+1)|<|f(xk)|满足,则称为"下山成功".
反之,如果在上述过程中找不到使下山条件|f(xk+1)|<|f(xk)|成立的下山因子λ,则称"下山失败",这时需另选初值x0重算.
YongChengComputingMethodsCh04方程求根的迭代法24/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT下山法例子例题:使用下山法求方程f(x)=x3–x–1=0的根,取x0=0.
6.
解:迭代公式如下:xk+1=xkλf(xk)f′(xk)=xkλx3kxk13x2k1牛顿下山法的计算结果:kλxk010.
611251.
14063211.
36681311.
32628411.
32472YongChengComputingMethodsCh04方程求根的迭代法25/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT弦截法牛顿法要计算f′(x),现用f(x)的值近似f′(x):认为切线斜率近似等于割线斜率.
f′(xk)≈f(xk)f(x0)xkx0xk+1=xkf(xk)f(xk)f(x0)(xkx0)迭代函数为:φ(x)=xf(x)f(x)f(x0)(xx0)单点弦截法为线性收敛.
YongChengComputingMethodsCh04方程求根的迭代法26/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT弦截法几何意义YongChengComputingMethodsCh04方程求根的迭代法27/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT快速弦截法快速弦截法也称为两点弦截法.
认为切线斜率近似等于割线斜率.
f′(xk)≈f(xk)f(xk1)xkxk1xk+1=xkf(xk)f(xk)f(xk1)(xkxk1)快速弦截法需要2个初值x0和x1,其收敛阶1.
618.
YongChengComputingMethodsCh04方程求根的迭代法28/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT快速弦截法例子例题:使用Newton法和快速弦截法求方程xex–1=0的根.
解:使用Newton法和快速弦截法,迭代公式分别如下:xk+1=xkf(xk)f′(xk)=xkxkexk1exk+xkexkxk+1=xkf(xk)(xkxk1)f(xk)f(xk1)kxkxk00.
50.
510.
571020.
620.
567160.
56531530.
567140.
56709440.
567140.
567143YongChengComputingMethodsCh04方程求根的迭代法29/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT埃特金迭代公式xk+1=φ(xk)xk+1=φ(xk+1)xk+1=xk+1(xk+1xk+1)2xk+12xk+1+xkYongChengComputingMethodsCh04方程求根的迭代法30/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT本章小结迭代法思想;开方法;Newton法;Newton法的改进;迭代过程的加速;YongChengComputingMethodsCh04方程求根的迭代法31/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT练习题1编程实现开方法.
2编程实现Newton法.
YongChengComputingMethodsCh04方程求根的迭代法32/33迭代法开方法Newton法改进的牛顿法埃特金方法BUCT谢谢!
Author:ChengYongAddress:Dept.
ofComputerBeijingUniversityofChemicalTechnologyBeijing,100029,ChinaEmail:buctcourse@163.
comYongChengComputingMethodsCh04方程求根的迭代法33/33
WordPress专业外贸企业网站搭建模版,特色专业外贸企业风格 + 自适应网站开发设计 通用流行的外贸企业网站模块 + 更好的SEO搜索优化和收录 自定义多模块的产品展示功能 + 高效实用的后台自定义模块设置!采用标准的HTML5+CSS3语言开发,兼容当下的各种主流浏览器: IE 6+(以及类似360、遨游等基于IE内核的)、Firefox、Google Chrome、Safari、Opera...
专心做抗投诉服务器的VirtVPS上线瑞士机房,看中的就是瑞士对隐私的保护,有需要欧洲抗投诉VPS的朋友不要错过了。VirtVPS这次上新的瑞士服务器采用E-2276G处理器,Windows/Linux操作系统可选。VirtVPS成立于2018年,主营荷兰、芬兰、德国、英国机房的离岸虚拟主机托管、VPS、独立服务器、游戏服务器和外汇服务器业务。VirtVPS 提供世界上最全面的安全、完全受保护和私...
hostkvm本月对香港国际线路的VPS、韩国CN2+bgp线路的VPS正在做7折终身优惠,对日本软银线路、美国CN2 GIA线路、新加坡直连线路的VPS进行8折终身优惠促销。所有VPS从4G内存开始支持Windows系统,当然主流Linux发行版是绝对不会缺席的!官方网站:https://hostkvm.com香港国际线路、韩国,7折优惠码:2021summer日本、美国、新加坡,8折优惠码:2...