二体问题之3：航天器轨道方程

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

1. 轨道方程

轨道方程可由积分常数拉普拉斯积分推导和定义：
$$\vec{L}=\vec{v}\times \vec{h}-\frac{\mu }{r}\vec{r}$$
在右侧点乘$\mathbf{r}$，已知二者都位于轨道面内，设夹角为$\theta$。
$$\begin{array}{rl}\vec{L}\cdot \vec{r}& =(\vec{v}\times \vec{h})\vec{r}-\frac{\mu }{r}\vec{r}\cdot \vec{r}\\ Lr\cos \theta & =-(\vec{h}\times \vec{v})\vec{r}-\frac{\mu }{r}{r}^2\\ & =-\vec{h}\cdot (\vec{v}\times \vec{r})-\mu r\\ & =h^2-\mu r\end{array}$$
整理可得矢径$\mathbf{r}$的大小$r$与真近点角$\theta$的关系
$$r=\frac{h^2}{\mu +L\cos \theta }=\frac{\frac{h^2}{\mu }}{1+\frac{L}{\mu }\cos \theta }$$
已知圆锥曲线方程为
$$r=\frac{p}{1+e\cos \theta }$$

则可定义极径$p=\frac{h^2}{\mu }$，偏心率$e=\frac{L}{\mu }$且$e^2={\left(\frac{\vec{L}}{\mu }\right)}^2=\frac{{\mu }^2+2E{h}^2}{{\mu }^2}=1+\frac{2E{h}^2}{{\mu }^2}$。此时可以给定2个积分常数（角动量和拉普拉斯常量）的几何表示。

2. 特殊几何参数

近地点p和远地点a处速度和位置矢量垂直，二者方向和大小为：
$$\{\begin{array}{l}\theta =0,r_p=\frac{P}{1+e},v_p=\frac{h}{r_p}=\sqrt{\frac{\mu }{p}}(1+e)\\ \theta =\pi ,r_a=\frac{p}{1-e},v_a=\frac{h}{r_a}=\sqrt{\frac{\mu }{p}}(1-e)\end{array}$$
其中，近地点p处$\theta$为0，为拉普拉斯常量相同方向。

进一步可以推导出，a为半长轴，p为极径，其关系如下
$$a=\frac{r_a+r_p}{2}=\frac{p}{2}\left(\frac{1-e+1+e}{1-e^2}\right)=\frac{P}{1-e^2}$$

将特殊点的位置大小和速度大小带入积分常数的能量积分可以得到其几何表示。

3. 圆与椭圆

当某一个圆半径与一个椭圆半长轴相等时，其相同x对应的y有如下关系
$$ye=\frac{b}{a}{y}_c$$
将圆与椭圆的方程变形可得。
$$\{\begin{array}{ll}\frac{x^2}{a^2}+\frac{y_2^2}{a^2}=1,& y_c=\sqrt{a^2-x^2}\\ \frac{x^2}{a^2}+\frac{y_p^2}{b^2}=1,& y_p=\frac{b}{a}\sqrt{a^2-x^2}\end{array}$$