三次样条插值

1.样条函数

有一函数满足$S(x)\in C^2[a,b]$（在$[a,b]$上该函数二阶导数连续）

2.三次样条插值

$x_0,x_1,\cdots,|x_n是$在区间$[a,b]$上给定的节点，若有三次样条函数$S(x)$满足$S(x_i)=y_i$，则$S(x)$为三次样条插值函数。

$$S(x)=\begin{cases} 不超过三次的多项式（共有4个系数）,[x_0,x_1]\\ 不超过三次的多项式（共有4个系数）,[x_1,x_2]\\ \vdots\\ 不超过三次的多项式（共有4个系数）,[x_{n-1},x_n]\\ \end{cases}$$
上式共有 $4n$个系数，所以要求所有系数，共需 $4n$个方程。
其中有，
$S(x_i)=y_i$, —- $(n+1)$个方程
$S(x_i-0)=S(x_i+0)$, —- $(n-1)$个方程(函数连续)
$S'(x_i-0)=S'(x_i+0)$, —- $(n-1)$个方程(一阶导数连续)
$S''(x_i-0)=S''(x_i+0)$, —- $(n-1)$个方程(二阶导数连续)
以上共有 $4n-2$个方程。
一般给出以下三种边界条件：
(1).已知两端的一阶导数值， $S'(x_0)=f'_0,S'(x_n)=f'_n$
(2).已知两端的二阶导数值， $S''(x_0)=f''_0,S''(x_n)=f''_n$
(3).周期边界
确定参数$a,b,c,d,e$的关系，使得函数$S(x)$是3次样条函数
$$S(x)=\begin{cases} a(x-2)^2+b(x-1)^3,x\in(-\infty,1)\\ c(x-2)^2,x\in[1,3)\\ d(x-2)^2+e(x-3)^3,x\in[3,+\infty)\\ \end{cases}$$
在 $x=1$处，有
$$\begin{cases} S(1-0)=S(1+0) \\ S'_1(1)=S'_2(1)\\ S''_1(1)=S''_2(1)\\ \end{cases}$$
解得，a=c;
在 $x=3$处也有三个式子，解得，c=d
所以a=c=d。