我不知道你的方程，我粗略得到一个R2 = 0.99以上的方程，代码及结果如下

高次多项式拟合代码 ：

load shuju.mat

B=fit([Dotx,Dotx],Dotz,'poly43')

结果输出如下：

Linear model Poly43:

B(x,y) = p00 + p10*x + p01*y + p20*x^2 + p11*x*y + p02*y^2 + p30*x^3 + p21*x^2*y

+ p12*x*y^2 + p03*y^3 + p40*x^4 + p31*x^3*y + p22*x^2*y^2

+ p13*x*y^3

Coefficients (with 95% confidence bounds):

p00 =    2.03e+27  (-2.866e+54, 2.866e+54)

p10 =   7.369e+39  (-3.292e+66, 3.292e+66)

p01 =  -7.369e+39  (-3.292e+66, 3.292e+66)

p20 =  -1.118e+28  (-2.253e+54, 2.253e+54)

p11 =  -3.039e+42  (-2.186e+69, 2.186e+69)

p02 =   3.039e+42  (-2.186e+69, 2.186e+69)

p30 =   -4.81e+14  (-1.609e+41, 1.609e+41)

p21 =   -3.18e+29  (-1.196e+56, 1.196e+56)

p12 =    3.18e+29  (-1.196e+56, 1.196e+56)

p03 =   4.851e+14  (-1.52e+41, 1.52e+41)

p40 =   3.179e+12  (-4.631e+39, 4.631e+39)

p31 =  -4.068e+12  (-4.657e+39, 4.657e+39)

p22 =  -3.629e+12  (-1.633e+39, 1.633e+39)

p13 =   4.518e+12  (-1.861e+39, 1.861e+39)，