CN109764879B - 一种卫星定轨方法、装置及电子设备 - Google Patents

Abstract

本申请提供一种卫星定轨方法、装置及电子设备，所述方法包括：获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

Description

一种卫星定轨方法、装置及电子设备

本申请要求在2018年12月28日提交中华人民共和国知识产权局、申请号为201811627358.1，发明名称为“一种卫星精密定轨方法及装置”的中国专利申请的优先权，其全部内容通过引用结合在本申请中。

技术领域

本申请涉及卫星导航技术领域，尤其涉及一种卫星定轨方法、装置及电子设备。

背景技术

卫星导航定位系统(GNSS)能够提供全球范围内实时的定位服务，并在世界各国的许多行业中得到了广泛的应用。目前标准的GNSS定位技术精度大约为5-10米。对于精度要求更高的应用，通常需要使用精密定位的方法。用户通过接收导航卫星发送的导航信号，并以导航卫星作为动态己知点，利用伪距等观测信息实时地测定运动载体的在航位置和速度，进而完成导航。

实时获取精确的卫星轨道，是实现精密卫星导航定位的关键技术。目前，为提高导航精度，卫星导航系统还可以包括卫星导航增强系统主要分为星基增强系统(SBAS)和地基增强系统(GBAS)两大类。星基增强系统如美国的广域增强系统(WAAS)、俄罗斯的差分校正和监测系统(SDCM)等，地基增强系统如美国的局域增强系统(LAAS)等。使用增强系统之后，卫星静态定位精度可以达到厘米级，动态精度可以达到米级(车道级)。

目前，卫星轨道的确定方法主要包括3种：地面监测站定轨、星上自主定轨和事后精密定轨。目前，GNSS系统，包括GPS，GLONASS，BDS等都使用少数的地面跟踪站实现轨道全弧段观测，然后计算和预报卫星轨道，上注到导航卫星。

地基增强通过提供差分修正信号，可达到提高卫星导航精度的目的；优化后的定位精度可以从毫米级至亚米级不等。基于连续运行的永久参考站，解算出的修正数包括区域信号(类似CORS信号)和广域差分信号(类似SBAS)，播发方式包括移动网络/UHF电台/同步卫星等。基于地基增强系统播发非差综合改正信息的方式改正流动站处相应误差，从而达到模糊度参数与位置参数的快速分离，可在几个历元内固定模糊度参数，实现实时的差分定位，如RTK等，具有高精度、高实时性等优势，目前地基增强系统中的精密定位的方法有差分定位方法和精密单点定位方法。其中，差分定位方法包括局域差分方法和广域差分方法。广域差分方法和精密单点定位都需要通过地面监测网的数据计算出各自精密的信号偏差，精密轨道等，然后在用户端使用这些精密轨道和钟差精细改正来提高定位精度。局域差分方法主要是将参考站的观测数据和坐标直接播放给用户，在用户端通过观测值差分的方式消除各自误差的影响，实现高精度相对定位。但是，广域差分方法缺乏精密电离层模型支持，需要20-30分钟收敛才能获得厘米级的定位结果；局域差分方法中，需要用户接收机与参考站之间的距离在一定范围内，对地基监测网络布站密度要求较高。

SBAS(Satellite-Based Augmentation System)星基增强系统，通过地球静止轨道(GEO)卫星搭载卫星导航增强信号转发器，可以向用户播发星历误差、卫星钟差、电离层延迟等多种修正信息，实现对于原有卫星导航系统定位精度的改进，从而成为各航天大国竞相发展的手段。

星上自主定轨是依靠星上的GNSS接收机或惯性测量单元进行轨道确定，其中，基于GNSS观测值的星上自主定轨方法能提供实时自主连续的卫星轨道，但是，由于受GNSS导航卫星广播星历的轨道误差和钟差的影响，基于GNSS观测值的星上自主定轨通常只能获得数米量级的定轨精度。为获得高精度的定轨，星基增强卫星主要分布在GEO轨道，分为面向民航等用户的广域差分完好性增强与面向测绘等高精度用户的广域精密定位增强，具备广域覆盖的特点，但是目前的广域精密定位技术收敛时间长，难以形成面实时性要求较高的PNT应用，需要解决其收敛时间较慢的问题。

发明内容

本发明实施例提供一种卫星定轨方法、装置及电子设备，用于有效提高精密定轨的收敛时间。

本发明实施例提供一种卫星定轨方法，该方法包括：

获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；

根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

一种可能的实现方式，所述第二观测数据包括载波相位观测数据；所述第二观测数据通过以下方式获得，包括：

