三次样条插值三弯矩matlab_三次样条（cubic spline）插值

当已知某些点而不知道具体方程时候，最经常遇到的场景就是做实验，采集到数据的时候，我们通常有两种做法：拟合或者插值。拟合不要求方程通过所有的已知点，讲究神似，就是整体趋势一致。插值则是形似，每个已知点都必会穿过，但是高阶会出现龙格库塔现象，所以一般采用分段插值。今天我们就来说说这个分段三次样条插值。

顾名思义，分段就是把区间[a,b]分成n个区间

,共有n+1个点，其中两个端点

。三次样条就是说每个小区间的曲线是一个三次方程，三次样条方程满足以下条件：

1，在每个分段小区间

上，

都是一个三次方程

2，满足插值条件，即

3, 曲线光滑，即

连续

则这个三次方程可以构造成如下形式：

这种形式,我们称这个方程为三次样条函数

。

从

可以看出每个小区间有四个未知数

，有n个小区间，则有4n个未知数，要解出这些未知数，则我们需要4n个方程来求解。

求解

我们要找出4n个方程来求解4n个未知数

首先，由于所有点必须满足插值条件，

，除了两个端点，所有n-1个内部点的每个点都满足

前后两个分段三次方程，则有2(n-1)个方程，再加上两个端点分别满足第一个和最后一个三次方程，则总共有2n个方程；

其次，n-1个内部点的一阶导数应该是连续的，即在第 i 区间的末点和第 i+1 区间的起点是同一个点，它们的一阶导数应该也相等，即

则有n-1个方程

另外，内部点的二阶导数也要连续，即