4.3 样条函数 插值

《4.3 样条函数 插值》由会员分享，可在线阅读，更多相关《4.3 样条函数 插值（42页珍藏版）》请在金锄头文库上搜索。

1、样条函数插值SPLINEINTERPOLATION信息与数学学院 内容提要引言样条函数的物理背景一般K次样条3次样条插值高次自然样条与B 样条基础 4 8样条函数插值4 8 1引言样条函数的物理背景回顾前面几节讲过的各种代数插值 它们有一个共同的弱点 那就是 它们都是相当刚性 stiff 的 也就是说 局部数据误差易向远处传播 放大 以Lagrange插值为例 设数据真值被代之以含有误差的 令是以为插值条件的插值多项式 于是最终的插值误差是由讲义第168页插值公式 6 知 上式右端第二项为这表明结点处的数据误差通过插值基函数放大和扩散 更何况 如果被插函数有奇点 甚至只要解析延拓到复平面有隐秘

2、奇点出现 则当为高次多项式时 误差放大和扩散还将助长很可怕的强振荡 展示Runge现象的著名例子就清楚地描述了这种振荡 右图 相同数据3次样条插值与Lagrangr插值效果比较CubicSplineInterpolationLagrangrInterpolation 如果采用分段多项式插值 则由于插值基函数只是局部活跃 它们的支集是局部紧致的 结点上的误差可以被控制在小的范围内 因而也带来了内在的高度稳定性 这是分段插值的一大优势 许多实际问题希望插值函数具有较高阶的整体光滑性 此时 高次Hermite插值或分段高次Hermite插值可以利用 注意 分段高次Lagrange插值和Newton插

3、值等是做不到的 在插值结点上它们只能保证插值函数连续 注 函数的支集Supp定义为Supp 但高次Hermite插值在许多场合中看不中用 提高Hermite插值多项式的次数就要增加约束条件 给出插值结点处被插函数及其直到足够高阶导数之值 作为约束条件的所有数据都是通过观测得到的 而观测总难免有误差 于是高次插值不仅增添了数据准备和计算的困难 也将导致更大的误差 还有许多应用不仅要求插值函数具有足够高阶的整体光滑性 还要求在某些结点处转折灵活 例如若干点处加载集中力的杆 梁或板弯曲 这就导致本节要讨论的样条函数 Spline 插值 数学里的样条 Spline 一词来源于它的直观几何背景 绘图员或

4、板金工人常用弹性木条或金属条加压铁 构成样条 来绘制或者放样成光顺曲线或者曲面 但它之所以成为数值分析的标志性成果之一并且在数学物理的广泛领域获得非常成功的应用 还在于它的明确的物理背景 请看下面的例子 例1 如图4 8 1 一均匀弹性弦两端固定于两点 在区间内取点列图4 8 1并在内结点集上分别给集中载荷则载荷分布可表为其中是集中于结点的点脉冲函数 事实上 在小变形和均匀分布外力假设下 上述弦的平衡问题的微分方程模型乃是两点边值问题现在是作用离散的集中力 此时弦达到平衡状态时位移函数应满足 因此有这意味着 在每一加载点处脉冲间断 是阶梯函数 是分段线性的连续函数 在每一内结点转折灵活 后面我

5、们将指出 如此的便是一次样条函数 例2 考察梁弯曲方程与例1类似地加载集中力 只是两端点除给零位移约束外还要加一阶或二阶导数约束 于是集中力作用下的梁弯曲方程成为此时我们得到 图4 8 2 在内结点上脉冲间断 为阶梯函数 在每个子区间上是三次多项式 和都是上的连续函数 这便是后面我们要着重讨论的三次样条 此例展示了三次样条的如下特征 它分段三次光滑 整体二次光滑 足够光滑 在内结点处三阶导数间断 转折灵活 上面两个例子分别涉及弹性弦和梁的小变形平衡 这就自然会想到弹性力学中联系平衡态与变形能的两个重要的极值原理 最小势能原理和虚功原理 既然一次和三次样条分别与弦和梁的小变形平衡问题联系着 那么

6、直觉告诉我们它们也应有相应的极值性质 后面我们将证明的确如此 作为伏笔 我们指出 上两例中弹性弦和梁的变形能分别表为和 4 8 2一般k次样条定义4 8 1 次样条函数 设是区间上的一个分划或分割 即称为定义在区间上关于分划的一个次样条函数 如果 在每一区间上是次数不超过的多项式 2 在区间上是次连续可微的 节点处的阶导数间断 因而转折灵活 次样条函数类记为为方便后面的讨论 我们将样条函数写成如下形式 例 根据上述定义 0次样条函数为分段常数 即阶梯函数 它可表为1次样条函数为分段线性函数 它可表为一般二次多项式不是严格意义下的二次样条 4 8 33次样条插值问题的提法 给定数据表构造3次样条

7、函数满足插值条件 构造方法 应具有如下形式并且满足条件 4 2 和 因是分段3次多项式 故在每个区间上都是3次多项式 从而共须个独立条件确定 和在个内结点连续 即满足条件 4 4 因而 4 4 给出了个条件 4 2 提供了个独立条件 还差2个条件 有多种给法 最常见的给法是 i 简支边界 导致三弯矩关系式 关系式 特别地 自然边界 三次自然样条 ii 固支边界 导致三转角关系式 关系式 注意 上述 给出的个条件是问题本身隐含的 和 共个独立条件须提供 故结点三次样插值问题只有个自由度 请与分段三次Hermite插值比较 三次自然样条插值关系式的构造方法记注意到于是连续的和分段线性的 从而在每个

