三次样条插值MATLAB实现.docxVIP

精品文档 精品文档 精品文档 在多中，往常根据、或得到的函数表或离散点上的信息， 去研究剖析函数的相关特性。其中插法是一种最基本的方法， 以下出最基本 的插——三次条插的基本提法： 插区 a,b 行区分：a x0 x1 xn b，函数y fx在点 xi上的yi f xi i 0,1,2, n ，并且如果函数Sx 在每个小区xi,xi1上 是三次多式，于 a,b上有二数，称Sx是a,b 上的三次条函数， 如果Sx在点xi 上足条件 Sxi yi i 0,1, n 称Sx三次条插函数。 三次条插提法： a,b上定的数表如下. x x0 1 ?? xn x y y0 y1 ?? yn 求一个分段三次多式函数Sx 足插条件Sxi yi i 0,1, n式，并在 插区 a,b上有二数。就需要推三次条插公式： f x在点xi的f xi mi(i 0,1, n)(不是定插数 表中的已知)。在每个小区 xi,xi1利用三次Hermite插公式，得三次插 公式： Sx  i xyi i 1 xyi1 i ximi i 1mi1，x xi,xi1。了得到个公 式需要4n个条件： (1).非端点的界点有 2n个；(2).一数有n 1个条件；(3).二 数有n 1个条件，其中界条件：○1 Sx0 m0Sxnmn 2 0S xn ○Sx0 3 4Sx0 0 5Sxn 6Sxn 1 ○3Sx0 ○y0 yn Sx0 0 Sxn 0 Sx0 0 Sxn 0 4 其中： xj 0, i j i xj 0 i xj 0且(i,j0,1)。 i 1, i j ixj 0, ij ，mi为对应变量的一阶导数。其推导过程如下： 1, i j 为了确定mi的值，把Sx展开为： Sx xxi1 2hi 2xxi yi xxi 2hi 2xi1 x yi1 hi 3 hi 3 + xxi1 2xximi xxi 2xxi1 mi1, hi 2 hi 2 这里hi xi1 xi，对Sx连续求两次导，得： Sx 6x2xi 2 4xi1mi 6x4xi 2 2xi1mi1 6xi xi 3 12x yi1yi。于是 hi hi hi 考虑S x在节点xi处的右极限值，得： Sxi 0 4 2 6 yi。 mi mi1 2yi1 hi hi hi 同理，在相邻小区间 xi1,xi 上可得S x的表达式为： 6x2xi1 4xi mi1 6x4xi1 2xi mi 6xi1xi 2x Sx 2 hi 2 2 yiyi1 hi 1 1 hi1 及S x在节点xi处的左极限值为： S xi 0 2 mi 4 mi 6 yi1 。利用Sx二阶导数于节点xi处 hi1 1 2yi hi 1 hi1 的连续性条件S xi 0 S xi 0，这里i 1,2, n 1，有下式建立： 1 1 1 1 mi1 yi1 yi yi yi1 ，用 1 1 除等式两 mi 1 2 hi mi 3 2 hi12 hi1 hi hi1 hi1 hi hi 边，并注意yi fi,yi1 yi fxi,xi1 ，上式可简记为： hi imi12mi imi1gi i 1,2,n 1, 且i hi 1 hi 1 gi 3ifxi1 ,xi ifxi,xi1 hi 1hi i i hi1 hi 最后求得 m1 mn的线性方程组为： 2 1 0 00 1 m1 g1 2 2 2 00 0 m2 g2 (**) 000 n12 n1 mn1 gn1 n00 0 n2 mn gn 经过以上复杂的求解和迭代， 就能够求解出插值函数的近似表达式。 得出来的表 达式就能够用MATLAB软件来求解。详细求解过程如下： 已知n对数据点 x1,y1,x2,y2,x3,y3, xn,yn,，假定函数关系为 y fx，但解析式不确定，数据插值就是结构函数关系式y gx，使 xii 1,2,3, ,n，知足关系gxi fxi。 例题：求知足下面函数表所给出的插值条件的三次自然样条函数。 x 1 2 4 5 yfx 1 3 4 2 剖析：表中所列出的是函数对点，首先要把对应的插值函数求出来，再用 MATLAB软件来求区间1,5上间隔为的各点的值。 求解过程如下： 因自然样条插值函数的边界条件为 S x0 Sxn 0, 这里n 3，故确定m0,m1,m2,m3的方程组形式形如上面的 (**)式，其中系数 i,i 和gi可按如下步骤进行： hi: h0 211, h1 42 h2541; i: 1 h0 h1 2, 2 h1 h2 1; h1 3 h2 3 i: 1 1 1 1, 2 1 2 2; 3 3 gi: g1 31fx1,x2 1fx0,x1 9 ; 2 g2 32fx2,x3 2fx1,x2 7; 2 g0 3fx0,x1 h0y