《计算方法课件第七章数值积分与数值微分》教学课件.pptVIP

《计算方法课件第七章数值积分与数值微分》教学课件.ppt

  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文档。上传文档
查看更多
由于高斯型求积公式的代数精度为2n+1,故 证毕 高斯型求积公式具有代数精度高、且总是收敛、稳定的优点。但当求积节点数增加时,前面的函数值不能在后面利用。因此,有时也可以将区间分化成若干个小区间,在每个小区间上应用低阶的Gauss型求积公式,即复化高斯求积公式。 §7.5 数值微分 以离散数据 近似表达 在节点 处的微分,通常称这类问题为数值微分。 一、Taylor展开法 根据Taylor展开式可得 则有: 类似地,由 可得下面的中点公式: 中点公式: 展开到3阶可得: 给出列表函数 ,可建立插值多项式 ,取 作为 的近似函数,则称为 插值型求导公式。 二、插值函数法 确定节点, 上的导数值,有余项 为讨论方便,假定所给节点是等距的。 1.一阶两点公式 2.一阶三点公式 3.二阶三点公式 Remark1:在数值微分计算中,并非步长越小精度越高。这是因为数值微分对舍入误差非常敏感,它随步长h的缩小而增大,导致计算不稳定。 Remark2:在数值微分计算中,当插值多项式收敛到函数f(x)时,P?n(x)不一定收敛到f ?(x)。 Remark3:为了避免上述问题,可以用样条插值函数的导函数来代替f(x)的导函数。 谢谢! 故由介值定理,一定在(a,b)有一点 使 Remark:若 ,则有误差估计式 证毕 2.复化辛浦生公式 将[a,b]等份成N个子区间[x2k,x2k+2](k=0,1,…,N-1),子区间长度 由 3.复化柯特斯公式 将[a,b]等份成N个子区间[x4k,x4k+4](k=0,1,…,N-1),子区间长度 由 例子1 若取9个节点,用复化梯形公式、复化辛浦生公式和复化柯特斯公式计算积分,其步长以及与9个节点所对应的求积系数分别是多少? 解:复化梯形公式:N=8,h=(b-a)/8,对应的求积系数为1、2、2、2、2、2、2、2、1。 复化辛浦生公式:N=4,h=(b-a)/4,对应的求积系数为1、4、2、4、2、4、2、4、1。 复化柯特斯公式:N=2,h=(b-a)/2,对应的求积系数为7、32、12、32、14、32、12、32、7。 # 用积分 计算ln2,要使所得积分近似值具有5位有效数字。问用复化梯形公式,复化Simpson公式时,至少要取多少个节点? 例子2 解:由 且 故,计算ln2时,要使误差不超过 也即计算2ln2,其误差不超过 。 即 其中 故区间应取671个,节点至少应取672。 其中 由 故区间N应取22,即45个节点 。 # 二、区间逐次分半求积法 ? 1.误差的事后估计法 复化求积公式是提高精度的一种有效方法,但在使用复化求积公式之前,必须根据复化求积公式的余项进行先验估计,以确定节点数目,从而确定合适的等分步长。因为余项表达式中包含了被积函数的导数,而估计各阶导数的最大值往往是很困难的,且估计的误差上界往往偏大。所以实际中,常常使用“事后估计误差”的方法,通过区间的逐次分半,在步长逐次分半的过程中,反复利用复化求积公式进行计算,查看相继两次计算结果的差值是否达到要求,直到所求得的积分值满足精度要求。 该法也称为步长自动选择的变步长求积法。 对梯形公式,假定区间分为N等份时,由公式 算出的积分近似值为TN ,因而有 再把各个小区间分别对分,得积分的近似值为T2N,则积分值为 其中 假定f??(x)有[a,b]上变化不大,即有 , 则 上式可改写为 计算时只需检验 是否满足?若不满足,则再把区间分半进行计算,直到满足要求为止 。 类似的,还可以得到下面的结论: 对于辛浦生公式,假定 在[a,b]上变化不大,则有 对于Cotes公式,假定 在[a,b]上变化不大,则有 ? 2.区间逐次分半的梯形公式 据此我们得到复化梯形公式区间逐次分半时的递推计算公式: 计算时只需检验 是否满足?若不满足,则再把区间分半进行计算,直到满足要求为止 。 §7.3 Romberg求积法 一、对低精度公式经过组合构造高精度公式 从SN及T2N,TN的计算公式可验证得到: 事实上, (1) (2) (3) 类似地,可以证明: 这个公式(3)称为Romberg公式。 由(1)(2)(3)组成的方法称为Romberg算法。

文档评论(0)

秦圈圈 + 关注
实名认证
文档贡献者

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

1亿VIP精品文档

相关文档