令：clearclfholdo
X11Y1125X2CD
ewpolyXYx11ypolyvalCxplotxyXYgrido
xp11z1125xp2plotxpzr得到插值函数和f（x）图形：
88
f百度文库让每个人平等地提升自我
下面再求三次样条插值函数，在MATLAB的Editor中建立一个多项式的Mfile输入下列程序代码：fu
ctio
ScsfitXYdx0dx
Nle
gthX1HdiffXDdiffYHAH2N1B2H1N1H2NCH2NU6diffDB1B1H12U1U13D1BN1BN1HN2UN1UN13DNfork2N1
tempAk1Bk1BkBktempCk1UkUktempUk1e
dMNUN1BN1forkN21199
f百度文库让每个人平等地提升自我
Mk1UkCkMk2Bke
dM13D1dx0H1M22MN13dx
DNHNMN2fork0N1
Sk11Mk2Mk16Hk1Sk12Mk12Sk13Dk1Hk12Mk1Mk26Sk14Yk1e
d当
10时，我们在Comma
dWi
dow中输入以下的命令：clearclcX11Y125X21dx0dx
ScsfitXYdx0dx
x11y1polyvalS1x1X1x20y2polyvalS2x2X2x30y3polyvalS3x3X3x41y4polyvalS4x4X4plotx1y1x2y2x3y3x4y4XY结果如图：
当
20时，我们在Comma
dWi
dow中输入以下的命令：1010
f百度文库让每个人平等地提升自我
clearclcX11Y125X21dx0dx
ScsfitXYdx0dx
x11y1polyvalS1x1X1x20y2polyvalS2x2X2x30y3polyvalS3x3X3x41y4polyvalS4x4X4plotx1y1x2y2x3y3x4y4XY结果如图：
实验三：
下列数据点的插值x01491625364y01234567
可以得到平方根函数的近似，在区间064上作图。（1）用这9各点作8次多项式插值L8x（2）用三次样条（自然边界条件）程序求S（x）。从结果看在064上，那个插值更精确；在区间01

L8x可由公式L
xyklkxj得出。k0
三次样条可以利用自然边界条件。写成矩阵：
1111
f百度文库让每个人平等地提升自我
20000M0
1
2
1
0
0

M1

0

0
20
23
22
03

MM
23


00042M4
其中jhj1，i
hj
，dj6fxj1xjxj1，
00d0d
0
hj1hj
hj1hj
l0xx1x4x9x16x25x36x49x64010409016025036049064
l1xx0x4x9x16x25x36x49x64101419116125136149164
l2xx0x1x9x16x25x36x49x64404149416425436449464
l3xx0x1x4x16x25x36x49x64909194916925936949964
l4x
x0x1x4x9x25x36x49x64
16016r