低轨卫星接收机若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；

将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；

其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；

所述确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程，包括：

根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

一种可能的实现方式，所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

本发明实施例提供一种卫星定轨装置，该方法包括：

收发单元，用于获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

处理单元，用于确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

一种可能的实现方式，所述处理单元，用于若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；

所述处理单元，具体用于根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

一种可能的实现方式，所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

本发明实施例提供一种电子设备，包括：

至少一个处理器；

以及与所述至少一个处理器通信连接的存储器；

所述存储器存储有可被所述至少一个收发机执行的指令，所述指令被所述至少一个收发机执行：

获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行：

确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

一种可能的实现方式，所述处理器，具体用于：

若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；

所述处理器，具体用于根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

一种可能的实现方式，所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

本发明实施例提供一种计算机存储介质，存储有计算机可执行指令，所述计算机可执行指令用于执行如本发明实施例提供任一项的方法。

本发明实施例提供一种计算机程序产品，所述计算机程序产品包括存储在计算机可读存储介质上的计算程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，使所述计算机执行上述实施例中任一种方法。

本发明实施例中，地面接收机和低轨卫星的接收机在当前定位时刻确定对导航卫星的第一观测值和第二观测值，分别对导航卫星的精密轨道和低轨卫星的精密轨道进行解算，由于低轨卫星相对地面运动速度较快(角速度大)，观测数据几何变化快，可加快载波相位整周模糊度估计收敛速度，大幅提高定位精度并缩短高精度定位收敛时间。因此，通过结合地面的观测数据和低轨卫星的观测数据，实现了对导航卫星的精密定轨的同时，减少了收敛时间。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简要介绍。

图1示出了本发明实施例提供的一种卫星定轨方法的系统架构示意图；

图2示出了本发明实施例提供的一种卫星定轨方法的流程示意图；

图3示出了本发明实施例提供的一种获取载波相位观测数据的方法的流程示意图；

图4示出了本发明实施例提供的一种卫星定轨方法的流程示意图；

图5示出了本发明实施例提供的一种卫星定轨方法的流程示意图；

图6示出了本发明实施例提供的一种卫星定轨装置的结构示意图；

图7示出了本发明实施例提供的一种电子设备的结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。

下面介绍本发明实施例的技术背景。

2015年，美国波音公司获得了美国军方约1.6亿美元合同，进一步开展基于铱星的GPS导航增强技术研究和应用实验。为此组建了新的研究团队，主要参与者更变为Satelles公司。新的低轨增强技术基于铱星系统新设置的GPS导航增强专用STL(Satelles Time andLocation)服务业务，在iGPS的基础上，进一步强化增强功能；国内低轨导航增强已经成为研究热点，多个星座计划或高校已经围绕低轨导航增强做了大量的体制论证工作，单高实时高精度服务尚处于空白；国外下一代铱星已经启动了低轨导航增强系统的建设，领先一步，有望实现低轨导航增强，并形成军民商应用；国内，低轨导航增强已经成为研究热点，多个星座计划或高校已经围绕低轨导航增强做了大量的体制论证工作，但都刚刚起步，高实时高精度服务尚处于空白。

如图1所示，本发明实施例提供的的系统架构示意图，该卫星定轨的系统架构中包括，导航卫星，低轨卫星，地面接收机和低轨卫星接收机，及地面的信号处理系统。地面接收机可以接收导航卫星发射的导航信号。低轨卫星相对地面运动速度较快，角速度大，观测数据几何变化快，可加快载波相位整周模糊度估计收敛速度，大幅提高定位精度并缩短高精度定位收敛时间，因此，相对传统的增强服务，由于LEO卫星对地面运动较快，星座的几何构型变化较快，通过实施精密定轨后的LEO卫星，向客户播发LEO星座卫星导航和增强信号有利于载波相位模糊度参数的快速收敛和固定，为解决精密单点定位PPP技术收敛时间过长的问题提供契机。本发明实施例采用LEO卫星上的星载接收机获取导航卫星的观测数据，以减少导航卫星定轨所需的收敛时间。

相对于传统的卫星导航设备，本发明实施例中的接收机由于安装在低轨卫星上，速度和动态都远大于车载、机载卫星导航设备，需要在卫星信号的捕获与跟踪算法上作相应的调整保证捕获灵敏度以及载波跟踪环的跟踪精度。

为了实现高实时高精度服务，作为核心载荷的多频星载接收机的硬件主要瓶颈在于制造卫星载荷复杂度高，软件算法则需要验证高精度轨道钟差改正数解算、卫星相位延迟解算、区域高精度电离层建模、区域观测数据与LEO观测数据联合定轨、实时区域站观测数据实时精密钟差解算模型、LEO精密钟差确定等关键技术。

