一种基于GNSS的高精度时间传递方法技术

本发明专利技术涉及一种基于GNSS的高精度时间传递方法，该方法将GNSS差分系统基站间观测值进行共视选星和数据预处理，通过载波相位观测值和伪距观测值组合得到观测矩阵，将载波相位单差模糊度和时间传递结果作为未知量，组成观测方程组。使用序贯最小二乘递归求得组合方程单差载波的模糊度，进而求解准确的钟差。更进一步地，当已知接收机时钟模型为线性漂移时，通过kalman滤波的方式求解相应未知量，模型更加准确，得到的钟差精度也更高。该方法能够同时实现低成本高精度等要求，能够达到纳秒级以下精度，是一种能够广泛使用的时间传递方法。

【技术实现步骤摘要】
一种基于GNSS的高精度时间传递方法
本专利技术属于时间同步
，具体涉及一种基于GNSS的高精度时间传递方法。
技术介绍
目前多个领域对高精度时间和频率传递的要求增加，例如电力系统、无线定位系统、分布式雷达、协同数据链、海上舰艇编队等领域。这些系统中，除了对时间传递的精度有着较高要求外，在应用于军事、空天等领域对设备的可靠性和适应性也有很高的要求。全球卫星导航系统(GNSS)的全天候、全覆盖、高精度等特点，使其成为了高精度时间传递中的重要方法。目前常用的时间传递方法包括卫星双向时间传递、卫星共视法、精密单点定位等技术，其中卫星双向时间传递方法采用卫星双向链路进行时间比对，是目前GNSS时间传递中精度最高的方法，但是其需要使用转发式卫星，应用条件十分受限。GNSS共视法通过消除卫星信号传递中的共同误差提高时间传递精度，但是其采用伪距作为观测值，精度不够高。PPP借助国际GNSS服务组织发布的精密星历和钟差产品可以达到亚纳秒级的时间传递精度，由于其性能不受距离的影响，应用广泛。但是该方法依赖精密星历和钟差产品，需要长时间观测，不能实时计算，只能用于事后分析。这些时间传递方法或多或少存在算法复杂、成本较高、使用场景受限，实时性较差等缺点，不能兼顾低成本、实时性与高精度。因此，急需一种低成本高精度适用范围广的实时时间传递方法。
技术实现思路
为解决所述问题，本专利技术提供了一种基于GNSS的时间传递方法，该方法能够同时实现低成本高精度等要求，能够达到纳秒级以下精度，是一种能够广泛使用的时间传递方法。本专利技术一种基于GNSS的高精度时间传递方法，包括如下步骤：步骤1，对GNSS差分系统中的移动站和基准站两个基站进行共视卫星选星，然后进行观测量数据预处理；所述观测量为基站到共视卫星的伪距和载波相位观测量；所述数据预处理为针对原始观测信息，对所述伪距观测量的粗差和所述载波相位观测量的周跳进行修正，以及对卫星信号传播过程中的误差进行修正；步骤2，将所述伪距和载波观测量进行基站间差分，组合成当前时刻观测方程：其中，和为伪距和载波的单差观测向量，σr为两者的噪声权重比，t为接收机相对钟差，τ＝c·t为将接收机钟差化为单位为米的量，为卫星到两个基站的几何距离差值，为各个共视卫星的载波相位单差整周模糊度，μ，ν为相应的码和载波观测噪声，λ为载波波长；步骤3，对所述时间传递过程中的每个历元建立观测方程，并将当前历元的观测方程与上个历元的观测方程进行联立建立历元间的观测方程；通过矩阵变换对所述历元间的观测方程进行降维运算，使其维数和未知量维数一致；步骤4，通过所述观测方程在历元间的迭代运算，计算出所述GNSS差分系统基站间的接收机钟差，得到时间传递结果；将所述时间传递结果补偿到移动站接收机钟差中，修正移动站的本地时间和输出的秒脉冲信号。进一步地，所述步骤3中，当所述共视卫星出现卫星升起或者降落的时候，对所述观测方程依据升降星进行矩阵变换，使方程迭代过程中保持矩阵可联立，具体包括以下步骤：步骤3.1，将所述当前时刻的共视卫星和上一时刻的共视卫星进行比较，分成上一时刻落下的卫星、保持不变的卫星和当前时刻升起的卫星，并将其变换位置在观测方程中通过矩阵的形式表现出来，完成迭代过程中的卫星置换；步骤3.2，将所述当前历元和上个历元的观测方程进行联立，对所述联立观测方程系数矩阵进行QR分解得到矩阵QT和上三角矩阵R；将矩阵QT分别与联立观测方程左式和系数矩阵相乘，得到降维后的观测方程，以此方法进行多历元观测方程的迭代；步骤3.3，通过所述迭代过程得到当前时刻k的时间传递结果τk。进一步地，接收机时钟为线性漂移模型时，将钟漂作为未知量加入到所述未知量矩阵中，所述观测方程使用钟差和钟漂估计,并使用kalman滤波模型进行未知量解算。本专利技术有益效果在于相比于卫星双向时间传递算法、PPP等传统时间传递方法，能够不依赖转发式卫星和精密星历等条件，实现高精度实时时间传递，同时能够应用于大多数场景中，适应性强可靠性高成本较低，能够被广泛采用。附图说明图1为本专利技术的基于GNSS的高精度时间传递原理图；图2为本专利技术中所述接收机钟差模型图；图3为本专利技术中GNSS时间传递方法系统框图。具体实施方式下面参照附图，对专利技术优选实施例详细进行描述。本专利技术提供了一种实时解算高精度时间传递结果的方法，利用序贯最小二乘实时解算模糊度和钟差。由于单差观测值之间不存在强相关性，序贯最小二乘方法将解算模糊度和钟差通过历元之间递归求解，初始化时间较短，解算准确率高，满足实时性要求，能达到实时时间传递亚纳秒级精度。为了进一步提高时间传递的精度，在时钟模型为线性漂移的情况下引入了Kalman滤波的方法，在序贯最小二乘基础上增加了钟漂估计，提高了解算的精度和准确度。基于GNSS的高精度时间传递原理图如图1所示。其中伪距和载波的测量值都包含了时钟偏差和其他不可忽略的误差，做站间单差后卫星端的误差基本消除，大气延迟等误差也可以通过模型削弱或消除，因此在去除误差项后差分测量值主要包含了接收机钟差和几何距离差值。对星站间观测量的精确估计得到两站除几何距离和误差外的差值即时间传递的结果。时间传递的精度和可靠性除了受上述误差的影响，还受接收机自身时钟的性能的影响，由于其他几种误差经过预处理以后已经基本消除，因此本专利技术关注如何针对接收机钟差模型进一步提高时间传递精度。接收机时钟频标由接收机内部的石英振荡器提供，通常不同型号的接收机内部运转的方式有三种：(a)线性漂移。(b)实时调整。(c)将钟差控制在1ms以内。图2为常见接收机钟差模型图。由于序贯最小二乘方法不受钟差模型的影响，同时可以满足实时性的要求，在接收机时钟模型未知的情况下，我们采用序贯最小二乘的方法进行时间传递的实时求解，并且取得可亚纳秒级的精度。当已知钟差模型为(a)时，在序贯最小二乘求解模糊度和钟差的基础上，可以采用kalman滤波算法对时钟钟漂进行估计，该方法由于充分考虑了模型时钟变化特点，精度较序贯最小二乘方法有进一步的提升。图3给出了高精度实时时间传递的系统框图。本专利技术基于GNSS的高精度时间传递方法包括如下步骤：步骤1，对GNSS差分系统中的移动站和基准站两个基站进行共视卫星选星，然后进行观测量数据预处理；观测量为基站到共视卫星的伪距和载波相位观测量；数据预处理为针对原始观测信息，对伪距观测量的粗差和载波相位观测量的周跳进行修正，以及对卫星信号传播过程中的误差进行修正。步骤2，针对GNSS差分系统基站观测到的卫星，将伪距和载波观测量进行基站间差分，组合成当前时刻观测方程：其中，和为伪距和载波的单差观测向量，σr为两者的噪声权重比，t为接收机相对钟差，τ＝c·t为将接收机钟差化为单位为米的量，为卫星到两个基站的几何距离差值，为各个共视卫星的载波相位单差整周模糊度，μ，v为相应的码和载波观测噪声，λ为载波波长；...

