第五章对一次三维旋转的几何结构进行了分析。给定一条有方向的单位轴 u \mathbf u u 和一个角度 θ \theta θ ,任意向量 p \mathbf p p 可分解为沿轴方向的分量与垂直于轴的分量。在旋转过程中,轴向分量保持静止,而垂直分量则在一个二维平面内发生转动。这种分解方式揭示了单次旋转的基本结构,并由此引出一个具体问题:如何在独立于临时平面坐标系的情况下,构造出表示旋转后结果的向量 p ′ \mathbf p' p ′ ?罗德里格斯公式正是为解决这一问题而提出的,其主要目的在于将几何分解过程转化为一个统一的数学表达式,使其能够作用于任意向量,从而在一个等式中呈现旋转的作用机制。
向量正交分解
设 u \mathbf u u 为满足条件的单位轴:
∣ u ∣ = 1 |\mathbf u|=1 ∣ u ∣ = 1
任意向量 p \mathbf p p 可向该轴进行投影,其轴向分量表示为:
p ∥ = ( p ⋅ u ) u \mathbf p_\parallel=(\mathbf p\cdot\mathbf u)\mathbf u p ∥ = ( p ⋅ u ) u
其中内积 p ⋅ u \mathbf p\cdot\mathbf u p ⋅ u 确定了 p \mathbf p p 在 u \mathbf u u 方向上的带符号长度,该标量与单位向量 u \mathbf u u 相乘即可得到沿轴方向的向量。通过将总向量减去该轴向投影,可以得到对应的垂直分量:
p ⊥ = p − p ∥ = p − ( p ⋅ u ) u \mathbf p_\perp=\mathbf p-\mathbf p_\parallel
=\mathbf p-(\mathbf p\cdot\mathbf u)\mathbf u p ⊥ = p − p ∥ = p − ( p ⋅ u ) u
对其进行正交性检验,计算两者内积:
u ⋅ p ⊥ = u ⋅ p − ( p ⋅ u ) ( u ⋅ u ) \mathbf u\cdot\mathbf p_\perp
=\mathbf u\cdot\mathbf p-(\mathbf p\cdot\mathbf u)(\mathbf u\cdot\mathbf u) u ⋅ p ⊥ = u ⋅ p − ( p ⋅ u ) ( u ⋅ u )
由单位向量的模长性质可知:
u ⋅ u = 1 \mathbf u\cdot\mathbf u=1 u ⋅ u = 1
于是该内积结果化简为:
u ⋅ p ⊥ = u ⋅ p − p ⋅ u = 0 \mathbf u\cdot\mathbf p_\perp
=\mathbf u\cdot\mathbf p-\mathbf p\cdot\mathbf u=0 u ⋅ p ⊥ = u ⋅ p − p ⋅ u = 0
这一结果从代数上证明了向量 p ⊥ \mathbf p_\perp p ⊥ 位于与轴垂直的平面内。
轴角旋转机制对上述两个分量采取了分离的处理方式:
p ∥ ↦ p ∥ \mathbf p_\parallel\mapsto \mathbf p_\parallel p ∥ ↦ p ∥
p ⊥ ↦ p ⊥ 在垂直平面里转 θ \mathbf p_\perp\mapsto
\mathbf p_\perp\text{ 在垂直平面里转 }\theta p ⊥ ↦ p ⊥ 在垂直平面里转 θ
这种分解处理蕴含着一种几何视角的转换。在传统物理描述中,“绕轴旋转”将分析重心置于旋转轴本身;而在罗德里格斯公式的推导中,主要运算聚焦于垂直平面内的二维旋转。在此过程中,旋转轴的作用体现为确立方向基准与提供投影参考,而几何变换发生在正交平面的分量上。
平面旋转
根据前述章节的推导,在给定的二维平面内,假设存在一组正交单位基底:
a , b \mathbf a,\quad \mathbf b a , b
且正向旋转定义为由 a \mathbf a a 转向 b \mathbf b b ,则任意平面向量:
x a + y b x\mathbf a+y\mathbf b x a + y b
在旋转角度 θ \theta θ 后将变换为:
( x cos θ − y sin θ ) a + ( x sin θ + y cos θ ) b (x\cos\theta-y\sin\theta)\mathbf a
+(x\sin\theta+y\cos\theta)\mathbf b ( x cos θ − y sin θ ) a + ( x sin θ + y cos θ ) b
这种代数表达可以进一步转换为几何描述。若将初始平面向量记为 q \mathbf q q ,并令其逆时针旋转 90 ∘ 90^\circ 9 0 ∘ 后的正交方向向量为 J q J\mathbf q J q ,则该向量旋转 θ \theta θ 的结果可等效表示为:
q cos θ + ( J q ) sin θ \mathbf q\cos\theta+(J\mathbf q)\sin\theta q cos θ + ( J q ) sin θ
这一恒等式表明,二维旋转能够通过两个相互正交的方向向量进行线性组合来实现:原方向分量对应 cos θ \cos\theta cos θ 的权重,而正交方向分量则对应 sin θ \sin\theta sin θ 的权重。
将上述结论推广至垂直于 u \mathbf u u 的三维空间截面中,此时平面向量即为 p ⊥ \mathbf p_\perp p ⊥ 。为了应用该几何结构,需要在此三维语境下构造出 p ⊥ \mathbf p_\perp p ⊥ 在其所在平面内沿正方向旋转 90 ∘ 90^\circ 9 0 ∘ 后的正交向量。向量叉乘满足这一条件:
u × p ⊥ \mathbf u\times\mathbf p_\perp u × p ⊥
由于叉积结果同时垂直于 u \mathbf u u 和 p ⊥ \mathbf p_\perp p ⊥ ,该向量仍位于与 u \mathbf u u 垂直的平面内。同时,根据右手定则,该叉乘结果给出了由旋转轴 u \mathbf u u 所定义的平面正向 90 ∘ 90^\circ 9 0 ∘ 偏转方向。因此,垂直分量 p ⊥ \mathbf p_\perp p ⊥ 旋转 θ \theta θ 后的状态可表达为:
p ⊥ cos θ + ( u × p ⊥ ) sin θ \mathbf p_\perp\cos\theta
+(\mathbf u\times\mathbf p_\perp)\sin\theta p ⊥ cos θ + ( u × p ⊥ ) sin θ
结合保持静止的轴向分量,整体旋转变换可写为:
p ′ = p ∥ + p ⊥ cos θ + ( u × p ⊥ ) sin θ \mathbf p'
=\mathbf p_\parallel
+\mathbf p_\perp\cos\theta
+(\mathbf u\times\mathbf p_\perp)\sin\theta p ′ = p ∥ + p ⊥ cos θ + ( u × p ⊥ ) sin θ
至此,基础的旋转公式框架已经得到:整体变换由保持不变的轴向投影、平面内原方向的余弦收缩以及平面内正交方向的正弦延伸三个独立部分共同叠加而成。
叉乘算符
在上述公式框架中,仍然包含依赖分解的中间变量 p ⊥ \mathbf p_\perp p ⊥ 。为了使旋转算子能够作用于初始向量 p \mathbf p p ,需要将该分解项化简。首先处理叉乘项:
u × p ⊥ = u × ( p − ( p ⋅ u ) u ) \mathbf u\times\mathbf p_\perp
=\mathbf u\times\bigl(\mathbf p-(\mathbf p\cdot\mathbf u)\mathbf u\bigr) u × p ⊥ = u × ( p − ( p ⋅ u ) u )
根据线性运算法则展开,得到:
u × p − ( p ⋅ u ) ( u × u ) \mathbf u\times\mathbf p
-(\mathbf p\cdot\mathbf u)(\mathbf u\times\mathbf u) u × p − ( p ⋅ u ) ( u × u )
由同一向量的自叉乘性质可知其等于零向量:
u × u = 0 \mathbf u\times\mathbf u=\mathbf 0 u × u = 0
所以该项简化为:
u × p ⊥ = u × p \mathbf u\times\mathbf p_\perp
=\mathbf u\times\mathbf p u × p ⊥ = u × p
这一代数结果与几何直觉一致:由于向量 p \mathbf p p 的轴向分量与 u \mathbf u u 共线,其在叉积运算中产生的正交扩张效应为零。因此,叉乘算子滤除了沿轴方向的投影,分离出参与垂直平面转动的有效几何部分。
随后,将定义式:
p ∥ = ( p ⋅ u ) u \mathbf p_\parallel=(\mathbf p\cdot\mathbf u)\mathbf u p ∥ = ( p ⋅ u ) u
与
p ⊥ = p − ( p ⋅ u ) u \mathbf p_\perp=\mathbf p-(\mathbf p\cdot\mathbf u)\mathbf u p ⊥ = p − ( p ⋅ u ) u
代回总旋转公式:
p ′ = ( p ⋅ u ) u + ( p − ( p ⋅ u ) u ) cos θ + ( u × p ) sin θ \mathbf p'
=(\mathbf p\cdot\mathbf u)\mathbf u
+\bigl(\mathbf p-(\mathbf p\cdot\mathbf u)\mathbf u\bigr)\cos\theta
+(\mathbf u\times\mathbf p)\sin\theta p ′ = ( p ⋅ u ) u + ( p − ( p ⋅ u ) u ) cos θ + ( u × p ) sin θ
对该等式进行同类项重组,便推导出标准的罗德里格斯旋转公式:
p ′ = p cos θ + ( u × p ) sin θ + u ( u ⋅ p ) ( 1 − cos θ ) \mathbf p'
=\mathbf p\cos\theta
+(\mathbf u\times\mathbf p)\sin\theta
+\mathbf u(\mathbf u\cdot\mathbf p)(1-\cos\theta) p ′ = p cos θ + ( u × p ) sin θ + u ( u ⋅ p ) ( 1 − cos θ )
等价表达式为:
p ′ = p cos θ + ( u × p ) sin θ + u ( u ⋅ p ) ( 1 − cos θ ) \boxed{
\mathbf p'
=\mathbf p\cos\theta
+(\mathbf u\times\mathbf p)\sin\theta
+\mathbf u(\mathbf u\cdot\mathbf p)(1-\cos\theta)
} p ′ = p cos θ + ( u × p ) sin θ + u ( u ⋅ p ) ( 1 − cos θ )
该公式的每一项都具有明确的几何来源。第一项 p cos θ \mathbf p\cos\theta p cos θ 将总体向量进行缩放,同时提取了轴向与垂直分量的余弦投影。第二项 ( u × p ) sin θ (\mathbf u\times\mathbf p)\sin\theta ( u × p ) sin θ 显式表达了垂直平面内发生的 90 ∘ 90^\circ 9 0 ∘ 正交偏转。第三项 u ( u ⋅ p ) ( 1 − cos θ ) \mathbf u(\mathbf u\cdot\mathbf p)(1-\cos\theta) u ( u ⋅ p ) ( 1 − cos θ ) 则是对轴向分量的补偿。由于旋转要求轴向分量保持不变,而第一项仅提取了该分量的 cos θ \cos\theta cos θ 倍,因此需要通过叠加 1 − cos θ 1-\cos\theta 1 − cos θ 倍的方式恢复其原长度。这种项数结构源于对“垂直平面二维旋转”与“轴向长度守恒”两种几何约束的代数综合过程。
边界检验
在确立公式代数形式后,可通过代入特定边界条件来验证其几何一致性。考虑向量 p \mathbf p p 平行于旋转轴的特殊情形,例如:
p = a u \mathbf p=a\mathbf u p = a u
此时存在叉积:
u × p = 0 \mathbf u\times\mathbf p=\mathbf 0 u × p = 0
以及内积:
u ⋅ p = a \mathbf u\cdot\mathbf p=a u ⋅ p = a
将这些条件代入罗德里格斯公式,计算可得:
p ′ = a u cos θ + 0 + a u ( 1 − cos θ ) = a u \mathbf p'
=a\mathbf u\cos\theta
+\mathbf 0
+a\mathbf u(1-\cos\theta)
=a\mathbf u p ′ = a u cos θ + 0 + a u ( 1 − cos θ ) = a u
该结果反映了轴上点在旋转变换下保持位置不变的刚体性质。
若向量 p \mathbf p p 垂直于旋转轴,此时内积满足:
u ⋅ p = 0 \mathbf u\cdot\mathbf p=0 u ⋅ p = 0
导致公式中的第三项消失,变换等式退化为:
p ′ = p cos θ + ( u × p ) sin θ \mathbf p'
=\mathbf p\cos\theta+(\mathbf u\times\mathbf p)\sin\theta p ′ = p cos θ + ( u × p ) sin θ
这一退化形式对应了纯二维平面旋转的代数结构。
进一步考察旋转角度的边界值,当 θ = 0 \theta=0 θ = 0 时,取三角函数值为:
cos 0 = 1 , sin 0 = 0 \cos0=1,\quad \sin0=0 cos 0 = 1 , sin 0 = 0
代入公式得到:
p ′ = p \mathbf p'=\mathbf p p ′ = p
当 θ = π \theta=\pi θ = π 时,存在:
cos π = − 1 , sin π = 0 \cos\pi=-1,\quad \sin\pi=0 cos π = − 1 , sin π = 0
公式化简为:
p ′ = − p + 2 u ( u ⋅ p ) \mathbf p'
=-\mathbf p+2\mathbf u(\mathbf u\cdot\mathbf p) p ′ = − p + 2 u ( u ⋅ p )
通过将正交分解式:
p = p ∥ + p ⊥ \mathbf p=\mathbf p_\parallel+\mathbf p_\perp p = p ∥ + p ⊥
代入该结果进行重组,可得:
− p + 2 p ∥ = p ∥ − p ⊥ -\mathbf p+2\mathbf p_\parallel
=\mathbf p_\parallel-\mathbf p_\perp − p + 2 p ∥ = p ∥ − p ⊥
该表达式说明,半圈旋转的效果表现为轴向分量守恒而垂直分量反向。上述系列验证表明,该代数公式保持了轴角旋转的几何物理意义。
反对称矩阵
尽管向量形式的罗德里格斯公式具有直观的几何解释,许多实际工程计算与理论推导要求将旋转变换表达为标准的矩阵乘法形式。这就要求将其转换为一个等效的 3 × 3 3\times3 3 × 3 旋转矩阵,而实现这一转换的主要途径在于引入反对称矩阵。
给定单位轴向量:
u = ( u x , u y , u z ) T \mathbf u=(u_x,u_y,u_z)^T u = ( u x , u y , u z ) T
可构造如下结构矩阵:
[ u ] × = [ 0 − u z u y u z 0 − u x − u y u x 0 ] [\mathbf u]_\times=
\begin{bmatrix}
0&-u_z&u_y\\
u_z&0&-u_x\\
-u_y&u_x&0
\end{bmatrix} [ u ] × = 0 u z − u y − u z 0 u x u y − u x 0
其代数作用定义为叉乘操作:
[ u ] × p = u × p [\mathbf u]_\times\mathbf p=\mathbf u\times\mathbf p [ u ] × p = u × p
该矩阵的一个基本性质在于其转置矩阵满足如下关系:
[ u ] × T = − [ u ] × [\mathbf u]_\times^T=-[\mathbf u]_\times [ u ] × T = − [ u ] ×
具备该特征的矩阵统称为反对称矩阵。
在此引入反对称矩阵来自线性代数的结构。由于叉积操作对于输入向量 p \mathbf p p 具有线性特征,在固定参量 u \mathbf u u 的条件下,与 u \mathbf u u 的叉乘运算构成了一个确定的线性算子,因而拥有与其对应的矩阵表示。这里的 [ u ] × [\mathbf u]_\times [ u ] × 正是这一线性算子的代数表示。
为进一步化简公式,考察该算子的平方项运算。对于任意向量 p \mathbf p p ,根据向量三重积的恒等式,可以展开得到:
u × ( u × p ) = u ( u ⋅ p ) − p ( u ⋅ u ) \mathbf u\times(\mathbf u\times\mathbf p)
=\mathbf u(\mathbf u\cdot\mathbf p)-\mathbf p(\mathbf u\cdot\mathbf u) u × ( u × p ) = u ( u ⋅ p ) − p ( u ⋅ u )
引入单位轴的归一化条件 u ⋅ u = 1 \mathbf u\cdot\mathbf u=1 u ⋅ u = 1 ,该恒等式即缩减为:
u × ( u × p ) = u ( u ⋅ p ) − p \mathbf u\times(\mathbf u\times\mathbf p)
=\mathbf u(\mathbf u\cdot\mathbf p)-\mathbf p u × ( u × p ) = u ( u ⋅ p ) − p
将其转换为矩阵乘法语言,便得到等效关系式:
[ u ] × 2 p = u ( u ⋅ p ) − p [\mathbf u]_\times^2\mathbf p
=\mathbf u(\mathbf u\cdot\mathbf p)-\mathbf p [ u ] × 2 p = u ( u ⋅ p ) − p
利用上述代数结果,可以对罗德里格斯公式进行重写:
p ′ = p cos θ + [ u ] × p sin θ + u ( u ⋅ p ) ( 1 − cos θ ) = p + [ u ] × p sin θ + ( u ( u ⋅ p ) − p ) ( 1 − cos θ ) = ( I + sin θ [ u ] × + ( 1 − cos θ ) [ u ] × 2 ) p \begin{aligned}
\mathbf p'
&=\mathbf p\cos\theta
+[\mathbf u]_\times\mathbf p\sin\theta
+\mathbf u(\mathbf u\cdot\mathbf p)(1-\cos\theta)\\
&=\mathbf p
+[\mathbf u]_\times\mathbf p\sin\theta
+\bigl(\mathbf u(\mathbf u\cdot\mathbf p)-\mathbf p\bigr)(1-\cos\theta)\\
&=\bigl(I+\sin\theta[\mathbf u]_\times
+(1-\cos\theta)[\mathbf u]_\times^2\bigr)\mathbf p
\end{aligned} p ′ = p cos θ + [ u ] × p sin θ + u ( u ⋅ p ) ( 1 − cos θ ) = p + [ u ] × p sin θ + ( u ( u ⋅ p ) − p ) ( 1 − cos θ ) = ( I + sin θ [ u ] × + ( 1 − cos θ ) [ u ] × 2 ) p
由此,提取右侧的算子部分,可以写出标准的三维旋转矩阵形式:
R = I + sin θ [ u ] × + ( 1 − cos θ ) [ u ] × 2 \boxed{
R
=I+\sin\theta[\mathbf u]_\times
+(1-\cos\theta)[\mathbf u]_\times^2
} R = I + sin θ [ u ] × + ( 1 − cos θ ) [ u ] × 2
上述表达式给出了罗德里格斯公式的矩阵形式,其数学作用在于将空间中的轴角结构转化为特殊正交群 S O ( 3 ) SO(3) S O ( 3 ) 的矩阵表达。在这一结构中,轴角参数提供了几何图像,矩阵形式确立了标准的线性变换框架,而反对称矩阵则将空间叉乘运算嵌入到规范的矩阵乘法体系之中。
从矩阵反求轴角参数时,常用两条由罗德里格斯公式导出的关系:
cos θ = tr R − 1 2 , R − R T = 2 sin θ [ u ] × . \cos\theta=\frac{\operatorname{tr}R-1}{2},
\qquad
R-R^T=2\sin\theta[\mathbf u]_\times. cos θ = 2 tr R − 1 , R − R T = 2 sin θ [ u ] × .
这些关系在 0 < θ < π 0<\theta<\pi 0 < θ < π 的常规区间内可以直接给出角度与带符号轴向。当 θ = 0 \theta=0 θ = 0 时,旋转轴没有唯一性;当 θ = π \theta=\pi θ = π 时,sin θ = 0 \sin\theta=0 sin θ = 0 ,反对称部分无法恢复轴向符号,需要转而利用 R + I R+I R + I 的列空间或等价的特征向量方法。由此可见,轴角反求公式本身也带有边界退化条件。
叉乘局限
尽管罗德里格斯公式的推导展现了代数对称性,且叉乘运算在其中具有作用,但有必要明确这种数学表述的几何适用边界。叉乘在此处发挥效用,依赖于三维空间的一项特定属性:任意二维平面都可由一条唯一的正交法向轴标记,且运算 u × p ⊥ \mathbf u\times\mathbf p_\perp u × p ⊥ 诱导了平面内的 90 ∘ 90^\circ 9 0 ∘ 正交偏转。
当维度发生变化时,该代数机制的适用前提随之改变。在二维空间中,旋转所在平面即覆盖了整个欧几里得空间,此时传统的法向轴需要向外扩展出额外维度或退化为纯标量。
在四维及更高维度空间中,二维旋转平面的正交补空间同样展现为一个多维子空间(例如在四维空间中表现为另一个二维平面)。此时,通过单一方向轴来标记二维平面的传统方法被正交多重结构所取代。
综上所述,三维叉积可看作一种特定的降维代换工具,它将二维截面内的 90 ∘ 90^\circ 9 0 ∘ 旋转操作等效转化为依托法向轴的外积运算。这种表述依赖于三维空间的对偶同构性质。在更广泛的维度框架下,分析旋转变换通常需要回归到张量代数与外代数的普适体系中。
这里的“对偶同构”在更标准的术语中常称为 Hodge 对偶。它把三维空间中的有向平面信息转写为法向轴信息。正文继续使用叉积,是因为它在三维旋转计算中最直接;第十三章会回到有向平面本身,说明这一记法背后的统一结构。
维度的特殊性也从侧面说明,旋转操作的主要几何载体是有向平面,而旋转轴是三维空间特有的一种等效几何标签。这一观点将在后续章节中继续展开。
小角度近似
罗德里格斯公式的另一项理论作用在于,它给出了微小旋转变换的一阶近似结构。
在旋转角 θ \theta θ 趋于微小的极限条件下,三角函数可通过泰勒展开进行线性化近似:
sin θ ≈ θ \sin\theta\approx\theta sin θ ≈ θ
并且:
1 − cos θ ≈ θ 2 2 1-\cos\theta\approx \frac{\theta^2}{2} 1 − cos θ ≈ 2 θ 2
在仅保留与 θ \theta θ 成正比的一阶主导项,并在截断高阶无穷小量后,旋转矩阵的表达式演变为:
R ≈ I + θ [ u ] × R
\approx I+\theta[\mathbf u]_\times R ≈ I + θ [ u ] ×
该近似式的代数结构表明,一次微小旋转在一阶近似下等效于在单位恒等变换 I I I 的基础上,沿着反对称矩阵所定义的方向进行微小的线性叠加。
若引入旋转向量:
ω = θ u \boldsymbol\omega=\theta\mathbf u ω = θ u
作为变量,则相应的反对称矩阵项可改写为:
θ [ u ] × = [ ω ] × \theta[\mathbf u]_\times=[\boldsymbol\omega]_\times θ [ u ] × = [ ω ] ×
由此推导出微小旋转的一阶线性化表达形式为:
R ≈ I + [ ω ] × R\approx I+[\boldsymbol\omega]_\times R ≈ I + [ ω ] ×
这一推论说明了旋转群在单位元附近的局部结构。对于大角度的有限旋转,三维旋转的复合过程保留非交换性,其代数运算呈现非线性行为。当旋转幅度收缩至微小邻域时,非线性效应只在高阶项中出现,旋转空间的局部结构可以近似为平直的线性空间。此时,系统脱离初始恒等姿态的微小演化方向由对应的反对称矩阵刻画。
这种局部线性化的视角引发了对系统连续演化过程的理论探讨。当物理旋转被视为随时间连续推进的动态过程时,其演化轨迹上的每一个瞬时状态均可等效为一次极小旋转作用。这种瞬时微分变化可以通过相应的微小旋转量进行数学描述,而反对称矩阵演化为表征旋转角速度的代数工具。
围绕这一连续微积分框架的构建,下一章将分析连续旋转动力学,并引入李代数的理论体系。
专题导航