In orbital mechanics, the eccentric anomaly is an angular parameter that defines the position of a body that is moving along an elliptic Kepler orbit. The eccentric anomaly is one of three angular parameters ("anomalies") that define a position along an orbit, the other two being the true anomaly and the mean anomaly.
Consider the ellipse with equation given by:
| x2 | |
| a2 |
+
| y2 | |
| b2 |
=1,
where a is the semi-major axis and b is the semi-minor axis.
For a point on the ellipse, P = P(x, y), representing the position of an orbiting body in an elliptical orbit, the eccentric anomaly is the angle E in the figure. The eccentric anomaly E is one of the angles of a right triangle with one vertex at the center of the ellipse, its adjacent side lying on the major axis, having hypotenuse a (equal to the semi-major axis of the ellipse), and opposite side (perpendicular to the major axis and touching the point P′ on the auxiliary circle of radius a) that passes through the point P. The eccentric anomaly is measured in the same direction as the true anomaly, shown in the figure as
\theta
\cosE=
| x | |
| a |
,
and
\sinE=
| y | |
| b |
The second equation is established using the relationship
| \left( | y |
| b |
\right)2=1-\cos2E=\sin2E
which implies that . The equation is immediately able to be ruled out since it traverses the ellipse in the wrong direction. It can also be noted that the second equation can be viewed as coming from a similar triangle with its opposite side having the same length y as the distance from P to the major axis, and its hypotenuse b equal to the semi-minor axis of the ellipse.
The eccentricity e is defined as:
e=\sqrt{1-\left(
| b | |
| a |
\right)2} .
From Pythagoras's theorem applied to the triangle with r (a distance FP) as hypotenuse:
\begin{align} r2&=b2\sin2E+(ae-a\cosE)2\\ &=a2\left(1-e2\right)\left(1-\cos2E\right)+a2\left(e2-2e\cosE+\cos2E\right)\\ &=a2-2a2e\cosE+a2e2\cos2E\\ &=a2\left(1-e\cosE\right)2\\ \end{align}
Thus, the radius (distance from the focus to point P) is related to the eccentric anomaly by the formula
r=a\left(1-e\cos{E}\right) .
With this result the eccentric anomaly can be determined from the true anomaly as shown next.
The true anomaly is the angle labeled
\theta
Using the formula for above, the sine and cosine of are found in terms of :
\begin{align} \cosE&=
| x | |
| a |
=
| ae+r\cosf | |
| a |
=e+(1-e\cosE)\cosf\\ ⇒ \cosE&=
| e+\cosf | |
| 1+e\cosf |
\\ \sinE&=\sqrt{1-\cos2E }=
| \sqrt{1-e2 | |
\sinf}{1+e\cosf}~. \end{align}
Hence,
\tanE=
| \sinE | |
| \cosE |
=
| \sqrt{1-e2 | |
\sinf}{e+\cosf}~.
where the correct quadrant for is given by the signs of numerator and denominator, so that can be most easily found using an atan2 function.
Angle is therefore the adjacent angle of a right triangle with hypotenuse
1+e\cosf ,
e+\cosf ,
\sqrt{1-e2 }\sinf .
Also,
| \tan | f |
| 2 |
=\sqrt{
| 1+e | |
| 1-e |
}\tan
| E | |
| 2 |
Substituting as found above into the expression for, the radial distance from the focal point to the point, can be found in terms of the true anomaly as well:[2]
r=
| a\left(1-e2\right) | |
| 1+e\cosf |
=
| p | |
| 1+e\cosf |
p\equiva\left(1-e2\right)
The eccentric anomaly E is related to the mean anomaly M by Kepler's equation:[3]
M=E-e\sinE
This equation does not have a closed-form solution for E given M. It is usually solved by numerical methods, e.g. the Newton–Raphson method. It may be expressed in a Fourier series as
E=M+
| infty | |
| 2\sum | |
| n=1 |
| Jn(ne) | |
| n |
\sin(nM)
Jn(x)
. Elements of analytic geometry . George Albert Wentworth . George A. Wentworth . 141 . The ellipse §126 . 2nd . Ginn & Co. . 1914.