| 1 |
//
|
| 2 |
// cOrbit.h
|
| 3 |
//
|
| 4 |
// This is the header file for the class cOrbit. This class accepts a
|
| 5 |
// single satellite's NORAD two-line element set and provides information
|
| 6 |
// regarding the satellite's orbit such as period, axis length,
|
| 7 |
// ECI coordinates/velocity, etc., using the SGP4/SDP4 orbital models.
|
| 8 |
//
|
| 9 |
// Copyright (c) 2002-2003 Michael F. Henry
|
| 10 |
//
|
| 11 |
#pragma once
|
| 12 |
|
| 13 |
#include "cTle.h"
|
| 14 |
#include "cJulian.h"
|
| 15 |
#include "cNoradBase.h"
|
| 16 |
#include "math.h"
|
| 17 |
|
| 18 |
using namespace std;
|
| 19 |
//////////////////////////////////////////////////////////////////////////////
|
| 20 |
|
| 21 |
class cVector;
|
| 22 |
class cGeoCoord;
|
| 23 |
class cEci;
|
| 24 |
|
| 25 |
//////////////////////////////////////////////////////////////////////////////
|
| 26 |
class cOrbit
|
| 27 |
{
|
| 28 |
public:
|
| 29 |
cOrbit(const cTle &tle);
|
| 30 |
virtual ~cOrbit();
|
| 31 |
|
| 32 |
// Return satellite ECI data at given minutes since element's epoch.
|
| 33 |
bool getPosition(double tsince, cEci *pEci) const;
|
| 34 |
|
| 35 |
double Inclination() const { return radGet(cTle::FLD_I); }
|
| 36 |
double Eccentricity() const { return m_tle.getField(cTle::FLD_E); }
|
| 37 |
double RAAN() const { return radGet(cTle::FLD_RAAN); }
|
| 38 |
double ArgPerigee() const { return radGet(cTle::FLD_ARGPER); }
|
| 39 |
double BStar() const { return m_tle.getField(cTle::FLD_BSTAR) / AE;}
|
| 40 |
double Drag() const { return m_tle.getField(cTle::FLD_MMOTIONDT); }
|
| 41 |
double mnMotion() const { return m_tle.getField(cTle::FLD_MMOTION); }
|
| 42 |
double mnAnomaly() const { return radGet(cTle::FLD_M); }
|
| 43 |
double mnAnomaly(cJulian t) const; // mean anomaly (in radians) at time t
|
| 44 |
|
| 45 |
cJulian Epoch() const { return m_jdEpoch; }
|
| 46 |
|
| 47 |
double TPlusEpoch(const cJulian &t) const; // time span [t - epoch] in secs
|
| 48 |
|
| 49 |
string SatName(bool fAppendId = false) const;
|
| 50 |
|
| 51 |
// "Recovered" from the input elements
|
| 52 |
double SemiMajor() const { return m_aeAxisSemiMajorRec; }
|
| 53 |
double SemiMinor() const { return m_aeAxisSemiMinorRec; }
|
| 54 |
double mnMotionRec() const { return m_mnMotionRec; } // mn motion, rads/min
|
| 55 |
double Major() const { return 2.0 * SemiMajor(); } // major axis in AE
|
| 56 |
double Minor() const { return 2.0 * SemiMinor(); } // minor axis in AE
|
| 57 |
double Perigee() const { return m_kmPerigeeRec; } // perigee in km
|
| 58 |
double Apogee() const { return m_kmApogeeRec; } // apogee in km
|
| 59 |
double Period() const; // period in seconds
|
| 60 |
|
| 61 |
protected:
|
| 62 |
double radGet(cTle::eField fld) const
|
| 63 |
{ return m_tle.getField(fld, cTle::U_RAD); }
|
| 64 |
|
| 65 |
double degGet(cTle::eField fld) const
|
| 66 |
{ return m_tle.getField(fld, cTle::U_DEG); }
|
| 67 |
|
| 68 |
private:
|
| 69 |
cTle m_tle;
|
| 70 |
cJulian m_jdEpoch;
|
| 71 |
cNoradBase *m_pNoradModel;
|
| 72 |
|
| 73 |
// Caching variables; note units are not necessarily the same as tle units
|
| 74 |
mutable double m_secPeriod;
|
| 75 |
|
| 76 |
// Caching variables recovered from the input TLE elements
|
| 77 |
double m_aeAxisSemiMinorRec; // semi-minor axis, in AE units
|
| 78 |
double m_aeAxisSemiMajorRec; // semi-major axis, in AE units
|
| 79 |
double m_mnMotionRec; // radians per minute
|
| 80 |
double m_kmPerigeeRec; // perigee, in km
|
| 81 |
double m_kmApogeeRec; // apogee, in km
|
| 82 |
};
|
| 83 |
|