第1章1-08样条插值.ppt

* §1.8 三次样条插值 1.8.1 引言 高次插值出现龙格现象 L-插值（牛顿插值） Hermite插值 分段插值 但分段线性插值在节点处不一定光滑 分段Hermite插值 但导数值不容易提取（找到） 三次样条插值（先由函数值确定导数值，再由分段 Hermite插值解决问题） 一、背景 举例： 1 汽车、船的外形设计，流体力学等要求流线型（光滑）； 2 木样条的来源。 定义 （3次样条函数） 在每一个小区间 上是次数 多项式。 若（1）中3次样条函数 还满足插值条件： 关于分划 称 为 的3次样条插值函数。 ，即具有连续的一阶，二阶导数。 满足下述条件： 如果函数 (1) 设有对[a,b]的分划 的一个3次样条函数。 为关于分划 则称 提出问题： 3次样条插值函数 是否存在？是否唯一？ 如何计算？误差估计？ 函数表 （2）设给定 二 、 样条函数的定义 分析： 因 上是3次多项式，即为 已有条件： 4n个待定系数: 个条件 内部条件： (内节点) 个条件 共有 个条件, 要唯一确定 ，还必须附加２个条件(边界条件)。 常见边界条件有三种： 第1种边界条件： 第2种边界条件： 若 ，称为自然边界条件。 已知,即 已知,即 为周期函数， 第3种边界条件（周期边界条件）： 此时称 为周期样条函数。 亦是周期函数，周期为 ，即取 要求 注： 一般不取一端是一阶导数而另一端是二阶导数。 三次样条插值函数的表达式 基本思路： 以分段三次Hermite插值为基础，由(1)函数表 ；(3)三种边界条件中 的某一种推导3次样条插值函数。 1、先确定插值函数 在节点处的一阶导数，记为 该方法即为3次样条插值函数的一阶导数表示。 2、先确定插值函数 在节点处的二阶导数，记为 该方法即为3次样条插值函数的二阶导数表示。 推导方法： （一阶导数表示） 1.8.2 三次样条插值函数的表达式 一、推导公式： 回忆： 分段3次Hermite插值 已知 问题： 求3次样条插值函数 。 函数表 设给定 不固定,是待定参数，共(n+1)个 若要 则要求满足： 个条件 加某一边界条件（2个） 个条件 令 对(8. 3)式求导: 由条件 及（8.5)、(8.6)得 把(8.4)代入(8.7)得到 所满足线性方程组 两边同除以 ，得 令 得 说明： （b）（8.8）式有n-1个方程,要确定n+1个未知量 还少两个方程,由边界条件补足. 方程组成的方程组.mj( j=0,1,…,n)在力学上叫做细梁xj( j=0,1,…,n) 的 （a）(8.8)式是关于n+1个未知量 个 三种边界条件 处的转角，数学上叫做变化率。方程(8. 8)反映了mj与mj-1,mj+1的 关系,因此(8.8)叫做三转角方程。 方程组(8. 8)为关于 所满足的方程组： (1)增加第1种边界条件： 则方程组(8. 8)为关于 所满足的方程组可写为： (1)增加第1种边界条件： 则方程组(8. 8)为关于 所满足的方程组： 再由 (2)增加第2种边界条件： 则由(8.5)式取 及(8.6)式取 得到2个方程 (利用(8.4)式中 的值）: 由下式 得( j = 0, j=n-1) 把(8.4)/式分别代入（**），得 由(8.4)式得 整理得两个方程: 于是,得到 所满足的线性方程组： 上式简记为 整理得两个方程: (3)增加第3种边界条件： 即得 又由条件 利用（*）式得到 令 得 两边同除以 ，得 所满足线性方程组: 此式加到（8.8）式即得 异矩阵，则方程组（8. 9）、（8. 10）或(8. 11)有唯一解 可用追赶法求解,从而由(8. 3)给出 表达式,且 具有连续的一阶,二阶导数(即 为3次样条插值函数)。 说明： 方程组（8. 9）、（8.10）或（8.11）系数矩阵都是严格 对角占优矩阵，这些方程组的系数阵为非奇 , 二、存在唯一性 且 (1)如果 是定义在 上函数且已知 函数表 定理 （3 次样条插值函数存在唯一） (2)给定边界条件 ，则 于 存在 唯一3次样条插值函数 ,且满足 三、计算步骤 : 先计算(8. 8)式中的 （若是第二类 或第三类边界条件，要计算 ）。 （2）用追赶法求解方程组(8. 9)(或