一种高轨卫星厘米级定轨系统与方法技术方案

本发明专利技术公开了一种高轨卫星厘米级精密定轨系统与方法，可以实现高轨卫星厘米级精密轨道确定，满足当前高轨卫星应用对厘米级高精度轨道的需求。从观测方程层面直接将低轨卫星观测方程和地面监测站观测方程进行融合，通过建立合理正确的函数关系，联合估计北斗/GNSS轨道参数、低轨卫星轨道参数和高轨卫星轨道参数，进而通过残差检验、迭代循环的方式得到北斗/GNSS卫星、低轨卫星、高轨卫星的最优估值，最后再通过轨道积分，得到北斗/GNSS卫星、低轨卫星、高轨卫星的厘米级精密轨道。高轨卫星的厘米级精密轨道。高轨卫星的厘米级精密轨道。

【技术实现步骤摘要】
一种高轨卫星厘米级定轨系统与方法


[0001]本专利技术属于卫星导航
，特别涉及一种高轨卫星厘米级精密定轨系统与方法。

技术介绍

[0002]近三十年来我国航天发展日新月异，各类人造卫星为我国国民经济发展、社会生产生活、国防安全建设等方面发挥了重大作用，也是我国迈向世界大国强国的重要标志。随着航天技术的不断发展与进步，卫星应用在各领域，如通信导航、遥感探测、气象研究、军事侦察、深空探测等方面的作用日益突出，成为各行业不可或缺的信息工具。人造卫星根据轨道高度，区分为低轨卫星、中轨卫星和高轨卫星。卫星的轨道精度作为各类卫星业务的“空间基准”，直接影响到卫星服务的连续性、可用性与应用水平，也在一定程度上决定了潜在用户和市场的开发与拓展潜力。
[0003]一般将轨道高度低于2000km的卫星称为低轨卫星，轨道高度介于2000km和35786km的为中轨卫星，而轨道高度大于35786km的为高轨卫星。当前低轨卫星大多应用于对地观测卫星、空间站、以及以“星链”为代表的下一代通信星座；中轨卫星一般应用于导航卫星星座，如我国的北斗系统、美国的GPS系统、欧盟的GALILEO系统等；而高轨卫星，由于其独特的高轨特性、静地特性，一颗卫星就几乎可以覆盖三分之一地球，越来越广泛的应用于通信、监测、气象、天文、深空探测等各个领域，具有中低轨不可替代的作用，成为地球空间段最重要的轨道资源。
[0004]相较于高轨卫星，当前大部分中低轨卫星均配备有星载GNSS接收机，借助于星载GNSS接收机的伪距和相位观测数据，可以实现中低轨卫星实时或事后厘米级精密定轨，满足其高精度轨道精度需求。对于高轨卫星定轨问题，传统定轨所采用的地面站跟踪方式，设施复杂，成本高，精度低，而且还受限于我国国土疆域限制，无法在国外区域布设跟踪站。如果采用类似与中低轨卫星搭载星载GNSS接收机的方式，由于高轨卫星运行轨道高于GNSS卫星轨道，而GNSS导航卫星信号发射方向是指向地心，会造成高轨卫星的星载GNSS接收机接收到的是从地球另一侧导航卫星发射的边缘信号，造成高轨卫星可见GNSS卫星数量过少和信噪比过低的问题，严重影响基于星载GNSS接收机的高轨卫星轨道确定精度，制约高轨卫星在通信、导航、遥感、侦察、科研等方面的应用潜力。

技术实现思路

