|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
写的MATLAB代码,最小二乘求距离出现错误,麻烦大神帮忙看看,非常感谢!2 W8 U; r' m! ?6 m! K3 W$ s* I
( Q3 q( Y8 X7 Zclc;
V. Y% B1 s2 v$ F$ kclear all;' C4 H& y$ M/ P; }
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];. ^) m4 M$ D5 [5 ?* R' X
P=[1 2 3 4 5 6 7 9 11 81 2 9 32 21 5 6 7 9 11 81 2];
8 x+ }/ R6 ~- {' G; c$ [! S) @( `0 ]$ p" {$ l
c=3*10^8;
* G2 G, h$ B2 a8 f# x( [X=zeros(4,1);
3 z" H- f/ s7 ^2 }+ } \x=zeros(4,1);2 x& Q! Q) s7 e8 N [
for i=1:1:21
. q+ l; [' e. p' t6 W Ri=sqrt((A(i,1)-x(1,1)).^2+(A(i,2)-x(2,1)).^2+(A(i,3)-x(3,1)).^2);# p6 E/ G9 G B4 l4 x+ n
p(i,1)=Ri+x(4,1)-P(i,1);
2 M" j( k, _) ?9 [. U2 A
7 ~0 x$ J) N5 ~1 |& a fxi=(A(i,1)-x(1,1))/Ri;
9 v3 f/ f0 i% h8 ^$ J fyi=(A(i,2)-x(2,1))/Ri;
; l' i! ]( J) |$ ~" q fzi=(A(i,3)-x(3,1))/Ri;& N% ?" s5 C8 }# Q4 ]
F(i,: )=[fxi fyi fzi -c];
/ O" |" r6 ?) c# G4 L( E6 b: [3 e& M" g5 J \
syms x0 y0 z0 t0 real;
( k+ D8 R4 X1 S) c4 f3 ? X=[x0 y0 z0 t0]';: ^; d( V, ~: z k5 x
# w) t2 b; U: }. l2 F1 ?' Y p(i,1)=F(i,: )*X;% ~ j9 \4 _) a" {) }
E1 t) m4 ]" x9 h) \ x(1,1)=x(1,1)+x0;, a0 v0 J" L, L7 d
x(2,1)=x(2,1)+y0;* G" K1 S) M; ~* U& Z
x(3,1)=x(3,1)+z0;
+ Z" k; |. _7 O: f! H x(4,1)=x(4,1)+t0;% v* c. a8 O' [% v* ?) @: V( S
end
5 i' g V5 b8 i' ]8 A/ |! G6 P+ p9 F# Q" P7 ^, E1 `
figure(1);
0 w0 ]* _( r6 A# uplot(x0,'r-');
. }/ o4 @3 A% ^; }* z; V! l3 Uhold on;
$ [# V% c: ~# A5 F/ b7 ofigure(2);- ]+ q: U B: F; }; {
plot(y0,'g*');
' m, t0 T' X& u8 [7 Rhold on;7 x! @9 c2 W. i, D+ S2 |8 h- x9 p
figure(3);; w+ [* {/ M/ l! J. B8 G
plot(z0,'b-');* Z- Q0 P! D, [( j: x( X" H( ^
hold on;
- c$ R& j* m6 Mxlabel('X');0 z9 s: d; L L# X8 O# K6 y
xlabel('Y');
" ^0 z m4 R8 @; m' Dlegend('measurement');
1 y% R2 _ H5 K8 o
4 n& K6 Y1 X4 N. `: u7 c6 r* d4 P P* G8 }( W8 p
% [' w% ?! y: r* O
编译出现错误:3 N+ o: o9 p1 O) A
错误使用 symengine (line 58)# v1 z+ |3 X1 N2 k
DOUBLE cannot convert the input expression into a double array.6 B% S( K; s0 ]1 Z
If the input expression contains a symbolic variable, use VPA.
' P4 l! t: ], \, t2 z! O- W, J
5 w* ~- o$ N) h" \8 E0 ]; L+ \出错 zxec (line 21) A: g' `; V; T, s% h8 G- L
p(i,1)=F(i,: )*X; l8 g, C, P4 O, O( C
|
|