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)
[mean/osc_oe_t2, mean/osc_oe_t2] = CL_ex_propagate(mod, type_oe, t1, mean_oe_t1, t2, res [, er, mu, j1jn])
Propagates orbital elements from one initial time to one final time.
This function is a wrap-up for all available analytical models.
The available 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.
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" or "cir" (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") | | |