5.6hermite插值.pptVIP

5.5 埃尔米特(Hermite) 插值 ★ Hermite插值描述： 设 f (x)具有一阶连续导数，已知节点上的函数值和导数值，即 (xi, f (xi)), (xi, f ’(xi)), i = 0, 1, 2, …, n, 若存在 2n+1次多项式 H2n+1(x) 满足 则称 H2n+1(x) 为 f (x) 关于节点{xi} (i = 0,1,2,…,n)的Hermite插值多项式。 记 f (xi) = yi, f (xi) = mi, i = 0,1,2,…,n . ★ 三次Hermite插值的构造 存在性 给定 f (xi) = yi, f (xi) = mi, i = 0, 1. 设 代入插值条件: H3(xi) = f(xi), H3(xi) = f (xi), i =0,1,得 基函数法 每个节点上对应两套基函数：x0: h0(x), g0(x); x1: h1(x), g1(x),满足 或简记为 先构造 h0(x), 设 h0(x) = (a + bx)(x ? x1)2 ,为方便计算，可设 由h0(x0) = 1, 得 a =1；由 所以， 同理 设 由 g0(x0) =1, 得 a =1。所以 注：我们知道，过 x0, x1 两点的Lagrange插值基函数为 显然， 于是，三次Hermite插值的基函数可表为 三次Hermite插值多项式为 容易验证，当 f(x)∈C4[a, b]时, 三次Hermite插值的误差为 提示：设 作辅助函数 固定x∈ [a, b]，则φ(t) 有三个零点 x0, x1, x, 且 x0, x1为二重零 点。反复应用Rolle定理可证。 ★ 高次Hermite插值的构造——插值基函数法 给定 n+1个节点 x0, x1,…,xn 上的函数值 f (xi)和导数值 f(xi) ，可以构造 2n+1 次Hermite插值多项式 满足 其中，hi(x), gi(x)分别为对应于函数值和导数值的不高于 2n+1 次插值基函数，它们满足 完全仿照三次Hermite插值基函数的求法，可得 容易验证， Hermite插值的误差(定理)： 当 f(x)∈C2n+2[a, b]时, 则存在ζ∈(a, b), 使 例5.8 给定 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. 例5.9 给定 f(0) = 1, f(1) = 2, f (0) = 2, 构造二次插值函数。 解 (1) 公式法 设 f ’(1) = m1，由三次Hermite插值公式得， 令 m1 = 0，得到二次Hermite插值函数 P2(x) = ?x2 + 2x + 1. 5.7 分段低次插值 n 次Lagrange插值多项式的误差： 插值多项式与被插函数的逼近程度同分点的数目和位置有关。一般地，分点越多，逼近程度越好，但也有例外。 例如 将[?1,1]10等分，步长 h = 2/10 = 0.2, 取节点 xi = ?1 + 0.2h, i = 0,1,2,…,10。以 (xi, f(xi))为插值点，构造L10(x)： Runge现象 插值多项式在插值区间上发生剧烈的震荡。它揭示了高次插值多项式存在的缺陷。 产生的原因 误差有截断误差和舍入误差两部分组成，而在插值的计算过程中，舍入误差可能会扩散或放大。 5.7.2 分段线性插值 区间[a, b]上的线性插值 I(x) 就是将每个小区间[xi , xi+1]上的线性插值 Ii (x) 连接起来, I(x) 为[xi , xi+1]上不高于一次的多项式。即 I(x)的图形是一条以 (xi, f(xi))为折点的折线。 用分段线性插值逼近上述例子的效果