Delta V for a bi-elliptic transfer
[delta_v,dv1,dv2,dv3,anv1,anv2,anv3]=CL_man_biElliptic(ai,af,rt [,mu])
Computes the maneuvers of a bi-elliptical transfer from a circular orbit with semi-major axis ai to a circular orbit with semi-major axis af.
The apogee radius of the elliptical transfer orbit is rt.
delta-v is the sum of the norms of the velocity increments.
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)
Radius at the position of the second maneuver [m] (1xN)
(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)
Total |delta-V| (=|dv1|+|dv2|+|dv3]) [m/s] (1xN)
First delta-V, in spherical coordinates in the QSW frame [lambda;phi;dv] [rad,rad,m/s] (3xN)
Second delta-V, in spherical coordinates in the QSW frame [lambda;phi;dv] [rad,rad,m/s] (3xN)
Third delta-V, in spherical coordinates in the QSW frame [lambda;phi;dv] [rad,rad,m/s] (3xN)
True anomaly at the position of the 1st maneuver : initial orbit is circular so this is an arbitrary value of 0 (1xN)
True anomaly at the position of the 2nd maneuver (either 0 or %pi) [rad] (1xN)
True anomaly at the position of the 3rd maneuver (either 0 or %pi) [rad] (1xN)
CNES - DCT/SB
1) Orbital Mechanics for engineering students, H D Curtis, Chapter 6 (equation 6.4a)