016四元数微分方程
果然,我還是沒有看平臺式慣導的初始對準和高度通道,而是直接來到了捷聯慣導的部分。之前看過四元數部分的內容,所以現在就當是復習了。回顧一下我之前做的筆記,好像沒有四元數微分方程,所以還是在這記錄一下吧,也可以說搬運一下。
在捷聯慣導中,表征從導航系n系到機體系b系的四元數為:
Q ? = c o s θ 2 + u ? n s i n θ 2 \vec{Q}=cos\frac{\theta}{2}+\vec{u}^nsin\frac{\theta}{2} Q?=cos2θ?+unsin2θ?
式中,剛體繞瞬時軸 u ? n \vec{u}^n un(單位向量)旋轉。
兩邊求導可得:
d Q ? d t = ? θ ˙ 2 s i n θ 2 + u ? n θ ˙ 2 c o s θ 2 + s i n θ 2 d u ? n d t \frac{d\vec{Q}}{dt}=-\frac{\dot{\theta} } {2}sin\frac{\theta}{2} + \vec{u}^n \frac{ \dot{\theta} }{2} cos\frac{\theta}{2} + sin\frac{\theta}{2} \frac{d\vec{u}^n}{dt} dtdQ??=?2θ˙?sin2θ?+un2θ˙?cos2θ?+sin2θ?dtdun?
可以這樣理解, u ? n \vec{u}^n un為某一瞬時固定軸,是一個常量,所以其導數為0,則:
d Q ? d t = ? θ ˙ 2 s i n θ 2 + u ? n θ ˙ 2 c o s θ 2 \frac{d\vec{Q}}{dt}=-\frac{\dot{\theta} } {2}sin\frac{\theta}{2} + \vec{u}^n \frac{ \dot{\theta} }{2} cos\frac{\theta}{2} dtdQ??=?2θ˙?sin2θ?+un2θ˙?cos2θ?
對于 u ? n θ ˙ 2 \vec{u}^n \frac{ \dot{\theta} }{2} un2θ˙?和四元數 Q ? \vec{Q} Q?的乘積有:
u ? n θ ˙ 2 ? Q ? = u ? n θ ˙ 2 ? ( c o s θ 2 + u ? n θ ˙ 2 ) = u ? n θ ˙ 2 c o s θ 2 + u ? n ? u ? n θ ˙ 2 s i n θ 2 \vec{u}^n \frac{ \dot{\theta} }{2} \otimes \vec{Q} =\vec{u}^n \frac{ \dot{\theta} }{2} \otimes (cos\frac{\theta}{2} + \vec{u}^n \frac{ \dot{\theta} }{2}) = \vec{u}^n \frac{ \dot{\theta} }{2} cos\frac{\theta}{2} + \vec{u}^n \otimes \vec{u}^n \frac{\dot{\theta} } {2}sin\frac{\theta}{2} un2θ˙??Q?=un2θ˙??(cos2θ?+un2θ˙?)=un2θ˙?cos2θ?+un?un2θ˙?sin2θ?
而 u ? n ? u ? n = ? 1 \vec{u}^n \otimes \vec{u}^n = -1 un?un=?1,所以:
u ? n θ ˙ 2 ? Q ? = u ? n θ ˙ 2 c o s θ 2 ? θ ˙ 2 s i n θ 2 = d Q ? d t \vec{u}^n \frac{ \dot{\theta} }{2} \otimes \vec{Q} = \vec{u}^n \frac{ \dot{\theta} }{2} cos\frac{\theta}{2} - \frac{\dot{\theta} } {2}sin\frac{\theta}{2} = \frac{d\vec{Q}}{dt} un2θ˙??Q?=un2θ˙?cos2θ??2θ˙?sin2θ?=dtdQ??
即:
d Q ? d t = θ ˙ 2 u ? n ? Q ? \frac{d\vec{Q}}{dt}=\frac{ \dot{\theta} }{2} \vec{u}^n \otimes \vec{Q} dtdQ??=2θ˙?un?Q?
其中,$
\dot{\theta}\vec{u}^n
$表示角速度標量與角速度方向的乘積即為角速度矢量。這里旋轉方向從n系到b系,在n系上的投影,如下:
θ ˙ u ? n = ω ? n b n \dot{\theta}\vec{u}^n = \vec{\omega}_{nb}^{n} θ˙un=ωnbn?
所以:
d Q ? d t = 1 2 ω ? n b n ? Q ? \frac{d\vec{Q}}{dt}=\frac{1}{2} \vec{\omega}_{nb}^n \otimes \vec{Q} dtdQ??=21?ωnbn??Q?
$ \vec{\omega}{nb}^n $ 需要通過 $ \vec{\omega}{nb}^b $ 求得:
ω ? n b n = Q ? ? ω ? n b b ? Q ? ? \vec{\omega}_{nb}^n = \vec{Q} \otimes \vec{\omega}_{nb}^b \otimes \vec{Q}^* ωnbn?=Q??ωnbb??Q??
即:
d Q ? d t = 1 2 ω ? n b n ? Q ? = 1 2 Q ? ? ω ? n b b ? Q ? ? ? Q ? = 1 2 Q ? ? ω ? n b b \frac{d\vec{Q}}{dt}=\frac{1}{2} \vec{\omega}_{nb}^n \otimes \vec{Q} =\frac{1}{2} \vec{Q} \otimes \vec{\omega}_{nb}^b \otimes \vec{Q}^* \otimes \vec{Q} =\frac{1}{2} \vec{Q} \otimes \vec{\omega}_{nb}^b dtdQ??=21?ωnbn??Q?=21?Q??ωnbb??Q???Q?=21?Q??ωnbb?
d Q ? d t = 1 2 M ′ ( ω ? n b b ) Q ? \frac{d\vec{Q}}{dt}=\frac{1}{2} M'(\vec{\omega}_{nb}^b) \vec{Q} dtdQ??=21?M′(ωnbb?)Q?
這就是四元數微分方程
式中:
ω ? n b b = ω ? i b b ? ω ? i n b = ω ? i b b ? C n b ω ? i n n = ω ? i b b ? C n b ( ω ? i e n + ω ? e n n ) \vec{\omega}_{nb}^b = \vec{\omega}_{ib}^b - \vec{\omega}_{in}^b = \vec{\omega}_{ib}^b - C_n^b \vec{\omega}_{in}^n = \vec{\omega}_{ib}^b - C_n^b( \vec{\omega}_{ie}^n + \vec{\omega}_{en}^n ) ωnbb?=ωibb??ωinb?=ωibb??Cnb?ωinn?=ωibb??Cnb?(ωien?+ωenn?)
ω ? i b b \vec{\omega}_{ib}^b ωibb? 為捷聯陀螺的輸出;
C n b C_n^b Cnb?由姿態更新的最新值確定;
ω ? i e n \vec{\omega}_{ie}^n ωien? 和 ω ? e n n \vec{\omega}_{en}^n ωenn?分別為導航系相對于地球的角速度和地球自轉角速度在導航系上的投影,對于導航系取地理坐標系的情況:
ω ? i e n + ω ? e n n = [ ? V N R M ω i e c o s L + V E R N ω i e s i n L + V E R N t a n L ] \vec{\omega}_{ie}^n + \vec{\omega}_{en}^n = \left[ \begin{matrix} -\frac{V_N}{R_M}\\ \\ \omega_{ie}cosL+\frac{V_E}{R_N}\\ \\ \omega_{ie}sinL+\frac{V_E}{R_N}tanL\\ \end{matrix} \right] ωien?+ωenn?=????????RM?VN??ωie?cosL+RN?VE??ωie?sinL+RN?VE??tanL????????
這在指北方位慣導系統的平臺指令角速度中已經證明過。
感覺最近推導的公式有點多,思緒有時候很混亂,所以打算捋一捋思緒,不知道又要有多長時間。
總結
以上是生活随笔為你收集整理的016四元数微分方程的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: Java中常见的对象类型简述(DO、BO
- 下一篇: 什么是哈希表?什么又是哈希冲突?哈希冲突