[0005]本专利技术解决的技术问题：针对高轨卫星对厘米级高精度轨道的需求，结合当前低轨卫星星座发展前景，充分利用高中低地多系统多层次观测数据，提出了一种适用于高轨卫星的厘米级精密定轨系统与方法，满足当前高轨卫星应用对厘米级高精度轨道的需求。
[0006]为实现上述目的，本专利技术公开了一种高轨卫星厘米级精密定轨系统与方法。
[0007]一种高轨卫星厘米级精密定轨系统，该定轨系统包括北斗/GNSS导航卫星、高轨卫星、低轨卫星星座、地面监测站、地面控制中心；
[0008]所述的北斗/GNSS导航卫星在星载原子钟的控制下生成下行导航信号，并向近地用户连续播发以提供定位测速授时服务，北斗/GNSS导航卫星包括中国北斗系统、美国GPS系统、俄罗斯GLONASS系统、欧盟GALILEO系统中的至少一个；
[0009]所述的近地用户是指低轨卫星星座和地面监测站；
[0010]所述的高轨卫星配置了高稳原子钟，用于提供高精度的时频基准，高稳原子钟秒稳～百秒稳优于1e
‑
12量级，万秒稳优于1e
‑
14量级；高轨卫星还配置了导航信号生成单元，导航信号生成单元能够在高稳原子钟提供的时频基准控制下，生成下行导航信号，并将生成的下行导航信号向近地用户连续播发；
[0011]所述的低轨卫星星座包括两个以上的低轨卫星，能够实现对高轨卫星的一重及以上覆盖，即任一时刻对高轨卫星的可见卫星数(仰角大于7
°
)不少于1颗；所有的低轨卫星均配置了高精度星载GNSS接收机，GNSS接收机能够接收北斗/GNSS导航卫星下行导航信号和高轨卫星下行导航信号；低轨卫星之间具有星间链路，通过星间链路能够将低轨卫星观测数据在低轨卫星之间进行实时数据传输；低轨卫星具有数据下传链路，能够将低轨卫星观测数据实时传输至地面控制中心；低轨卫星观测数据包括两部分内容，一部分内容为低轨卫星星载GNSS接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；另一部分内容为低轨卫星星载GNSS接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；并将由低轨卫星星载GNSS接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为低轨卫星观测数据A，将由低轨卫星星载GNSS接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为低轨卫星观测数据B，即低轨卫星观测数据A是由低轨卫星星载GNSS接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；低轨卫星观测数据B是由低轨卫星星载GNSS接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；
[0012]所述的地面监测站配备有高精度GNSS监测接收机，能够接收北斗/GNSS下行导航信号和高轨卫星播发的下行导航信号，生成地面监测站观测数据，并将生成的地面监测站观测数据实时传输至地面控制中心；地面监测站观测数据包括两部分内容，一部分为由地面监测站高精度GNSS监测接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；另一部分为由地面监测站高精度GNSS监测接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；将由地面监测站高精度GNSS监测接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为地面监测站观测数据A，将由地面监测站高精度GNSS监测接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为地面监测站观测数据B，即地面监测站观测数据A是指由地面监测站高精度GNSS监测接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；地面监测站观测数据B是指由地面监测站高精度GNSS监测接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；
[0013]所述的地面控制中心包括数据管理子系统和数据处理子系统，数据管理子系统用
于接收、存储、管理低轨卫星观测数据；数据管理子系统用于接收、存储、管理地面监测站观测数据；数据管理子系统用于接收、存储、管理高轨卫星精密轨道和精密钟差；数据处理子系统利用地面监测站观测数据、低轨卫星观测数据进行数据处理，生成高轨卫星厘米级精密轨道和精密钟差；
[0014]所述的数据处理子系统利用地面监测站观测数据、低轨卫星本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种高轨卫星厘米级定轨系统，其特征在于：该定轨系统包括北斗/GNSS导航卫星、高轨卫星、低轨卫星星座、地面监测站、地面控制中心；所述的北斗/GNSS导航卫星在星载原子钟的控制下生成下行导航信号，并向近地用户连续播发；所述的高轨卫星配置了高稳原子钟和导航信号生成单元，导航信号生成单元能够在高稳原子钟提供的时频基准控制下，生成下行导航信号；所述的低轨卫星星座包括两个以上的低轨卫星，能够实现对高轨卫星的一重及以上覆盖，低轨卫星均配置了星载GNSS接收机，GNSS接收机能够接收北斗/GNSS导航卫星下行导航信号和高轨卫星下行导航信号；低轨卫星之间具有星间链路，通过星间链路能够将低轨卫星观测数据在低轨卫星之间进行实时数据传输；低轨卫星具有数据下传链路，能够将低轨卫星观测数据实时传输至地面控制中心；所述的地面监测站配备有高精度GNSS监测接收机，能够接收北斗/GNSS下行导航信号和高轨卫星播发的下行导航信号，生成地面监测站观测数据，并将生成的地面监测站观测数据实时传输至地面控制中心；所述的地面控制中心包括数据管理子系统和数据处理子系统，数据管理子系统用于接收、存储、管理低轨卫星观测数据和地面监测站观测数据；数据处理子系统用于对地面监测站观测数据、低轨卫星观测数据进行数据处理，生成高轨卫星厘米级精密轨道。2.根据权利要求1所述的一种高轨卫星厘米级定轨系统，其特征在于：北斗/GNSS导航卫星包括中国北斗系统、美国GPS系统、俄罗斯GLONASS系统、欧盟GALILEO系统中的至少一个；所述的近地用户是指低轨卫星星座和地面监测站。3.根据权利要求1所述的一种高轨卫星厘米级定轨系统，其特征在于：所述的高轨卫星配置的高稳原子钟秒稳～百秒稳优于1e
‑
12量级，万秒稳优于1e
‑
14量级；高轨卫星将生成的下行导航信号向近地用户连续播发；所述的低轨卫星星座中任一时刻对高轨卫星的可见卫星数(仰角大于7
°
)不少于1颗。4.根据权利要求1所述的一种高轨卫星厘米级定轨系统，其特征在于：低轨卫星观测数据包括两部分内容，一部分内容为低轨卫星星载GNSS接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；另一部分内容为低轨卫星星载GNSS接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；并将由低轨卫星星载GNSS接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为低轨卫星观测数据A，将由低轨卫星星载GNSS接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为低轨卫星观测数据B。5.根据权利要求1所述的一种高轨卫星厘米级定轨系统，其特征在于：地面监测站观测数据包括两部分内容，一部分为由地面监测站高精度GNSS监测接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；另一部分为由地面监测站高精度GNSS监测接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据；将由地面监测站高精度GNSS监测接收机接收北斗/GNSS导航卫星下行导航信号，然后生成的伪距观测数据、载波
相位观测数据和多普勒观测数据称为地面监测站观测数据A，将由地面监测站高精度GNSS监测接收机接收高轨卫星下行导航信号，然后生成的伪距观测数据、载波相位观测数据和多普勒观测数据称为地面监测站观测数据B。6.根据权利要求1所述的一种高轨卫星厘米级定轨系统，其特征在于：所述的数据处理子系统利用地面监测站观测数据、低轨卫星观测数据进行数据处理，生成高轨卫星厘米级精密轨道和精密钟差的方法有两种，第一种为高中低地一体化精密定轨与时间同步方法，另一种为高轨卫星单系统精密定轨方法，其中高中低地一体化精密定轨与时间同步方法的的步骤包括：第一步，建立地面监测站对北斗/GNSS导航卫星和高轨卫星观测模型的误差方程；地面监测站对北斗/GNSS导航卫星和高轨卫星观测模型的误差方程为：v
sta
＝F(X
GNSS
，X
HEO
，X
sta_y
，X
sta_n
，t
sta
)
‑
obs
sta
式中，v
sta
表示地面监测站相对于北斗/GNSS导航卫星观测方程的残差和相对于高轨卫星观测方程的残差；F(X
GNSS
，X
HEo
，X
sta_y
，X
sta_n
，t
sta
)表示X
GNSS
，X
HEO
，X
sta_y
，X
sta_n
，t
sta
与obs
sta
的函数关系；obs
sta
表示地面监测站对北斗/GNSS导航卫星和高轨卫星的观测数据，包括伪距观测数据、载波相位观测数据和多普勒观测数据；X
GNSS
表示北斗/GNSS卫星的动力学参数，包括北斗/GNSS卫星初始状态向量、大气阻力参数、光压参数、经验力参数；X
HEO
表示高轨卫星的动力学参数，包括高轨卫星初始状态向量、大气阻力参数、光压参数、经验力参数；X
sta_y
表示与地面监测站坐标相关参数，包括测站坐标、地球固体潮改正、UT1变化率；X
sta_n
表示与地面监测站坐标无关的参数，包括地面监测站接收机钟差、北斗/GNSS卫星钟差、高轨卫星钟差、模糊度参数、对流层参数和电离层参数；t
sta
表示地面监测站观测历元；第二步，建立低轨卫星对北斗/GNSS导航卫星和高轨卫星观测模型的误差方程；低轨卫星对北斗/GNSS导航卫星和高轨卫星观测模型的误差方程表示为：v
leo
＝G(X
GNSS
，X
HEO
，X
LEO
，X
leo_y
，X
leo_n
，t
leo
)
‑
obs
leo
式中，v
leo
表示低轨卫星相对于北斗/GNSS导航卫星观测方程的残差和相对于高轨卫星观测方程的残差；G(X
GNSS
，X
HEO
，X
LEO
，X
leo_y
，X
leo_n
，t
leo
)表示X
GNSS
，X
HEO
，X
LEO
，X
leo_y
，X
leo_n
，t
leo
与obs
leo
的函数关系；obs
leo
表示低轨卫星对北斗/GNSS导航卫星和高轨卫星的观测数据，包括伪距观测数据、载波相位观测数据和多普勒观测数据，X
GNSS
表示北斗/GNSS卫星动力学参数，包括北斗/GNSS卫星初始状态向量、大气阻力参数、光压参数、经验力参数；X
HEO
表示高轨卫星的动力学参数，包括高轨卫星初始状态向量、大气阻力参数、光压参数、经验力参数；X
LEO
表示低轨卫星的动力学参数，包括低轨卫星初始状态向量、大气阻力参数、光压参数、经验力参数；X
leo_y
表示低轨卫星坐标参数，X
leo_n
表示与低轨卫星坐标无关的参数，包括低轨卫星星载接收机钟差、北斗/GNSS卫星钟差、高轨卫星钟差、模糊度参数和电离层参数；t
leo
表示低轨卫星观测历元；第三步，将第一步建立的误差方程在待估参数近似值第三步，将第一步建立的误差方程在待估参数近似值处进行泰勒公式展开，得到：
第四步，将第二步建立的误差方程在待估参数近似值第四步，将第二步建立的误差方程在待估参数近似值处进行泰勒公式展开，得到：式中，dX
GNSS
、dX
HEO
、dX
sta_y
、dX
sta_n
、dX
LEO
、dX
leo_y
、dX
leo_n
表示待估参数的改正数；第五步，采用最小二乘估计，得到待估参数改正数的最优估值第五步，采用最小二乘估计，得到待估参数改正数的最优估值第六步，将待估参数初始值第六步，将待估参数初始值和待估参数改正数dX
GNSS
、dX
HEO
、dX
sta_y
、dX
sta_n
、dX
LEO
、dX
leo_y
、dX
leo_n
相加，得到待估参数的最优估值。7.根据权利要求6所述的一种高轨卫星厘米级定轨系统，其特征在于：其中高轨卫星单系统精密定轨的方法步骤包括：第一步，从国际IGS数据中心或全球科研机构，获取北斗/GNSS精密轨道、精密钟差和地面监测站精确坐标文件；从低轨卫星星座运行控制机构获取低轨卫星精密轨道文件；第二步，建立地面监测站对高轨卫星观测误差方程，地面监测站对高轨卫星观测误差方程表示为：v
staGEO
＝F
GEO
(X
HEO
，t
sta
)
‑
obs
staGEO
式中，v
staGEO
表示地面监测站对高轨卫星观测方程的残差；F
GEO
(X
HEO
，t
sta
)表示X
HEO
与obs
staGEO
的函数关系；obs
staGEO
表示地面监测站对高轨卫星的观测数据，包括伪距观测数据、载波相位观测数据和多普勒观测数据；X
HEO
表示高轨卫星的动力学参数，包括高轨卫星初始状态向量、大气阻力参数、光压参数、经验力参数；t
sta
表示地面监测站观测历元；第三步，建立低轨卫星对高轨卫星观测误差方程，低轨卫星对高轨卫星观测误差方程表示为：v
leoGEO
＝G
GEO
(X
HEO

【专利技术属性】
技术研发人员：蒙艳松，周泉，边朗，王瑛，张蓬，严涛，田野，
申请(专利权)人：西安空间无线电技术研究所，
类型：发明
国别省市：