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

1. 轨道方程

轨道方程可由积分常数拉普拉斯积分推导和定义:
L ⃗ = v ⃗ × h ⃗ − μ r r ⃗ \vec{L}=\vec{v} \times \vec{h}-\frac{\mu}{r} \vec{r} L =v ×h −rμ​r
在右侧点乘 r \mathbf{r} r,已知二者都位于轨道面内,设夹角为 θ \theta θ。
L ⃗ ⋅ r ⃗ = ( v ⃗ × h ⃗ ) r ⃗ − μ r r ⃗ ⋅ r ⃗ L r cos ⁡ θ = − ( h ⃗ × v ⃗ ) r ⃗ − μ r r 2 = − h ⃗ ⋅ ( v ⃗ × r ⃗ ) − μ r = h 2 − μ r \begin{aligned} \vec{L} \cdot \vec{r} &=(\vec{v} \times \vec{h}) \vec{r}-\frac{\mu}{r} \vec{r} \cdot \vec{r} \\ L r \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{aligned} L ⋅r Lrcosθ​=(v ×h )r −rμ​r ⋅r =−(h ×v )r −rμ​r2=−h ⋅(v ×r )−μr=h2−μr​
整理可得矢径 r \mathbf{r} r的大小 r r r与真近点角 θ \theta θ的关系
r = h 2 μ + L cos ⁡ θ = h 2 μ 1 + L μ cos ⁡ θ r=\frac{h^{2}}{\mu+L \cos \theta}=\frac{\frac{h^{2}}{\mu}}{1+\frac{L}{\mu} \cos \theta} r=μ+Lcosθh2​=1+μL​cosθμh2​​
已知圆锥曲线方程为
r = p 1 + e cos ⁡ θ r=\frac{p}{1+e \cos \theta} r=1+ecosθp​
二体问题之3:航天器轨道方程

则可定义极径 p = h 2 μ p=\frac{h^{2}}{\mu} p=μh2​,偏心率 e = L μ e=\frac{L}{\mu} e=μL​且 e 2 = ( L ⃗ μ ) 2 = μ 2 + 2 E h 2 μ 2 = 1 + 2 E h 2 μ 2 e^{2}=\left(\frac{\vec{L}}{\mu}\right)^{2}=\frac{\mu^{2}+2 E h^{2}}{\mu^{2}}=1+\frac{2 E h^{2}}{\mu^{2}} e2=(μL ​)2=μ2μ2+2Eh2​=1+μ22Eh2​。此时可以给定2个[[积分常数]](角动量和拉普拉斯常量)的几何表示。

2. 特殊几何参数

近地点p和远地点a处速度和位置矢量垂直,二者方向和大小为:
{ θ = 0 , r p = P 1 + e , v p = h r p = μ p ( 1 + e ) θ = π , r a = p 1 − e , v a = h r a = μ p ( 1 − e ) \left\{\begin{array}{l} \theta=0, r_{p}=\frac{P}{1+e}, \quad v_{p}=\frac{h}{r_{p}}=\sqrt{\frac{\mu}{p}}(1+e) \\ \theta=\pi, r_{a}=\frac{p}{1-e}, \quad v_{a}=\frac{h}{r_{a}}=\sqrt{\frac{\mu}{p}}(1-e) \end{array}\right. ⎩⎨⎧​θ=0,rp​=1+eP​,vp​=rp​h​=pμ​ ​(1+e)θ=π,ra​=1−ep​,va​=ra​h​=pμ​ ​(1−e)​
其中,近地点p处 θ \theta θ为0,为拉普拉斯常量相同方向。

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

进一步可以推导出,a为半长轴,p为极径,其关系如下
a = r a + r p 2 = p 2 ( 1 − e + 1 + e 1 − e 2 ) = P 1 − e 2 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}} a=2ra​+rp​​=2p​(1−e21−e+1+e​)=1−e2P​

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

3. 圆与椭圆

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

当某一个圆半径与一个椭圆半长轴相等时,其相同x对应的y有如下关系
y e = b a y c y e=\frac{b}{a} y_{c} ye=ab​yc​
{ x 2 a 2 + y 2 2 a 2 = 1 , y c = a 2 − x 2 x 2 a 2 + y p 2 b 2 = 1 , y p = b a a 2 − x 2 \begin{cases}\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{cases} {a2x2​+a2y22​​=1,a2x2​+b2yp2​​=1,​yc​=a2−x2 ​yp​=ab​a2−x2 ​​

上一篇:opencv+python识别猪


下一篇:P3901 数列找不同 题解