<< CL_op_frozenOrbitDerCir Orbit properties CL_op_incdotSsoCirc >>

celestlab >> Orbit properties > CL_op_gauss

CL_op_gauss

Gauss equations

Calling Sequence

[M, N] = CL_op_gauss(type_oe, oe, frame [[, mu]])

Description

Parameters

type_oe:

(string) Type of input orbital elements ("kep", "cir", "cireq" or "equin""). (1x1)

oe:

Input orbital elements. (6xN)

frame:

(string) Coordinates frame for acceleration : "qsw", "tnw" or "inertial". (1x1)

mu:

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

M:

Effects of acceleration on orbital elements. (6x3xN)

N:

Mean motion. [rad/s] (6xN)

Authors

Bibliography

See also

Examples

// Keplerian orbital elements
kep = [7000.e3; 0.01; 1.8; 0.1; 0.2; 0.3];
[M, N] = CL_op_gauss("kep", kep, "qsw");

// d(kep)/dt due to some acceleration:
acc = [1e-3; 1e-4; -2e-3]; // in qsw frame
dkepdt = M * acc + N;

// Same orbit, circular orbital elements
[cir, dcirdkep] = CL_oe_kep2cir(kep);
[M2, N2] = CL_op_gauss("cir", cir, "qsw");

// Check consistency using jacobian:
M2 - dcirdkep * M   // --> zero

Report an issue
<< CL_op_frozenOrbitDerCir Orbit properties CL_op_incdotSsoCirc >>