找回密码
 注册
关于网站域名变更的通知
查看: 392|回复: 1
打印 上一主题 下一主题

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

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2020-10-23 19:21 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

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

该用户从未签到

2#
发表于 2020-10-23 19:35 | 只看该作者
你的代码有问题。
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-7-30 01:11 , Processed in 0.125000 second(s), 24 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表