CN108181633A - 一种gnss时间频率传递接收机及接收方法 - Google Patents

一种gnss时间频率传递接收机及接收方法 Download PDF

Info

Publication number
CN108181633A
CN108181633A CN201711392705.2A CN201711392705A CN108181633A CN 108181633 A CN108181633 A CN 108181633A CN 201711392705 A CN201711392705 A CN 201711392705A CN 108181633 A CN108181633 A CN 108181633A
Authority
CN
China
Prior art keywords
gps
time
positioning system
galileo
coordinate
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.)
Pending
Application number
CN201711392705.2A
Other languages
English (en)
Inventor
程禹
吴限德
白文彬
韩华
田光远
谢亚恩
孙楚琦
路胜卓
陈卫东
倪蜂棋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201711392705.2A priority Critical patent/CN108181633A/zh
Publication of CN108181633A publication Critical patent/CN108181633A/zh
Pending legal-status Critical Current

Links

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/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
    • 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/30—Acquisition or tracking or demodulation of signals transmitted by the system code related
    • 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时间频率传递接收机及接收方法，属于卫星导航***接收机时间比对研究领域。其特征在于：时频单元产生射频单元和信号处理单元所需要的时频信号，同时输出PPS信号；射频单元为信号处理单元提供中频信号；信号接收处理单元完成导航信号的采样、捕获、跟踪、导航电文解调、伪距测量并上报各种观测量；监控单元对观测数据进行传输和对各单元下发控制指令。多***融合共视时间比对方法统一了现有的几种定位***的时间和空间坐标；计算GNSS卫星的位置参数；进而计算出时间差值，最后对时间比对结果进行卡尔曼滤波和RTS事后处理。

Description

