CL_3b_lissajous

Description

• Computes a Lissajous orbit around a Lagrangian point.

• Notes:

  - Before using this function, it is needed to create an "environment" (env) for the chosen libration point and the chosen system (see CL_3b_environment).

  - For the definition of adimensional quantities, see CL_3b_environment.

Parameters

env:

(struct) Lagrangian point structure

Ax:

Amplitude in the "x" direction (adimensional)

phix:

Phase in the "x" direction [rad]

Az:

Amplitude in the "z" direction (adimensional)

phiz:

Phase in the "z" direction [rad]

epsilon:

Accuracy for convergence

t_orb:

Time instants at which the orbit is computed (adimensional) (1xN)

orb:

Position and velocity [x;y;z;vx;vy;vz] defining the orbit (adimensional) (6xN)

omega:

Orbit's angular rate in "xy" plane (adimensional)

nu:

Orbit's angular rate along "z" axis (adimensional)

See also

• CL_3b_environment

• CL_3b_halo

Authors

• CNES - DCT/SB

Examples

env = CL_3b_environment("S-EM", "L2"); Ax = 30.e6 / env.D; // adimensional phix = 0; Az = 100.e6 / env.D; // adimensional phiz = 0; epsilon = 1.e-10; t_orb = linspace(0,360*86400,100) * env.OMEGA; // 360 days -> adimensional orb = CL_3b_lissajous(env, Ax, phix, Az, phiz, epsilon, t_orb); // Plot orbit (normalized coordinates) scf(); param3d((orb(1,:)-env.gl), orb(2,:), orb(3,:)); xtitle("Orbit (3D), origin = Lagrange point");