Examples
t1 = 21915; mean_t1 = [42166712 ; 7.9e-3 ; CL_deg2rad([97.2; 89; 125; 0])]; t2 = [21916 21918 21920]; tman = [21915.5 21916.5]; dvman = [ [ 1 ; 1 ; 1] , [ 1 ; 0 ; 1]]; [mean_t2,osc_t2] = CL_ex_lyddaneMan(t1,mean_t1,t2,tman,dvman,er=6378138) |  |  |
Lyddane orbit propagation analytical model with maneuvers - DEPRECATED
[mean_t2,osc_t2] = CL_ex_lyddaneMan(t1,mean_t1,t2 [,tman,dvman,er,mu,j1jn])
This function is deprecated.
Replacement function: CL_ex_propagateMan
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 type of orbital elements is keplerian.
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
Warning :
- The input argument "zonals" is deprecated as of CelestLab v3.0.0. It has been replaced by "j1jn".
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
t1 = 21915; mean_t1 = [42166712 ; 7.9e-3 ; CL_deg2rad([97.2; 89; 125; 0])]; t2 = [21916 21918 21920]; tman = [21915.5 21916.5]; dvman = [ [ 1 ; 1 ; 1] , [ 1 ; 0 ; 1]]; [mean_t2,osc_t2] = CL_ex_lyddaneMan(t1,mean_t1,t2,tman,dvman,er=6378138) | | |