コード例 #1
0
def test_flight_ang():
    r""" Test for Radius of Apoapsis Function
    """

    e = .5
    theta = 1.2
    actual_out = RV2COE.flight_ang(e, theta)
    expected_out = 0.37578860669992659

    np.testing.assert_allclose(actual_out, expected_out)
コード例 #2
0
        theta = RV2COE.true_anom(r_input, h_input, e_vector)

        i = RV2COE.inclination(RV2COE.unit_vector(h_input))

        a = RV2COE.semi_major_axis(p, np.linalg.norm(e_vector))

        raan = RV2COE.R_A_A_N(n_vector)

        w = RV2COE.arg_of_periapsis(n_vector, e_vector, h_input)

        r_p = RV2COE.rad_peri(p, np.linalg.norm(e_vector))

        r_a = RV2COE.rad_apo(p, np.linalg.norm(e_vector))

        gamma = RV2COE.flight_ang(np.linalg.norm(e_vector), theta)

        period = RV2COE.Orbit_Period(mu, a)

        if (theta < 0):
            theta = 2 * np.pi + theta
        else:
            theta = theta
        if (raan < 0):
            raan = 2 * np.pi + raan
        else:
            raan = raan
        if (w < 0):
            w = 2 * np.pi + w
        else:
            w = w