CN112946693B - 一种卫星导航***的***时间偏差的确定方法及*** - Google Patents
一种卫星导航***的***时间偏差的确定方法及*** Download PDFInfo
- Publication number
- CN112946693B CN112946693B CN202110149454.5A CN202110149454A CN112946693B CN 112946693 B CN112946693 B CN 112946693B CN 202110149454 A CN202110149454 A CN 202110149454A CN 112946693 B CN112946693 B CN 112946693B
- Authority
- CN
- China
- Prior art keywords
- satellite navigation
- clock
- navigation system
- deviation
- difference
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000005070 sampling Methods 0.000 claims description 57
- 238000012937 correction Methods 0.000 claims description 39
- 230000009897 systematic effect Effects 0.000 claims description 22
- 238000009499 grossing Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 3
- 239000005433 ionosphere Substances 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 239000005436 troposphere Substances 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
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/20—Integrity monitoring, fault detection or fault isolation of space segment
-
- 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
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Computer Security & Cryptography (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Electric Clocks (AREA)
Abstract
本发明涉及一种卫星导航***的***时间偏差的确定方法。所述方法包括:获取设定时间段内各卫星导航***的广播钟差参数和精密钟差;基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差;基于所述广播钟差和对应的精密钟差计算对应卫星的钟差之差;根据所述卫星导航***内所有卫星的钟差之差计算***偏差;对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差。本发明能提高***时间偏差测量的可及性、便捷性和精度。
Description
技术领域
本发明涉及卫星导航领域,特别是涉及一种卫星导航***的***时间偏差的确定方法及***。
背景技术
卫星导航***已经成为最重要的定位和授时基础设施。目前主要的卫星导航***包括美国的GPS、俄罗斯的GLONASS、欧盟的Galileo和中国的BDS。卫星导航***的接收机是通过测量电磁波从发射到被接收所经历的时间实现定位和授时的,因此,每一个卫星导航***维持着专门的***时间。例如,GPS的***时间由主控站、监测站和卫星的原子钟组共同维持和实现,并溯源到美国海军天文台维持的协调世界时UTC。BDS的***时间由位于主控站的原子钟维持和实现,并溯源到中国科学院国家授时中心维持的协调世界时UTC。GPS的***时间和BDS的***时间都是连续的原子时***,根据定义,前者的时间起点是UTC1980年1月6日0时0分0秒,后者的时间起点是UTC 2006年1月1日0时0分0秒,即GPS的***时间和BDS的***时间分别在这两个时刻与UTC对齐。已知的,UTC存在跳秒问题,从1980年1月6日0时0分0秒到2006年1月1日0时0分0秒期间,UTC发生了14次跳秒,导致名义上GPS的***时间比BDS的***时间快了14秒。由于GPS的***时间和BDS的***时间都是通过具体的原子钟维持和实现的,原子钟输出的时间信号总是或多或少的存在误差,尽管原子钟已经是目前人类所能获得的最高精度的守时装置。实际的GPS的***时间与BDS的***时间之间的偏差(不含名义偏差)就是GPS和BDS的***时间偏差,其量级可能达到十几、甚至几十纳秒。如果不考虑不同卫星导航***的***时间之间存在的偏差,直接将不同卫星导航***的测量数据联合用于定位或授时,将带来额外的误差。为了使用户在定位授时过程中可以有更多的选择,即可以将属于不同卫星导航***的卫星联合起来(如2颗GPS卫星+2颗BDS卫星)用于定位和授时,卫星导航***的主控站需要测量出不同导航卫星***的***时间偏差并作为改正参数,同轨道和钟差参数一起播发给用户接收机。例如,BDS***的主控站会事先测量出BDS的***时间相对于GPS的***时间的偏差,编入导航电文(如参数:A0GPS、A1GPS)发播给用户,以便用户能够联合BDS和GPS卫星进行定位、测速和授时。
已知的,测量不同卫星导航***的***时间偏差的方法有①基于用户端接收机秒脉冲信号的测量方法:用不同的接收机接收不同卫星导航***的导航信号并输出1PPS时间信号,通过精密时间间隔计数器对不同卫星导航***接收机输出的1PPS时间信号进行直接或间接比对。②直接比对法:采用双向时间频率传递、卫星共视比对链路对各卫星导航***的时间基准站的时间进行直接或间接的比对。③接收机钟差比对法:将多卫星导航***接收机架设在具有已知坐标的站点上,分别用不同的卫星导航***数据进行接收机钟差解算,将不同卫星导航***数据解算得到的接收机钟差进行求差就得到不同卫星导航***的***时间之差。
已知的测量不同卫星导航***的***时间偏差的方法具有明显的不足。所述方法①和②需要专门的设备才能进行卫星导航******时间的比对,如方法①要求接收机能够输出1PPS秒脉冲信号,且需要有专门的时间间隔计数器才能进行***时间偏差测量;方法②需要有专门的双向时间频率传递或卫星共视设备和链路。方法③需要有多卫星导航***接收机,且受用户接收机端的观测误差、多路径效应、对流程延迟、电离层延迟等环境因素影响较大,监测精度不高。
发明内容
本发明的目的是提供一种卫星导航***的***时间偏差的确定方法及***,提高了***时间偏差确定的可及性、便捷性和精度。
为实现上述目的,本发明提供了如下方案:
一种卫星导航***的***时间偏差的确定方法,包括:
获取设定时间段内各卫星导航***的钟差数据;所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和在同一时间基准下的所述卫星导航***内所有卫星的精密钟差;
基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差;
基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差;
对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差;
对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差。
可选的,所述基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差,具体为:
根据公式计算卫星的广播钟差;其中,/>为第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。
可选的,所述基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差,具体为:
根据公式计算卫星的钟差之差;其中,/>表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,/>表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项。
可选的,所述根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差,具体为:
根据公式计算卫星导航***g的***偏差;其中,/>表示第i个采样点对应的卫星导航***g的***偏差,MED{}表示求中位数,表示第i个采样点对应的卫星导航***g内的卫星s的钟差之差。
可选的,所述根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差,具体为:
根据公式计算卫星导航***g和卫星导航***f之间的***时间偏差;其中,/>表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,/>表示第i个采样点对应的卫星导航***g的***偏差,/>表示第i个采样点对应的卫星导航***f的***偏差。
可选的,在所述对于任意两个所述卫星导航***,根据相应的两个所述***偏差计算两个所述卫星导航***之间的***时间偏差之后,还包括:
对所述***时间偏差进行平滑,得到平滑后的***时间偏差值;
采用最小二乘法对所述平滑后的***时间偏差进行多项式拟合,得到多项式系数;
利用所述多项式系数预测未来某一时刻的***时间偏差。
一种卫星导航***的***时间偏差的确定***,包括:
获取模块,用于获取设定时间段内各卫星导航***的钟差数据;所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和在同一时间基准下的所述卫星导航***内所有卫星的精密钟差;
广播钟差确定模块,用于基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差;
钟差之差确定模块,用于基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差;
***偏差确定模块,用于对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差;
***时间偏差确定模块,用于对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差。
可选的,所述广播钟差确定模块,具体包括:
广播钟差确定单元,用于根据公式
计算卫星的广播钟差;其中,/>为第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。
可选的,所述钟差之差确定模块,具体包括:
钟差之差确定单元,用于根据公式计算卫星的钟差之差;其中,/>表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,/>表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项。
可选的,所述***时间偏差确定模块,具体包括:
***时间偏差确定单元,用于根据公式计算卫星导航***g和卫星导航***f之间的***时间偏差;其中,/>表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,/>表示第i个采样点对应的卫星导航***g的***偏差,表示第i个采样点对应的卫星导航***f的***偏差。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明在计算***时间偏差时考虑到了***中所有卫星的钟差之差,而且不需要专门的时间比对或时间同步设备,也不需要处理接收机数据,不受测站环境、接收机噪声、多路径效应、对流层和电离层延迟误差的影响,提高了***时间偏差确定的可及性、便捷性和精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例1提供的一种卫星导航***的***时间偏差的确定方法流程图;
图2为本发明实施例2提供的一种卫星导航***的***时间偏差的确定***的结果示意图;
图3为本发明实施例1步骤4得到的GPS***时间和GLONASS***时间相对于精密钟差基准时间的***偏差序列示意图;
图4为本发明实施例1步骤5得到的GPS和GLONASS***时间偏差及步骤6平滑后的GPS和GLONASS***时间偏差示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示,本实施例提供了一种卫星导航***的***时间偏差的确定方法,所述方法包括:
101:获取设定时间段内各卫星导航***的钟差数据。所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和精密钟差,所述精密钟差是指在同一时间基准下得到的包含不同卫星导航***内卫星的精密钟差的数据;其中,广播钟差参数,以(二次)多项式系数的形式给出,还包括对应的参考时间参数。例如国际卫星导航***服务组织(IGS)多卫星导航***实验项目(MGEX)提供的同一时间基准的精密钟差产品或全球连续监测评估***(iGMAS)提供的多卫星导航***精密钟差产品。广播钟差参数可以通过多卫星导航***接收机获得,也可以从公开的卫星导航***数据中心,如IGS数据中心获得,典型地,广播钟差参数以RINEX格式存放。
102:基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差。
103:基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差。
104:对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差。
105:对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差。
在实际应用中,102具体为:
用广播钟差参数计算对应于精密钟差采样时刻的卫星导航***的广播钟差,即根据公式
计算卫星的广播钟差;其中,为第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。假设5分钟计算一个广播钟差,则在1小时内就会有12个采样点。
在实际应用中,103具体为:
计算每颗卫星的广播钟差与精密钟差之差;需要考虑天线相位中心偏差改正,码间偏差改正和相对论效应改正,即根据公式
计算卫星的钟差之差;其中,表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,/>表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项。
其中天线相位中心偏差改正项是由于广播钟差和精密钟差解算时采用的天线相位偏差(PCO,Phase CenterOffset)参数的不同引起的,主要考虑卫星体Z方向PCO差异的影响,公式为:
Δpco=Zbrd-Zpre (3)
其中Zbrd和Zpre分别表示广播钟差和精密钟差采用的Z方向PCO值。
码间偏差改正项Δdcb是由于广播钟差和精密钟差解算时采用的测距码的不同引起。例如对于GPS,广播钟差和精密钟差都采用P1和P2码的消电离层组合,因此此项改正为0;BDS民用信号的广播钟差对应于B3I,而精密钟差通常采用B1I和B3I或B2a和B1C的消电离层组合,需要采用已知的方法计算该项改正。
对于GPS、Galileo和BDS的广播钟差不包含相对论效应,GLONASS的广播钟差则包含有相对论效应,而精密钟差产品,如IGS MGEX项目提供的多卫星导航***精密钟差不包含相对论效应,因此对于GPS、Galileo和BDS***,相对论该正项Δrel为0;对于GLONASS,相对论改正项需要采用已知的方法计算。
在实际应用中,104具体为:
根据公式
计算卫星导航***g的***偏差;其中,表示第i个采样点对应的卫星导航***g的***偏差,MED{}表示求中位数,/>表示第i个采样点对应的卫星导航***g内的卫星s的钟差之差。
在实际应用中,104具体还可以为:
根据公式
计算卫星导航***g的***偏差;其中,表示第i个采样点对应的卫星导航***g的***偏差,MEAN{}表示求平均值,/>表示第i个采样点对应的卫星导航***g的卫星s的钟差之差。
在上述步骤中,采用平均值或中位值,是因为全星座所有卫星的钟差才能代表***时间,其中中位值比平均值具有更好的抗差性能。
在实际应用中,计算两个卫星导航***之间的***时间偏差,采用如下公式:
其中,表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,/>表示第i个采样点对应的卫星导航***g的***偏差,/>表示第i个采样点对应的卫星导航***f的***偏差。
在实际应用中,在所述对于任意两个所述卫星导航***,根据相应的两个所述***偏差计算两个所述卫星导航***之间的***时间偏差之后,还包括:
对所述***时间偏差进行平滑,得到平滑后的***时间偏差值,具体为对序列,采用已知的方法(如Vondrak方法)进行平滑,得到平滑后的***时间偏差值/>采用Vondrak的优势是该方法不需要事先知道序列的变化规律,可以通过调节参数兼顾平滑结果的符合度和平滑度,是经典的数据平滑方法。
对***时间偏差进行预报;采用最小二乘法对所述平滑后的***时间偏差进行多项式拟合,得到多项式系数;利用所述多项式系数预测未来某一时刻的***时间偏差,多项式可以为一次多项式。
本实施例还采用一个具体的实施案例对上述方法的技术方案及效果进行详细阐述。该实施例中,确定2020年10月01日GLONASS和GPS的***时间偏差。
步骤1、收集数据:收集2020年10月01日的卫星导航***广播钟差参数和多卫星导航***精密钟差。从IGS网站下载GPS和GLONASS RINEX格式的导航文件和IGS MGEX分析中心GBM提供的多***精密钟差文件。
步骤2、采用公式(1)用广播钟差参数计算对应于精密钟差采样时刻的GPS和GLONASS的广播钟差序列,由于精密钟差的时间基准名义上对齐到GPS时,所以对于GPS,ΔTGPS=0;对于GLONASS***ΔTGLO=-18s。
步骤3、采用公式(2)计算钟差之差,其中Zbrd的值采用卫星厂家提供的值,Zpre由IGS网站ANTEX文件提供。
码间偏差改正项Δdcb是由于广播钟差和精密钟差解算采用的测距码的不同引起。对于GPS和GLONASS,此项改正都为0。对于GPS卫星,Δrel为0;对于GLONASS卫星,其中r和v分别为卫星的位置和速度向量,c为光在真空中的传播速度。
步骤4、如图3所示,采用公式(4)计算***偏差,得到GPS***时间与精密钟差基准时间的***偏差以及GLONASS***时间与精密钟差基准时间的***偏差/>对于GPS***nGPS=32,对于GLONASS***nGLO=24。
步骤5、计算GPS和GLONASS之间的***时间偏差。
采用如下公式:
其中,表示Ti时刻GPS和GLONASS的***时间偏差,如图4中的正十字图例所示。
步骤6、对GPS和GLONASS之间的***时间偏差进行平滑;对序列,采用已知的方法(如Vondrak方法)进行平滑,得到平滑后的***时间偏差值/>结果如图4中的斜十字所示。
实施例2
本实施例提供了一种与实施例1对应的卫星导航***的***时间偏差的确定***,如图2所示,所述***包括:
获取模块A1,用于获取设定时间段内各卫星导航***的钟差数据;所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和在同一时间基准下的所述卫星导航***内所有卫星的精密钟差。
广播钟差确定模块A2,用于基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差。
钟差之差确定模块A3,用于基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差。
***偏差确定模块A4,用于对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差。
***时间偏差确定模块A5,用于对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差。
作为一种可选的实施方式,所述广播钟差确定模块,具体包括:
广播钟差确定单元,用于根据公式
计算卫星的广播钟差;其中,/>为第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。
作为一种可选的实施方式,所述钟差之差确定模块,具体包括:
钟差之差确定单元,用于根据公式计算卫星的钟差之差;其中,/>表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,/>表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项。
作为一种可选的实施方式,所述***偏差确定模块,具体包括:
***偏差确定单元,用于根据公式计算卫星导航***g的***偏差;其中,/>表示第i个采样点对应的卫星导航***g的***偏差,MED{}表示求中位数,/>表示第i个采样点对应的卫星导航***g内的卫星s的钟差之差。
作为一种可选的实施方式,所述***时间偏差确定模块,具体包括:***时间偏差确定单元,用于根据公式计算卫星导航***g和卫星导航***f之间的***时间偏差;其中,/>表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,/>表示第i个采样点对应的卫星导航***g的***偏差,/>表示第i个采样点对应的卫星导航***f的***偏差。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
1)不需要专门的比对设备就可以对不同卫星导航***的***时间偏差进行测量。
2)在测量***时间偏差时考虑到星座中所有卫星的钟差差异的中位值或平均值,理论上更准确,具有更好的稳健性,因为***时间是由全星座所有卫星共同代表,而不是某一颗或者某一部分卫星代表。
3)不需要处理接收机数据,不受测站环境,接收机噪声,多路径效应,对流层和电离层延迟误差的影响。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种卫星导航***的***时间偏差的确定方法,其特征在于,包括:
获取设定时间段内各卫星导航***的钟差数据;所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和在同一时间基准下的所述卫星导航***内所有卫星的精密钟差;
基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差;
基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差;
对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差,具体为:
根据公式计算卫星导航***g的***偏差;其中,表示第i个采样点对应的卫星导航***g的***偏差,MED{}表示求中位数,表示第i个采样点对应的卫星导航***g内的卫星s的钟差之差;
对于任意两个所述卫星导航***,根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差,所述基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差,具体为:
根据公式计算卫星的钟差之差;其中,表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,/>表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项。
2.根据权利要求1所述的卫星导航***的***时间偏差的确定方法,其特征在于,所述基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差,具体为:
根据公式计算卫星的广播钟差;其中,/>为第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。
3.根据权利要求1所述的卫星导航***的***时间偏差的确定方法,其特征在于,所述根据相应的两个所述卫星导航***的***偏差计算两个所述卫星导航***之间的***时间偏差,具体为:
根据公式计算卫星导航***g和卫星导航***f之间的***时间偏差;其中,/>表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,/>表示第i个采样点对应的卫星导航***g的***偏差,/>表示第i个采样点对应的卫星导航***f的***偏差。
4.根据权利要求1所述的卫星导航***的***时间偏差的确定方法,其特征在于,在所述对于任意两个所述卫星导航***,根据相应的两个所述***偏差计算两个所述卫星导航***之间的***时间偏差之后,还包括:
对所述***时间偏差进行平滑,得到平滑后的***时间偏差值;
采用最小二乘法对所述平滑后的***时间偏差进行多项式拟合,得到多项式系数;
利用所述多项式系数预测未来某一时刻的***时间偏差。
5.一种卫星导航***的***时间偏差的确定***,其特征在于,包括:
获取模块,用于获取设定时间段内各卫星导航***的钟差数据;所述钟差数据包括所述卫星导航***内所有卫星的广播钟差参数和在同一时间基准下的所述卫星导航***内所有卫星的精密钟差;
广播钟差确定模块,用于基于所述广播钟差参数计算所述卫星导航***内各所述卫星的广播钟差;
钟差之差确定模块,用于基于所述广播钟差和对应卫星的精密钟差计算对应卫星的钟差之差;
***偏差确定模块,用于对于任意一个卫星导航***,根据所述卫星导航***内所有卫星的钟差之差计算所述卫星导航***的***偏差,具体为:根据公式
计算卫星导航***g的***偏差;其中,/>表示第i个采样点对应的卫星导航***g的***偏差,MED{}表示求中位数,/>表示第i个采样点对应的卫星导航***g内的卫星s的钟差之差;
***时间偏差确定模块,用于对于任意两个所述卫星导航***,根据相应的两个所述卫星导,所述钟差之差确定模块,具体包括:
钟差之差确定单元,用于根据公式计算卫星的钟差之差;其中,/>表示第i个采样点对应的卫星导航***g内卫星s的钟差之差,表示第i个采样点对应的卫星导航***g内卫星s的广播钟差,/>表示第i个采样点对应的卫星导航***g内卫星s的精密钟差,Δpco表示天线相位中心偏差改正项,Δdcb表示码间偏差改正项,Δrel表示相对论改正项航***的***偏差计算两个所述卫星导航***之间的***时间偏差。
6.根据权利要求5所述的卫星导航***的***时间偏差的确定***,其特征在于,所述广播钟差确定模块,具体包括:
广播钟差确定单元,用于根据公式
星的广播钟差;其中,为第i个采样点对应的卫星导航***g内卫星s的广播钟差,为广播钟差参数的0次项系数,/>为广播钟差参数的一次项系数,/>为广播钟差参数的二次项系数,/>为卫星导航***g内卫星s的广播钟差参数的参考时刻;Ti为第i个采样点对应的时间,ΔTg为卫星导航***g的***时间与精密钟差所用的时间基准的名义偏差。
7.根据权利要求5所述的卫星导航***的***时间偏差的确定***,其特征在于,所述***时间偏差确定模块,具体包括:
***时间偏差确定单元,用于根据公式计算卫星导航***g和卫星导航***f之间的***时间偏差;其中,/>表示第i个采样点对应的卫星导航***g和卫星导航***f的***时间偏差,Bi g表示第i个采样点对应的卫星导航***g的***偏差,/>表示第i个采样点对应的卫星导航***f的***偏差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110149454.5A CN112946693B (zh) | 2021-02-03 | 2021-02-03 | 一种卫星导航***的***时间偏差的确定方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110149454.5A CN112946693B (zh) | 2021-02-03 | 2021-02-03 | 一种卫星导航***的***时间偏差的确定方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112946693A CN112946693A (zh) | 2021-06-11 |
CN112946693B true CN112946693B (zh) | 2024-01-23 |
Family
ID=76242328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110149454.5A Active CN112946693B (zh) | 2021-02-03 | 2021-02-03 | 一种卫星导航***的***时间偏差的确定方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112946693B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115993617B (zh) * | 2023-02-13 | 2024-01-26 | 中国科学院国家授时中心 | 一种gnss***时间偏差的监测方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702030A (zh) * | 2009-10-29 | 2010-05-05 | 中国科学院国家授时中心 | 一种站间钟差辅助的导航卫星精密定轨方法 |
CN104297764A (zh) * | 2013-07-16 | 2015-01-21 | 凹凸电子(武汉)有限公司 | 一种提高导航***时间的pps准确度的方法和接收机 |
CN105738927A (zh) * | 2016-02-01 | 2016-07-06 | 中国人民解放军61540部队 | 顾及卫星导航***间偏差的融合定位导航授时计算模型构建方法 |
CN108241158A (zh) * | 2017-12-25 | 2018-07-03 | 航天恒星科技有限公司 | 一种gnss卫星异常状态快速监测技术 |
CN109061696A (zh) * | 2018-09-28 | 2018-12-21 | 中国人民解放军61540部队 | 一种确定导航卫星轨道和钟差的方法 |
CN110488328A (zh) * | 2019-07-18 | 2019-11-22 | 北京未来导航科技有限公司 | 低轨卫星导航增强平台的电文收发方法及*** |
CN110851773A (zh) * | 2019-10-29 | 2020-02-28 | 中海北斗(深圳)导航技术有限公司 | Gnss实时钟差评估算法 |
CN110988951A (zh) * | 2019-12-06 | 2020-04-10 | 中国地质大学(北京) | 多源数据融合实时导航定位方法及*** |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
WO2020228754A1 (zh) * | 2019-05-16 | 2020-11-19 | 北京合众思壮科技股份有限公司 | 一种低轨卫星定轨方法、装置及*** |
-
2021
- 2021-02-03 CN CN202110149454.5A patent/CN112946693B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702030A (zh) * | 2009-10-29 | 2010-05-05 | 中国科学院国家授时中心 | 一种站间钟差辅助的导航卫星精密定轨方法 |
CN104297764A (zh) * | 2013-07-16 | 2015-01-21 | 凹凸电子(武汉)有限公司 | 一种提高导航***时间的pps准确度的方法和接收机 |
CN105738927A (zh) * | 2016-02-01 | 2016-07-06 | 中国人民解放军61540部队 | 顾及卫星导航***间偏差的融合定位导航授时计算模型构建方法 |
CN108241158A (zh) * | 2017-12-25 | 2018-07-03 | 航天恒星科技有限公司 | 一种gnss卫星异常状态快速监测技术 |
CN109061696A (zh) * | 2018-09-28 | 2018-12-21 | 中国人民解放军61540部队 | 一种确定导航卫星轨道和钟差的方法 |
WO2020228754A1 (zh) * | 2019-05-16 | 2020-11-19 | 北京合众思壮科技股份有限公司 | 一种低轨卫星定轨方法、装置及*** |
CN110488328A (zh) * | 2019-07-18 | 2019-11-22 | 北京未来导航科技有限公司 | 低轨卫星导航增强平台的电文收发方法及*** |
CN110851773A (zh) * | 2019-10-29 | 2020-02-28 | 中海北斗(深圳)导航技术有限公司 | Gnss实时钟差评估算法 |
CN111025346A (zh) * | 2019-11-18 | 2020-04-17 | 广州南方卫星导航仪器有限公司 | 一种快速估计gnss精密卫星钟差的方法及存储介质 |
CN110988951A (zh) * | 2019-12-06 | 2020-04-10 | 中国地质大学(北京) | 多源数据融合实时导航定位方法及*** |
Non-Patent Citations (17)
Title |
---|
Atomic Clock Performance Assessment of BeiDou-3 Basic System with the Noise Analysis of Orbit Determination and Time Synchronization;Xiaolin Jia;《remote sensing》;全文 * |
BDS卫星钟差半参数平差模型异常数据探测与处理;李文涛;颜雄;夏磊;任青阳;陈洁;潘雄;;测绘学报(第01期);全文 * |
GNSS开放服务的监测评估;焦文海;丁群;李建文;卢晓春;冯来平;马加庆;陈罡;;中国科学:物理学 力学 天文学(第05期);全文 * |
GPS/Galileo组合定位和Galileo卫星钟评估;于碧云;张慧君;李孝辉;;时间频率学报(第02期);全文 * |
Time Synchronization Algorithm for Wireless Sensor Networks Based on Frequency Offset Estimation;Shaojun Yu;《iJOE》;全文 * |
利用SLR及伪距数据进行卫星预报钟差精度验证;毛悦;贾小林;秦显平;;武汉大学学报(信息科学版)(第12期);全文 * |
北斗三号基本***空间信号性能分析;毛悦;《测绘科学技术学报》;全文 * |
北斗三号空间信号测距误差评估与对比分析;***;《测绘学报》;1213-1221 * |
北斗卫星导航***误差分析与评估;翟显;《现代导航》;10-16 * |
周江文.《抗差最小二乘法》.武汉:华中理工大学出版社,1997,11. * |
基于GLONASS频间偏差的GNSS时差监测方法研究;孔思嘉;谢维华;;全球定位***(第04期);全文 * |
基于非同时双向星间链路的自主时间同步仿真分析;方琳;杨旭海;孙保琪;钦伟瑾;孔垚;;天文学报(第05期);全文 * |
星间单程伪距联合监测站数据确定北斗三号卫星轨道和钟差;阮仁桂;《测绘学报》;全文 * |
***时间偏差应用于GNSS位置解算及完好性改善研究;王欢;《中国优秀硕士学位论文全文数据库 信息科技辑》;6-8 * |
***时间偏差应用于GNSS位置解算及完好性改善研究;王欢;《中国优秀硕士学位论文全文数据库信息科技辑》;6-9 * |
翟显.北斗卫星导航***误差分析与评估.《现代导航》.2018,10-16. * |
顾及不同天线相位中心改正模型的北斗空间信号精度评估方法;张柔;胡志刚;陶钧;王晨;赵齐乐;;武汉大学学报(信息科学版)(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112946693A (zh) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110031881B (zh) | 高精度星间激光测距辅助精密单点定位的方法 | |
WO2020228754A1 (zh) | 一种低轨卫星定轨方法、装置及*** | |
Lombardi et al. | Time and frequency measurements using the global positioning system | |
CN109001771B (zh) | 导航卫星和低轨卫星实时钟差确定及预报方法和*** | |
CN102998681B (zh) | 一种卫星导航***的高频钟差估计方法 | |
JP2005521863A (ja) | 衛星測位システム受信機での時間の決定とそのための方法 | |
CN114994727B (zh) | 一种实现高精度时间校准及卫星定位的设备 | |
CN112543078A (zh) | 一种基于卫星共视的网络时间服务器校准方法与装置 | |
CN117388881B (zh) | 一种低轨卫星的星载原子钟向UTC(k)的溯源方法及*** | |
CN113325446B (zh) | 一种多模共频gnss载波相位时间传递方法及*** | |
CN115993617B (zh) | 一种gnss***时间偏差的监测方法 | |
CN110082801B (zh) | 一种芯片原子钟辅助bds/gps接收机的定位方法 | |
CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测***及方法 | |
CN112946693B (zh) | 一种卫星导航***的***时间偏差的确定方法及*** | |
CN111970002A (zh) | 基于北斗geo卫星的原子钟远程频率传递与比对方法 | |
WO2022001368A1 (zh) | 雷达及雷达*** | |
GB2459334A (en) | Reducing the time required to establish a position solution in a satellite positioning system | |
CN210742507U (zh) | 一种基于全球导航卫星***的标准时间频率源装置 | |
CN116224746A (zh) | 一种星地原子钟融合的高稳定时间基准建立方法 | |
EP3748401A1 (en) | System and method for position determination of a stationary gnss receiver using a distributed time signal | |
Zhang et al. | Study on PPP time comparison based on BeiDou-3 new signal | |
Gogoi et al. | NavIC receiver clock offsets estimation with common view master clock method | |
CN118091718B (zh) | 通过低轨卫星下行导航信号改善ut1解算精度的方法 | |
CN110764123A (zh) | 一种基于glonass广播星历的伪距定位改进方法 | |
CN114755702B (zh) | 基于ppp的精密卫星星历时间基准稳定度评估方法及*** |
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 |