2021-08-05十四讲第三章

这章主要描述了何描述三维刚体运动,主要围绕两个主题展开:旋转和位移。在书中提到了有四种方法可以描述旋转和位移:旋转矩阵、矩阵、旋转向量、欧拉角、四元数。
旋转矩阵用矩阵的基坐标变换来描述旋转,加上一个向量,就可以完整地表示旋转和平移,但由于旋转矩阵的运算不是线性的,因此改进为由于旋转向量方向是旋转轴的方向,长度为旋转的角度,因此改进为变换矩阵。但这还没完,矩阵的表示法数据冗余度大,并且矩阵必须是正交阵,我们又提出了旋转向量,他的方向与旋转轴一致,长度和旋转角度相同,因此也可以表示旋转和位移。欧拉角比较直观,它一个旋转拆分成分别绕三个坐标轴旋转,可以用于判断和想象旋转是什么样,但是会出现万向锁问题。旋转向量和欧拉角这种三个自由度的表示法,会出现奇异性(例如万向锁),但矩阵表示又不够紧凑,而四元数则兼顾不带奇异性和紧凑型,用一个标量和一个向量组成的四元数字来描述旋转和平移,唯一的缺点可能是计算比较复杂。

slambook第三章.png

1旋转矩阵

1.1公式推导

