<< CL_man_dvBiElliptic Trajectory and maneuvers CL_man_dvHohmannG >>

celestlab >> Trajectory and maneuvers > CL_man_dvHohmann

CL_man_dvHohmann

Hohmann transfer

Calling Sequence

[deltav, dv1, dv2, anv1, anv2] = CL_man_dvHohmann(ai, af [[, mu, res="d"]])
man = CL_man_dvHohmann(ai, af [[, mu, res="s"]])

Description

Parameters

ai :

Semi-major axis of initial circular orbit. [m] (1xN or 1x1)

af :

Semi-major axis of final circular orbit. [m] (1xN or 1x1)

mu :

(optional) Gravitational constant. Default is %CL_mu. [m^3/s^2] (1x1)

res :

(string, optional) Type of output: "d" for multiple outputs or "s" for a single structure output. Default is "d".

deltav :

Som of norms of velocity increments. [m/s] (1xN)

dv1:

First velocity increment in cartesian coordinates in the "qsw" local orbital frame. [m/s]. (3xN)

dv2:

Second velocity increment in cartesian coordinates in the "qsw" local orbital frame. [m/s] (3xN)

anv1:

True anomaly at the location of the first velocity increment (in the initial orbit): as the initial orbit is circular, anv1 is arbitrarily set to 0. [rad] (1xN)

anv2:

True anomaly at the location of the second velocity increment (in the intermediate orbit). [rad] (1xN)

man:

Structure containing all the output data.

Authors

See also

Examples

// 7200km to 7000km :
ai = 7200.e3;
af = 7000.e3;
[deltav, dv1, dv2, anv1, anv2] = CL_man_dvHohmann(ai, af);

// Check results:
kep0 = [ai; 0; 1; 0; 0; 0];
kep0(6) = CL_kp_v2M(kep0(2), anv1);
kep1 = CL_man_applyDvKep(kep0, dv1, dv_frame = "qsw");
kep1(6) = CL_kp_v2M(kep1(2), anv2);;
kep2 = CL_man_applyDvKep(kep1, dv2, dv_frame = "qsw")

Report an issue
<< CL_man_dvBiElliptic Trajectory and maneuvers CL_man_dvHohmannG >>