自编的三次样条插值matlab程序(含多种边界条件)txt.DOC

数值计算第二次大作业——验证三次样条函数插值是否有几何不变性（1）给定的插值条件如下： i 0 1 2 3 4 5 6 7 Xi 8.125 8.4 9.0 9.485 9.6 9.959 10.17 10.2 Yi 0.0774 0.099 0.28 0.60 0.708 1.200 1.800 2.177端点边界条件为第一类边界条件（给定一阶导数）： .Y . 0.01087 .0. Y .100 .7.三次样条函数的构造过程如下：设 x1 . x2 . x3 . xn.1 .xn 共 n个插值节点，则经过数据点 . x1, y1 . ,. x2, y2 . ,.,.xn , yn . 的三次样条 S . x. 是一组三次多项式： .S .........231111111112,,xabxxcxxdxxxxx.........,. . 23. . . a .. x . x . c .. x . d .x . x , x ..x , x ,.S2 x 2 b22 . 2 x 2 . 22 . 23 . . (1.1).. .S .x . x .3, x ..x , x .. n. .....nnnnnxabxxcxxd............. n .n.1 n.1. 由节点处的连续性可知： S . x .. y , S . x .. y ,i .1, 2, .n .1. ii iii.1 i.1 .a . yi , .1, 2, .n .1. ii . .y y ......2321121121121bxxcxxdxx......., (1.2).... .y . y . b . x . x .. c .x . x .2 . d .x . x .3 . nn.1 n.1 nn.1 n.1 nn.1 n.1 nn.1 由节点处的一阶与二阶光滑性可知： Si 11,iiiiiixSxSxS..... .. . . ..xi . ,i .1, 2, ., n. . 20 . S .x .. S .x .. b . 2c .x . x .. 3d .x . x .. b12 221121 121 2. .. . 2 (1.3) .0 . Sn.2 . xn.1 .. Sn.1 . xn.1 .. bn.2 . 2cn.2 .xn.1 . xn.2 .. 3dn.2 .xn.1 . xn.2 .. bn.1.. .0 . S . x .. S . x .. 2c . 6d .x . x .. 2c1222 1121 2 ... ..0 . Sn 211122nnnnxSxc........ .. .. 6dn.2 . xn.1 . xn.2 .. 2cn.1 又设 cn .Sn 1.. xn . 2 ，记 .. x . x , .. y . yi , .1, 2, ., n .1 ，则由(1.3)可得： ii.1 ii i.1 i c . ci.1 idi . , i . 1, 2, ., n .1. (1.4)3.i 从（1.2）解得： bi ..i . ci.i . di.i 212iiiicc..... . . ., i .1, 2, ., n .1. (1.5).i .i 3 将(1.4)与(1.5)代入(1.3)得： . ..2 .1 . ..c . 2..1 .. . 2 ..2c3 . 311 2 c ..., ...2 .1 .... (1.6). ..n.1 .n.2 ...n.2cn.2 . 2..n.2 ..n.1 .cn.1 ..n.1cn . 3..... ..n.1 .n.2 . 增加两个端点边界条件，因为 2c 1111,2nnSxcS.... . . xn . ，故有： 1.第零类边界条件：自然样条， c1 . 0, cn . 0. 2.第一类边界条件：给定端点一阶导数，设 S 1111,nxvS... .. xn ..vn ，则有： .2.c ..c . 3.. .. v ..... ..n .11121111111123nnnnnnccv........... 3.第二类边界条件：给定端点二阶导数，设 S1111,nxvS... .. xn ..vn ，则有： .2c1 . v1.2c . v.nn ci . c ..结合(1.6)及所给的边界条件即可解出 .ci .,而 di ..311,2iiiiiibcc.....i.i3 . . ，故可得到最终各个子区间上的三次样条函数。根据以上过程进行 matlab编程，编写三次样条 spline3函数，具体见附录。因此所编函数可对第一题求解： clear;clc;format short g;x1=[8.125 8.4 9 9