File Exchange

image thumbnail

Closest Approach Between the Earth and Heliocentric Objects

version (351 KB) by David Eagle
MATLAB script that predicts closest approach between the Earth and heliocentric objects.


Updated 02 Dec 2013

View License

The cae2ho.m script uses a Runge-Kutta-Fehlberg 7(8) numerical method to numerically integrate the first-order form of the orbital equations of motion. This is a variable step size method of order 7 with an 8th order error estimate which is used to dynamically change the integration step size during the simulation. This software also uses a one-dimensional minimization algorithm due to Richard Brent to solve the close approach problem. Additional information about this numerical method can be found in the book, Algorithms for Minimization Without Derivatives, R. Brent, Prentice-Hall, 1972. As the title indicates, this algorithm does not require derivatives of the objective function. This feature is important because the analytic first derivative of many objective functions may be difficult to derive. The objective function for this program is the scalar geocentric distance of the celestial body or spacecraft.

The lunar and planetary ephemeris used in this computer program is based on JPL DE421. A Windows compatible DE421 binary ephemeris file is available at

Cite As

David Eagle (2021). Closest Approach Between the Earth and Heliocentric Objects (, MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R2012a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!