第3章微分运动和速度学习内容:1微分关系2坐标系以及关节的微分运动3雅克比矩阵的相关运算及其与速度之间的关系学习重点:雅克比矩阵的计算1微分关系的概念微分运动就是指机器人的微小运动,而微分关系是指微分运动与速度之间的关系。2微分关系的理论推导下面这幅图是具有两个自由度的简单机构。其中每个连杆都能独立旋转,表示第一个连杆相对于参考坐标系的旋转角度,表示第二个连杆相对于第一个连杆的旋转角度。§3.1微分关系12微分关系让我们计算一下B点的速度ABABVVV/根据物理学中的相关公式,可以得到212122121121221211)cos()cos(cos)sin()sin(sinllllllVVYXBB接下来让我们对B点的位置方程求微分)sin(sin)cos(cos2121121211llYllXBB方程两边对和求微分,可得到12212122121121221211)cos()cos(cos)sin()sin(sinddlllllldydxBB可以看到,微分方程与速度方程极为相似,只不过二者表达的物理含义不同,如果在微分方程的两边同时除以dt,则两方程就完全相同了。3微分方程的结构微分关系212122121121221211)cos()cos(cos)sin()sin(sinddlllllldydxBBB点的微分运动方程雅克比矩阵关节的微分运动3.6§3.2雅克比矩阵1研究雅克比矩阵的意义由式3.6可以看到,雅克比矩阵将单个关节的微分运动或速度转换为感兴趣点的微分运动或速度,也可以将单个关节的运动与整个机构的运动联系起来。2雅克比矩阵的计算由式3.6可以看到,由于角度是时变的,所以雅克比矩阵也是时变的。所以我们可以通过对位置方程中的时间变量求导的方法来计算雅克比矩阵假设有一组变量为的方程:则变量和函数间的微分关系可以表示为jxiYjjiijjiijixxfYxxxxfxfxfxfxfxfYYY或211121211121根据上述关系,我们可以建立机器人的关节微分运动和机器人手坐标系微分运动之间的关系雅克比矩阵654321ddddddzyxdzdydx雅克比机器人机器人手沿x,y,z轴的微分运动机器人手绕x,y,z轴的微分旋转关节的微分运动DDJ雅克比矩阵例3.1给定某一时刻的机器人雅克比矩阵,给定关节的微分运动,求机器人手坐标系的线位移微分运动和角位移微分运动。100000000100002000000010000101010002J2.0001.01.00DzyxdzdydxJDD2.01.001.01.002.0001.01.00100000002000100000000000010010101002解:由例题可知:刚体或坐标系的微分运动包含微分移动矢量和微分转动矢量。前者由沿三个坐标轴的微分移动组成;后者由绕三轴的微分转动组成。雅克比矩阵的构造:一、矢量方程微分法;二、位姿方程微分法。§3.3坐标系的微分运动假设坐标系相对于参考坐标系做微量运动。一种情况是可以不考虑产生微分运动的原因来观察坐标系的微分运动,另一种情况是通过引起微分运动的机构来考察该微分运动。前一种情况只研究坐标系的运动以及坐标系表示的变化(如图a所示)。后一种情况则将研究产生该运动的机构的微分运动以及它与坐标系运动的联系(如图b所示)。此时,机器人关节的微量运动会使得手坐标系也产生微量运动。假设有一个机器人要将两片工作焊接在一起,为了获得最好的焊接质量,要求机器人以恒速运动,也就是说要求指定的手坐标系的微分运动能表示按特定姿态的恒速运动。这就涉及到坐标系的微分运动,而该运动是由机器人产生的。因此,应计算每一时刻各关节的速度,以使得由机器人产生的总运动就等于坐标系的期望速度。坐标系微分运动可以分为如下三个运动:●微分平移●微分旋转●微分变换(平移与旋转)我们首先研究坐标系的微分运动,然后研究机器人机构的微分运动,最后建立两者之间的联系。1微分平移微分平移就是坐标系平移一个为分量,因此它可以用Trans(dx,dy,dz)来表示,其含义是坐标系沿3条坐标轴做了微小量的运动。2微分旋转微分旋转是坐标系的小量旋转,它通常用来描述,即坐标系轴转动角度。绕三轴的转动分别定义为因为转动很小,所以),(dkRotkˆd1cos(sinxxx用弧度)因此,表示绕x,y,z轴的微分旋转矩阵为:10000100001001),(10000100100001),(zzzzRotxxxxRot10000100010001),(yyyyRot注意:上述矩阵违反了每个向量长度为1的规定。例如。然而由于微分值很小,在计算时高阶微分可以忽略不计。并且这样的近似并不影响计算结果。所以我们采用这样的向量长度。1122x1000010001000110000100100001),(),(yyxxyyRotxxRot另外我们再来看看矩阵乘法的顺序,不同的顺序是否会产生不同的结果。10000101001xyxzxy100001010011000010010000110000100010001),(),(xyxyyxxxyyxxRotyyRot比较两式,如果忽略掉所有的高阶微分变换,上述两式的结果是相同的。因此乘法的顺序并不影响最终计算结果。也就是说在微分运动分析中满足交换律。1000010101),(),(),(),(zyxzxyxzyxyxzyzzzRotyyRotxxRotdkRot绕三条坐标轴的三个微分运动可以表示为1000010101xyxzyz02.0,05.0,1.0zyx习题:求绕三个坐标轴作微分旋转所产生的总微分变换1000011.005.001.0102.0005.002.011000010101),(xyxzyzkRot解:3坐标系的微分变换坐标系的微分变换是微分平移和微分旋转运动的合成。如果用T表示原始坐标系,并假定由于微分变换所引起的坐标系T的变化量用dT表示,则有:][),(),,(),(),,(TIdkRotdzdydxTransdTTdkRotdzdydxTransdTT或可令:]),(),,([][]][[][IdkRotdzdydxTransTdT我们称为微分算子,用它乘以一个坐标系将导致坐标系的变化进一步求得0000000100001000010000110000101011000100010001),(),,(dzxydyxzdxyzxyxzyzdzdydxIdkRotdzdydxTrans习题对如下的坐标系B,绕y轴做0.1弧度的微分转动,然后微分平移[0.1,0,0.2],求微分变换的结果10003010500110100B]][[][BdB提示:解:00008.01.00000004.001.00]][[][0,1.0,0,2.0,0,1.0BdBzyxdzdydx其中,dB矩阵表示坐标系B的变化,该矩阵的每个元素表示坐标系中相应元素的变化。如,本例中dB意味着该坐标系沿x轴移动了0.4个单位的微小量,沿y轴无运动,沿z轴移动了-0.8个单位的微小量。它也意味着坐标系的旋转使得向量没有改变,而在向量的分量上改变了0.1,在向量的分量上改变了-0.1。—微分变化的理解naoxoza10002.21.01050014.1011.0000008.01.00000004.001.0010003010500110100dBBBoriginalnew由此,我们可求上例中坐标系B运动后的位姿,如下3.6坐标系之间的微分变化前面介绍的微分算子是相对于固定参考坐标系来说的,同样的,我们可以定义另外一个微分算子,是相对于当前坐标系的,这样使得可以在该坐标系(当前)中计算同样的变换。由于是相对于当前坐标系的,必须用右乘该坐标系的。如下式所示:]][[][][]][[][]][[][]][[]][[111TTTTTTTTdTTTT因此,上式可以用来计算相对于本身坐标系的微分算子。将上式矩阵相乘并加以简化,得到的结果如下:T10001paaaapoooopnnnnTzyxzyxzyx0000000dzxydyxzdxyz00000001dzxydyxzdxyzTTTTTTTTTTTTdpadzdpodydpndxazoynxTTTTTT应注意,看上去如同矩阵,但所有元素都是相对于当前坐标系的,这些元素可从以上矩阵相乘的结果求得,结果归纳如下:T例对如下的坐标系B,绕y轴做0.1弧度的微分转动,然后微分平移[0.1,0,0.2],求微分变换的结果。10003010500110100B解:00008.01.00000004.001.001000301050011010000002.0001.000001.01.000]][[][0,1.0,0,2.0,0,1.0BdBzyxdzdydx举例说明如何求得相对于本身坐标系的微分算子,回忆下面的例题(上节课出现过):]8.004.0[]2.001.0[]103.0[]103.0[351001.00ˆˆˆ2.001.001.00]3510[]001[100010,,,,,,,,,,,,,,,,,,,,,,,,dpkjipdpaon现在求出相对于本身坐标系的微分算子:由给定的信息中可以得到以下向量,用来计算向量BdBB001.04.08.00][