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

FindPosition SGP4 changes

feature/19
Daniel Warner 2011-03-29 20:30:46 +01:00
parent e2bffd6fb3
commit 204cc4d4b9
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
SGDP4.cpp
View File

@ -309,21 +309,27 @@ void SGDP4::FindPosition(double tsince) {
final_aycof = 0.25 * i_a3ovk2_ * final_sinio; final_aycof = 0.25 * i_a3ovk2_ * final_sinio;
#endif #endif
} else { } else {
xincl = Inclination(); xincl = Inclination();
omega = omgadf; omega = omgadf;
double xmp = xmdf; double xmp = xmdf;
if (!i_use_simple_model_) { if (!i_use_simple_model_) {
const double delomg = i_omgcof_ * tsince; const double delomg = i_omgcof_ * tsince;
const double delm = i_xmcof_ * (pow(1.0 + i_eta_ * cos(xmdf), 3.0) - i_delmo_); const double delm = i_xmcof_ * (pow(1.0 + i_eta_ * cos(xmdf), 3.0) - i_delmo_);
const double temp = delomg + delm; const double temp = delomg + delm;
xmp = xmdf + temp;
xmp += temp;
omega -= temp; omega -= temp;
const double tcube = tsq * tsince; const double tcube = tsq * tsince;
const double tfour = tsince * tcube; const double tfour = tsince * tcube;
tempa -= i_d2_ * tsq - i_d3_ * tcube - i_d4_ * tfour; tempa -= i_d2_ * tsq - i_d3_ * tcube - i_d4_ * tfour;
tempe += BStar() * i_c5_ * (sin(xmp) - i_sinmo_); tempe += BStar() * i_c5_ * (sin(xmp) - i_sinmo_);
templ += i_t3cof_ * tcube + tfour * (i_t4cof_ + tsince * i_t5cof_); templ += i_t3cof_ * tcube + tfour * (i_t4cof_ + tsince * i_t5cof_);
} }
a = RecoveredSemiMajorAxis() * pow(tempa, 2.0); a = RecoveredSemiMajorAxis() * pow(tempa, 2.0);
e = Eccentricity() - tempe; e = Eccentricity() - tempe;
xl = xmp + omega + xnode + RecoveredMeanMotion() * templ; xl = xmp + omega + xnode + RecoveredMeanMotion() * templ;