[发明专利]一种基于GNSS的高精度时间传递方法

权利要求书

1.一种基于GNSS的高精度时间传递方法，其特征在于，包括如下步骤：

步骤1，对GNSS差分系统中的移动站和基准站两个基站进行共视卫星选星，然后进行观测量数据预处理；所述观测量为基站到共视卫星的伪距和载波相位观测量；所述数据预处理为针对原始观测信息，对所述伪距观测量的粗差和所述载波相位观测量的周跳进行修正，以及对卫星信号传播过程中的误差进行修正；

步骤2，将所述伪距和载波观测量进行基站间差分，组合成当前时刻观测方程：

其中，和为伪距和载波的单差观测向量，σ_r为两者的噪声权重比，t为接收机相对钟差，τ＝c·t为将接收机钟差化为单位为米的量，为卫星到两个基站的几何距离差值，为各个共视卫星的载波相位单差整周模糊度，μ，ν为相应的码和载波观测噪声，λ为载波波长；

步骤3，对所述时间传递过程中的每个历元建立观测方程，并将当前历元的观测方程与上个历元的观测方程进行联立建立历元间的观测方程；通过矩阵变换对所述历元间的观测方程进行降维运算，使其维数和未知量维数一致；

步骤4，通过所述观测方程在历元间的迭代运算，计算出所述GNSS差分系统基站间的接收机钟差，得到时间传递结果；将所述时间传递结果补偿到移动站接收机钟差中，修正移动站的本地时间和输出的秒脉冲信号。

2.如权利要求1所述的一种基于GNSS的高精度时间传递方法，其特征在于，所述步骤3中，当所述共视卫星出现卫星升起或者降落的时候，对所述观测方程依据升降星进行矩阵变换，使方程迭代过程中保持矩阵可联立，具体包括以下步骤：

步骤3.1，将所述当前时刻的共视卫星和上一时刻的共视卫星进行比较，分成上一时刻落下的卫星、保持不变的卫星和当前时刻升起的卫星，并将其变换位置在观测方程中通过矩阵的形式表现出来，完成迭代过程中的卫星置换；

步骤3.2，将所述当前历元和上个历元的观测方程进行联立，对所述联立观测方程系数矩阵进行QR分解得到矩阵Q^T和上三角矩阵R；将矩阵Q^T分别与联立观测方程左式和系数矩阵相乘，得到降维后的观测方程，以此方法进行多历元观测方程的迭代；

步骤3.3，通过所述迭代过程得到当前时刻k的时间传递结果τ_k。

3.如权利要求1或2所述的一种基于GNSS的高精度时间传递方法，其特征在于，接收机时钟为线性漂移模型时，将钟漂作为未知量加入到所述未知量矩阵中，所述观测方程使用钟差和钟漂估计,并使用kalman滤波模型进行未知量解算。