Examples
mean_kep0 = [7.e6; 1.e-3; 1; %pi/2; 0.1; 0.2]; t0 = 0; t = 1:4; osc_kep = CL_ex_propagate("eckhech","kep",t0,mean_kep0,t,"o") |  |  |
Orbit propagation (all analytical models)
[result1, result2] = CL_ex_propagate(mod, type_oe, t1, mean_oe_t1, t2, res [, er, mu, j1jn])
Propagates orbital elements using one analytical model.
The available propagation models are:
"central": Central force (osculating elements = mean elements)
"j2sec": Secular effects of J2 (osculating elements = mean elements by convention)
"lydsec": Lyddane (mean elements include secular effects only)
"lydlp": Lyddane (mean elements include secular and long period effects)
"eckhech": Eckstein-Hechler (mean elements include secular and long period effects)
Notes:
- There can be 1 or N initial times, and 1 or N final times.
- Conversions take place if the type of orbital elements is not the "natural" type for the model. There is no conversion for the "central" and "j2sec" models.
See Propagation models for more details.
(string) Model name: "central", "j2sec", "lydsec", "lydlp", "eckhech". (1x1)
(string) Type of orbital elements used for input/output: "kep", "cir", "cireq" or "equin" (1x1)
Initial time [days] (1x1 or 1xN)
Mean orbital elements at time t1 (6x1 or 6xN)
Final time [days] (1xN or 1x1)
(string) Type of output (mean or osculating): "m", "o", "mo", "om" (1x1)
(optional) Equatorial radius [m]. Default is %CL_eqRad
(optional) Gravitational constant [m^3/s^2]. Default is %CL_mu
(optional) Vector of zonal harmonics. Default is %CL_j1jn (Nz x 1)
Mean or osculating orbital elements at t2 (6xN)
CNES - DCT/SB
mean_kep0 = [7.e6; 1.e-3; 1; %pi/2; 0.1; 0.2]; t0 = 0; t = 1:4; osc_kep = CL_ex_propagate("eckhech","kep",t0,mean_kep0,t,"o") | | |