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

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

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
写的MATLAB代码,最小二乘求距离出现错误,麻烦大神帮忙看看,非常感谢!  ^. m! O, [6 g. @
3 ^& c7 D: t) X1 X  X
clc;
& V2 W( p* o0 p0 [( Zclear all;' E' E' c3 q$ H$ z( p. y2 G% x
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];
* n, v" _" I3 ^P=[1 2 3 4 5 6 7 9 11 81 2 9 32 21 5 6 7 9 11 81 2];. R  u% i7 b) D$ w

  C, O0 c& J+ \4 {) S" ?5 `c=3*10^8;
8 a+ H6 J8 l8 V  J* [2 fX=zeros(4,1);0 `5 x% P; t" E' f" N
x=zeros(4,1);
, A5 `7 v: r: R/ Zfor i=1:1:21
: m+ a5 k& Q$ N' e. _* V) h. p/ C    Ri=sqrt((A(i,1)-x(1,1)).^2+(A(i,2)-x(2,1)).^2+(A(i,3)-x(3,1)).^2);1 _3 A4 q( }& [- O1 V: c
    p(i,1)=Ri+x(4,1)-P(i,1);
, r( ]  _# g+ x& ]3 O/ F2 v! z: j8 G6 f) g
    fxi=(A(i,1)-x(1,1))/Ri;! S1 Z' Q4 d& z# m; A
    fyi=(A(i,2)-x(2,1))/Ri;
- T( L5 j) `- d: ~. C/ c" P2 ~2 C    fzi=(A(i,3)-x(3,1))/Ri;! [/ M5 x% @4 I. Q& M: n5 N" [: y# ^
    F(i,: )=[fxi fyi fzi -c];
$ u3 V/ @2 W1 ?0 f( D* j
! T( A" v/ s4 d- g0 C    syms x0 y0 z0 t0 real;( A* `' X. W, k1 K. ]5 K; j
    X=[x0 y0 z0 t0]';9 k) M& Z# g. R) h) w

+ X% m( J1 g# |# d1 A/ t0 f* F    p(i,1)=F(i,: )*X;' `, K. k. j& Q1 b" h) z. ?
! E. h1 O4 @: c+ t4 ^1 |7 p6 I
    x(1,1)=x(1,1)+x0;$ V8 b9 @! c" ^( X8 q8 e% C
    x(2,1)=x(2,1)+y0;, J3 g2 N  J& v# W* Q
    x(3,1)=x(3,1)+z0;
" ]0 C2 M$ e! C/ x* u5 W* j, |    x(4,1)=x(4,1)+t0;
, w4 I  @) I" [* F# fend
/ E: v5 D$ s" w% A
4 z1 p/ O2 `* N% ?+ Y& C$ Jfigure(1);
* c" V3 R# ~& ~! y) ^; Z( wplot(x0,'r-');+ |) s1 {4 D" b- `, T
hold on;3 m# y! ~( a5 s- v4 }
figure(2);
% \0 F$ p2 ?& O" d! ^8 G* bplot(y0,'g*');
# b4 S0 F8 [- Q6 J# _4 bhold on;( G4 V  a7 o9 ~
figure(3);
3 V% [  m4 \! O3 jplot(z0,'b-');% E$ O7 ~, G. b2 }9 E& z% l
hold on;
0 _, q6 q1 ?6 ]# {( K$ f/ Oxlabel('X');; g2 k) s3 K6 G( @. r' V5 q  G
xlabel('Y');
" E7 x6 o1 v! Rlegend('measurement');, ?+ j1 X1 W* k$ p8 X+ b

% u$ D* c1 e! ^1 }* A
8 r) v# ?' c# u  s  m% z: L, b% ^0 p7 i. Z  K
编译出现错误:
! N( ^+ v! z1 |; r/ g错误使用 symengine (line 58)1 ]' Y% Z3 S& n2 g
DOUBLE cannot convert the input expression into a double array.' z8 z4 g; h8 s% {' |' \0 R
If the input expression contains a symbolic variable, use VPA.
! \: ~8 H) S- w* T  A+ L- Y) U. O" Z' P9 N, \& j6 J4 N! q
出错 zxec (line 21)
0 S, S  `" N1 {$ [    p(i,1)=F(i,: )*X;, s" Z) I; Y! h, I

该用户从未签到

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

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 10:19 , Processed in 0.156250 second(s), 23 queries , Gzip On.

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

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

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