第四章探讨了欧拉角的机制。欧拉角将一次姿态转换为三次基本旋转,这套表示方法引入了顺序约定、固定轴与自身轴的区分,同时会引发万向节锁等坐标退化现象。为了寻找更直接的表示方式,有必要探究单次 SO(3)SO(3) 旋转自身的几何结构。轴角表示法正是基于这一需求建立的,其基本思路在于确定旋转中保持不动的方向,并计算垂直于该方向的平面所旋转的角度。这种几何构想在直觉上具有合理性,而要确立其理论基础,首先要解决的问题便是证明三维旋转存在一条不动轴。

不动轴

三维旋转在代数上由 SO(3)SO(3) 群中的矩阵 RR 描述:

RTR=I,detR=1R^TR=I,\quad \det R=1

其中 RTR=IR^TR=I 表明该变换保持空间的长度与角度,而 detR=1\det R=1 则表示三维取向保持不变。几何上的旋转轴表现为一条穿过原点的直线,且该直线上的点在旋转操作后保持在原位置。使用向量语言表述,即要求存在一个非零向量 u\mathbf u,满足方程:

Ru=uR\mathbf u=\mathbf u

上述方程指出,轴角模型所指的“不动轴”具有逐点固定特性。如果某条直线仅作为整体集合被映射回自身(例如发生方向反转的反向映射 Ru=uR\mathbf u=-\mathbf u),其上的点已发生位移。因此,寻找不动轴的目标,是找出变换矩阵 RR 对应于特征值为 11 的特征方向:

Ru=uR\mathbf u=\mathbf u

这一几何要求转化为代数问题:我们需要证明,任意 SO(3)SO(3) 矩阵都具有特征值 11

旋转轴

从直观的三维空间结构来看,二维旋转平面占据了两个正交方向。三维空间拥有三个独立的维度,当其中两个维度构成旋转平面发生转动时,会留下一个独立方向不参与该平面内的转动。欧拉旋转定理确立了这一事实:任何固定原点、保持长度与取向的三维旋转,都会保留一条过原点的不动轴。

这一几何直觉可以建立在代数推导之上。该论证包含两个主要环节:首先定位一个对应于实特征值的特征方向,随后利用长度与取向的守恒性确定该特征值为正数。考察特征方程,若存在非零向量 v\mathbf v 满足

Rv=λvR\mathbf v=\lambda\mathbf v

则可将其等价变形为

(RλI)v=0(R-\lambda I)\mathbf v=\mathbf 0

由于 v\mathbf v 为非零向量,矩阵 RλIR-\lambda I 为奇异矩阵,因此行列式为零:

det(RλI)=0\det(R-\lambda I)=0

对于 3×33\times3 实矩阵 RR,上述行列式展开为关于 λ\lambda 的三次实系数特征多项式。奇数次实多项式在全实数域内至少有一个实根。因此,RR 拥有至少一个实特征值 λ\lambda,并对应于非零特征向量 v\mathbf v,满足:

Rv=λvR\mathbf v=\lambda\mathbf v

根据正交矩阵的长度保持性质,存在如下等式:

Rv=v|R\mathbf v|=|\mathbf v|

另一方面,结合特征值的定义,可得:

λv=λv|\lambda\mathbf v|=|\lambda|\,|\mathbf v|

两相比较,特征值的绝对值满足:

λ=1|\lambda|=1

在实数范围内,满足该条件的特征值仅可能为

λ=1\lambda=1

λ=1\lambda=-1

λ=1\lambda=1 时,所求的不动方向即告确立。当首先确定的实特征值为 1-1 时,需要进一步考察矩阵的行列式。行列式的值恒等于所有特征值的乘积。基于 SO(3)SO(3) 矩阵的属性:

detR=1\det R=1

这意味着剩余两个特征值的乘积等于

1-1

根据实系数多项式的代数性质,复数特征值总是以共轭对的形式出现。一对共轭复数 a+bia+biabia-bi 的乘积为正实数:

(a+bi)(abi)=a2+b2(a+bi)(a-bi)=a^2+b^2

乘积等于 1-1 的要求说明剩余两个特征值需要为实数,且各自的绝对值同样受到正交性约束,为 +1+11-1。要使其乘积达到 1-1,这两个实特征值的组合为 +1+11-1。因此,在任何特征值组合下,矩阵包含至少一个 +1+1 的特征值。

代数推导由此得出结论:任意 SO(3)SO(3) 矩阵均拥有满足如下条件的非零向量 u\mathbf u

Ru=uR\mathbf u=\mathbf u

