The orbit in (any) time: connecting the mean anomaly (linear in time) to the true anomaly (geometry) through the eccentric anomaly. See also ENAE 441: Kepler’s Equation.

Mean Anomaly as a Function of Time

  • From any epoch 𝑡0 (not only perigee):
    • 𝑀−𝑀0=𝑛(𝑡−𝑡0)
    • ̇𝑀=𝑛=√𝜇𝑎3

Kepler’s Equation

  • The area swept out since perigee, in terms of the eccentric anomaly 𝐸:
    • 𝐴(𝑡)=12𝑎𝑏(𝐸−𝑒sin𝐸)
  • Substitute into the areal definition of the mean anomaly:
    • 𝑀=2𝐴(𝑡)𝑎𝑏=𝐸−𝑒sin𝐸

Kepler's Equation

  • 𝑀=𝐸−𝑒sin𝐸
  • Easy direction: 𝐸→𝑀. Given 𝑀 (from time), solving for 𝐸 requires iteration

Eccentric and True Anomaly

  • Position in terms of 𝐸 (Cartesian, origin at the focus):
    • 𝑟2=𝑎2(1−𝑒2)sin2𝐸+𝑎2(𝑒−cos𝐸)2
  • Expand:
    • 𝑟2=𝑎2(1+𝑒2−𝑒2sin2𝐸⏟𝑒2cos2𝐸−2𝑒cos𝐸)=𝑎2(1−𝑒cos𝐸)2
  • Take the square root:

Radius-Eccentric Anomaly Equation

  • 𝑟=𝑎(1−𝑒cos𝐸)
  • Equate with the orbit equation (𝑝=𝑎(1−𝑒2)):
    • 𝑎(1−𝑒cos𝐸)=𝑎(1−𝑒2)1+𝑒cos𝜃
  • Relations between 𝐸 and 𝜃:
    • cos𝜃=cos𝐸−𝑒1−𝑒cos𝐸Curtis (3.10b)cos𝐸=𝑒+cos𝜃1+𝑒cos𝜃Curtis (3.10a)sin𝜃=√1−𝑒2sin𝐸1−𝑒cos𝐸sin𝐸=√1−𝑒2sin𝜃1+𝑒cos𝜃Curtis (3.11)
  • Quadrant-safe forms with the two-argument arctangent:
    • 𝐸=arctan(√1−𝑒2sin𝜃,𝑒+cos𝜃)𝜃=arctan(√1−𝑒2sin𝐸,cos𝐸−𝑒)

Half-Angle Relations

  • Half-angle identity:
    • tan2𝑥2=1−cos𝑥1+cos𝑥
  • Applied to 𝜃 with (3.10b):
    • tan2𝜃2=1+𝑒1−𝑒tan2𝐸2
  • 𝜃2 and 𝐸2 are always in the same quadrant, so

Half-Angle Tangent Relations

  • tan𝜃2=√1+𝑒1−𝑒tan𝐸2
  • tan𝐸2=√1−𝑒1+𝑒tan𝜃2

From True to Mean Anomaly

AnomalyNatureStep to the next
True𝜃Geometrichalf-angle tangent relations: tan𝐸2=√1−𝑒1+𝑒tan𝜃2
Eccentric𝐸IntermediateKepler’s equation: 𝑀=𝐸−𝑒sin𝐸
Mean𝑀Dynamiċ𝑀=𝑛=√𝜇𝑎3
  • Mean anomaly vs true anomaly for ellipses of various eccentricities, Curtis Fig. 3.2 (from Dr. Healy's notes)
    • Curtis Fig. 3.2: 𝑀 vs 𝜃 for 𝑒=0,0.2,0.5,0.8,0.9
    • 𝑒=0: a straight line (𝑀=𝜃)
    • Larger 𝑒: 𝑀 changes slowly near perigee (𝜃≈0, fast motion) and very quickly near apogee (𝜃≈𝜋, slow motion)

Hyperbolic Orbits

  • 𝐹 (or 𝐻) = hyperbolic anomaly, −∞≤𝐹≤+∞
  • Radius:
    • 𝑟=|𝑎|(𝑒cosh𝐹−1)
  • Relation to the true anomaly:
    • tan𝜃2=√𝑒+1𝑒−1tanh𝐹2
  • Mean motion:
    • ̇𝑀=𝑛=√𝜇|𝑎|3
  • See Hyperbolic Orbits

Parabolic Motion

  • Parabolic anomaly:
    • 𝐷=tan𝜃2

Barker's Equation

  • 𝑀=𝐷+13𝐷3
  • ̇𝑀=1𝑟32𝑝√𝜇2
  • 𝑟𝑝: radius of perigee