The drawing on the right gives a geometric construction of those angles (no, don't try to puzzle out the details). The orbital ellipse is enclosed in a circle of radius a, and given a position P of the satellite, a corresponding point Q on the circle can be drawn, sharing the same line perpendicular to the ellipse's axis. Then E is the angle between the long axis of the ellipse and the line drawn from the center of the circle to Q ("eccentric" might mean here "from the center").
Suppose the elements a, e and M(0) at time t=0 are given, and we need to find the value of φ at some different time t. With f known, the above equation gives r, and (r ,φ) together pin-point the satellite's position in its orbital plane. The first step is to derive
M = M(0) + 360°(t/T)
We assume the period T is known (this requires the 3rd law and is discussed for circular orbits in sections 20 and 20a). It can then be shown that the angle E satisfies "Kepler's equation"
M = E – (180°/π)e sinE
where π = 3.14159256... is the ratio between the circumference of a circle and its diameter. How did that number suddenly crop up, you may ask? The fact is, the division of the circle into 360 degrees may be convenient to use (we inherited it from the ancient Babylonians) but the number 360 has no particular place in mathematics. It is probably related to the number of days in a year. The "natural" division of angles which arises in calculus and other branches of math is into radians, with 360 degrees equal to 2π = 6.2831... radians (making each radian equal to about 57.3 degrees). With angles measured in radians, Kepler's equation simplifies to
M = E – e sinE
No matter which form is used, mathematics knows no formula which gives E in terms of M. However, solutions can often be approximated to any degree of accuracy by iteration--by starting with an approximate solution, then improving it again and again by an appropriate procedure ("algorithm"--more about that word, here). If the eccentricity e is not too big--the ellipse not much different from a circle--then M and E are not too different. So an initial guess
E' = M
may not be too far off. Putting this guess into the term sinE gives an improved guess E"
E" = M + (180°/π)e sinE'
One can now insert E" in the sinE term and get an even closer guess, and so on and so forth... until the first (say) ten decimal digits of the value of E no longer change, at which point we may decide we have E to sufficient accuracy and stop the process. Computer handle such a process of continuous improvement ("iteration of the solution"--one form of an algorithm) very rapidly, and other methods also exist, with sufficient speed even when e is not very small.
Given E, a number of formulas will give the true anomaly φ. For instance, one can first derive
r = a(1 – e cos E)
And then cos φ can be found from
r = a(1 – e2)/(1 + e cos φ)
and sin φ follows from cos φ. All this is easily and automatically computed nowadays, but must have been a real hassle in the days before computers.
The Orbit in Space
The 3 remaining orbital elements are all angles giving the position of the orbit in 3 dimensions. They are described below, but their actual use belongs to a university course in orbital mechanics and will be omitted. The angles:
- Inclination i.
The argument of perigee ω (small Greek omega).
The longitude of the ascending node Ω (capital omega).
To orient the orbit in 3 dimensions requires a reference plane and a reference direction. For satellite orbits, the reference plane--the horizontal plane in the drawing--is usually the Earth's equatorial plane (sometimes it is the plane of the ecliptic). The reference direction in either case is the direction from the center of the Earth to the vernal equinox (which belongs to both above planes). We will call it the x direction, since that is its role in (x,y,z) coordinates used in orbital calculations.
Two non-parallel planes always intersect along a line--the way the plane of a door intersects the plane of the wall along the door's hinge. The orbital plane and the equatorial plane (used for reference) do so too, and their intersection is called the line of nodes N. Let the origin O of our coordinates be the center of the Earth, which is also the focus of the ellipse; this point belongs to both the equatorial plane and the orbital plane, and is therefore also on their intersection line N (drawing).
- The inclination i is the opening angle of the "hinge" along N. It is best defined by erecting at O lines perpendicular to each plane and measuring the angle between them (drawing).
- The angle Ω is measured in the equatorial plane between N and the reference direction x. One can imagine rotating the "hinge" N around point O, without changing the inclination: the orbital plane then covers all possible values of Ω.
But what is this "ascending node" bit? The above definition contains some ambiguity: N defines two lines coming out of O, in opposite directions. From which of them should Ω be measured? To resolve this we note that the plane of the equator divides space into two parts, one north of it and one south of it. Specifying "the ascending node" selects the branch on which the satellite crosses as it enters the northern half-space, rather than the one crossed when leaving it.
- Finally, ω is the angle measured in the orbital plane between N and the direction from O to the perigee point P. If perigee lies on the "hinge", on the side of positive x, then ω = 0; rotating the orbit by 90o until the line OP is perpendicular to N gives ω = 90o, rotating it further until it reaches the negative side of the x-direction gives ω = 180o.
Suppose you have the orbital elements of some satellite, e.g. the space shuttle (you can often get them off the world-wide web). The first three (a, e, M), with M given at some particular time, enable you to calculate where the satellite will be at any time in its orbit. With (i, ω, Ω) you can then find where it would be in the sky.
[For those who would like to know more about the way such rotations in space are calculated, see problem 8 in the "trigonometry proficiency drill." That problem only deals with rotations in two dimensions (e.g. on a sheet of paper), but an additional short section linked to its end explains how this is extended to 3D.]
Questions from Users:
Launch window of the Space Shuttle
*** Motion of the Moon