CN114755702A - 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 - Google Patents
基于ppp的精密卫星星历时间基准稳定度评估方法及系统 Download PDFInfo
- Publication number
- CN114755702A CN114755702A CN202210599663.4A CN202210599663A CN114755702A CN 114755702 A CN114755702 A CN 114755702A CN 202210599663 A CN202210599663 A CN 202210599663A CN 114755702 A CN114755702 A CN 114755702A
- Authority
- CN
- China
- Prior art keywords
- ppp
- time
- station
- stability
- noise
- 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.)
- Granted
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims abstract description 22
- 238000012360 testing method Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000003089 Pariser Parr Pople method Methods 0.000 claims description 33
- 238000000034 method Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 16
- 230000000875 corresponding effect Effects 0.000 claims description 11
- 230000007774 longterm Effects 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000012805 post-processing Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000002596 correlated effect Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 claims description 2
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000007619 statistical method Methods 0.000 abstract description 2
- 238000004590 computer program Methods 0.000 description 4
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 3
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 239000013078 crystal Substances 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 238000013097 stability assessment Methods 0.000 description 2
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- 101100522111 Oryza sativa subsp. japonica PHT1-11 gene Proteins 0.000 description 1
- 101100545228 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) ZDS1 gene Proteins 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明提供一种基于PPP的精密卫星星历时间基准稳定度评估方法及系统,包括使用外接原子钟组频率源的IGS观测站的观测数据与不同分析中心的实时产品进行事后PPP处理,得到不同分析中心产品对应的观测站单站授时结果;选取一个测站作为中心测站,进行站间时间同步,得到时间同步结果;计算单站授时与站间时间同步结果的重叠Allan方差,将Allan方差结果线性组合,分离产品基准稳定度与PPP估计噪声稳定度。本发明可仅使用PPP技术通过数据处理统计分析实现,无需硬件设备;通过PPP技术与Allan方差结果的线性组合,即可得到PPP产品的时间基准稳定度,为授时接收机的PPP产品选择提供参考。
Description
技术领域
本发明属于全球导航卫星系统(Global Navigation Satellite System,GNSS)技术领域,特别涉及一种基于精密单点定位技术(Precise Point Positioning,PPP)的精密卫星星历时间基准稳定度的评估方法及系统。
背景技术
全球导航卫星系统(GNSS)是一种星基无线电定位系统,主要包括美国全球定位系统(Global Positioning System,GPS)、中国的北斗(BDS)、俄罗斯的格洛纳斯(GLONASS),以及欧洲的伽利略(Galileo)四大系统。在授时领域,利用GNSS接收机进行授时具有精度高、成本低和稳定等特点,授时接收机将在授时领域得到越来越广泛的运用。
GNSS授时接收机进行单向授时,其精度取决于卫星星历产品的精度。使用伪距单点定位技术,结合广播星历,其钟差的估计精度在20ns左右,单向授时精度在20ns~50ns;使用PPP技术,其钟差的估计精度在0.2ns左右,授时精度小于2ns。GNSS授时接收机使用实时PPP数据进行实时高精度单向授时,避免了传统PPP产品发布延迟的缺点。这种授时接收机从网络端接收到的实时卫星钟差改正数信息,其输出的时频信号短期稳定度取决于本地晶振或者原子钟,长期稳定度取决于接收的卫星改正数产品时间基准稳定度。不同分析中心播发的实时PPP改正数具有不同的时间基准稳定度,对不同分析中心的产品时间基准稳定度进行评估,可以为GNSS授时接收机运行时对实时PPP产品的选择提供参考,从而输出长期稳定度更好的时频信号。但是,本领域当前还缺乏有效的评估技术手段。
发明内容
针对基于PPP技术的GNSS授时接收机,其接收到的实时PPP产品的时间基准稳定度会影响到其输出的时频信号稳定度的情况。本发明提供了一种适用于基于PPP技术的精密卫星星历时间基准稳定度的评估方案,通过对不同分析中心产品的后处理,进行模拟单向授时与站间时间同步,进而计算出不同分析中心产品的时间基准稳定度。与传统的基于PPP单向授时的时间基准稳定度评估相比,本发明进一步的剔除了其中的PPP估计噪声稳定度,使得到的结果更具参考性。
为解决上述技术问题,本发明采用如下的技术方案:
一种基于PPP的精密卫星星历时间基准稳定度评估方法,包括以下步骤:
步骤S1,使用外接原子钟组频率源的IGS观测站的观测数据与不同分析中心的实时产品进行事后PPP处理,得到不同分析中心产品对应的观测站单站授时结果;选取一个测站作为中心测站,进行站间时间同步,得到时间同步结果;
步骤S2,计算单站授时与站间时间同步结果的重叠Allan方差,将Allan方差结果线性组合,分离产品基准稳定度与PPP估计噪声稳定度,实现方式包括以下子步骤,
步骤S2.1,计算钟差序列的Allan方差,包括将此Allan方差看作产品时间基准的稳定度与PPP估计噪声的线性组合;
步骤S2.2,计算时间同步结果序列的Allan方差,包括将每个观测站的PPP估计噪声视为具有相同噪声水平系数与成分,则将此Allan方差看作是系数为2的PPP估计噪声Allan方差的线性组合;
步骤S2.3,将S2.1与S2.2分别得到的Allan方差所表征的线性组合组成线性多项式,其待估参数分别为产品时间基准稳定度与PPP估计噪声,求解出两个待估参数,得到产品时间基准稳定度。
而且,步骤S1的实现过程包括以下子步骤,
步骤S1.1,选择IGS中GNSS接收机外接原子钟组的观测站,得到观测数据;
步骤S1.2,存储不同分析中心的改正数产品,组合广播星历与改正数进行精密星历恢复,得到精密星历;
步骤S1.3,使用观测数据与恢复后的精密星历进行PPP事后处理,得到的钟差序列看作包括产品时间基准、接收机钟差和PPP估计噪声的线性组合;
步骤S1.4,以S1.1中选择的观测站中的某一个站作为中心站,建立时间链路,进行站间时间同步,计算得到站间时间同步结果。
而且,步骤S1.4中,进行站间时间同步时,将钟差序列对应时标作差,计算得到站间时间同步结果,实现方式如下,
其中,ta为测站a的接收机真实钟差,tb为测站b的接收机真实钟差,εa为测站a的接收机钟差估计噪声,εb为测站b的接收机钟差估计噪声,为测站a与测站b的之差,长时间尺度上可视为定值。Δtab即为时间同步结果。
而且,步骤S2.1实现方式如下,
S(x)=xi+2m-2xi+m+xi (五)
其中,τ=mτ0表示平均时间,m表示平均因子,τ0表示测量间隔,N表示观测值数量,x表示N个以τ0为间隔的相位观测值;xi表示表示测量时间为i的相位数据,xi+2m、xi+m分别表示在xi之后,与xi时间间隔为m与2m的相位数据;S(x)表示相位数据的三差算子;表示平均时间为τ,计算是使用的相位数据为x的Allan方差函数;
其中,表示以接收机估计钟差序列计算的Allan方差,表示卫星产品时间基准的Allan方差,表示PPP估计噪声的Allan方差,S(t0)表示产品时间基准t0的三差算子结果,S(εr)表示PPP估计噪声εr的三差算子结果,表示产品时间基准与PPP估计噪声的Allan方差交叉项。
而且,步骤S2.2实现方式如下,
将式三代入式四中,当tr具有高稳定性,在长时间尺度上视为定值;考虑使用相同分析中心产品时,不同测站观测数据进行PPP解算的估计噪声具有相同稳定度,记为r为测站标号;相应地,S(ε)表示PPP估计噪声ε的三差算子结果,S(εa)表示测站a的PPP估计噪声εa的三差算子结果,S(εb)表示测站b的PPP估计噪声εb的三差算子结果;通过模拟调频白噪声与调相白噪声序列,得知 所以有:
而且,步骤S2.3实现方式如下,
其中,uk=(1 1 … 1)T为k×1的1值向量,zk=(0,0 ... 0)T为k×1的0值向量,uk-1=(1 1 … 1)T为(k-1)×1的1值向量,zk-1=(0,0 ... 0)T为(k-1)×1的0值向量;
表示测站r1的PPP接收机估计钟差序列的Allan方差,...,表示测站rk的PPP接收机估计钟差序列的Allan方差,表示测站r1与r2的时间同步结果序列的Allan方差,...,表示表示测站r1与rk的时间同步结果序列的Allan方差通过式八,求解PPP产品的时间基准稳定度,同时得到PPP钟差估算噪声稳定度。
而且,用于GNSS授时接收机运行时,根据产品时间基准稳定度选择实时PPP产品,从而输出长期稳定度更好的时频信号。
另一方面,本发明提供一种基于PPP的精密卫星星历时间基准稳定度评估系统,用于实现如上所述的一种基于PPP的精密卫星星历时间基准稳定度评估方法。
本发明公开了一种基于精密单点定位(Precise Point Positioning,PPP)技术的精密卫星星历时间基准稳定度的评估方案。当GNSS授时接收机进行授时时,使用的实时PPP产品的时间基准稳定度会间接影响其时频输出的长期稳定度。本发明首先利用公开发布的观测数据,使用恢复的实时PPP产品进行PPP解算,得到测站钟差序列,即为单向授时结果。可以对多个不同测站重复此过程,选取其中一个测站为中心,建立时间同步链路,获取时间同步结果。对单向授时结果与时间同步结果求取重叠Allan方差,将Allan方差结果进行线性组合,即可分别求取PPP接收机钟差估计噪声稳定度与PPP产品时间基准稳定度。本发明可仅使用PPP解算通过数据处理统计分析实现,无需硬件设备。
通过PPP解算与Allan方差结果的线性组合,即可得到PPP产品的时间基准稳定度,为授时接收机的PPP产品选择提供参考。
附图说明
图1是本发明实施例原理的示意图。
具体实施方式
以下通过实施例结合附图说明本发明的技术方案。
本发明实施例使用IGS发布的公开观测数据与公开播发的实时PPP改正数,利用基于PPP的单向授时与站间时间同步,通过重叠Allan方差表征频率稳定度,对求取的重叠Allan方差进行线性近似,解算线性方程组,得到实时PPP产品的时间基准稳定度。通过对比不同产品的时间基准稳定度,可以为授时接收机提供产品选择的参考,提高授时接收机输出时频信号的长期稳定度。
参见图1,实施例提供一种基于PPP的精密卫星星历时间基准稳定度评估方法,所述方法中包含PPP解算、建立时间链路进行时间同步、Allan方差计算。
下面将以基于IGS公开发布的观测数据与精密星历为例,对本发明的实施例做进一步说明。
实施例提供的一种基于PPP的精密卫星星历时间基准稳定度评估方法,包括以下步骤:
步骤S1使用外接原子钟组频率源的IGS观测站的观测数据与不同分析中心的实时产品进行事后PPP处理,得到不同分析中心产品对应的观测站单站授时结果。选取一个测站作为中心测站,进行站间时间同步,得到时间同步结果。
进一步的,步骤1包括:
步骤S1.1获取GNSS观测数据。选择外接氢原子频标或原子钟组的IGS观测站,具体实施时也可使用自采集的外接氢原子频标或原子钟组的接收机。使用的观测数据,在采集时需要让接收机外接高稳定度频标,其目的是在GNSS观测数据中引用高稳定度参考频率源。目前外接氢原子频标并提供稳定观测数据的IGS观测站有BRUX,IENG,NIST,NRC1,OP71,OPMT,PTBB,PT11,ROAG,SPT0,TWTF,USN7,USN8,WAB2等;数据可通过https://cddis.nasa.gov/archive/下载获取。
步骤S1.2存储不同分析中心的改正数产品。不同分析中心播发的是基于广播星历的改正数数据,组合广播星历与改正数从而进行精密星历恢复,得到精密星历。
实施例中,通过BKG提供的BKG Ntrip Client软件(BNC,https://igs.BKG.bund.de/)进行实时改正数接收存储。通过参考《北斗卫星导航系统地基增强服务接口控制文件1.0》(http://www.csno-tarc.cn/support/documents)附录1与附录2将改正数恢复为精密星历。
步骤S1.3使用观测数据与恢复后的精密星历进行PPP事后处理,得到的钟差序列可以看作产品时间基准、接收机实际钟差、PPP估计噪声等的线性组合。
下载GNSS数据进行PPP解算,PPP为现有技术,本发明不予赘述。
步骤S1.4以步骤S1.1中选择的观测站中的某一个站作为中心站,建立时间链路,进行站间时间同步,即将钟差序列对应时标作差,计算得到站间时间同步结果。设测站a与测站b之间建立时间链路a-b,以时间链路a-b为例进行说明,将测站a的估计钟差与测站b的估计钟差对应时标作差:
其中,ta为测站a的接收机真实钟差,tb为测站b的接收机真实钟差,εa为测站a的接收机钟差估计噪声,εb为测站b的接收机钟差估计噪声,为测站a与测站b的之差,长时间尺度上可视为定值。Δtab即为时间同步结果。
步骤S2计算单站授时与站间时间同步结果的重叠Allan方差,将Allan方差结果线性组合,分离产品基准稳定度与PPP估计噪声稳定度。
进一步的,步骤2包括:
步骤S2.1计算钟差序列的Allan方差(阿伦方差),将此Allan方差可以看作产品时间基准的稳定度与PPP估计噪声的线性组合;
S(x)=xi+2m-2xi+m+xi (五)
其中,τ=mτ0表示平均时间,m表示平均因子,τ0表示测量间隔,N表示观测值数量,x表示N个以τ0为间隔的相位观测值;xi表示表示测量时间为i的相位数据,xi+2m、xi+m分别表示在xi之后,与xi时间间隔为m与2m的相位数据;S(x)表示相位数据的三差算子;表示平均时间为τ,计算是使用的相位数据为x的Allan方差函数。
其中,表示以接收机估计钟差序列计算的Allan方差,表示卫星产品时间基准的Allan方差,表示PPP估计噪声的Allan方差,S(t0)表示产品时间基准t0的三差算子结果,S(εr)表示PPP估计噪声εr的三差算子结果,表示产品时间基准与PPP估计噪声的Allan方差交叉项。
步骤S2.2计算时间同步结果序列的Allan方差,将每个观测站的PPP估计噪声视为具有相同噪声水平系数与成分,则此Allan方差可以看作是系数为2的PPP估计噪声Allan方差的线性组合;
将式三代入式四,考虑测站外接主动式氢钟,即tr具有高稳定性,在长时间尺度上可视为定值;考虑使用相同分析中心产品时,不同测站观测数据进行PPP解算的估计噪声具有相同稳定度,εr对不同测站观测数据可视为具有相同水平系数与成分的噪声ε,即有(测站标号r=a,b,…),相应地,S(ε)表示PPP估计噪声ε的三差算子结果,S(εa)表示测站a的PPP估计噪声εa的三差算子结果,S(εb)表示测站b的PPP估计噪声εb的三差算子结果;通过模拟调频白噪声与调相白噪声序列,可得知即测站a与测站b的PPP估计噪声交叉项的稳定度远小于估计噪声的稳定度,所以有:
步骤S2.3将步骤S2.1与步骤S2.2得到的Allan方差所表征的线性组合组成线性多项式,其待估参数分别为产品时间基准稳定度与PPP估计噪声,系数矩阵秩为2,能够求解出两个待估参数,得到产品时间基准稳定度。
其中,uk=(1 1 … 1)T为k×1的1值向量,zk=(0,0 ... 0)T为k×1的0值向量,uk-1=(1 1 … 1)T为(k-1)×1的1值向量,zk-1=(0,0 ... 0)T为(k-1)×1的0值向量;
表示测站r1的PPP接收机估计钟差序列的Allan方差,...,表示测站rk的PPP接收机估计钟差序列的Allan方差,表示测站r1与r2的时间同步结果序列的Allan方差,...,表示表示测站r1与rk的时间同步结果序列的Allan方差。
通过式八即可求解PPP产品的时间基准稳定度,同时也可得到PPP钟差估算噪声稳定度。
本发明强调的是一种基于PPP技术的PPP产品时间基准稳定度评估方法,也可用于PPP钟差估算噪声稳定度的评估。
实时卫星钟差改正数产品由不同分析中心计算并发布,由于不同分析中心对卫星钟差解算策略的不同,会导致产品间中时间基准稳定度的差异。进行GNSS授时接收机进行PPP单向授时,授时接收机输出时频信号的短期稳定度由接收机本地的晶振或原子钟决定,长期稳定度则取决于使用的实时卫星钟差产品的时间基准稳定度。对不同分析中心发布的实时卫星钟差产品时间基准稳定度进行评估,可以预测GNSS授时接收机输出时频信号的长期稳定度,选取时间基准稳定度更高的实时卫星钟差产品可以获取更稳定的PPP单向授时时频信号输出。
具体实施时,本发明技术方案提出的方法可由本领域技术人员采用计算机软件技术实现自动运行流程,实现方法的系统装置例如存储本发明技术方案相应计算机程序的计算机可读存储介质以及包括运行相应计算机程序的计算机设备,也应当在本发明的保护范围内。
在一些可能的实施例中,提供一种基于PPP的精密卫星星历时间基准稳定度评估系统,包括处理器和存储器,存储器用于存储程序指令,处理器用于调用存储器中的存储指令执行如上所述的一种基于PPP的精密卫星星历时间基准稳定度评估方法。
在一些可能的实施例中,提供一种基于PPP的精密卫星星历时间基准稳定度评估系统,包括可读存储介质,所述可读存储介质上存储有计算机程序,所述计算机程序执行时,实现如上所述的一种基于PPP的精密卫星星历时间基准稳定度评估方法。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (9)
1.一种基于PPP的精密卫星星历时间基准稳定度评估方法,其特征在于,包括以下步骤:
步骤S1,使用外接原子钟组频率源的IGS观测站的观测数据与不同分析中心的实时产品进行事后PPP处理,得到不同分析中心产品对应的观测站单站授时结果;选取一个测站作为中心测站,进行站间时间同步,得到时间同步结果;
步骤S2,计算单站授时与站间时间同步结果的重叠Allan方差,将Allan方差结果线性组合,分离产品基准稳定度与PPP估计噪声稳定度,实现方式包括以下子步骤,
步骤S2.1,计算钟差序列的Allan方差,包括将此Allan方差看作产品时间基准的稳定度与PPP估计噪声的线性组合;
步骤S2.2,计算时间同步结果序列的Allan方差,包括将每个观测站的PPP估计噪声视为具有相同噪声水平系数与成分,则将此Allan方差看作是系数为2的PPP估计噪声Allan方差的线性组合;
步骤S2.3,将S2.1与S2.2分别得到的Allan方差所表征的线性组合组成线性多项式,其待估参数分别为产品时间基准稳定度与PPP估计噪声,求解出两个待估参数,得到产品时间基准稳定度。
2.根据权利要求1所述一种基于PPP的精密卫星星历时间基准稳定度评估方法,其特征在于:步骤S1的实现过程包括以下子步骤,
步骤S1.1,选择IGS中GNSS接收机外接原子钟组的观测站,得到观测数据;
步骤S1.2,存储不同分析中心的改正数产品,组合广播星历与改正数进行精密星历恢复,得到精密星历;
步骤S1.3,使用观测数据与恢复后的精密星历进行PPP事后处理,得到的钟差序列看作包括产品时间基准、接收机钟差和PPP估计噪声的线性组合;
步骤S1.4,以S1.1中选择的观测站中的某一个站作为中心站,建立时间链路,进行站间时间同步,计算得到站间时间同步结果。
5.根据权利要求4所述一种基于PPP的精密卫星星历时间基准稳定度评估方法,其特征在于:步骤S2.1实现方式如下,
S(x)=xi+2m-2xi+m+xi (5)
其中,τ=mτ0表示平均时间,m表示平均因子,τ0表示测量间隔,N表示观测值数量,x表示N个以τ0为间隔的相位观测值;xi表示表示测量时间为i的相位数据,xi+2m、xi+m分别表示在xi之后,与xi时间间隔为m与2m的相位数据;S(x)表示相位数据的三差算子;表示平均时间为τ,计算是使用的相位数据为x的Allan方差函数;
7.根据权利要求6所述一种基于PPP的精密卫星星历时间基准稳定度评估方法,其特征在于:步骤S2.3实现方式如下,
其中,uk=(1 1…1)T为k×1的1值向量,zk=(0,0…0)T为k×1的0值向量,uk-1=(1 1…1)T为(k-1)×1的1值向量,zk-1=(0,0…0)T为(k-1)×1的0值向量;
表示测站r1的PPP接收机估计钟差序列的Allan方差,…,表示测站rk的PPP接收机估计钟差序列的Allan方差,表示测站r1与r2的时间同步结果序列的Allan方差,…,表示表示测站r1与rk的时间同步结果序列的Allan方差
通过式8,求解PPP产品的时间基准稳定度,同时得到PPP钟差估算噪声稳定度。
8.根据权利要求1或2或3或4或5或6或7所述一种基于PPP的精密卫星星历时间基准稳定度评估方法,其特征在于:用于GNSS授时接收机运行时,根据产品时间基准稳定度选择实时PPP产品,从而输出长期稳定度更好的时频信号。
9.一种基于PPP的精密卫星星历时间基准稳定度评估系统,其特征在于:用于实现如权利要求1-8任一项所述的一种基于PPP的精密卫星星历时间基准稳定度评估方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210599663.4A CN114755702B (zh) | 2022-05-30 | 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210599663.4A CN114755702B (zh) | 2022-05-30 | 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114755702A true CN114755702A (zh) | 2022-07-15 |
CN114755702B CN114755702B (zh) | 2024-08-02 |
Family
ID=
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115453582A (zh) * | 2022-08-17 | 2022-12-09 | 武汉大学 | 基于钟模型的卫星钟差估计方法及系统 |
CN116224746A (zh) * | 2023-01-29 | 2023-06-06 | 北京航空航天大学 | 一种星地原子钟融合的高稳定时间基准建立方法 |
CN117970767A (zh) * | 2024-03-29 | 2024-05-03 | 中国科学院国家授时中心 | 一种基于空间站高性能频率源的低轨星载钟噪声确定方法 |
CN115453582B (zh) * | 2022-08-17 | 2024-07-30 | 武汉大学 | 基于钟模型的卫星钟差估计方法及系统 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115453582A (zh) * | 2022-08-17 | 2022-12-09 | 武汉大学 | 基于钟模型的卫星钟差估计方法及系统 |
CN115453582B (zh) * | 2022-08-17 | 2024-07-30 | 武汉大学 | 基于钟模型的卫星钟差估计方法及系统 |
CN116224746A (zh) * | 2023-01-29 | 2023-06-06 | 北京航空航天大学 | 一种星地原子钟融合的高稳定时间基准建立方法 |
CN117970767A (zh) * | 2024-03-29 | 2024-05-03 | 中国科学院国家授时中心 | 一种基于空间站高性能频率源的低轨星载钟噪声确定方法 |
CN117970767B (zh) * | 2024-03-29 | 2024-06-11 | 中国科学院国家授时中心 | 一种基于空间站高性能频率源的低轨星载钟噪声确定方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111045034B (zh) | 基于广播星历的gnss多系统实时精密时间传递方法及系统 | |
CN110058287B (zh) | 一种低轨卫星定轨方法、装置及系统 | |
US7498984B2 (en) | Positioning system, information supply device, terminal device, control method of information supply device, control program of information supply device, and computer readable recording medium recording control program of information supply device | |
CN107728171B (zh) | 基于粒子滤波的gnss相位系统间偏差实时追踪和精密估计方法 | |
JP2017151113A (ja) | エフェメリス拡張システムとgnssでの使用方法 | |
JP2005524073A (ja) | ハイブリッド位置決定システムにおける無線端末のための位置決定 | |
US8358241B2 (en) | Position calculating method and position calculating device | |
CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及系统 | |
CN115993617B (zh) | 一种gnss系统时间偏差的监测方法 | |
CN111650614A (zh) | 卫星定位方法、芯片、模组和电子设备 | |
CN114721012B (zh) | 一种电离层延迟效应改正方法及装置 | |
CN114002712B (zh) | 基于数据自适应选择的gnss实时卫星钟差确定系统及方法 | |
CN116449400B (zh) | 一种北斗三号ppp服务实时卫星钟差评价方法及系统 | |
CN116224746A (zh) | 一种星地原子钟融合的高稳定时间基准建立方法 | |
CN114755702A (zh) | 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 | |
RU2486544C2 (ru) | Способ определения координат потребителя в системах глонасс/gps и устройство для его реализации | |
CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
CN114755702B (zh) | 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 | |
CN112946693B (zh) | 一种卫星导航系统的系统时间偏差的确定方法及系统 | |
CN115603847A (zh) | 基于北斗短报文和PPP-B2b服务的亚纳秒级时间同步方法及装置 | |
Piriz et al. | Performance of the NeQuick G iono model for single-frequency GNSS timing applications | |
EP3748401A1 (en) | System and method for position determination of a stationary gnss receiver using a distributed time signal | |
CN118033696B (zh) | 实时多gnss时间尺度保持方法及系统 | |
CN116893433B (zh) | 一种实现跟踪站观测值预测的方法及装置 | |
Tseng et al. | Introduction of GNSS Up-Sampled All-in-View Time Transfers |
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 |