数值分析第2章3-5节.pptVIP

  1. 1、本文档共24页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
2.3 均差与牛顿插值公式 2.3.2 牛顿插值公式 2.4 差分与等距节点插值 2.5 埃尔米特插值 * 2.3.1 均差及其性质 问题:利用插值基函数得到的拉格朗日插值多项式有何优缺点? 优点:结构紧凑,便于理论分析,易于编程求解。 缺点:当插值节点增减时全部插值基函数均要随之变化,整个公式也将发生变化. 问题:如何改进? (3.1) 其中 为待定系数, 确定 . 为了克服这一缺点,可把插值多项式表示为如下便于 计算的形式: 可由 个插值条件 由 , 当 时, 当 时, 依此递推可得到 . 当 时, 推得 推得 由 称 为函数 关 于点 的一阶均差. 称为 的二阶均差. 定义2 一般地,称 为 的 阶均差 (均差也称为差商). 均差有如下的基本性质: 这个性质可用归纳法证明. 1° 阶均差可表示为函数值 的线 性组合, 这性质也表明均差与节点的排列次序无关,称为均差 的对称性. 即 3° 若 在 上存在 阶导数,且节点 这公式可直接用罗尔定理证明. 2° 由性质1°及均差定义可得 即 则 阶均差与导数关系如下: 均差计算可列均差表如下(表2-1). 根据均差定义,把 看成 上一点, 可得 只要把后一式代入前一式,就得到 其中 显然, 由前式确定的多项式 满足插值条件, 且次数不超过 , 称 为牛顿(Newton)均差插值多项式. 系数 就是均差表2-1中加横线的各阶均差,它比拉格朗日插值计算量省,且便于程序设计. 其系数为 它就是形如(3.1)的多项式, 但其更有一般性,它对 是由离散点给出的情形或导数不存在时也是适用的. 由插值多项式惟一性知,此处的插值余项,与拉格朗日插值多项式的余项应该是等价的. 牛顿插值多项式的优点还在于它的递进性,当增加 插值节点时,只要在原来插值多项式的基础上增加一项 即可. 解 首先根据给定函数表造出均差表 给出 的函数表(见表2-2),求3次牛顿插 值多项式,并由此计算 的近似值. 例 27 9 3 1 3 2 1 0 2 6 9 2 4/3 6 18 27 3 2 3 1 1 0 三阶均差 二阶均差 一阶均差 实际应用时经常遇到等距节点的情形,这时插值公式 可以进一步简化,计算也简单得多. 2.4.1 差分及其性质 设函数 在等距节点 上 的值 为已知,这里 为常数,称为步长. 为了得到等距节点的插值公式,先介绍差分的概念. 记号 定义3 分别称为 在 处以 为步长的向前差分,向后差分 及中心差分. 符号 , , 分别称为向前差分算子,向后差分算子 及中心差分算子. 利用一阶差分可定义二阶差分为 一般地可定义 阶差分为 由上面定义的差分算子,易推出差分可用函数值表示,函数值也可用差分表示,均差也可用差分表示。由此可推导出各种形式的等距节点插值公式。 课堂练习:利用下表数据的Newton向前插值公式以及Newton向后插值公式计算 近似值。 1.82212 1.49182 1.22140 1.00000 f(x) 0.6 0.4 0.2 0.0 x 有些实际的插值问题不但要求在节点上函数值相等, 下面只讨论函数值与导数值个数相等的情况. 满足这种要求的插值多项式就是埃尔米特插值多项式. 而且还要求对应的导数值也相等,甚至要求高阶导数也相等.

文档评论(0)

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

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

1亿VIP精品文档

相关文档