20201205 旋转矩阵导数的推导过程
本文不講旋轉(zhuǎn)矩陣導(dǎo)數(shù)的證明,直接講其中一種推導(dǎo)過程。
對象:姿態(tài)旋轉(zhuǎn)矩陣
坐標(biāo)系定義:
歐拉旋轉(zhuǎn)定理:
FB\mathcal F_{B}FB? 相對于 FR\mathcal F_{R}FR? 的旋轉(zhuǎn)可以表示成繞某一個單位軸 e\boldsymbol ee 旋轉(zhuǎn) φ\varphiφ
相關(guān)定義:
旋轉(zhuǎn)矩陣 RRR:從 FR\mathcal F_{R}FR? 到 FB\mathcal F_{B}FB? 的姿態(tài)旋轉(zhuǎn)矩陣,其大小為R=cos?φI3+(1?cos?φ)eeT?sin?φe×R = \cos \varphi I_{3}+\left(1-\cos \varphi \right) \boldsymbol e {\boldsymbol e}^{\mathrm{T}}-\sin \varphi \boldsymbol e^{\times}R=cosφI3?+(1?cosφ)eeT?sinφe×。
旋轉(zhuǎn)角速度 ω\boldsymbol \omegaω:FB\mathcal F_{B}FB? 相對于 FR\mathcal F_{R}FR? 的旋轉(zhuǎn)角速度在 FB\mathcal F_{B}FB? 上的投影,其大小為 ω=ωe\boldsymbol \omega = \omega \boldsymbol eω=ωe, 其中 ω\omegaω 為 ω\boldsymbol \omegaω 的轉(zhuǎn)速大小。
e×=[0?e3e2e30?e1?e2e10]\boldsymbol{e}^{\times}=\left[\begin{array}{ccc}0 & -e_{3} & e_{2} \\ e_{3} & 0 & -e_{1} \\ -e_{2} & e_{1} & 0\end{array}\right]e×=???0e3??e2???e3?0e1??e2??e1?0????。
旋轉(zhuǎn)矩陣導(dǎo)數(shù)推導(dǎo)過程:
在 ttt 時刻,存在極小的時間段 Δt\Delta tΔt 的旋轉(zhuǎn)可以滿足:
R(t+Δt)=R(Δt)R(t)\boldsymbol R(t+\Delta t) = \boldsymbol R(\Delta t) \boldsymbol R(t)R(t+Δt)=R(Δt)R(t)
其中,根據(jù)旋轉(zhuǎn)矩陣的定義,可以得:
R(Δt)=cos?ΔφI3+(1?cos?Δφ)eeT?sin?Δφe×\boldsymbol R(\Delta t) = \cos \Delta \varphi I_{3}+\left(1-\cos \Delta \varphi \right) \boldsymbol e {\boldsymbol e}^{\mathrm{T}}-\sin \Delta \varphi \boldsymbol e^{\times}R(Δt)=cosΔφI3?+(1?cosΔφ)eeT?sinΔφe×
根據(jù)極限的思想,在極小的 Δt\Delta tΔt 中產(chǎn)生的 Δφ\Delta \varphiΔφ 也是極小的,那么也就是說
cos?Δφ=1andsin?Δφ=Δφ\cos \Delta \varphi = 1 \quad \text{and} \quad \sin \Delta \varphi = \Delta \varphicosΔφ=1andsinΔφ=Δφ
將其帶入 R(Δt)\boldsymbol R(\Delta t)R(Δt) 的表達(dá)式,可以得到
R(Δt)=I3?Δφe×\boldsymbol R(\Delta t) = I_{3}- \Delta \varphi \boldsymbol e^{\times}R(Δt)=I3??Δφe×
dRdt=lim?Δt→0R(t+Δt)?R(t)Δt=lim?Δt→0R(t)?Δφe×R(t)?R(t)Δt=lim?Δt→0?Δφe×R(t)Δt=lim?Δt→0?ωe×R(t)=?ω×R\begin{aligned} \frac{d \boldsymbol R}{d t}&=\lim _{\Delta t \rightarrow 0} \frac{\boldsymbol R(t+\Delta t)-\boldsymbol R(t)}{\Delta t} \\&=\lim _{\Delta t \rightarrow 0} \frac{\boldsymbol R(t)- \Delta \varphi \boldsymbol e^{\times} \boldsymbol R(t)-\boldsymbol R(t)}{\Delta t} \\&=\lim _{\Delta t \rightarrow 0} \frac{- \Delta \varphi \boldsymbol e^{\times} \boldsymbol R(t)}{\Delta t} \\&=\lim _{\Delta t \rightarrow 0} {- \omega \boldsymbol e^{\times} \boldsymbol R(t)} \\&= {- \boldsymbol \omega^{\times} \boldsymbol R} \end{aligned}dtdR??=Δt→0lim?ΔtR(t+Δt)?R(t)?=Δt→0lim?ΔtR(t)?Δφe×R(t)?R(t)?=Δt→0lim?Δt?Δφe×R(t)?=Δt→0lim??ωe×R(t)=?ω×R?
結(jié)論:R˙=?ω×R\dot {\boldsymbol R} = - \boldsymbol \omega^{\times} \boldsymbol RR˙=?ω×R
參考文獻(xiàn):章仁為. 衛(wèi)星軌道姿態(tài)動力學(xué)與控制. 5.3節(jié)姿態(tài)運(yùn)動學(xué)方程
總結(jié)
以上是生活随笔為你收集整理的20201205 旋转矩阵导数的推导过程的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 2020-12-03 The Geome
- 下一篇: 2020-12-15 IEEE写作格式