CN109655846A - 一种多站差分后处理高精度时间同步方法及系统 - Google Patents
一种多站差分后处理高精度时间同步方法及系统 Download PDFInfo
- Publication number
- CN109655846A CN109655846A CN201910091296.5A CN201910091296A CN109655846A CN 109655846 A CN109655846 A CN 109655846A CN 201910091296 A CN201910091296 A CN 201910091296A CN 109655846 A CN109655846 A CN 109655846A
- Authority
- CN
- China
- Prior art keywords
- satellite navigation
- satellite
- common
- receiver
- navigation receiver
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000012805 post-processing Methods 0.000 title claims abstract description 18
- 230000005540 biological transmission Effects 0.000 claims abstract description 29
- 238000012545 processing Methods 0.000 claims abstract description 23
- 238000012937 correction Methods 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 8
- 230000001960 triggered effect Effects 0.000 claims description 7
- 239000005433 ionosphere Substances 0.000 claims description 5
- 230000006854 communication Effects 0.000 claims description 4
- 238000004891 communication Methods 0.000 claims description 4
- 230000001360 synchronised effect Effects 0.000 claims description 4
- 239000005436 troposphere Substances 0.000 claims description 4
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 239000000969 carrier Substances 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000004069 differentiation Effects 0.000 description 3
- 230000002457 bidirectional effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000007175 bidirectional communication Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000003068 static effect Effects 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/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
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
本发明涉及一种多站差分后处理高精度时间同步方法及系统,其包括安装在每个飞行器上的卫星导航接收机,卫星导航接收机接收卫星导航信号,产生PPS秒脉冲信号,产生时刻为TA0、TB0,将卫星导航接收机基带原始观测量信息以及TA0、TB0都传输至数据传输系统,由数据传输系统将接收到的基带原始观测量信息发送到地面接收机系统,飞行器之间不需要进行通信;地面接收机系统同时接收多个飞行器上的卫星导航接收机的原始观测量信息,并传输至数据处理系统,数据处理系统同时解算多个卫星导航接收机的原始观测量信息,得到共视卫星,再经过相对差分计算,修正每个接收机触发PPS的时标,实现多站时间高度同步。本发明能在卫星导航高精度授时领域中应用。
Description
技术领域
本发明涉及一种卫星导航高精度授时领域,特别是关于一种多站差分后处理高精度时间同步方法及系统。
背景技术
卫星导航技术在飞行器上使用非常普遍,提供实时位置和时间。由于卫星接收机各自有各自的时钟源,虽然都接收导航卫星的扩频信息得到同步时间,与世协调时(UTC)保持同步,但是受卫星自身时钟源抖动、空间电离层误差干扰、以及卫星导航接收机自身晶振漂移等多种因素影响,一般接收机触发的PPS精度相对于UTC时间在±30ns量级,卫星接收机两两之间的时间精度可能达到60ns,不满足高度时间同步要求的使用场合应用。
目前比较成熟的卫星共视技术可以使两地间时间高度统一,达到1ns级别,再结合双向比对技术可以更高精度。但是该技术实施的前提是两地卫星导航接收机是静态的,且预先精确测量卫星接收机天线位置,共视同步过程需要10分钟到30分钟之久,无法直接使用在飞行器场合上使用。
发明内容
针对上述多站之间无站间数据链时的高精度时间同步技术难题,本发明的目的是提供一种多站差分后处理高精度时间同步方法及系统,其能使相互间距离为20千米以内的飞行器,快速获取5ns高精度时间同步精度。
为实现上述目的,本发明采取以下技术方案:一种多站差分后处理高精度时间同步系统,其包括安装在每个飞行器上的卫星导航接收机、卫星导航天线、数据传输系统和数据传输系统天线,以及设置在地面接收站的地面接收机系统和数据处理系统;所述卫星导航接收机通过所述卫星导航天线接收卫星导航信号,产生PPS秒脉冲信号,产生时刻为TA0、TB0,将所述卫星导航接收机基带原始观测量信息以及TA0、TB0都传输至所述数据传输系统,由所述数据传输系统将接收到的基带原始观测量信息经所述数据传输系统天线发送所述到地面接收机系统,飞行器之间不需要进行通信;所述地面接收机系统同时接收多个飞行器上的所述卫星导航接收机的原始观测量信息,并传输至所述数据处理系统,所述数据处理系统同时解算多个所述卫星导航接收机的原始观测量信息,得到共视卫星,再经过相对差分计算,修正每个接收机触发PPS的时标为TA1,TB1,实现多站时间高度同步。
一种基于上述系统的多站差分后处理高精度时间同步方法,其包括以下步骤:1)卫星导航接收机接收导航卫星信号,产生PPS秒脉冲信号,卫星导航接收机实时计算出PPS的时标分别为TA0和TB0,并输出基带原始观测量信息,通过数据传输系统发送到地面数据接收系统;飞行器上的数据传输系统负责向地面接收系统发送卫星观测量信息;2)数据接收系统接收多个飞行器上的原始观测量信息,并由地面接收站的数据处理系统同时解算多个卫星接收机原始观测量信息,找出共视卫星,再经过相对差分计算,修正每个卫星导航接收机触发PPS的时标。
进一步,所述步骤1)中,PPS的时标TA0和TB0计算公式如下:
TA0=tAt+tAk+tAR+ΔtA,TB0=tBt+tBk+tBR+ΔtB
式中,tAt,tBt分别代表A与B卫星导航接收机自身时钟引入的误差;tAK,tBK分别代表A与B卫星导航接收机计算时标时由多颗卫星空间公有误差引入的误差量;tAR,tBR代表A与B卫星导航接收机触发PPS时刻的理想时标;ΔtA,ΔtB代表系统噪声。
进一步,所述步骤1)中,卫星导航接收机观测量误差包含接收机公有误差和自有误差,卫星导航接收机公有误差通过伪距差分消除,其方法如下:
卫星导航接收机距离基站在200km以内情况下,利用基站伪距改正数修正卫星导航接收机n的伪距测量值:消除公有误差,即:
式中,ΔD表示差分后的系统误差;(xn,yn,zn)表示卫星导航接收机n的坐标;表示伪距改正数;表示卫星导航接收机n接收到卫星i的伪距;卫星i的坐标为(xi,yi,zi);卫星导航接收机n的观测量误差经过上述修正,公有误差基本被消除。
进一步,所述步骤2)中,每个卫星导航接收机触发PPS的时标修正包括以下步骤:2.1)以其中任意一个卫星导航接收机为主站,其它卫星导航接收机为从站,然后对主站进行定位解算,获得主站当前时刻位置、速度,以及PPS触发时刻的时标,并得到该时刻每个卫星差分修正量信息;2.2)采用主站生成的差分和共视数据参数,处理其它从站卫星导航接收机数据,根据相对差分原理以及卫星共视授时基本原理,重新解算出各个从站差分定位结果,并重新计算出各从站PPS秒脉冲触发时刻的时标TA1与TB1,即:
TA1=tAt+tAk+tAR+Δt0,TB1=tBt+tAk+tBR+Δt1
式中,tAt,tBt分别代表A与B卫星导航接收机自身时钟引入的误差;tAK,tBK分别代表A与B卫星导航接收机计算时标时由多颗卫星空间公有误差引入的误差量;tAR,tBR代表A与B卫星导航接收机触发PPS时刻的理想时标;ΔtA,ΔtB代表系统噪声;Δt0和Δt1表示差分和共视处理后卫星导航接收机A与B的系统噪声;卫星接收机两个卫星导航接收机在短基线情况下空间公有误差接近相等,都用主站的tAk修正,又采取共视信息,都用同一个颗卫星钟差信息计算各个PPS触发的时标,计算出tAt与tBt的差值,达到A与B卫星导航接收机时间同步,或时间基准同源。
进一步,所述步骤2.2)中,采用伪距差分相对定位原理,针对飞行器间空间距离小于10km,以其中任意一个载体为基准站,载体上的卫星导航接收机高速运动,基准站基准位置由该飞行器上卫星导航接收机自身单点定为获得,后处理时将计算出的伪距误差存盘,用作其它从站的修正量,即是基站和从站基于同一个误差基准即公有误差,得到从站差分相对位置。
进一步,所述步骤2.2)中,导航卫星共视法包括以下步骤:2.2.1)假设A、B两个卫星导航接收机在两地,同时测量卫星i的信号,那么两卫星导航接收机的伪距观测方程式改写为:
其中,分别为A、B卫星导航接收机对卫星i的伪距测量值,两卫星导航接收机到卫星的几何距离差分处理后是已知的,而两地的接收机距离较近时,两地的电离层和对流层基本相互抵消掉;δtA、δtB表示A、B卫星导航接收机的钟差引起的伪距误差量;δti表示卫星i的钟差引起的伪距误差量;和分别表示A与B卫星导航接收机测量噪声;
2.2.2)将上述步骤中的两式相减可表示为:
其中,(δtA-δtB)为A与B两卫星接收机的钟差之差,表示A与B两卫星接收机对卫星i的观测量噪声误差;通过对多个卫星形成观测量方程,计算出A与B卫星导航接收机时钟之差,然后对A与B卫星导航接收机的时钟进行修正,得到同步修正时间。
进一步,所述步骤2)中,共视卫星的寻找方法包括以下步骤:(1)共视基准站发起共视命令报文;(2)共视卫星导航接收机接到共视报文后,发送共视查询报文1类或者2类;(3)共视基准站接收到共视查询报文后,快速匹配共视卫星号;(4)计算当前共视卫星号时差模型,并按照预先设定的频度发送共视数据报文,循环执行;(5)共视卫星导航接收机接收到共视数据报文,对当前共视卫星号PPS进行补偿输出,循环执行;(6)如果共视基准站的当前共视卫星失效,回到步骤(1);(7)如果共视接收机的当前共视卫星失效,回到步骤(2)。
本发明由于采取以上技术方案,其具有以下优点:1、本发明基于卫星共视基本原理,飞行器不需要接收对方信号,只需要将自身信息发送到地面数据处理中心,由数据处理中心计算飞行器间时间同步误差,使相互间距离20千米以内的飞行器,快速获取5ns高精度时间同步精度。2、本发明不依赖站间数据链进行双向比对时间同步,可简化系统设计难度,降低技术复杂程度。
附图说明
图1站间数据链双向比对时间同步示意图
图2是本发明的系统结构示意图
图3是本发明的技术原理分析示意图
图4是本发明的差分后处理的预处理流程图
图5是本发明的差分后处理流程图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1、图2所示,本发明提供一种多站差分后处理高精度时间同步系统,其包括安装在每个飞行器上的卫星导航接收机、卫星导航天线、数据传输系统(可以只具有单向发射功能)和数据传输系统天线,以及设置在地面接收站的地面接收机系统和数据处理系统。其中,飞行器上的数据传输系统可以采用只具有单向发射功能的数据传输系统。卫星导航接收机通过卫星导航天线接收卫星导航信号,产生PPS秒脉冲信号,产生时刻为TA0、TB0。将卫星导航接收机基带原始观测量信息以及TA0、TB0都传输至数据传输系统,由数据传输系统将接收到的基带原始观测量信息经数据传输系统天线发送到地面接收机系统,飞行器之间不需要进行通信。地面接收机系统同时接收多个飞行器上的卫星导航接收机的原始观测量信息,并传输至数据处理系统,数据处理系统同时解算多个卫星导航接收机的原始观测量信息,得到共视卫星(即多个卫星导航接收机同时使用共同的卫星信号),再经过相对差分计算,修正每个接收机触发PPS的时标为TA1,TB1,达到多站时间高度同步目的。
基于上述系统,本发明提供一种多站差分后处理高精度时间同步方法,其包括以下步骤:
1)卫星导航接收机接收导航卫星信号,产生PPS秒脉冲信号,PPS秒脉冲信号触发时刻与原始观测量采集时刻一致,卫星导航接收机实时计算出PPS的时标分别为TA0和TB0,并输出基带原始观测量信息,通过数据传输系统发送到地面数据接收系统;飞行器上的数据传输系统负责向地面接收系统发送卫星观测量信息,飞行器之间不需要通信。由于A与B是独立工作,TA0和TB0计算时没有消除各种公有误差和自身误差。因此TA0和TB0同步精度很差。
其中,PPS的时标TA0和TB0计算公式如下:
TA0=tAt+tAk+tAR+ΔtA,TB0=tBt+tBk+tBR+ΔtB (1)
式中,tAt,tBt分别代表A与B卫星导航接收机自身时钟引入的误差;tAK,tBK分别代表A与B卫星导航接收机计算时标时由多颗卫星空间公有误差(含卫星钟差)引入的误差量;tAR,tBR代表A与B卫星导航接收机触发PPS时刻的理想时标;ΔtA,ΔtB代表系统噪声。
2)数据接收系统接收多个飞行器上的原始观测量信息,并由地面接收站的数据处理系统同时解算多个卫星接收机原始观测量信息,找出共视卫星,再经过相对差分计算,修正每个卫星导航接收机触发PPS的时标;
其中,如图4所示,具体包括以下步骤:
2.1)以其中任意一个卫星导航接收机为主站,其它卫星导航接收机为从站。然后对主站进行定位解算,获得主站当前时刻位置、速度,以及PPS触发时刻的时标,并得到该时刻每个卫星差分修正量信息。
2.2)如图5所示,采用主站生成的差分和共视数据参数,处理其它从站卫星导航接收机数据,根据相对差分原理以及卫星共视授时基本原理,重新解算出各个从站差分定位结果,并重新计算出各从站PPS秒脉冲触发时刻的时标TA1与TB1。即:
TA1=tAt+tAk+tAR+Δt0,TB1=tBt+tAk+tBR+Δt1 (2)
式中,Δt0和Δt1表示差分和共视处理后卫星导航接收机A与B的系统噪声。卫星接收机两个卫星导航接收机在短基线情况下空间公有误差接近相等,都用主站的tAk修正,又采取共视信息,都用同一个颗卫星钟差信息计算各个PPS触发的时标,计算出tAt与tBt的差值,达到A与B卫星导航接收机时间同步,或时间基准同源。
由于A与B卫星导航接收机采用各自独立时钟工作,触发PPS的时刻并不在同一个时刻,有前后之分,即:经过共视和后差分处理后TA1值不是直接等于TB1值,而是经共视和差分处理,A与B卫星导航接收机使用共同的卫星信息,确保所用卫星钟差一致,相当于A与B的时间基准一致,再借助短基线下空间误差几乎相等,确保了A与B虽然不在同一时刻采样,但是几乎确保A与B使用的是同一个时钟源,确保A与B时间同步。
可以通过(TA0-TBO)-(TA1-TB1)计算验证,该值小于5ns,验证飞行器间相互时钟得到了高度同步。
上述步骤1)中,卫星导航接收机观测量误差包含接收机公有误差和自有误差两部分:
(1)卫星导航接收机公有误差:星历误差、卫星钟误差、电离层误差和对流层误差等;
(2)卫星导航接收机自身误差:内部噪声、通道延迟、多径效应;
卫星导航接收机公有误差可以通过差分方法完全予以消除,自身误差无法消除。
卫星导航接收机公有误差的消除方法如下:
公有误差可以通过伪距差分或载波相位差分消除。由于本发明应用于飞行器场景,飞行器高速飞行,姿态不稳定,导航卫星载波相位难以保证持续稳定锁定,实现载波相位差分难度大,但是完善卫星导航接收机环路跟踪可以稳定锁定扩频码相位,得到稳定伪距,可以通过伪距差分消除公有误差。
通常伪距差分实现过程是在基准站上测量采集导航卫星信息,根据基准站m(也是一个台卫星导航接收机)已知坐标(x0,y0,z0)和卫星i的坐标(xi,yi,zi),求出各卫星每一时刻与基准站的真实距离ri,如:基准站m到卫星i,则表示为
基准站伪距公式:
式中,为卫星导航接收机到第i颗卫星的伪距;为卫星导航接收机到第i颗卫星的真实距;Δtm为卫星导航接收机m的钟差;Δti为第i颗卫星钟差;为卫星i的星历误差;为卫星i传输到m的电离层误差;是卫星i传输到m的大气层误差;c表示光速。
根据卫星星历和基准站m坐标已知可求出卫星i到基准站真实距离则伪距改正数可以表示为:
卫星导航接收机n接收到卫星i的伪距为:
卫星导航接收机距离基站在200km以内情况下,利用基站伪距改正数修正卫星导航接收机n的伪距测量值:可以消除星历误差、电离层误差和对流层误差。即:
即:
式中,ΔD表示差分后的系统误差。(xn,yn,zn)表示卫星导航接收机n的坐标。
卫星坐标通过星历解出,基站坐标是已知的,差分基站可以给出每颗卫星改正数,卫星导航接收机n的观测量误差经过上述修正,第一类误差基本被消除。
上述步骤2.2)中,采用伪距差分相对定位原理,针对飞行器间空间距离小于10km,以其中任意一个载体为基准站,载体上的卫星导航接收机高速运动,基准站基准位置由该飞行器上卫星导航接收机自身单点定为获得,后处理时将计算出的伪距误差存盘,用作其它从站的修正量,即是基站和从站基于同一个误差基准(第一类误差),得到从站差分相对位置。
该思路与常规的差分定位原理一样,只不过此处的基站是运动的,基站位置时刻变化。短距离范围内,空间误差一致,以其中一个为基站,修改其它运动的卫星导航接收机,修正后的相对伪距误差减小到1m以内,考虑dop,最终相对定位精度可达1~2m。
上述步骤2.2)中,除采用相对差分定位原理获得高精度相对差分位置外,还采用导航卫星共视法提高时间同步精度。
导航卫星共视法是远距离实现数纳秒级时间同步的一种技术手段。基本原理是任意两地(或多地)的卫星导航接收机对同一颗卫星的时间信号同时测量,从而比较两地的时钟或者振荡器频率。主要步骤如下:
2.2.1)假设A、B两个卫星导航接收机在两地,同时测量卫星i的信号,那么两卫星导航接收机的伪距观测方程式改写如下:
其中,分别为A、B卫星导航接收机对卫星i的伪距测量值,两卫星导航接收机到卫星的几何距离差分处理后是已知的,而两地的接收机距离较近时,两地的电离层和对流层基本相互抵消掉;δtA、δtB表示A、B卫星导航接收机的钟差引起的伪距误差量;δti表示卫星i的钟差引起的伪距误差量;和分别表示A与B卫星导航接收机测量噪声。
2.2.2)将上述步骤中的两式相减可表示为:
其中,(δtA-δtB)为A与B两卫星接收机的钟差之差,表示A与B两卫星接收机对卫星i的观测量噪声误差,目前卫星导航接收机噪声都能控制在1m以内。通过对多个卫星形成观测量方程,可计算出A与B卫星导航接收机时钟之差,然后即可对A与B卫星导航接收机的时钟按式(2)进行修正,得到同步修正时间。
综上所述,本发明中数据传输系统单向向地面接收机系统传输数据,而不必进行飞行器间的双向通信,可以避免恶意干扰情况下,飞行器之间无法正常通信,无法高精度时间同步的可能。一般情况下飞行器远离地面接收站,容易被恶意干扰,而地面数据接收系统在地面接收站,不容易被干扰,所以即便飞行器所在位置被恶意干扰,也不影响数据发送,地面收据接收系统任然可以正常接收数据,任然能够通过地面后处理软件计算出每个时刻的高精度相对同步时间。
上述步骤2)中,共视卫星的寻找方法包括以下步骤:
(1)共视基准站发起共视命令报文。
(2)共视卫星导航接收机接到共视报文后,发送共视查询报文1类或者2类。
(3)共视基准站接收到共视查询报文后,快速匹配共视卫星号。
(4)计算当前共视卫星号时差模型,并按照预先设定的频度发送共视数据报文,循环执行。
(5)共视卫星导航接收机接收到共视数据报文,对当前共视卫星号PPS进行补偿输出,循环执行。
(6)如果共视基准站的当前共视卫星失效,回到步骤(1)。
(7)如果共视接收机的当前共视卫星失效,回到步骤(2)。
上述各实施例仅用于说明本发明,各个步骤都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别步骤进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (8)
1.一种多站差分后处理高精度时间同步系统,其特征在于:包括安装在每个飞行器上的卫星导航接收机、卫星导航天线、数据传输系统和数据传输系统天线,以及设置在地面接收站的地面接收机系统和数据处理系统;所述卫星导航接收机通过所述卫星导航天线接收卫星导航信号,产生PPS秒脉冲信号,产生时刻为TA0、TB0,将所述卫星导航接收机基带原始观测量信息以及TA0、TB0都传输至所述数据传输系统,由所述数据传输系统将接收到的基带原始观测量信息经所述数据传输系统天线发送所述到地面接收机系统,飞行器之间不需要进行通信;所述地面接收机系统同时接收多个飞行器上的所述卫星导航接收机的原始观测量信息,并传输至所述数据处理系统,所述数据处理系统同时解算多个所述卫星导航接收机的原始观测量信息,得到共视卫星,再经过相对差分计算,修正每个接收机触发PPS的时标为TA1,TB1,实现多站时间高度同步。
2.一种基于如权利要求1所述系统的多站差分后处理高精度时间同步方法,其特征在于包括以下步骤:
1)卫星导航接收机接收导航卫星信号,产生PPS秒脉冲信号,卫星导航接收机实时计算出PPS的时标分别为TA0和TB0,并输出基带原始观测量信息,通过数据传输系统发送到地面数据接收系统;飞行器上的数据传输系统负责向地面接收系统发送卫星观测量信息;
2)数据接收系统接收多个飞行器上的原始观测量信息,并由地面接收站的数据处理系统同时解算多个卫星接收机原始观测量信息,找出共视卫星,再经过相对差分计算,修正每个卫星导航接收机触发PPS的时标。
3.如权利要求2所述同步方法,其特征在于:所述步骤1)中,PPS的时标TA0和TB0计算公式如下:
TA0=tat+tAk+tAR+ΔtA,TB0=tBt+tBk+tBR+ΔtB
式中,tat,tBt分别代表A与B卫星导航接收机自身时钟引入的误差;tAK,tBK分别代表A与B卫星导航接收机计算时标时由多颗卫星空间公有误差引入的误差量;tAR,tBR代表A与B卫星导航接收机触发PPS时刻的理想时标;ΔtA,ΔtB代表系统噪声。
4.如权利要求2或3所述同步方法,其特征在于:所述步骤1)中,卫星导航接收机观测量误差包含接收机公有误差和自有误差,卫星导航接收机公有误差通过伪距差分消除,其方法如下:
卫星导航接收机距离基站在200km以内情况下,利用基站伪距改正数修正卫星导航接收机n的伪距测量值:消除公有误差,即:
式中,ΔD表示差分后的系统误差;(xn,yn,zn)表示卫星导航接收机n的坐标;表示伪距改正数;表示卫星导航接收机n接收到卫星i的伪距;卫星i的坐标为(xi,yi,zi);卫星导航接收机n的观测量误差经过上述修正,公有误差基本被消除。
5.如权利要求2所述同步方法,其特征在于:所述步骤2)中,每个卫星导航接收机触发PPS的时标修正包括以下步骤:
2.1)以其中任意一个卫星导航接收机为主站,其它卫星导航接收机为从站,然后对主站进行定位解算,获得主站当前时刻位置、速度,以及PPS触发时刻的时标,并得到该时刻每个卫星差分修正量信息;
2.2)采用主站生成的差分和共视数据参数,处理其它从站卫星导航接收机数据,根据相对差分原理以及卫星共视授时基本原理,重新解算出各个从站差分定位结果,并重新计算出各从站PPS秒脉冲触发时刻的时标TA1与TB1,即:
TA1=tAt+tAk+tAR+Δt0,TB1=tBt+tAk+tBR+Δt1
式中,tat,tBt分别代表A与B卫星导航接收机自身时钟引入的误差;tAK,tBK分别代表A与B卫星导航接收机计算时标时由多颗卫星空间公有误差引入的误差量;tAR,tBR代表A与B卫星导航接收机触发PPS时刻的理想时标;ΔtA,ΔtB代表系统噪声;Δt0和Δt1表示差分和共视处理后卫星导航接收机A与B的系统噪声;卫星接收机两个卫星导航接收机在短基线情况下空间公有误差接近相等,都用主站的tAk修正,又采取共视信息,都用同一个颗卫星钟差信息计算各个PPS触发的时标,计算出tAt与tBt的差值,达到A与B卫星导航接收机时间同步,或时间基准同源。
6.如权利要求5所述同步方法,其特征在于:所述步骤2.2)中,采用伪距差分相对定位原理,针对飞行器间空间距离小于10km,以其中任意一个载体为基准站,载体上的卫星导航接收机高速运动,基准站基准位置由该飞行器上卫星导航接收机自身单点定为获得,后处理时将计算出的伪距误差存盘,用作其它从站的修正量,即是基站和从站基于同一个误差基准即公有误差,得到从站差分相对位置。
7.如权利要求5所述同步方法,其特征在于:所述步骤2.2)中,导航卫星共视法包括以下步骤:
2.2.1)假设A、B两个卫星导航接收机在两地,同时测量卫星i的信号,那么两卫星导航接收机的伪距观测方程式改写为:
其中,分别为A、B卫星导航接收机对卫星i的伪距测量值,两卫星导航接收机到卫星的几何距离差分处理后是已知的,而两地的接收机距离较近时,两地的电离层和对流层基本相互抵消掉;δtA、δtB表示A、B卫星导航接收机的钟差引起的伪距误差量;δti表示卫星i的钟差引起的伪距误差量;和分别表示A与B卫星导航接收机测量噪声;
2.2.2)将上述步骤中的两式相减可表示为:
其中,(δtA-δtB)为A与B两卫星接收机的钟差之差,表示A与B两卫星接收机对卫星i的观测量噪声误差;通过对多个卫星形成观测量方程,计算出A与B卫星导航接收机时钟之差,然后对A与B卫星导航接收机的时钟进行修正,得到同步修正时间。
8.如权利要求1所述同步方法,其特征在于:所述步骤2)中,共视卫星的寻找方法包括以下步骤:
(1)共视基准站发起共视命令报文;
(2)共视卫星导航接收机接到共视报文后,发送共视查询报文1类或者2类;
(3)共视基准站接收到共视查询报文后,快速匹配共视卫星号;
(4)计算当前共视卫星号时差模型,并按照预先设定的频度发送共视数据报文,循环执行;
(5)共视卫星导航接收机接收到共视数据报文,对当前共视卫星号PPS进行补偿输出,循环执行;
(6)如果共视基准站的当前共视卫星失效,回到步骤(1);
(7)如果共视接收机的当前共视卫星失效,回到步骤(2)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910091296.5A CN109655846B (zh) | 2019-01-30 | 2019-01-30 | 一种多站差分后处理高精度时间同步方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910091296.5A CN109655846B (zh) | 2019-01-30 | 2019-01-30 | 一种多站差分后处理高精度时间同步方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109655846A true CN109655846A (zh) | 2019-04-19 |
CN109655846B CN109655846B (zh) | 2021-03-26 |
Family
ID=66121800
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910091296.5A Expired - Fee Related CN109655846B (zh) | 2019-01-30 | 2019-01-30 | 一种多站差分后处理高精度时间同步方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109655846B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045395A (zh) * | 2019-05-17 | 2019-07-23 | 北京和协航电信息科技有限公司 | 一种基于导航卫星信号的检靶方法、系统和设备 |
CN110275134A (zh) * | 2019-06-27 | 2019-09-24 | 清华大学 | 一种基于虚拟到达频差的非共视连续信号无源定位方法 |
CN110289903A (zh) * | 2019-06-26 | 2019-09-27 | 清华大学 | 一种基于链路测量的卫星网络自主时间同步方法和系统 |
CN110456397A (zh) * | 2019-07-29 | 2019-11-15 | 桂林电子科技大学 | 一种多天线超短基线定位监测方法、装置及存储介质 |
CN110764401A (zh) * | 2019-10-29 | 2020-02-07 | 北京无线电计量测试研究所 | 一种船载时间同步校准设备 |
CN110907971A (zh) * | 2019-12-03 | 2020-03-24 | 深圳华大北斗科技有限公司 | 高空设备的卫星定位方法、装置、计算机设备和存储介质 |
CN111044972A (zh) * | 2019-12-12 | 2020-04-21 | 北京航空航天大学 | 基于gnss精密时间同步的飞行器时差定位方法及系统 |
CN111142122A (zh) * | 2020-01-08 | 2020-05-12 | 中国电子科技集团公司第五十四研究所 | 一种基于bds相对定位的高精度授时方法 |
CN111175797A (zh) * | 2020-02-02 | 2020-05-19 | 上海交通大学 | 基于虚拟质心的多gnss接收机协同导航方法 |
CN111337961A (zh) * | 2020-05-21 | 2020-06-26 | 深圳市西博泰科电子有限公司 | 提高车辆在城市中定位精度的方法、装置、设备及介质 |
CN111464397A (zh) * | 2020-04-08 | 2020-07-28 | 清华大学 | 一种双向距离和钟差测量的方法和系统 |
CN112422222A (zh) * | 2020-11-04 | 2021-02-26 | 深圳市三七智联科技有限公司 | 授时模组与基站的接口方法、卫星导航系统及存储介质 |
CN112884342A (zh) * | 2021-03-10 | 2021-06-01 | 陕西九州遥感信息技术有限公司 | 一种水色卫星大气层顶辐射产品质量评价与交叉定标方法 |
CN113489509A (zh) * | 2021-07-06 | 2021-10-08 | 阎镜予 | 大规模gnss伪卫星间时频同步方法、装置 |
CN113721277A (zh) * | 2021-09-24 | 2021-11-30 | 苏州常棣华电子科技有限公司 | 一种高精度卫星定位方法(grp)和装置 |
CN114563804A (zh) * | 2021-12-27 | 2022-05-31 | 中国人民解放军空军工程大学 | 一种gnss/ins紧组合导航系统的自适应容错方法 |
CN117148405A (zh) * | 2023-09-03 | 2023-12-01 | 南京林业大学 | 基于时间同步算法的环岛下自动驾驶车辆高精度定位算法 |
CN118226467A (zh) * | 2024-05-27 | 2024-06-21 | 江苏深蓝航天有限公司 | 一种时间同步方法、装置、电子设备、介质以及产品 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6983160B2 (en) * | 2001-10-25 | 2006-01-03 | Motorola, Inc. | Base site and method for GPS-based regional time synchronization |
CN102890279A (zh) * | 2011-07-20 | 2013-01-23 | 中国科学院电子学研究所 | 一种利用全球定位系统进行纳秒级精确时间同步的方法 |
CN104125029A (zh) * | 2013-04-23 | 2014-10-29 | 广州南方卫星导航仪器有限公司 | 北斗精密授时技术 |
CN104199280A (zh) * | 2014-09-23 | 2014-12-10 | 中国电子科技集团公司第二十九研究所 | 一种基于差分gps的时间同步误差测量方法 |
CN105607089A (zh) * | 2016-01-18 | 2016-05-25 | 西安测绘研究所 | 一种基于伪距辅助载波相位的北斗共视时间传递的方法 |
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN108333604A (zh) * | 2017-12-27 | 2018-07-27 | 和芯星通科技(北京)有限公司 | 一种利用卫星定位的方法和装置、卫星授时方法和装置 |
CN108931915A (zh) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | 利用导航卫星的授时方法和装置、计算机可读存储介质 |
-
2019
- 2019-01-30 CN CN201910091296.5A patent/CN109655846B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6983160B2 (en) * | 2001-10-25 | 2006-01-03 | Motorola, Inc. | Base site and method for GPS-based regional time synchronization |
CN102890279A (zh) * | 2011-07-20 | 2013-01-23 | 中国科学院电子学研究所 | 一种利用全球定位系统进行纳秒级精确时间同步的方法 |
CN104125029A (zh) * | 2013-04-23 | 2014-10-29 | 广州南方卫星导航仪器有限公司 | 北斗精密授时技术 |
CN104199280A (zh) * | 2014-09-23 | 2014-12-10 | 中国电子科技集团公司第二十九研究所 | 一种基于差分gps的时间同步误差测量方法 |
CN105607089A (zh) * | 2016-01-18 | 2016-05-25 | 西安测绘研究所 | 一种基于伪距辅助载波相位的北斗共视时间传递的方法 |
CN107579794A (zh) * | 2017-08-21 | 2018-01-12 | 中国科学院国家授时中心 | 一种基于北斗geo导航卫星的精密共视时间频率传递方法 |
CN108333604A (zh) * | 2017-12-27 | 2018-07-27 | 和芯星通科技(北京)有限公司 | 一种利用卫星定位的方法和装置、卫星授时方法和装置 |
CN108931915A (zh) * | 2018-05-08 | 2018-12-04 | 和芯星通科技(北京)有限公司 | 利用导航卫星的授时方法和装置、计算机可读存储介质 |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110045395A (zh) * | 2019-05-17 | 2019-07-23 | 北京和协航电信息科技有限公司 | 一种基于导航卫星信号的检靶方法、系统和设备 |
CN110289903B (zh) * | 2019-06-26 | 2020-04-14 | 清华大学 | 一种基于链路测量的卫星网络自主时间同步方法和系统 |
CN110289903A (zh) * | 2019-06-26 | 2019-09-27 | 清华大学 | 一种基于链路测量的卫星网络自主时间同步方法和系统 |
CN110275134A (zh) * | 2019-06-27 | 2019-09-24 | 清华大学 | 一种基于虚拟到达频差的非共视连续信号无源定位方法 |
CN110275134B (zh) * | 2019-06-27 | 2021-03-05 | 清华大学 | 一种基于虚拟到达频差的非共视连续信号无源定位方法 |
CN110456397B (zh) * | 2019-07-29 | 2021-02-09 | 桂林电子科技大学 | 一种多天线超短基线定位监测方法、装置及存储介质 |
CN110456397A (zh) * | 2019-07-29 | 2019-11-15 | 桂林电子科技大学 | 一种多天线超短基线定位监测方法、装置及存储介质 |
CN110764401B (zh) * | 2019-10-29 | 2021-11-16 | 北京无线电计量测试研究所 | 一种船载时间同步校准设备 |
CN110764401A (zh) * | 2019-10-29 | 2020-02-07 | 北京无线电计量测试研究所 | 一种船载时间同步校准设备 |
CN110907971A (zh) * | 2019-12-03 | 2020-03-24 | 深圳华大北斗科技有限公司 | 高空设备的卫星定位方法、装置、计算机设备和存储介质 |
CN110907971B (zh) * | 2019-12-03 | 2022-07-12 | 深圳华大北斗科技股份有限公司 | 高空设备的卫星定位方法、装置、计算机设备和存储介质 |
CN111044972A (zh) * | 2019-12-12 | 2020-04-21 | 北京航空航天大学 | 基于gnss精密时间同步的飞行器时差定位方法及系统 |
CN111142122A (zh) * | 2020-01-08 | 2020-05-12 | 中国电子科技集团公司第五十四研究所 | 一种基于bds相对定位的高精度授时方法 |
CN111175797A (zh) * | 2020-02-02 | 2020-05-19 | 上海交通大学 | 基于虚拟质心的多gnss接收机协同导航方法 |
CN111175797B (zh) * | 2020-02-02 | 2023-11-03 | 上海交通大学 | 基于虚拟质心的多gnss接收机协同导航方法 |
CN111464397A (zh) * | 2020-04-08 | 2020-07-28 | 清华大学 | 一种双向距离和钟差测量的方法和系统 |
CN111337961A (zh) * | 2020-05-21 | 2020-06-26 | 深圳市西博泰科电子有限公司 | 提高车辆在城市中定位精度的方法、装置、设备及介质 |
CN112422222B (zh) * | 2020-11-04 | 2023-01-10 | 深圳市三七智联科技有限公司 | 授时模组与基站的接口方法、卫星导航系统及存储介质 |
CN112422222A (zh) * | 2020-11-04 | 2021-02-26 | 深圳市三七智联科技有限公司 | 授时模组与基站的接口方法、卫星导航系统及存储介质 |
CN112884342A (zh) * | 2021-03-10 | 2021-06-01 | 陕西九州遥感信息技术有限公司 | 一种水色卫星大气层顶辐射产品质量评价与交叉定标方法 |
CN112884342B (zh) * | 2021-03-10 | 2024-03-12 | 陕西九州遥感信息技术有限公司 | 一种水色卫星大气层顶辐射产品质量评价与交叉定标方法 |
CN113489509A (zh) * | 2021-07-06 | 2021-10-08 | 阎镜予 | 大规模gnss伪卫星间时频同步方法、装置 |
CN113721277A (zh) * | 2021-09-24 | 2021-11-30 | 苏州常棣华电子科技有限公司 | 一种高精度卫星定位方法(grp)和装置 |
CN114563804A (zh) * | 2021-12-27 | 2022-05-31 | 中国人民解放军空军工程大学 | 一种gnss/ins紧组合导航系统的自适应容错方法 |
CN114563804B (zh) * | 2021-12-27 | 2024-05-24 | 中国人民解放军空军工程大学 | 一种gnss/ins紧组合导航系统的自适应容错方法 |
CN117148405A (zh) * | 2023-09-03 | 2023-12-01 | 南京林业大学 | 基于时间同步算法的环岛下自动驾驶车辆高精度定位算法 |
CN118226467A (zh) * | 2024-05-27 | 2024-06-21 | 江苏深蓝航天有限公司 | 一种时间同步方法、装置、电子设备、介质以及产品 |
Also Published As
Publication number | Publication date |
---|---|
CN109655846B (zh) | 2021-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109655846B (zh) | 一种多站差分后处理高精度时间同步方法及系统 | |
CN110031881B (zh) | 高精度星间激光测距辅助精密单点定位的方法 | |
CN108931915B (zh) | 利用导航卫星的授时方法和装置、计算机可读存储介质 | |
EP2673658B1 (en) | Method and system for determining clock corrections | |
EP2661933B1 (en) | System and method for time synchronizing wireless network access points | |
US7616682B2 (en) | Method and device for chronologically synchronizing a location network | |
AU2011325913B2 (en) | Wide area positioning system | |
CN108196272A (zh) | 一种基于实时精密单点定位的卫星导航定位装置及方法 | |
WO2014171999A2 (en) | Method and system for high-accuracy differential tracking of global positioning system (gps) receivers | |
EP1709460A1 (en) | Transfer of calibrated time information in a mobile terminal | |
US7679554B1 (en) | Communications link time transfer to improve navigation system accuracy | |
US7605749B2 (en) | Multifrequency line biases for multifrequency GNSS receivers | |
CN111044972B (zh) | 基于gnss精密时间同步的飞行器时差定位方法及系统 | |
CN110749907A (zh) | 一种基于北斗动定位中接收机的钟差补偿方法及其系统 | |
CN111998849A (zh) | 一种基于惯性导航系统的差分动动定位的方法 | |
CN112630811B (zh) | 一种实时ppp-rtk组合定位方法 | |
WO2009130305A1 (en) | Method of positioning using satellites | |
CN114895327A (zh) | 一种多系统upd产品估计与检核方法 | |
JP2010060421A (ja) | 移動体用測位システム及びgnss受信装置 | |
CN108549094A (zh) | 基于中继卫星辅助的地基导航系统 | |
JP3524018B2 (ja) | 測位システム及び測位方法 | |
RU2825248C1 (ru) | Способ определения местоположения абонентских терминалов, перемещающихся в зоне покрытия локальной системы навигации | |
KR20130111603A (ko) | 무선 네트워크 액세스 포인트들을 시간 동기화하는 시스템 및 방법 | |
Aloi et al. | A relative technique for characterization of PCV error of large aperture antennas using GPS data | |
CN118393531A (zh) | 集成gnss/uwb定位系统的校准方法、校准装置 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210326 Termination date: 20220130 |
|
CF01 | Termination of patent right due to non-payment of annual fee |