[发明专利]面向实时应用的GPS频间钟差新预报方法

权利要求书

1.面向实时应用的GPS频间钟差新预报方法，其特征在于，包括以下步骤：

步骤一、获取当前时间弧段的GPS三频载波相位观测值，并估计得到该时间弧段t时刻的IFCB值，即IFCB时间序列，其中，仅对初始日期的初始时刻的IFCB值置零，初始时刻以后的IFCB值估计均不置零；

步骤二、建立IFCB预报模型，表达式如式(1)所示：

其中，IFCB(t)表示步骤一中估计得到的t时刻的IFCB值，m、h均为线性趋势项系数，n为高阶谐波函数的阶数，λ_i为周期项振幅，θ_i为周期项相位，T_i为时间周期；

步骤三、建立IFCB的变化率模型，并根据步骤一中IFCB时间序列计算IFCB变化率，拟合得到系数h、λ_i、θ_i的值；

步骤四、基于预报时间弧段的初始时刻的IFCB值和步骤三中的系数h、λ_i、θ_i，计算得到系数m的值；

步骤五、将系数h、λ_i、θ_i、系数m代入式(1)，即得更新后IFCB预报模型，对下一时间弧段的每一时刻进行预报；

其中，步骤三计算系数h、λ_i、θ_i的方法具体为：将公式(1)展开得式(2)：

其中，

采用最小二乘估计对式(2)求导，即得到IFCB变化率模型，表达如式(4)所示：

其中，为t时刻的IFCB变化率，基于步骤一中估计的IFCB序列，采用一阶中心差分法估计t时刻的IFCB变化率：

其中，IFCB(t+1)和IFCB(t-1)分别为与t时刻两侧相邻历元的IFCB估计值，Δt为IFCB采样间隔；

基于IFCB时间序列，根据式(5)计算得到对应历元的IFCB变化率序列，基于IFCB变化率序列，根据式(4)，计算得到系数h、a_i、b_i的值，a_i、b_i为式(2)中周期项系数。

2.如权利要求1所述的面向实时应用的GPS频间钟差新预报方法，其特征在于，步骤一中采用历元间差分方法基于GPS三频载波相位观测值估计IFCB时间序列。

3.如权利要求1所述的面向实时应用的GPS频间钟差新预报方法，其特征在于，针对IFCB预报模型，谐波阶数n为六阶，且周期T_i(i＝1～6)分别为12h、8h、6h、4.8h、4h和3h。

4.如权利要求1所述的面向实时应用的GPS频间钟差新预报方法，其特征在于，步骤四中系数m的计算公式如式(6)所示：

其中，t₀为预报时间弧段的初始时刻。

5.如权利要求1所述的面向实时应用的GPS频间钟差新预报方法，其特征在于，时间弧段设置为24h。

6.电子设备，其特征在于，包括：

一个或多个处理器；

存储器；以及

一个或多个程序，其中所述一个或多个程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，所述程序包括用于执行权利要求1～5任一项所述方法的指令。

7.计算机可读存储介质，其特征在于，存储于具有一种电子设备结合使用的计算机程序，所述计算机程序可被处理器执行以实现权利要求1～5任一项所述的方法。