A long-term evolution rapid
analysis method for an inclined
geosynchronous orbit comprises the following steps: firstly, deducing a primary or secondary average perturbation
force function of each perturbation item, including 1:1
resonance parts of earth non-spherical perturbation field
harmonic items J22, J31, J32, J33, J41, J42, J43 and J44 with
harmonic items, sun-moon
gravitation perturbationLegendre expansion intercepting 4-order items and
solar light pressure perturbation; secondly, establishing a secondary average semi-analytical
orbit recursor of the IGSO
orbit through a Lagrange typeorbit perturbation equation in combination with an average perturbation
force function; comparing and analyzing IGSO
orbit long-term evolution conditions under consideration of different perturbationsources and orders, so that the model is simplified, and the orbit
recursion efficiency is further improved; and finally, by means of efficient orbit
recursion, drawing a series of dynamic grid diagrams of complete orbit elements and initial epoch moment combinations, wherein the divisible orbit elements are (e, i, [
Omega], [
omega]), all binary combinations of the divisible orbit elements are e-i, [
omega]-[
Omega];e-[
omega], i-[
Omega], e-[Omega], i-[omega], all binary combinations of the initial epoch moment are InitialEpoch-e, InitialEpoch-i, InitialEpoch-[Omega], and InitialEpoch-[omega], and rapid and comprehensive analysis of IGSO orbit long-term evolution is completed according to the dynamic grid chart.