三维空间曲线拟合
 

%function hsyc
x=[19369427.05 19365589.18 19395396.98 19378215.94 19396290.16 19372408.39];
y=[3960436.062 3959204.844 3956360.172 3965447.232 3970705.408 3964632.235];
z=[52.1399 43.6719 37.4569 50;56.3249 85.4918];
Z=interp2(x,y,z,19374774,3956111,'cubic')