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

机器人在关节空间圆弧的轨迹规划仿真出现错误

[复制链接]

该用户从未签到

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

EDA365欢迎您登录!

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

x
机器人在关节空间圆弧的轨迹规划仿真出现错误,仿真结果显示
* M; B: M( A+ I( j! s7 t- P??? Error using ==> plot
$ L7 X2 k) ], B4 p# `6 c& g( ^( {Conversion to double from cell is not possible.- S6 i* T8 O5 N+ J
Error in ==> sy2 at 29- E- ?. D3 c( m) F3 ^2 y$ M, C/ C3 r) C
plot(t{i},q(:,1))
3 l4 o* M/ b# c6 U" S1 ?麻烦帮忙解答,谢谢!7 R2 w3 `5 c  N6 \7 Q8 @' x

2 j8 \) S( q7 L& ^% g6 h6 l% ^程序如下:
! ]# h# }3 Q* q. eL{1}=link([0 0 -pi/2 0 0],'mod');
" N- R& Z* Z+ U" X0 Z; N3 f9 S( D  j0 CL{2}=link([pi/2 0 pi/2 0 0],'mod');
3 J. p8 Z& w( fL{3}=link([0 1150 0 0 0],'mod');9 M+ A5 a( W8 j: V5 c
L{4}=link([pi/2 0 0 1466 0],'mod');
% K$ j1 N! C9 M$ YL{5}=link([pi/2 -188 0 0 0],'mod');& v4 @8 `! G# Y
hithd=robot(L,'hithdrobot');5 B! o- L; Z# ~7 k( E+ y
clear L
5 x  h. N1 z, j! B% talpha=0:pi/80:2*pi;%定义圆弧角
3 R+ f8 r1 K9 {" iR=30;%给定半径
# {( i  T7 o: g" m: q4 ta1=2022.5+R*sin(alpha);%定义x变量
3 |$ ~' t$ j9 C+ Z, i) p) v9 c- ~b1=30-R*cos(alpha);%定义y变量* H" y3 W- N: h8 X. U. r
for i=1:1:41
$ E! C: |; u6 B' ]T{i}=[0 1 0 a1(:,i);0 0 1 b1(:,i);1 0 0 -230;0 0 0 1];%定义半圆弧上41个点的齐次变换矩阵" b0 T: X: y' O! ~% X* c# G& y- t$ O
end& `; N( U* g% k
for i=1:1:40
$ M! F6 N; `2 {2 v2 E   t{i}=[i:.025:i+1]';* l) G6 Y0 \. ~9 `
Ta{i}=ctraj(T{i},T{i+1},length(t{i}));%生成40个线段的笛卡尔轨迹
/ F$ ^. T4 n; Hend
! b+ B, }2 j& ~9 c- N$ ^for i=1:1:407 I9 V2 k# g- R. j' `  Z3 t
tic;0 R3 V% E8 F+ d9 M0 t5 d7 G$ Q
q{i}=ikine(hithd,Ta{i});%运动学逆解40个线段的关节坐标
: R! q: Z; @+ `% \toc
7 ?+ P  }6 s3 X! ?: iend  D: }. A" @) C' p$ U
subplot(3,1,1)/ O8 X' E9 ~0 ]  S
plot(t{i},q(:,1))
" A# D; G+ E* i" e* e6 pxlabel('Time(s)');
* ^/ {5 u# p3 i. i4 o9 s9 Qylabel('Joint 1(rad)'): X  ~9 O$ w; |4 A2 Q& ^4 Z) l  ~; ]
subplot(3,1,2)$ L, H+ ~' k: s: I( v6 D( t; [
plot(t{i},q(:,2))
& ?9 v- q. ]7 H, M% m) g) p. h- Dxlabel('Time(s)');8 K* J9 X4 f: L+ a
ylabel('Joint 2(rad)')
1 w% H2 e6 E% ?; [8 |" A4 j0 Rsubplot(3,1,3). x5 p; }. ?$ j4 i1 M
plot(t{i},q(:,3))2 a& [! f: a$ s6 M1 {
xlabel('Time(s)');
2 q. a9 O* q' Dylabel('Joint 3(rad)')5 A0 R9 B9 [9 a3 r  v9 [
hold on' ^6 A8 I/ ?4 B! g# N

该用户从未签到

2#
发表于 2020-11-18 17:13 | 只看该作者
帮你顶一下

该用户从未签到

4#
发表于 2020-11-19 14:22 | 只看该作者
楼主解决了吗?我也遇到这个问题了
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

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

EDA365公众号

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

GMT+8, 2025-11-24 05:08 , Processed in 0.125000 second(s), 23 queries , Gzip On.

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

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

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