数值计算方法-插值法.pptxVIP

  1. 1、本文档共52页,可阅读全部内容。
  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文档。上传文档
查看更多
数值计算方法-插值法

;;定理1 n次代数插值问题的解是存在且惟一的 ; 这是一个关于待定参数 的n+1阶线性方 程组,其系数矩阵行列式为 ;;定理2 设f(x)在?a, b?有n+1阶导数, x0, x1,…, xn 为 ?a, b?上n+1个互异的节点, p(x)为满足 p(xi) = f(xi) (i=1,2, …, n) 的n 次插值多项式,那么对于任何x ? ?a, b?有 插值余项;;拉格朗日插值多项式 两个插值点可求出一次插值多项式,而三 个插值点可求出二次插值多项式。插值点增加到n+1 个时,也就是通过n+1个不同的已知点 ,来构造一个次数为n的代数多项式P(x)。与推导抛物插值的基函数类似,先构造一个特殊n次多项式 的插值问题,使其在各节点 上满足 ;其中 为待定常数。由条件 ,可求得 ;以n+1个n次基本插值多项式 为基础,就能直接写出满足插值条件 的n次代数插值多项式。 事实上,由于每个插值基函数 都是n次值多项式,所以他们的线性组合;§3 均差与牛顿插值多项式 拉格朗日插值多项式结构对称,使用方便。但由于是用基函数构成的插值,这样要增加一个节点时,所有的基函数必须全部重新计算,不具备承袭性,还造成计算量的浪费。这就启发我们去构造一种具有承袭性的插值多项式来克服这个缺点,也就是说,每增加一个节点时,只需增加相应的一项即可。这就是牛顿插值多项式。 ;由线性代数知,任何一个不高于n次的多项式, 都可以表示成函数;;3.1差商及其性质;f[xi,xj,xk]是指;差商表;xi;在n+1个节点处各阶差商的计算方法;这个性质可用数学归纳法证明(用Lagrange插值多项式比较最高项系数来得到);f[x0 , x1]=;性质3 若f[x, x0, x1 , …, xk ]是 x 的 m 次多项式, 则 f[x, x0, x1 ,…, xk , xk+1]是 x 的 m-1 次多项式 证:由差商定义 ;4.4 .1 差商及其性质;性质5 k阶差商 和k阶导数之间有下 列关系 这个性质可直接用罗尔(Rolle)定理证明(或以下方法即余项方法) ;牛顿(Newton)插值多项式 ;;为牛顿插值多项式的误差。由插值多项式的存在惟一性定理知,满足同一组插值条件的拉格朗日插值多项式P(x)与牛顿插值多项式Nn(x)实际上是同一???多项式,仅是同一插值多项式的不同表达形式而已,因此得到牛顿插值多项式的误差与拉格朗日插值多项式的误差也完全相等。故有 ;;f(x)=f(x0)+(x- x0)f[x1,x0]+(x- x0)(x-x1)f[x1,x0,x];;xi;xi;牛顿插值余项; Rn(n)(x) =f (n)(x)- Pn(n)(x);增加新节点x,并且f(x)为(n+1)阶可导时,有;4.4 .1 差商及其性质;4.4 .1 差商及其性质;例2.15 求 并估计其误差;f 3(x) =;差分与等距节点插值;x;;;例2.16 计算 f (x) = x3在等距节点0,1,2,3, 4上的各 阶差分值;牛顿前插公式;;x;向后差分;同样利用向后差分可以得到牛顿向后插值公式 其中 ,公式 称之为牛顿向后插值公式余项。

文档评论(0)

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

分享好文档!

1亿VIP精品文档

相关文档