1.地球同步卫星小推力南北保持控制解析计算方法，其特征在于，具体计算过程包括如下步骤：

步骤1、解析计算地球同步卫星倾角矢量的长期摄动变化量；

步骤2、计算电推力器点火序列，解析计算出整个规划周期内每一次电推点火的开关机时刻。

2.如权利要求1所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，步骤1的具体步骤为：确定倾角矢量为(ix iy)，其中ix＝icosΩ，ix＝isinΩ，i和Ω分别为地球同步卫星轨道要素的倾角和升交点赤经；计算规划时段t0到tf内，倾角矢量的长期摄动变化量(Δix)长期和(Δiy)长期。

3.如权利要求2所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤1具体按照以下步骤：

步骤1.1、计算地球非球形引力对倾角矢量的摄动量；

步骤1.2、计算太阳引力对倾角矢量的摄动量；

步骤1.3、计算月球引力对倾角矢量的摄动量；

步骤1.4、将步骤1.1至步骤1.3分别求得的地球非球形引力、太阳引力、月球引力对倾角矢量的摄动量求和，即为地球同步卫星倾角矢量的长期项摄动变化量。

4.如权利要求3所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤1.1中地球非球形引力对倾角矢量的摄动量计算公式为：其中，n为卫星轨道角速度，a为卫星轨道半长轴，μe为地球引力常数，Re为地球半径，J2、J22、J31、J33、J4、J43、λ22、λ31、λ33、λ42、λ43是地球非球形四阶主项摄动系数，ix、iy为倾角矢量初值。

5.如权利要求4所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤1.2中太阳引力对倾角矢量的摄动量计算公式为：其中，μs为太阳引力常数，rs为日地距离，Ωs为地心惯性系下的太阳轨道升交点赤经，is为太阳轨道倾角。

6.如权利要求5所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤1.3中月球引力对倾角矢量的摄动量计算公式为：其中，μl为月球引力常数，rl为地月距离，Ωl为月球轨道升交点赤经，il为月球轨道倾角。

7.如权利要求6所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤1.4中地球同步卫星倾角矢量的长期项摄动变化量由上述三个部分组成，即(Δix)长期＝(Δix)非球形，长期+(Δix)太阳，长期+(Δix)月球，长期                       (7)(Δiy)长期＝(Δiy)非球形，长期+(Δiy)太阳，长期+(Δiy)月球，长期                       (8)。

8.如权利要求7所述的地球同步卫星小推力南北保持控制解析计算方法，其特征在于，所述步骤2的具体过程如下：

步骤2.1、根据步骤1中地球同步卫星倾角矢量的长期项摄动变化量得到一个周期南北保持的倾角矢量目标控制量；一个周期南北保持的倾角矢量目标控制量能够抵消周期内的长期项摄动量，即

(Δiy)control＝‑(Δiy)长期                                    (9)(Δix)control＝‑(Δix)长期                                      (10)；

步骤2.2、计算每次小推力点火时长，得到规划周期内的倾角控制量；

根据步骤2.1中一个周期南北保持的倾角矢量目标控制量得到控制方向角，控制方向角αi＝arctan2[(Δiy)control,(Δix)control]，南侧电推力器的点火中心赤经为InnerTwoPi(αi)，北侧电推力器的点火中心赤经为InnerTwoPi(αi+π)；

点火序列规划周期为N天，其中的M天进行电推点火，M≤N，为轨道确定、载荷控制预留空间，每天有2次南北控制的时机，电推点火总次数为2M，计算每次小推力点火时长为Tfire＝2arcsin[(Δi)controlVsωe/acc/M/4]/ωe                  (11)其中，(Δi)control为规划周期内的倾角控制量，acc为电推力器法向合加速度，Vs为地球同步卫星标称轨道速度，ωe为地球自转角速度；

电推力器一天内的两次点火时刻t1,on、t1,off和t2,on、t2,off分别由南侧和北侧点火中心赤经和点火时长计算出来，设当天00:00:00.00时刻的格林尼治恒星时角为Θ0，卫星平经度为λ；

t1,on＝InnerTwoPi(αi‑λ‑Θ0)/ωe‑Tfire/2                        (13)t1,off＝t1,on+Tfire                                  (14)t2,on＝InnerTwoPi(αi+π‑λ‑Θ0)/ωe‑Tfire/2                      (15)t2,off＝t2,on+Tfire                                  (16)计算出来的t1,on、t1,off和t2,on、t2,off是从当天00:00:00.00时刻开始的积秒，以此类推，解析计算出整个规划周期内每一次电推点火的开关机时刻。