三次样条插值计算步骤

三次样条插值在实罩举袜际中有着广泛的应用，在计算机上也容易实现。下面介绍用计算机求取三样条插值函数S（x）的算法步骤：

（1）输入初始节点离散数据xi，yi（i＝0，1，…，n）；

（2）依据式（6－46），计算hi＝xi－xi－1，λi和Ri（i＝1，…，n－1）；

（3）根据实际问题，从式（6－49）、式（6－51）物激和式（6－53）中选择一类对应的边界条件，求取v0，w0，u0，R0，un，vn，wn，Rn；

（4）根据形成的方程组（6－54）的特点，选用追赶法、高斯法等解方程组，求出Mi（i＝0，1，2，…，n）；

（5）依据式答禅（6－41）、式（6－42），计算插值点的三样条插值函数值和该点的导数值。

void SPL(int n, double *x, double *y, int ni, double *xi, double *yi)； 是你所要。

已知 n 个点 x,y x 必须已按袭搭顺序排好。要插值 ni 点，横坐标 xi[], 输出 yi[]。

程序里用double 型，拍卖拿保证计算精度。

SPL调用现成的程序。

现成的程序很多。端点处理方法配谈不同，结果会有不同。想同matlab比较，你需 尝试 调用 spline（）函数 时，令 end1 为 1, 设 slope1 的值，令 end2 为 1 设 slope2 的值。

---