中学课件数值分析三次样条插值.pptVIP

华长生制作 * 2.3.3 三次样条插值 高次插值出现龙格现象 L-插值（牛顿插值） Hermite插值 分段插值 但分段线性插值在节点处不一定光滑 分段Hermite插值 但导数值不容易提取（找到） 三次样条插值（先由函数值确定导数值，再由分段 Hermite插值解决问题） 举例： 1 汽车、船的外形设计，流体力学等要求流线型（光滑）； 2 木样条的来源。 背景 什么是样条: 样条本质上是一段一段的三次多项式拼合而成的曲线 在拼接处,不仅函数是连续的,且一阶和二阶导数也是连续的 1946年,Schoenberg将样条引入数学,即所谓的样条函数 一、三次样条插值函数 定义1. ------(1) 二、三次样条插值多项式 ------(2) ------(3) ------(4) 少两个条件 并且我们不能只对插值函数在中间节点的状态进行限制 也要对插值多项式在两端点的状态加以要求 也就是所谓的边界条件: 第一类(转角)边界条件: 第二类(弯矩)边界条件 第三类(周期)边界条件 ------(5) ------(6) ------(7) 加上任何一类边界条件(至少两个)后 一般使用第一、二类边界条件, 即 ------(8) 或 常用第二类边界条件 2、 三弯矩构造法 　　三次样条插值函数 可以有多种表达式，有时用二阶导数值 　表示时，使用更方便。 在力学上解释为细梁在 处的弯矩，并且得到的弯矩与相邻两个弯矩有关，故称用 表示 　的算法为三弯矩算法。 　对 　 积分两次,并利用插值条件定出积分常数,可以得到 　　由于 在区间 　上是三次多项式，故 在 　上是线性函数，可表示为 其中 这是三次样条插值函数的表达式，当求出 后, 就由上式完全确定. 对 求导得 由此可得 当 　　时, 　 的表达式可得,因此有 利用条件 　得 其中 上述方程组是关于 　的方程组,有 个未知数,但只有 个方程.可由前面给出的任一种边界条件补充两个方程。 对于转角边界条件，可以导出两个方程 这样可解出 ，从而得 的表达式,若令 　　　　　　　　　　　　　　 则方程组可以写成矩阵形式 对于弯矩边界条件,直接得 然后再求解矩阵方程。 对于周期性边界条件，有 可解出 ，方程组的矩阵形式为 其中 例1. 对于给定的节点及函数值 解: 得方程组 * * * * *