下面结合附图介绍本发明实施例提供的技术方案。

结合图1，本发明实施例提供一种卫星定轨方法，如图2所示，该方法包括：

步骤201：获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

步骤202：确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；

步骤203：根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

需要说明的是，根据不同应用场合，星载接收机的功能不同，一般来说多数星载接收机用于卫星定轨、授时和定姿等。定轨精度一般在1-2米，根据频点和星座规划不同，一般分单星座单频，多星座多频接收机，本发明实施例中，为了实现模糊度搜索效率提高，可以更快获取模糊度整数解，可以采用双星座双频星载接收机。具体的，该接收机可以同时接收GPS的L1、L2和/或北斗卫星B1、B2频点的信号，以提供双频的高精度观测数据。

本发明实施例中，地面接收机和低轨卫星的接收机在当前定位时刻确定对导航卫星的第一观测值和第二观测值，分别对导航卫星的精密轨道和低轨卫星的精密轨道进行解算，由于低轨卫星相对地面运动速度较快(角速度大)，观测数据几何变化快，可加快载波相位整周模糊度估计收敛速度，大幅提高定位精度并缩短高精度定位收敛时间。因此，通过结合地面的观测数据和低轨卫星的观测数据，实现了对导航卫星的精密定轨的同时，减少了收敛时间。

本发明的技术解决问题首先是观测数据生成问题是LEO卫星观测弧段短、运行速度快、大气阻力影响大，导致地面站所接收到的LEO卫星观测数据中周跳较多、粗差大；本发明实施例中的星载接收机环路使用考虑到角速度比较大，速度很快，载波环路机制需要保证跟踪的准确性，常见的接收机的载波相位观测数据一般只用来平滑伪距和多普勒的计算，计算时只需要载波相位前后时刻的增量数据减去中频累加部分。

而如果用于非差的PPP精密定位，则需要连续的载波相位累加值输出，需要考虑到钟差调整带来的影响。

一种可能的实现方式，所述第二观测数据包括载波相位观测数据；所述第二观测数据通过以下方式获得，包括：

低轨卫星接收机若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；

将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

本发明实施例中，载波环路控制主要包括锁频环FLL、二阶锁相环PLL2以及2阶锁频辅助3阶锁相环的机制下切换。

扩频信号的捕获方法主要有串行捕获、并行捕获和基于FFT的快速捕获。串行捕获方法是一个频域和时域的顺序二维搜索过程，其捕获时间较长；并行捕获方法采用多个捕获通道，各个通道分别并行完成接收信号与不同码相位和不同多普勒频率的本地再生信号的相关计算，相对于串行捕获，捕获速度大。

直扩信号的捕获过程就是通过在时间域和频域进行二维搜索来检测y的模平方是否超过由噪声统计决定的门限值，基于FFT的快速捕获方法可以在给定ω下，在同一时间段内搜索伪码所有相位，因此捕获速度很高。

在高动态环境下，由于多普勒频移的不确定性，直接捕获载波相位有较大的难度；另外，为了提高动态跟踪能力，势必增加环路带宽，而增加环路带宽就会引入较大的跟踪误差。在初始捕获时，采用FFT可以对信号进行快速捕获，为了解决高动态捕获能力和提高跟踪精度的矛盾，在需要提高精度时，载波环还可以采用FLL+PLL混合载波跟踪算法。FLL环直接跟踪载波频率，通过载波鉴频器输出多普勒频率估计误差，具有较好的动态性能，但跟踪精度比PLL环的低。

简单的PLL由频率基准、相位检波器、电荷泵、环路滤波器和压控振荡器(VCO)组成。基于PLL技术的频率合成器将增加两个分频器：一个用于降低基准频率，另一个则用于对VCO进行分频。PLL是作为闭环控制系统工作，用于比较基准信号与VCO的相位。增设基准和反馈分频器的频率合成器负责比较两个由分频器的设定值调节相位。该相位比较在相位检波器中完成，该相位-频率检波器生成一个误差电压，此误差电压在±2π的相位误差范围内近似为线性，并在误差大于±2π的情况下保持恒定。相位-频率比较器所采用的这种双模式操作可生成针对大频率误差(比如，当PLL在上电期间起动时)的较快的PLL锁定时间，并避免被锁定于谐波之上。

