CL_ex_lyddane — Lydanne orbit extrapolation analytic model
[mean_kep_t2,osc_kep_t2] = CL_ex_lyddane(t1,mean_kep_t1,t2[,er[,mu[,zonals]]])
Zonals coefficients up to J5 are taken into account
This function is vectorized for mean_kep_t1 or t2; but not for both at same time. When size(mean_kep_t1)=6xN and size(t2)=1xM, M or N must be 1.
- This function does not work for inclinations close (0.0022918 deg) to the critical inclinations (63.43494882 deg and 116.5650512 deg)
- This function does not work for eccentricities higher than 0.9
CNES julian date of initial parameter (1xN)
keplerian mean parameters at date t1 [sma,ecc,inc,pom,gom,anm]' (6xN)
CNES julian date (1xM)
(optional) equatorial radius [m] (default is %CL_eqRad)
(optional) geocentric gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) vector of zonals coefficients J1 to Jn (troncated to J5) to be used (default is %CL_j1jn(1:5)) (1 x N)
keplerian mean parameters at date t2 [sma,ecc,inc,pom,gom,anm]' (6 x max(N,M))
keplerian osculatory parameters at date t2 [sma,ecc,inc,pom,gom,anm]' (6 x max(N,M))