一种GNSS时间频率传递接收机及接收方法
技术领域
本发明涉及一种GNSS接收机及多***融合共视时间比对方法，属于卫星导航***接收机时间比对研究领域。
背景技术
在定位导航精度不断提高的背景下，高精度的时间同步已经成为各个GNSS卫星定位导航***中的关键技术，是定位导航精度的重要保证。时间比对算法是实现时间同步技术的关键，是确定不同地面接收机时间差的重要方法。目前现有的时间比对算法包括共视法时间比对，全视法时间比对，载波相位法时间比对等算法。Allan和Weiss等人提出共视时间比对算法；Jung与Petit提出了全视时间比对算法。全视法时间比对需要多个接收机同时观测多颗卫星，这同样对接收机和比对算法有较高要求；载波相位法虽然精度较高，但观测器材的成本也相对较高，且算法复杂，不易实现。
发明内容
本发明的目的是这样实现的：
一种GNSS时间频率传递接收机，包括天线单元、时频单元、射频单元、信号处理单元和监控单元，其特征在于：时频单元产生射频单元和信号处理单元所需要的时频信号，同时输出PPS信号；射频单元把经过天线单元处理的射频信号进行分离，经下变频及低噪声放大、AGC控制后，为信号处理单元提供中频信号；信号接收处理单元完成导航信号的采样、捕获、跟踪、导航电文解调、伪距测量并上报各种观测量；监控单元对观测数据进行传输和对各单元下发控制指令。
信号处理单元包括信号捕获模块、通道跟踪模块和信息处理模块；信号捕获模块接收用户配置参数；跟踪通道组包含伪码发生器、载波NCO、码NCO和累加器；信息处理模块的工作过程包含以下步骤：
步骤一分别通过GPS定位***、COMPASS定位***、Galileo定位***和GLONASS定位***观测各自***内的卫星，各获得一组时间比对数据；
步骤二统一时间***，包括GPS定位***与COMPASS定位***时间转换、GPS定位***与Galileo定位***时间转换、GPS定位***与GLONASS定位***时间转换；
步骤三对步骤二得到的统一后的时间数据进行解码和预处理获得GNSS卫星轨道位置，并对获得的GNSS卫星位置坐标进行地球自转效应修正；
步骤四统一空间***；
步骤五根据统一的空间***，计算GNSS卫星与两个GNSS接收机的几何距离；计算卫星与接收机然后计算信号传播过程中的对流层时延；将这些参数代入到由伪距观测方程推导出的地面接收机时间差观测方程中，进而计算出时间差值；
步骤六对时间比对结果进行Kalman滤波和RTS事后处理。
步骤二中所述的GPS定位***与COMPASS定位***时间转换过程为：
COMPASS定位***和GPS定位***时间的换算公式为：
tGPS＝tCOM-ΔtCOM,GPS
ΔtCOM,GPS＝A0GPS+A1GPS×tCOM
式中，tCOM为指定用户计算的COMPASS时间，A0GPS是COMPASS相对于GPS定位***时间的钟差，A1GPS是COMPASS定位***相对于GPS定位***时间的钟速，取周内秒计数部分；
步骤二中所述的GPS定位***与Galileo定位***时间转换过程为：
Galileo导航电文给出Galileo***时间与GPS时间之间实现转换所需的A0G、A1G、toG、WNoG四个参数，其中，toG、WNoG指出这套参数的参考时间，A0G、A1G为关于GGTO的一线性模型中的两个系数；给定一个GPS时间tGPS及其所属的星期数WNGPS，得此时的Galileo/GPS***时间偏差ΔtGAL,GPS为：
ΔtGAL,GPS＝A0G+A1G(tGPS-toG+604800|WNGPS-WN0G|64)[s]
式中，偏差ΔtGAL,GPS定义为相同时刻的Galileo***时间tGAL减去GPS时间tGPS，即
tGPS＝tGAL-ΔtGAL,GPS
|i|64代表i对64的模；GGTO的准确度2σ预计为5ns，ΔtGAL,GPS作为一个未知量在Galileo/GPS联合定位、定时计算中；
步骤二中所述的GPS定位***与GLONASS定位***时间转换过程为：
GLONASS-M型卫星在其播发的每一帧导航电文的第5串中提供GLONASS时间tGLO相对于GPS时间tGPS的秒内偏差值τGPS；GPS NAV电文提供tGPS与UTC之间由跳秒引入的整数秒差异值ΔtLS；再加上考虑tGLO相对于UTC的3小时差异，得如下tGLO和tGPS之间的转换关系式：
tGLO＝(tGPS-ΔtLS)-τGPS+3×3600[s]
GPS和GLONASS两***时间之间的秒内偏差值τGPS从GLONASS导航电文中解调获得。
步骤四中所述的统一空间***为：
GPS定位***与GLONASS定位***坐标系转换：
设空间一点P在一个直角坐标系S中的位置坐标为(x,y,z)，则该点在另一个不同却几乎平行的直角坐标系S′中的位置坐标(x′,y′,z′)为：
式中，坐标平移量(Δx,Δy,Δz)为坐标系S的原点在坐标系S′中的坐标，δε、δψ和δω为坐标系S分别绕X、Y和Z坐标轴旋转而能得到S′相应坐标轴的所需旋转角度，δs体现为坐标尺度因子；式中等号右边(x,y,z)为一个PZ-90坐标值，等号左边(x′,y′,z′)为相应的WGS-84坐标值。
步骤六包含以下步骤：
步骤一状态变量Xk与其Kalman估计之间的均方误差矩阵称为估计误差的协方差矩阵，用Ck表示；给定C0一个初值，根据下式
计算出P1；其中，Pk为状态变量Xk与其在无观测噪声与模型噪声条件下的估计之间的均方误差阵，Qk为1×1阶模型噪声叭的协方差阵；
步骤二得到P1后，根据Kalman增益矩阵Gk的表示式
求得G1，其中Rk为1×1阶观测噪声Vk的协方差阵；
步骤三根据下式
得k＝1时刻的状态变量估计值即k＝1时刻钟差Kalman估计值；
步骤四将P1带入下式
Ck＝(I-GkHk)Pk
求得k＝1时刻的1×1阶估计误差的协方差阵C1；然后，进入下次循环；
在近实时共视中，对异地钟差X0有初步的估计，用这个值对滤波器进行初始化，使滤波器的收敛速度加快；此时，估计误差的协方差矩阵初值取为：
步骤五在卡尔曼滤波的基础上，利用整个时间间隔内所有量测数据得到状态的最小方差估计，对结果进行RTS固定区间最优平滑，平滑解算需要在滤波过程中实时存储数据，所存储的数据为4个矩阵，分别为估计值一步转移阵估计的均方差阵以及一步预测均方差阵
平滑公式为：
K＝N-1，N-2，……2，1，0；
其中
式中，为RTS平滑值；Ks,k为滤波增益；Ps,k为协方差阵。
一种嵌入多***融合共视时间比对算法的GNSS时间频率传递接收方法，包含以下步骤：
步骤一分别通过GPS定位***、COMPASS定位***、Galileo定位***和GLONASS定位***观测各自***内的卫星，各获得一组时间比对数据；
步骤二统一时间***，包括GPS定位***与COMPASS定位***时间转换、GPS定位***与Galileo定位***时间转换、GPS定位***与GLONASS定位***时间转换；
步骤三对步骤二得到的统一后的时间数据进行解码和预处理获得GNSS卫星轨道位置，并对获得的GNSS卫星位置坐标进行地球自转效应修正；
步骤四统一空间***；
步骤五根据统一的空间***，计算GNSS卫星与两个GNSS接收机的几何距离；计算卫星与接收机然后计算信号传播过程中的对流层时延；将这些参数代入到由伪距观测方程推导出的地面接收机时间差观测方程中，进而计算出时间差值；
步骤六对时间比对结果进行Kalman滤波和RTS事后处理。
步骤二中所述的GPS定位***与COMPASS定位***时间转换过程为：
COMPASS定位***和GPS定位***时间的换算公式为：
tGPS＝tCOM-ΔtCOM,GPS
ΔtCOM,GPS＝A0GPS+A1GPS×tCOM
式中，tCOM为指定用户计算的COMPASS时间，A0GPS是COMPASS相对于GPS定位***时间的钟差，A1GPS是COMPASS定位***相对于GPS定位***时间的钟速，取周内秒计数部分；
步骤二中所述的GPS定位***与Galileo定位***时间转换过程为：
Galileo导航电文给出Galileo***时间与GPS时间之间实现转换所需的A0G、A1G、toG、WNoG四个参数，其中，toG、WNoG指出这套参数的参考时间，A0G、A1G为关于GGTO的一线性模型中的两个系数；给定一个GPS时间tGPS及其所属的星期数WNGPS，得此时的Galileo/GPS***时间偏差ΔtGAL,GPS为：
ΔtGAL,GPS＝A0G+A1G(tGPS-toG+604800|WNGPS-WN0G|64)[s]
式中，偏差ΔtGAL,GPS定义为相同时刻的Galileo***时间tGAL减去GPS时间tGPS，即
tGPS＝tGAL-ΔtGAL,GPS
|i|64代表i对64的模；GGTO的准确度2σ预计为5ns，ΔtGAL,GPS作为一个未知量在Galileo/GPS联合定位、定时计算中；
步骤二中所述的GPS定位***与GLONASS定位***时间转换过程为：
GLONASS-M型卫星在其播发的每一帧导航电文的第5串中提供GLONASS时间tGLO相对于GPS时间tGPS的秒内偏差值τGPS；GPS NAV电文提供tGPS与UTC之间由跳秒引入的整数秒差异值ΔtLS；再加上考虑tGLO相对于UTC的3小时差异，得如下tGLO和tGPS之间的转换关系式：
tGLO＝(tGPS-ΔtLS)-τGPS+3×3600[s]
GPS和GLONASS两***时间之间的秒内偏差值τGPS从GLONASS导航电文中解调获得。
步骤四中所述的统一空间***为：
GPS定位***与GLONASS定位***坐标系转换：
设空间一点P在一个直角坐标系S中的位置坐标为(x,y,z)，则该点在另一个不同却几乎平行的直角坐标系S′中的位置坐标(x′,y′,z′)为：
式中，坐标平移量(Δx,Δy,Δz)为坐标系S的原点在坐标系S′中的坐标，δε、δψ和δω为坐标系S分别绕X、Y和Z坐标轴旋转而能得到S′相应坐标轴的所需旋转角度，δs体现为坐标尺度因子；式中等号右边(x,y,z)为一个PZ-90坐标值，等号左边(x′,y′,z′)为相应的WGS-84坐标值。
步骤六包含以下步骤：
步骤一状态变量Xk与其Kalman估计之间的均方误差矩阵称为估计误差的协方差矩阵，用Ck表示；给定C0一个初值，根据下式
计算出P1；其中，Pk为状态变量Xk与其在无观测噪声与模型噪声条件下的估计之间的均方误差阵，Qk为1×1阶模型噪声叭的协方差阵；
步骤二得到P1后，根据Kalman增益矩阵Gk的表示式
求得G1，其中Rk为1×1阶观测噪声Vk的协方差阵；
步骤三根据下式
得k＝1时刻的状态变量估计值即k＝1时刻钟差Kalman估计值；
步骤四将P1带入下式
Ck＝(I-GkHk)Pk
求得k＝1时刻的1×1阶估计误差的协方差阵C1；然后，进入下次循环；
在近实时共视中，对异地钟差X0有初步的估计，用这个值对滤波器进行初始化，使滤波器的收敛速度加快；此时，估计误差的协方差矩阵初值取为：
步骤五在卡尔曼滤波的基础上，利用整个时间间隔内所有量测数据得到状态的最小方差估计，对结果进行RTS固定区间最优平滑，平滑解算需要在滤波过程中实时存储数据，所存储的数据为4个矩阵，分别为估计值一步转移阵估计的均方差阵以及一步预测均方差阵
平滑公式为：
K＝N-1，N-2，……2，1，0；
其中
式中，为RTS平滑值；Ks,k为滤波增益；Ps,k为协方差阵。
本发明的有益效果：第一，接收机无需连接计算机便可以独立进行时间比对工作。第二，采用多***融合共视时间比对算法，在易于实现的前提下，保证了单***观测卫星数目较少时的时间比对精度。第三，加入Kalman滤波和RTS事后处理的时间比对相比于单纯的共视时间比对算法精度更高。
附图说明
附图1GNSS接收机原理图。
附图2是多***融合时间比对方法流程图。
具体实施方式
下面结合附图对本发明做详细地描述：
本发明属于卫星导航***接收机时间比对研究领域，具体涉及到接收机硬件设计：时频单元，射频单元，信号处理单元，监控单元；时间空间***转换算法和数据的后期处理。多***融合共视时间比对算法的关键是时间和空间***转换算法，获得实时时间差后对结果进行滤波和事后处理。
本发明包括GNSS接收机硬件设计，图1为共视接收机原理图。
在GNSS共视接收机中嵌入多***融合共视时间比对算法模块。融合了GPS，GLONASS，Galileo和COMPASS四个卫星导航***的相关接收数据，融合解算出接收机间的时间差。
GNSS共视接收机主要包括时频单元，射频单元，信号处理单元和监控单元等部分。
时频单元
产生射频单元和基带信号处理单元所需要的时频信号，同时输出PPS信号。
射频单元
射频单元把经过天线模块处理的射频信号进行分离，经下变频及低噪声放大、AGC控制后，为信号处理模块提供满足一定信噪比(S/N)和幅度要求的中频信号。
由于入口信号电平低(比白噪声电平低很多)，为了给后端处理部分提供足够幅度的中频信号，变频器单元与天线单元通道净增益应大于110dB，其中天线部分增益大约40dB，如果不考虑传输电缆的影响，两路信号的增益都应该大于70dB。故要求射频芯片配置成低中频及零中频输出接口方式。
信号处理单元
信号接收处理单元主要是完成导航信号的采样、捕获、跟踪、导航电文解调、伪距测量并上报各种观测量等任务。
接收机信号处理单元包括信号捕获模块、通道跟踪模块(含跟踪通道组和流程控制单元)和信息处理模块。
捕获模块接收用户配置参数，根据不同的配置组合，可以完成所有频点的伪码捕获工作，为了扩展方便，该模块把所有频点的码速率和载波频率都设置成可配置的，这样该模块更容易实现各频点互操作功能，而且也容易适应不同的射频信道。
跟踪通道组包含伪码发生器、载波NCO、码NCO和累加器功能。
信息处理模块的主要工作是接收信号处理单元提供的卫星观测量和导航电文，进行自主完好性检测，并对多***不同的电文参数、时间***、坐标***进行统一，然后完成定位解算并输出结果等。
监控单元
监控单元主要对观测数据进行传输和对各单元下发控制指令。
本发明包括伪距方程的建立，延迟修正值的算法建模以及滤波和事后处理。图2为时间比对算法流程图。
多***融合共视时间比对算法，包括时间和空间***统一算法，并对实时时间比对就结果进行滤波和事后处理。多***共视时间比对算法是基于单***双频接收机时间比对算法设计的，首先，各***先通过观测***内卫星，各获得一组时间比对数据，并根据数据精度分配权重，加权得出时间比对的最终结果。各***的时间比对算法依然是以导航电文、广播星历和接收机位置为输入，通过对星历电文的解码和预处理获得卫星轨道位置。对获得的卫星位置坐标进行地球自转效应修正并计算GNSS卫星与两个接收机的几何距离。计算卫星与接收机然后计算信号传播过程中的对流层时延，最终将这些参数代入到由伪距观测方程推导出的地面接收机时间差观测方程中，进而计算出时间差值最后，对时间比对结果进行卡尔曼滤波和RTS事后处理。
首先，统一时间和空间***。
GPS***与COMPASS***时间转换算法
COMPASS***采用的时间***，与GPS时间类似，COMPASS电文提供了其***时间同GPS***时间的同步参数，如表1所示。
表1 COMPASS与GPS的时间同步参数
表1中，A0GPS是COMPASS相对于GPS***时间的钟差，A1GPS是COMPASS相对于GPS***时间的钟速，COMPASS和GPS***时间的换算公式为：
(1)tGPS＝tCOM-ΔtCOM,GPS (1-1)
(2)ΔtCOM,GPS＝A0GPS+A1GPS×tCOM (1-2)
式中，tCOM为指定用户计算的COMPASS时间，取周内秒计数部分。
GPS***与Galileo***时间转换
Galileo***时间(GST)基于一个连续运行的原子时，它通过对一系列原子频率标准的整合来维持，其中氢原子钟被作为主钟。Galileo***时间与国际原子时(TAI)之间存在一个整数秒的恒定差异，两者之间的秒内偏差(2σ)被控制在28ns以内。因为Galileo***时间不跳秒，所以它与UTC之间整数秒的差异随着UTC的跳秒而变化。可见，Galileo***时间与GPS时间在运行上极为相似。
Galileo***时间的原点定义为UT时间的1999年8月22日(星期日)零时零点，即从8月21日星期六午夜转变到8月22日星期日凌晨的零时零点，在这一原点起始时刻，Galileo***时间比UTC超前13s，而这整数秒差异值会逐渐随着UTC的跳秒变大。GPS时间在这Galileo***时间原点处刚好也比UTC超前13s，因而Galileo***时间与GPS时间之间的整数秒差异等于0。
Galileo***时间在导航电文中由32位数据比特给出，它由以下两部分组成：
1)高12比特的星期数(WN)是从Galileo***时间原点开始算起的星期计数值对4096的模，这相当于每约78年重新置零一次；
2)低20比特的周内时(TOW)是从每个星期交界处开始算起的、以秒为单位的计数值，其值范围0～604799s，在每星期结束(亦即开始)时重返至零。
Galileo导航电文给出了Galileo***时间与GPS时间之间实现转换所需的A0G、A1G、toG、WNoG四个参数，它们一起称为Galileo/GPS时间差异参数(GGTO)，其中，toG、WNoG指出这套参数的参考时间，A0G、A1G为关于GGTO的一线性模型中的两个系数。给定一个GPS时间tGPS及其所属的星期数WNGPS，可得此时的Galileo/GPS***时间偏差ΔtGAL,GPS为
ΔtGAL,GPS＝A0G+A1G(tGPS-toG+604800|WNGPS-WN0G|64)[s](1-3)
式中，偏差ΔtGAL,GPS定义为相同时刻的Galileo***时间tGAL减去GPS时间tGPS，即
tGPS＝tGAL-ΔtGAL,GPS (1-4)
而|i|64代表i对64的模。GGTO的准确度(2σ)预计为5ns，而ΔtGAL,GPS也可以作为一个未知量在Galileo/GPS联合定位、定时计算中。
GPS***与GLONASS***时间转换
我们可从GPS和GLONASS两方面导航电文中解调出这两个***时间的关系参数：一方面，GLONASS-M型卫星在其播发的每一帧导航电文的第5串中提供GLONASS时间tGLO相对于GPS时间tGPS的秒内偏差值τGPS；另一方面，GPS NAV电文提供tGPS与UTC之间由跳秒引入的整数秒差异值ΔtLS。这样，再加上考虑tGLO相对于UTC的3小时差异，可得如下tGLO和tGPS之间的转换关系式
tGLO＝(tGPS-ΔtLS)-τGPS+3×3600[s] (1-5)
GPS和GLONASS两***时间之间的秒内偏差值τGPS可以从GLONASS导航电文中解调获得。在较短的一个时间段内和对非精密定位而言，值τGPS可以被认为是保持不变或者说是可以预测的。
GPS***与GLONASS***坐标系转换
WGS-84和PZ-90两坐标系有着相同的定义，可是因为这两个坐标系是各自通过一些不同的观测基准站点加以实现的，而这些站点的位置坐标事实上决定着坐标系，所以这两个如此得以实现的坐标系最后也不可能完全一致。WGS-84与起先的PZ-90在地球表面的坐标差异可达20m。对GPS、GLONASS、Galileo和COMPASS四大GNSS而言，GNSS联合定位通常将不同GNSS的卫星位置坐标值转换到同一个坐标系中去，而WGS-84通常被选择作为此共同的坐标转换目标，因此只需要考虑WGS-84与PZ-90坐标系之间的差异及其变换。
假设空间一点P在一个直角坐标系S中的位置坐标为(x,y,z)，那么该点在另一个不同却几乎平行的直角坐标系S′中的位置坐标(x′,y′,z′)为
式中，坐标平移量(Δx,Δy,Δz)为坐标系S的原点在坐标系S′中的坐标，δε、δψ和δω为坐标系S分别绕X、Y和Z坐标轴旋转而能得到S′相应坐标轴的所需旋转角度，δs体现为坐标尺度因子。如式5-20所示的这一坐标变换模型称为赫尔模特(Helmert)变换。它一共包含着7个参数，其中绕Z轴旋转的角度参量δω对整个变换而言显得最为敏感与重要。因为δε、δψ和δω这三个坐标轴旋转角度值都很小，即都很接近于0，所以如式5-20所显示的那样，这些值都分别约等于他们相应的正弦值，而它们的余弦值都约等于1。我们一般总是将PZ-90坐标值通过这种赫尔模特变换成WGS-84坐标值，即式5-20中的等号右边(x,y,z)为一个PZ-90坐标值，而等号左边(x′,y′,z′)为相应的WGS-84坐标值。
虽然GLONASS接口控制文件和俄罗斯政府方面均没有给出一套官方式的从PZ-90到WGS-84的赫尔模特变换参数，但是多个不用的民间组织都曾奋力解决这一技术问题。很多文献讨论了关于从PZ-90到WGS-84的赫尔模特变换及其参数，表5-2给出了其中的四套不同的变换参数值，表中每一行对应于一套独立的变换参数值。顺便提提下，表中各个角度的单位“″”代表秒，而60″为1′，60′为1°。例如，第3套参数中的δω值为-0.4″，这相当于可直接用来代入式5.15中的δω值为-1.94×10-6rad。
表5-1几套不同的从PZ-90到WGS-84的赫尔模特变换参数值
从测量数据的时间跨度和空间跨度方面讲，表5-2中的第3套和第4套赫尔模特变换参数值相对可靠，他们所对应的坐标系的原点之间相差2.5m，在绕Z轴的旋转角度δω上相差0.07″。如果对地球表面赤道上的一点进行坐标变换，那么利用这两套变换参数所得到的WGS-84坐标值结果之间最多相差4.6m，而该差异量一般来说甚至小于GNSS单点定位的误差量；如果对处在2万千米高空的卫星进行坐标变换，那么根据这两套变换参数所得的WGS-84坐标值结果之间最多相差11m，大致相当于将GNSS卫星的星历误差量加倍，而这通常仍能满足一般性导航应用对GNSS单点定位的精度要求。我们可以这么说，对精度要求为米级的导航定位应用***来讲，表5-2中的任意一套变换参数均能胜任。
本算法采用第三套赫尔模特变换参数值。
总而言之，PZ-90与WGS-84两坐标系之间的差异很小，所以一些不严格、非精密的民用导航应用***可以对此不做任何坐标变换，相应的坐标***差异和坐标变换误差均可视为只是卫星星历误差的一部分，我们还需要特别指出的是，俄罗斯在2007年9月更新了GLONASS的地球参考框架，这使得更新后的PZ-90.02与ITRF基本一致，他们两者之间已被建议不做任何坐标变换。
将GLONASS卫星位置从PZ-90转换到WGS-84之后的坐标值一般指的是卫星在信号发射时刻的位置，因而在定位计算时，我们还需要再进行地球自转校正对流层，将GLONASS卫星位置坐标值最终转换到信号接收时刻的WGS-84坐标值。
采用多***融合共视时间比对算法最终可以的到不同接收机时间差数据，但由于观测噪声的影响，这组数据的精度很难达到要求。卡尔曼滤波是一个对信号或数据进行处理和变换的过程，其最主要目的是去掉或削弱不想要的成分对估计值的影响，并增强想要的成分的权重。并且卡尔曼滤波采用递推处理，用上个采样时刻的异地钟差估计值和当前时刻的异地钟差观测值，来估计当前时刻的异地钟差估计值，当前时刻以后的观测值不会对当前时刻的估计值产生任何影响，因而适合于实时的共视观测资料处理。因此对时间差数据进行卡尔曼滤波是提高时间比对精度的常用方法。对含噪声的钟差数据(观测量)进行Kalman滤波，可估计出准确的钟差。假设在k时刻钟差真值用xk表示，它构成了状态变量Xk，这里
Xk＝(xk)
Kalman滤波的状态方程为
Xk＝Φk,k-1Xk-1+Wk-1
Φk,k-1为状态转移矩阵，Wk-1为模型噪声。
Kalman滤波器的动态***维数n、观测***维数m均为1。
考虑接收机未能按照共视表规定时刻及时锁定卫星，甚至在共视表所规定的整个跟踪时间段内，一直未能锁定卫星，致使共视数据中缺少该记录。我们采用等间隔的Kalman滤波器，对于共视间歇和未能锁星成功的时间段，我们用前3个时刻点上钟差估计值进行二阶多项式外插，作为当前时刻的观测值x(k)，使Kalman滤波继续进行。二阶多项式外插考虑到了两站原子钟的频差和相对漂移的影响，可以对间歇时间段内两站钟差变化量进行较为准确的估计，从而保证了Kalman滤波器性能。
在共视资料预处理后，对含噪声的异地钟差数据序列进行Kalman滤波，其算法过程为：
第一步，状态变量Xk(只含钟差真值一个分量)与其Kalman估计之间的均方误差矩阵称为估计误差的协方差矩阵，用Ck表示。给定C0一个初值，根据下式
可以计算出P1。其中，Pk为状态变量Xk与其在无观测噪声与模型噪声条件下的估计之间的均方误差阵，Qk为1×1阶模型噪声叭的协方差阵。
第二步，得到P1后，根据Kalman增益矩阵Gk的表示式
求得G1，其中Rk为1×1阶观测噪声Vk的协方差阵。
第三步，根据下式
可得k＝1时刻的状态变量估计值(即k＝1时刻钟差Kalman估计值)。
第四步，将P1带入下式
Ck＝(I-GkHk)Pk
可求得k＝1时刻的1×1阶估计误差的协方差阵C1。然后，进入下次循环。
在近实时共视中，应该对异地钟差X0有初步的估计，用这个值对滤波器进行初始化，会使滤波器的收敛速度加快。此时，估计误差的协方差矩阵初值取为
RTS固定区间最优平滑算法是在卡尔曼滤波的基础上，利用整个时间间隔内所有量测数据得到状态的最小方差估计，可以获得比卡尔曼滤波精度更高的融合结果。平滑解算过程相对于滤波过程是逆向的。因此，RTS固定区间平滑算法在传递对准精度评估等侧重于初始状态获取的应用中，最终平滑值的读取方式与普通前向滤波器估计值的读取方式相反。平滑过程首先是进行前向卡尔曼滤波，获得滤波估计值，然后经过一个反向平滑过程，进而的到平滑估计值。因此，平滑解算需要在滤波过程中实时存储数据，所存储的数据为4个矩阵，分别为估计值一步转移阵估计的均方差阵以及一步预测均方差阵
平滑公式为：
其中
式中，为RTS平滑值；Ks,k为滤波增益；Ps,k为协方差阵。
经过RTS处理后的时间差数据即为最终数据。