鉴相器有两个输入，分别是输入信号和压控振荡器的输出信号，在二者相位差和频率差不是很大的情况下，鉴相器的输出与两输入信号之差成正比，鉴相器的输出为模拟信号，其通过低通滤波器虑除高频杂波，后进入压控振荡器，压控振荡器的输出频率随其输入电压的改变而改变。从原理图上看，PLL实际上是一负反馈系统，只要输入信号在正常范围内，输出信号在“一定时间内”都能跟上。输入信号发生变化后，输出信号跟踪输入信号的过程称之为捕获；输出信号跟踪完毕时称之为锁定；输入信号变化过快导致输出信号无法跟踪时称为失锁。通过PLL可以方便实现N倍频。

为进一步提高精度，可以采用二阶锁相环PLL2以及2阶锁频辅助3阶锁相环等方式。

具体的，PLL2和FLL2+PLL3中的时域和频域的公式如下:

2阶PLL环路滤波器时域核心公式如下：

其中x_n表示输入信号，y_n表示输出信号，u_n表示输出信号，a₂为环路参数，T表示调整周期，ω₀表示环路特征值；

变换到Z域可以得到

U表示鉴相输入，z表示Z变换；

则有

利用Z变换变换到时域，可以得到输出信号

2阶锁频环辅助3阶锁相环环路滤波器核心公式如下：

其中，v_n表示一级滤波器锁相输入信号，ω_f表示鉴相特征值，ω_p表示鉴频特征值，f_n表示频率；

转换到Z域可以得到

本发明实施例中，载波环路的调整和载波相位观测数据的生成机制和普通接收机不同；在获取载波相位观测数据时，根据失锁捕获以及跟踪精度，确定跟踪环路，以提高载波相位数据的精度。

在具体实施过程中，观测数据的质量和环路状态、环路精度评估机制和bit反转与否都相关，因此，在输出载波相位观测数据之前，还可以包括环路状态、环路精度评估和bit反转等参数的判断，以提高观测数据精度。

如图3所示，本发明实施例提供一种获取载波相位观测数据的方法，包括：

步骤301：接收导航卫星的直发信号，并将接收到的直发信号作为接收机中载波环路的输入信号；

步骤302：判断当前载波环路的锁定指示是否高于跟踪门限；若是，则执行复位，对所述直发信号重新跟踪；若否，则执行步骤303；

步骤303：根据当前环路滤波的需要，确定滤波器类型；若确定当前跟踪状态为初始状态或者当前跟踪精度较低，则执行步骤304；若确定当前跟踪状态为稳定跟踪或者当前跟踪精度较高，可以执行步骤305或步骤306，可以根据跟踪精度确定选择的环路；

步骤304：将输入信号依次输入至鉴频器，和FLL环路滤波器，获得输出信号；

步骤305：将输入信号依次输入至鉴相器，和PLL2环路滤波器，获得输出信号；

步骤306：将输入信号依次输入至鉴相器，和FLL2+PLL3环路滤波器，获得输出信号；

步骤307：将环路滤波器的输出信号输入至载波环路的数字控制振荡器NCO，输出载波相位累加值；

步骤308：判断载波环路状态，若确定载波环路状态失锁，则执行步骤302；若确定载波环路状态为稳定跟踪，则执行步骤309；

步骤309：判断Bt帧同步是否有效；若是，则执行步骤310；若否，则执行步骤313；

步骤310：判断环路跟踪精度是否满足要求；若是，则执行步骤311；若否，则执行步骤313；

步骤311：判断载波Bit是否需要反转；若否，则执行步骤312；若是则执行步骤313；

步骤312：确定当前时刻的载波相位累加值，并将载波相位累加值作为输出的载波相位观测数据；

步骤313：确定载波相位观测数据无效。

通过上述方法，可以实现一个完整的环路控制和切换机制，以获得满足精度需要的载波相位观测数据。

为了实现非差的精密定位，本发明实施例中，由于卫星和接收机未校准的初始相位延迟，载波相位中所包含的模糊度信息要完成的保留下来，所以整数部分和小数部分的处理必须分开，以保留载波相位观测值任意一部分，提高定轨的精度。

如图4所示，载波相位观测数据计算流程中所有的载波相位计算都是整数部分和小数部分分开，包括：

步骤401、将上一时刻的载波相位累加值、当前时刻的载波相位累加值、中频累加值和本地时钟调整量，输入至载波相位计算单元，确定出载波相位观测数据的整数部分和小数部分；

具体的，确定出载波相位观测数据的整数部分和小数部分为当前时刻载波相位累加值的整数部分和小数部分减去上一时刻的载波相位累加值整数部分和小数部分，再减去数字中频频率*测量时间间隔。

步骤402、对载波相位观测数据进行半周补偿，输出最终的载波相位观测数据。

