第二章：插值法.ppt

在区间 上 （3） 为分段Hermite插值 §5 三次样条 定义 　　　设 。三次样条函数 , 且在每个 上为三次多项式 。若它同时还满足 则称为 f 的三次样条插值函数 注：三次样条与分段埃尔米特插值的根本区别在于S(x)自身光滑，不需要知道 f 的导数值（除了在2个端点可能需要）；而埃尔米特插值依赖于f 在所有插值点的导数值。 f(x) H(x) S(x) 2、样条函数的求解 分析： 1、区间 上的三次多项式，共需待定系数 个。 2、已知条件有 个 个 个 个 共计 个 （1）补充条件 1、 边界条件 2、 特别 称自然边界条件 3、 周期性条件（略） （2）三转角方程 假定 在区间 上的Hermite插值函数为： 则 Hermite插值 ？ 所以 P36 （2.8.7） 所以 进而 整理得 P36 同理 第一步：在区间 上做Hermite插值函数 第二步：求 第三步：求 三转角方程： 的矩阵形式 第 一 类 边 界 条 件 如果已知 特别 则在区间 上 即 整理得 同理，在可得区间 上 所以 整理得 由 得 对于自然条件 有 因为 是三次多项式，所以 是线性函数。 （3）三弯矩方程 在区间 上令 积分两次得 事实上 因此 所以 类似用同样的方法在区间 做 经过两次积分，一次微分，并求左极限得 由 并整理得 三弯矩 其中 类似用同样的方法在区间 做 经过两次积分，一次微分，并求左极限得 由 并整理得 三弯矩 其中 已知等距节点 三、Newton等距插值 1、差分 定义 简记为 简记为 简记为 向前差分 向后差分 中心差分 前差算子 后差算子 高阶向前差分 高阶向后差分 如 2、高阶差分 又如 3、前差与后差的关系 一般有 再定义 前移算子 不变算子 后移算子 则有 因此 4 差商与差分的关系 m阶向前差商与m阶向前差分的关系 查看全部 m阶向后差商与m阶向后差分的关系 又 所以 5、差分的计算 6、等距节点的Newton插值 已知等距节点 得 令 由Newton插值公式 其中 参照 即 前插公式 同理可得后插公式（P27) §3 厄尔米特插值 不仅要求函数值重合，而且要求若干阶导数也重合。 即：要求插值函数 ? (x) 满足? (xi) = f (xi), ?’ (xi) = f ’ (xi), …, ?(mi) (xi) = f (mi) (xi). 注：? N 个条件可以确定 阶多项式。 N ? 1 ?要求在1个节点 x0 处直到m0 阶导数都重合的插值多项式即为Taylor多项式 其余项为 ?一般只考虑 f 与f ’的值。 当 较大时用待定系数法求 是困难的 令 为 次多项式 且满足 其中 且满足 令 为 次多项式 所以 为Hermite插值多项式。 Kronecker（克罗内克）符号 柏林科学院院士，巴黎科学院通讯院士，伦敦皇家学会外籍会员。 主张分析学应奠基于算术，而算术的基础是整数。克罗内克名言：“上帝创造了整数，其余都是人做的工作” 令 则 其中 又 则 由（1）（2） 得 所以 其中 则 所以 令 则 又 由 得 所以 例 已知 求三次多项式 满足 解 所以 验证： 练习：（ P43（19））求四次多项式 满足 令 解： 在点0，1，2上做Lagrange插值函数 则 因此 所以 §4 分段低次插值 Remember what I have said? Increasing the degree of interpolating polynomial will NOT guarantee a good result, since high-degree polynomials are