matlab code
y=[432441.33 433215.293 433482.888 434025.39 433350 434170 430160 431320 432020 432375 435113 434825 436045 435145 434760 435310 435650 431820 432285 433525 434400];
z=[197.689 231.15 224.16 252.58 196 194 215 215 181.5 182.2 242 185 218.2 185 240 233 226 188 194.5 242 242];
[X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),50)',linspace(min(y),max(y),50),'v4');%插值
sfun2d = Bspline(Z,3);%拟合曲面 需安装 bsspline tools
surf(X,Y,double(sfun2d))%作图
hold on;plot3(x,y,z,'ko')
最小二乘法逼近用sftool就可以多项式3次就可拟合的差不多