上海交大数值分析课件数值分析2-6(三次样条插值).ppt

一、三次样条的产生和背景 第二章 插值法 §5 三次样条插值 二、三次样条函数的定义 三、三转角方程 四、三弯矩方程 预备知识 已知：4个条件 xi xk xk+1 yi = f(xi) yk yk+1 求:一个次数不超过3的多项式H3(x) Hermite插值： 结论： 其中 实际中有许多计算问题对插值函数的光滑性有较高的要求，例如飞机机翼外形、发动机进、排气口都要求有连续的二阶导数。 一、 三次样条的产生和背景 1.问题的产生 显然我们前面介绍的方法已不能解决这个问题。 2.样条的概念(Spline) 样条是工程设计中使用的一种绘图工具，它是富有弹性的细木条或细金属条。绘图员利用它把一些已知的点连接成一条光滑曲线称为样条曲线，样条曲线在连接点处有连续的曲率（即连续的二阶导数），它实际上是分段三次曲线拼接而成，在连接点上要求二阶导数连续。 二、 三次样条函数的定义 若函数S(x)∈C2[a,b],且在每个小区间[xj,xj+1]上是三次多项式，其中 a =x0＜x1 ＜… ＜xn=b 是给定节点，则称S(x)是节点x0，x1， … ，xn上的三次样条函数。 1.三次样条的定义 a.S(x)∈C2[a,b] b.S(x)在[xj,xj+1]上是三次多项式 即： 三次样条函数 2.三次样条插值函数的定义 + S(xi) = yi 3.求解三次样条插值函数的已知条件数和未知条件数 未知参数个数 4n 已知条件个数 插值条件： n+1 S(x)∈C2[a,b] ：3(n-1) 共 计： 4n-2 缺少条件，通常在插值区间的端点给出，称为边界条件。 4.常用的三种边界条件 1°已知两端的一阶导数值，即： 2°已知两端的二阶导数值，即： 3°当f(x)是以xn-x0为周期的周期函数时，则要求S(x)也是周期函数，即 周期样条 三、 求解方法之一：三转角方程 设在[a,b]上给出插值条件： 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 思路: (1)首先要补条件：每个区间上构造三次多项式需要四个条件，但现在最多有三个，故要补充条件，形成四个； (2)补什么条件：或函数值，或一阶导数值，或二阶导数值。这里选一阶导数较合适； (3)如何补？若随意给，则只能保证构造出的插值函数的函数值和一阶导数值连续，但不一定能保证二阶导数值连续，故只能选那组使二阶导数连续的一阶导数值。 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn x1处: 得到与m0,m1,m2有关的等式 x2处: 得到与m1,m2,m3有关的等式 共n-1个等式 设在[a,b]上给出插值条件： 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 设法求出 求解过程具体如下： 2.求解 mj 的思路 由内部节点上的二阶导数连续求出 考虑S(x)在[xj , xj+1]上的表达式 hj=xj+1-xj 对S(x)求二阶导数得： 于是 同理可得S(x)在区间[xj-1 , xj]上的二阶导数： 于是 由条件 可得 进一步简化为 写成矩阵形式为 四、求解方法之二：三弯矩方程 设在[a,b]上给出插值条件： 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 思路: (1)首先要补条件：每个区间上构造三次多项式需要四个条件，但现在最多有三个，故要补充条件，形成四个； (2)补什么条件：或函数值，或一阶导数值，或二阶导数值。这里选二阶导数较合适； (3)如何补？若随意给，则只能保证所构造出的函数的函数值和相邻两段在公共点的二阶导数的极限值连续，但不能保证一阶导数连续，故只能选那组使一阶导数连续的二阶导数值。 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn x1处: 得到与M0,M1,M2有关的等式 x2处: 得到与M1,M2,M3有关的等式 共n-1个等式 求解方法之二：三弯矩方程 设在[a,b]上给出插值条件： 1.条件 xj x0 x1 x2 … xn f(xj) f0 f1 f2 … fn 求三次样条插值函数 S(x) 2.求解S(x)的思路及求解 1)首先确定S(x)与二阶导数值的关系 2）求出中间节点上的一阶导数值 1)首先确定S(x)与二阶导数值的关系 由于S(x)在区间[xj ,xj+1]上是三次多项式， 故S〃(x)在[xj ,xj+1]上是线性函数，可表示为 对S〃(x)积分两次并利用S(xj)=yj 及S(xj+1)=yj+1