样条(Spline)函数.docVIP

样条（Spline）函数 I、三次Spline插值函数的定义 给定区间上的一个分划，且 … 和一组函数值，，，…，，如果具有下列性质： ； 在每个子区间上，是一个三次多项式； ，=0，1，…，。（称为插值条件） 则称是关于分划的分段三次样条函数，简称为Spline函数。 将表示成如下分段形式： 其中是一个三次多项式，且满足插值条件： ， 为了得到的具体表达式，根据可知，对每个内部结点有： 下面利用上述条件导出的表达式。 II、三次Spline插值函数的表示 设，则。利用分段Hermite插值公式得 （*） 其中和称为Hermite插值的基函数（在下面附录中将详细介绍），这里的一阶导数是未知待定的。由上式不难看出，只要确定了各插值点的一阶导数，Spline插值函数就相应确定了。对（*）两边同时求二阶导得： 同理： 利用上面两式再结合二阶导连续的条件可得： 利用和的性质（在附录中介绍）上式等价于： （**） 其中，若令，，则（**）可划为： ， (***) 其中， 这里的(***)构成了关于个未知数的方程组，但是此处只有个方程，要确定这些未知数的值，还需要再加两个条件，一般在的端点上各附加一个条件，称为边界条件，不同的边界条件对应着不同的边值问题，下面是常见的三种边界条件： 第一类边界条件： ， 即是，。记，，，。 第二类边界条件： ， 如果，称对应的Spline函数为自然Spline函数。 即有： 化简之后为： 记，，，。 第三类边界条件： ， 这类函数一般要求，它对应的Spline函数称为周期Spline函数。 针对第一类边界条件和第二类边界条件，他们构成的方程组有如下统一形式： 因为该系数矩阵为三对角阵，且严格对角占优，因此该方程组有唯一解，由此可见，第一、二类型的边值问题的Spline插值函数是存在唯一的，利用追赶法可求得，这样就得到了，从而确定了Spline插值函数。 附录： 前面学过分段三次Hermite插值多项式的表达式为： 其中 其中。 分段三次Hermite插值多项式也可表示成以结点的函数值和导数值为系数、关于基函数和的线性组合： 其中 基函数和有如下性质： ①、 ②、 ③、 ①、且 ②、 ③、 当时，， 这些性质的证明只需要对相应的函数求导即可得出，我就不那个了……！ 引入上述基函数以后，此时有： ， 这正是上面定义样条函数是用到的式子。