Claims (9)

1.一种GNSS时间频率传递接收机，包括天线单元、时频单元、射频单元、信号处理单元和监控单元，其特征在于：时频单元产生射频单元和信号处理单元所需要的时频信号，同时输出PPS信号；射频单元把经过天线单元处理的射频信号进行分离，经下变频及低噪声放大、AGC控制后，为信号处理单元提供中频信号；信号接收处理单元完成导航信号的采样、捕获、跟踪、导航电文解调、伪距测量并上报各种观测量；监控单元对观测数据进行传输和对各单元下发控制指令。
2.根据权利要求1所述的一种GNSS时间频率传递接收机，其特征在于，信号处理单元包括信号捕获模块、通道跟踪模块和信息处理模块；信号捕获模块接收用户配置参数；跟踪通道组包含伪码发生器、载波NCO、码NCO和累加器；信息处理模块的工作过程包含以下步骤：
步骤一 分别通过GPS定位***、COMPASS定位***、Galileo定位***和GLONASS定位***观测各自***内的卫星，各获得一组时间比对数据；
步骤二 统一时间***，包括GPS定位***与COMPASS定位***时间转换、GPS定位***与Galileo定位***时间转换、GPS定位***与GLONASS定位***时间转换；
步骤三 对步骤二得到的统一后的时间数据进行解码和预处理获得GNSS卫星轨道位置，并对获得的GNSS卫星位置坐标进行地球自转效应修正；
步骤四 统一空间***；
步骤五 根据统一的空间***，计算GNSS卫星与两个GNSS接收机的几何距离；计算卫星与接收机然后计算信号传播过程中的对流层时延；将这些参数代入到由伪距观测方程推导出的地面接收机时间差观测方程中，进而计算出时间差值；
步骤六 对时间比对结果进行Kalman滤波和RTS事后处理。
3.根据权利要求2所述的一种GNSS时间频率传递接收机，其特征在于，步骤二中所述的GPS定位***与COMPASS定位***时间转换过程为：
COMPASS定位***和GPS定位***时间的换算公式为：
tGPS＝tCOM-ΔtCOM,GPS
ΔtCOM,GPS＝A0GPS+A1GPS×tCOM
式中，tCOM为指定用户计算的COMPASS时间，A0GPS是COMPASS相对于GPS定位***时间的钟差，A1GPS是COMPASS定位***相对于GPS定位***时间的钟速，取周内秒计数部分；
步骤二中所述的GPS定位***与Galileo定位***时间转换过程为：
Galileo导航电文给出Galileo***时间与GPS时间之间实现转换所需的A0G、A1G、toG、WNoG四个参数，其中，toG、WNoG指出这套参数的参考时间，A0G、A1G为关于GGTO的一线性模型中的两个系数；给定一个GPS时间tGPS及其所属的星期数WNGPS，得此时的Galileo/GPS***时间偏差ΔtGAL,GPS为：
ΔtGAL,GPS＝A0G+A1G(tGPS-toG+604800|WNGPS-WN0G|64)[s]
式中，偏差ΔtGAL,GPS定义为相同时刻的Galileo***时间tGAL减去GPS时间tGPS，即
tGPS＝tGAL-ΔtGAL,GPS
|i|64代表i对64的模；GGTO的准确度2σ预计为5ns，ΔtGAL,GPS作为一个未知量在Galileo/GPS联合定位、定时计算中；
步骤二中所述的GPS定位***与GLONASS定位***时间转换过程为：
GLONASS-M型卫星在其播发的每一帧导航电文的第5串中提供GLONASS时间tGLO相对于GPS时间tGPS的秒内偏差值τGPS；GPS NAV电文提供tGPS与UTC之间由跳秒引入的整数秒差异值ΔtLS；再加上考虑tGLO相对于UTC的3小时差异，得如下tGLO和tGPS之间的转换关系式：
tGLO＝(tGPS-ΔtLS)-τGPS+3×3600[s]
GPS和GLONASS两***时间之间的秒内偏差值τGPS从GLONASS导航电文中解调获得。
4.根据权利要求2所述的一种GNSS时间频率传递接收机，其特征在于，步骤四中所述的统一空间***为：
GPS定位***与GLONASS定位***坐标系转换：
设空间一点P在一个直角坐标系S中的位置坐标为(x,y,z)，则该点在另一个不同却几乎平行的直角坐标系S′中的位置坐标(x′,y′,z′)为：
式中，坐标平移量(Δx,Δy,Δz)为坐标系S的原点在坐标系S′中的坐标，δε、δψ和δω为坐标系S分别绕X、Y和Z坐标轴旋转而能得到S′相应坐标轴的所需旋转角度，δs体现为坐标尺度因子；式中等号右边(x,y,z)为一个PZ-90坐标值，等号左边(x′,y′,z′)为相应的WGS-84坐标值。
5.根据权利要求2所述的一种GNSS时间频率传递接收机，其特征在于，步骤六包含以下步骤：
步骤一状态变量Xk与其Kalman估计之间的均方误差矩阵称为估计误差的协方差矩阵，用Ck表示；给定C0一个初值，根据下式
计算出P1；其中，Pk为状态变量Xk与其在无观测噪声与模型噪声条件下的估计之间的均方误差阵，Qk为1×1阶模型噪声叭的协方差阵；
步骤二得到P1后，根据Kalman增益矩阵Gk的表示式
求得G1，其中Rk为1×1阶观测噪声Vk的协方差阵；
步骤三根据下式
得k＝1时刻的状态变量估计值即k＝1时刻钟差Kalman估计值；
步骤四将P1带入下式
Ck＝(I-GkHk)Pk
求得k＝1时刻的1×1阶估计误差的协方差阵C1；然后，进入下次循环；
在近实时共视中，对异地钟差X0有初步的估计，用这个值对滤波器进行初始化，使滤波器的收敛速度加快；此时，估计误差的协方差矩阵初值取为：
步骤五在卡尔曼滤波的基础上，利用整个时间间隔内所有量测数据得到状态的最小方差估计，对结果进行RTS固定区间最优平滑，平滑解算需要在滤波过程中实时存储数据，所存储的数据为4个矩阵，分别为估计值一步转移阵估计的均方差阵以及一步预测均方差阵
平滑公式为：
K＝N-1，N-2，……2，1，0；
其中
式中，为RTS平滑值；Ks,k为滤波增益；Ps,k为协方差阵。
6.一种GNSS时间频率传递接收方法，其特征在于，包含以下步骤：
步骤一 分别通过GPS定位***、COMPASS定位***、Galileo定位***和GLONASS定位***观测各自***内的卫星，各获得一组时间比对数据；
步骤二 统一时间***，包括GPS定位***与COMPASS定位***时间转换、GPS定位***与Galileo定位***时间转换、GPS定位***与GLONASS定位***时间转换；
步骤三 对步骤二得到的统一后的时间数据进行解码和预处理获得GNSS卫星轨道位置，并对获得的GNSS卫星位置坐标进行地球自转效应修正；
步骤四 统一空间***；
步骤五 根据统一的空间***，计算GNSS卫星与两个GNSS接收机的几何距离；计算卫星与接收机然后计算信号传播过程中的对流层时延；将这些参数代入到由伪距观测方程推导出的地面接收机时间差观测方程中，进而计算出时间差值；
步骤六 对时间比对结果进行Kalman滤波和RTS事后处理。
7.根据权利要求6所述的一种GNSS时间频率传递接收方法，其特征在于，步骤二中所述的GPS定位***与COMPASS定位***时间转换过程为：
COMPASS定位***和GPS定位***时间的换算公式为：
tGPS＝tCOM-ΔtCOM,GPS
ΔtCOM,GPS＝A0GPS+A1GPS×tCOM
式中，tCOM为指定用户计算的COMPASS时间，A0GPS是COMPASS相对于GPS定位***时间的钟差，A1GPS是COMPASS定位***相对于GPS定位***时间的钟速，取周内秒计数部分；
步骤二中所述的GPS定位***与Galileo定位***时间转换过程为：
Galileo导航电文给出Galileo***时间与GPS时间之间实现转换所需的A0G、A1G、toG、WNoG四个参数，其中，toG、WNoG指出这套参数的参考时间，A0G、A1G为关于GGTO的一线性模型中的两个系数；给定一个GPS时间tGPS及其所属的星期数WNGPS，得此时的Galileo/GPS***时间偏差ΔtGAL,GPS为：
ΔtGAL,GPS＝A0G+A1G(tGPS-toG+604800|WNGPS-WN0G|64)[s]
式中，偏差ΔtGAL,GPS定义为相同时刻的Galileo***时间tGAL减去GPS时间tGPS，即
tGPS＝tGAL-ΔtGAL,GPS
|i|64代表i对64的模；GGTO的准确度2σ预计为5ns，ΔtGAL,GPS作为一个未知量在Galileo/GPS联合定位、定时计算中；
步骤二中所述的GPS定位***与GLONASS定位***时间转换过程为：
GLONASS-M型卫星在其播发的每一帧导航电文的第5串中提供GLONASS时间tGLO相对于GPS时间tGPS的秒内偏差值τGPS；GPS NAV电文提供tGPS与UTC之间由跳秒引入的整数秒差异值ΔtLS；再加上考虑tGLO相对于UTC的3小时差异，得如下tGLO和tGPS之间的转换关系式：
tGLO＝(tGPS-ΔtLS)-τGPS+3×3600[s]
GPS和GLONASS两***时间之间的秒内偏差值τGPS从GLONASS导航电文中解调获得。
8.根据权利要求6所述的一种GNSS时间频率传递接收方法，其特征在于，步骤四中所述的统一空间***为：
GPS定位***与GLONASS定位***坐标系转换：
设空间一点P在一个直角坐标系S中的位置坐标为(x,y,z)，则该点在另一个不同却几乎平行的直角坐标系S′中的位置坐标(x′,y′,z′)为：
式中，坐标平移量(Δx,Δy,Δz)为坐标系S的原点在坐标系S′中的坐标，δε、δψ和δω为坐标系S分别绕X、Y和Z坐标轴旋转而能得到S′相应坐标轴的所需旋转角度，δs体现为坐标尺度因子；式中等号右边(x,y,z)为一个PZ-90坐标值，等号左边(x′,y′,z′)为相应的WGS-84坐标值。
9.根据权利要求6所述的一种GNSS时间频率传递接收方法，其特征在于，步骤六包含以下步骤：
步骤一 状态变量Xk与其Kalman估计之间的均方误差矩阵称为估计误差的协方差矩阵，用Ck表示；给定C0一个初值，根据下式
计算出P1；其中，Pk为状态变量Xk与其在无观测噪声与模型噪声条件下的估计之间的均方误差阵，Qk为1×1阶模型噪声叭的协方差阵；
步骤二 得到P1后，根据Kalman增益矩阵Gk的表示式
求得G1，其中Rk为1×1阶观测噪声Vk的协方差阵；
步骤三 根据下式
得k＝1时刻的状态变量估计值即k＝1时刻钟差Kalman估计值；
步骤四 将P1带入下式
Ck＝(I-GkHk)Pk
求得k＝1时刻的1×1阶估计误差的协方差阵C1；然后，进入下次循环；
在近实时共视中，对异地钟差X0有初步的估计，用这个值对滤波器进行初始化，使滤波器的收敛速度加快；此时，估计误差的协方差矩阵初值取为：
步骤五 在卡尔曼滤波的基础上，利用整个时间间隔内所有量测数据得到状态的最小方差估计，对结果进行RTS固定区间最优平滑，平滑解算需要在滤波过程中实时存储数据，所存储的数据为4个矩阵，分别为估计值一步转移阵估计的均方差阵以及一步预测均方差阵
平滑公式为：
K＝N-1，N-2，……2，1，0；
其中
式中，为RTS平滑值；Ks,k为滤波增益；Ps,k为协方差阵。
CN201711392705.2A 2017-12-21 2017-12-21 一种gnss时间频率传递接收机及接收方法 Pending CN108181633A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711392705.2A CN108181633A (zh) 2017-12-21 2017-12-21 一种gnss时间频率传递接收机及接收方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711392705.2A CN108181633A (zh) 2017-12-21 2017-12-21 一种gnss时间频率传递接收机及接收方法

