Coordinate SystemsASEN 3200 NotesThe transformation matrix between these frames isWhere ,,and are unit vectors along the P, Q, W and X, Y, Z axes respectively and the etc. are direction cosines.Thus, we have directly, i.e.,Coordinate SystemsASEN 3200 NotesGeorge H. BornSeveral Coordinate Systems are used extensively in Orbit Mechanics. They will be reviewed here.1. Earth Centered Interval (ECI) and Earth Centered Fixed (ECF)The ECI frame has its origin at the center of mass of the Earth but has a fixed inertial direction along the intersection of the Earth equatorial plane and the ecliptic plane. Although this frame is referred to as inertial it is actually only pseudo inertial because the center of mass of the Earth accelerates due to perturbations from the Moon and other planets. The ECF frame has the same origin but is fixed in the Earth with its X- axis through the Greenwich meridian (zero longitude). g Greenwich Siderial Time g = )(00ttg(1) 0g Greenwich hour angle (Greenwich Siderial Time) at the epoch,0t.The transformation matrix between these frames is 10000ggggECICSSCZYXECFZYX. (2)1ECIZZ,ECFYYgXX2/6/20012. Another useful frame is the P,Q,W frame where Pand Q are unit vectors in the orbit plane with Pdirected to perigee, Walong the angular momentum vector and Qcompleting a right hand triad, i.e., WQP To go from XYZ (ECI) to P,Q,WWe first rotate about Z through , the right ascention of the ascending node. Next rotate about the new X-axis through I, the inclination,then rotate about the resulting Z-axis through , the argument of perigee.Here P, Q, W are the components of a vector expressed in the PQW frame, i.e., PQWr = WWQQPP . (3)Remember P,Q,Ware unit vectors. These three rotations will yield the desired transformation matrix as the product of three matrices, i.e.,WQP = 10000CSSCIIIICSSC0000110000CSSCZYX(4)or,WQP = IIIIIIICCSSSSCCCCSSSCCCSSSCCSSCSCSCCIIZYX(5)Hence,2S/C Orbit PlaneYZQWXPjikEquatorIECIPQW = IIIIIIICCSSSSCCCCSSSCCCSSSCCSSCSCSCCII.Then,WQP = ECIPQWZYXandZYX = TECIPQWWQP, (6)and PQWECI = TECIPQW.If we have unit vectors P,Q,Wexpressed in the ECI frame (as we obtain using Gibbs Method for example), thenkZjYiXPPPP,kZjYiXQQQQ,kZjYiXWWWW,and WQP = WWWQQQPPPZYXZYXZYXkji. (7)Where P,Q,Wand kji ,, are unit vectors along the P, Q, W and X, Y, Z axes respectively and the PPPZYX ,,etc. are direction cosines.Thus, we have ECIPQW directly, i.e.,ECIPQW=WWWQQQPPPZYXZYXZYX. (8)Also, note that P lies along the eccentricity vector and Wlies along the angular momentum vector. Consequently, if the position and velocity vectors are available in the ECI frame. We may obtain P,Q,Wdirectly as follows3rrhreECI1, (9)eeP , (10)rrhECI, (11)hhW , (12)andPWQ since P,Q,Wform a right hand triad.(13)3. Another frame often used is the RIC frame with unit vectors IR,and C.In this frame Rlies along the instantaneous radius vector, Ilies in the orbit plane normal to Rand in the direction of motion of the spacecraft. Cis normal to the orbitplane and lies along the angular momentum vector. Hence, WC of the PQW frame.Rand Ialign with Pand Qwhen the spacecraft is at perigee. Although Rand Irotate with the spacecraft radius vector, at each instant in time the frame is consideredfixed, so we do not differentiate these unit vectors when transforming velocity to this frame i.e., velocity magnitude in this frame has the same value as in the inertial, ECI, frame. This frame is useful for displaying the difference between two orbits in the radial, in-track, and cross-track directions. To obtain the transformation between the ECI and RIC frames, assume we are given the position and velocity vectors in the ECI frame. ThenkZjYiXr , (14)kZjYiXrECI, (15)rrhECI, (16) and krZjrYirXrrR , (17)rrrrhhCECIECI, (18)RCI . (19)Or in matrix notation4CIR = zyxzyxzyxCCCIIIRRRkji. (20)DefineECIRIC = zyxzyxzyxCCCIIIRRR,where the elements of ECIRICare the direction cosines between the ECI and RIC frames and are given by Eqns (17), (18), and (19).Hence, the coordinates in one frame can be computed from coordinate in the other, i.e.,CIR = ECIRICZYX andZYX = TECIRICCIR. (21)If we wish to examine the difference between two orbits in the IR,and Cdirections, we do the following. Designate one orbit as the reference orbit. Call the position and velocity in this orbit r and rECI. Then use equation (17), (18), (19) and (20) to compute the unit vectors CIR ,,and ECIRICusing r and rECI. Next compute position and velocity difference in the ECI frame, ECIECIECIECIECIrrrZYXrrrZYX,. (22)Then use equation (21) to compute differences in the radial, in-track and cross-track directions in position and velocity
View Full Document