写的MATLAB代码，最小二乘求距离出现错误？

SallyDai2

写的MATLAB代码，最小二乘求距离出现错误，麻烦大神帮忙看看，非常感谢！
  @" k: C" [; i* F! }% q1 o
. U, _: p) J, n0 f' M- c3 zclc;
" i6 ?  @" C2 Mclear all;
/ \) U1 [3 x# K6 B8 |5 N# lA=[1 2 3;3 4 5;6 7 8;6 7 8;6 7 8;11 12 14;22 23 43;1 2 3;3 4 5;6 7 8;6 7 8;6 7 8;11 12 14;22 23 43;1 2 3;3 4 5;6 7 8;6 7 8;6 7 8;11 12 14;22 23 43];
/ A$ Z& q  y8 lP=[1 2 3 4 5 6 7 9 11 81 2 9 32 21 5 6 7 9 11 81 2];

" }: \) Y" T5 N4 o5 {, B0 ac=3*10^8;
# {1 v4 f/ p: U7 @! z; UX=zeros(4,1);
x=zeros(4,1);
for i=1:1:21
    Ri=sqrt((A(i,1)-x(1,1)).^2+(A(i,2)-x(2,1)).^2+(A(i,3)-x(3,1)).^2);
    p(i,1)=Ri+x(4,1)-P(i,1);

    fxi=(A(i,1)-x(1,1))/Ri;
8 j4 X6 k7 Y/ T+ P+ x    fyi=(A(i,2)-x(2,1))/Ri;
: @. f% D$ T4 J4 ]' `  V    fzi=(A(i,3)-x(3,1))/Ri;
    F(i,: )=[fxi fyi fzi -c];

    syms x0 y0 z0 t0 real;
    X=[x0 y0 z0 t0]';
8 }  t2 g" }, ]  u# s8 w: f
. _# [! i6 B/ ]    p(i,1)=F(i,: )*X;
' Z# \3 h9 u5 H5 N4 s& N) g
) O3 x5 K$ Z& N0 i2 u# e    x(1,1)=x(1,1)+x0;
: j0 w$ y7 k  f- ^6 z* g    x(2,1)=x(2,1)+y0;
* |1 [: {4 g. ?4 e3 O% D7 ^    x(3,1)=x(3,1)+z0;
    x(4,1)=x(4,1)+t0;
end
& x. B( Z* M- P7 q
+ D- [/ r5 ?- m/ H4 b1 a! J! xfigure(1);
plot(x0,'r-');
: w" u* q0 f: k6 J8 n2 Ahold on;
figure(2);
plot(y0,'g*');
hold on;
  f; h! L6 K) w! x) h; b+ Nfigure(3);
9 x' g- q0 C! Uplot(z0,'b-');
4 D4 I8 H  C0 N+ c$ lhold on;
( a+ x6 ~9 p0 J# r, j0 x/ qxlabel('X');
xlabel('Y');
. o9 \1 @$ e( g8 @2 u" L$ W# tlegend('measurement');


6 V4 [1 [8 m! _
编译出现错误：
& ?# y5 q8 c' N; g, V& K) w错误使用 symengine (line 58)
7 q# T5 O( v/ X; @DOUBLE cannot convert the input expression into a double array.
If the input expression contains a symbolic variable, use VPA.
) P, c: r5 e1 K4 I- z
/ y! ]) T8 b+ n1 f6 B7 j出错 zxec (line 21)
, ^& M) Y: `( ~$ i  s/ N4 q$ E    p(i,1)=F(i,: )*X;
6 z% `8 s! t- w/ B

wu68aq

你的代码有问题。