Be the first to rate this file! 29 Downloads (last 30 days) File Size: 263 KB File ID: #35566
image thumbnail

Vectorized Analytic Two Body Propagator (Kepler Universal Variables)

by

 

Analytic propagation routine uses universal variables to solve a single formula for all orbit types

| Watch this File

File Information
Description

The Keplerian universal variable two body propagation method, using the Newton-Raphson iteration technique, provides a highly efficient propagation routine which doesn't have to transform coordinate systems from cartesian to Keplerian and back again. This method will work on all orbit types (circular, eccentric, parabolic, and hyperbolic) while providing a positional error which is less than 0.5 km (this is obtained from propagation comparisons with the STK implementation of the analytic two body propagation routine).

Example function call is shown below:
>> [r,v] = keplerUniversal(r0,v0,t,mu);

Where:

Inputs:
r0 is the initial position vectors in the ECI coordinate frame of reference (km) [3xN]
v0 is the initial velocity vectors in ECI coordinate frame of reference (km/s) [3xN]
t is the time vector in seconds [1xN]
mu is the Gravitational Constant of the central body which defaults to Earth if not specified (double)

Outputs:
r is the propagated position vector in the ECI coordinate frame of reference (km) [3xN]
v is the propagated velocity vector in ECI coordinate frame of reference (km/s) [3xN]

MATLAB release MATLAB 7.11 (R2010b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
11 Apr 2012 Darin Koblick

John,

The zip file seems to unzip without problems. Try using either the default windows zip utility, the MATLAB command unzip(), or the free 7zip utility. If all else fails, you can view the source code here:

http://www.mathworks.com/matlabcentral/fileexchange/35566-vectorized-analytic-two-body-propagator-kepler-universal-variables/content/keplerUniversal.m

11 Apr 2012 John Fuller

Downloaded but can't seem to unzip. Not sure if the file is corrupted, may want to check.

Contact us