第3章插值法.docVIP

§6 三次样条插值 分段低次差值都有一致收敛性，但对某些问题，光滑性达不到要求，譬如外观设计，像机翼形线，早期的电视机荧光屏，船体放样等问题，要求型值线过插值节点且具有二阶连续导数。这样的插值函数叫做样条函数。 样条函数还在数值微分、数值积分、微分方程和积分方程的数值解等方面有重要应用。 一、三次样条函数 定义 设是给定的节点，若函数满足： （1）； （2）在每个子区间上是三次多项式； （3）， 则称为函数的三次样条插值函数。 由定义，在每个小区间是三次多项式，需确定4个待定系数。共个子区间，故应确定个参数。 根据样条插值函数的定义，应满足： （1）在个内节点处连续，故有 （2）在个内节点处一阶导数连续，故有 （3）在个内节点处二阶导数连续，故有 （4）在个节点处满足插值条件： 以上共个条件，要确定个参数，尚需附加另外两个条件。应用中，在区间的端点处各加一个条件，称为边界条件。 边界条件通常有以下三类。 第一类：已知（指定）两端点的一阶导数值，即 该类边界条件用于样条曲线在两端有确定的倾斜角的情况。 第二类：已知（指定）两端点的二阶导数值，即 特别，时称为自然边界条件。此时样条曲线在两个端点不施加外力，样条在两端自然伸直，曲率。该类边界条件一般用于样条曲线在两端自然伸直的情况。 第三类：周期条件，即 这样确定的样条函数称为周期样条函数，用于样条曲线具有周期性的情况。 二、三弯矩方程 从力学的角度，样条曲线的二阶导数称为弯矩。 如果用待定系数法求样条函数，需解一个阶的线性方程组，当较大时，计算量相当大。而用三弯矩方程，只需解一个不超过阶的线性方程组即可。特别，对于常用的边界条件，相应的，方程组得阶数仅为。 设插值节点为，记 注意到三次样条插值函数在内的二阶导数存在，且为折线，若记在节点处的待定值为 则可用分段线性插值表示为 （1） 其中。 因在上连续，即在内节点处的二阶导数连续，故由（1）式可知 显然成立。 对（1）式两边作两次不定积分，并写成下述形式 （2） 其中，为待定系数，即（2）式的最后两项为待定线性函数。 根据插值条件，得 代入（2）式，可得 （3） 其中。 下面利用一阶导数连续条件，建立方程组。 对（3）式两边求一阶导数，可得 （4） 其中 ， 为一阶均差。 由（4）式可得 于是，由一阶导数在内节点处的连续条件 得 （5） 记 其中为二阶均差，则（5）式可写为 （6） 这样得到的（6）式，就是待定值满足的线性方程组，每个方程含有三个弯矩值，故称之为三弯矩方程。 方程祖（6）含有个未知数，但只有个方程，另外两个方程由边界条件获得，它们分别为： 第一类： 第二类： 即为已知。 第三类： 利用边界条件 得： 由于满足前面得到的个方程，将其中的第1个方程与合二为一，即得 其中 将三者分别与（6）联立，得到对应的线性方程组如下： 第一种： (7) 第二种（含个方程）： (8) 第三种： (9) 把上面的三个方程组统一写成，注意到以及，可知系数矩阵是严格对角占优阵，因此方程组的解存在且唯一，可用追赶法求解（具体方法见教材PP58-60）。把求出的弯矩值代入，即得到三次样条插值多项式。 例1：给定函数表 0 1 2 3 1 2 4 5 1 3 4 2 求满足自然边界条件得三次样条插值函数，并计算的近似值。 解：（1）该问题的边界条件是第2类，三弯矩方程组为 其中 由于，故三弯矩方程组为 （2）计算系数矩阵和常数项。 0 1 1 1 1 2 3 2 1/3 2/3 -3 2 4 4 1 2/3 1/3 -5 3 5 2 三弯矩方程组为 解得 代入 （3） 其中 得样条插值函数为 在区间上， 在区间上， 在区间上， 作业 习题三：1，2，3，7，17，21 6