KJ135
/
sgp4
mirror of http://172.16.0.12:5000/LAMPSARToolSoftware/sgp4.git
4
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

Changes to mean motion / semi major axis recovery

feature/19
Daniel Warner 2011-03-29 17:52:34 +01:00
parent 7eb83ea9e9
commit 09f9da3bf9
1 changed files with 6 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
SGDP4.cpp
View File

@ -91,12 +91,13 @@ void SGDP4::SetTle(const Tle& tle) {
double eosq = Eccentricity() * Eccentricity();
double betao2 = 1.0 - eosq;
double betao = sqrt(betao2);
double del1 = 1.5 * constants_.CK2 * i_x3thm1_ / (a1 * a1 * betao * betao2);
double ao = a1 * (1.0 - del1 * (0.5 * constants_.TWOTHRD + del1 * (1.0 + 134.0 / 81.0 * del1)));
double delo = 1.5 * constants_.CK2 * i_x3thm1_ / (ao * ao * betao * betao2);
double temp = (1.5 * constants_.CK2) * i_x3thm1_ / (betao * betao2);
double del1 = temp / (a1 * a1);
double a0 = a1 * (1.0 - del1 * (1.0 / 3.0 + del1 * (1.0 + del1 * 134.0 / 81.0)));
double del0 = temp / (a0 * a0);
recovered_mean_motion_ = MeanMotion() / (1.0 + delo);
recovered_semi_major_axis_ = ao / (1.0 - delo);
recovered_mean_motion_ = MeanMotion() / (1.0 + del0);
recovered_semi_major_axis_ = a0 / (1.0 - del0);
/*
* find perigee and period