Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

M002-Cartesian State Vectors To Keplerian Orbit Elements

Download as pdf or txt
Download as pdf or txt
You are on page 1of 2

M.Eng.

René Schwarz
Master of Engineering in Computer Science and Communication Systems (M.Eng.)
Bachelor of Engineering in Mechatronics, Industrial and Physics Technology (B.Eng.)

Memorandum № 2

Cartesian State Vectors −→ Keplerian Orbit Elements


This document is subject to the conditions of the Creative Commons Attribution-NonCommercial 4.0 International

Inputs Outputs
License (CC BY-NC 4.0), which can be examined at https://creativecommons.org/licenses/by-nc/4.0/. Other

 cartesian state vectors  a traditional set of Keplerian Orbit Elements


– position vector r(t) [m] – Semi-major axis a [m]
– velocity vector ṙ(t) [ ms ] – Eccentricity e [1]
 standard gravitational parameter µ = GM of the central – Argument of periapsis ω [rad]
body, if different from Sun (G…Newtonian constant of – Longitude of ascending node (LAN) Ω [rad]
m3 – Inclination i [rad]
gravitation [ kg·s 2 ], M …central body mass [kg])
– Mean anomaly M [rad]
licenses may apply for particular contents or supplemental documents/files.
©2017 M.Eng. René Schwarz (rene-schwarz.com), unless otherwise stated.

1 Algorithm
1. Preparations:
[ ]
m2
a) Calculate orbital momentum vector h s :
h = r × ṙ (1)

b) Obtain the eccentricity vector e [1] from


