Summary

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

Part of a series on
Astrodynamics
Orbital mechanics
Orbital parameters Apsis Argument of periapsis Azimuth Eccentricity Inclination Mean anomaly Orbital nodes Semi-major axis True anomaly
Types of two-body orbits, by eccentricity Circular orbit Elliptic orbit Transfer orbit (Hohmann transfer orbit Bi-elliptic transfer orbit) Parabolic orbit Hyperbolic orbit Radial orbit Decaying orbit
Equations Dynamical friction Escape velocity Kepler's equation Kepler's laws of planetary motion Orbital period Orbital velocity Surface gravity Specific orbital energy Vis-viva equation
Celestial mechanics
Gravitational influences Barycenter Hill sphere Perturbations Sphere of influence
N-body orbits Lagrangian points (Halo orbits) Lissajous orbits Lyapunov orbits
Engineering and efficiency
Preflight engineering Mass ratio Payload fraction Propellant mass fraction Tsiolkovsky rocket equation
Efficiency measures Gravity assist Oberth effect
v t e

In celestial mechanics, true anomaly is an angular parameter that defines the position of a body moving along a Keplerian orbit. It is the angle between the direction of periapsis and the current position of the body, as seen from the main focus of the ellipse (the point around which the object orbits).

The true anomaly is usually denoted by the Greek letters ν or θ, or the Latin letter f, and is usually restricted to the range 0–360° (0–2π^c).

As shown in the image, the true anomaly f is one of three angular parameters (anomalies) that defines a position along an orbit, the other two being the eccentric anomaly and the mean anomaly.

Formulas

From state vectors

For elliptic orbits, the true anomaly ν can be calculated from orbital state vectors as:

${\displaystyle \nu =\arccos {{\mathbf {e} \cdot \mathbf {r} } \over {\mathbf {\left|e\right|} \mathbf {\left|r\right|} }}}$

(if r ⋅ v < 0 then replace ν by 2π − ν)

where:

• v is the orbital velocity vector of the orbiting body,
• e is the eccentricity vector,
• r is the orbital position vector (segment FP in the figure) of the orbiting body.

Circular orbit

For circular orbits the true anomaly is undefined, because circular orbits do not have a uniquely determined periapsis. Instead the argument of latitude u is used:

${\displaystyle u=\arccos {{\mathbf {n} \cdot \mathbf {r} } \over {\mathbf {\left|n\right|} \mathbf {\left|r\right|} }}}$

(if r_z < 0 then replace u by 2π − u)

where:

• n is a vector pointing towards the ascending node (i.e. the z-component of n is zero).
• r_z is the z-component of the orbital position vector r

Circular orbit with zero inclination

For circular orbits with zero inclination the argument of latitude is also undefined, because there is no uniquely determined line of nodes. One uses the true longitude instead:

${\displaystyle l=\arccos {r_{x} \over {\mathbf {\left|r\right|} }}}$

(if v_x > 0 then replace l by 2π − l)

where:

• r_x is the x-component of the orbital position vector r
• v_x is the x-component of the orbital velocity vector v.

From the eccentric anomaly

The relation between the true anomaly ν and the eccentric anomaly E is:

${\displaystyle \cos {\nu }={{\cos {E}-e} \over {1-e\cos {E}}}}$

or using the sine^[1] and tangent:

${\displaystyle {\begin{aligned}\sin {\nu }&={{{\sqrt {1-e^{2}\,}}\sin {E}} \over {1-e\cos {E}}}\\[4pt]\tan {\nu }={{\sin {\nu }} \over {\cos {\nu }}}&={{{\sqrt {1-e^{2}\,}}\sin {E}} \over {\cos {E}-e}}\end{aligned}}}$

or equivalently:

${\displaystyle \tan {\nu \over 2}={\sqrt {{1+e\,} \over {1-e\,}}}\tan {E \over 2}}$

so

${\displaystyle \nu =2\,\operatorname {arctan} \left(\,{\sqrt {{1+e\,} \over {1-e\,}}}\tan {E \over 2}\,\right)}$

An equivalent form avoids the singularity as e → 1, however it does not produce the correct value for ${\displaystyle \nu \,}$:

${\displaystyle \nu =2\,\operatorname {arg} \left({\sqrt {1-e\,}}\,\cos {\frac {E}{2}},\;{\sqrt {1+e\,}}\sin {\frac {E}{2}}\right)}$

or, with the same problem as e → 1 ,

${\displaystyle \nu =\operatorname {arg} \left(\cos {E}-e,\;{\sqrt {1-e^{2}\,}}\sin {E}\right)}$.

In both of the above, the function arg(x, y) is the polar argument of the vector (x y), available in many programming languages as the library function named atan2(y,x) (note the reversed order of x and y).

From the mean anomaly

The true anomaly can be calculated directly from the mean anomaly via a Fourier expansion:^[2]

${\displaystyle \nu =M+\left(2e-{\frac {1}{4}}e^{3}\right)\sin {M}+{\frac {5}{4}}e^{2}\sin {2M}+{\frac {13}{12}}e^{3}\sin {3M}+\operatorname {O} \left(e^{4}\right)}$

where the ${\displaystyle \operatorname {O} \left(e^{4}\right)}$ means that the omitted terms are all of order e⁴ or higher. Note that for reasons of accuracy this approximation is usually limited to orbits where the eccentricity (e) is small.

The expression ${\displaystyle \nu -M}$ is known as the equation of the center.

Radius from true anomaly

The radius (distance between the focus of attraction and the orbiting body) is related to the true anomaly by the formula

${\displaystyle r=a\,{1-e^{2} \over 1+e\cos \nu }\,\!}$

where a is the orbit's semi-major axis.

See also

• Kepler's laws of planetary motion
• Eccentric anomaly
• Mean anomaly
• Ellipse
• Hyperbola

References

Further reading

• Murray, C. D. & Dermott, S. F., 1999, Solar System Dynamics, Cambridge University Press, Cambridge. ISBN 0-521-57597-4
• Plummer, H. C., 1960, An Introductory Treatise on Dynamical Astronomy, Dover Publications, New York. OCLC 1311887 (Reprint of the 1918 Cambridge University Press edition.)

External links

• Federal Aviation Administration - Describing Orbits