最小二乘拟合椭球系数

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];
已知测量数据x y z：
       76.561      -125.92       187.94
       321.38       193.05       133.07
$ _8 f1 B) n# @7 ?       2.3395       439.42        94.38
9 ]# w1 Q1 {4 W4 C5 n7 W      -201.23       119.15       138.61
# U. {/ U  A3 D4 q# v, p4 L       352.77       422.41      -86.162
        371.4       119.72        167.8
       361.08       -106.7      -136.33
       320.15       195.71      -414.88
         66.5       -105.1      -423.75
       283.17      -214.49      -45.219
      -14.378      -109.17       175.69
      -215.19      -79.986      -163.08
1 b- o% y- V9 Q  `0 ^9 l; ^- V
& y$ C+ `1 K% ~0 s! s' k  @* T& b程序如下：
1 X1 V2 O' u$ u7 C; @2 L- g: z9 Exdata = [m(:,2) m(:,3) m(:,4)];
! _" [. `5 {6 q$ s0 a% Xydata = [0 0 0 0 0 0 0 0 0 0 0 0]';
k=[0 0 0 0 0 0 0 0 0 0];
( Q* k* i/ R7 g7 {! k/ o- P+ AF = @(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);
( @3 a' \5 n) g0 m8 B5 }k0=ones(1,10);
[k,resnorm]=lsqcurvefit(F,k0,xdata,ydata);

最后结果趋近于零，为什么？求指教
# r% _1 p% t  q' @" T/ o

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
B: -8.16614311365315E-6
) j# q5 y: b: G' E. AC: -8.62656710048033E-6
D: -9.30535264389442E-7
E: 1.3203421356307E-6
5 q- g4 e* B: x" a; hF: -2.39934410974343E-6
2 }1 _7 z- K0 r6 o" ]G: 0.00165699153341434
7 D/ [% K' h$ ]/ cH: 0.0022664316097256
1 A9 L& W+ ^7 l" I2 dK: -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
/ f5 w$ I) q# y4 i8 i/ IA: -9.52309128703986E-6
B: -8.166143 ...

您的意思是我的代码没有问题？是方程模型错了？
- j$ A# V- l9 m" y