第七章 非线性方程的数值解法.pptVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
上一页 下一页 返回 上一页 下一页 返回 例8 现在用函数 构造Steffensen迭代法, k 0 1 … 5 6 xk 1.5 1.4162930 … 1.3247180 1.3247180 可见,Steffensen迭代法对这种不收敛的情形同样有效。 §3 一元方程的常用迭代法 一、Newton迭代法 设x*是方程f (x)=0的实根。取 x0 ? x*,将 f (x)在 x0 做一阶 Taylor展开: ,? 在 x0 和 x 之间。 将 (x* ? x0)2 看成高阶小量,则有: x y x* x0 上一页 下一页 返回 由图示,可见xk+1为函数f (x)在点xk处的切线与横坐标轴的交点,所以,Newton迭代法也称切线法。 Newton迭代格式 上一页 下一页 返回 与上一节例7相比,牛顿法的收敛速度快很多。 例9 上一页 下一页 返回 牛顿迭代法的流程图 注:牛顿迭代法的收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 上一页 下一页 返回 定理 上一页 下一页 返回 (收敛的充分条件)设 f ?C2[a, b],若 (1) f (a) f (b) 0; (2) 在整个[a, b]上 f ’’ 不变号且 f ’(x) ? 0; (3) 选取 x0 ? [a, b] 使得 f (x0) f ”(x0) 0; 则牛顿迭代法产生的序列{ xk } 收敛到f (x) 在 [a, b] 的唯一根。 定理 有根 根唯一 产生的序列单调有界,保证收敛。 保证f (x) 上凸或下凸 上一页 下一页 返回 10 Q1: 若    ,牛顿迭代法是否仍收敛? 设 x* 是 f 的 m 重零点, 则: 且 A1: 有局部收敛性,但仅为线性收敛. 下面介绍计算重根的牛顿迭代法 上一页 下一页 返回 * 第七章 非线性方程的数值解法 求 f (x) = 0 的根 第一节 方程求根的二分法 第二节 一元方程的不动点迭代法 第三节 一元方程的常用迭代法 上一页 下一页 返回 本章要讨论的问题: (非线性)方程f (x)=0的数值解法 首先引入定义: 1. f (x)=0的解x*称为方程f (x)=0的根或函数f (x) 的零点. 若f (x)=(x-x*)m g(x) , g (x*) ≠0 ,其中m为正整数, 则称 x*为方程f (x)=0的m重根,或称 x*为函数 f (x)的m重零点. 上一页 下一页 返回 §1 方程求根的二分法 方程求根步骤: ①根的隔离 确定根所在的区间[a, b], 使在[a, b]内至少有方程的一个根. 有根区间 ②近似根的精确化 已知根的一个近似值后, 构造某种算法, 将此近似值精确化, 使其满足给定的精度要求. 越小越好 上一页 下一页 返回 下面介绍方程求根的二分法 在确立了有根区间[a, b]后, 需要逐步缩小有根区间. 取[a, b]的中点x0=(a+b)/2,将区间一分为二.若 f(x0)=0,则x0 就是方程的根,否则判别根 x* 在x0 的左侧还是右侧. 不论出现哪种情况, (a1, b1)均为新的有根区间, 它的长度只有原有根区间长度的一半, 达到了压缩有根区间的目的. 上一页 下一页 返回 重复以上过程, 继续进行二分, 可得一系列有根区间 由于每个小区间都是有根区间, 所以这个点就是所求方程的根. 同时, 在每次二分时, 所求出的中点 形成一个无穷数列 这个数列必定收敛于 x*. 上一页 下一页 返回 a b x0 x1 a1 b2 When to stop? x* b1 上一页 下一页 返回 a2 误差 分析: 第1步产生的 有误差 第 k 步产生的 xk-1 有误差 对于给定的精度 ? ,可估计二分法所需的步数 k : ① 算法简单; ② 对f (x) 要求不高(只要连续即可) ,收敛

文档评论(0)

考试资料文库 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档