CL_op_driftJ2 — Drift of orbital elements due to J2
[drift_pom,drift_gom,drift_M[,dpomdaei[,dgomdaei[,dMdaei]]]] = CL_op_driftJ2(sma,ecc,inc[,er[,mu[,j2]]])
Results are obtained computing the following expressions. On this formula, a is semi majors axis (sma), e is eccentricity (ecc) and i is inclination (inc). drpom, drgom and drM are, respectively, drift of the perigee argument, ascendant node longitude and mean anomaly due to j2. dpomdaei, dgomdaei and dMdaei are derivatives of, respectively, drpom, drgom and drM with respect to a, e and i. n is the mean motion.
semi major axis [m] (1xN)
eccentricity (1xN)
inclination [rad] (1xN)
(optional) equatorial radius [m] (default is %CL_eqRad)
(optional) geocentric gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) geopotentiel elliptic factor (second zonal harmonic) (default is %CL_j2)
drift of the perigee argument [rad/s] (1xN)
drift of the ascendant node longitude [rad/s] (1xN)
drift of the mean anomaly [rad/s] (1xN)
derivative of drift_pom with respect to sma, ecc and inc (3xN)
derivative of drift_gom with respect to sma, ecc and inc (3xN)
derivative of drift_M with respect to sma, ecc and inc (3xN)