三次样条插值.docxVIP

三次样条插值 算法原理 由于在许多实际问题中，要求函数的二阶导数连续，人们便提出了三次样条插值函数，三次样条插值函数是由分段三次函数拼接而成的，在连接点处二阶导数连续。 设S(x)在节点处的二阶导数，其中为待定参数。由S（x）是分段三次多项式可知，是分段线性函数，在子区间上可以表示为 其中，对S（x）两端积分两次得 其中和为积分常数。由插值条件得 由此解得 代入得： 求导得： 令得在处的左导数 又令得在处右导数 ， 从而有，由在节点处一阶导数的连续性知，即 两端同乘得， 记，，则关于的方程组写成。 三种边界条件的三弯矩方程： （1）第一种边界条件：已知。取，这时方程组减少了两个未知量，变成只含n-1个未知量的n-1个方程的方程组，用矩阵表示为 可用追赶法求解出后，即得三条样条插值函数。 第二种边界条件，已知，记，则有，得 ，即 ，其中，得到方程组，用矩阵表示为 ，该方程组的系数矩阵是严格三对角占优矩阵，可用追赶法求解。 （3）第三种边界条件：周期型边界条件。已知是以为周期的周期函数，则由周期性可知，，这时将点看成是内节点，则有，也即 ，其中，方程组第1个方程为：，所以方程组为 ，用矩阵表示为 ，显然系数矩阵为严格对角占优矩阵，可用LU分解法求解。  程序框图 源程序 function x=mchase(A,d) %追赶法 n=length(d); u=zer