[发明专利]无溯源参考条件下时间统一系统时间自主同步方法

权利要求书

1.一种无溯源参考条件下时间统一系统时间自主同步方法，其特征在于包括下述步骤：

步骤1，将搭建时间统一系统所需要使用的频率标准送到保持有UTC(k)的时间实验室或国家一级时间计量单位，由实验室或计量单位准确计量频率标准的频率稳定度σ_y(τ)和频偏f_d；

步骤2，由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差数据；判断当前系统运行状态，若系统正常运行，则转入步骤3；若系统异常，则转入步骤7；

步骤3，对步骤2得到的时间差数据进行预处理，剔除奇异值；

步骤4，对时间差数据中已经剔除的数据位进行内插，保证时间差数据的连续性；

步骤5，对步骤4得到的时间差数据作最小二乘拟合，预测下一小时整点时刻时间统一系统输出时间与溯源参考时间之差；如果预测值超过设定的阈值，则进入步骤6；否则，返回步骤2；

步骤6，根据时间差数据的拟合函数y(t)＝a₀+a₁t+a₂t²计算得到频率控制量f_adj，即根据步骤5中得到整点时刻时间统一系统输出时间与溯源参考时间之差的预报值y(t)，计算得到频率控制量f_adj，a₀为常数项，a₁为一次项系数，a₂为二次项系数；系统初始时刻频率控制量f’_adj＝—f_d，则当前时刻的频率控制量为f’_adj+f_adj；通过频率控制量对频率标准的频率进行调整；然后返回步骤2；

步骤7，以最小二乘法求出时间差数据序列的拟合函数，拟合函数y(t)＝a₀+a₁t+a₂t²，设定频率调整周期为一小时，以此函数预测下一调整周期内，时间统一系统输出时间与溯源参考时间的时差，每一个调整周期起始时刻的预测时差值记为y₀，

y₀＝y(86400+N×3600)＝a₀+a₁×(86400+N×3600)+a₂×(86400+N×3600)²；

每一调整周期终止时刻的预测时差值记为y₁，

y₁＝y(86400+(N+1)×3600)＝a₀+a₁×(86400+(N+1)×3600)+a₂×(86400+(N+1)×3600)²；

其中N＝0,1,2,3……是时间统一系统进入自主同步工作模式后经历的小时数；

步骤8，若|y₀||y₁|，则不调整频率标准的频率，维持当前时刻系统状态；若|y₀|≤|y₁|，转入步骤9；

步骤9，以频率标准的万秒稳定度σ_y(10k)作为频率调整增量，在系统异常后的每小时起始时刻调整频率标准输出信号的频率，并修正方程中的a₁项，修正后的a₁记为a₁′，预测时差的拟合函数修正为y(t)＝a₀+a₁′t+a₂t²，当y₁0时，f_adj的增量为-σ_y(10k)，a₁′＝a₁-σ_y(10k)；当y₁≤0时，f_adj的增量为+σ_y(10k)，a₁′＝a₁+σ_y(10k)；重复本步骤，直至|y₀||y₁|时，转入步骤10；

步骤10，若y₀×y₁0，系统维持当前状态，若y₀×y₁≤0，则将频率控制量设置为f_adj＝-f_d，直至系统恢复正常，转入步骤2。