math lessons - Eccentric anomaly

The eccentric anomaly is the angle between the direction of periapsis and the current position of an object on its orbit, projected onto the ellipse's circumscribing circle perpendicularly to the major axis, measured at the centre of the ellipse. In the diagram below, it is E (the angle zcx).

Calculation

In astrodynamics eccentric anomaly E can be calculated as follows:

$E=\arccos {{1-\left [ \mathbf{r} \right ] / a} \over e}$

where:

$\mathbf{r}\,\!$ is the orbiting body's position vector (segment sp),
$a\,\!$ is the orbit's semi-major axis (segment cz), and
$e\,\!$ is the orbit's eccentricity.

The relation between E and M, the mean anomaly, is:

$M = E - e \cdot \sin{E}.\,\!$

For small values of $e$ ( $e < 0.6627434$ ) this equation can be solved iteratively, starting from $E 0 = M$ and using the relation $E_{i+1} = M + e\,\sin E_i$ . The first few terms of the expansion in powers of $e$ are:

$E_1 = M + e\,\sin M$
$E_2 = M + e\,\sin M + \frac{1}{2} e^2 \sin 2M$
$E_3 = M + e\,\sin M + \frac{1}{2} e^2 \sin 2M + \frac{1}{8} e^3 (3\sin 3M - \sin M)$ .