Cartesian to keplerian orbital elements
[kep[,jacob]] = CL_oe_car2kep(pos_car,vel_car[,mu])
Converts position and velocity to classical keplerian orbital elements for elliptic (non circular and non equatorial), hyperbolic or parabolic orbit.
The transformation jacobian is optionally computed.
See Orbital elements for more details
Notes:
- for eccentricites close to 0, ω and M are not defined independently. Only ω + M is.
- for inclinations close to 0, ω and Ω are not defined independently. Only ω + Ω is.
- for inclinations and eccentricities close to 0, ω, Ω and M are not defined independently. Only ω + Ω + M is.
Jacobian transformation is not exact for any of the above cases
position [X;Y;Z] [m] (3xN)
velocity [Vx;Vy;Vz] [m/s] (3xN)
(optional) gravitational constant. [m^3/s^2] (default value is %CL_mu)
classical keplerian orbital elements [sma;ecc;inc;pom;raan;anm] [m,rad] (6xN)
(optional) transformation jacobian d(sma,ecc,inc,pom,raan,anm)/d(X,Y,Z,Vx,Vy,Vz) (6x6xN)
1) CNES - MSLIB FORTRAN 90, Volume V (mv_car_kep)
CNES - DCT/SB