第2章+插值法4.ppt

* §6 样条插值 问题的提出: 上面讨论的分段低次插值函数都有一致收敛性， 但光滑性较差，对于像高速飞机的机翼形线，船 体放样等型值线往往要求有二阶光滑度，即有二 阶连续导数，早期工程师制图时，把富有弹性的 细长木条（所谓样条）用压铁固定在样点上，在 其它地方让它自由弯曲，然后画下长条的曲线， 称为样条曲线。它实际上是由分段三次曲线并接 而成，在连接点即样点上要求二阶导数连续，从 数学上加以概括就得到数学样条这一概念。下面 我们讨论最常用的三次样条函数。 * 三次样条函数: 定义:函数 ，且在每个小区间 上是三次多项式，其中 是给定节点，则称 是 节点 上的三次样条函数。 若在节点 上给定函数值 ，并成立 则称 为三次样条插值函数。 * 从定义知要求出 在每个小区间 上是三次多项式，要确定4个待定系数，共有n个小区间，故应确定4n个参数。 根据 在 上二阶导数连续，在节点 处应满足连续性条件 共有3n-3个条件，再加上 满足插值条件 ，共有4n-2个条件，因此 还需要2个条件才能确定 。 * 通常可在区间 端点 上各加一 个条件（称为边界条件），可根据实际问题的 要求给定。常见的有以下三种： 1° 已知两端的一阶导数值，即 ，其特殊情况 ，样条曲线在端点 处呈水平状态。 2° 两端的二阶导数已知，即 其特殊情况 ，称为自然 边界条件。 * 3°当 是以 为周期的周期函数时， 则要求 也是周期函数。这时边界条件应 满足 而此时 。这样确定的样条函 数 ，称为周期样条函数。 * 三转角方程: 现在构造满足条件 及加上相应边界条件的三次样条函数 的表 达式。 若假定 在节点 处的值为 ，则由分段三次埃尔米特插值 公式可得 其中 是插值基函数。 * 显然，表达式中 及 在整个区间 上连续，且满足 ；现需 确定 ，可利用 及某一边界条件来确定。为了求出 ，我们考 虑 在 上的表达式 * 这里 。对 求二次导数得 于是 同理，可得 在区间 上的表达式 * 及 由条件 ， 可得 用 除全式，并注意 ，上面方程可简化为 * 其中 此方程是关于未知数 的n-1个方 程，若加上边界条件： ，则方程 变为只含