机器人的数学基础齐次变换矩阵及其运算

机器人的数学基础齐次变换矩阵及其运算齐次变换矩阵及其运算由于各种原因，变换矩阵应写成方型形式，3*3或4*4均可.为保证所表示的矩阵为方阵，如果在同一矩阵中既表示姿态又表示位置，那么可在矩阵中加入比例因子使之成为4*4矩阵。第2页,共35页，星期六，2024年，5月变换可定义为空间的一个运动。已知一直角坐标系中的某点坐标，那么该点在另一直角坐标系中的坐标可通过齐次坐标变换来求得。变换可分为如下形式：纯平移纯旋转平移与旋转的结合第3页,共35页，星期六，2024年，5月1.平移的齐次变换空间某一点在直角坐标系中的平移,由A(x,y,z)平移至A′(x′,y′,z′),即

a′=Trans(Δx,Δy,Δz)a

平移算子第4页,共35页，星期六，2024年，5月①算子左乘:表示点的平移是相对固定坐标系进行的坐标变换。②算子右乘:表示点的平移是相对动坐标系进行的坐标变换。③该公式亦适用于坐标系的平移变换、物体的平移变换,如机器人手部的平移变换。第5页,共35页，星期六，2024年，5月例动坐标系{A}相对于固定坐标系的X0、Y0、Z0轴作(-1,2,2)平移后到{A’}；动坐标系{A}相对于自身坐标系(即动系)的X、Y、Z轴分别作(-1,2,2)平移后到{A’’}。已知A,写出坐标系{A’}、{A’’}第6页,共35页，星期六，2024年，5月2．旋转的齐次变换点在空间直角坐标系中的旋转如图所示。A(x,y,z)绕Z轴旋转θ角后至A’(x’,y’,z’),则A与A’之间的关系为：记为:a′=Rot(z,θ)a

旋转算子第7页,共35页，星期六，2024年，5月同理，绕x轴、Y轴旋转算子内容为：绕Z轴旋转算子内容为：第8页,共35页，星期六，2024年，5月如图所示单操作手臂，并且手腕也具有一个旋转自由度。已知手部的起始位姿矩阵为G1.若手臂绕Z0轴旋转90°，则手臂到达G2；若手臂不动，仅手部绕手腕Z1轴转90°，则手部到达G3.写出手部坐标系G2、G3表达式。第9页,共35页，星期六，2024年，5月第10页,共35页，星期六，2024年，5月3．复合齐次变换复合变换是由固定参考坐标系或当前运动坐标系的一系列沿轴平移和绕轴旋转变换所组成的。任何变换都可以分解为按一定顺序的一组平移和旋转变换。相对于固定坐标系相对于动坐标系算子左乘算子右乘第11页,共35页，星期六，2024年，5月已知坐标系中点U的位置矢量，将此点绕Z轴旋转90°，再绕Y轴旋转90°，如图所示，求旋转变换后所得的点W。第12页,共35页，星期六，2024年，5月平移变换和旋转变换可以组合在一个齐次变换中。上例中点U若还要作4i-3j+7k的平移，则只要左乘上平移变换算子即可得到最后的列阵表达式。第13页,共35页，星期六，2024年，5月第14页,共35页，星期六，2024年，5月

齐次变换矩阵的数学意义：

（1）同一点在不同坐标系{B}和{A}中的变换；（2）描述坐标系{B}相对于坐标系{A}的位置和方位；（3）点的运动算子。第15页,共35页，星期六，2024年，5月4．变换矩阵相乘对于给定的坐标系｛A｝、｛B｝、｛C｝，已知｛B｝相对｛A｝的描述为，｛C｝相对｛B｝的描述为，则。从而定义复合变换表示｛C｝相对于｛A｝的描述，是两变换矩阵的乘积。注意：变换矩阵相乘不满足“交换律”，变换矩阵的左乘和右乘的运动解释不同。第16页,共35页，星期六，2024年，5月复合变换可解释为：(1)和分别代表同一坐标系{C}相对于{A}和{B}的描述。则表示坐标系{C}从映射为的变换。

