最小二乘拟合椭球系数

IBB-EUT

椭球方程为一般二次曲面方程：f=a*x*x+b*y*y+c*z*z+d*2*x*y+e*2*x*z+f*2*y*z+p*2*x+q*2*y+*r*2*z+g;
需要拟合出系数[a b c d e f p q r g];
: ?& a7 J. ?0 j; X3 R已知测量数据x y z：
       76.561      -125.92       187.94
       321.38       193.05       133.07
       2.3395       439.42        94.38
" @* X5 E$ |* {5 m      -201.23       119.15       138.61
       352.77       422.41      -86.162
        371.4       119.72        167.8
       361.08       -106.7      -136.33
9 R1 [2 s& Q5 O) t% {3 f$ {5 K       320.15       195.71      -414.88
% N* y6 U& b) d  _" ]( m         66.5       -105.1      -423.75
       283.17      -214.49      -45.219
' ^! F. q0 S7 ~9 ?9 ~      -14.378      -109.17       175.69
' O* Y8 Q) E1 Z. M: [      -215.19      -79.986      -163.08

, g4 o& P; m+ S4 z: Q程序如下：
5 m3 e) P) O  `" Sxdata = [m(:,2) m(:,3) m(:,4)];
( O6 X& z5 F: F; sydata = [0 0 0 0 0 0 0 0 0 0 0 0]';
7 A. G" u2 f* B) u$ Pk=[0 0 0 0 0 0 0 0 0 0];
! y. _$ N0 B4 ~: mF = @(k,xdata) k(1)*xdata(:,1).*xdata(:,1) + k(2)*xdata(:,2).*xdata(:,2) + k(3)*xdata(:,3).*xdata(:,3) + k(4)*2*xdata(:,1).*xdata(:,2) + k(5)*2*xdata(:,1).*xdata(:,3) + k(6)*2*xdata(:,2).*xdata(:,3) + k(7)*2*xdata(:,1) + k(8)*2*xdata(:,2) + k(9)*2*xdata(:,3) + k(10);
  Z, g% n+ J5 s2 u7 sk0=ones(1,10);
9 F1 a0 _& U3 J9 w# P2 j[k,resnorm]=lsqcurvefit(F,k0,xdata,ydata);
1 O/ G, B. z/ |, S; ?2 K! _7 A/ s
! |$ c% b( W0 y1 y6 ?( R( y最后结果趋近于零，为什么？求指教
0 B$ R) w% J% }& {0 E

zaiyiaaaa

帮你顶一下

IBB-EUT

等大神帮忙

nkkopd

参考下：A*x^2+B*y^2+C*z^2+D*x*y+E*x*z+F*y*z+G*x+H*y+K*z+1 = 0
- r. j0 j$ \5 D/ h7 K" n, LA: -9.52309128703986E-6
8 r" P( f0 V+ \) uB: -8.16614311365315E-6
" t1 b# N" ~0 \0 Q1 A5 rC: -8.62656710048033E-6
3 e7 _, e9 ?0 x! G/ U1 S6 bD: -9.30535264389442E-7
E: 1.3203421356307E-6
F: -2.39934410974343E-6
G: 0.00165699153341434
& e$ s; V3 e9 H5 QH: 0.0022664316097256
! y$ D; B1 [1 V! \K: -0.00210652935033001

IBB-EUT

nkkopd 发表于 2020-7-15 14:21
参考下：A*x^2+B*y^2+C*z^2+D*x*y+E*x*z+F*y*z+G*x+H*y+K*z+1 = 0
0 V6 F/ F9 o; o' K, C1 I/ OA: -9.52309128703986E-6
B: -8.166143 ...

您的意思是我的代码没有问题？是方程模型错了？