三维三次样条插值计算三坐标

 clear;
x=[2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3];                      %测量数据
y=[0 5 15 25 35 45 55 65 75];
z=[0.984 0.983 1.000 1.007 1.007 1.007 1.007 1.006 1.005 1.003 1.001 0.998;
1.137 1.137 1.154 1.161 1.161 1.161 1.161 1.160 1.158 1.157 1.154 1.151;
1.444 1.444 1.463 1.469 1.470 1.470 1.469 1.468 1.466 1.464 1.461 1.457;
1.751 1.752 1.771 1.778 1.778 1.779 1.778 1.776 1.744 1.771 1.768 1.763;
2.057 2.059 2.079 2.086 2.087 2.087 2.086 2.084 2.081 2.077 2.074 2.069;
2.365 2.366 2.387 2.394 2.395 2.395 2.394 2.392 2.389 2.384 2.380 2.373;
2.670 2.673 2.696 2.702 2.702 2.703 2.702 2.700 2.695 2.691 2.686 2.678;
2.977 2.980 3.004 3.010 3.011 3.011 3.010 3.007 3.002 2.997 2.991 2.984;
3.284 3.287 3.312 3.318 3.319 3.319 3.317 3.314 3.309 3.303 3.297 3.289];     
x1=2:2:120;                                                                             %生成插值数据矩阵
y1=0:2:75;
[x2y2]=meshgrid（x1y1）;
x3=[2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3
    2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3
    2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3
    2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3
    2.7 13.9 26.1 36.9 47.3 56.5 66.8 76.9 87.9 97.4 106.9 117.3