| 1 |
kusanagi |
1.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 |
|
|
|
Parent Directory
| 
Revision Log