Path: Orbit/RHSOrbit
% High fidelity orbit model right-hand-side.
Includes higher harmonics of the earth's gravitational field and sun and
moon accelerations as options.
- The drag option = drag coefficient * area / mass where the area is in m^2
- The sun option = area / mass
- The moon option is 0 or 1
The gravity model should be unnormalized.
--------------------------------------------------------------------------
Form:
xDot = FOrbHF( x, t, options, jD, aExt, mu, s, c, j, a )
xDot = FOrbHF( x, t, options, jD, aExt, d )
xDot = FOrbHF( t, x, options, jD, aExt, mu, s, c, j, a ) (for ode113)
--------------------------------------------------------------------------
------
Inputs
------
x (6,1) The state vector [x;y;z;dx/dt;dy/dt;dz/dt]
t (1,1) Time
options (1,5) [Highest zonal harmonic,Highest tesseral harmonic,
Sun,Moon,Drag]
jD Julian date
aExt (3,1) External acceleration vector [x;y;z] ECI
d (.) Gravity model struct, unnormalized
- or -
mu Spherical gravitational potential
s (36,36) S terms
c (36,36) C terms
j (36) m = 0 terms
a Earth radius
-------
Outputs
-------
xDot (6,1) The derivative of the state state vector
--------------------------------------------------------------------------
See also AGravity, APlanet, AtmDens2, SolarF, and FOrb.
--------------------------------------------------------------------------