[发明专利]有势位力作用下的导航卫星轨道确定改进的无奇点技术

权利要求书

1.导航卫星轨道确定改进的无奇点技术，其特征在于，包括如下步骤：

步骤1 ：定义和求取摄动位相对于六个开普勒根数之一（平近点角）的偏导数，即六个中间轨道变量之一的时间导数，即

                                                                                                                     (1)

步骤2 ：定义和求取摄动位相对于六个开普勒根数之二（近地点角）的偏导数，即六个中间轨道变量之二的时间导数，即

                                                                      (2)

步骤3 ：定义和求取摄动位相对于六个开普勒根数之三（扁率）的偏导数，即六个中间轨道变量之三的时间导数，即

                                                                      (3)

步骤4 ：定义和求取摄动位相对于六个开普勒根数之四（升交点赤经）的偏导数，即六个中间轨道变量之四的时间导数，即

                                                                     (4)

步骤5 ：定义和求取摄动位相对于六个开普勒根数之五（倾角）的偏导数，即六个中间轨道变量之五的时间导数，即

                                                                    (5)

步骤6 ：定义和求取摄动位相对于六个开普勒根数之六（半长轴）的偏导数，即六个中间轨道变量之六的时间导数，即

                                                                    (6)

其中R 是摄动力的位函数，方程右侧为位函数相对于六个开普勒根数 (a, e, ω, i, ?, M，即半长轴，扁率，近地点角，倾角，升交点赤经，平近点角) 的偏导数，左侧是六个中间变量（用开普勒根数带下标1表示）对时间的导数

步骤7 ：根据步骤5中的推导，判断偏导数是否含有因子 sini

记为

步骤8 ：根据步骤3中的推导，判断偏导数是否含有因子 e 

记为

步骤9 ：根据初始轨道根数，判断，即轨道是否为近赤道轨道

步骤10 ：根据初始轨道根数，判断，即轨道是否为近圆形轨道

步骤11 ：根据步骤1组成如下第一个无奇点方程，即

     ,                                                                       (7)

步骤12 ：根据步骤1和步骤2，组成如下第二个无奇点方程，即

  ,                                                     (8)

步骤13 ：根据步骤3和步骤5，以及步骤7-10，组成如下第三个无奇点方程，即

  ,    (9)

步骤14 ：根据步骤2和步骤4，组成如下第四个无奇点方程，即

  ,                                           (10)  

步骤15 ：根据步骤5，以及步骤7和步骤9，组成如下第五个无奇点方程，即

  ,                                 (11)

步骤16 ：根据步骤3和步骤6，以及步骤8和步骤10，组成如下第六个无奇点方程，即

  .                            (12)

步骤17 ：根据积分中值定理积分上述（7-12）式，即得到卫星受有势位力摄轨道的理论表达式

步骤18 ：根据卫星观测数据，对上述结果进行拟合平差滤波，即获得有势位力作用下的轨道确定或轨道改进结果。