基础-四元数基础知识-用于航天器姿态动力学

四元数的历史

推荐观看3blue1brown的视频和交互操作动画，链接:
  四元数的可视化
  交互操作动画
待续。。。
计算坐标系与描述坐标系。
采用四元数进行计算的好处（姿态解算）：

1. 无奇点，方程不会退化；
2. 比方向余弦法需解算9个联立微分方程相比，仅需要4个联立微分方程；
3. 四元数微分方程的约束条件只有1个，比方向余弦法中的6个要少得多。

四元数的运算法则

定义1： 称量$q$
$$q=q_0+q_{1}i+q_{2}j+q_{3}k\text{\hspace{1em}(1)}$$
为四元数，其中$q_0$、$q_1$、$q_2$、$q_3$为任意实数。$q_0$称为四元数的标部，$q_{1}i+q_{2}j+q_{3}k$称为四元数的矢部，并记作$\mathbf{q}$。因此$q$也可以记作：
$$q=q_0+\mathbf{q}\text{\hspace{1em}(2)}$$
定义2： $i$、$j$、$k$满足如下运算法则：
$$\{\begin{array}{rl}{i}^2& =j^2=k^2=-1\\ jk& =-kj=i\\ ki& =-ik=j\\ ij& =-ji=k\end{array}\text{\hspace{1em}(3)}$$
定义3： 称$q^{\ast }$为$q$的共轭四元数：
$$q^{\ast }=q_0-q_{1}i-q_{2}j-q_{3}k=q_0-\mathbf{q}\text{\hspace{1em}(4)}$$
定义4： 设$p$、$q$为任意两个四元数：
$$\begin{gathered} p=p_0+p_{1}i+p_{2}j+p_{3}k \\ q=q_0+q_{1}i+q_{2}j+q_{3}k \end{gathered}$$
则$p$与$q$的相加定义为：
$$\begin{array}{rl}p+q& =(p_0+q_0)+(p_1+q_1)i\\ & +(p_2+q_2)j+(p_3+q_3)k\end{array}\text{\hspace{1em}(5)}$$
称$p+q$ 为四元数$p$与$q$之和。
定义5： 设$p$、$q$为任意两个四元数，则$p$与$q$的相乘定义为：
$$\begin{array}{rl}pq& =(p_0+p_{1}i+p_{2}j+p_{3}k)(q_0+q_{1}i+q_{2}j+q_{3}k)\\ & =p_0{q}_0-(p_1{q}_1+p_2{q}_2+p_3{q}_3)+(p_{1}i+p_{2}j+p_{3}k)q_0\\ & +p_0(q_{1}i+q_{2}j+q_{3}k)+[(p_2{p}_3-p_3{q}_2)i\\ & +(p_3{q}_1-p_1{q}_3)j+(p_1{q}_2-p_2{q}_1)k]\end{array}\text{\hspace{1em}(6)}$$
称$pq$为四元数$p$与$q$的积。

定义6：

定义7：

定义8：

四元数的对数运算

一般情况下，如无特殊说明，我们所指的四元数就是单位四元数。
单位四元数$q=[\eta ,\mathbf{\xi }]$的对数定义为：
$$\ln q=\left[0,\frac{\arccos \eta }{2\sqrt{1-{\eta }^2}}\mathbf{\xi }\right]$$

单位四元数表示旋转，一个绕单位轴$\mathbf{n}$，转过角度$\varphi$的坐标系转动，可以表示为：
$$q=\left[\cos \left(\frac{\varphi }{2}\right),\sin \left(\frac{\varphi }{2}\right)\mathbf{n}\right]$$

对上式所示的使用欧拉轴和角表示的单位四元数，其等价形式为：
$$\ln q=\left[0,\frac{\varphi }{2}\mathbf{n}\right]$$

简写为
$$\ln q=\frac{\varphi }{2}\mathbf{n}$$

定理1： 已知随时间变化的四元数$q(t)$和角速度$\mathbf{\omega }(t)$，有如下等式成立：
$$\frac{1}{2}\frac{\text{d}}{\text{d}t}({\Gamma }^T\Gamma )={\Gamma }^T\mathbf{\omega }$$

其中，$\Gamma =2\ln q$

四元数表示旋转

aaa

四元数运动学方程

欧拉运动学方程是用欧拉角与欧拉角速率表示刚体的角速度。刚体的角速度同样可以用四元数及其导数进行表示，称为四元数运动学方程。

姿态动力学方程

$$\begin{array}{c}{I}_x{\dot{\omega }}_x+(I_z-I_y){\omega }_y{\omega }_z=L_x\\ I_y{\dot{\omega }}_y+(I_x-I_z){\omega }_x{\omega }_z=L_y\\ I_z{\dot{\omega }}_z+(I_y-I_x){\omega }_x{\omega }_y=L_z\end{array}$$

四元数运动学方程

$$\dot{q}=\frac{1}{2}q\mathbf{\omega }$$

用四元数表示的欧拉角

$$\{\begin{array}{c}\varphi =\varphi (q_0,q_1,q_2,q_3)\\ \psi =\psi (q_0,q_1,q_2,q_3)\\ \theta =\theta (q_0,q_1,q_2,q_3)\end{array}\text{\hspace{1em}(xx)}$$
利用上式可消除姿态动力学方程右端的欧拉角，以及四元数方程求解后再变换为欧拉角。

用欧拉角表示的四元数

$$\{\begin{array}{l}{q}_0=q_0(\varphi ,\psi ,\theta )\\ q_1=q_1(\varphi ,\psi ,\theta )\\ q_2=q_2(\varphi ,\psi ,\theta )\\ q_3=q_3(\varphi ,\psi ,\theta )\\ q_0^2+q_1^2+q_2^2+q_3^2=1\end{array}$$
  利用上式可把欧拉角表示的运动初始条件变换为四元数表示的运动初始条件。
  应用上述方程可以把姿态动力学问题变换为以$({\omega }_x,{\omega }_y,{\omega }_z,q_0,q_1,q_2,q_3{)}^T$为状态的动力学方程。

导航计算图

关于万向节死锁（Gimbal Lock）

万向死节锁是由于动态欧拉角导致的。物体绕世界坐标系的轴所旋转的角度就是静态欧拉角；物体绕自身坐标系的轴所旋转的角度就是动态欧拉角。世界坐标系就是我们常说的惯性系。
描述的时候采用的是世界坐标系，旋转的时候用的是物体坐标系，所以导致认知上的错误。