这在代数层面印证了欧拉旋转定理的主要内容,确立了三维空间中保持长度与取向的旋转操作包含一条穿过原点的不动方向。

垂直平面

确立不动轴之后,需要进一步考察轴外空间点的运动规律。设 u\mathbf u 为单位不动轴向量:

u=1,Ru=u|\mathbf u|=1,\quad R\mathbf u=\mathbf u

选取任意一个垂直于该轴的向量 v\mathbf v

uv=0\mathbf u\cdot\mathbf v=0

在经历旋转矩阵 RR 作用后,向量 v\mathbf v 映射为 RvR\mathbf v。通过内积的守恒性与轴的不动性,可以证明旋转后的向量仍然垂直于 u\mathbf u

uRv=RuRv=uv=0\mathbf u\cdot R\mathbf v =R\mathbf u\cdot R\mathbf v =\mathbf u\cdot\mathbf v =0

这一内积结果表明,与轴正交的向量在变换后仍位于初始的垂直平面内部,该垂直平面作为一个整体对该变换封闭。在确立了平面的不变性之后,需进一步界定平面内的变换类型。在该垂直平面内构建一组单位正交基 a,b\mathbf a,\mathbf b,并使其与不动轴 u\mathbf u 共同构成右手坐标系:

a,b,u\mathbf a,\mathbf b,\mathbf u

基于 u\mathbf u 的不变性,整个三维变换可分解为两部分操作:轴方向的恒等映射,以及垂直平面内的二维变换。平面内的变换继承了原矩阵的长度与角度保持性质,属于二维正交变换。由于 RSO(3)R\in SO(3) 具有行列式为 +1+1 的特性,它保持三维坐标架的取向。考虑到轴向分量固定,平面内的二维部分也保持自身的二维取向,表现为二维旋转 R(θ)R(\theta)

至此,三维旋转展现出双重结构:轴方向的坐标值保持不变,而垂直于轴的平面则经历一次特定的二维旋转。“绕轴转”这一术语指示了运动的空间基准。在此过程中,轴维持静止状态,承载角位移的实体是该垂直平面。不动轴在三维空间中同时承担着双重功能:它既标定了变换中的固定方向,又以法向量的身份界定了发生二维旋转的参考平面。

旋转角符号

基于不动轴与垂直平面内的旋转特性,可以建立轴角参数化表示。该表示法需要两个参量:其一为单位轴向量:

u\mathbf u

其二为旋转角度:

θ\theta

这组参数的几何意义明确:空间中任意向量沿 u\mathbf u 的分量维持恒定,而其垂直于 u\mathbf u 的分量在正交平面内经历角度为 θ\theta 的二维旋转。几何轴作为一条直线,在轴角参数化中被赋予了特定的方向向量 u\mathbf u。引入方向性旨在为旋转角度的正负提供参照基准,从而确立垂直平面内“正向旋转”的数学定义。

在垂直平面内任选一单位向量 a\mathbf a,并通过叉积构造向量:

b=u×a\mathbf b=\mathbf u\times\mathbf a

这组基向量与轴向量共同确立了一个标准的右手坐标系:

a,b,u\mathbf a,\mathbf b,\mathbf u

在此坐标框架下,正角度 θ\theta 的旋转遵循平面正交矩阵 R(θ)R(\theta) 的演化规则:

acosθa+sinθb\mathbf a\mapsto \cos\theta\,\mathbf a+\sin\theta\,\mathbf b

当轴向量取向反转为 u-\mathbf u 时,由右手定则生成的正交基向量将相应变为 b-\mathbf b,这导致垂直平面内的旋转正方向随之反转。因此,以下两组轴角参数:

(u,θ)(\mathbf u,\theta)

(u,θ)(-\mathbf u,-\theta)

映射到同一个三维旋转变换。这一等价性说明了轴角表示中的符号绑定关系:轴的方向向量与角度的正负号存在耦合关系,同步翻转两者的符号将对应同一旋转。

向量分解

针对三维空间中的任意向量 p\mathbf p,可依据选定的单位轴进行正交分解,将其表达为沿轴分量与垂直分量之和:

p=p+p\mathbf p=\mathbf p_\parallel+\mathbf p_\perp

其中,沿轴分量由点积形式给出:

p=(pu)u\mathbf p_\parallel=(\mathbf p\cdot\mathbf u)\mathbf u

而垂直于轴的分量则为余下的残差向量:

p=p(pu)u\mathbf p_\perp=\mathbf p-(\mathbf p\cdot\mathbf u)\mathbf u

轴角模型对这两种正交分量施加了不同的几何操作:

