NA002b插值.ppt

§3 牛顿插值 /* Newton’s Interpolation */ Lagrange 插值虽然易算，但若要增加一个节点时，全部基函数 li(x) 都需重新算过。 能否将 Ln(x) 改写成 的形式，希望每加一个节点时， 只附加一项上去即可？ ? ? ? ? ? 差商(亦称均差) /* divided difference */ 1阶差商 /* the 1st divided difference of f w.r.t. xi and xj */ 2阶差商 I. Newton (1642-1726) 英国数学、物理学家 微积分创立者之一 Newton插值多项式(1687) §3 Newton’s Interpolation 1 1 1 0 1 0 1 0 1 1 1 0 1 0 ] , , ... , [ ] , , ... , [ ] , , ... , [ ] , ... , , [ ] , ... , [ + + - - + + + - - = - - = k k k k k k k k k k k x x x x x f x x x f x x x x x f x x x f x x f (k+1)阶差商： Warning: my head is exploding… What is the point of this formula? 差商的值与 xi 的顺序无关！ 注：? 阶差商必须由 个节点构成， 个节点是构造不出 阶差商的。 ? 为统一起见,补充定义函数 为零阶差商。 2).差商具有对称性，即 3).若 在 上有 阶导数,且节点 则 阶差商与 阶导数有如下关系式： 4).若 是 次多项式，则其 阶差商 当 时是一个 次多项式，而当 时恒为零. ?差商性质： §3 Newton’s Interpolation 1). 阶差商可表为函数值 的线性组合，即 §3 Newton’s Interpolation 三阶差商 二阶差商 一阶差商 差商计算可列差商表如下 Newton插值是通过选取特殊的基函数来实现的，这时，取 作为Newton插值的以 为节点的基函数，而次数 不超过 的多项式 可表示为 其中 是待定系数，由插值条件 决定。 ? 牛顿插值 /* Newton’s Interpolation */ §3 Newton’s Interpolation 通过插值条件运用数学归纳法可以求得 因此就得到下列的满足插值条件(2.1)的 次插值多项式 ? §3 Newton’s Interpolation §3 Newton’s Interpolation 1 2 … … … … n?1 1 + (x ? x0) ? 2 + … … + (x ? x0)…(x ? xn?1) ? n?1 Nn(x) Rn(x) ci = f [ x0, …, xi ] ? 插值余项 /* Remainder */ 注：? 由唯一性可知 Nn(x) ? Ln(x)， 只是算法不同，故其余项也相同，即 HW p.50 #6,#8,#9 0.41075 0.57815 0.69675 0.88811 0.40 0.55 0.65 0.80 例 已知函数 在各节点处的函数值如下,用Newton插值 的值. 法求 §3 Newton’s Interpolation 解: 0.19733 0.28000 0.35893 1.11600 1.18600 1.27573 0.41075 0.57815 0.69675 0.88811 0.40 0.55 0.65 0.80 三阶 二阶 一阶 §3 Newton’s Interpolation §3 Newton’s Interpolation ? 等距节点公式 /* Formulae with Equal Spacing */ 向前差分 /* forward difference */ i i i f f f - = ? + 1 向后差分 /* backward difference */ 1 1 1 - - - ? - ? = ? i k i k i k f f f i?1 i i f f f - = ? 中心差分 /* centered difference */ 其中 当节点等距分布时: More given