Long period perturbations of earth satellite orbits

All the equations involved in extending the PS phi solution to include the long periodic and second order secular effects of the zonal harmonics are presented. Topics covered include DSphi elements and relations for their conconical transformation into the PS phi elements; the solution algorithm based on the Von Zeipel method; and the elimination of long periodic terms and analytical integration of primed variables. The equations were entered into the ASOP program, checked out, and verified. Comparisons with numerical integrations show the long period theory to be accurate within several meters after 800 revolutions.