Name

CL_oe_car2kep — Cartesien to keplerian orbital elements

Calling Sequence

   [kep[,jacob]] = CL_oe_car2kep(pos_car,vel_car[,mu])
   
   

Description

  • Given the cartesien parameters (position and velocity), the function computes the keplerien parameters [a,e,i,w,omega,M] for a non circular and non equatorial orbit (elliptic, hyperbolic or parabolic orbits).

    The transformation jacobian is optionally computed.

  • Warning : transformation is not exact for eccentricites and/or inclinations close to 0
  • Last update : 29/11/2007

Parameters

pos_car:

position in cartesien coordinates [x;y;z] [m] (3xN)

vel_car:

velocity in cartesien coordinates [vx;vy;vz] [m/s] (3xN)

mu :

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

kep:

keplerien parameters [sma,ecc,inc,pom,gom,anm] [m],[rad] (6xN)

jacob:

(optional) transformation jacobian d(sma,ecc,inc,pom,gom,anm)/d(x,y,z,vx,vy,vz) (6x6xN)

Bibliography

1 CNES - MSLIB FORTRAN 90, Volume V (mv_car_kep)

See also

CL_oe_kep2car, CL_oe_car2cir, CL_oe_car2cirEqua, CL_co_car2sph, CL_co_car2ell

Authors

CNES - DCT/SB

Examples

// Example 1
kep=[7000000,0.001,0.9,%pi/3,%pi/4,0]'
[pos_car,vel_car,jacob1]=CL_oe_kep2car(kep)
[kep2,jacob2]=CL_oe_car2kep(pos_car,vel_car)

// Example 2
pos = [-26655470,29881667,-113657]';
vel = [-1125,-1122,195]';
[kep,jacob1] = CL_oe_car2kep(pos,vel);
[pos2,vel2,jacob2] = CL_oe_kep2car(kep);