p保持不动\mathbf p_\parallel \quad\text{保持不动} p在垂直平面内转 θ\mathbf p_\perp \quad\text{在垂直平面内转 }\theta

这种分解方式呈现了轴角框架的主要优势。它提供了一种整体解析方式,将单次旋转解构为一维空间的恒等映射与二维子空间的正交旋转。基于这种几何分工,建立 RpR\mathbf p 的显式代数表达式成为下一步演算内容。在三维向量代数中,利用叉乘运算可以构造出垂直平面内的正交基准。“将 p\mathbf p_\perp 在垂直平面内旋转 9090^\circ”的向量可由如下叉积给出:

u×p\mathbf u\times\mathbf p_\perp

该结果向量既垂直于轴 u\mathbf u,也垂直于初始分量 p\mathbf p_\perp,从而位于给定的垂直平面内。遵循由 u\mathbf u 定义的右手正向准则,该叉积结果表征了 p\mathbf p_\perp 旋转 9090^\circ 后的状态。这一代数结构奠定了罗德里格斯公式的理论基础。通过融合轴向不变性与平面内旋转律,该几何结构将导向一个能够作用于任意向量 p\mathbf p 的简洁代数表达式。

旋转向量

轴角表示依赖于单位轴 u\mathbf u 与角度 θ\theta 两个参数。其中,单位轴向量包含三个空间分量,并受限于长度约束:

u=1|\mathbf u|=1

这一约束将轴方向的轨迹限定于二维的单位球面之上,从而提供 2 个独立的自由度。结合角度 θ\theta 贡献的 1 个自由度,整个轴角系统对应三维旋转所需的 3 个自由度。在工程与代数实践中,这两种量常被融合成单一的旋转向量:

ω=θu\boldsymbol\omega=\theta\mathbf u

旋转向量通过其方向同时编码了轴线方位与角度的正负号,并以其模长指示角位移大小。当 θ\theta 取负值时,ω\boldsymbol\omega 指向 u-\mathbf u 侧,其模长为 θ|\theta|。这种紧凑的记法将空间旋转转化为欧氏空间中的三维坐标。由此需要考察其在代数运算中的行为:旋转向量的线性加法与实际的旋转复合之间存在特定的对应关系。

针对有限角度旋转,群乘法的非交换性构成了主要约束。常规的向量加法满足交换律:

ω1+ω2=ω2+ω1\boldsymbol\omega_1+\boldsymbol\omega_2 = \boldsymbol\omega_2+\boldsymbol\omega_1

而三维空间中依次绕不同坐标轴旋转的结果依赖于执行序列。因此,旋转向量的定位明确体现为单次旋转姿态的坐标映射,有穷旋转的序列复合需要借助非线性代数结构来处理。

在微小旋转的极限框架下,旋转向量的线性叠加特性得以显现。当参与复合的旋转角度趋近于无穷小时,序列差异所引发的效应属于二阶小量。首次微小转动所产生的坐标系偏移量正比于第一个转角;随后的微小转动作用于此偏斜状态上,由此产生的耦合项量级等于两个微小角度的乘积。在只保留一阶线性主项的近似条件下,旋转向量表现出近似加法特性。这一线性化特质使其在描述瞬时动力学、连续演化以及角速度积分时发挥着重要作用,明确区分了无穷小生成元与大角度位移的代数边界。

零度退化

在轴角参数空间中,零旋转对应着一种退化状态。零旋转等价于单位恒等变换:

R=IR=I

在这种状态下,空间中的任意方向向量都满足特征方程:

Iu=uI\mathbf u=\mathbf u

这导致整个三维空间内的任意过原点直线都可作为不动轴。这一现象反映了零旋转所具备的各向同性对称特征,各个空间方向在该变换下表现出等价的几何地位。在使用轴角系统进行参数化时,零旋转展现为无穷多种坐标表示:

(u,0)(\mathbf u,0)

其中 u\mathbf u 可遍历单位球面上的任意一点。这些形式在代数上共同映射至同一个恒等旋转矩阵。引入旋转向量体系后:

ω=θu\boldsymbol\omega=\theta\mathbf u

零旋转状态被统一映射到原点:

ω=0\boldsymbol\omega=\mathbf 0

这种原点处的合并特性为微积分运算与无穷小分析提供了结构上的便利,同时也保留了零旋转的各向同性特征。在此情境下,代数映射的层次结构显现出来。物理旋转操作由单位矩阵 II 定义,而其相应的数学表达形式对应一个连续集合 (u,0)(\mathbf u,0)。单一群元素对应连续多组参数的现象构成了参数化系统中的多重覆盖特性。这一现象指明,轴角虽然符合三维空间的几何直觉,仍是一种带有拓扑退化的参数化方案。

