三次养条插值计算.doc

所谓样条函数，从数学上而言，就是按一定光滑性要求“装配”起来的分段多项式。具体而言，称具有分划 的分段k次式Sk (x)为k次样条函数.如果它在每个内节点 上具有直到(k-1)阶连续导数，则称内节点 为样条函数Sk (x)的节点。 样条函数常简称为样条，其特点是，它既具有充分的光滑性（Sk (x)具有直到(k-1)阶连续导数），又保留有一定的间断性（Sk (x)的k阶导数在节点处可能有间断）。光滑性保证了外形曲线的美观；间断性则使它能转折自如地被灵活运用。 样条（spline）这个词本来是在飞机和船舶的制造过程中为了描绘出光滑的外形曲线所用的一种绘图工具。它是一种富有弹性的细长条，使用时，用压铁固定在一些给定的点（称为样点）上，而在其他地方任它自由弯曲，然后依样画下光滑的曲线，就称为样条曲线。从数学上看，在小挠度的假定下，它实际上是由一段一段的三次多项式曲线拼接而成的，在拼接处，函数自身不仅是连续的，而且它的一阶和二阶导数也是连续的，但是三阶导数一般是不连续的。因此，样条曲线具有非常好的光滑性。从数学上加以概括就得到了三次样条插值函数这一概念。 近几十年来，函数逼近在理论研究和实际应用中均获得了重大的发展，其中非常流行和十分活跃的分支就是样条函数。它已经以各种方式应用到数据拟合、数值微分、数值积分、微分方程和积分方程的数值求解中，而且三次样条插值函数的计算也比较简单。 2.三次样条插值 样条插值实际上就是一种改进的分段插值。 定义 对于函数f (x)在(n+1)个互异点 处给定函数值 其中， 若函数S(x)满足条件： (1) S(x)在每一个子区间 上是一个三次多项式； (2) S(x)在每一个内节点 上具有直到二阶的连续导数即 则称S(x)为节点 上的三次样条函数。若S(x)在所有节点上满足插值条件： 则称S(x)为三次样条插值函数。 由以上定义可知，三次样条插值函数就是全部通过样点的二阶连续可微的分段三次多项式函数。 由定义的条件（1）知，S(x)在每一个子区间[xj-1,xj]上是一个三次多项式，若记为Sj(x)则可设 从而有 这是一个分段函数，其中Sj(x)由 唯一确定，因此要确定整个区间[a,b]上的三次样条插值函数S(x)，必须确定4n个未知系数 下面分析这4n个未知系数的确定方法。 在[xi, xi+1]上构造一个三次多项式 共有 4n 个未知量，因此需要 4n 个条件。 由定义中的 (3)知，有n+1个条件； 由定义中的 (2)知，有3(n ?1)个条件，即 上面共有4n-2 个条件，因此还需要2个两个条件才能确定4n个系数，这就要用到边界条件。 3.三次样条插值函数的求法 为了构造三次样条函数，需要补充边界条件，一共有三种边界条件，它们分别是： (1) 第一种边界条件 （1） (2) 第二种边界条件 （2） 若 ，则称之为自然边界，满足自然边界条件的样条函数称为自然样条。它是通过所有数据点的插值函数中，总曲率最小的唯一函数，因此自然三次样条是插值所有数据点的最光滑的函数。 (3)第三种边界条件 当被插值函数是以b-a为周期的周期函数时，则要求S(x)也是周期函数，这时边界条件为 （3） 设对区间[a,b]的分划 是任意的。记 ，由于S(x)在每个子区间[xj-1,xj]上是一个三次多项式Sj(x)，所以S’’j(x)在区间[xj-1,xj]上是关于x的线性函数，且 因此S’’j(x)的表达式可表示为： （4） 对上式积分两次得：