二体问题之6：轨道根数及其转化

注：笔记，恳请批评指正。

1. 轨道坐标系

1. 黄道坐标系
   原点位于太阳质心，轨道面位于黄道面（地球绕太阳运行平面），x轴指向春分点。

2. 地心黄道坐标系
   原点位于地球质心，轨道面位于黄道面（地球绕太阳运行平面），x轴指向春分点。

3. 地心赤道
   原点位于地球质心，轨道面位于赤道面（地球绕太阳运行平面），x轴指向春分点。

4. 地心拱线
   原点位于地球质心，轨道面位于航天器运行轨道面，x轴指向近地点。

5. 地心轨道
   原点位于地球质心，轨道面位于航天器运行轨道面，x轴指向轨道矢径。

2. 轨道根数$\vec{b}$

![[Pasted image 20211025114018.png]]
形状：半长轴、偏心率
面指向：轨道倾角、升交点赤经
拱线指向：近地点辐角

3. 特殊轨道根数

1. 地球赤道轨道
   轨道面位于地球赤道面。升交点赤经和近地点辐角共面。定义近地点经度为赤道面上$I$轴和$\vec{e}$近地点方向夹角
   $${\tilde{\omega }}_{\text{true }}=\Omega +\omega$$
   对于赤道平面轨道的二维情况，近地点$J$分量在负轴的，逆时针减。
   $$\cos \left({\tilde{\omega }}_{\text{the }}\right)=\frac{\hat{I}\cdot \vec{e}}{|\hat{I}||\vec{e}|},e_j<0,{\tilde{\omega }}_{\text{the }}=360^{\circ }-{\tilde{\omega }}_t$$
2. 圆倾角轨道
   不存在近地点，近地点辐角和偏移角无法定义。定义近地点角度为$n$（轨道面与赤道面交线方向）和航天器矢径方向夹角。
   $$\cos (u)=\frac{\vec{n}\cdot \vec{r}}{|\vec{n}|\cdot |\vec{r}|},r_k<0,u=360^{\circ }-u$$
   三维情况，u$K$分量在负轴的，逆时针减。
3. 圆赤道轨道
   地球赤道的延申，只有一个真经度，定义为$I$轴和矢径的夹角。
   $$\cos \left({\lambda }_{\mathrm{tm}}\mathbf{e}\right)=\frac{\hat{I}\vec{r}}{|\hat{I}||\vec{r}|}$$

4. RV2COE

通过计算[[积分常数]]，根据定义推导求解公式。

1. 偏心率
   $$\begin{array}{rl}\vec{e}=\frac{\vec{L}}{\mu }& =\frac{\vec{v}\times \vec{h}-\frac{\mu }{r}\vec{r}}{\mu }\\ & =\frac{\vec{v}((\vec{r}\times \vec{v})}{\mu }\cdot \frac{\vec{r}}{r}\\ & =\frac{(\vec{v}\cdot \vec{v})\vec{r}-(\vec{r}\cdot \vec{v})\vec{v}}{\mu }-\frac{\vec{r}}{r}\\ & =\frac{1}{\mu }\left[\left(v^2\cdot \frac{\mu }{r}\right)\vec{r}-(\vec{r}\cdot \vec{v})\vec{v}\right]\end{array}$$
2. 半长轴
   $$\begin{array}{rl}& \xi =\frac{v^2}{2}-\frac{\mu }{r}\\ & a=-\frac{\mu }{2\xi }\end{array}$$
3. 倾角
   $$\cos (\bar{ı})=\frac{\hat{k}\cdot \vec{h}}{|\hat{k}\Vert \vec{h}|}$$
4. 升交点赤经
   $$\vec{n}=\hat{k}\times \vec{h}$$
   $$\cos (\Omega )=\frac{\hat{I}\cdot \vec{n}}{|\hat{I}||\vec{n}|}{n}_j<0,\Omega =360^{\circ }-\Omega$$
5. 近地点辐角
   $$\cos (\omega )=\frac{\vec{n}\cdot \vec{e}}{|\tilde{n}||\vec{e}|}{e}_k<0,\omega =360^{\circ }-\omega$$
6. 真近点角
   $$\cos (v)=\frac{\vec{e}\cdot \vec{r}}{|\vec{e}||\vec{F}|}\vec{r}\cdot \vec{v}<0,v=360^{\circ }-v$$

5. COE2RV

基本思路是首先用偏移角表示轨道坐标系$PQW$速度和矢径矢量。再通过坐标变换转得到赤道地心坐标系$IJK$。

在$PQW$中
$${\vec{r}}_{PQW}=\left[\begin{array}{c}r\cos (v)\\ r\sin (v)\\ 0\end{array}\right],{\vec{v}}_{\text{PQW}}=\left[\begin{array}{c}r\cos (v)-r\dot{v}\sin (v)\\ \dot{r}\sin (v)+r\dot{c}\cos (v)\\ 0\end{array}\right]$$

根据[[速度描述]]
$$\begin{array}{rl}& r=\frac{p}{1+e\cos (v)}\\ & r\dot{v}=\frac{h}{r}=\frac{\sqrt{\mu P}}{P}(1+e\cos (v))\\ & r=\sqrt{\frac{\mu }{p}}(\mathrm{estn}(v))\end{array}$$

得到
$${\vec{v}}_{PQW}=\left[\begin{array}{c}-\sqrt{\frac{\mu }{P}}\sin (v)\\ \sqrt{\frac{\mu }{p}}(e+\cos (v))\\ 0\end{array}\right]$$

接下来坐标变换

想得到$IJK$中坐标表示，续从右向左将$PQW$绕Z逆时针转近地点辐角，再绕X逆时针转轨道倾角，再绕Z逆时针转升交点赤经。
$$\left[\frac{IJK}{\text{ POW }}\right]=ROT3(-\Omega )ROT1(-i)\mathrm{ROT}3(-\omega )$$
$$\begin{array}{rl}& r_{\text{IJk }}=\left[\frac{IJK}{PQW}\right]r_{PQW}\\ & v_{\text{IJk }}=\left[\frac{IJK}{PQW}\right]v_{\text{PQW }}\end{array}$$