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

SallyDai2

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

, }$ x  i) k" j3 cclc;
clear 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];
P=[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;
  M, l' @9 d& t- Y( JX=zeros(4,1);
x=zeros(4,1);
$ }3 _. ^7 I7 n+ j. _for i=1:1:21
% a  o% Z! t, W    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);

' i7 M; I4 b3 a# n    fxi=(A(i,1)-x(1,1))/Ri;
    fyi=(A(i,2)-x(2,1))/Ri;
    fzi=(A(i,3)-x(3,1))/Ri;
6 r' Q9 b. d: u    F(i,: )=[fxi fyi fzi -c];

    syms x0 y0 z0 t0 real;
; O1 G' f* g' O" M% X. w    X=[x0 y0 z0 t0]';

    p(i,1)=F(i,: )*X;
5 @+ e  q+ R/ a" Y  ^1 F. o
    x(1,1)=x(1,1)+x0;
    x(2,1)=x(2,1)+y0;
    x(3,1)=x(3,1)+z0;
    x(4,1)=x(4,1)+t0;
  F2 l$ d8 m) A) o6 v- i/ M8 hend

figure(1);
plot(x0,'r-');
; a& _! ?2 L9 F6 Rhold on;
. x. x: P9 c; k+ |' k6 Nfigure(2);
2 f# |7 D+ i& y! k- O" B5 d# [" F' H. `' Mplot(y0,'g*');
hold on;
figure(3);
1 A8 x$ ^0 A% G/ ?% fplot(z0,'b-');
hold on;
' h/ J7 A3 l' Q6 Y: `0 Pxlabel('X');
& i" X- J( D2 ~& e5 b5 wxlabel('Y');
+ z! y9 j* E/ o, dlegend('measurement');
; Q; A" @' z# Y% h9 h7 d
9 u4 v: T) y7 D0 V* Y- v, Y4 C
- ?& G! U! r) K- v5 T6 r; G
, D& v* b% C. }8 U编译出现错误：
$ c! g  H7 s) \5 R3 h7 l错误使用 symengine (line 58)
DOUBLE cannot convert the input expression into a double array.
If the input expression contains a symbolic variable, use VPA.
" a, o5 _7 K/ k2 q$ f! R0 ?
3 n  `0 _4 I, b4 Z出错 zxec (line 21)
    p(i,1)=F(i,: )*X;
6 U; k; X6 i3 c6 c$ C6 F! V

wu68aq

你的代码有问题。