Formula & Calculator
Kepler's Equation
Transcendental relation between mean anomaly, eccentric anomaly, and eccentricity, solved iteratively to locate a body in its orbit at a given time.
Interpretation
Kepler's equation: M = E − e·sin(E), where M is mean anomaly, E is eccentric anomaly, e is eccentricity. It relates time to position in an elliptic orbit. Example: e=0.1, E=1 rad → M = 1 − 0.1×0.8415 = 0.9158 rad.
Variables
| Symbol | Quantity | Unit |
|---|---|---|
| M | Mean anomaly | rad |
| E | Eccentric anomaly | rad |
| e | Eccentricity |
What it means
Kepler’s equation is a fundamental relation in celestial mechanics that links the mean anomaly (which increases uniformly with time) to the eccentric anomaly (which describes the position on the ellipse). This equation is transcendental and must be solved iteratively. It is used to compute the position of a body in an orbit at a given time, which is essential for ephemeris prediction and orbit propagation. Understanding Kepler’s equation is crucial for astrodynamics and for mission planning.
Worked example
Kepler's Equation – Two Examples
Real‑World| Parameter | Value |
|---|---|
| E | 1.0 rad |
| e | 0.1 |
| Parameter | Value |
|---|---|
| E | 1.5 |
| e | 0.3 |
Common mistakes
- Kepler’s equation: M = E − e·sin(E).
- M: Mean anomaly (radians).
- E: Eccentric anomaly (radians).
- e: Orbital eccentricity.
- Transcendental equation – solve iteratively (e.g., Newton’s method).
Applications
Kepler's equation, M = E − e·sinE, relates the mean anomaly (time) to the eccentric anomaly in an elliptical orbit. It is essential for predicting satellite positions over time. Engineers use this to compute true anomaly from time, to design orbital events (e.g., eclipse, ground track), and to plan manoeuvres. Solving Kepler's equation is a core task in astrodynamics, and it underpins orbit propagation software. By using this equation, aerospace engineers can accurately determine the position of a spacecraft at any given time.
- Orbit propagation and prediction of satellite positions
- Mission timing and event planning (eclipses, ground station visibility)
- Rendezvous and flyby calculations
- Orbital mechanics education and software implementation
- Numerical integration and error analysis
Frequently Asked Questions
It is a transcendental relation between mean anomaly, eccentric anomaly, and eccentricity, solved iteratively to locate a body in its orbit at a given time.
M = mean anomaly (radians)
E = eccentric anomaly (radians)
e = eccentricity
It allows the conversion from time (mean anomaly) to position (true anomaly) in an orbit. It is essential for orbit determination and ephemeris generation.
- Attempting to solve Kepler's equation for E directly in closed form; it requires iterative (e.g. Newton‑Raphson) numerical solution.
- Using degrees instead of radians.
- Forgetting that M = n·(t − tp), where n is the mean motion.
For e = 0.5 and M = 1.0 rad, solve E = M + e·sin(E) iteratively. Initial guess E = M = 1.0. Iteration: E_new = 1.0 + 0.5·sin(1.0) = 1.0 + 0.4207 = 1.4207. Continue: E_new = 1.0 + 0.5·sin(1.4207) = 1.0 + 0.490 = 1.490. Next: 1.0 + 0.5·sin(1.49) = 1.0 + 0.497 = 1.497. Converges to E ≈ 1.5 rad.
Using tan(θ/2) = √((1+e)/(1−e))·tan(E/2).
n = √(μ/a³), the average angular rate of the orbiting body.
For small e, the equation is nearly linear. For large e, the convergence may be slower, and the solution is more sensitive.
Kepler's equation becomes Barker's equation, which has a closed‑form solution.