The orbital variation of a lunar satellite is somewhat different from an artificial Earth satellite, due to slow rotation of the Moon and the features of lunar gravity field. In this paper, the above problems are studied thoroughly ; all disturbing factors are investigated in detail; and an analytical solution, on the premise of certain accuracy, is given by the method of quasi-average elements. Comparing with numerical integration shows that the analytical solution is able to describe the motion of lunar satellite with high precesion.
展开▼