Code covered by the BSD License  

Highlights from
A MATLAB Script for Predicting the Evolution of Lunar Orbits

A MATLAB Script for Predicting the Evolution of Lunar Orbits

by

 

Script for propagating lunar orbits subject to non-spherical lunar gravity and third-body gravity.

moon_pa(jdate)
function tmatrix = moon_pa(jdate)

% transformation matrix from lunar mean equator and IAU node of j2000
% to the lunar principal axes system using JPL binary ephemeris

% input

%  jdate = TDB julian date

% output

%  tmatrix = transformation matrix

% Orbital Mechanics with MATLAB

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

global tmatrix2
      
% compute lunar libration angles (radians)

sv = jplephem (jdate, 15, 0);

phi = sv(1);

theta = sv(2);

psi = sv(3);

% compute lunar libration matrix

tmatrix1 = matran (phi, 3, theta, 1, psi, 3, 0.0d0, 0);

% create moon_j2000 to lunar principal axes transformation matrix

tmatrix = tmatrix1 * tmatrix2;





Contact us