Eccentricity and argument of periapsis of a frozen orbit
[ecc,pom] = CL_op_frozenOrbit(sma,inc,[,er,zonals])
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 zonals are limited to J2 and J3, the frozen eccentricity is given by:
eccg = -0.5*sin(inc)*(Req/sma)*J3/J2
NB: The relative accuracy decreases as the number of zonal terms increases. The result can be very inaccurate for 40 terms or more.
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