分段三次样条插值机算方法.pdfVIP

苎三±=星±里!墨望：：!竺：!：：婴竺：堡璺苎垄：皇三：堡璺些至型堑兰查墨望 一 分段三次样条插值的机算方法 李迎春 (天津市电子计算机研究所300060 摘要：利用分段三次样条插值方法可以构造一个整体上具有二阶连续导数的插值函数，上机找到求解的通用 方法和步骤。 在工程计算中，为满足实际问题的要求，需要将一些指定点连接成一条光滑曲线。利用样条插值的方 法即可达到目的。由于分段三次样条插值的方法计算简单，所构造的函数的整体光措性好。而且能满足一 般实际问题的要求。所以有很高的实用价值，本文给出用分段三次样条插值的方法由一些指定点构造插值 函数的求解过程及上机计算源程序。 1、预备知识 1．1、将一些指定点连接成一条光滑曲线的实际问题可数学化为： f(xi)=y。(f：0,1，2，…．。)，构造一个三次样条插值函数s(工)，使其满足下列条件： 2)在每个小区间kf，‘+l】上是一个三次多项式。f-0，1，2，…，咒一1； 1)s∽)：y。。i_0．I，2，…。n； 2i=1，2，…，H一1。 3)s(工)在k，6】上具有二阶连续导数，即j‘‘’(百)=s‘‘’(xj)，k=o…1 1．2、构造分段三次样条插值函数的过程： 、’ 函数工(曲的计算公式： JHI—xj工j—jf十1 工f一。f+I工f+i一，‘ 2Ⅲ，+- J 一^i +(z—zt)(篇1)2m·+(x—zXI—ZJ+ r+-)(意I+1 2 m ‰一 =p 0 ％mH 一 =P+风孤‰ 盯= n +1三“ m卜 一n 其中， 第二十一届中国(天津)’2007 IT、闻络、信息技术、电子、仪器仪表创新学术会议 a。2-，a·=i：；晕i，cr=-，：一．n一-，，。。=。， hf=JJ“一zf，(i=1,2，．．^一1)． 2mB+口。mf：芦o (1一口I)mo+2ml+alm 即 2；13 (1一a．一1)ntn一2+2m。一I+口。一lmn=卢n一1 5 l 13一 (1一at)’h爿+；h器曲，一in 方程组系数矩阵为强三角矩阵，其行列式不为0，所以方程组解存在且唯一。 对方程组求解，可导出递推公式：mi=4；m。+以(f=n，n—l，．．…IO) 其中： 铲面羔石^=譬云等鲁撕叱…峨no-一誓^=譬 代人s(j)即得所求三次样条插值函数。 2、求分段三次样条插值函数的函数值的机算步骤 当z≠而(f=0,1…2．珊)时，计算s(曲函数值的算法归纳如下： 时计算：包=k，一薯“=0,1…2∽一1)，并赋值给对应数组元素日(f)(f=0,1…2，Jl—1)： 2．计算含，‰，砚，．．Jn。的线性方程组的系数％，常数项晟(i=0,1，2，…，n) a。=t·a·=j：；气言‘t=-，：．．．．一一t)，a。；。， 将讼i，成赋值给对应数组元素A(f)，丑(D“=0,1…2．．，以)： ’ 24-(1一tl J)aI-I