Lyddane orbit propagation analytical model with maneuvers
[mean_t2[,osc_t2]] = CL_ex_lyddaneMan(t1,mean_t1,t2 [,tman,dvman,er,mu,zonals])
Computes the mean and osculating orbital elements mean_t2 and osc_t2 at time t2 given the mean elements mean_t1 at time t1.
Maneuvers (i.e. velocity increments defined by tman,dvman) are taken into account.
tman defines the time of the maneuver.
dvman = [delta-Vx;delta-Vy;delta-Vz] defines the delta-V in the QSW local frame at the time of maneuver.
If one maneuver occurs at one of the final times, the output orbital elements mean_t2 and osc_t2 include the effect of the maneuver at that time.
The final times as well as the maneuver time should be sorted (increasing order).
Zonals coefficients up to J5 are taken into account.
The orbital elements are the following:
Warnings :
- This function does not work for inclinations close to the critical inclinations (63.43494882 deg and 116.5650512 deg)
- This function does not work for eccentricities higher than 0.9
Initial time [days] (1x1)
Mean Keplerian elements at time t1 [sma;ecc;inc;pom;raan;anm] (6x1)
Final time [days] (1xN)
(optional) Maneuver times [days] (1xP)
(optional) Delta-V in QSW local frame [dvx;dvy;dvz] [m/s] (3xP)
(optional) Equatorial radius [m] (default is %CL_eqRad)
(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) Vector of zonal coefficients J1 to Jn (troncated to J5) to be used (default is %CL_j1jn(1:5)) (1xM)
Mean Keplerian elements at t2 [sma;ecc;inc;pom;raan;anm] (6 x N)
Osculating Keplerian elements at t2 [sma;ecc;inc;pom;raan;anm] (6 x N)
CNES - DCT/SB