机器人学数学基础

机器人学数学基础杜志江科学园C1栋机器人研究所512室86414462-25duzj01@研究方向：医疗机器人智能移动机器人参考教材：付京逊《机器人学》蔡自兴《机器人学》理查德·鲍尔《机器人操作手·数学·编程与控制》第2页,共44页，星期六，2024年，5月机器人运动学第3页,共44页，星期六，2024年，5月引言

机器人位置和姿态的描述机器人可以用一个开环关节链来建模由数个驱动器驱动的转动或移动关节串联而成一端固定在基座上，另一端是自由的，安装工具，用以操纵物体人们感兴趣的是操作机末端执行器相对于固定参考坐标数的空间几何描述，也就是机器人的运动学问题机器人的运动学即是研究机器人手臂末端执行器位置和姿态与关节变量空间之间的关系第4页,共44页，星期六，2024年，5月运动学研究的问题Whereismyhand?DirectKinematicsHERE!HowdoIputmyhandhere?InverseKinematics:Choosetheseangles!运动学正问题运动学逆问题第5页,共44页，星期六，2024年，5月机器人示例第6页,共44页，星期六，2024年，5月动画示例第7页,共44页，星期六，2024年，5月

哈佛大学RogerBrockett建立的指数积公式运动学滚动接触非完整控制数学基础-刚体运动参考文献：机器人操作的数学导论作者：理查德·摩雷李泽湘夏卡恩·萨斯特里翻译：徐卫良钱瑞明（东南大学）研究运动学的方法第8页,共44页，星期六，2024年，5月丹纳维特（Denavit）和哈顿伯格（Hartenberg）于1955年提出了一种矩阵代数方法解决机器人的运动学问题—D-H方法具有直观的几何意义能表达动力学、计算机视觉和比例变换问题其数学基础即是齐次变换第9页,共44页，星期六，2024年，5月第二章数学基础—齐次坐标和齐次变换第10页,共44页，星期六，2024年，5月2.1点和面的齐次坐标

2.1.1点的齐次坐标

一般来说，n维空间的齐次坐标表示是一个（n+1）维空间实体。有一个特定的投影附加于n维空间，也可以把它看作一个附加于每个矢量的特定坐标—比例系数

式中i,j,k为x,y,z轴上的单位矢量，a=,b=,c=，w为比例系数

显然，齐次坐标表达并不是唯一的，随w值的不同而不同。在计算机图学中，w作为通用比例因子，它可取任意正值，但在机器人的运动分析中，总是取w=1

。列矩阵第11页,共44页，星期六，2024年，5月[例]：可以表示为：

V=[3451]T

或V=[68102]T

或V=[-12-16-20-4]T

第12页,共44页，星期六，2024年，5月

齐次坐标与三维直角坐标的区别V点在ΣOXYZ坐标系中表示是唯一的（x、y、z）而在齐次坐标中表示可以是多值的。不同的表示方法代表的V点在空间位置上不变。

第13页,共44页，星期六，2024年，5月

几个特定意义的齐次坐标：[000n]T—坐标原点矢量的齐次坐标，n为任意非零比例系数[1000]T—

指向无穷远处的OX轴[0100]T—

指向无穷远处的OY轴[0010]T—

指向无穷远处的OZ轴

2个常用的公式：第14页,共44页，星期六，2024年，5月2.1.2平面的齐次坐标平面齐次坐标由行矩阵P=[abcd]来表示当点v=[xyzw]T处于平面P内时，矩阵乘积PV=O，或记为

如果定义一个常数m=，则有：=可以把矢量解释为某个平面的外法线，此平面沿着法线方向与坐标原点的距离为。=第15页,共44页，星期六，2024年，5月[例]：一个平行于x、y轴，且在z轴上的坐标为单位距离的平面P可以表示为：或

有：PV=例如：点V=[102011]T

必定处于此平面内，而点V=[0021]T处于平P的上方点V=[0001]T处于P平面下方。因为：

与点矢相仿，平面也没有意义第16页,共44页，星期六，2024年，5月2.2旋转矩阵及旋转齐次变换

2.2.1旋转矩阵

设固定参考坐标系直角坐标为ΣOxyz，动坐标系为ΣO´uvw，研究旋转变换情况。①

初始位置时，动静坐标系重合，O、O´重合，如图。各轴对应重合，设P点是动坐标系ΣO´uvw中的一点，且固定不变。则P点在ΣO´uvw中可表示为：

、、为坐标系ΣO´uvw的单位矢量，则P点在Σoxyz中可表示为：第17页,共44页，星期六，2024年，5月当动坐标系ΣO´uvw绕O点回转时，求P点在固定坐标系Σoxyz中的位置已知：P点在ΣO´uvw中是不变的仍然成立，由于ΣO´uvw回转，则：用矩阵表示为:（2-7）第18页,共44页，星期六，2024年，5月反过来：2.2.2旋转齐次变换用齐次坐标变换来表示式（2-7）第19页,共44页，星期六，2024年，5月2.2.3三个基本旋转矩阵和合成旋转矩阵

三个基本旋转矩阵

即动坐标系求的旋转矩阵，也就是求出坐标系中各轴单位矢量在固定坐标系中各轴的投影分量，很容易得到在重合时，有：第20页,共44页，星期六，2024年，5月由图2-5可知，在y轴上的投影为，在z轴上的投影为,在y轴上的投影为，在z轴上的投影为，所以有：方向余弦阵第21页,共44页，星期六，2024年，5月同理：

