0% found this document useful (0 votes)
177 views

M002-Cartesian State Vectors To Keplerian Orbit Elements

This document describes an algorithm to convert cartesian state vectors describing the position and velocity of an object in orbit to Keplerian orbital elements. The algorithm uses vectors related to the orbital momentum, eccentricity and position to calculate the semi-major axis, eccentricity, inclination, longitude of ascending node, argument of periapsis and mean anomaly.

Uploaded by

Faisal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
177 views

M002-Cartesian State Vectors To Keplerian Orbit Elements

This document describes an algorithm to convert cartesian state vectors describing the position and velocity of an object in orbit to Keplerian orbital elements. The algorithm uses vectors related to the orbital momentum, eccentricity and position to calculate the semi-major axis, eccentricity, inclination, longitude of ascending node, argument of periapsis and mean anomaly.

Uploaded by

Faisal
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 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

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy