[发明专利]一种基于北斗三号PPP-B2b信息的连续时间传递方法

权利要求书

1.一种基于北斗三号PPP-B2b信息的连续时间传递方法，其特征在于，包括如下步骤：

S1：若观测数据连续，利用北斗三号PPP-B2b信息进行PPP解算，获得时间传递结果；

所述PPP时间传递采用双频消电离层组合，数学模型如下：

式中：j表示卫星号，Φ和P表示载波相位和伪距消电离层组合观测量，ρ为卫星到接收机的几何距离，c为光速，δt^j卫星钟差，δt为接收机钟相对于北斗时的时间传递结果，λ、N分别为消电离层组合载波长和模糊度，T是对流层延迟，为载波和伪距测量噪声；

S2：若观测数据异常而造成观测数据瞬时中断，利用接收机历史钟差数据预报的钟差结果，作为当前历元的时间传递结果；

当观测数据恢复时，将预报的当前历元时间传递结果作为约束条件，来约束待估计的钟差参数，辅助PPP时间传递算法，加快数据的重新收敛进行解算，进而实现连续的时间传递。

2.根据权利要求1所述的一种基于北斗三号PPP-B2b信息的连续时间传递方法，其特征在于：所述步骤S1包含如下步骤：

S11.获取观测数据、广播星历和PPP-B2b信息，并利用PPP-B2b信号提供的实时卫星轨道、钟差和码间偏差的改正信息对卫星轨道、卫星钟差和卫星码间偏差进行改正，得到最终的卫星轨道、钟差结果和修正后的观测值；

S12.数据预处理是剔除观测数据中的粗差数据，并探测周跳，探测周跳采用TurboEdit和电离层残差法，对周跳进行标记；

S13.数据处理是对固体潮、海潮、相位缠绕、对流层延迟、地球自转、相对论效应和接收机天线相位中心的误差进行改正；

S14.参数估计是求解接收机位置、接收机钟差、对流层天顶延迟湿分量、消电离层组合模糊度的未知参数；

S15.对卡尔曼滤波输出的残差结果进行检验，如果最大残差大于阈值，则判断为粗差，剔除后重新进行卡尔曼滤波解算，直到残差检验通过；

S16.输出并存储解算得到的时间传递结果。

3.根据权利要求1所述的一种基于北斗三号PPP-B2b信息的连续时间传递方法，其特征在于：所述步骤S2包含如下步骤：

S21.在进行时间传递解算前，判断观测数据是否连续；

S22.如果连续则采用步骤S1的数据处理方法；

S23.如果数据仅当前历元异常中断，进行标记，并利用存储的接收机钟差数据进行建模，建立接收机钟差模型，对当前历元钟差结果进行预报；

S24.若之前观测数据中断，当前历元观测数据恢复正常，则可在原始观测方程基础上，根据步骤S22利用建立的接收机钟差模型预报当前历元的时间传递结果，将预报的当前历元时间传递结果作为约束条件，来约束待估计的钟差参数，辅助PPP时间传递算法，加快数据的重新收敛进行解算，实现连续的时间传递。

4.根据权利要求3所述的一种基于北斗三号PPP-B2b信息的连续时间传递方法，其特征在于：所述步骤S23中当前历元钟差结果进行预报的步骤如下：

a.建立接收机钟差模型：

x(t)＝a₀+a₁(t-t₀)+a₂(t-t₀)²+ε_x(t)

式中，a₀、a₁、a₂表示钟差、钟速、钟漂，ε_x(t)为随机变化的不确定性分量；

b.建立钟差模型矩阵：

X＝Ha+e

c.引入遗忘因子矩阵W＝diag(λ^n-1,λ^n-2,...,λ⁰),0＜λ≤1，并计算和T_m(m≥3)

b)迭代计算k(k＝m+1,m+2,…,n)时刻的T_k：

c)预测当前历元t_k时刻的钟差，也即是时间传递结果：

采用递推遗忘因子最小二乘算法进行预报，将预报得到的时间传递结果存储起来，并可以将当前预报结果输出，以确保时间传递结果的连续性。