CN108732597B - 一种多卫星导航***的时间基准建立方法及*** - Google Patents
一种多卫星导航***的时间基准建立方法及*** Download PDFInfo
- Publication number
- CN108732597B CN108732597B CN201810564951.XA CN201810564951A CN108732597B CN 108732597 B CN108732597 B CN 108732597B CN 201810564951 A CN201810564951 A CN 201810564951A CN 108732597 B CN108732597 B CN 108732597B
- Authority
- CN
- China
- Prior art keywords
- satellite
- clock
- navigation system
- time
- satellite navigation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
- G01S19/235—Calibration of receiver components
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种多卫星导航***的时间基准建立方法及***,所述方法包括:获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据地面站GNSS观测数据和星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;根据地面站观测模型和星载GNSS观测模型,对多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;根据各个卫星导航***的钟差构建多卫星导航***的时间尺度;将时间尺度约束到某一卫星导航***的***时间下,生成多卫星导航***中各个卫星的统一时间基准。本发明能够基于多卫星导航***星载原子钟,实现多卫星导航***的连续、稳定的时间基准的建立和维护。
Description
技术领域
本发明涉及卫星导航技术领域,尤其涉及一种基于多卫星导航***星载原子钟的多卫星导航***的时间基准建立方法及***。
背景技术
时间作为组成卫星导航***的三大基础要素(时间、轨道、信号)之一,对卫星导航***起支撑作用,是维护***正常运行秩序和工作流程的重要保障,在很大程度上决定了卫星导航***提供定位、授时、测速的功能及性能。
导航***的时间基准作为卫星导航***的统一时间参考,是根据守时钟组的对比数据产生的一种时间尺度,导航***需要将每个星载原子钟的时间统一到时间基准上,发射给用户使用。每个导航***都建立了专门的时间基准(***时间)。目前,卫星导航***时间基准包括如下两种定义方式:
1)主钟:由主控站的主钟定义,如GLONASS的***时间GLONASST,由地面监控部分配置的若干台高精度的原子钟,经过数据对比到GLONASST,溯源至俄罗斯国家标准时间UTC(SU),GLONASS卫星时间由星载原子钟保持,综合控制中心将其与GLONASST每天对比两次,将星钟修正值传至卫星;北斗的***时间BDT也是一种主钟定义的方式,由位于北斗二代地面运控***主控站的时频***建立并维持,溯源至中科院国家授时中心UTC(NTSC),采用卫星双向时间传递技术传给每一监测站钟,用无线电双向时间比对技术传递给每一卫星钟,保证***的时间同步。
2)纸面钟:是一种为合成钟方式,由所有地面钟和卫星钟组成的钟组定义,***时间尺度由各个钟的加权平均得到合成钟。如GPS的***时间为GPST,通过Kalman滤波算法,由地面主控站对监测站的高精度原子钟和卫星的星载原子钟的组合钟进行共同加权,建立和维护GPST,溯源至美国海军天文台UTC(USNO)。Galileo的时间***GST同样采用组合钟定义,由地面控制中心接受来自测站的观测数据,通过共视法获得UTC,滤波产生钟差改正数和平均频率,平均频率作用于Galileo***主钟,用于***基准GST的建立。而且国际GNSS服务组织(IGS)的时间尺度建立也是一种合成钟的方式,利用几个稳定的测站和GPS卫星钟通过加权平均得到一个内部的时间基准,并约束到GPST上,使得最终综合的钟差产品具有相同的时间基准。
由于不同的卫星导航***在设计理念、时间、轨道、频率等方面存在着差异,导致现有建立时间基准的方法存在以下几个问题:
1)各个***确定的***时间存在时间偏差,导致最终求解出的卫星钟差也存在时间偏差;
2)以上各***在建立和维护各自***时间基准时并未利用低轨卫星的星载原子钟,在卫星端只用到导航卫星星载原子钟,无法适用于包括各个卫星导航***和低轨卫星的多卫星导航***。
发明内容
鉴于上述问题,本发明提出了一种多卫星导航***的时间基准建立方法及***,能够基于多卫星导航***星载原子钟,实现连续、稳定的时间基准的建立和维护。
本发明的一个方面,提供了一种多卫星导航***的时间基准建立方法,包括:
获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;
根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;
根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度;
将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差。
可选地,所述根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差,包括:
根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理;
采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数;
将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
可选地,所述根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度,包括:
对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值;
构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程;
将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程;
根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程;
根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列;
根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。
可选地,所述根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列包括:
对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列。
可选地,所述将所述时间尺度约束到某一卫星导航***的***时间下,形成所述多卫星导航***中各个卫星的统一时间基准,包括:
构建约束算法的状态方程和观测方程;
根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
本发明的另一个方面,提供了一种多卫星导航***的时间基准建立***,包括:
模型构建单元,用于获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;
钟差获取单元,用于根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;
时间尺度建立单元,用于根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度;
时间基准生成单元,用于将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差。
可选地,所述钟差获取单元,包括:
线性化处理模块,用于根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理;
第一计算模块,用于采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数;
约束处理模块,用于将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
可选地,所述时间尺度建立单元,包括:
相位转换模块,用于对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值;
第一模型构建模块,用于构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程;
第二模型构建模块,用于将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程;
第三模型构建模块,用于根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程;
第二计算模块,用于根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列;
时间尺度构建模块,用于根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。
可选地,所述时间尺度构建模块,具体用于对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列。
可选地,所述时间基准生成单元,包括:
第四模型构建模块,用于构建约束算法的状态方程和观测方程;
第三计算模块,用于根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
本发明实施例提供的多卫星导航***的时间基准建立方法及***,利用涵盖所有***的地面站网GNSS和星载GNSS观测数据建立时间基准时,保证了结果的稳定性和可靠性;在时间基准的建立过程中,利用所有星载原子钟形成原子钟组,并通过加权来维持,保证时间基准短周期项更高的稳定度和精度;通过将时间基准最终约束到某一个卫星导航***时间下,保证了时间基准长周期项的稳定,而且最终获取的各个卫星的钟差具有统一的时间基准,方便了多卫星导航***数据的融合处理及应用。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1为本发明实施例的多卫星导航***的时间基准建立方法的流程图;
图2为本发明实施例的多卫星导航***的时间基准建立方法中步骤S12的细分流程图;
图3为本发明实施例的多卫星导航***的时间基准建立方法中步骤S13的细分流程图;
图4为本发明实施例的多卫星导航***的时间基准建立***的结构示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
本技术领域技术人员可以理解,除非特意声明,这里使用的单数形式“一”、“一个”、“所述”和“该”也可包括复数形式。应该进一步理解的是,本发明的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件,但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。
本技术领域技术人员可以理解,除非另外定义,这里使用的所有术语(包括技术术语和科学术语),具有与本发明所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是,诸如通用字典中定义的那些术语,应该被理解为具有与现有技术的上下文中的意义一致的意义,并且除非被特定定义,否则不会用理想化或过于正式的含义来解释。
图1示意性示出了本发明一个实施例的多卫星导航***的时间基准建立方法的流程图。参照图1,本发明实施例的多卫星导航***的时间基准建立方法具体包括以下步骤:
S11、获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型。
地面站的GNSS观测方程为:
星载GNSS观测方程为:
式中,f为频率,为载波相位观测值(以距离表示),s表示导航卫星,k表示测站或低轨卫星,为伪距观测值,为模糊度,为地面站或低轨卫星与导航卫星几何距离,dtk为测站或低轨卫星钟差,·dts为导航卫星钟差,为电离层对相位和伪距的影响,为对流层延迟, 为其他各种改正(包括多路径、天线相位中心、地球固体潮、海洋负荷潮、相对论效应等)对相位和伪距的影响。
S12、根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差。
S13、根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度。
S14、将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差,实现多卫星导航***的统一时间基准的卫星钟差。
本发明实施例提供的多卫星导航***的时间基准建立方法,利用涵盖所有***的地面站网GNSS和星载GNSS观测数据建立时间基准时,保证了结果的稳定性和可靠性;在时间基准的建立过程中,利用所有星载原子钟形成原子钟组,并通过加权来维持,保证时间基准短周期项更高的稳定度和精度;通过将时间基准最终约束到某一个卫星导航***时间下,保证了时间基准长周期项的稳定,而且最终获取的各个卫星的钟差具有统一的时间基准,方便了多卫星导航***数据的融合处理及应用。
本发明实施例中,如图2所示,步骤S12中的根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差,具体包括以下步骤:
S121、根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理。具体包括:根据预处理后的中高轨GNSS卫星地面跟踪数据和低轨卫星的GNSS卫星观测数据分别生成中高轨GNSS卫星到低轨卫星及地面站的观测距离O;根据预设的参数初始值计算中高轨GNSS卫星到低轨卫星及地面站的几何距离;对计算出来的几何距离加以修正,得到中高轨GNSS卫星到低轨卫星及地面站的计算距离C;将中高轨GNSS卫星到地面站的观测距离O与计算距离C进行差值运算,生成地面站观测模型的先验残差,并将中高轨GNSS卫星到低轨卫星的所述观测距离O与所述计算距离C进行差值运算,生成星载GNSS观测模型的先验残差;根据各个观测模型中参数向量的近似值,采用观测模型对相应的参数向量计算偏导数,得到地面站观测模型的第一信息矩阵和星载GNSS观测模型的第二信息矩阵;根据所述地面站观测模型的先验残差和所述第一信息矩阵构建第一观测模型对应的线性化后的观测方程,以及根据所述星载GNSS观测模型的先验残差和所述第二信息矩阵构建第二观测模型对应的线性化后的观测方程。
S122、采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数。
S123、将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
在多GNSS导航卫星和低轨卫星融合精密定轨和时间同步中采用简化动力学定轨,估计得出各卫星轨道参数和力模型参数等状态向量、以及地面站、低轨星载钟差参数、导航卫星钟差参数等,同时将各导航卫星的钟差参数约束到各自的导航星历上,即
本发明实施例中,如图3所示,步骤S13中根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度,具体包括以下步骤:
S131、对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值。
S132、构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程。
单个卫星钟的状态方程可以写成如下形式:
其中
单个卫星钟的观测方程可以写成如下形式:
式中,n表示测量噪声,一般为白噪声,其方差r为:
S133、将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程。
以所有***(包括各个导航***和低轨卫星)的N个卫星钟构成钟组,则钟组的状态方程为:
权重的计算公式为:
S134、根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程。
S135、根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列。
S136、根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。进一步地,所述根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列包括:对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列其中,GPST为某一导航卫星***的卫星钟时间,NEWT为所述***时间。
本发明实施例中,步骤S14中的将所述时间尺度约束到某一卫星导航***的***时间下,形成所述多卫星导航***中各个卫星的统一时间基准,具体包括以下步骤:
构建约束算法的状态方程和观测方程;
根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
约束到某一***时间可以利用该***时间基准稳定的长期项,使新的基准的短周期项更加稳定,同样以约束到ORIT为例,约束算法的状态方程和观测方程可以表示为:
得到多卫星导航***的统一时间基准之后,每一卫星导航***将统一的时间基准附加到原始钟差上,获取具有相同时间基准的卫星钟差,实现多***的统一时空基准的卫星钟差;
将原始的钟差归算到NEWT,即
即可得到具有统一时间基准的卫星钟差。
本发明实施例建立的多卫星导航***时间基准,还通过实时的计算更新,保证时间基准长短周期项的稳定性和可靠性。
对于方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明实施例并不受所描述的动作顺序的限制,因为依据本发明实施例,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本发明实施例所必须的。
图4示意性示出了本发明一个实施例的多卫星导航***的时间基准建立***的结构示意图。参照图4,本发明实施例的多卫星导航***的时间基准建立***具体包括模型构建单元401、钟差获取单元402、时间尺度建立单元403以及时间基准生成单元404,其中:
模型构建单元401,用于获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;
钟差获取单元402,用于根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;
时间尺度建立单元403,用于根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度;
时间基准生成单元404,用于将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差。
在本发明的一个可选实施例中,所述钟差获取单元402,包括线性化处理模块、第一计算模块以及约束处理模块,其中:
线性化处理模块,用于根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理;
第一计算模块,用于采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数;
约束处理模块,用于将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
在本发明的一个可选实施例中,所述时间尺度建立单元403,包括相位转换模块、第一模型构建模块、第二模型构建模块、第三模型构建模块、第二计算模块以及时间尺度构建模块,其中:
相位转换模块,用于对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值;
第一模型构建模块,用于构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程;
第二模型构建模块,用于将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程;
第三模型构建模块,用于根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程;
第二计算模块,用于根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列;
时间尺度构建模块,用于根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。
进一步地,所述时间尺度构建模块,具体用于对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列。
在本发明的一个可选实施例中,所述时间基准生成单元404,包括第四模型构建模块和第三计算模块,其中:
第四模型构建模块,用于构建约束算法的状态方程和观测方程;
第三计算模块,用于根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
对于***实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本发明实施例提供的多卫星导航***的时间基准建立方法及***,利用涵盖所有***的地面站网GNSS和星载GNSS观测数据建立时间基准时,保证了结果的稳定性和可靠性;在时间基准的建立过程中,利用所有星载原子钟形成原子钟组,并通过加权来维持,保证时间基准短周期项更高的稳定度和精度;通过将时间基准最终约束到某一个卫星导航***时间下,保证了时间基准长周期项的稳定,而且最终获取的各个卫星的钟差具有统一的时间基准,方便了多卫星导航***数据的融合处理及应用。
以上所描述的装置实施例仅仅是示意性的,其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下,即可以理解并实施。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件。基于这样的理解,上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在计算机可读存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
此外,本领域的技术人员能够理解,尽管在此的一些实施例包括其它实施例中所包括的某些特征而不是其它特征,但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如,在下面的权利要求书中,所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (8)
1.一种多卫星导航***的时间基准建立方法,其特征在于,包括:
获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;
根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;
根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度;
将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差;
其中,所述根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度,包括:
对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值;
构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程;
将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程;
根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程;
根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列;
根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。
2.根据权利要求1所述的方法,其特征在于,所述根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差,包括:
根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理;
采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数;
将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
3.根据权利要求1所述的方法,其特征在于,所述根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列包括:
对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列。
4.根据权利要求1所述的方法,其特征在于,所述将所述时间尺度约束到某一卫星导航***的***时间下,形成所述多卫星导航***中各个卫星的统一时间基准,包括:
构建约束算法的状态方程和观测方程;
根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
5.一种多卫星导航***的时间基准建立***,其特征在于,包括:
模型构建单元,用于获取多卫星导航***的地面站GNSS观测数据和星载GNSS观测数据,分别根据所述地面站GNSS观测数据和所述星载GNSS观测数据构建对应的地面站观测模型和星载GNSS观测模型;
钟差获取单元,用于根据所述地面站观测模型和星载GNSS观测模型,对所述多卫星导航***中的导航卫星以及低轨卫星进行联合精密定轨和时间同步,获取各个卫星导航***的钟差;
时间尺度建立单元,用于根据所述各个卫星导航***的钟差构建多卫星导航***的时间尺度;
时间基准生成单元,用于将所述时间尺度约束到某一卫星导航***的***时间下,生成所述多卫星导航***中各个卫星的统一时间基准,以供每一卫星导航***将所述统一时间基准附加到各自对应的原始钟差上,以获取各自在所述统一时间基准下的卫星钟差;
所述时间尺度建立单元,包括:
相位转换模块,用于对所述各个卫星导航***的钟差进行相位转换,得到各个卫星导航***的频率观测值;
第一模型构建模块,用于构建每一卫星钟的时钟状态方程,基于所述时钟状态方程构建时钟观测方程;
第二模型构建模块,用于将各个导航卫星以及低轨卫星对应的卫星钟组合成钟组,根据每一卫星钟的时钟状态方程构建所述钟组的状态方程;
第三模型构建模块,用于根据所述钟组的状态方程以及所述各个卫星导航***的频率观测值构建所述钟组的观测方程;
第二计算模块,用于根据所述钟组的观测方程,计算各个卫星钟相对于欲约束到的卫星导航***的***时间的频率序列;
时间尺度构建模块,用于根据所述频率序列计算各个卫星钟相对于所述***时间的时间序列,并根据所述时间序列构建所述时间尺度。
6.根据权利要求5所述的***,其特征在于,所述钟差获取单元,包括:
线性化处理模块,用于根据预设的参数初始值分别对所述地面站观测模型和星载GNSS观测模型进行线性化处理;
第一计算模块,用于采用最小二乘法计算线性化处理后的所述地面站观测模型和星载GNSS观测模型,得到各个卫星导航***的钟差参数;
约束处理模块,用于将所述各个卫星导航***的钟差参数约束到各自对应的导航星历。
7.根据权利要求5所述的***,其特征在于,所述时间尺度构建模块,具体用于对所述频率序列进行积分运算得到相应的卫星钟时间相对于所述***时间的时间序列。
8.根据权利要求5所述的***,其特征在于,所述时间基准生成单元,包括:
第四模型构建模块,用于构建约束算法的状态方程和观测方程;
第三计算模块,用于根据所述约束算法的状态方程和观测方程,采用卡尔曼滤波法计算约束到某一卫星导航***的***时间下的时间序列,将该时间序列作为所述多卫星导航***中各个卫星的统一时间基准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810564951.XA CN108732597B (zh) | 2018-06-04 | 2018-06-04 | 一种多卫星导航***的时间基准建立方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810564951.XA CN108732597B (zh) | 2018-06-04 | 2018-06-04 | 一种多卫星导航***的时间基准建立方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108732597A CN108732597A (zh) | 2018-11-02 |
CN108732597B true CN108732597B (zh) | 2020-10-02 |
Family
ID=63931903
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810564951.XA Active CN108732597B (zh) | 2018-06-04 | 2018-06-04 | 一种多卫星导航***的时间基准建立方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108732597B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109520512A (zh) * | 2018-12-28 | 2019-03-26 | 上海海积信息科技股份有限公司 | 一种卫星精密定轨方法及装置 |
CN109870898B (zh) * | 2019-02-27 | 2020-11-17 | 武汉大学 | 基于ppp的gnss授时接收机时钟组合调控方法及*** |
CN110045395B (zh) * | 2019-05-17 | 2021-10-08 | 北京和协航电信息科技有限公司 | 一种基于导航卫星信号的检靶方法、***和设备 |
CN111948686B (zh) * | 2020-08-05 | 2024-04-26 | 航天恒星科技有限公司 | 时间同步方法及装置 |
CN112325879A (zh) * | 2020-11-03 | 2021-02-05 | 中国电子科技集团公司信息科学研究院 | 一种基于多源传感器集成的仿生复合导航授时微*** |
CN112445120B (zh) * | 2020-11-27 | 2022-04-22 | 西安空间无线电技术研究所 | 一种分布式无中心天基时间基准建立与保持*** |
CN113253311B (zh) * | 2021-04-06 | 2024-02-06 | 探索数据科技(深圳)有限公司 | 联合卫星导航方法、***、电子设备及存储介质 |
CN113253314B (zh) * | 2021-05-11 | 2023-06-23 | 中国科学院国家授时中心 | 一种低轨卫星间时间同步方法及*** |
CN113341685A (zh) * | 2021-06-30 | 2021-09-03 | 北京微纳星空科技有限公司 | 一种卫星时间管理方法、电子设备及介质 |
CN113608427B (zh) * | 2021-07-09 | 2022-07-05 | 中国科学院国家授时中心 | 一种集中式天基时间基准建立方法 |
CN113640838B (zh) * | 2021-08-19 | 2024-06-18 | 中国电力科学研究院有限公司 | 一种多***联合卫星导航定位授时装置及方法 |
CN115639743B (zh) * | 2022-10-19 | 2023-05-30 | 中国科学院国家授时中心 | 一种基于整网时间比对的天基时间基准建立方法及*** |
CN116299617B (zh) * | 2023-02-16 | 2023-10-31 | 北京航空航天大学 | 一种gnss星地原子钟实时动态模型的建立方法 |
CN117377057B (zh) * | 2023-12-08 | 2024-02-13 | 中国科学院国家授时中心 | 一种低轨卫星和地面站的星地时间同步方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5963167A (en) * | 1996-03-13 | 1999-10-05 | California Institute Of Technology | Analyzing system for global positioning system and general satellite tracking |
CN102937710A (zh) * | 2012-11-13 | 2013-02-20 | 北京空间飞行器总体设计部 | 一种确定原子钟非理想特性引起用户测距误差的方法 |
CN103235321A (zh) * | 2013-03-29 | 2013-08-07 | 东南大学 | Gps单站伪距定位精密授时方法 |
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
CN105867108A (zh) * | 2016-06-20 | 2016-08-17 | 中国科学院国家授时中心 | 一种国家标准时间远程复现方法 |
CN106443739A (zh) * | 2016-09-09 | 2017-02-22 | 清华大学 | 辅助增强导航方法及设备 |
CN107153209A (zh) * | 2017-07-06 | 2017-09-12 | 武汉大学 | 一种短弧段低轨导航卫星实时精密定轨方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7312747B2 (en) * | 2005-10-03 | 2007-12-25 | Trimble Navigation Limited | Multiple-GNSS and FDMA high precision carrier-phase based positioning |
US20150212211A1 (en) * | 2014-01-24 | 2015-07-30 | Qualcomm Incorporated | Methods and systems for multi-gnss operation |
-
2018
- 2018-06-04 CN CN201810564951.XA patent/CN108732597B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5963167A (en) * | 1996-03-13 | 1999-10-05 | California Institute Of Technology | Analyzing system for global positioning system and general satellite tracking |
CN102937710A (zh) * | 2012-11-13 | 2013-02-20 | 北京空间飞行器总体设计部 | 一种确定原子钟非理想特性引起用户测距误差的方法 |
CN103235321A (zh) * | 2013-03-29 | 2013-08-07 | 东南大学 | Gps单站伪距定位精密授时方法 |
CN105158780A (zh) * | 2015-07-24 | 2015-12-16 | 北京跟踪与通信技术研究所 | 一种基于多种导航卫星可互换的导航定位方法 |
CN105867108A (zh) * | 2016-06-20 | 2016-08-17 | 中国科学院国家授时中心 | 一种国家标准时间远程复现方法 |
CN106443739A (zh) * | 2016-09-09 | 2017-02-22 | 清华大学 | 辅助增强导航方法及设备 |
CN107153209A (zh) * | 2017-07-06 | 2017-09-12 | 武汉大学 | 一种短弧段低轨导航卫星实时精密定轨方法 |
Non-Patent Citations (6)
Title |
---|
A control theory approach to clock steering techniques;Marcello Farina等;《IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control》;20101231;第57卷(第10期);全文 * |
Local representations of UTC in national laboratories;Peter B Whibberley等;《Metrologia》;20111231;第48卷(第4期);全文 * |
伍贻威等.原子钟两级驾驭算法及在建立GNSS基准中的应用.《电子学报》.2016,第44卷(第7期), * |
伽利略***GSTB-V1的时间尺度算法与实现;郑子扬等;《舰船电子工程》;20081231;第28卷(第10期);全文 * |
原子钟两级驾驭算法及在建立GNSS基准中的应用;伍贻威等;《电子学报》;20160731;第44卷(第7期);第1节引言 * |
原子钟时间尺度与Kalman滤波器;卫国;《陕西天文台台刊》;19901231;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108732597A (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108732597B (zh) | 一种多卫星导航***的时间基准建立方法及*** | |
CN111045034B (zh) | 基于广播星历的gnss多***实时精密时间传递方法及*** | |
Strasser et al. | Processing of GNSS constellations and ground station networks using the raw observation approach | |
Haines et al. | One-centimeter orbit determination for Jason-1: new GPS-based strategies | |
CN103221839B (zh) | 使用区域增强消息的gnss信号处理 | |
CN102498414B (zh) | 用以估计轨道的gnss信号处理 | |
Wilson et al. | New and improved-the broadcast interfrequency biases | |
US8120530B2 (en) | Method and device for using GNSS satellite trajectory extension data in mobile apparatus | |
CN108120994B (zh) | 一种基于星载gnss的geo卫星实时定轨方法 | |
WO2017160619A1 (en) | Navigation satellite orbit and low latency clock determination with wide-lane and narrow-lane bias corrections | |
CN103502843A (zh) | 具有三角相的gnss信号处理 | |
CN117388881B (zh) | 一种低轨卫星的星载原子钟向UTC(k)的溯源方法及*** | |
JP2007529010A (ja) | 2周波数の一方で測定データが利用できない場合に短期間バックアップ2周波数ナビゲーションを行なう方法 | |
WO2017160618A1 (en) | Navigation satellite wide-lane bias determination and over-range adjustment system and method | |
Li et al. | Precise orbit determination for the FY-3C satellite using onboard BDS and GPS observations from 2013, 2015, and 2017 | |
CN113204042B (zh) | 一种基于精密单点定位的多星座联合列车定位方法 | |
Hwang et al. | GPS‐Based Orbit Determination for KOMPSAT‐5 Satellite | |
Kuang et al. | Real-time GPS satellite orbit and clock estimation based on OpenMP | |
CN110554443A (zh) | 基于载波相位观测值和点加速度法确定地球重力场的方法 | |
Gong et al. | Precise real-time navigation of the small TJU-1 satellite using GPS, GLONASS and BDS | |
Jin et al. | Analysis of a federal Kalman filter-based tracking loop for GPS signals | |
CN116299586B (zh) | 基于广播星历的精密单点定位方法、接收机、设备和介质 | |
WO2020231641A1 (en) | Methods and enhancing non-global navigation satellite system location and timing pseudorange positioning calculations and systems thereof | |
CN116859420A (zh) | 一种低轨卫星增强全球卫星导航***(LeGNSS)数据仿真的方法 | |
CN110988932A (zh) | 一种提高星载gps接收机实时钟差解算精度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |