2.5 的角和轴表示虽然直观,却留下两个缺陷:一是不唯一 ——( ϑ , r ) (\vartheta,\boldsymbol r) ( ϑ , r ) 与 ( − ϑ , − r ) (-\vartheta,-\boldsymbol r) ( − ϑ , − r ) 同姿(定理 2.5.A );二是 sin ϑ = 0 \sin\vartheta=0 sin ϑ = 0 处逆解奇异 (定理 2.5.B )。本节要讲的单位四元数 (unit quaternion,又称欧拉参数)用 4 个参数就把这两点一并修掉——它至今仍是全姿态机动平台(卫星等)与计算机图形学的标准姿态表示。
本节定位(教材补充 · 课程未授)
第 2 周录播未展开本节,仅在讲 RPY 奇点时提及「卫星等全姿态机动飞行器改用四元数」(录播 0:21:26)。为使本编完整覆盖原书 §2.1–§2.12,这里按原书 §2.6 收录本节;考试不作要求 ,但理解它与 2.5 的关系,能解释「为什么工程上会放弃直观的欧拉角」。
定义 2.6.1 · 单位四元数
单位四元数 (欧拉参数)定义为四元组 Q = { η , ϵ } Q=\{\eta,\boldsymbol\epsilon\} Q = { η , ϵ } ,其中
η = cos ϑ 2 , ϵ = sin ϑ 2 r ( 原书式 2.30 – 2.31 ) \eta=\cos\frac{\vartheta}{2},\qquad \boldsymbol\epsilon=\sin\frac{\vartheta}{2}\,\boldsymbol r\qquad(\text{原书式 }2.30\text{–}2.31) η = cos 2 ϑ , ϵ = sin 2 ϑ r ( 原书式 2.30 – 2.31 )
η \eta η 称标量部分 ,ϵ = [ ϵ x ϵ y ϵ z ] T \boldsymbol\epsilon=[\epsilon_x\ \epsilon_y\ \epsilon_z]^T ϵ = [ ϵ x ϵ y ϵ z ] T 称向量部分 ,二者受单位约束
η 2 + ϵ x 2 + ϵ y 2 + ϵ z 2 = 1 ( 原书式 2.32 ) \eta^2+\epsilon_x^2+\epsilon_y^2+\epsilon_z^2=1\qquad(\text{原书式 }2.32) η 2 + ϵ x 2 + ϵ y 2 + ϵ z 2 = 1 ( 原书式 2.32 )
故独立参数仍为 3 个(4 参数 − 1 约束),与最简表示的要求不矛盾。
原理 2.6.A · 四元数修掉了哪些缺陷
关键差别在半角打包 :代入 − ϑ -\vartheta − ϑ 与 − r -\boldsymbol r − r ,得 cos − ϑ 2 = cos ϑ 2 \cos\frac{-\vartheta}{2}=\cos\frac{\vartheta}{2} cos 2 − ϑ = cos 2 ϑ 、sin − ϑ 2 ( − r ) = sin ϑ 2 r \sin\frac{-\vartheta}{2}(-\boldsymbol r)=\sin\frac{\vartheta}{2}\boldsymbol r sin 2 − ϑ ( − r ) = sin 2 ϑ r ——同一结果 。
消除了角和轴的表示不唯一 ——( ϑ , r ) (\vartheta,\boldsymbol r) ( ϑ , r ) 与 ( − ϑ , − r ) (-\vartheta,-\boldsymbol r) ( − ϑ , − r ) 是同一旋转的两组不同参数(定理 2.5.A );半角打包后这两组参数映射到同一个 四元数,歧义坍缩为一个表示(原书:a rotation by − ϑ -\vartheta − ϑ about − r -\boldsymbol r − r gives the same quaternion as a rotation by ϑ \vartheta ϑ about r \boldsymbol r r ; this solves the nonuniqueness problem);
剩下的双重覆盖 Q ∼ − Q Q\sim-Q Q ∼ − Q 来自 ϑ ∼ ϑ + 2 π \vartheta\sim\vartheta+2\pi ϑ ∼ ϑ + 2 π ——同一个物理旋转,转角再多加满一整圈(ϑ + 2 π \vartheta+2\pi ϑ + 2 π )时 η → − η \eta\to-\eta η → − η 、ϵ → − ϵ \boldsymbol\epsilon\to-\boldsymbol\epsilon ϵ → − ϵ ,故 Q Q Q 与 − Q -Q − Q 表示同一旋转。逆解式 2.34 的 η ≥ 0 \eta\ge0 η ≥ 0 约定(ϑ ∈ [ − π , π ] \vartheta\in[-\pi,\pi] ϑ ∈ [ − π , π ] )为每个旋转矩阵选定唯一 一个四元数(η = 0 \eta=0 η = 0 边界处由 sgn \operatorname{sgn} sgn 规则仍给出确定答案);
消除了逆解奇点 ——式 2.34–2.35 对任何旋转矩阵都有效,没有 sin ϑ = 0 \sin\vartheta=0 sin ϑ = 0 那样的失效点(对照 定理 2.5.B )。
把原理 2.6.A 的前两条结论画成一张「参数 → 旋转 → 四元数」的映射图,关系会立刻清楚:半角打包把角和轴的两组 参数坍缩成一个 四元数;而 ϑ ∼ ϑ + 2 π \vartheta\sim\vartheta+2\pi ϑ ∼ ϑ + 2 π 留下的 Q ∼ − Q Q\sim-Q Q ∼ − Q 双覆盖,则由 η ≥ 0 \eta\ge0 η ≥ 0 约定裁决。
四元数如何修复角和轴的两个缺陷
左侧映射图:参数组 (ϑ,r) 与 (−ϑ,−r) 对应同一旋转 R,经半角打包映射到同一个四元数 Q;ϑ 与 ϑ+2π 对应同一旋转但映射到 Q 与 −Q 两个对径点。右侧为单位球面示意,Q 与 −Q 为对径两点。
修复 ①:表示不唯一 → 半角打包坍缩为一个 Q
(ϑ, r)
(−ϑ, −r)
同一旋转 R
同一个 Q
{cos(ϑ/2), sin(ϑ/2)·r}
仍保留:ϑ ∼ ϑ+2π → Q ∼ −Q(双覆盖)
ϑ
ϑ + 2π
同一旋转 R
Q
−Q
Q
−Q
ϑ 增加 2π
单位约束 η² + ‖ε‖² = 1(单位球面示意)
取 η ≥ 0 → 每个旋转矩阵选定唯一一个 Q
四元数对角和轴缺陷的修复(本教程绘):左上——半角打包使 ( ϑ , r ) (\vartheta,\boldsymbol r) ( ϑ , r ) 与 ( − ϑ , − r ) (-\vartheta,-\boldsymbol r) ( − ϑ , − r ) 映到同一个 Q Q Q (定理 2.5.A 的不唯一性消失);左下——ϑ ∼ ϑ + 2 π \vartheta\sim\vartheta+2\pi ϑ ∼ ϑ + 2 π 仍给出 Q ∼ − Q Q\sim-Q Q ∼ − Q 的双覆盖,η ≥ 0 \eta\ge0 η ≥ 0 约定把它也裁决掉。
给定四元数,旋转矩阵为(原书式 2.33)
R ( η , ϵ ) = [ 2 ( η 2 + ϵ x 2 ) − 1 2 ( ϵ x ϵ y − η ϵ z ) 2 ( ϵ x ϵ z + η ϵ y ) 2 ( ϵ x ϵ y + η ϵ z ) 2 ( η 2 + ϵ y 2 ) − 1 2 ( ϵ y ϵ z − η ϵ x ) 2 ( ϵ x ϵ z − η ϵ y ) 2 ( ϵ y ϵ z + η ϵ x ) 2 ( η 2 + ϵ z 2 ) − 1 ] R(\eta,\boldsymbol\epsilon)=\begin{bmatrix} 2(\eta^2+\epsilon_x^2)-1 & 2(\epsilon_x\epsilon_y-\eta\epsilon_z) & 2(\epsilon_x\epsilon_z+\eta\epsilon_y)\ 2(\epsilon_x\epsilon_y+\eta\epsilon_z) & 2(\eta^2+\epsilon_y^2)-1 & 2(\epsilon_y\epsilon_z-\eta\epsilon_x)\ 2(\epsilon_x\epsilon_z-\eta\epsilon_y) & 2(\epsilon_y\epsilon_z+\eta\epsilon_x) & 2(\eta^2+\epsilon_z^2)-1 \end{bmatrix} R ( η , ϵ ) = 2 ( η 2 + ϵ x 2 ) − 1 2 ( ϵ x ϵ y + η ϵ z ) 2 ( ϵ x ϵ z − η ϵ y ) 2 ( ϵ x ϵ y − η ϵ z ) 2 ( η 2 + ϵ y 2 ) − 1 2 ( ϵ y ϵ z + η ϵ x ) 2 ( ϵ x ϵ z + η ϵ y ) 2 ( ϵ y ϵ z − η ϵ x ) 2 ( η 2 + ϵ z 2 ) − 1
证明 · 式 2.33 是 Rodrigues 紧凑形式的半角改写(本编补充)
把 η = cos ϑ 2 \eta=\cos\frac\vartheta2 η = cos 2 ϑ 、ϵ = sin ϑ 2 r \boldsymbol\epsilon=\sin\frac\vartheta2\,\boldsymbol r ϵ = sin 2 ϑ r 代入 2.5.1 的紧凑形式 R = c ϑ I + ( 1 − c ϑ ) r r T + s ϑ S ( r ) R=c_\vartheta I+(1-c_\vartheta)\boldsymbol r\boldsymbol r^{\mathrm T}+s_\vartheta S(\boldsymbol r) R = c ϑ I + ( 1 − c ϑ ) r r T + s ϑ S ( r ) 。用倍角公式 c ϑ = 2 η 2 − 1 = 1 − 2 sin 2 ϑ 2 c_\vartheta=2\eta^2-1=1-2\sin^2\frac\vartheta2 c ϑ = 2 η 2 − 1 = 1 − 2 sin 2 2 ϑ 、s ϑ = 2 sin ϑ 2 cos ϑ 2 s_\vartheta=2\sin\frac\vartheta2\cos\frac\vartheta2 s ϑ = 2 sin 2 ϑ cos 2 ϑ ,三项逐个改写:
c ϑ I = ( 2 η 2 − 1 ) I , ( 1 − c ϑ ) r r T = 2 sin 2 ϑ 2 r r T = 2 ϵ ϵ T , s ϑ S ( r ) = 2 sin ϑ 2 cos ϑ 2 S ( r ) = 2 η S ( ϵ ) . c_\vartheta I=(2\eta^2-1)I,\qquad (1-c_\vartheta)\boldsymbol r\boldsymbol r^{\mathrm T}=2\sin^2\tfrac\vartheta2\ \boldsymbol r\boldsymbol r^{\mathrm T}=2\,\boldsymbol\epsilon\,\boldsymbol\epsilon^{\mathrm T},\qquad s_\vartheta S(\boldsymbol r)=2\sin\tfrac\vartheta2\cos\tfrac\vartheta2\ S(\boldsymbol r)=2\eta\,S(\boldsymbol\epsilon). c ϑ I = ( 2 η 2 − 1 ) I , ( 1 − c ϑ ) r r T = 2 sin 2 2 ϑ r r T = 2 ϵ ϵ T , s ϑ S ( r ) = 2 sin 2 ϑ cos 2 ϑ S ( r ) = 2 η S ( ϵ ) .
(S ( ϵ ) = sin ϑ 2 S ( r ) S(\boldsymbol\epsilon)=\sin\frac\vartheta2\,S(\boldsymbol r) S ( ϵ ) = sin 2 ϑ S ( r ) ——S ( ⋅ ) S(\cdot) S ( ⋅ ) 对参数线性。)合并得四元数紧凑形式
R = ( 2 η 2 − 1 ) I + 2 ϵ ϵ T + 2 η S ( ϵ ) , R=(2\eta^2-1)I+2\,\boldsymbol\epsilon\,\boldsymbol\epsilon^{\mathrm T}+2\eta\,S(\boldsymbol\epsilon), R = ( 2 η 2 − 1 ) I + 2 ϵ ϵ T + 2 η S ( ϵ ) ,
逐元素展开:对角线 = 2 η 2 − 1 + 2 ϵ i 2 = 2 ( η 2 + ϵ i 2 ) − 1 =2\eta^2-1+2\epsilon_i^2=2(\eta^2+\epsilon_i^2)-1 = 2 η 2 − 1 + 2 ϵ i 2 = 2 ( η 2 + ϵ i 2 ) − 1 ;非对角 ( i , j ) = 2 ϵ i ϵ j + 2 η S ( ϵ ) i j (i,j)=2\epsilon_i\epsilon_j+2\eta\,S(\boldsymbol\epsilon)_{ij} ( i , j ) = 2 ϵ i ϵ j + 2 η S ( ϵ ) ij ,而 S ( ϵ ) x y = − ϵ z S(\boldsymbol\epsilon)_{xy}=-\epsilon_z S ( ϵ ) x y = − ϵ z 等——即式 2.33 的九个元素。□
□
逆问题(由 R = [ r i j ] R=[r_{ij}] R = [ r ij ] 提取四元数,原书式 2.34–2.35)——取 η ≥ 0 \eta\ge0 η ≥ 0 (对应 ϑ ∈ [ − π , π ] \vartheta\in[-\pi,\pi] ϑ ∈ [ − π , π ] ,可覆盖所有旋转):
η = 1 2 r 11 + r 22 + r 33 + 1 , ϵ = 1 2 [ sgn ( r 32 − r 23 ) r 11 − r 22 − r 33 + 1 sgn ( r 13 − r 31 ) r 22 − r 33 − r 11 + 1 sgn ( r 21 − r 12 ) r 33 − r 11 − r 22 + 1 ] \eta=\tfrac12\sqrt{r_{11}+r_{22}+r_{33}+1},\qquad \boldsymbol\epsilon=\tfrac12\begin{bmatrix}\operatorname{sgn}(r_{32}-r_{23})\sqrt{r_{11}-r_{22}-r_{33}+1}\ \operatorname{sgn}(r_{13}-r_{31})\sqrt{r_{22}-r_{33}-r_{11}+1}\ \operatorname{sgn}(r_{21}-r_{12})\sqrt{r_{33}-r_{11}-r_{22}+1}\end{bmatrix} η = 2 1 r 11 + r 22 + r 33 + 1 , ϵ = 2 1 sgn ( r 32 − r 23 ) r 11 − r 22 − r 33 + 1 sgn ( r 13 − r 31 ) r 22 − r 33 − r 11 + 1 sgn ( r 21 − r 12 ) r 33 − r 11 − r 22 + 1
(sgn ( x ) = 1 ( x ≥ 0 ) \operatorname{sgn}(x)=1\ (x\ge0) sgn ( x ) = 1 ( x ≥ 0 ) 、− 1 ( x < 0 ) -1\ (x<0) − 1 ( x < 0 ) )
证明 · 式 2.34–2.35 的三条线索(本编补充)
求 η \eta η ——取迹。 对四元数紧凑形式取迹(S ( ϵ ) S(\boldsymbol\epsilon) S ( ϵ ) 对角为零、tr ( ϵ ϵ T ) = ∣ ϵ ∣ 2 = 1 − η 2 \operatorname{tr}(\boldsymbol\epsilon\boldsymbol\epsilon^{\mathrm T})=|\boldsymbol\epsilon|^2=1-\eta^2 tr ( ϵ ϵ T ) = ∣ ϵ ∣ 2 = 1 − η 2 ):
tr R = 3 ( 2 η 2 − 1 ) + 2 ∣ ϵ ∣ 2 = 6 η 2 − 3 + 2 ( 1 − η 2 ) = 4 η 2 − 1 ⟹ η = 1 2 tr R + 1 \operatorname{tr}R=3(2\eta^2-1)+2|\boldsymbol\epsilon|^2=6\eta^2-3+2(1-\eta^2)=4\eta^2-1 \ \Longrightarrow\ \eta=\tfrac12\sqrt{\operatorname{tr}R+1} tr R = 3 ( 2 η 2 − 1 ) + 2∣ ϵ ∣ 2 = 6 η 2 − 3 + 2 ( 1 − η 2 ) = 4 η 2 − 1 ⟹ η = 2 1 tr R + 1
(开方取正根,即 η ≥ 0 \eta\ge0 η ≥ 0 约定。)
求 ∣ ϵ i ∣ |\epsilon_i| ∣ ϵ i ∣ ——用对角线。 对角元素 r i i = 2 η 2 − 1 + 2 ϵ i 2 r_{ii}=2\eta^2-1+2\epsilon_i^2 r ii = 2 η 2 − 1 + 2 ϵ i 2 解出 ϵ i 2 = r i i + 1 2 − η 2 = r i i − r j j − r k k + 1 4 \epsilon_i^2=\dfrac{r_{ii}+1}{2}-\eta^2=\dfrac{r_{ii}-r_{jj}-r_{kk}+1}{4} ϵ i 2 = 2 r ii + 1 − η 2 = 4 r ii − r jj − r kk + 1 (i , j , k i,j,k i , j , k 互异)——即式 2.35 根号内的三个组合。
定符号——用反对称部分。 紧凑形式中前两项对称、末项反对称,故 R − R T = 4 η S ( ϵ ) R-R^{\mathrm T}=4\eta\,S(\boldsymbol\epsilon) R − R T = 4 η S ( ϵ ) :ϵ x \epsilon_x ϵ x 的符号 = sgn ( r 32 − r 23 ) =\operatorname{sgn}(r_{32}-r_{23}) = sgn ( r 32 − r 23 ) (η > 0 \eta>0 η > 0 时直接读出;η = 0 \eta=0 η = 0 即 ϑ = π \vartheta=\pi ϑ = π 的边界由 sgn \operatorname{sgn} sgn 规则约定)。□
□
逆四元数由转置关系给出(原书式 2.36):
Q − 1 = { η , − ϵ } ⟺ R − 1 = R T . Q^{-1}={\eta,-\boldsymbol\epsilon}\quad\Longleftrightarrow\quad R^{-1}=R^T . Q − 1 = { η , − ϵ } ⟺ R − 1 = R T .
(把 { η , − ϵ } {\eta,-\boldsymbol\epsilon} { η , − ϵ } 代入式 2.33:ϵ → − ϵ \boldsymbol\epsilon\to-\boldsymbol\epsilon ϵ → − ϵ 使反对称件反号、对称件不变,恰好逐元素转置。)
两个旋转依次施加(R 1 R_1 R 1 后 R 2 R_2 R 2 )对应四元数乘积(原书式 2.37,⊗ \otimes ⊗ 为四元数乘算子):
Q 1 ⊗ Q 2 = { η 1 η 2 − ϵ 1 T ϵ 2 , η 1 ϵ 2 + η 2 ϵ 1 + ϵ 1 × ϵ 2 } . Q_1\otimes Q_2={\eta_1\eta_2-\boldsymbol\epsilon_1^T\boldsymbol\epsilon_2,\ \ \eta_1\boldsymbol\epsilon_2+\eta_2\boldsymbol\epsilon_1+\boldsymbol\epsilon_1\times\boldsymbol\epsilon_2}. Q 1 ⊗ Q 2 = { η 1 η 2 − ϵ 1 T ϵ 2 , η 1 ϵ 2 + η 2 ϵ 1 + ϵ 1 × ϵ 2 } .
乘法公式从哪来——「标量 + 向量」的数(本编补充)
四元数是形如 η + ϵ x i + ϵ y j + ϵ z k \eta+\epsilon_x i+\epsilon_y j+\epsilon_z k η + ϵ x i + ϵ y j + ϵ z k 的数(Hamilton 1843):标量 + 向量 ,虚单位满足 i 2 = j 2 = k 2 = i j k = − 1 i^2=j^2=k^2=ijk=-1 i 2 = j 2 = k 2 = ijk = − 1 (等价地 i j = k , j k = i , k i = j ij=k,\ jk=i,\ ki=j ij = k , jk = i , ki = j ,反序差一负号)。按分配律逐项相乘,再把标量部分与向量部分分别收集,乘法规则就浓缩成一条:
( a , u ) ⊗ ( b , v ) = ( a b − u ⋅ v , a v + b u + u × v ) . (a,\boldsymbol u)\otimes(b,\boldsymbol v)=\big(ab-\boldsymbol u\cdot\boldsymbol v,\ \ a\boldsymbol v+b\boldsymbol u+\boldsymbol u\times\boldsymbol v\big). ( a , u ) ⊗ ( b , v ) = ( ab − u ⋅ v , a v + b u + u × v ) .
式 2.37 就是它在 ( η 1 , ϵ 1 ) (\eta_1,\boldsymbol\epsilon_1) ( η 1 , ϵ 1 ) 、( η 2 , ϵ 2 ) (\eta_2,\boldsymbol\epsilon_2) ( η 2 , ϵ 2 ) 上的展开。两个直接推论:乘法含叉积项 → 不可交换 (Q 1 ⊗ Q 2 ≠ Q 2 ⊗ Q 1 Q_1\otimes Q_2\neq Q_2\otimes Q_1 Q 1 ⊗ Q 2 = Q 2 ⊗ Q 1 ,与旋转合成不可交换同源,2.3 );「共轭」Q ˉ = { η , − ϵ } \bar Q=\{\eta,-\boldsymbol\epsilon\} Q ˉ = { η , − ϵ } 满足 Q ⊗ Q ˉ = { η 2 + ∣ ϵ ∣ 2 , 0 } = { 1 , 0 } Q\otimes\bar Q=\{\eta^2+|\boldsymbol\epsilon|^2,\boldsymbol 0\}=\{1,\boldsymbol 0\} Q ⊗ Q ˉ = { η 2 + ∣ ϵ ∣ 2 , 0 } = { 1 , 0 } ——单位四元数的共轭就是逆 ,即式 2.36。
用 2.5.3 的语言读四元数(本编补充)
四元数把指数映射「搬到一半角上」:对纯向量四元数 ( 0 , u ) (0,\boldsymbol u) ( 0 , u ) 定义指数(泰勒级数,与 定理 2.5.D 同法,用 u × ( u × ⋅ ) \boldsymbol u\times(\boldsymbol u\times\cdot) u × ( u × ⋅ ) 的三循环可证)
exp ( ( 0 , u ) ) = ( cos ∣ u ∣ , sin ∣ u ∣ u ∣ u ∣ ) , \exp\big((0,\boldsymbol u)\big)=\Big(\cos|\boldsymbol u|,\ \ \sin|\boldsymbol u|\,\frac{\boldsymbol u}{|\boldsymbol u|}\Big), exp ( ( 0 , u ) ) = ( cos ∣ u ∣ , sin ∣ u ∣ ∣ u ∣ u ) ,
则角和轴表示正好是 Q = exp ( ( 0 , ϑ 2 r ) ) = ( cos ϑ 2 , sin ϑ 2 r ) Q=\exp\big((0,\tfrac\vartheta2\boldsymbol r)\big)=\big(\cos\frac\vartheta2,\sin\frac\vartheta2\boldsymbol r\big) Q = exp ( ( 0 , 2 ϑ r ) ) = ( cos 2 ϑ , sin 2 ϑ r ) 。取了一半的角 ,所以 ϑ → ϑ + 2 π \vartheta\to\vartheta+2\pi ϑ → ϑ + 2 π (S O ( 3 ) SO(3) SO ( 3 ) 上仍是同一旋转)时 Q → − Q Q\to-Q Q → − Q ——这正是双覆盖 Q ∼ − Q Q\sim-Q Q ∼ − Q 的代数根源,也是式 2.33 里处处「半角」的来历。
为什么控制领域偏爱四元数
单位四元数满足的单位约束把参数限制在一个光滑曲面上 (R 4 \mathbb{R}^4 R 4 中的单位球面),姿态随时间的演化因而是球面上的一条连续轨迹——既不撞奇点、也没有「同一姿态两套参数」的跳变,插值与数值积分都稳定。代价:4 个参数不再能直观读成角度,且乘法含叉积项(定理 2.2.B 的「三重身份」在四元数里退化为纯粹的代数运算)。欧拉角给人看,四元数给机器算 ——这是姿态表示选型最实用的一句总结。
默写点 · 合上书再写(主动回忆)
不看正文与速查卡,先写出下面几条,再展开核对。写不出来就说明这一节还没「拿走」。
写出 ( ϑ , r ) → ( η , ϵ ) (\vartheta,\boldsymbol r)\to(\eta,\boldsymbol\epsilon) ( ϑ , r ) → ( η , ϵ ) 的换算与单位约束。
四元数的独立参数是几个?为什么?
四元数修掉了角和轴的哪两个缺陷?
核对
η = cos ϑ 2 \eta=\cos\frac{\vartheta}{2} η = cos 2 ϑ 、ϵ = sin ϑ 2 r \boldsymbol\epsilon=\sin\frac{\vartheta}{2}\boldsymbol r ϵ = sin 2 ϑ r ;约束 η 2 + ϵ x 2 + ϵ y 2 + ϵ z 2 = 1 \eta^2+\epsilon_x^2+\epsilon_y^2+\epsilon_z^2=1 η 2 + ϵ x 2 + ϵ y 2 + ϵ z 2 = 1 。
3 个(4 参数 − 1 约束),与最简表示的要求不矛盾。
① 表示不唯一 :半角打包使 ( ϑ , r ) (\vartheta,\boldsymbol r) ( ϑ , r ) 与 ( − ϑ , − r ) (-\vartheta,-\boldsymbol r) ( − ϑ , − r ) 映到同一四元数;② 逆解奇点 :式 2.34–2.35 对任何 R R R 都有效,没有 sin ϑ = 0 \sin\vartheta=0 sin ϑ = 0 那样的失效点。
一条主线 :半角打包 ( ϑ , r ) ↦ { cos ϑ 2 , sin ϑ 2 r } (\vartheta,\boldsymbol r)\mapsto{\cos\frac{\vartheta}{2},\sin\frac{\vartheta}{2}\boldsymbol r} ( ϑ , r ) ↦ { cos 2 ϑ , sin 2 ϑ r } → 单位约束(4 参数 − 1 = 3)→ 正解式 2.33 / 逆解式 2.34–2.35 → 逆 Q − 1 = { η , − ϵ } Q^{-1}={\eta,-\boldsymbol\epsilon} Q − 1 = { η , − ϵ } 、乘 Q 1 ⊗ Q 2 Q_1\otimes Q_2 Q 1 ⊗ Q 2 。
三个修复 :表示不唯一(定理 2.5.A ,半角打包使 ( ϑ , r ) (\vartheta,\boldsymbol r) ( ϑ , r ) 与 ( − ϑ , − r ) (-\vartheta,-\boldsymbol r) ( − ϑ , − r ) 映到同一四元数)、Q ∼ − Q Q\sim-Q Q ∼ − Q 双重覆盖(源自 ϑ ∼ ϑ + 2 π \vartheta\sim\vartheta+2\pi ϑ ∼ ϑ + 2 π ,由 η ≥ 0 \eta\ge0 η ≥ 0 约定消除)、逆解奇点(定理 2.5.B )。
一个取舍 :4 参数非最简 × 无奇点 × 乘法含叉积 = 控制 / 图形学首选,人机交互界面仍常用欧拉角。
常见错误清单 · 自查两条
忘记单位约束——四元数必须归一化才能当旋转用,长期数值迭代后需重新归一化;
把「Q Q Q 与 − Q -Q − Q 同旋转」当成角和轴式的不唯一。前者源自 ϑ \vartheta ϑ 与 ϑ + 2 π \vartheta+2\pi ϑ + 2 π 等价(同一旋转、两个四元数),但逆解的 η ≥ 0 \eta\ge0 η ≥ 0 约定已为每个 R R R 选定唯一一个 Q Q Q ——锁定符号是人为却合法的约定;角和轴的不唯一则不同,它是两组参数 都留在表示里、无人裁决。