最小二乘拟合椭球系数

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];
8 T2 ?5 c8 w, _& V  o已知测量数据x y z：
- H+ N! M% F! [4 l9 e- M' R       76.561      -125.92       187.94
       321.38       193.05       133.07
       2.3395       439.42        94.38
      -201.23       119.15       138.61
' D) Q0 G' ~8 e' W. [       352.77       422.41      -86.162
+ k$ O; s1 C# `7 k9 h! H        371.4       119.72        167.8
5 G% A! T) m5 l; J' `       361.08       -106.7      -136.33
       320.15       195.71      -414.88
         66.5       -105.1      -423.75
7 V8 S7 G/ {& K# F+ L. K3 B       283.17      -214.49      -45.219
: D% d' W( {" T" g7 m3 c* Q      -14.378      -109.17       175.69
      -215.19      -79.986      -163.08

程序如下：
6 l* E% A: G0 C- txdata = [m(:,2) m(:,3) m(:,4)];
ydata = [0 0 0 0 0 0 0 0 0 0 0 0]';
$ M1 {5 E" J6 u% M! r4 u3 kk=[0 0 0 0 0 0 0 0 0 0];
' D. z7 _9 c' V- P8 x. cF = @(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);
4 E% x% |/ |! @5 _5 Rk0=ones(1,10);
[k,resnorm]=lsqcurvefit(F,k0,xdata,ydata);

8 S9 y& h- X" ^9 k8 z最后结果趋近于零，为什么？求指教

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
A: -9.52309128703986E-6
9 d  w% ?$ X( V1 K7 ~+ l* u# PB: -8.16614311365315E-6
; {; p8 A5 w* uC: -8.62656710048033E-6
D: -9.30535264389442E-7
E: 1.3203421356307E-6
: ?6 l- ~3 u6 ?" [7 V$ y, GF: -2.39934410974343E-6
( C7 \1 U+ h# O- }G: 0.00165699153341434
H: 0.0022664316097256
4 J0 h! q0 R# O* F5 I- i2 n. ?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
; x- l# R8 F7 m" VA: -9.52309128703986E-6
B: -8.166143 ...

您的意思是我的代码没有问题？是方程模型错了？
7 V% J& P3 l2 ~+ `4 s3 C9 O