三次样条插值的MAAB实现.docxVIP

MATLAB 程序设计期中考查 在 多 中，通常根据 、 或 得到的函数表或离散点上的信息， 去研究分析函数的有关特性。 其中插 法是一种最基本的方法， 以下 出最基本 的插 ——三次 条插 的基本提法： 插 区 a,b 行划分： a x0 x1 xn b ，函数 y f x 在 点 xi 上的 yi f xi i 0,1,2, n ，并且如果函数 S x 在每个小区 xi , xi 1 上 是三次多 式，于 a,b 上有二 数， 称 S x 是 a,b 上的三次 条函数， 如果 S x 在 点 xi 上 足条件 称 S x 三次 条插 函数。 三次 条插 提法： a, b 上 定的数表如下 . x0 x1 ?? xn y0 y1 ?? yn 求一个分段三次多 式函数 S x 足插 条件 S xi yi i 0,1, n 式，并在 插 区 a,b 上有二 数。 就需要推 三次 条插 公式： f x 在 点 xi 的 f xi mi ( i 0,1, n )( 不是 定插 数 表中的已知 )。在每个小区 xi , xi 1 利用三次 Hermite 插 公式，得三次插 公式： S x  i x yi i 1 x yi 1 i xi mi i 1mi 1 ， x xi , xi 1 。 了得到 个公 式需要 4n 个条件： (1).非端点 的界点有 2n 个； (2).一 数 有 n 1 个条件； (3). 二 数 有 n 1 个条件，其中 界条件：○1 S x0 m0 S xn mn 其中： i x j 0 , i j i x j0 i x j0且 ( i , j 0,1)。 1, i j i x j 0, i j ， mi 为对应变量的一阶导数。其推导过程如下： 1, i j 为了确定 mi 的值，把 S x 展开为： + xxi 1 2 x xi mi x xi 2 x xi 1 mi 1 , hi 2 hi 2 这里 hi xi 1 xi ，对 S x 连续求两次导，得： S x 6x 2xi 2 4xi 1 mi 6x 4xi 2 2xi 1 mi 1 6 xi xi 3 1 2x yi 1 yi 。于是 hi hi hi 考虑 S x 在节点 xi 处的右极限值，得： S xi 4 2 6 yi 1 yi 。 0 mi mi 12 hi hi hi 同理，在相邻小区间 xi 1, xi 上可得 S x 的表达式为： 及 S x 在节点 xi 处的左极限值为： S xi 0 2 mi 1 4 mi 6 2 yi yi 1 。利用 S x 二阶导数于节点 xi 处 hi 1 hi 1 hi 1 的连续性条件 S xi 0 S xi 0 ，这里 i 1,2, n 1 ，有下式成立： 1 1 1 1 yi 1 yi yi yi 1 1 1 i 1 mi 1 2 i 1 i mi i mi 1 3 2 2 ，用 hi 1 hi 除等式两 h h h h hi hi 1 边，并注意 yi fi , yi 1 yi f xi , xi 1 ，上式可简记为： hi 且 i hi i 1 i hi 1 gi 3 i f xi 1, xi i f xi , xi 1 hi hi 1 hi hi 1 最后求得 m1 mn 的线性方程组为： 2 1 0 0 0 1 m1 g1 2 2 2 0 0 0 m2 g 2 (**) 0 0 0 n 1 2 n 1 mn 1 g n 1 n 0 0 0 n 2 mn g n 通过以上复杂的求解和迭代， 就可以求解出插值函数的近似表达式。 得出来的表 达式就可以用 MATLAB 软件来求解。具体求解过程如下： 已知 n 对数据点  x1 , y1  , x2  , y2  , x3 , y3  ,  xn  , yn  , ，假设函数关系为 y f x ，但解析式不确定，数据插值就是构造函数关系式  y g x ，使 xi  i  1,2,3,  , n  ，满足关系  g xi  f xi  。 例题：求满足下面函数表所给出的插值条件的三次自然样条函数。 1 2 4 5 1 3 4 2 分析：表中所列出的是函数对点，首先要把对应的插值函数求出来，再用 MATLAB 软件来求区间 1,5 上间隔为 0.5 的各点的值。 求解过程如下： 因自然样条插值函数的边界条件为 这里 n 3，故确定 m0 , m1