Publications (1)

Publication Number Publication Date
CN108181633A true CN108181633A (zh) 2018-06-19

Family

ID=62546862

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711392705.2A Pending CN108181633A (zh) 2017-12-21 2017-12-21 一种gnss时间频率传递接收机及接收方法

Country Status (1)

Country Link
CN (1) CN108181633A (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109634092A (zh) * 2019-01-16 2019-04-16 和芯星通科技(北京)有限公司 一种基于gnss接收机的授时方法及gnss接收机
CN109631883A (zh) * 2018-12-17 2019-04-16 西安理工大学 一种基于节点信息共享的载机局部姿态精确估计方法
CN111025359A (zh) * 2019-12-30 2020-04-17 中国人民解放军陆军工程大学 一种用于卫星导航的跳频测距方法和***
CN111141279A (zh) * 2019-12-20 2020-05-12 北京小马慧行科技有限公司 行车轨迹的处理方法及装置
CN112083462A (zh) * 2020-09-14 2020-12-15 哈尔滨工程大学 一种海洋精密单点定位方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102023290A (zh) * 2010-11-04 2011-04-20 中国民用航空总局第二研究所 高精度分布式脉冲信号到达时间差检测***
CN105044747A (zh) * 2015-08-27 2015-11-11 交通信息通信技术研究发展中心 一种基于多星共视和滤波的时间同步装置及其方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102023290A (zh) * 2010-11-04 2011-04-20 中国民用航空总局第二研究所 高精度分布式脉冲信号到达时间差检测***
CN105044747A (zh) * 2015-08-27 2015-11-11 交通信息通信技术研究发展中心 一种基于多星共视和滤波的时间同步装置及其方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
DEFRAIGNE PASCALE 等: "Time and frequency transfer combining GLONASS and GPS data", 《2011 JOINT CONFERENCE OF THE IEEE INTERNATIONAL FREQUENCY CONTROL SYMPOSIUM/EUROPEAN FREQUENCY AND TIME FORUM PROCEEDINGS》 *
XIANDE WU 等: "Accuracy analysis of a combined Common-view time Comparison algorithm of GNSS navigation system based on Kalman filtering and RTS filtering", 《2017 CHINESE AUTOMATION CONGRESS (CAC)》 *
ZHANG VS 等: "Multi-channel GPS/GLONASS common-view between NIST and USNO", 《PROCEEDINGS OF THE 2000 IEEE/EIA INTERNATIONAL FREQUENCY CONTROL SYMPOSIUM & EXHIBITION》 *
戚素娟: "GPS/GLONASS组合共视技术研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
李征航 等: "《GPS测量与数据处理 第三版》", 31 May 2016, 武汉大学出版社 *
赵树强 等: "《箭载GNSS测量数据处理》", 31 May 2015, 国防工业出版社 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109631883A (zh) * 2018-12-17 2019-04-16 西安理工大学 一种基于节点信息共享的载机局部姿态精确估计方法
CN109631883B (zh) * 2018-12-17 2022-12-09 西安理工大学 一种基于节点信息共享的载机局部姿态精确估计方法
CN109634092A (zh) * 2019-01-16 2019-04-16 和芯星通科技(北京)有限公司 一种基于gnss接收机的授时方法及gnss接收机
CN111141279A (zh) * 2019-12-20 2020-05-12 北京小马慧行科技有限公司 行车轨迹的处理方法及装置
CN111141279B (zh) * 2019-12-20 2022-07-01 北京小马慧行科技有限公司 行车轨迹的处理方法及装置
CN111025359A (zh) * 2019-12-30 2020-04-17 中国人民解放军陆军工程大学 一种用于卫星导航的跳频测距方法和***
CN112083462A (zh) * 2020-09-14 2020-12-15 哈尔滨工程大学 一种海洋精密单点定位方法
CN112083462B (zh) * 2020-09-14 2022-12-13 哈尔滨工程大学 一种海洋精密单点定位方法

Similar Documents

Publication Publication Date Title
CN108181633A (zh) 一种gnss时间频率传递接收机及接收方法
CN110187364B (zh) 一种低轨导航增强精密改正数据生成、上注***及方法
Luthcke et al. The 1-centimeter orbit: Jason-1 precision orbit determination using GPS, SLR, DORIS, and altimeter data special issue: Jason-1 calibration/validation
AU2012205602B2 (en) Method and system for determining clock corrections
Arikan et al. Regularized estimation of vertical total electron content from GPS data for a desired time period
Bertiger et al. Sub-centimeter precision orbit determination with GPS for ocean altimetry
CN103728647B (zh) 一种基于卫星载波信号调制的弹体滚转角测量方法
CN103033188B (zh) 基于综合孔径观测的导航卫星自主时间同步方法
CN102033236B (zh) 一种卫星导航位置速度联合估计方法
CN107607971B (zh) 基于gnss共视时间比对算法的时间频率传递方法及接收机
CN111596321B (zh) 利用非差改正数的多gnss多路径误差恒星日滤波方法及***
KR20160147775A (ko) 대기 관측을 위한 무선항법 신호 처리를 위한 방법 및 장치
CN106850177A (zh) 一种基于精密单点定位的实时时间传递方法
Capuano et al. Orbital filter aiding of a high sensitivity GPS receiver for lunar missions
CN114994727B (zh) 一种实现高精度时间校准及卫星定位的设备
CN112748449A (zh) 一种卫星导航接收机锁相环与锁频环结合的矢量跟踪方法
CN112230254B (zh) 一种gps载波相位多径误差的校正方法及装置
Natras et al. Regional ionosphere delay models based on CORS data and machine learning
CN111812681B (zh) 大气区域建模方法、装置、电子设备及存储介质
CN110988932B (zh) 一种提高星载gps接收机实时钟差解算精度的方法
CN107643527B (zh) Gps共视信号模拟生成方法及gps共视信号模拟器
CN116699671A (zh) 一种基于随机森林回归的电离层幅度闪烁指数计算方法
CN108051834A (zh) 一种包含gnss共视时间比对算法的时间频率传递接收机
Li et al. Precise orbit determination for LEO satellites: single-receiver ambiguity resolution using GREAT products
CN105510936A (zh) 星载gnss联合定轨方法

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20180619