三次样条插值MAAB实现.docVIP

三次样条插值MAAB实现 三次样条插值MAAB实现 PAGE / NUMPAGES 三次样条插值MAAB实现 MATLAB程序设计期中考察 在很多问题中，往常依据实验、观察或经验获得的函数表或失散点上的信息， 去研究剖析函数的相关特征。此中插值法是一种最基本的方法，以下给出最基本 的插值问题——三次样条插值的基本提法： 对插值区间a,b进行区分：ax0x1xnb，函数yfx在节点 x上的值yifxii0,1,2,n，而且假如函数Sx在每个小区间xi,xi1上 i 是三次多项式，于a,b上有二阶连续导数，则称Sx是a,b上的三次样条函数， 假如Sx在节点xi上还知足条件 则称Sx为三次样条插值函数。 三次样条插值问题提法：对a,b上给定的数表以下. xx1??xn 0 yy1??yn 0 求一个分段三次多项式函数Sx知足插值条件Sxiyii0,1,n式，并在 插值区间a,b上有二阶连续导数。这就需要推导三次样条插值公式： 记fx在节点xi处的值为fximi(i0,1,n)(这不是给定插值问题数 表中的已知值)。在每个小区间xi,xi1利用三次Hermite插值公式，得三次插值 公式： Sxixyxyxmm，xxi,xi1。为了获得这个公 ii1i1iiii1i1 式需要4n个条件： (1).非端点处的界点有2n个；(2).一阶导数连续有n1个条件；(3).二阶导数 连续有n1个条件，此中界限条件：○1Sx0m0Sxnmn 此中： 0, x j i 1, i j i j i x 0 i xj 0 且( i, j 0,1)。 j 0, x j i 1, i j i j ，mi 为对应变量的一阶导数。其推导过程以下： 为了确立mi的值，把Sx睁开为： 22 xxxxxxxx im 1iii1 +m, ii1 22 hh ii 这里hixi1xi，对Sx连续求两次导，得： 6x 2x 4x 1 6x 4x 2x 6 x x 2x i i y y i i 1 i i 1 S x m m 1 i i 1 i i 2 2 3 h h h i i i 。于是 考虑Sx在节点xi处的右极限值，得： 426 Sxi0mmy1y。 ii1ii 2 hh h ii i 同理，在相邻小区间xi1,xi上可得Sx的表达式为： 及Sx在节点xi处的左极限值为： 246 Sxi0mmyiyi。利用Sx二阶导数于节点xi处 i1i1 2 hh h i1i1i1 的连续性条件Sxi0Sxi0，这里i1,2,n1，有下式建立： 1 1 1 1 y y y i 1 i i m 2 m m 3 i 1 i i 1 2 hi h h h h h 1 i 1 i i i i y i 2 1 1 ，用 1 1 hi h 1 i 除等式两 yy i1i, 边，并注意yif,fxixi1，上式可简记为： i h i hh ii1 且11 ii hhhh i1ii1i 最后求得m1mn的线性方程组为： 2 0 0 0 1 1 m 1 g 1 2 2 0 0 0 2 m 2 g 2 (**) 0 0 0 n 2 1 n 1 m n 1 g n 1 n 0 0 0 2 n m n g n 经过以上复杂的求解和迭代，就能够求解出插值函数的近似表达式。得出来的表 达式就能够用MATLAB 软件来求解。详细求解过程以下： 已知n对数据点x1,y1,x2,y2,x3,y3,xn,yn,，假定函数关系为 yfx，但分析式不确立，数据插值就是结构函数关系式ygx，使 xii1,2,3,,n，知足关系gxifxi。 例题：求知足下边函数表所给出的插值条件的三次自然样条函数。 1245 1342 剖析：表中所列出的是函数对点，第一要把对应的插值函数求出来，再用 MATLAB软件来求区间1,5上间隔为0.5的各点的值。 求解过程以下： 因自然样条插值函数的界限条件为 这里n3，故确立m0,m1,m2,m3的方程组形式形如上边的(**) 式，此中系数i,i 和gi可按以下步骤进行： 将上述参数带入(**)式，获得以下方程组： 解得： 由公式 22 xxxxxxxx i1m iii1 +m1, ii 22 hh ii 可知， 由所求出的表达式可知区间1,5可分为1,44,5，对两个区间分别用MATLAB 命令即可： 针对第一个区间： 13227 yxxx1;其图像以下 884 命令以下： x=1:4;y=(-1/8)*x.^3+(2/8)*x.^2+(7/4)*x-1;xi=1:0.5:4; y1=interp1(x,y,xi,spline) 其运转结果以下： y1 = Columns 1 through 6 0.87501.76562.50002.98443.12502.8281 Column 7 2