Code covered by the BSD License  

Highlights from
The Gravity Perturbed Hohmann Transfer

The Gravity Perturbed Hohmann Transfer

by

 

28 Feb 2013 (Updated )

MATLAB script for solving the Hohmann transfer problem perturbed by non-spherical Earth gravity.

ceqm1 (t, y)
function ydot = ceqm1 (t, y)

% first order form of Cowell's equations of orbital motion

% version for ode45

% input

%  t = current simulation time
%  y = current eci state vector

% output

%  ydot = eci acceleration vector

% Orbital Mechanics with MATLAB

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% compute non-spherical gravity perturbations

agrav = gravity(t, y);

% total acceleration vector

ydot = [ y(4)
         y(5)
         y(6)
         agrav(1)
         agrav(2)
         agrav(3)];

Contact us