机器人运动中存在两个坐标系,一个是惯性坐标系(或世界坐标系)p_w其基坐标为\begin{bmatrix}{e1},e2,e3\end{bmatrix}
一个是机器人的移动坐标系p_c,其及坐标为
\begin{bmatrix}e1^{\prime},e2^{\prime},e3^{\prime}\end{bmatrix}
对于同一个三维世界中的某个向量可以分别用这两个坐标系表示:
p= \begin{bmatrix}e_1,e_2,e_3\end{bmatrix}\begin{bmatrix}a_1\\a_2\\a_3\end{bmatrix}=\begin{bmatrix}e_1^{\prime},e_2^{\prime}{,e_3^{\prime}}\end{bmatrix}\begin{bmatrix}{a_1^{\prime}}\\a_2^{\prime}\\a_3^{\prime}\end{bmatrix}
两边同乘\begin{bmatrix}{e_1^{T}}\\{e_2^{T}}\\e_3^{T}\end{bmatrix},于是得到旋转矩阵R=\begin{bmatrix}{e_1^{T}}\\{e_2^{T}}\\e_3^{T}\end{bmatrix}\begin{bmatrix}e_1,e_2,e_3\end{bmatrix}
知道了旋转阵,也就描述了这个旋转,若已知\vec{e1},则可知\vec{e1^{'}},反之亦然。R^{-1}是R的逆变换。

1.2旋转矩阵定义

旋转矩阵等价于正交阵:
SO(n)={\{}{R{\in}R^{n{\times}n}}|{RR^T=I},{det(R)=1}{\}}
SO(n)也称为特殊正交群。

1.3齐次坐标和变换矩阵

由旋转矩阵R和一个平移向量t,可以完整表示一个欧式空间所有的变换:
a^{\prime}=Ra+t
假设我们有多次运算,比如b=R_1a+t_1,c=R_2b+t_2,那么a到c的变换是c=R_2(R_1a+t_1)+t_2
这样的形式在多次变换后十分复杂,于是引入了变换矩阵。

1.3.1齐次坐标

在一般向量末尾加上一维,规定为1,成为齐次坐标,一个齐次坐标数乘一个k,仍然表示同一个点,因此,多个齐次向量可以是相同的:[1,1,1,1]和[2,2,2,2]。但是我们习惯把最后一维强制变为1,于是齐次坐标表示唯一。

1.3.2变换矩阵

加入齐次坐标以后,
\begin{bmatrix}{a^{\prime}}\\1\end{bmatrix}=\begin{bmatrix}R,t\\0^{T}{,}1\end{bmatrix}\begin{bmatrix}a\\1\end{bmatrix}\triangleq{T}\begin{bmatrix}a\\1\end{bmatrix}
那么,多次变换就方便了:
\widetilde{b}=T_1\widetilde{a},\widetilde{c}=T_2\widetilde{b}\Rightarrow{T_2}{T_1}\widetilde{a}
定义:
SE(3)={\{T}=\begin{bmatrix}{R,t}\\{0,1}\end{bmatrix}\in{R^{4{\times}{4}}}|R\in{SO(3)},t{\in}{R^3}{\}}
SE(3)也被称作特殊欧式群。
类似的,变换矩阵的逆也是一个相反的变换。

2.旋转向量

旋转矩阵和变换矩阵用多余3维的数据表示3维的旋转数据,有数据冗余,并且旋转矩阵必须是正交阵,所以提出旋转向量,旋转向量的大小是旋转角度,方向是旋转轴的方向。

2.1旋转向量变为旋转矩阵R=cos{\theta}I+(1-cos{\theta})n{n^T}+sin{\theta}{n^{\land}}

2.2由矩阵矩阵得到旋转向量

两边取迹trace,得到旋转角\theta{\theta}={arccos(tr(R)-1)\over 2}
再求旋转向量\vec{n},由于旋转向量\vec{n}绕旋转向量不变,所以R\vec{n}=\vec{n}.
因此解出向量和旋转角。

3.欧拉角

欧拉角是提供了直观的方式考察旋转,把旋转分解到三个坐标轴上。有多种转角方式包括ZYX,XYZ等,字母的次序代表旋转轴的旋转次序。

3.1ZYX旋转

ZYX旋转(yaw-pitch-roll,偏航-俯仰-滚转),依次绕ZYX旋转。

3.2万向锁

当俯仰角为正负90度时,第一次和第三次旋转的旋转轴会相同因此出现万向锁。

4.四元数

4.1四元数定义

定义一
q=q_0+q_1i+q_2j+q_3k
并且i,j,k满足:
\begin{cases}{i^2}={j^2}={k^2}=-1\\{ij=k},{ji=-k}\\{jk=i},{kj=-i}\\{ki=j},{ik={-j}}\end{cases}
可以看到,按照ijk邻接的循环次序为正,否则为负。
定义二
q=[s,\vec{v}],s=q_0\in\mathbb{R},v=[q_1,q_2,q_3]^T\in\mathbb{R^3}
四元数与旋转向量的关系
旋转向量到四元数
q=[{cos\theta\over2},n_x{sin\theta\over2},n_y{sin\theta\over2},n_z{sin\theta\over2}]
四元数到旋转向量
\begin{cases}\theta=2arccosq_0\\\begin{bmatrix}n_x,n_y,n_z\end{bmatrix}^T={\begin{bmatrix}q_1,q_2,{q_3}\end{bmatrix}^T/{sin{\theta\over2}}}\end{cases}

从式子中而可以看出,旋转角\theta加上2pi,q变为相反数,这说明互为相反数的两个四元数可以表示同一个旋转。

4.2四元数运算

4.2.1加减

q_a{\pm}q_b=[s_a{\pm}s_b,v_a{\pm}v_b]

4.2.2乘法

四元数的乘法依然得到四元数
q_aq_b=[s_as_b-v_a^Tv_b,s_av_b+s_bv_a+v_a{\times}v_b]

4.2.3共轭

q_a^*=[s_a,-v_b]
共轭相乘,得到一个实四元数
q^*q=qq^*=[s_a^2+v^Tv,0]

4.2.4模长

||q_a||=\sqrt{s_a^2+x_a^2+x_b^2+x_c^2}

4.2.5逆

q^{-1}=q^*/||q^2||
四元数的逆和四元数相乘可以得到四元数1
qq^{-1}=q^{-1}q=1
并且具有这样的性质
(q_aq_b)^{-1}=q_b^{-1}q_a^{-1}

4.2.6数乘和点乘

数乘
kq=[ks,kv]
点乘
q_aq_b=s_as_b+x_ax_b+y_ay_b+z_az_b

4.3用四元数表示旋转

单位四元数可以表示任意一个旋转
这里将描述如何将一个三维点旋转,p是一个三维空间点,q为四元数,n,\theta是指定旋转。
第一步用四元数表示三维点:
p=[0,x,y,z]=[0,v]
第二步:用四元数表示这个旋转
q=[\cos{\theta\over2},n\sin{\theta\over2}]
第三步:旋转
p^{\prime}=qpq^{-1}
同理逆变换为
q^{-1}p^{\prime}q=p,
他的几何意义是把坐标转到四维空间,再转回来

4.4四元数到旋转矩阵的转换

四元数到旋转矩阵R
设四元数q=q_0+q_1i+q_2j+q_3k,对应的旋转矩阵R
R=\begin{bmatrix}1-2q_2^2-2q_3^2&2q_1q_2-2q_0q_3&{2q_1q_3+2q_0q_2}\\2q_1q_2+2q_0q_3&{1-2q_1^2-2q_2^2}&{2q_2q_3-2q_0q_1}\\{2q_1q_3-2q_0q_2}&{2q_2q_3+2q_0q_1}&{1-2q_1^2-2q_2^2}\end{bmatrix}
如何记住这个矩阵?
首先主对角线上那个没有减哪个,减二倍方。对于其他项,矩阵位置对应的q_iq_j在前,剩下的在后,符号联系行列式的计算,主对角线系列加号,副对角线系列是减号。
旋转矩阵R到四元数
注意到表示旋转的四元数都是单位四元数,因此,
q_0={\sqrt{tr(R)+1}\over2},q_1=\frac{m_{23}-m_{32}}{4q_0},q_1=\frac{m_{31}-m_{13}}{4q_0},q_1=\frac{m_{12}-m_{21}}{4q_0}

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 220,548评论 6 513
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 94,069评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 166,985评论 0 357
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,305评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,324评论 6 397
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 52,030评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,639评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,552评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 46,081评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,194评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,327评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 36,004评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,688评论 3 332
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,188评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,307评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,667评论 3 375
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,337评论 2 358

推荐阅读更多精彩内容