湖南科技大学数值计算方法课件.pptVIP

下载本文档

8
0
约 39页
2017-08-21 发布于重庆
举报
版权申诉

湖南科技大学数值计算方法课件.ppt

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

湖南科技大学数值计算方法课件.ppt

上一页下一页返回 ? 差分与差商的关系又上一页下一页返回差分表如下： * 第二章插值法第一节拉格朗日插值多项式第三节牛顿插值多项式上一页下一页返回本章要讨论的基本问题是：研究用便于计算的简单函数近似代替一个复杂函数在某些点上的值的方法。当精确函数 y = f (x) 非常复杂或未知时，在一系列点 x0 … xn 处测得函数值 y0 = f (x0), …， y n = f (xn)，由此构造一个简单易算的近似函数 ? f (x)，满足条件 = yi (i = 0, … n)。称为f (x) 的插值函数, 点 x0 … xn称为插值节点. 称f (x)为被插值函数。将求的方法称为插值法。上一页下一页返回第一节 Lagrange插值多项式这就是拉格朗日（ Lagrange ）插值问题. 插值条件 x0 x1 x2 x3 x4 x L(x) ? f(x) 一、问题的提出上一页下一页返回定理 (唯一性) 满足的 n 阶插值多项式是存在且唯一的。证明：设所求的插值多项式为二、插值多项式的存在唯一性其系数行列式为范德蒙行列式上一页下一页返回由Crame法则，方程组有唯一解，即Ln(x)存在且唯一。三、 Lagrange 插值多项式上一页下一页返回下面用基函数的方法来构造Lagrange插值多项式. 首先考虑最简单的插值多项式：在[a,b] 上有两个插值节点x0,x1，且已知f (x)在节点上的函数值y0, y1。现在要求一个多项式L1(x)，使得：若能够找到这样的函数li(x)，即且次数不能超过1。则恰好满足插值条件。问题：怎样求出这样的li(x) , i=0,1 上一页下一页返回不妨先求l0(x) ，由于l0(x)在x1处函数值为0，显然应包括 x-x1这个因子；又因它的次数不能超过1，则l0(x)=A(x-x1) ，而l0(x)在x0处函数值为1，故 A=1/(x0-x1) ，即得：同理可得则称之为线性插值多项式。同理称之为Lagrange插值基函数. 上一页下一页返回利用拉格朗日基函数,可以构造多项式插值多项式为：线性插值多项式：n=1 y=f (x) x y x0 x1 y=L1(x) 几何意义：上一页下一页返回抛物插值多项式：n=2 插值多项式为： x y y=L2 (x) x0 y=f(x) x1 x2 几何意义：上一页下一页返回上一页下一页返回例1 解：例2 解：用线性插值. 上一页下一页返回用抛物插值计算上一页下一页返回上一页下一页返回 Lagrange插值法的流程图四、插值余项设节点在(a , b)内存在, 考察截断误差，且 f 满足条件 , 罗尔定理: 若充分光滑，，则存在使得。推广：若使得使得存在使得 Rn(x) 至少有个零点 n+1 ? = - = n i i n x x x K x R 0 ) ( ) ( ) ( 任意固定 x ? xi (i = 0, …, n), 考察 ? = - ? = n i i x t x K t Rn t 0 ) ( ) ( ) ( ) ( j ?(t)有 n+2 个不同的零点 x0 … xn x ! ) 1 ( ) ( ) ( ) 1 ( + - + n x K f n x 注意这里是对 t 求导 = + - - + + ! ) 1 )( ( ) ( ) ( ) 1 ( ) 1 ( n x K L f n n n x x ! ) 1 ( ) ( ) ( ) 1 ( + = + n f x K n x 上一页下一页返回注：? 若记，则插值余项为