Eccentricity and argument of periapsis of a frozen orbit
[ecc,pom] = CL_op_frozenOrbit(sma,inc [,er,j1jn])
Computes the eccentricity and argument of periapsis of a frozen orbit, so that the mean value of eccentricity and argument of periapsis remain constant over time.
Freezing the orbit is possible by balancing the effects of harmonic even and odd terms of the potential (J1=0, J2, J3, ...).
The argument of periapis returned is 90 degrees or -90 degrees. The eccentricity is usually small: of the order of 1.e-3 for the Earth.
If j1jn are limited to J2 and J3, the frozen orbit eccentricity is given by:
eccg = -0.5 * sin(inc) * (Req/sma) * (J3/J2)
Note:
The relative accuracy decreases as the number of zonal terms increases because of numerical errors. The result can be very inaccurate for 40 terms or more.
Warning :
- The input argument "zonals" is deprecated as of CelestLab v3.0.0. It has been replaced by "j1jn".
Semi major axis [m] (1xN or 1x1)
Inclination [rad] (1xN or 1x1)
(optional) Equatorial radius [m]. Default is %CL_eqRad.
(optional) Vector of zonal coefficients J1 to Jn (1xNz or Nzx1). Default is %CL_j1jn(1:3).
Mean eccentricity (1xN)
Mean argument of periapsis [rad] (1xN)
CNES - DCT/SB