ṙ × h r
e= − (2)
µ ∥r∥
3
with standard gravitational parameter µ = µÀ = 1.327 124 400 18 · 1020 (±8 · 109 ) ms2 for the Sun as
central body.
[ 2]
c) Determine the vector n ms pointing towards the ascending node and the true anomaly ν [rad] with
{ ⟨e,r⟩
T T
arccos ∥e∥∥r∥ for ⟨r, ṙ⟩ ≥ 0
n = (0, 0, 1) × h = (−hy , hx , 0) ν= ⟨e,r⟩ (3)
2π − arccos ∥e∥∥r∥ otherwise.
cbn

2. Calculate the orbit inclination i by using the orbital momentum vector h, where hz is the third component of h:
hz
i = arccos (4)
∥h∥
3. Determine the orbit eccentricity e [1], which is simply the magnitude of the eccentricity vector e, and the eccen-
tric anomaly E [1]:
tan ν
e = ∥e∥ E = 2 arctan √ 2 (5)
Errors, comments or ideas regarding this paper?

1+e
1−e

4. Obtain the longitude of the ascending node Ω and the argument of periapsis ω:
{ { ⟨n,e⟩
nx
arccos ∥n∥ for ny ≥ 0 arccos ∥n∥∥e∥ for ez ≥ 0
Ω= ω = ⟨n,e⟩ (6)
→ https://goo.gl/XB85LZ

2π − arccos ∥n∥
nx
for ny < 0 2π − arccos ∥n∥∥e∥ for ez < 0

5. Compute the mean anomaly M with help of Kepler’s Equation from the eccentric anomaly E and the eccen-
tricity e:
M = E − e sin E (7)
6. Finally, the semi-major axis a is found from the expression
1
a= ∥ṙ∥2
. (8)
2
∥r∥ − µ

Page 1 of 2
version: 2017/10/05 14:55
M.Eng. René Schwarz (rene-schwarz.com): Memorandum Series
Cartesian State Vectors −→ Keplerian Orbit Elements (Memorandum № 2)

2 Constants and Conversion Factors


Universal Constants
Symbol Description Value Source
1 −11 m3
G Newtonian constant of gravitation G = 6.67428(67) · 10 kg·s2
[1, pp. 686–689]

Conversion Factors
Conversion Source
Astronomical Units → Meters 1 AU = 1.495 978 707 00 · 1011 (±3) m [4, p. 370 f.]
Julian Days → Seconds 1 d = 86 400 s [3]
Degrees → Radians 1◦ = 1 ◦ · π
180◦
rad ≈ 0,017453293 rad

3 References
Equations 1 and 8: [2, p. 28]; Eq. 2: [8]; Eq. 3: [9, 12]; Eq. 5: [7, 10]; Eq. 4: [11]; Eq. 6: [6, 9]; Eq. 7: [5, p. 26]; Value for
µÀ : [3].
[1] Mohr, Peter J.; Taylor, Barry N.; Newell, David B.: CODATA recommended values of the fundamental physical constants: 2006. In: Review of
Modern Physics 80 (2): 633–730. American Physical Society, 2008. ISSN: 1539-0756. DOI: 10.1103/RevModPhys.80.633. [→ cited on page 2]
[2] Montenbruck, Oliver; Gill, Eberhard: Satellite Orbits: Models, Methods, Applications. Corrected 3rd printing. Springer, Heidelberg, 2005.
ISBN: 9783540672807. [→ cited on page 2]
[3] NASA/JPL: Astrodynamic Constants. Online available at http://ssd.jpl.nasa.gov/?constants (retrieved 2011/04/26). [→ cited on page 2]
[4] Pitjeva, E.; Standish, E.: Proposals for the masses of the three largest asteroids, the Moon-Earth mass ratio and the Astronomical Unit. In: Celestial
Mechanics and Dynamical Astronomy 103 (4): 365–372. Springer Netherlands, 2009. ISSN: 0923-2958. DOI: 10.1007/s10569- 009- 9203- 8. [→
cited on page 2]
[5] Standish, E. Myles; Williams, James G.: Orbital Ephemerides of the Sun, Moon and Planets. Online available at ftp://ssd.jpl.nasa.gov/pub/
eph/planets/ioms/ExplSupplChap8.pdf (retrieved 2011/05/15). [→ cited on page 2]
[6] Wikipedia (ed.): Argument of periapsis. Online available at http://en.wikipedia.org/w/index.php?title=Argument%5C_of%5C_periapsis&
oldid=450370870 (retrieved 2011/10/29). [→ cited on page 2]
[7] Wikipedia (ed.): Eccentric anomaly. Online available at http://en.wikipedia.org/w/index.php?title=Eccentric%5C_anomaly&oldid=416453169
(retrieved 2011/10/29). [→ cited on page 2]
[8] Wikipedia (ed.): Eccentricity vector. Online available at http://en.wikipedia.org/w/index.php?title=Eccentricity%5C_vector&oldid=443226923
(retrieved 2011/10/29). [→ cited on page 2]
[9] Wikipedia (ed.): Longitude of the ascending node. Online available at http://en.wikipedia.org/w/index.php?title=Longitude%5C_of%5C_the%
5C_ascending%5C_node&oldid=453399307 (retrieved 2011/10/29). [→ cited on page 2]
[10] Wikipedia (ed.): Orbital eccentricity. Online available at http : / / en . wikipedia . org / w / index . php ? title = Orbital % 5C _ eccentricity & oldid =
453223480 (retrieved 2011/10/29). [→ cited on page 2]
[11] Wikipedia (ed.): Orbital inclination. Online available at http://en.wikipedia.org/w/index.php?title=Orbital%5C_inclination&oldid=452407742
(retrieved 2011/10/29). [→ cited on page 2]
[12] Wikipedia (ed.): True anomaly. Online available at http : / / en . wikipedia . org / w / index . php ? title = True % 5C _ anomaly & oldid = 427251318
(retrieved 2011/05/15). [→ cited on page 2]

1 The numbers in parentheses in 6.67428(67)·10−11 are a common way to state the uncertainty; short notation for (6.67428±0.0000067)·10−11 .

Page 2 of 2
version: 2017/10/05 14:55

You might also like