样条函数，B-样条，三次样条插值，光滑余因子.docVIP

样条逼近方法 教学目的及要求： 掌握样条函数及性质、B-样条及性质、三次样条插值。 借助于多项式来逼近，虽然有很多优点，但由于多项式乃幂级数的特例，其在一点附近的性质足以决定它的整体性质。然而自然界较大范围内的许多现象，如物理或生物现象间的关系往往呈现互不关联、互相割裂的本性。亦即在不同区域中，它们的性状可以完全不相关。另一方面，从数学上讲，例如在多项式插值理论中，具有n个插值点的一元插值多项式是一个n-1次的多项式，它可能有n-3个拐点。这对于比较平滑的函数来说就不是那么理想了。 本章介绍的样条（函数）是一种分段多项式，各相邻段上的多项式之间又具有某中连接性质。因而它既保持了多项式的简单性和逼近的可行性，又在各段之间保持了相对独立的局部性质。数十年来的理论和实践表明，样条是一类特别有效的逼近工具。 §1. 样条函数及其基本性质 设给定一组结点 （1.1） 又设分段函数S(x)满足条件： 于每个区间上，S(x)是一个次数不超过n的实系数代数多项式； S(x)于上具有一直到n-1阶的连续导数。 则称为n次样条函数。常把以（1.1）为结点的n次样条函数的总体记为称为样条结点。 一个（奇次）2n-1次样条函数，如果起在区间上的表达式都是n-1次多项式（并不要求该两个n-1次多项式相同），则特别称之为2n-1次的自然样条函数。以（1.1）为结点的2n-1次自然样条函数的总体记为显然 （1.2） 下面来给出样条函数类中任一样条函数的一般表达式。 对于任意给定的以（1.1）为结点的n次样条函数根据定义，其在每个子区间上均为一n次多项式。特别地，于子区间内是一n次多项式。不妨设该多项式为。 今考虑S(x)于上的表达式。由定义，S(x)于上的表达式仍为一个n次多项式。若设该n次多项式为，并考虑下述n次多项式的性质： 按n次样条函数的定义，与于点处的值以及1阶、2阶、一直到n-1阶导数值皆相等： 亦即 是故是含这个因子。由于是一n次多项式，所以存在某常数使得 （1.3） 亦即 （1.4） 它说明S(x)于区间上的表达式恰为其前一区间上的表达式加上的某一常数倍。这样一来，S(x)于上的统一表达式应为 （1.5） 为把（1.5）写成一个统一的表达式，引入记号 （1.6） 则（1.5）所示的S(x)又可紧凑地表示为 继续采用这种分析方法，可得S(x)于整个实轴上的表达式为 （1.7） 此即为下述定理所叙述的事实。 定理 1 任一均可唯一地表现为 （1.8） 其中为实数。 显然，由（1.8）式所给出的任一函数S(x)必然满足n次样条函数的定义，亦即因而定理1可进一步写成 定理 2 为使，必须且只须存在和N个实数，使得（1.8）成立： 定理1和定理2说明函数系 （1.9） 构成n次样条函数类的一组基底。 由（1.2）和定理2可知，任一均可表示为 （1.10） 其中 当然，一函数S(x)只是满足（1.10）还不足以保证它一定是一个自然样条函数。因为它在上是否为一个n-1次的多项式尚不能保证，为保证这点，便须要求S(x)于中的表达式 亦为一n-1次多项式。即要求上述求和号这一项中n次以上的各方幂项之系数为0。但 要求上式中的系数为0，即得 （1.11） 定理 3 为使，必须且只须存在和满足线性约束（1.11）的实数，使得 （1.12） 下面讨论样条函数的积分关系式。先给出 定理 4 设S(x)由（1.8）所给出，其中n=2k-1,且 （1.13） 又设f(x)满足下述三性质： 1． 且于每个开区间内连续； 2. ; 3. ,则 （1.14） 证明 逐次采用分部积分法，有