需要说明的是，在载波相位观测数据计算流程中，同时需要把上一时刻载波相位观测值未调整部分减去，保证载波相位累加值的变化和伪距观测值保持一致。

本发明实施例提供的一种卫星定轨方法，具体处理过程可以包括：

步骤一、地面接收机和LEO星载接收机接收导航卫星播发的导航直发信号，对直发信号进行捕获、跟踪；

步骤二、在每一历元，地面接收机对导航直发信号进行测量，产生伪距、载波相位的第一观测数据；LEO星载接收机对导航直发信号进行测量，产生伪距、载波相位的第二观测数据；

步骤三、利用第一观测数据和第二观测数据建立第一观测方程和第二观测方程；

其中，地面接收机的第一观测方程和星载北斗接收机的第二观测方程可以分别表示:

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi

式中,Y_GROUNDi和Y_LEOi分别表示地面和星载北斗接收机在t_i时刻的观测值；F_GB和F_LB分别表示地面和星载北斗接收机在t_i时刻的观测函数；X_oi和X_Oi分别表示除了X_BDi和X_LEOi之外的其它待估参数,如钟差、模糊度、大气参数等；ξ_GROUNDi和ξ_LEOi是相应观测值的观测误差。

其中，若确定存在多种卫星导航系统，则以一种卫星导航系统为基准，对其他卫星导航系统和低轨卫星观测数据进行归一化处理，得到统一时间基准观测方程；

步骤四、对第一观测数据和第二观测数据，进行周跳粗差的预处理。

步骤五、利用第一观测数据，第二观测数据，卫星精密钟差，以及精密轨道模型，进行观测方程的解算，得到导航卫星和leo卫星的精密定位和接收机钟差；

步骤六、通过通信链路接收地基监测网播发的非差综合改正信息，或者，leo卫星播发非差综合改正信息；

步骤七、根据接收非差综合改正信息计算用户概略位置相对每颗导航卫星和低轨卫星的误差改正参数；

步骤八、采用精密单点定位模式进行定位处理，授时和测速结果及载波相位模糊度参数等。

通过本申请提供的一种卫星定轨方法，可在全球获得近实时的精密定位、测速和授时结果。该星载接收机提供的双频观测数据结合GNSS卫星精密轨道和钟差及硬件延迟建模信息以及地面站观测数据可实现快速实时的LEO/GNSS联合PPP定轨。

本发明的技术解决问题涉及到观测数据的处理问题，特别是关于双频载波相位处理的问题，关于PPP的研究主要基于双频观测数据，已经实现了静态mm-cm、动态cm-dm级的定位精度。

与单频精密单点定位不同的是，双频精密单点定位对双频伪距和相位观测值进行组合来消除电离层一阶项的影响。双频精密单点定位通常组成Melbourne-Wübbena和Geometry-Free组合，采用Turbo-Edit方法进行周跳探测。双频精密单点定位的观测模型及双频载波与伪距的无电离层组合观测方程为：下面通过一个具体实施例对本发明技术方案进行详细说明。

如图5所示，基于非差的PPP定轨的流程示意图，可以包括：

步骤501、对第一观测数据和第二观测数据周跳粗差探测数据预处理。

步骤502、利用第一观测数据和第二观测数据建立第一观测方程和第二观测方程；

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；

其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

其中，地面接收机接收导航卫星的导航信号产生的第一观测数据包含多星座多频点伪距、载波相位和多普勒观测数据；LEO星载接收机接收导航卫星的导航信号产生的第二观测数据。

第一观测方程可以包括载波相位观测方程和伪距方程；第二观测方程也可以包括载波相位观测方程和伪距方程。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；所述确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程，包括：

根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

在具体实施过程中，可以根据双频的第一观测数据和第二观测数据，构造无电离层组合的第一观测数据和无电离层组合的第二观测数据，消除一阶电离层延迟影响，减少未知参数，具体的，载波相位观测方程和伪距方程的模型为：

其中：λ_lc为无电离层组合相位波长，

为无电离层载波相位观测值(以距离表示)；p表示卫星，k表示测站；

为无电离层伪距观测值；

为站星几何距离；

为模糊度；dt_k为接收机钟差，dt^p为卫星钟差；

对流层延迟；

为多路径效应；

相位和伪距观测的噪声。

在具体实施过程中，第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

其中，对流层延迟通常可以分为干分量和湿分量两部分。干分量可以通过模型进行改正，湿分量作为待估参数进行估计。为了减少待估参数的数量，可以使用映射函数将斜延迟投影到天顶方向，只估计一个天顶湿延迟。

