数值计算方法第05章插值法介绍.ppt

* * ? 插值余项与误差估计 截断误差或插值余项 定理 若 则存在? ? (a, b), 使得 证明 故 其中 K (x)是与 x有关的待定函数. 如何求 K (x) ? * 现把 x看成是[a, b]上的固定点, 作辅助函数 即 F(t )在[a, b]上有 n+2 个零点. 根据Rolle定理, F ?(t )在 F(t )的两个零点之间至少有一个零点, 故 F ?(t )在(a, b)内至少有 (n+1)+(n+1)个零点. 对F ?(t )再应用Rolle 定理， 可知F ??(t )在(a, b)内至少有(2 n+1) 个零点. 依此类推, F(2n+2) (t )在(a, b)内至少有一个零点, 记之为??(a, b), 使得 则 * 因此 若 则 * ? 两个节点的三次Hermite插值多项式的截断误差 * 定理 满足 的 2n+1 阶Hermite插值多项式是唯一存在的. 因为H(x)为至多2n+1次多项式，故H(2n+2)(x)=0. 从而 ? Hermite插值多项式的唯一性 证明 假设 H(x)与 H(x) 是满足相同插值条件的 2n+1次Hermite多项式， H(x)也是 H(x) 的 (2n+1) 次Hermite插值多项式. 由余项公式 H(x)= H(x) * ? 分段三次Hermite插值定义 给定函数表 求分段三次Hermite插值函数H(x), 使其满足 (1) (2) 在每个小区间[xi, xi+1] ( i=0, 1,…, n－1)上，H(x)是三次多项式. * ? 分段三次Hermite插值函数H(x)的分段表达式 * ? 分段三次Hermite插值函数的误差估计 其中 ? 分段三次Hermite插值函数的导数在整个区间 [a, b]上是连续，但二阶导数在内节点处不连续. ? 为了使二阶导数也在内节点处连续，可用三次样条插值函数. * ? Hermite插值的一般形式 求一至多 n+m+1 次的多项式 H(x), 使得 已知函数 f (x)在(n+1)个互异节点 处的函数值； 以及某些节点上的导数值 * ? 插值多项式H(x)存在唯一 ? 插值余项 ? 如何求？ ? 待定系数法 ? 求基函数方法 ? 利用Newton插值或Lagrange插值方法. * 例 按下表求 f (x)的三次Hermite插值多项式H(x), 并写出截断误差R (x)=f (x)－H(x)的表达式. 0 1 2 3 2 －1 1 解 设 则由已知条件得 故 * 例 求满足 H( xj )=f ( xj ) ( j=0,1,2 )及 H ?(x1)=f ?(x1)的插值多项式及其余项表达式. 其中A为待定常数， 解 设 由 * 余项表达式 * 解2 用插值基函数方法 * 数学上三次样条插值中的“样条”来源于工程师的做法。放样：要求曲线光滑地通过型值点。 工程师的做法：使用一种有弹性的细长木条(或金属条)，称之为样条(Spline)，强迫它弯曲通过已知点。 §5 样条插值 样条: 是 指飞机或轮船等的制造过程中为描绘 出光滑的外形曲线(放样)所用的工具. 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的 1946年,Schoenberg将样条引入数学,即所谓的样条函数 * ? 三次样条插值 定义已知函数 y=f (x)在区间[a, b]上的 (n+1)个节点 上的函数值 yi =f (xi ) (i=0,1,…,n), 求插值函数 S(x) 在每一个小区间 上是三次多项式; (1) (2) (3) 函数S(x)称为f (x)的三次样条插值函数. * x0 x1 x2 xn S0 (x) S1 (x) Sn?1 (x) Si(x) = ai x3 + bi x2 + cixi + di 4n 个未知数 连续斜率与曲率 xn?1 三次样条插值 * ? 从定义知要求出S(x), 在每个小区间[xj , xj+1]上要确定4个待定系数, 故应确定 个参数. ? 根据S(x)在区间[a, b]上二阶连续可导, 在(n?1）个内节点xj (j=1, 2, …, n?1)处应满足连续性条件 共有3n?3个条件, 再加上S(x)满足插值条件 yi =f (xi ) (i=0,1,…,n), 共有 个条件, 因此还需要 个条件才能确定S(x). 通常在区间[a, b]端点a=x0, b=xn 上各加一个条件(边界条件). 4n 4n?2 2 * ? 三类