果然,我还是没有看平台式惯导的初始对准和高度通道,而是直接来到了捷联惯导的部分。之前看过四元数部分的内容,所以现在就当是复习了。回顾一下我之前做的笔记,好像没有四元数微分方程,所以还是在这记录一下吧,也可以说搬运一下。
在捷联惯导中,表征从导航系n系到机体系b系的四元数为:
Q⃗ =cosθ2+u⃗ nsinθ2
式中,刚体绕瞬时轴
u⃗ n
(单位向量)旋转。
两边求导可得:
dQ⃗ dt=−θ˙2sinθ2+u⃗ nθ˙2cosθ2+sinθ2du⃗ ndt
可以这样理解,
u⃗ n
为某一瞬时固定轴,是一个常量,所以其导数为0,则:
dQ⃗ dt=−θ˙2sinθ2+u⃗ nθ˙2cosθ2
对于
u⃗ nθ˙2
和四元数
Q⃗
的乘积有:
u⃗ nθ˙2⨂Q⃗ =u⃗ nθ˙2⨂(cosθ2+u⃗ nθ˙2)=u⃗ nθ˙2cosθ2+u⃗ n⨂u⃗ nθ˙2sinθ2
而
u⃗ n⨂u⃗ n=−1
,所以:
u⃗ nθ˙2⨂Q⃗ =u⃗ nθ˙2cosθ2−θ˙2sinθ2=dQ⃗ dt
即:
dQ⃗ dt=θ˙2u⃗ n⨂Q⃗
其中,
θ˙u⃗ n
表示角速度标量与角速度方向的乘积即为角速度矢量。这里旋转方向从n系到b系,在n系上的投影,如下:
θ˙u⃗ n=ω⃗ nnb
所以:
dQ⃗ dt=12ω⃗ nnb⨂Q⃗
ω⃗ nnb
需要通过
ω⃗ bnb
求得:
ω⃗ nnb=Q⃗ ⨂ω⃗ bnb⨂Q⃗ ∗
即:
dQ⃗ dt=12ω⃗ nnb⨂Q⃗ =12Q⃗ ⨂ω⃗ bnb⨂Q⃗ ∗⨂Q⃗ =12Q⃗ ⨂ω⃗ bnb
dQ⃗ dt=12M′(ω⃗ bnb)Q⃗
这就是四元数微分方程
式中:
ω⃗ bnb=ω⃗ bib−ω⃗ bin=ω⃗ bib−Cbnω⃗ nin=ω⃗ bib−Cbn(ω⃗ nie+ω⃗ nen)
ω⃗ bib
为捷联陀螺的输出;
Cbn
由姿态更新的最新值确定;
ω⃗ nie
和
ω⃗ nen
分别为导航系相对于地球的角速度和地球自转角速度在导航系上的投影,对于导航系取地理坐标系的情况:
ω⃗ nie+ω⃗ nen=⎡⎣⎢⎢⎢⎢⎢⎢⎢⎢⎢−VNRMωiecosL+VERNωiesinL+VERNtanL⎤⎦⎥⎥⎥⎥⎥⎥⎥⎥⎥
这在指北方位惯导系统的平台指令角速度中已经证明过。
感觉最近推导的公式有点多,思绪有时候很混乱,所以打算捋一捋思绪,不知道又要有多长时间。