[发明专利]一种铷钟频率修正值的预估方法

权利要求书

1.一种铷钟频率修正值的预估方法，其特征在于：所述方法包括：

S1，建立铷钟相对UTC(NIM)的频差的预估模型：

y＝D*(t-t₀)+y₀(1)

其中，y表示在未来时刻t，铷钟相对UTC(NIM)的预估频差；D是铷钟的频率漂移；y₀表示在时刻t₀铷钟相对UTC(NIM)的频差；

S2，分别计算铷钟的频率漂移D和在时刻t₀铷钟相对UTC(NIM)的频差y₀；

S3，用预估方程预测未来m天数据；

S4，用补偿量对S3得到的预测值进行补偿；

其中，所述S2是这样实现的：

使用公式(2)获得铷钟的频率漂移D：

D_n＝a×(s_n-D_n-1)+D_n-1(2)

其中，D_n表示在第n次计算出的铷钟频率漂移；D_n-1表示在第n-1次计算出的铷钟频率漂移，用最近60天铷钟相对UTC(NIM)频差实测值做直线拟合的斜率来赋初值；s_n表示最近60天铷钟相对UTC(NIM)频差数据作一次线性拟合计算出的斜率；系数a取值范围是0.5～1.5；

使用公式(3)获得在时刻t₀铷钟相对UTC(NIM)的频差y₀：

$y_{0n}=b\×\frac{1}{m}\underset{i=1}{\overset{m}{\mathrm{\Σ}}}(f_i-y_i)+y_{0n-1}---\left(3\right)$

其中，y_0n表示在第n次计算出的铷钟相对UTC(NIM)的频差；y_0n-1表示在第n-1次计算出的铷钟相对UTC(NIM)的频差，用铷钟实测值减去UTC(NIM)的频差赋初值；f_i表示第i天实际测量的铷钟相对UTC(NIM)的频差；y_i表示第i天预估的铷钟相对UTC(NIM)的频差；m表示使用第n次计算前多少天来计算本次y_0n；系数b取值范围是0.1～2.0。

2.根据权利要求1所述的铷钟频率修正值的预估方法，其特征在于：所述a取1。

3.根据权利要求1所述的铷钟频率修正值的预估方法，其特征在于：所述b取0.8。

4.根据权利要求1所述的铷钟频率修正值的预估方法，其特征在于：所述S3是这样实现的：将S2得到的铷钟的频率漂移D和在时刻t₀铷钟相对UTC(NIM)的频差y0代入到公式(1)中得到直线方程，将直线外推即得到未来m天数据。

5.根据权利要求4所述的铷钟频率修正值的预估方法，其特征在于：所述S4是这样实现的：

补偿量由公式(4)计算：

$Y_n=\underset{i=1}{\overset{m}{\mathrm{\Σ}}}(f_i-y_i)---\left(4\right)$

式中，Y_n表示在第n次计算的补偿量；m表示第n-1次计算到第n次计算间的天数；

将补偿量取均值加入到预测的未来m天数据中。