利用相对论效应、地球自转、天线相位中心等模型对观测方程进行改正，消去部分参数，同时忽略残留的卫星轨道和钟差误差。

根据导航卫星和低轨卫星的运动方程和变分方程积分,就可以分别得到它们的初始参考轨道和状态转移矩阵；

根据卫星精密定轨问题的一般性描述可知，对导航卫星和低轨卫星的运动方程和变分方程积分,就可以分别得到它们的初始参考轨道和状态转移矩阵

φ_LEO(t_i,t₀)，其中状态转移矩阵应满足方程：

式中,

和

分别表示北斗卫星和低轨卫星在时刻t_i的状态改正数向量,

和

分别表示它们在初始时刻t₀(参考历元)的状态改正数向量。上式用于将在其它时刻的状态改正数映射到初始时刻，以便参与最后的最优参数估计。

根据初始参考轨道，对观测方程进行线性化。

具体的，在初始参考轨道的接收机近似位置进行泰勒展开，舍弃二阶项，得到：

其中，(x,y,z)为低轨卫星或导航卫星精密轨道坐标，(xr,0，yr,0，zr,0)为接收机近似位置。

于是观测方程可以简化写为：

V＝AΔX+L

其中V为观测残差，A为系数矩阵，ΔX为包括接收机坐标改正、接收机钟差、对流层天顶湿延迟、载波相位模糊度在内的未知向量，L为计算向量。

步骤503、结合参数最优估计方法对线性化的观测方程进行解算，获得导航卫星和低轨卫星的精密定轨位置。

具体的，在进行参数估计和模糊度固定处理的具体实施过程中，可以采用最小二乘法或Kalman滤波进行综合PPP处理。

星地联合定轨中存在观测数据量大、估计参数多的问题，所以在数据处理中需要使用实时预消除参数的方法消除法方程中失效的参数(包括历元参数、时段参数)，使之有效地减小法方程的大小，加快法方程处理的时间。而且，如果在进行消参数的同时也把消除的参数与其它参数的联系方程保存起来，则可以在求解法方程后通过回代的方式恢复这些参数。

对于频繁的历元参数的消除，本项目将采取极为有效的索引策略来加快法方程的更新，极大地缩短了消参数的时间消耗；对于模糊度参数，本专利将定义好模糊度的有效时间段，在该模糊度参数刚消失时即将其从法方程中消去。从总体上优化软件的数据结构，减少函数间大型数组的传递和继承的时间，提高软件的运行效率。

下面将通过数学描述来说明参数预消除与恢复的基本流程。

设有误差方程：

将参数向量

分解成

则其法方程和相应的二次型为：

现假设

为当前历元中“失效”的参数向量，为减少法方程的维数，需对

实行预消除。

对公式中的第二式两端同时乘以

得：

再两式相加，得：

令

则消除参数向量

后的法方程为：

法方程的解为：

则可以恢复预消除参数向量

的解：

结合有：

由此可见，在消除参数向量

后，相应的二次型

需在原有的二次型上加改正

在Kalman滤波中，需要提供合适的观测值随机模型以及状态向量动态模型。随机模型描述的是观测值的统计特性，通常用观测值的方差协方差阵表示。从观测方程可知，消电离层组合观测值是原始观测值的线性组合，假设不同频率上的观测值不相关，消电离层组合观测值的初始方差可以通过误差传播定律计算得到。具体的方差可以定义为初始方差和卫星高度角的函数。假设不同卫星、不同系统的观测值不相关，以及不同类型的观测值，即伪距和相位观测值不相关，就可以得到观测值的方差协方差阵。

针对状态向量的动态模型，静态接收机坐标可以表示为常数，动态接收机坐标和接收机钟差可以表示为随机游走或者一阶高斯马尔科夫过程，对流层天顶湿延迟可以表示为随机游走过程，载波相位模糊度参数可以表示为常数，于是得到状态方程。

X_k＝Φ(t_k,t_k-1)X_k-1+w_k-1

式中，X为待估计的接收机坐标改正、接收机钟差等参数，Φ为状态转移矩阵，w_k-1为状态转移噪声。综合观测方程和状态方程，可应用标准Kalman滤波过程进行参数估计。这里由于未进行卫星相位小数偏差改正，所以仅获得载波相位模糊度浮点解结果。若进一步利用低轨卫星增强信息中包含的卫星相位小数偏差改正进行观测方程改正，则可恢复模糊度的整数特性，实现模糊度固定，得到载波相位模糊度固定解结果，进一步缩短初始化时间，提高定位、测速和授时精度。

