文档介绍：一、三次样条的产生和背景
第二章插值法
§6 三次样条插值
二、三次样条函数的定义
三、三转角方程
四、三弯矩方程
实际中有许多计算问题对插值函数的光滑性有较高的要求,例如飞机机翼外形、发动机进、排气口都要求有连续的二阶导数。
一、三次样条的产生
显然我们前面介绍的方法已不能解决这个问题。
二、三次样条函数的定义
若函数S(x)∈C2[a,b],且在每个小区间[xj,xj+1]上是三次多项式,其中
a =x0<x1 <…<xn=b
是给定节点,则称S(x)是节点x0,x1, …,xn上的三次样条函数。

(x)∈C2[a,b]
(x)在[xj,xj+1]上是三次多项式
即:
三次样条函数

+
S(xi) = yi

未知参数个数
4n
已知条件个数
插值条件: n+1
S(x)∈C2[a,b] :3(n-1)
共计: 4n-2
缺少条件,通常在插值区间的端点给出,称为边界条件。

1°已知两端的一阶导数值,即:
2°已知两端的二阶导数值,即:
3°当f(x)是以xn-x0为周期的周期函数时,则要求S(x)也是周期函数,即
周期样条
三、求解方法之一:三转角方程
设在[a,b]上给出插值条件:

xj
x0
x1
x2
…
xn
f(xj)
f0
f1
f2
…
fn
求三次样条插值函数 S(x)
xj
f2
x2
fn
xn
…
f1
f0
f(xj)
…
x1
x0
思路: (1)首先要补条件:每个区间上构造三次多项式需要四个条件,但现在最多有三个,故要补充条件,形成四个;
x1处:
得到与m0,m1,m2有关的等式
x2处:
得到与m1,m2,m3有关的等式
共n-1个等式
(2)补什么条件:这里选一阶导数较合适;
(3)如何补?若随意给,则只能保证构造出的插值函数的函数值和一阶导数值连续,但不一定能保证二阶导数值连续,故只能选那组使二阶导数连续的一阶导数值。
设在[a,b]上给出插值条件:

xj
x0
x1
x2
…
xn
f(xj)
f0
f1
f2
…
fn
求三次样条插值函数 S(x)
设法求出
求解过程具体如下:
mj 的思路
由内部节点上的二阶导数连续求出
考虑S(x)在[xj , xj+1]上的表达式
hj=xj+1-xj
对S(x)求二阶导数得:
于是