[论文解读] Derivation of the Euler-Rodrigues formula for three-dimensional rotations from the general formula for four-dimensional rotations

Johan Ernest Mebius|ArXiv.org|Jan 26, 2007

Mathematics and Applications被引用 31

一句话总结

本文通过将一般4D旋转矩阵特化，推导出3D旋转的欧拉-罗德里格斯公式，其中条件 $a_{00} = 1$ 强制了四元数结构。证明了左、右四元数互为逆元，从而得到标准的3D旋转公式 $\mathbf{q} \to \mathbf{p} \mathbf{q} \mathbf{p}^{-1}$，并将推导扩展至3D旋转变换（rotoreflections），通过设定 $a_{00} = -1$ 实现。主要贡献是基于矩阵理论，从4D旋转原理严格推导出欧拉-罗德里格斯公式。

ABSTRACT

The general 4D rotation matrix is specialised to the general 3D rotation matrix by equating its leftmost top element (a00) to 1. Its associate matrix of products of the left-hand and right-hand quaternion components is specialised correspondingly. Inequalities involving the angles through which the coordinate axes in 3D space are displaced are used to prove that the left-hand and the right-hand quaternions are each other's inverses, thus proving the Euler-Rodrigues formula. A general procedure to determine the Euler parameters of a given 3D rotation matrix is sketched. By equating the leftmost top element to -1 instead of +1 in the general 4D rotation matrix, one proves the counterpart of the Euler-Rodrigues formula for 3D rotoreflections. Keywords: Euler--Rodrigues formula, Euler parameters, quaternions, four--dimensional rotations, three--dimensional rotations, rotoreflections

研究动机与目标

• 从一般4D旋转矩阵形式推导出3D旋转的欧拉-罗德里格斯公式。
• 通过四元数分解，在4D旋转矩阵与3D旋转之间建立数学联系。
• 在条件 $a_{00} = 1$ 下，证明4D旋转分解中的左、右四元数互为逆元，从而导出标准3D旋转公式。
• 通过设定 $a_{00} = -1$，将推导扩展至3D旋转变换，提供相应的公式。

提出的方法

• 通过强制 $a_{00} = 1$，将一般4D旋转矩阵 $A = M_L M_R$ 特化，使其简化为3D旋转矩阵。
• 利用乘积 $ap, aq, ar, as, \dots, dp, dq, dr, ds$ 的关联矩阵 $M$，将4D旋转参数与3D旋转四元数联系起来。
• 应用旋转角的不等式，证明 $ap \leq 0$ 且 $bq, cr, ds \geq 0$，从而约束四元数分量之间的符号关系。
• 由恒等式 $a^2 = p^2$ 及 $bq, cr, ds$ 的非负性，推导出 $a = -p$ 与 $q = b, r = c, s = d$，证明左、右四元数互为逆元。
• 利用矩阵 $M$ 通过 $a = (a_{11} + a_{22} + a_{33} + 1)/4$ 等公式，从给定的3D旋转矩阵重构欧拉参数 $a, b, c, d$。
• 通过设定 $a_{00} = -1$ 推导出旋转变换的对应公式，得到带符号特征 $(-a^2 - b^2 + c^2 + d^2, \dots)$ 的修正矩阵表达式。

实验结果

研究问题

• RQ1 如何从一般4D旋转矩阵形式严格推导出3D旋转的欧拉-罗德里格斯公式？
• RQ2 4D旋转矩阵元素的何种条件会强制左、右四元数互为逆元，从而简化为3D旋转？
• RQ3 如何从关联矩阵 $M$ 计算给定3D旋转矩阵的欧拉参数？
• RQ4 3D旋转变换公式的数学结构是什么？它与4D旋转框架有何关联？
• RQ5 3D旋转中的角度不等式在确定四元数分量之间符号关系时起到什么作用？

主要发现

• 在4D旋转矩阵中，条件 $a_{00} = 1$ 强制左、右四元数互为逆元，从而导出欧拉-罗德里格斯公式 $\mathbf{q} \to \mathbf{p} \mathbf{q} \mathbf{p}^{-1}$。
• 欧拉参数通过 $a = (a_{11} + a_{22} + a_{33} + 1)/4$，$b = (a_{11} - a_{22} + a_{33} - 1)/4$，$c = (a_{11} + a_{22} - a_{33} - 1)/4$，$d = (-a_{11} + a_{22} + a_{33} - 1)/4$ 从3D旋转矩阵恢复，且其平方和为1。
• 在约束 $ap \leq 0$ 和 $bq, cr, ds \geq 0$ 下，推导证明了 $a = -p$，$q = b$，$r = c$，$s = d$，确认了逆关系。
• 对于旋转变换，设定 $a_{00} = -1$ 得到矩阵 $\begin{bmatrix} -a^2-b^2+c^2+d^2 & 2ad-2bc & -2ac-2bd \\ -2ad-2bc & -a^2+b^2-c^2+d^2 & 2ab-2cd \\ 2ac-2bd & -2ab-2cd & -a^2+b^2+c^2-d^2 \end{bmatrix}$，且满足 $a^2 + b^2 + c^2 + d^2 = 1$。
• 证明中使用的角度不等式表明，每个轴位移角的余弦值至多为 $\cos\alpha$，这支持了对 $bq, cr, ds$ 的符号约束。
• 4D旋转矩阵的关联矩阵 $M$ 秩为1且范数为1，使得能够从矩阵元素重构四元数分量。