[发明专利]一种利用Algos和Kalman组合产生原子时标的方法

摘要

本发明公开一种利用Algos和Kalman组合产生原子时标的方法，首先对原始钟差数据进行离群值检测；计算前六个月的TA‑K哈德玛方差以确定每台钟的权重，作为计算本月TA‑K权重；其次用改进的Algos算法得到的结果，作为Kalman算法的噪声矩阵和初值的输入，利用哈德玛方差确定每台钟的权重；最后预测下一时刻的钟差值。采用本发明的技术方案，既可以降低Kalman算法的发散性，又能使时间尺度稳定性好。

主权项

1.一种利用Algos和Kalman组合产生原子时标的方法，其特征在于，包括以下步骤：步骤1、对实验室测得的原子钟钟差数据进行离群值检测；步骤2、根据步骤1得到的原子钟钟差数据计算每台原子钟分配权重；步骤3、计算实验室单个原子钟与实验室原子钟组平均时间尺度之差xi(t)，即TA‑K数据，若每台原子钟的权重大于最大权重上限，则该台原子钟的权重为最大权重上限；设h'i(t)为t时刻原子钟Hi时间改正量，则将权重、原子钟数据、h'i(t)通过如下公式计算：其中，xij(t)为原子钟Hj和原子钟Hi在t时刻的钟差值；xi(t)为实验室单个原子钟与实验室原子钟组平均时间尺度之差，xi(t)与xj(t)不相同，N为1…i…j…N，N为原子钟组中原子钟的数量；步骤4、根据TA‑K数据xi(t)和每台钟的权重计算钟组的原子时标；将x_i(t)数据作为Kalman滤波器模型的输入，计算t时刻原子钟频差值的预测值为：其中，A(t)为t时刻原子钟频差值，为A(t)的临时原子钟频差值预测值，K为Kalman滤波器的增益，H为测量矩阵，计算t时刻钟组的原子时标TA(t)：其中，ω_i为每台原子钟的权重，A_i(t)为每台原子钟的当前原子钟频差值，为每台原子钟的预测原子钟频差值。