由于增加低轨卫星的观测数据，低轨卫星的快速移动特性极大提升了解算效率，从而使PPP收敛时间大幅降低。

步骤504、利用低轨卫星播发的导航卫星增强信息和模型进行误差改正。

基于相同的发明构思，本发明实施例提供一种卫星定轨装置，如图6所示，该方法包括：

收发单元601，用于获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

处理单元602，用于确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

一种可能的实现方式，所述处理单元602，用于若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；

所述处理单元602，具体用于根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

一种可能的实现方式，所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

基于相同的发明构思，本申请提供一种电子设备，包括至少一个处理器；以及，与所述至少一个处理器通信连接的存储器；所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行，以使所述至少一个处理器能够执行上述实施例中的卫星定轨方法。

以一个处理器为例，图7为本申请提供的一种电子设备的结构示意图。如图7所示，该电子设备包括处理器701、存储器702和收发机703；其中，处理器701、存储器702和收发机703通过总线704相互连接。

其中，存储器702用于存储程序。具体地，程序可以包括程序代码，程序代码包括计算机操作指令。存储器702可以为易失性存储器(volatile memory)，例如随机存取存储器(random-access memory，简称RAM)；也可以为非易失性存储器(non-volatile memory)，例如快闪存储器(flash memory)，硬盘(hard disk drive，简称HDD)或固态硬盘(solid-state drive，简称SSD)；还可以为上述任一种或任多种易失性存储器和非易失性存储器的组合。

存储器702存储了如下的元素，可执行模块或者数据结构，或者它们的子集，或者它们的扩展集：

操作指令：包括各种操作指令，用于实现各种操作。

操作系统：包括各种系统程序，用于实现各种基础业务以及处理基于硬件的任务。

总线704可以是外设部件互连标准(peripheral component interconnect，简称PCI)总线或扩展工业标准结构(extended industry standard architecture，简称EISA)总线等。总线可以分为地址总线、数据总线、控制总线等。为便于表示，图7中仅用一条粗线表示，但并不表示仅有一根总线或一种类型的总线。

收发机703可以为通过通信接口进行通信，通信接口可以有线通信接入口，无线通信接口或其组合，其中，有线通信接口例如可以为以太网接口。以太网接口可以是光接口，电接口或其组合。无线通信接口可以为WLAN接口。

处理器701可以是中央处理器(central processing unit，简称CPU)，网络处理器(network processor，简称NP)或者CPU和NP的组合。还可以是硬件芯片。上述硬件芯片可以是专用集成电路(application-specific integrated circuit，简称ASIC)，可编程逻辑器件(programmable logic device，简称PLD)或其组合。上述PLD可以是复杂可编程逻辑器件(complex programmable logic device，简称CPLD)，现场可编程逻辑门阵列(field-programmable gate array，简称FPGA)，通用阵列逻辑(generic array logic，简称GAL)或其任意组合。一种可能的设计中，存储器702也可以和处理器701集成在一起。

所述存储器702，用于存储一个或多个可执行程序，可以存储所述处理器701在执行操作时所使用的数据。

收发机703，用于获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；

所述处理器701，用于：确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述第一观测方程和所述第二观测方程，从而确定所述导航卫星的轨道位置。

一种可能的实现方式，所述处理器701，用于：

若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

一种可能的实现方式，所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；

其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

一种可能的实现方式，所述第一观测数据和所述第二观测数据均为双频观测数据；

所述处理器701，具体用于根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程。

一种可能的实现方式，所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星接收机与导航卫星间的对流层延迟；

为地面接收机伪距方程的多路径效应；

为低轨卫星接收机伪距方程的多路径效应；

为地面接收机载波相位观测方程的多路径效应、

为低轨卫星接收机载波相位观测方程的多路径效应。

上述产品可执行本申请实施例所提供的方法，具备执行方法相应的功能模块和有益效果。未在本实施例中详尽描述的技术细节，可参见本申请实施例所提供的方法。

本申请实施例提供了一种计算机程序产品，其中，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，其中，当所述程序指令被计算机执行时，使所述计算机执行本申请上述方法实施例中任一项数据库同步延迟的确定方法。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例做出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。

显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。

Claims (9)

1.一种卫星定轨方法，其特征在于，该方法包括：

获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；所述第一观测数据和所述第二观测数据均为双频观测数据；所述第一观测数据和所述第二观测数据均包括载波相位观测数据；所述载波相位观测数据的整数部分和小数部分分开处理，以完整地保留载波相位中所包含的模糊度信息，包括：将上一时刻的载波相位累加值、当前时刻的载波相位累加值、中频累加值和本地时钟调整量，输入至载波相位计算单元，确定出载波相位观测数据的整数部分和小数部分；对载波相位观测数据进行半周补偿，输出最终的载波相位观测数据；

