第四章 插值法.ppt

设f x lnx，给定f 1 0, f 2 0.693147, f ’ 1 1, f ’ 2 0.5。用三次Hermite插值多项式H3 x 计算f 1.5 的近似值。 解 记x0 1,x1 2, 利用 2 可得 得三次Hermite插值多项式 由此得f 1.5 的近似值H3 1.5 0.409074 例3.11 练习 给定 f ? 1 0, f 1 4, f ? 1 2, f 1 0, 求H3 x , 并计算 f 0.5 . 解 x0 ? 1, x1 1, f 0.5 ≈H3 0.5 3.5625. 例3.12 给定 f 0 1, f 1 2, f 0 2, 构造二次插值函数。 解 1 公式法 设 f ’ 1 m1，由三次Hermite插值公式得， 令 m1 0，得到二次Hermite插值函数 P2 x ?x2 + 2x + 1. 解 2 扩展牛顿法--用牛顿差商表构造Hermite插值 写成差商表的形式，将带导数的节点X0及其上的函数值重复一遍，无导数的节点X1不重复，即 x f x 一阶 差商 x0 0 1 x1 0 1 2 X1 x2 1 2 1 ?1 练习 已知 计算 f 1.36 。 解 x f x 一阶 1.2 0.6 1.2 0.6 0.5 1.4 0.9 1.5 5 1.4 0.9 0.7 ?4 ?45 1.6 1.1 1.0 1.5 13 145 当x 0.9时, 用向后插值公式 2 , 这时t x-x3 /h 0.5. 将t - 0.5代入 2 , 得 f 0.9 ≈N3 0.9 2.46875. 说明：当构造向后牛顿差分的表达式时，如果采用的节点与向前差分的相同或是其中的一部分，则可以利用向前差分表求向后牛顿差分的表达式，例如 x y ?y ?2y ?3y ?4y x0 y0 x1 y1 ?y0 x2 y2 ?y1 ?2y0 x3 y3 ?y2 ?2y1 ?3y0 x4 y4 ?y3 ?2y2 ?3y1 ?4y0 Newton向后插值公式 所以可以用向前差分表中最后一行的数据计算Newton向后插值。 例 3.10 已知f x sinx的数值如表2-11的第2列，分别用二次Newton向前、向后插值公式求sin0.57891的近似值。 0.4 0.38942 0.5 0.47943 0.09001 0.6 0.56464 0.08521 0.00480 0.7 0.64422 0.07958 -0.00563 -0.00083 x sinx △ △2 △3 表2-11 0.57891 向前 向后 （只需构造向前差分表） 0.4 0.38942 0.5 0.47943 0.09001 0.6 0.56464 0.08521 0.00480 0.7 0.64422 0.07958 -0.00563 -0.00083 x sinx △ △2 △3 表2-6 解 作向前差分表如表2-6，使用Newton向前差分公式x0 0.5, x1 0.6 , x2 0.7 , x 0.57891, h 0.1,则 t x-x0 /h 0.7891, 即sin0.57891≈0.54714。误差为 向前 0.4 0.38942 0.5 0.47943 0.09001 0.6 0.56464 0.08521 0.00480 0.7 0.64422 0.07958 -0.00563 -0.00083 x sinx △ △2 △3 表2-6 向后 　　若用Newton向后插值公式（2.30），则可取x0 0.4, x1 0.5, x2 0.6, x 0.57891,h 0.1, t x-x2 /h - 0.2109。于是 即sin0.57891≈0.54707。误差为 最后，我们指出，在带有差商的差商型的插值公式中， 为了计算差商需要进行多次除法，因此当节点等距时，应 当用差分代替差商，若节点是可以随意取的，则自然应当 选为等距的。此外，利用差分作插值多项式，也要合理 地决定差分的阶数，以避免高阶差分的误差积累，而这种 积累有时是很严重的。请看下例。 例：求序列 0，0，0， ，0，0，0 的差分表 表中仅有一个微小误差 ，但在各阶差分中，随着差分的阶数的增高，误差影响增大很快，到了六阶差分时，误差积累已达 . 0 0 0 0 0 0 通常差分阶数选取的一般原则为，当差分表中 的k阶差分等于（或近似等于）0时，则插值多项式的差 分阶数只选取到k-1阶。 上面讨论差商时，节点彼此相异。 下面引入带重合节点的差商概念。若 存在。则定义 三.