计算方法第2章插值法与逼近.ppt

第2章 插值法与逼近 关于节点 的二重Hermite插值多项式为 2.5 分段低次插值 龙格现象 例：在[?5, 5]上考察 的Ln(x)。取 - 5 - 4 - 3 - 2 - 1 0 1 2 3 4 5 - 0.5 0 0.5 1 1.5 2 2.5 n 越大， 端点附近抖动 越大，称为 Runge 现象 Ln(x) ? f (x) ? 分段低次插值 例：在[?5, 5]上考察 的Ln(x)。取 在实际工作中，常用的是一、二、三次插值。 2.5.1 分段线性插值 2.5.2 分段二次插值 2.5.3 分段三次Hermite插值多项式 分段三次Hermite插值的余项公式为 2.6 三次样条插值 2.6.1 样条函数的概念 2.6.2 三次样条插值函数的定义 注：三次样条与分段 Hermite 插值的根本区别在于S(x)自身光滑，不需要知道 f 的导数值（除了在2个端点可能需要）；而Hermite插值依赖于f 在所有插值点的导数值。 2.6.3 边界条件和三次样条插值多项式 在 上，记 ] , [ ) ( ) ( 1 ] [ j j j x x x x S x S - ? = 对每个j, 此为3次多项式 则 S[j]”(x) 为 次多项式，需 个点的值确定之。 1 2 设 S[j]”(xj?1) = Mj?1， S[j]”(xj) = Mj 对于x ?[xj?1, xj ] 可得到 S[j]”(x) = j j j j j j h x x M h x x M 1 1 - - - + - 积分2次，可得 S[j]’(x) 和 S[j](x) ： j j j j j j j A h x x M h x x M + - + - - - - - 2 ) ( 2 ) ( 2 1 1 2 1 S[j]’(x) = j j j j j j j j B x A h x x M h x x M + + - + - - - 6 ) ( 6 ) ( 3 1 3 1 S[j](x) = 利用已知S[j](xj?1) = yj?1 S[j](xj) = yj 可解 j j j j j j j h M M h y y A 6 1 1 - - - - - = j j j j j j j j j j j j h x x h M y h x x h M y B x A 1 2 2 1 1 ) 6 ( ) 6 ( - - - - - + - - = + 下面解决 Mj ： 利用S’ 在 xj 的连续性 [xj?1, xj ]: S[j]’(x) = j j j j j j j j j j j h M M x x f h x x M h x x M 6 ] , [ 2 ) ( 2 ) ( 1 1 2 1 2 1 - - - - - - + - + - - 1 1 1 1 2 1 1 2 1 6 ] , [ 2 ) ( 2 ) ( + + + + + + + - - + - + - - j j j j j j j j j j j h M M x x f h x x M h x x M [xj , xj+1]: S[j+1]’(x) = 利用S[j]’(xj) = S[j+1]’(xj)，合并关于Mj?1、 Mj、 Mj+1的同类项，并记 , , , 整理后得到： 1 1 j j j j h h h + + + = l 1 j j - = l m ]) , [ ] , [ ( 6 1 1 1 j j j j j j j x x f x x f h h g - + + - + = 2 1 1 g M M M j j j j j j = + + + - l m ? j ? 1 n?1 即：有 个未知数， 个方程。 n?1 n+1 还需2个边界条件 上一页 下一页 返回首页 2.1 插值与逼近的基本概念 2.2 拉格朗日插值多项式 2.3 牛顿插值多项式 2.4 Hermite插值 2.5 分段低次插值 2.6 三次样条插值 2.7 最小二乘法 2.8 数值微分 2.1 插值与逼近的基本概念 2.1.1 插值与逼近的概念 逼近是指对函数类A中给定的函数 要求在另一类较简单便于计算的函数类 中求函数 属于