CL_gm_intersectCoplanOrb

Description

• Computes the argument of latitude (= true anomaly + arg of periapsis) of the intersection of 2 coplanar orbits.

  The orbit are defined by 3 parameters:

  sma: Semi major axis (positive),

  ecc: Excentricity,

  pom: Argument of periapsis = angle between some arbitrary direction in the orbit plane and the periapsis.

  The number of solutions (0, 1 or 2) is returned in numsol.

• Notes:

  - The orbits can be of any type (elliptical or hyperbolic).

  - If the number of intersections is infinite, no solution is returned (numsol == 0).

  - If there are no solutions, both psoa and psob are set to %nan.

  - If there is 1 solution, it is returned in psoa, and psob is set to %nan.

  - If there are 2 solutions, psoa contains the one with the smaller radius.

Parameters

sma1:

Semi major axis of orbit 1 [m] (1xN or 1x1)

ecc1:

Excentricity of orbit 1 (1xN or 1x1)

pom1:

Argument of periapsis of orbit 1 [rad] (1xN or 1x1)

sma2:

Semi major axis of orbit 2 [m] (1xN or 1x1)

ecc2:

Excentricity of orbit 2 (1xN or 1x1)

pom2:

Argument of periapsis of orbit 2 [rad] (1xN or 1x1)

psoa:

Argument of latitude of 1st intersection (or %nan) [rad] (1xN)

psob:

Argument of latitude of 2nd intersection (or %nan) [rad] (1xN)

numsol:

Number of intersections (1xN)

Authors

• CNES - DCT/SB

Examples

// Define 2 elliptical orbits with 2 intersections [sma1, ecc1] = CL_op_rarp2ae(10.e6, 8.e6); [sma2, ecc2] = CL_op_rarp2ae(11.e6, 7.e6); // lower periapsis, higher apoapsis) pom1 = 0; // rad pom2 = 1; // rad [psoa, psob, numsol] = CL_gm_intersectCoplanOrb(sma1, ecc1, pom1, sma2, ecc2, pom2); // Check result res1 = CL_kp_characteristics(sma1, ecc1, [psoa, psob] - pom1); res2 = CL_kp_characteristics(sma2, ecc2, [psoa, psob] - pom2); res1.r - res2.r // => 0