《计算方法》第四章插值方法概要.ppt

例题: 已知函数 f(x) 有如下数据: 插值的几何意义 从几何上看，插值就是求一条曲线 使其通过给定的 个点 ， 并且与已知曲线 有一定的近似度。 线性插值与其基函数示意图 （Lagrange）插值多项式 由于 是不能确定，因此我们并不能确定误差的大小 但如能求出 ，那么用 逼近 的截断误差限是： 当 时， 当 时 其截断误差为： 其中, 因为 可取 于是： 用抛物线插值时，取所有节点，得到 §4.3.1 差商的概念 从零阶差商出发，归纳地定义各阶差商: n 阶差商的概念 差商的基本性质 性质1：差商可表示为函数值的线性组合，即： 性质2：差商关于所含节点是对称的，即： 差商的基本性质 性质3： 性质4：设 在 存在 n 阶导数，且 则 ，使得: 差商的计算-差商表 §4.3.2 差分 在前面的讨论中，节点是任意分布的，但实际上经常遇 到等距节点的情况，这时插值公式可以得到简化，为此，我 们先介绍差分的概念。 设函数 在等距节点 上的值 为已知，这里 为常数，称为步长。 差分的定义 记号 分别称为 在 处以 为步长的 向前差分、向后差分、中心差分 符号 、 、 分别称为向前差分算子、向后差分算子、 中心差分算子. 高阶差分 用一阶差分可以定义二阶差分 一般地可定义 m 阶差分为： 中心差分定义为: 以此类推。 不变算子 I、移位算子 E 定义 从而可得： 于是得到： 同理，由于： 得到： 由于： 得到： 由差分的定义及不变算子和移位算子有如下性质: 差分的性质 性质1：各阶差分均可用函数值表示，如： 性质2：某点的函数可用各阶差分来表示： 性质3：差商与差分有如下关系： 性质4：差分与导数有如下关系： 差分的计算 4.4 牛顿插值公式 根据差商的定义，把 看成 上的一点， 可得： 4.4 牛顿插值公式 根据差商的定义，把 看成 上的一点， 可得： 其中 显然 满足插值条件，且次数不超过 ,它就 是插值多项式，其系数为： 我们称 为牛顿插值多项式. 已知 的函数表，求4 次牛顿插值多项式, 并求 已知 的函数表，求4 次牛顿插值多项式, 并求 已知 的函数表，求4 次牛顿插值多项式, 并求 牛顿插值公式和Lagrange插值公式比较 和 均是 n 次多项式,且均满足插值条件: 由多项式的唯一性, ,因而,两个公式 的余项是相等的,即 当插值多项式从 n-1 次增加到 n 次时, 拉格朗日型插值必须重新计算所有的基本插值多项式; 而对于牛顿型插值,只需用表格再计算一个 n 阶差商, 然后加上一项即可。 4.5 分段插值公式 在区间 [a, b] 上用插值多项式 P 逼近函数 f 时，f 和P 在每个节点上的差异(理论上)应该为零。自然，我们期望 在一切中间点上也能很好地逼近 f ，并且当插值点增加时 这种逼近效果应该越来越好。 但上述的期望不可能实现的。当认识到这一点时，在数 学界曾引起强烈的震动。 4.5.1 分段线性插值 分段线性插值特别简单，从几何上看，就是用折线逼近 曲线。 分段线性插值的数学定义 设 是区间 上的函数，在节点 上的函数值为 ， 求一分段折线函数 满足： （1） （2） 在 上， 是一次多项式。 （3） 则称 为 的分段线性插值函数。 4.5.1 分段线性插值的基函数 当 时, 显然 是 的线性组合： 在区间 上的值为： 4.5.2 埃尔米特 (Hermite