lecture06 - Computation of orbital elements from Cartesian...

Info iconThis preview shows pages 1–7. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: Computation of orbital elements from Cartesian Ref.: Curtis Sec. 4.4, Vallado Section 2.5 We know about orbital elements, but how do we find them? If were given Cartesian state vector (position vector r and velocity vector r ), we can compute them using the conserved quantities. L. Healy ENAE404 Spring 2007 Lecture 6 (Feb. 13) 1 IJK to Classical Elements To compute the classical elements a, e, i, , , ( t ) from the Cartesian state vectors r , r , we will need to know how to take dot and cross prod- ucts in Cartesian coordinates. Start with angular momentum vector h = r r semilatus rectum p = h 2 / eccentricity vector e = r h - r node direction n = K h . K is the unit vector along the polar axis. L. Healy ENAE404 Spring 2007 Lecture 6 (Feb. 13) 2 Semimajor axis a and eccentricity e Eccentricity can be computed as magnitude of the eccentricity vector, e = | e | . The semimajor axis can be computed in one of two ways: a = 2 r- v 2 !- 1 = p 1- e 2 . Now we need to find the angular elements i , , , . L. Healy ENAE404 Spring 2007 Lecture 6 (Feb. 13) 3 Finding angular elements After finding the unit vectors compute angles, RAAN: = arctan( n J , n I ) or cos = n I and if n J < 0, the angle is in the lower halfplane < < 2 . Inclination: angular momentum direction gives inclination, since its always perpen- dicular to the orbital plane: cos i = h K No need for arctangent here, 0 i 180 . Argument of perigee is the angle from the ascending node to the eccentricity vector, cos = n e If the third component of e is negative, e K < 0, then the argument of perigee is below the equatorial plane, 2 . L. Healy ENAE404 Spring 2007 Lecture 6 (Feb. 13) 4 Finding angular elements contd. True anomaly cos = r e = 1 e p r- 1 If r r < 0, then the true anomaly is below the major axis, < 2 , because (recall r r = r r ) from perigee to apogee 0 < < , dis- tance r is positive and increasing, so r r = r r > from apogee to perigee < 2 , dis- tance r is decreasing, so r r = r r < 0. L. Healy ENAE404 Spring 2007 Lecture 6 (Feb. 13) 5 Coordinate systems We have computed the orbital elements from the Cartesian IJK state vector. What about going the other direction? It is possible to compute the inverse transformation, but we shall need to discuss coordinate systems and transformations first. So far we have used the coordinate system IJK but have skirted what coordinate systems are about and how to get from one to another (transformation). This requires some atten- tion, because coordinates are an essential part of everything we do in astrodynamics....
View Full Document

Page1 / 24

lecture06 - Computation of orbital elements from Cartesian...

This preview shows document pages 1 - 7. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online