Name

CL_ex_lyddane — Lydanne orbit extrapolation analytic model

Calling Sequence

   [mean_kep_t2,osc_kep_t2] = CL_ex_lyddane(t1,mean_kep_t1,t2[,er[,mu[,zonals]]])
   
   

Description

  • Computes the mean parameters mean_kep_t2 and the osculatory parameters osc_kep_t2 at dates t2 given mean parameters mean_kep_t1 at dates t1.

    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.

  • Warnings :

    - 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

  • Last update : 28/3/2008

Parameters

t1:

CNES julian date of initial parameter (1xN)

mean_kep_t1:

keplerian mean parameters at date t1 [sma,ecc,inc,pom,gom,anm]' (6xN)

t2:

CNES julian date (1xM)

er:

(optional) equatorial radius [m] (default is %CL_eqRad)

mu:

(optional) geocentric gravitational constant [m^3/s^2] (default value is %CL_mu)

zonals:

(optional) vector of zonals coefficients J1 to Jn (troncated to J5) to be used (default is %CL_j1jn(1:5)) (1 x N)

mean_kep_t2:

keplerian mean parameters at date t2 [sma,ecc,inc,pom,gom,anm]' (6 x max(N,M))

osc_kep_t2:

keplerian osculatory parameters at date t2 [sma,ecc,inc,pom,gom,anm]' (6 x max(N,M))

Authors

CNES - DCT/SB

Bibliography

1 CNES - MSLIB FORTRAN 90, Volume E (me_lyddane)

See also

CL_ex_meanLyddane, CL_ex_lyddaneMan

Examples

t1 = 12584
mean_kep_t1 = [42166712,7.9e-3,CL_deg2rad(97.2),CL_deg2rad(89),CL_deg2rad(125),CL_deg2rad(0)]'
t2 = 12587:0.2:12588;
[mean_kep_t2,osc_kep_t2] = CL_ex_lyddane(t1,mean_kep_t1,t2,er=6378138)

See 'CelestLab > Demos > EXTRAPOLATION > lyddane' for more examples