(2)坐标系{C}相对于{A}的描述是这样得到的：最初{C}与{A}重合，首先相对于{A}作运动，到达{B},然后相对{B}作运动，到达最终位置{C}。第17页,共35页，星期六，2024年，5月5.变换矩阵求逆如果知道坐标系｛B｝相对于｛A｝的描述。希望得到｛A｝相对于｛B｝的描述，这是个齐次变换求逆问题。对4*4矩阵直接求逆；利用齐次变换矩阵的特点，简化矩阵求逆运算。求逆问题可以描述为：已知，求解。利用旋转矩阵正交性利用复合变换公式(2.13),求出在{B}中描述。第18页,共35页，星期六，2024年，5月第19页,共35页，星期六，2024年，5月下面我们写出变换矩阵的一般表达形式

nxoxaxpxnyoyaypyT=nzozazpz0001式中n,o,a

是旋转变换列向量，p

是平移向量，其逆是

nxnynz-p.noxoyoz-p.oT-1=axay

az-p.a

0001式中的“.”表示向量的点积。第20页,共35页，星期六，2024年，5月计算T矩阵的逆矩阵。-0.5第21页,共35页，星期六，2024年，5月6变换方程为了描述机器人的操作，必须建立机器人本身各连杆之间，机器人与周围环境之间的运算关系。为此要规定各种坐标系来描述机器人与环境的相对位姿关系。{B}代表基座坐标系；{W}代表腕部坐标系；{T}代表工具坐标系；{S}代表工作站坐标系；{G}代表目标坐标系；它们之间的位姿关系用相应的齐次变换来描述。描述工作站坐标系相对于基座的位姿；描述目标坐标系相对于{S}的位姿；描述腕部{W}相对于基座{B}的位姿；………………第22页,共35页，星期六，2024年，5月对物体进行操作时，工具坐标系{T}相对于目标坐标系{G}的位姿直接影响操作效果。它是机器人控制和规划的目标。 实际上，它与其他变换之间的关系类似于空间尺寸链，则是封闭环。如图所示，工具坐标系{T}相对于基座坐标系{B}的描述可用两种变换矩阵的乘积来表示：令上面两式相等，则得变换方程第23页,共35页，星期六，2024年，5月 变换方程中的任一变换矩阵都可用其余的变换矩阵来表示。例如，为了对目标物进行有效操作，工具坐标系{T}相对于目标坐标系{G}的位姿是预先规定的，需要改变以达到这一目的，即通常规定，求。根据变换方程，可以立即求出第24页,共35页，星期六，2024年，5月旋转变换通式令是过原点的单位矢量，求绕k旋转θ角的旋转矩阵R(k,θ)。问题描述：令即R(k,θ)表示坐标系{B}相对于参考系{A}的方位。

坐标系{B}由坐标系{A}绕轴旋转角得到。k{A}第25页,共35页，星期六，2024年，5月xAyAzAxByBzB旋转变换通式再定义两坐标系{A’}和{B’}，分别与{A}和{B}固接，但要求(1){A’}和{B’}的z轴与k重合。(2)旋转之前{A’}和{B’}重合，{A}和{B}也重合。第26页,共35页，星期六，2024年，5月又因为所以可以得到：坐标系{B}绕k轴相对于{A}旋转θ角相当于：坐标系{B’}相对于{A’}的z轴旋转θ角，保持其他关系不变。则xAyAzAxByBzB坐标系{A}经过如下变换到坐标系{B}:第27页,共35页，星期六，2024年，5月把上式右端三矩阵相乘，并运用旋转矩阵的正交性质：第28页,共35页，星期六，2024年，5月①该式为一般旋转齐次变换通式，概括了绕X、Y、Z轴进行旋转变换的情况。反之，当给出某个旋转齐次变换矩阵，则可求得k及转角θ。②变换算子公式不仅适用于点的旋转，也适用于矢量、坐标系、物体的旋转。③左乘是相对固定坐标系的变换；右乘是相对动坐标系的变换。当kx=1，ky=kz=0时…………当ky=1，kx=kz=0时…………当kz=1，kx=ky=0时…………第29页,共35页，星期六，2024年，5月反之，若给出某个旋转齐次矩阵则可根据求出其等效矢量k及等效转角θ第30页,共35页，星期六，2024年，5月等效转轴和等效转角

给定旋转矩阵，求对应的等效旋转轴和等效转角设，令第31页,共35页，星期六，2024年，5月得到：方程两边矩阵的非对角元素成对相减，得到：两边平方后相加，所以整理后得到：所以，第32页,共35页，星期六，2024年，5月所以：方程两边矩阵的非对角元素成对相减，整理得到：（1）多值性：和值并不唯一，一般选取。（2）病态情况：当很小时，转轴不能确定，需要其它方法。注意：第3