8、上是线性的 故可表为对此式积分两次并应用条件 4 2 可得到 微分 4 6 可得到和由 4 4 第二式 4 7 与 4 8 应相等 于是得到 利用条件得到关于的线性方程组其中解出代入 4 6 即得到 CubicSplineInterpolationLagrangrInterpolation 类似地可构造关系式 留作练习 下面的两个定理有重要意义 定理 4 1 设则定理 4 2 设则 我们只证明第二个定理 证明之前 先解释它的重要意义 前面曾提到 梁弯曲达到平衡态时变形能 内能 可表为于是 4 11 意味着 在满足端点约束和以及插值条件 的一切连续二次可微的函数中 三次自然样条插值函数使得变形能

9、达到最小的 当时 注意曲线的曲率 而积分达到最小意味着三次自然样条插值函数是各种可能的插值函数中使得均方曲率为最小的插值函数 即在一定意义下最为光顺的插值函数 这是三次自然样条插值函数的一个非常特别和有用的性质 下面证明定理 4 2 证明令 则并且下证 利用分部积分并注意到以及得到 4 8 4高次自然样条与B 样条高次自然样条确切地说 高次自然样条只对奇数阶定义 下面定义次自然样条 定义4 8 2 自然次样条函数 给定区间上的分划如定义4 8 1 函数称为自然样条 如果它 i 在每个子区间上都是次数不超过的多项式 在区间和上是次数不超过的多项式 分划下的次自然样条函数构成的集合记为 依此定义

10、三次自然样条在区间和上退化为一次多项式 为方便后面导出B 样条 我们考虑利用构造Lagrange插值或Newton插值的结点基函数那样的方法 使一般次自然样条能用适当的基函数作 元件 来 组合装配 出来 对次幂函数作一种 截肢 手术 截去的部分并以0接上 得到所谓半截幂函数 truncatedpowerfunction 这种手术导致处阶导数间断 有跃度 但直到阶导数仍连续 我们来观察跃度 阶梯函数的基函数此函数本身在处有间断 跃度为1 一次单项式的半截幂及其1阶导数 次单项式的半截幂及其阶导数 至此我们已经可以观察到 单项式半截幂可被视为一个结点的次样条 平移后得到的也是一个结点的次样条 平移

11、 压缩 组合后得到的是两个结点的次样条 尤为重要的是 利用半截幂这样简单的样条可构造出一系列新的样条 定理 4 3 前述分划上存在唯一自然样条且可表为其中证明 留作讨论题 思考题 半截幂函数的支集Supp为何 定性地分类函数数据误差的扩散情况 B 样条基础前已指出 我们可以利用半截幂这样简单的样条作为基装配成一系列新的样条 下面讨论通过提供基函数 BasesFunction 来构造样条函数空间的一般方法 首先将有限的结点集扩展为无穷集分划其次定义半截幂的一阶差分如此得到的样条称为0次B 样条 B spline 有如下性质 对所有的和所有的 对所有的于是0次样条函数可表为阶梯函数递归地定义次B

12、样条如下 记则 4 12 可写成例1 1次B 样条 图4 8 3 图4 8 3于是1次样条函数可表为分段线性函数 例2 2次B 样条 图4 8 4 图4 8 4例3 3次B 样条 图4 8 5 图4 8 5 次B 样条有如下性质 证明作为课外讨论题 B 样条组是上的线性无关组 在相同的分划之下 越大 的集Supp越大 振幅越小 因而误差随远离结点而衰减 记则有 若为一致分划 则B 样条提供了一类具有若干特殊性质的基 B 样条插值我们的兴趣在于利用B 样条作插值 为此 我们限制于区间 插值问题的提法为 求形如满足条件 和 某种边界约束 有关B 样条的更多理论与应用讨论留待下次课 小结 样条相对于

13、其它分段多项式插值的主要优点是它保证了较高阶的整体光滑性 实践中 曲率间断就须仔细打量才能察觉 三次样条能保证二阶导数连续 三阶导数仅在结点间断 如此的光滑通常是足够的 三阶导数在结点间断还带来了转折灵活的特点 B 样条的支集虽然随阶数的提高而扩大 但结点数据误差仍可控制于局部 且愈远离结点愈衰减 这使得样条插值的稳定性要优于其它高次多项式插值 次B 样条之性质 意味着越大越光滑 但它的支集也就越宽 误差将传播到愈远处才能消失为0 这就是高次样条也会有多余的波动的原因 样条插值计算是隐式的 需要形成并求解一个线性方程组 不如前面讲过的多项式插值计算简便 有关样条插值收敛性和误差估计方面的问题我们没有提及 建议阅读教材 样条函数还有其它一些有用性质 例如 次自然样条也有类似于 4 11 的最优性不等式这个不等式深刻揭示了自然样条的和谐之美 老师鼓励并且非常乐意帮助对样条函数方面的知识有兴趣的学生更深入地学习