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