University:
California State University, Northridge
Course:
ASTR 300 | Introductory Astrophysics
Academic year:

2024
Views:

486

Pages:

2
Author:

Shyla Hinton

Numerical Solution of Kepler's equation Since the mean anomaly is proportional to time it is an easily measured quantity for an orbiting body. However, in order to compute the location of this body we need to find either the eccentric or the true anomaly. For this one needs to solve Kepler's equation for . Since this is a transcendental equation, the solution can normally not be found in closed form. However, for small we can devise a simple iteration scheme. To start we rewrite Kepler's equation as (28) In this form, the equation has the form of a fixed point equation. In Calculus one encounters the following Fixed Point Theorem For a proof of this theorem we refer to a Calculus textbook. Kepler's equation perfectly fits mold of this theorem. First of all, the function in satisfies Moreover, we have that Therefore we can use this iteration to find an approximate solution of the equation. Example: A satellite moves on an elliptic orbit with ellipticity . The mean anomaly is , Find an approximate value of the eccentric anomaly. In this case . We set . Using the iteration we get We see that after only two iterations, we get and the approximate value of the eccentric anomaly is or . The method converges not always that fast, depending on the values of and . For larger values of , we can also use Newton's method. For this we rewrite the equation as Newton's method yields the iteration (starting with ): which converges faster that the simple fixed point iteration above.

Numerical Solution of Kepler’s Equation

Related Documents

Recommended Documents

Description

Related Documents

Free up your schedule!

Take 5 seconds to unlock