反向轴

在轴角参数化域中,另一个引发拓扑等价现象的特殊边界为:

θ=π\theta=\pi

此即通常所指的 180180^\circ 旋转。围绕单位轴 u\mathbf u 的半圈旋转与围绕 u-\mathbf u 的半圈旋转构成相同的空间算子:

(u,π)=(u,π)(\mathbf u,\pi)=(-\mathbf u,\pi)

该等效性可通过正交分解模型获得代数验证。任意向量 p\mathbf p 被分解为:

p+p\mathbf p_\parallel+\mathbf p_\perp

在经历半圈旋转操作后,沿轴向的分量保持不变:

pp\mathbf p_\parallel\mapsto \mathbf p_\parallel

处于正交平面的分量则映射到反向位置:

pp\mathbf p_\perp\mapsto -\mathbf p_\perp

将参考轴由 u\mathbf u 切换至 u-\mathbf u,轴向所定义的空间直线与正交法平面保持不变。此外,半圈旋转的轨迹在终态上呈现出方向上的等价性:顺时针转动与逆时针转动将到达该二维平面的对跖点。这使得两种相反的轴向定义在 π\pi 角度下对应同一几何变换。

引入角度参量的周期性特征:

(u,θ)(\mathbf u,\theta)

(u,θ+2kπ)(\mathbf u,\theta+2k\pi)

(其中 kk 取整数)表征相同的物理动作。结合符号翻转规则与周期律,整个角度定义域可限制到半圈区间之内。对于分布在如下范围的角度:

π<θ<2π\pi<\theta<2\pi

u\mathbf u 旋转 θ\theta 的动作,等效于绕该轴施加如下负向旋转:

θ2π\theta-2\pi

进一步对该负角度及对应轴向量应用符号联合反转,即可导出新参数系:

(u,2πθ)(-\mathbf u,2\pi-\theta)

此时,新的角度参数落入标准区间:

0<2πθ<π0<2\pi-\theta<\pi

借由这种代数转换,越过半圈的旋转都可改写为绕反向基准轴的、小于半周的正角位移。为减少重复,通常取如下的角度约束域:

0θπ0\le \theta\le \pi

在几何图景中,该规范域可表示为一个半径为 π\pi 的实心球体:球内向量的方向表示旋转轴 u\mathbf u,而从球心测量的径向距离对应旋转角度 θ\theta。在此流形结构下,原点融合了全部的零旋转状态,而位于球体表面的对跖点则通过识别机制粘合为同一个半圈旋转算子。这一特征说明,轴角系统既映射了单次旋转的几何规律,也保留了自身特定的流形拓扑识别规则。

奇点问题

本章的论述确立了轴角方法的理论作用:它提供了理解单次旋转几何结构的视角,并为代数计算公式的推导提供了基础。相较于将整体姿态拆解为多次基本运动的欧拉角系统,轴角模型专注于刻画旋转自身几何属性------通过一条固定不动的空间直轴,以及其正交平面内发生的二维角位移来完成参数化。该范式呈现出简洁的几何图景,使得三维旋转的整体映射关系较为直观。

在应用轴角体系时,有两个拓扑与代数特性需要被纳入考量:

复合机制:当执行多重轴角旋转的复合操作时,合成矩阵所对应的轴向量与角位移遵循非线性复合律。三维空间中群元素的乘法运算具备序列依赖性,轴角参数反映了 SO(3)SO(3) 群的非交换代数结构。

维度特化:轴角表示法依赖于三维特有结构。在三维欧氏空间中,任意二维旋转平面的法空间为一维直线,从而允许单一法向量充当几何标识。当向四维及更高维度泛化时,正交变换能够表现为多组相互正交的二维平面的同步旋转。高维体系需要采用诸如复合旋转面或二重向量等更广义的几何工具进行系统表达。

总体而言,轴角准确捕捉了三维旋转中的基本几何结构,显示了法平面作为转动载体的物理意义。为了将该几何构造转化为可供算法调用的解析形式,正交分解理论导向一个具体的代数综合问题:在已知单位轴 u\mathbf u、旋转角度 θ\theta 以及目标向量 p\mathbf p 的前提下,如何将如下两项映射:

p不动\mathbf p_\parallel\quad\text{不动}

以及

p在平面内转 θ\mathbf p_\perp\quad\text{在平面内转 }\theta

整合为一个作用于原始向量 p\mathbf p 的闭式函数方程?这一代数综合过程将引入叉积算符,并从几何分解中推导出罗德里格斯旋转公式。

专题导航