Tracking Earth-Orbiting Satellites

Coordinate Transformations Orbital elements offer an excellent means of keeping track of earth orbiting satellites. However, in order to find the satellite in the sky ant any given time, these elements need to be transformed into a Cartesian coordinate system. his transformation is relatively simple once we can understand the motion of the satellite geometrically. To start with the satellite moves in a plane which goes through the center of the earth. The direction of this plane is determined by three angles, the argument of thw perigee , the right ascension of the node , and the angle of inclination .We start with a Cartesian coordinate system with coordinates in which the coordinate is along the semi-major axis of the ellipse pointing toward the perigee, - the -coordinate is perpendicular to the axis in the plane of motion, and the - coordinate is normal to this plane. The origin is fixed at the center of the earth. Using Kepler's equation we get the eccentric anomaly from the mean anomaly and can give the Cartesian coordinates of the satellite in the coordinate system , : In the next step we transform these Cartesian coordinates into a new Cartesian coordinate system , where is along the ascending node, is perpendicular to in the plane perpendicular to the axis of the earth, and is along the axis. This is a simple rotation about the axis by the negative angle of the right ascension of the perigee . It is clear, that . For the and coordinates, a little trigonometry (see Figure x, below) gives: Figure 4: Rotation of Coordinate systems We can make this more efficient by describing it via the rotation matrix If introduce the notation the coordinate transformation as and , we can write In linear algebra, one studies these rotation matrices in detail. A rotation by an angle about the -axis is given by If the rotation is about the -axis the matrix is and a similar formula holds for a rotation about the -axis. To continue our coordinate transformation we need to rotate the about the -axis (which is also the line of nodes) by an angle of system , this will give us the coordinates in a Cartesian coordinate system in which the axis is along the axis of the earth and the axis is along the line of nodes. Following our discussion above this transformation is given by where The final step is to rotate the system about the -axis by an angle of , this will give us the desired coordinates in the Cartesian coordinate system, in which the -axis points to the vernal equinox. This transformation is again given by a matrix in particular it is given by where . We can now combine these three rotations into a single matrix , where the product on the right is the usual matrix multiplication. The coordinate transformation is now: While this is an elegant way of writing the process down mathematically, in practice it is easier doing the three steps one at a time.