Details

Type: Bug

Status: Open (View Workflow)

Priority: P2

Resolution: Unresolved

Affects Version/s: R2018A

Fix Version/s: Someday

Component/s: Nav

Labels:
Description
To help solve the Doppler noise issue, GMAT is performing many calculations with more precision.
Determine if changes need to be made to the matlab prototypes. Implement easy enhancements.
It may not be easy to get more matlab precision for a propagation  investigate. (The work estimate for this ticket can be reduced if we transfer the matlab propagation work to a different Sprint)
Gliffy Diagrams
Attachments
Issue Links
 relates to

GMT6195 Matlab routines for Precision Earth rotation  Part I
 Closed
Activity
 All
 Comments
 Work Log
 History
 Activity
 Subversion
Note from John on matlab routine sidereal.m: It looks like Vallado provides two different methods to calculate the rotation. I don’t have a copy of the book in front of me, but I believe that sidereal.m uses equation 343, while GMAT uses equation 345 (2nd edition). I could change sidereal.m to use equation 345 like we use in GMAT, but the epoch is only given in 1 word time, while the changes I made require 2 word time for the extra needed precision.
Note from John on what the GMAT code change:
In AxisSystem::ComputeSiderealTimeRotation() I changed the following:
Real ThetaGmst = ((67310.54841 * sec2deg) +
((876600 * hour2deg) + (8640184.812866 * sec2deg))*tUT1 +
(0.093104 * sec2deg)*tUT12  (6.2e06 * sec2deg)*tUT13 )
Became:
Real ThetaGmst = ((67310.54841 * sec2deg) +
(mjdUT1.GetSec() + mjdUT1.GetFracSec())/SECS_PER_DAY*TWO_PI_DEG +
(8640184.812866 * sec2deg)*tUT1 +
(0.093104 * sec2deg)*tUT12  (6.2e06 * sec2deg)*tUT13 )