You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
68 lines
2.1 KiB
68 lines
2.1 KiB
|
|
#include <catch2/catch.hpp> |
|
|
|
#include <glm/glm.hpp> |
|
|
|
#include "../src/orbits.cpp" |
|
|
|
// NOTE: using WithinAbs instead of Approx because: |
|
// https://github.com/catchorg/Catch2/issues/1962 |
|
// https://github.com/catchorg/Catch2/issues/1507 |
|
using Catch::Matchers::WithinAbs; |
|
|
|
|
|
TEST_CASE("orbit construction", "[orbits]") |
|
{ |
|
|
|
} |
|
|
|
// NOTE: example 4.6 in "Space Flight Dynamics" by Craig A. Kluever |
|
TEST_CASE("orbit propagation", "[orbits]") |
|
{ |
|
double a = 26564.5; // NOTE: semi-major axis in km |
|
double e = 0.7411; // NOTE: eccentricity |
|
double mu = 398601.68; // NOTE: gravitational parameter |
|
double r = 6378; // NOTE: body radius in km |
|
double initial_anom = 260 * M_PI / 180; // NOTE: radians |
|
system_2body sys = systemInit(gravBodyInit(mu, r), orbitInit(a, e)); |
|
sys.time_step = 60 * 50; // NOTE: seconds |
|
|
|
double E1 = getEccAnomFromTrueAnom(sys.ep.e, initial_anom); |
|
REQUIRE_THAT(E1, WithinAbs(-0.8615, 1e-4)); |
|
|
|
double M1 = getMeanAnomFromEccAnom(E1, sys.ep.e); |
|
REQUIRE_THAT(M1, WithinAbs(-0.2992, 1e-4)); |
|
|
|
double n = getMeanMotion(sys.body.mu, sys.ep.a); |
|
REQUIRE_THAT(n, WithinAbs(0.00014582, 1e-8)); |
|
|
|
double M2 = getPropagatedMeanAnom(M1, n, sys.time_step); |
|
REQUIRE_THAT(M2, WithinAbs(0.1383, 1e-5)); |
|
|
|
// TODO: could also test for the other trial values listed in table 4.1 |
|
double E2_1 = getInitialTrialValue(M2, sys.ep.e); |
|
REQUIRE_THAT(E2_1, WithinAbs(0.315452, 1e-5)); |
|
|
|
double ecc_anom = getPropagatedEccAnomaly(sys, initial_anom); |
|
REQUIRE_THAT(ecc_anom, WithinAbs(0.481518, 1e-5)); |
|
|
|
double true_anom = getPropagatedTrueAnomaly(sys, initial_anom); |
|
REQUIRE_THAT(true_anom, WithinAbs(1.1339, 1e-4)); |
|
|
|
double r2 = orbitGetRadialDistance(sys.ep.e, sys.ep.p, true_anom); |
|
REQUIRE_THAT(r2, WithinAbs(9116.1, 0.1)); |
|
|
|
glm::vec2 pos = polarToRect(true_anom, r2); |
|
REQUIRE_THAT(pos.x, WithinAbs(3856.9, 0.1)); |
|
REQUIRE_THAT(pos.y, WithinAbs(8259.9, 0.1)); |
|
} |
|
|
|
// NOTE: example 2.5 - c |
|
TEST_CASE("orbital period", "[orbits]") |
|
{ |
|
double a = 24371; // NOTE: semi-major axis in km |
|
double mu = 398601.68; // NOTE: gravitational parameter |
|
|
|
double T = orbitGetPeriod(a, mu); |
|
REQUIRE_THAT(T, WithinAbs(37863.5, 50e-3)); |
|
}
|
|
|