Ch2插值法..pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * §2.6 分段低次插值法 一、高次插值的病态性质 龙格(Runge)现象…. 二、分段线性插值 所谓分段线性插值就是用通过插值点的折线段逼近f(x). 二、分段三次埃尔米特插值 分段线性插值函数导数间断，若已知节点上函数值和导数，可构造一个导数连续的插值函数Ih(x)，满足 作业 P60, 15,16. §2.7 样条插值 问题背景… 一、样条插值的概念 二、三次样条插值函数的建立 三、误差界与收敛性 作业 P60, 20(1)只列出方程组; 20(2)应得结果. * * * * * * * * 差商的基本性质: 由上式得差商表: k xk f(xk) 一阶差商 二阶差商 三阶差商 … 0 1 2 3 4 ┆ x0 x1 x2 x3 x4 ┆ f(x0) f(x1) f(x2) f(x3) f(x4) ┆ f[x0, x1] f[x1, x2] f[x0,x1,x2] f[x2, x3] f[x1,x2,x3] f[x0,x1,x2,x3] f[x3, x4] f[x2, x3,x4] f[x1,x2,x3,x4] … ┆ ┆ ┆ 规定函数值为零阶差商 Newton插值多项式 差商表: k xk f(xk) 一阶差商 二阶差商 三阶差商 … 0 1 2 3 4 ┆ x0 x1 x2 x3 x4 ┆ f(x0) f(x1) f(x2) f(x3) f(x4) ┆ f[x0, x1] f[x1, x2] f[x0,x1,x2] f[x2, x3] f[x1,x2,x3] f[x0,x1,x2,x3] f[x3, x4] f[x2, x3,x4] f[x1,x2,x3,x4] … ┆ ┆ ┆ P26-27 例2.4 * 解：先造差商表 例1 k xk f(xk) 一阶差商 二阶差商 三阶差商 四阶差商 0 1 2 3 4 0.40 0.55 0.65 0.80 0.90 0.41075 0.57815 0.69675 0.88811 1.02652 1.1160 1.1860 0.2800 1.2757 0.3588 0.197 1.3841 0.4336 0.241 0.034 * 由Newton公式得四次插值多项式为： k xk f(xk) 一阶差商 二阶差商 三阶差商 四阶差商 0 1 2 3 4 1 2 3 4 5 1 4 7 8 6 3 3 0 1 -1 -1/3 -2 -3/2 -1/6 1/24 §2.4 差分与等距节点的Newton插值 上节讨论任意分布节点的插值公式，应用时常碰到等距节点的情形，此时插值公式可简化，为此先介绍差分． 一、差分及其性质 差分的基本性质: * 可以证明 如 * 差分表 差分表: ?2f0 ?2f1 … ┆ ?2f2 ┆ ┆ ?f0 ?f1 ?f2 ?f3 ┆ f0 f1 f2 f3 f4 ┆ 0 1 2 3 4 ┆ ?2 ?3 … ? fk k 二、Newton向前、向后插值公式 称为牛顿后插公式。 * Newton插值法的优点是计算较简单,尤其是增加节点时, 计算只要增加一项,这是Lagrange插值无法比的. 但是Newton插值仍然没有改变Lagrange插值的插值曲线 在节点处有尖点，不光滑，插值多项式在节点处不可导 等缺点 P31 例2.5 §2.5 埃尔米特（Hermite)插值