三次样条差值

1计算原理和程序流程图

设在区间[a, b]上给定n+1个节点xi（a≤x₀<x₁< …< x_n≤b）,在节点x_i处的函数值为y_i=f(x_i) (i= 0,1,…,n)。若函数S(x)满足以下三个条件：

(1) 在每个子区间[x_i-1, x_i] (i= 0,1,…,n)上，S(x)是三次多项式；

(2) S(x_i) = y_i(i=0,1,…,n)；

(3) 在区间[a, b]上，S(x)的二阶导数S”(x)连续。

则称S(x)为函数y_i= f(x)在区间[a,b]上的三次样条插值函数。

由定义可知S(x)共有4n个待定参数，根据条件(3)可得如下3n-3个方程，再加上条件(2)的n+1 个方程，共有4n-2 个方程，因此还需要增加两个条件才可以确定出4n 个待定参数。所增加的条件称为边界条件或端点条件。常用的三种边界条件为：①已知 f”(a), f”(b)；②已知 f”(a)= f”(b)；③已知 f(x)是以 T=b-a 为周期的周期函数。

三次样条插值计算流程图如下：