确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程；

根据所述第一观测数据和所述第二观测数据，解算所述消除后的第一观测方程和所述消除后的第二观测方程，从而确定所述导航卫星的轨道位置。

2.如权利要求1所述的卫星定轨方法，其特征在于，所述第二观测数据包括载波相位观测数据；所述第二观测数据通过以下方式获得，包括：

低轨卫星接收机若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；

将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

3.如权利要求1所述的卫星定轨方法，其特征在于，

所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；

其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

4.如权利要求3所述的卫星定轨方法，其特征在于，

所述第一观测方程中的观测函数包括伪距观测函数和载波相位观测函数；所述第二观测方程中的观测函数包括伪距观测函数和载波相位观测函数；

所述第一观测方程的伪距观测函数为：

所述第二观测方程的伪距观测函数为：

所述第一观测方程的载波相位观测函数为：

所述第二观测方程的载波相位观测函数为：

其中：λ_lc为无电离层组合相位波长，p表示导航卫星；

为导航卫星与地面接收机的几何距离；

为导航卫星与地面接收机的几何距离；

为第一观测方程中的整周模糊度和

为第二观测方程中的整周模糊度；dt_BD,i为地面接收机钟差，dt_LEO,i为低轨卫星接收机钟差；dt_i ^p为卫星钟差；

为地面与导航卫星间的对流层延迟，

为低轨卫星与导航卫星的多路径效应；c为真空中的光速。

5.一种卫星定轨装置，其特征在于，该装置包括：

收发单元，用于获取地面接收机确定的导航卫星的第一观测数据和低轨卫星接收机确定的所述导航卫星的第二观测数据；所述第一观测数据和所述第二观测数据均为双频观测数据；所述第一观测数据和所述第二观测数据均包括载波相位观测数据；所述载波相位观测数据的整数部分和小数部分分开处理，以完整地保留载波相位中所包含的模糊度信息，包括：将上一时刻的载波相位累加值、当前时刻的载波相位累加值、中频累加值和本地时钟调整量，输入至载波相位计算单元，确定出载波相位观测数据的整数部分和小数部分；对载波相位观测数据进行半周补偿，输出最终的载波相位观测数据；

处理单元，用于确定所述地面接收机的第一观测方程和所述低轨卫星接收机的第二观测方程；根据双频观测数据，消除所述第一观测方程和所述第二观测方程中的电离层延迟误差，得到消除后的第一观测方程和第二观测方程；根据所述第一观测数据和所述第二观测数据，解算所述消除后的第一观测方程和所述消除后的第二观测方程，从而确定所述导航卫星的轨道位置。

6.如权利要求5所述的卫星定轨装置，其特征在于，所述处理单元，用于若确定所述低轨卫星接收机的载波跟踪环路的跟踪状态为锁定状态，则根据测量精度确定所述载波跟踪环路的滤波器组合；将获得的导航卫星的信号输入至所述载波跟踪环路的滤波器组合，并将输出的载波相位数据作为所述第二观测数据中的载波相位观测数据。

7.如权利要求5所述的卫星定轨装置，其特征在于，

所述第一观测方程为：

Y_GROUNDi＝F_GB(X_BDi,X_oi,t_i)+ξ_GROUNDi；其中,Y_GROUNDi为在t_i时刻的第一观测数据；F_GB表示地面接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_oi表示第一观测方程中的观测模型的参数；ξ_GROUNDi为所述第一观测数据的观测误差；

所述第二观测方程为：

Y_LEOi＝F_LB(X_BDi,X_LEOi,X_Oi,t_i)+ξ_LEOi；

其中,Y_LEOi为在t_i时刻的第二观测数据；F_LB为低轨卫星接收机在t_i时刻的观测函数；X_BDi为导航卫星的轨道位置；X_LEOi为低轨卫星的轨道位置；X_Oi为第二观测方程中的观测模型的参数；ξ_LEOi为所述第二观测数据的观测误差。

8.一种电子设备，其特征在于，包括：

至少一个处理器；与所述至少一个处理器通信连接的存储器；

其中，所述存储器存储有可被所述至少一个处理器执行的指令，所述指令被所述至少一个处理器执行权利要求1～4任一权利要求所述方法。

9.一种计算机程序产品，其特征在于，所述计算机程序产品包括存储在计算机可读存储介质上的计算程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，使所述计算机执行权利要求1～4任一权利要求所述方法。

Images