计算方法教学课件第六章.数值积分和数值微分.ppt

计算方法教学课件第六章.数值积分和数值微分.ppt

  1. 1、本文档共109页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
例6.14 用复化梯形公式计算定积分 才能使误差不超过 解:取 ,则 ,又区间长度b-a=1,对 复化梯形公式有余项 即 ,n≥212.85,取n=213,即将区间 [0,1]分为213等份时,用复化梯形公式计算误差 不超过 。 问区间[0,1]应分多少等份 6.6 龙贝格(Romberg)求积法 复化求积方法对于提高计算精度是行之有效的方法,但复化公式的一个主要缺点在于要先估计出步长。若步长太大,则难以保证计算精度,若步长太小,则计算量太大,并且积累误差也会增大。在实际计算中通常采用变步长的方法,即把步长逐次分半,直至达到某种精度为止。 6.6.1变步长的梯形公式 变步长复化求积法的基本思想是在求积过程中,通过对计算结果精度的不断估计,逐步改变步长(逐次分半),直至满足精度要求为止。即按照给定的精度实现步长的自动选取。 设将积分区间[a,b]n等分,即分成n个子区间,一共有n+1个节点,即x=a+kh, k=0,1,…,n,步长 。对于某个子区间 ,利用梯形公式计算积分近似值有 对整个区间[a,b]有 将子区间 再二等份,取其中点 作新节点,此时区间数增加了一倍为2n,对某个子区间 ,利用复化梯形公式计算其积分近似值 。 对整个区间[a,b]有 比较 和 有 (6.7) (6.7)式称为变步长梯形公式 当把积分区间分成n等份,用复化梯形 公式计算积分I的近似值 时,截断误差为 若把区间再分半为2n等份,计算出定积分 的近似值 ,则截断误差为 当 在区间[a,b]上变化不大时,有 所以 可见,当步长二分后误差将减至 ,将 上式移项整理,可得验后误差估计式 上式说明,只要二等份前后两个积分值 和 相当接近,就可以保证计算结果 的误差很小,使 接近于积分值I。 6.6.2 变步长的梯形求积算法实现 (1)变步长的梯形求积法的计算步骤 ① 变步长梯形求积法。它是以梯形求积公式为基础,逐步减少步长,按如下递推公式求二分后的梯形值 其中Tn和T2n分别代表二等分前后的积分值 ② 如果 , (ε为给定的误差限 ) 则T2n作为积分的近似值, 否则继续进行二等分, 即 转 ①再计算,直到满足所要求的精度为止,最终取二分后的积分值T2n 作为所求的结果 (2)变步长梯形公式的流程图 (3)程序实现(见附录A A-13 变步长梯形求积法) 例6.15 用变步长梯形求积法计算定积分 解: 先对整个区间?0,1?用梯形公式,对于 所以有 然后将区间二等份,由于 ,故有 进一步二分求积区间,并计算新分点上的函数值 有 这样不断二分下去,计算结果如P139列表所示。积分的准确值为0.9460831,从表中可看出用变步长二分10次可得此结果。 6.6.3 龙贝格求积公式 变步长梯形求积法算法简单,但精度较差,收敛速度较慢,但可以利用梯形法算法简单的优点,形成一个新算法,这就是龙贝格求积公式。龙贝格公式又称逐次分半加速法。 根据积分区间分成n等份和2n等份时的误差估计式(6.8)可得 所以积分值 的误差大致等于 ,如果用 对 进行修正时, 与 之和比 更接近积分真值,所以可以将 看成是对 误差的一种补偿,因此可得到具有更好效果的式子. (6.9) 考察 与n等份辛卜生公式 之间的关系。将 复化梯形公式 梯形变步长公式 代入(6.9) 表达式得 故 这就是说,用梯形法二分前后两个积分值 和 作线性组合,结果却得到复化辛卜生公式计算得到 的积分值 。 再考察辛卜生法。其截断误差与 成正比,因此,如果将步长折半,则误差减至 ,即有 由此可得 可以验证,上式右端的值其实等于Cn,就是说,用辛卜生公式二等份前后的两个积分值Sn和S2n 作线性组合后,可得到柯特斯公式求得的积分值Cn,即有 (6.11) 用同样的方法,根据柯特斯公式的误差公式,可进一步导出龙贝格公式 (6.12) 在变步长的过程中运用(6.10)、(6.11)和(6.12),就能将粗糙的梯形值Tn逐步加工成精度较高的辛卜生值Sn、柯特斯值Cn和龙贝格值Rn或者说,将收敛缓慢的梯形值序列Tn加工成收敛迅速的龙贝格值序列Rn,这种加速方法称为龙贝格

文档评论(0)

autohhh + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档