- 1、本文档共18页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
上海交大数值分析课件数值分析2-3(牛顿插值法)
二、差商的计算 §3 差商与牛顿插值 一、差商及其性质 三、牛顿插值公式 四、牛顿插值法举例 二、差商的计算 五、牛顿插值法特点 一、差商及其性质 1. 差商的定义 二阶差商 函数关于 一阶差商 xi, xj 一般的k阶差商定义为 特别地,f(x)关于一个点xi的零阶差商定义为函数值本身,即 1.差商与节点的排列次序无关,称为差商的对称性; 性质: 2.高阶差商可由低阶差商反复作一阶差商得到,计算具有递推性 3.若f(x)在[a, b]上存在n阶导数,则 思考题:设f(x)=x3,则 f[x0, x1, x2, x3 ]= ? , f[x0, x1, x2, x3, x4 ]= ? 答:1,0 f[x0,x1,x2,x3] ∶ ∶ ∶ 三阶差商 f[x0,x1,x2] f[x1,x2,x3] ∶ ∶ ∶ f[x0,x1] f[x1,x2] f[x2,x3] ∶ ∶ ∶ f(x0) f(x1) f(x2) f(x3) ∶ ∶ ∶ x0 x1 x2 x3 ∶ ∶ ∶ 二阶差商 一阶差商 f(xi) xi 例 已知函数y= f (x)的观测数据如下,试构造差商表,并求 f [2,4,5,6]的值 1 5 9 -4 13 f(x) 0 2 4 5 6 x 解 构造差商表如下 1 四阶 -1 5 三阶 0 -5 15 2 2 -13 17 1 5 9 -4 13 0 2 4 5 6 二阶 一阶 f(xi) xi 由表可知 f[2,4,5,6] =5 二、牛顿差商插值多项式 由差商定义 可得 又 代入得 插值余项 牛顿插值 多项式 n=1 n=2 牛顿插值与拉氏插值的比较: 注:牛顿插值只需增加一项, 拉氏插值需要重新计算! 1.由插值多项式的唯一性,两余项是等价的,即 2.牛顿插值法由差商表中的差商值可判断出插值多项式的次数,而拉氏插值法则要计算到最后。 两点说明: 例 已知函数y=f(x)的观测数据如下,试用全部节点构造牛顿差商插值多项式,并用二次插值求f(3)的近似值。 1 5 9 -4 13 f(x) 0 2 4 5 6 x 四、牛顿插值法举例 解 构造差商表如下 1 四阶 -1 5 三阶 0 -5 15 2 2 -13 17 1 5 9 -4 13 0 2 4 5 6 二阶 一阶 f(xi) xi 由表可知 用二次插值求f (3)时,取 x0=2, x1=4, x2=5, 得 思考:若本题只给出前三个点,结果如何?请你总结牛顿插值法何时停止? * * *
文档评论(0)