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

SallyDai2

写的MATLAB代码，最小二乘求距离出现错误，麻烦大神帮忙看看，非常感谢！

clc;
3 A" g, R& [0 F; U: O+ T+ O9 qclear all;
A=[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];
; v. M0 q* X4 x: hP=[1 2 3 4 5 6 7 9 11 81 2 9 32 21 5 6 7 9 11 81 2];

c=3*10^8;
X=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);

' }5 L/ c: _) n/ u6 L    fxi=(A(i,1)-x(1,1))/Ri;
    fyi=(A(i,2)-x(2,1))/Ri;
    fzi=(A(i,3)-x(3,1))/Ri;
/ w8 {) i) \" B" Z' ~9 L7 J2 n; b    F(i,: )=[fxi fyi fzi -c];
& |* v6 F  v4 ^* C5 B0 j6 W- Y
2 M: i4 U/ {# p9 V, F" v    syms x0 y0 z0 t0 real;
    X=[x0 y0 z0 t0]';
# ^" D3 \& M( l
9 Y$ ^) C, E: q2 v- R# x6 F    p(i,1)=F(i,: )*X;
' Y9 d* D0 I3 A6 ~% t; r) d
% Q1 X3 X$ s5 z    x(1,1)=x(1,1)+x0;
9 S7 k' V+ G3 m2 u    x(2,1)=x(2,1)+y0;
    x(3,1)=x(3,1)+z0;
    x(4,1)=x(4,1)+t0;
5 q1 S# n: I5 _9 `$ @3 Pend
7 A; W6 i* u) k1 M) e6 b* z) j2 v
& g( t' R" v' S3 \. K+ W' |3 Yfigure(1);
6 z" V( z! W! |9 C& `. r) iplot(x0,'r-');
hold on;
% v5 {2 @1 L" ~figure(2);
plot(y0,'g*');
hold on;
( y& E+ ?' J$ a; m; r. |  s1 m2 Ofigure(3);
7 ^) Z# L9 y" _# f$ t2 s# zplot(z0,'b-');
: F# K$ y$ E9 V3 phold on;
3 G  K) G' O! d, S* gxlabel('X');
xlabel('Y');
1 ^5 G2 |7 ]" ]legend('measurement');
+ Z8 p% Y/ l) Q9 m& W( C
1 r6 I1 q+ A$ R! D8 c" ]

编译出现错误：
( q4 k7 j3 B) s5 l. J# z+ h错误使用 symengine (line 58)
DOUBLE cannot convert the input expression into a double array.
If the input expression contains a symbolic variable, use VPA.

& a: Q, s+ T" d3 I! F, p出错 zxec (line 21)
7 [4 u: C+ k+ C3 Z; u7 A    p(i,1)=F(i,: )*X;
( @( j2 }) X0 t. M# \

wu68aq

你的代码有问题。