三个基本旋转矩阵:

第22页,共44页，星期六，2024年，5月

合成旋转矩阵:例1：在动坐标中有一固定点，相对固定参考坐标系做如下运动：①R（x,90°）；②R(z,90°)；③R(y,90°)。求点在固定参考坐标系下的位置。解1：用画图的简单方法第23页,共44页，星期六，2024年，5月解2：用分步计算的方法①R（x,90°）②R（z,90°）③R（y,90°）（2-14）（2-15）（2-16）第24页,共44页，星期六，2024年，5月

上述计算方法非常繁琐，可以通过一系列计算得到上述结果。将式（2-14）（2-15）（2-16）联写为如下形式：R4x4为二者之间的关系矩阵，我们令：定义1：当动坐标系绕固定坐标系各坐标轴顺序有限次转动时，其合成旋转矩阵为各基本旋转矩阵依旋转顺序左乘。注意：旋转矩阵间不可以交换

第25页,共44页，星期六，2024年，5月

平移齐次变换矩阵注意：平移矩阵间可以交换，平移和旋转矩阵间不可以交换第26页,共44页，星期六，2024年，5月2.2.4相对变换

举例说明：例1：动坐标系∑0′起始位置与固定参考坐标系∑0重合,动坐标系∑0′做如下运动：①R(Z,90º)②R（y,90º）③Trans(4，-3,7)，求合成矩阵解1：用画图的方法：第27页,共44页，星期六，2024年，5月解2：用计算的方法根据定义1，我们有：

以上均以固定坐标系多轴为变换基准，因此矩阵左乘。如果我们做如下变换，也可以得到相同的结果：例2：①先平移Trans(4,-3,7)；②绕当前轴转动90º；③绕当前轴转动90º；求合成旋转矩阵。（2-20）第28页,共44页，星期六，2024年，5月解1：用画图的方法解2：用计算的方法（2-21）第29页,共44页，星期六，2024年，5月式（2-20）和式（2-21）无论在形式上，还是在结果上都是一致的。因此我们有如下的结论：动坐标系在固定坐标系中的齐次变换有2种情况：定义1：如果所有的变换都是相对于固定坐标系中各坐标轴旋转或平移，则依次左乘，称为绝对变换。定义2：如果动坐标系相对于自身坐标系的当前坐标轴旋转或平移，则齐次变换为依次右乘，称为相对变换。结果均为为动坐标系在固定坐标中的位姿（位置+姿态）。相对于固定坐标系，

也就是说，动坐标系绕自身坐标轴做齐次变换，要达到绕固定坐标系相等的结果，就应该用相反的顺序。

第30页,共44页，星期六，2024年，5月右乘的意义：机器人用到相对变换的时候比较多例如机械手抓一个杯子，如右图所示，手爪需要转动一个角度才抓的牢，相对于固定坐标系表达太麻烦，可以直接根据手爪的坐标系表示但也要知道在∑O中的位姿，就用右乘的概念。

oH第31页,共44页，星期六，2024年，5月2.2.5绕通过原点的任意轴旋转的齐次变换有时动坐标系∑O´可能绕过原点O的而分量分别为rx、ry、rz的任意单位矢量r转动φ角。研究这种转动的好处是可用∑O´绕某轴r的一次转动代替绕∑O各坐标轴的数次转动为推导此旋转矩阵，可作下述变换：绕X轴转α角，使r轴处于XZ平面内绕Y轴转-β角，使r

轴与OZ轴重合绕OZ轴转动φ角绕Y轴转β角绕X轴转-α角第32页,共44页，星期六，2024年，5月由上图容易求出：由定义1和定义2，上述5次旋转的合成旋转矩阵为：（2-25）第33页,共44页，星期六，2024年，5月带入式（2-25），得由该式可以推出3个基本旋转矩阵第34页,共44页，星期六，2024年，5月2.2.6齐次变换矩阵的几何意义

设T=，有一个手爪，已知其在∑O的位置，设一个该坐标系∑O´，已知，，那么∑O´在∑O中的齐次坐标变换为，如果手爪转了一个角度，

则：第35页,共44页，星期六，2024年，5月T反映了∑O´在∑O中的位置和姿态，即表示了该坐标系原点和各坐标轴单位矢量在固定坐标系中的位置和姿态。该矩阵可以由4个子矩阵组成，写成如下形式：为姿态矩阵（旋转矩阵），表示动坐标系∑O´在固定参考坐标系∑O中的姿态，即表示∑O´各坐标轴单位矢量在∑O各轴上的投影为位置矢量矩阵，代表动坐标系∑O´坐标原点在固定参考坐标系∑O中的位置为透视变换矩阵，在视觉中进行图像计算，一般置为0为比例系数

第36页,共44页，星期六，2024年，5月如果需要求解∑O在∑O´中的位置和姿态，此时的齐次变换矩阵为，即求逆矩阵：

其中：这些式子以后经常遇到，在机器人计算中，所要求的就是齐次变换矩阵第37页,共44页，星期六，2024年，5月2.2.7透镜成像的齐次变换

第38页,共44页，星期六，2024年，5月因此，进行机器人运动学计算时，不能省略透视矩阵，有摄像头时，透视矩阵为

[0-0]，没有摄像头时为[000]。第39页,共44页，星期六，2024年，5月知识点：

点和面的齐次坐标和齐次变换三个基本旋转矩阵绝对变换：如果所有的变换都是相对于固定坐标系中各坐标轴