11641691625163616491664
l5x
x0x1x4x9x16x36x49x64
2502512542592516253625492564
l6x
x0x1x4x9x16x25x49x64
3603613643693616362536493664
l7x
（x0）x1x4x9x16x25x36x64
（490）49149449949164925（49364964
l8x
x0x1x4x9x16x25x36x49
6406416446496416642564366449
L8xl1x2l2x3l3x4l4x5l5x6l6x7l7x8l8x
求三次样条插值函数由MATLAB计算：
可得矩阵形式的线性方程组为：
20000
02500
0

0
00

0
0
0
02000003750
000000
0075002000004167
00000
000625020000043750000
000058332000004500000
000005625200000458300
000000550020000046340
00000005417200000
0M00
0


M
1


10

00

MM
23

0100286
0M400119
0

M
5



00061
0
0

5357
MM
67


00035
00022
20000M80
在MATLAB中的Editor中输入程序代码，以下是三次样条函数的程序代码：fu
ctio
tgsa
ci
st
代表元素数，st代表端点的一阶导。
1212
f1313
百度文库让每个人平等地提升自我
y012345678x01491625364964
9forj11
1
hjxj1xje
d
forj21
1rjhjhjhj1
e
dforj11
1uj1rje
dforj11
1fjyj1yjhje
d
forj21
1dj6fjfj1hj1hj
e
dd10d
0azeros
forj11
ajj2e
dr10u
0forj11
1aj1juj1ajj1rje
dbi
vambdtapzeros
14p矩阵为Sx函数的系数矩阵
f百度文库让每个人平等地提升自我
forj11
1
pj1mj6hj
pj2mj16hj
pj3yjmjhj26hj
pj4yj1mj1hj26hj
e
d
p
解得
M00M1M2M3M4M5M6M7
M80，则三次样条函数：
（01x）300868x03（01x）10868x0，x01

0028（94

x）3

0003（1x
1）3

0593（84

x）
0638（8x

1），x

14
0001（99x）30000（9x4）30353（59x）0627（1x4），x49
Sx
0（0020500（x6）1360x0）030（1（0xx
9）30459（016x）05708x9，x91616）30443（625x）05600x16，x1625
（036x）3（0x25）30459（936x）05470x25，x2536
（048x）3（0x36）30463（348x）05404x36，x3648
（064x）30x4830468（964x）05333x48，x4864
下面进行画图，在Comma
dWi
dow中输入画图的程序代码：
画图形比较那个插值更精确的函数：
x001491625364964
y001r