【技术保护点】
1.一种基于GNSS的高精度时间传递方法，其特征在于，包括如下步骤：/n步骤1，对GNSS差分系统中的移动站和基准站两个基站进行共视卫星选星，然后进行观测量数据预处理；所述观测量为基站到共视卫星的伪距和载波相位观测量；所述数据预处理为针对原始观测信息，对所述伪距观测量的粗差和所述载波相位观测量的周跳进行修正，以及对卫星信号传播过程中的误差进行修正；/n步骤2，将所述伪距和载波观测量进行基站间差分，组合成当前时刻观测方程：/n

【技术特征摘要】
1.一种基于GNSS的高精度时间传递方法，其特征在于，包括如下步骤：
步骤1，对GNSS差分系统中的移动站和基准站两个基站进行共视卫星选星，然后进行观测量数据预处理；所述观测量为基站到共视卫星的伪距和载波相位观测量；所述数据预处理为针对原始观测信息，对所述伪距观测量的粗差和所述载波相位观测量的周跳进行修正，以及对卫星信号传播过程中的误差进行修正；
步骤2，将所述伪距和载波观测量进行基站间差分，组合成当前时刻观测方程：



其中，和为伪距和载波的单差观测向量，σr为两者的噪声权重比，t为接收机相对钟差，τ＝c·t为将接收机钟差化为单位为米的量，为卫星到两个基站的几何距离差值，为各个共视卫星的载波相位单差整周模糊度，μ，ν为相应的码和载波观测噪声，λ为载波波长；
步骤3，对所述时间传递过程中的每个历元建立观测方程，并将当前历元的观测方程与上个历元的观测方程进行联立建立历元间的观测方程；通过矩阵变换对所述历元间的观测方程进行降维运算，使其维数和未知量维数一致；
步骤4，通过所述观测方程在历元间的迭代运算，计算出所述GNSS差分系统基站间的接收机钟差，得到时间传递结果；将所述时间传递结果补偿到移...

【专利技术属性】
技术研发人员：张航，
申请(专利权)人：北京自动化控制设备研究所，
类型：发明
国别省市：北京;11