Eccentric anomaly

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

In orbital mechanics, 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.

Graphical Representation

The eccentric anomaly of point P is the angle E. The center of the ellipse is point C, and the focus is point F.

Consider the ellipse with equation given by:

\frac{x^2}{a^2} + \frac{y^2}{b^2} = 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 to the right. The eccentric anomaly, E, is observed by drawing a right triangle with one vertex at the center of the ellipse, 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 crosses through the point P. The eccentric anomaly is measured in the same direction as the true anomaly, shown in the figure as f. The eccentric anomaly E in terms of these coordinates is given by:[1]

\cos E = \frac{x}{a} ,

and

\sin E = \frac{y}{b}

The second equation is established using the relationship {\left( \frac{y}{b}\right)}^2 = 1-{\cos}^2 E={\sin}^2 E, which implies that \sin E = \pm \frac{y}{b}. The equation \sin E = - \frac{y}{b} 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 being the similar triangle with adjacent side through P and the minor auxiliary circle, hypotenuse b, and whose opposite side is y.

Formulas

Radius and eccentric anomaly

The eccentricity e is defined as:

e=\sqrt{1 - \left(\frac{b}{a}\right)^2 } \ .

From Pythagoras' theorem applied to the triangle with r (a distance FP) as hypotenuse:

\begin{align} r^2 &= b^2 \sin^2E + (ae-a\cos E)^2 \\
&=a^2(1-e^2)(1-\cos^2 E)+a^2 (e^2 -2e \cos E +\cos^2 E)\\
&=a^2  -2a^2e \cos E +a^2e^2 \cos^2 E \\
&=a^2 (1-e \cos E )^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 \cdot \cos{E} \right ) \ .

With this result the eccentric anomaly can be determined from the true anomaly as shown next.

From the true anomaly

The true anomaly is the angle labeled f in the figure, located at the focus of the ellipse; it is often referred to as θ as in the calculations below. The true anomaly and the eccentric anomaly are related as follows.[2]

Using the formula for r above, the sine and cosine of E are found in terms of θ:

  \cos E = \frac{x}{a} =  \frac{ae +r \cos \theta}{a} = e+ (1-e \cos E) \cos \theta \ \to \cos E = \frac{ e + \cos \theta }{1 + e \cos \theta }
  \sin  E = \sqrt{1 - \cos^2 E}  =  \frac{ \sqrt{1 - e^2} \, \sin \theta }{1 +  e \cos \theta } \ .

Hence,

 \tan E =\frac{\sin E}{\cos E} = \frac{ \sqrt{1-e^2} \sin \theta }{e + \cos \theta} \ .

Angle E is therefore the adjacent angle of a right triangle with hypotenuse 1 + e cosθ, adjacent side e + cosθ, and opposite side √(1-e2) sinθ.

Also,

\tan \frac{\theta}{2} = \sqrt{\frac{1+e}{1-e}} \cdot \tan \frac{E}{2}

Substituting cosE as found above into the expression for r, the radial distance from the focal point to the point P, can be found in terms of the true anomaly as well:[2]

r = \frac{a \left( 1-e^2\right)}{1+e\cos \theta} \ .

From the mean anomaly

The eccentric anomaly E is related to the mean anomaly M by Kepler's equation:[3]

M =  E - e \cdot \sin E

This equation does not have a closed-form solution for E given M. It is usually solved by numerical methods, e.g. Newton-Raphson method.

In-line references and notes

  1. Lua error in package.lua at line 80: module 'strict' not found.
  2. 2.0 2.1 Lua error in package.lua at line 80: module 'strict' not found.
  3. Lua error in package.lua at line 80: module 'strict' not found.

Background references

  • Murray, Carl D.; & Dermott, Stanley F. (1999); Solar System Dynamics, Cambridge University Press, Cambridge, GB
  • Plummer, Henry C. K. (1960); An Introductory Treatise on Dynamical Astronomy, Dover Publications, New York, NY (Reprint of the 1918 Cambridge University Press edition)

See also

de:Exzentrische Anomalie