很久以前我在Mathematica写了一个三次样条包.这是我将该包翻译成Matlab的.注意我在大约7年内没有看过三次样条函数,所以我基于我自己的文档.你应该检查我说的一切.
基本问题是我们给出了n个数据点(x(1),y(1)),…,(x(n),y(n)),我们希望计算分段三次插值.插值定义为
S(x) = { Sk(x) when x(k) <= x <= x(k+1)
{ 0 otherwise
这里Sk(x)是形式的三次多项式
Sk(x) = sk0 + sk1*(x-x(k)) + sk2*(x-x(k))^2 + sk3*(x-x(k))^3
样条曲线的属性为:
>样条曲线通过数据点Sk(x(k))= y(k)
>样条在端点处是连续的,因此在插值间隔Sk(x(k 1))= Sk 1(x(k 1))中的任何地方都是连续的
>样条曲线具有连续的一阶导数Sk'(x(k 1))= Sk 1′(x(k 1))
>样条具有连续的二阶导数Sk”(x(k 1))= Sk 1”(x(k 1))
要从一组数据点构造三次样条,我们需要求解系数
对于每个n-1个三次多项式,sk0,sk1,sk2和sk3.这是总共4 *(n-1)= 4 * n-4个未知数.属性1提供n个约束,属性2,3,4每个提供额外的n-2个约束.因此,我们有n 3 *(n-2)= 4 * n – 6个约束和4 * n – 4个未知数.这留下了两个自由度.我们通过在起始节点和结束节点处设置二阶导数等于零来修复这些自由度.
设m(k)= Sk”(x(k)),h(k)= x(k 1) – x(k)和d(k)=(y(k 1) – y(k))/ H(K).下列
三期复发关系成立
h(k-1)*m