Delta V for a Hohmann transfer
[delta_v,dv1,dv2,anv1,anv2] = CL_man_hohmann(ai,af [,mu])
Computes the maneuvers of a Hohmann transfer from an initial circular orbit with semi major axis ai to a final circular orbit with semi major axis af.
The output argument delta_v is the sum of the norms of the velocity increments required (|dv1| + |dv2|).
Velocity increments are expressed in spherical coordinates in the QSW frame: [lambda; phi; dv], where lambda is the in-plane angle (+%pi: towards planet and +%pi/2: ~along velocity), phi is the out-of-plane angle, positive towards the angular momentum vector (the angular momentum vector is perpendicular to the orbit plane and oriented according to right hand rule), dv is the norm of the velocity increment.
Semi-major axis of initial circular orbit. [m] (1xN)
Semi-major axis of final circular orbit. [m] (1xN)
(optional) Gravitational constant. [m^3/s^2] (default value is %CL_mu)
Total delta-v required = |dv1| + |dv2|. [m/s] (1xN)
First velocity increment, in spherical coordinates in the QSW frame [lambda;phi;|dv|] [rad,rad,m/s]. (3xN)
Second velocity increment, in spherical coordinates in the QSW frame [lambda;phi;|dv|]. [rad,rad,m/s] (3xN)
True anomaly at the location of the first velocity increment (in the initial orbit): as the initial orbit is circular, anv1 is set to 0 arbitrarily. (1xN)
True anomaly at the location of the second velocity increment (in the intermediate orbit): %pi if ai < af, 0 otherwise. [rad] (1xN)
CNES - DCT/SB