This preview shows pages 1–7. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full 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
 Spring '11

Click to edit the document details