Given parameters |
Symbol |
radial distance vector |
$\vec{r}$ |
veloicity vector |
$\vec{v}$ |
standard gravitational parameter |
$\mu$ |
semi-major axis |
a |
1. Calculate the specific angular momentum vector, $\vec{h}$
Specific angular momentum is the cross product of the radial distance and velocity vectors, and is useful for finding many of the other Keplerian paramenters, both when modeling 2d and 3d orbits.
$$\vec{h} = \vec{r} \times \vec{v}$$
2. Calculate the eccentricity vector, $\vec{e}$
The eccentricity vector is a vector with the magnitude of the orbital eccentricity, that points from the center of the object being orbited, towards the pericenter.
$$\vec{e}=\frac{\vec{v} \times\vec{h}}{\mu} - \frac{\vec{r}}{|\vec{r}|}$$
It can also be calculated without using the angular momentum vector:
$$\vec{e} = \left(\frac{|\vec{v}|^2}{\mu} - \frac{1}{|\vec{r}|}\right)\vec{r} - \frac{\vec{r} \cdot\vec{v}}{\mu} \vec{v}$$
3. Calculate the pericenter distance, $q$
For elliptical orbits ($e<=1 ,a > 0$), and hyperbolic trajectories ($ e > 1 ,a < 0$)
$$q = a(1-|\vec{e}|)$$
For parabolic trajectories ($e = 1, a = \pm \infty$):
$$q = \frac{|\vec{h}|^2}{2\mu}$$
Circular orbits don't have a well-defined pericenter.
4. Calculate the pericenter position vector, $\vec{q}$
The pericenter position vector, $\vec{q}$, has magnitude of the pericenter distance $q$, in the direction of the eccentricity vector.
$$\vec{q} = q \frac{\vec{e}}{|\vec{e}|}$$