非保守力作用下的导航卫星轨道确定改进的无奇点技术的制作方法
【专利摘要】本发明涉及一种非保守力下的导航星轨道确定改进的无奇点技术。技术特征在于：通过高斯方程获得非保守力摄动的理论表达。根据高斯--许方程，组成受摄轨道无奇点的理论表达，根据观测进行理论拟合，从而确定导航星轨道或轨道改进。因北斗导航星轨道都是近圆和赤道的，经典轨道理论均存在奇点，从而定轨或轨道改进失败。本发明技术可供北斗导航星轨道定轨和改进，而没有奇点。另本发明技术简单，可减少计算量，适合星上自行定轨。因卫星摄动力包括有势位力与非保守力两种，故本发明与发明权利登记号201310545077.2所述发明在有势位力下的导航星轨道确定改进的无奇点技术联合应用，构成了一个完整的无奇点轨道确定改进技术。
【专利说明】 非保守力作用下的导航卫星轨道确定改进的无奇点技术
【技术领域】
[0001]本发明涉及一种利用无奇点理论和技术确定非保守力作用下的卫星轨道和对卫星轨道进行改进的方法和技术，涉及卫星导航定位【技术领域】，涉及导航卫星测量定轨领域，属于大地测量与卫星导航定位学科，属于天体测量与天体力学学科。
【背景技术】
[0002]非保守力作用下的导航卫星的轨道确定和改进的关键问题在于卫星运动方程应用。
[0003]目前主要描述导航卫星运动的方程有两类:1)直角坐标系下的运动方程，该方法无法获得理论解即理论的轨道描述，所以只能用数值积分求解，需要积分器，需要较大的计算能力，很难实时，无法对运动进行定性分析，无法星上自行定轨应用。2)开普勒变量下的高斯方程，应用该方程可以求得卫星轨道的理论描述，定轨和轨道改进变成拟合平差滤波，计算量大大减小，适合于星上自行定轨应用，可以实时，但是轨道理论描述存在奇点问题。
[0004]GPS导航卫星的轨道基本上都是近圆形的轨道，所以卫星轨道有近圆奇点问题。北斗卫星导航系统除了近圆轨道卫星外，还有赤道轨道卫星，既存在近圆奇点问题，也存在赤道轨道奇点问题。
[0005]上述现有理论方法技术要么存在计算量大，需要积分器，很难实时，无法定性分析的缺点，要么存在奇点问题，无法用于实时导航卫星星上的自行定轨。
【发明内容】
[0006]要解决的技术问题
为了使得能够自行星上定轨和轨道改进，本发明提出一种无奇点的技术能使非保守力作用下的导航卫星轨道的理论描述不含奇点，从而使得到导航卫星轨道的确定和改进都可以获得理论描述而不会发生奇异，从而通过卫星观测拟合平差滤波，实现实时星上自行定轨。利用本发明进行卫星定轨和轨道改进可以实现北斗系统的自主能力，减少对地面控制站的依赖，从而提高北斗系统在实战中的生存能力，从而提高北斗应用的范围，并具有计算简单、精度高，成本降低等优点。综合本发明与发明登记号201310545077.2的发明，全部受摄运动的卫星轨道均无奇点问题，从而形成了一整套可实时星上自行定轨的无奇点技术。
[0007]技术方案
本发明的技术特征在于:根据不含奇点的六个中间变量的时间导数，作出判断，决定如何组成在非保守力作用下的无奇点的运动方程，然后利用中值积分定理求得轨道的理论描述，并通过对卫星的观测拟合平差滤波获得轨道确定和改进，具体步骤如下:
步骤1:求取非保守力作用下第一个中间变量对时间的导数；
步骤2:求取非保守力作用下第二个中间变量对时间的导数；
步骤3:求取非保守力作用下第三个中间变量对时间的导数；
步骤4:求取非保守力作用下第四个中间变量对时间的导数； 步骤5:求取非保守力作用下第五个中间变量对时间的导数；
步骤6:求取非保守力作用下第六个中间变量对时间的导数；
步骤7:判断(5)式是否含有因子sini ;
步骤8:判断(3)式是否含有因子
步骤9:根据初始轨道根数判断轨道倾角是否为零，即轨道是否为赤道轨道；
步骤10:根据初始轨道根数判断轨道椭圆扁率是否为零，即轨道是否为圆形轨道； 步骤11:根据(7)式组成开普勒根数轨道半长轴相对于时间的导数的微分方程；
步骤12:根据(8)式组成开普勒根数轨道轨道椭圆扁率相对于时间的导数的微分方
程;
步骤13:根据(9)式组成开普勒根数近地点角相对于时间的导数的微分方程；
步骤14:根据(10)式组成开普勒根数轨道轨道面倾角相对于时间的导数的微分方
程;
步骤15:根据(11)式组成开普勒根数轨道升交点赤经相对于时间的导数的微分方
程;
步骤16:根据(12)式组成开普勒根数轨道平近点角相对于时间的导数的微分方程；步骤17:根据中值积分定理积分方程(7-11)从而获得非保守力作用下无奇点的轨道理论描述；
步骤18:根据对卫星的观测和对轨道理论描述进行拟合平差滤波，从而确定非保守力作用下的卫星轨道或其改进；
步骤19:根据发明权利登记号201310545077.2所发明的保守位力作用下的导航卫星轨道确定改进的无奇点技术，可获得在保守力位作用下轨道确定或轨道改进的结果；
步骤20:综合上述步骤18和步骤19所获得的结果，即本发明与发明权利登记号201310545077.2所述发明保守力作用下的导航卫星轨道确定改进的无奇点技术的联合应用，构成了一个完整的无奇点轨道确定改进技术。
[0008]有益效果
本发明提出的利用无奇点理论推导非保守力作用下轨道理论描述并利用对卫星观测的拟合平差滤波从而确定卫星轨道或对卫星轨道进行改进的技术和方法，解决了天体力学和轨道理论中的奇点问题，使卫星定轨变成了用理论模型对卫星观测值进行拟合的简单过程，从而不需要积分器，不需要求解变分方程，适合于星上自行实时定轨的应用。可提高导航卫星的定轨精度从而提高导航卫星的定位精度。另外本发明技术计算简化，可以有效用于实时星载应用。
【专利附图】
【附图说明】
[0009]图1为四个角度开普勒根数的几何意义示意图。(其中equator为赤道，orbit为轨道，node为升交点，perigee为近地点，satellite为卫星，pole为极，vernalequinox为春分点经圈，Y为春分点，Ω, ω , i, f分别为升交点赤经,近地点角，轨道倾角，真近点角)。
[0010]图2为真近点角/扁近点角万之间的关系图。(其中a为轨道椭圆长半轴，b为轨道短半轴，S为卫星$力卫星S在以a为半径的圆上的投影，O为轨道椭圆的焦点，灯为椭圆中心)。
[0011]图3为轨道面内的非保守力两分量示意图。(其中O是轨道椭圆的焦点，f是真近点角，fr, fa, fh是摄动是力的三个分量，前两个分量是在轨道面内，第一个分量是径向分量，第二个分量与第一个垂直且指向卫星运动方向，第三个分量与前两个分量构成右手系)。
[0012]图4为发明技术的实施步骤的示意框图。
【具体实施方式】
[0013]现结合实施例对本发明作进一步描述。
[0014]本实施例的定轨轨道改进技术包括五部分:中间变量对时间导数的推导，奇点与否的判断，组成无奇点的卫星受摄运动方程，利用中值定理积分求解受摄运动方程，根据卫星观测进行拟合平差滤波确定轨道和对轨道进行改进，根据发明权利登记号201310545077.2所述也获得保守力位作用下的轨道解或改进，本发明与发明权利登记号201310545077.2所述的技术联在一起，组成全部的无奇点基础的知识。
[0015]中间变量对时间偏导数的求法原理与高等数学中微分方法原理相同。
[0016]分别求得非保守力摄动相对于六个中间变量的时间导数。
[0017]分别对非保守力摄动相对于开普勒轨道根数轨道面倾角和轨道椭圆扁率的时间导数作出是否奇异的判断。
[0018]分别对卫星轨道是否为圆形轨道或赤道轨道作出判断。
[0019]由上述实施例可以看出，非保守力作用下的导航卫星轨道确定改进的无奇点技术，主要利用无奇点理论，依据摄动位对轨道根数的偏导数的推导及其是否奇异的判断，决定如何组成无奇点的运动方程，从而整个问题没有奇点，可以积分获得理论描述，从而可以根据观测，通过拟合平差滤波确定卫星轨道。
[0020]由于不需要数值积分，所以不需要积分器，使得计算量大幅度减低，从而使得本发明在星上的实时应用成为可能。
[0021]利用本发明进行卫星定轨和轨道改进可以实现北斗系统的自主能力，减少对地面控制站的依赖，从而提高北斗系统在实战中的生存能力，从而提高北斗应用的范围，并具有计算简单、精度高，成本降低等优点。
[0022]可见，利用非保守力作用下的导航卫星轨道确定改进的无奇点技术可以使导航卫星系统的作用精度有很大的提高。
[0023]利用本发明可以解决非保守力作用下的卫星定轨和轨道改进，利用发明权利登记号201310545077.2所发明的保守位力作用下的导航卫星轨道确定改进的无奇点技术，可获得在保守力位作用下轨道确定或轨道改进的结果。
[0024]综合上述所获结果，即本发明与发明权利登记号201310545077.2所述发明保守力作用下的导航卫星轨道确定改进的无奇点技术的联合应用，构成了一个完整的无奇点轨道确定改进技术。
【权利要求】
1.非保守力作用下的导航卫星轨道确定改进的无奇点技术，其特征在于，包括如下步骤: 步骤1:通过非保守力函数，定义和求取六个中间轨道变量之一的时间导数，即
2.根据权利要求1所述，其特征在于:根据根据非保守力函数和定义(1-6)求取六个中间变量(开普勒根数加下标I)的时间导数(步骤1-6)。
3.根据权利要求1所述，其特征在于:根据(3，5)两式以及初始轨道根数对奇异情况作出判断(步骤7-10)。
4.根据权利要求1所述，其特征在于:根据判据获得非保守力下的无奇点运动方程(步骤 11-16)。
5.根据权利要求1所述，其特征在于:根据积分中值定理获得非保守力作用下无奇点的运动方程的理论解表达式(步骤17)。
6.根据权利要求1所述，其特征在于:根据对卫星的观测数据拟合平差滤波获得确定的非保守力下的轨道或轨道改进(步骤18)。
7.根据发明权利登记号201310545077.2所发明的保守力位作用下的导航卫星轨道确定改进的无奇点技术，其特征在于:可获得在保守力位作用下轨道确定或轨道改进的结果(步骤19)。
8.根据权利要求1所述，根据发明权利登记号201310545077.2所发明的保守力作用下的导航卫星轨道确定改进的无奇点技术，其特征在于:综合上述步骤18和步骤19所获得的结果，即本发明与发明权利登记号201310545077.2所发明的保守力作用下的导航卫星轨道确定改进的无奇点技术的联合应用，构成了一个完整的无奇点轨道确定改进技术(步骤 20)。
【文档编号】G01S19/01GK103744089SQ201310579068
【公开日】2014年4月23日 申请日期:2013年11月15日 优先权日:2013年11月15日
【发明者】许国昌, 许佳, 曲云英, 李忠华, 李鹏 申请人:福建纳威导航科技有限责任公司