mission02.py

import numpy as np
from toolbox import *

"""
Task 3 - Sun-Pointing Reference Frame Orientation

Determine an analytic expressions for the sun pointing reference frame Rs by defining the DCM [RsN]

As the axis are fixed, there is no angular velocity.
"""
matrix_RsN = np.round(EAtoDCM313(SUN_FRAME), 6)

print(f"Matrix RsN at 0s:\n{matrix_RsN}")


"""
Task 4 - Nadir-Pointing Reference Frame Orientation

Determine an analytic expression for the nadir pointing reference frame Rn by defining the DCM [RnN]
"""
## Matrix RnH
matrix_RnH = np.round(EAtoDCM313(NADIR_FRAME), 6)

## Matrix HN

## Find HN at 330s

matrix_HN = orbit_integrator('LMO', 330, 1)

## Matrix RnN at 330s
matrix_RnN = np.round(np.matmul(matrix_RnH, matrix_HN), 6)

print(f"Matrix RnN at 330s:\n{matrix_RnN}")


## Angular velocity in inertial frame

matrix_NH = matrix_HN.T

w_inertial_nadir_330 = np.matmul(matrix_NH, w_NADIR.T)

print(f"Angular velocity in N frame at 330s:\n{w_inertial_nadir_330}")


"""
Task 5 - GMO-Pointing Reference Frame Orientation

Determine an analytic expression for the communication mode reference frame Rc by defining DCM [RcN]
"""