第四章 多项式插值与函数逼近1.pptVIP

写成方程组的形式： n-1个方程 n+1个未知数 三转角方程 上述方程组称为 的m连续方程 M、m连续方程的求解：需要补充附加条件 3、边界条件 ?已知端点的斜率： ?已知端点的二阶导数： ?设 是以 为周期的周期函数，对 附加周期性条件： 即要求三次样条插值函数在端点处函数值、一阶导数值和二阶导数值相同。 M连续方程在各类边界条件下的求解方法 ?对于第一类边界条件 由 得 从而得到方程组（三对角）： 可用追赶法求解 ?对于第二类边界条件 类似地可以得到方程组（三对角）： 上述两种情况得到的方程组，可以写成统一形式： 其中 时为第二类边界条件 时为第一类边界条件 ?对于第三类边界条件 得到方程 其中 第三类边界条件对应的方程组： 对三对角算法经过修改后可以求解上述方程组 不是三对角方程组 注：三次样条与分段 Hermite 插值的根本区别在于S x 自身光滑，不需要知道 f 的导数值（除了在2个端点处的函数值）；而Hermite插值依赖于f 在许多插值节点的导数值。 f x H x S x 0 1 2 3 0 2 3 6 1 0 例1：已知函数 在的 数据表： 解： 求 在区间 上的三次样条插值函数。 三弯矩方程组 用追赶法求解方程组得 三次样条插值函数 4、三次样条函数的性质 性质1（极小模性质） 设 是任一被插函数， 是自然三次样条插值函数，则成立 式中等号当且仅当 时成立。 证明： 因为 故结论成立 （积分可加性质） 下面证明 分步积分 分步积分 等号成立 为线性函数 （由插值条件） 性质2（最佳逼近性质） 设 是任一被插函数， 是带有斜率边界条件的三次插值样条函数， 是与 有相同分割的任一三次样条函数，则有 证明： 证明方法与性质1类似 （自己证明） 性质3（误差估计） 设函数 ， 是区间 的一个分割， 是关于 的带有Ⅰ型 斜率边界 或Ⅱ型 二阶导数边界 边界条件的插值函数，则有误差估计 其中 是分割比，并且系数 与 是最优估计。 性质说明：三次样条插值函数本身连同它的一、二、三 阶导数分别收敛到 及其相应导数，具有强收敛性。 * 全导数的Hermite插值多项式 其中 x0 x1 x2 x3 x4 x H9 x ? f x 全导数的Hermite插值多项式的几何意义 如n 1时Hermite插值多项式 为 ? Quiz: 给定 xi i +1, i 0, 1, 2. 下面哪个与 的图像形状相似？ 斜率 1 x 0 - - 1 0.5 1 2 3 4 5 6 y x y 0 - - - 1 0.5 1 2 3 4 5 6 二、 全导数的Hermite插值多项式的存在唯一性和余项 满足前述插值条件的不超过2n+1次的插值多项式 是唯一存在的。 证明： 存在性由前面的构造过程已证； 也是满足前述插值条件的不超过2n+1次的插值多项式 设 令 因此，多项式 有n+1个二重零点 矛盾！ 设被插值函数 在区间[a,b]上有2n+1阶连续导 数， 在区间（a,b）内存在, 是 关于互异节点 的满足前述插值条 件的不超过2n+1次的插值多项式，则对 成立 证明方法同以前多项式插值余项完全类似 0 1 2 1 2 3 -1 0 1 应用Hermite插值计算 的近似值。 例1：已知函数 在点 数据表： 解： 三、 不完全导数的Hermite插值多项式举例 例2：已知函数 在 的数据表： 构造Hermite插值多项式 满足 并推导其插值余项（已知 具有4阶连续导数）。 解： 首先构造满足插值条件 的多项式 令 由 插值余项的推导： 设 构造辅助函数 有零点 且 由罗尔 Roll 定理： 至少有4个零点 反复利用Roll定理 至少有1个零点 即 ？ §4 分段插值 一、高次插值评述 从插值余项角度分析 为了提高插值精度，一般来说应该增加插值节点的个数，这 从插值余项的表达式也可以看出，但不能简单地这样认为，主要 原因有： ?余项公式成立的前提条件是 有足够阶连续导数（即函数足够光滑），但随着节点个数的增加，这个条件一般很难成立；； ?随着节点个数的增加， 可能会增大。 随着节点个数增加到某个值，误差反而会增加。 注意下面图中 曲线的变化情况！ 例3：在[?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 ? 二、 分段插值的构造方法 将插值区间划分为若干个小区间（通常取等距划分） 采用低次插值 在区间 上