CN109100927B - 无溯源参考条件下时间统一系统时间自主同步方法 - Google Patents
无溯源参考条件下时间统一系统时间自主同步方法 Download PDFInfo
- Publication number
- CN109100927B CN109100927B CN201810825290.1A CN201810825290A CN109100927B CN 109100927 B CN109100927 B CN 109100927B CN 201810825290 A CN201810825290 A CN 201810825290A CN 109100927 B CN109100927 B CN 109100927B
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- time difference
- difference data
- adj
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000002159 abnormal effect Effects 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims description 8
- 230000001360 synchronised effect Effects 0.000 claims description 6
- 230000005856 abnormality Effects 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 2
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04F—TIME-INTERVAL MEASURING
- G04F5/00—Apparatus for producing preselected time intervals for use as timing standards
- G04F5/14—Apparatus for producing preselected time intervals for use as timing standards using atomic clocks
-
- G—PHYSICS
- G04—HOROLOGY
- G04G—ELECTRONIC TIME-PIECES
- G04G7/00—Synchronisation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Electric Clocks (AREA)
Abstract
本发明提供了一种无溯源参考条件下时间统一系统时间自主同步方法,测量时间统一系统输出时间与溯源参考之间的时间差数据;若系统正常,对时间差数据作最小二乘拟合,预测下一整点时刻时间统一系统输出时间与溯源参考时间之差,根据时间差数据的拟合函数计算得到频率控制量,对频率标准的频率进行调整;若系统异常,以最小二乘法求出时间差数据序列的拟合函数,预测下一调整周期内时间统一系统输出时间与溯源参考时间的时差,进而计算频率调整增量。本发明能够在时间溯源设备异常的条件下,依然与溯源参考时间保持高精度同步。
Description
技术领域
本发明涉及一种在无时间溯源参考条件下,时间统一系统时频信号与溯源参考时间自主同步的方法。
背景技术
时间统一系统是提供溯源参考时间信号和标准频率信号以实现整个系统的时间和频率的高精度同步,由各种电子设备组成的一套完整系统,简称“时统”。时统一般由频率标准、相位微调器、时间传递接收机、时间间隔计数器、数字钟、时频信号分配器和控制计算机等组成。高精度的时间统一系统一般配备原子钟作为频率标准。当原子钟自由振荡时,由于其固有频率老化特性或受其他环境因素的影响,原子钟输出的时间、频率信号具有明显的漂移特性。如果不对它们进行实时调整,时统输出的时间将逐渐偏离溯源参考时间。因此,时间统一系统为了提供溯源参考时间,必须对其输出的时频信号进行控制,使时统时间与国家溯源参考时间保持同步。在计量学的量值传递中,这种做法称为溯源。
现在的时统大多以授时型接收机、GNSS卫星共视、卫星双向时间传递或光纤时间传递等技术手段实现高精度的时间溯源。
为了实现时间溯源,时间统一系统首先用以上技术手段之一测量出时统时间与溯源参考时间的时差,根据时差结果,对原子钟输出的频率和相位进行精密的调整与控制,使时统时间与溯源参考时间同步。实现时间同步的调整与控制技术称为钟驾驭技术。实现对原子钟输出的频率和相位进行调整与控制的设备称为相位微调器。在高精度的时间应用场合,如守时与授时系统、卫星导航、航天和导弹试验、高速通信系统等,一般不采用“简单粗暴”的直接调整信号相位的方法,而是以不破坏原子钟稳定度为原则,结合原子钟的稳定度指标,使用最小化控制调整策略,保持输出的实时物理信号的准确度和稳定度,“缓步微调”信号频率,达到时间同步的目的。
在时间溯源设备工作稳定时,可以实时得到时统输出时间与溯源参考时间的时差,通过对原子钟的频率和相位进行精密的调整与控制,时间统一系统可以实现较高的时间同步精度。但当溯源参考时间异常,无法提供时差数据时,由于时统内频率标准固有的频率漂移特性,如果不对其进行调整与控制,时间统一系统输出的时间将逐渐与溯源参考时间偏离,造成系统失效。
发明内容
为了克服现有技术的不足,本发明提供一种无溯源参考条件下时间统一系统时间自主同步方法,能够使时间统一系统在时间溯源设备异常,无法确切获取时间统一系统时间与溯源参考时间时差的条件下,依然可以与溯源参考时间保持高精度同步,提高了无溯源参考时,时间统一系统的输出时间与溯源参考时间的同步精度。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1,将搭建时间统一系统所需要使用的频率标准送到保持有UTC(k)的时间实验室或国家一级时间计量单位,由实验室或计量单位准确计量频率标准的频率稳定度σy(τ)和频偏fd;
步骤2,由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差数据;判断当前系统运行状态,若系统正常运行,则转入步骤3;若系统异常,则转入步骤7;
步骤3,对步骤2得到的时间差数据进行预处理,剔除奇异值;
步骤4,对时间差数据中已经剔除的数据位进行内插,保证时间差数据的连续性;
步骤5,对步骤4得到的时间差数据作最小二乘拟合,预测下一小时整点时刻时间统一系统输出时间与溯源参考时间之差;如果预测值超过设定的阈值,则进入步骤6;否则,返回步骤2;
步骤6,根据时间差数据的拟合函数y(t)=a0+a1t+a2t2计算得到频率控制量fadj;系统初始时刻频率控制量f’adj=—fd,则当前时刻的频率控制量为f’adj+fadj;通过频率控制量对频率标准的频率进行调整;然后返回步骤2;
步骤7,以最小二乘法求出时间差数据序列的拟合函数,拟合函数y(t)=a0+a1t+a2t2,设定频率调整周期为一小时,以此函数预测下一调整周期内,时间统一系统输出时间与溯源参考时间的时差,每一个调整周期起始时刻的预测时差值记为y0,
y0=y(86400+N×3600)=a0+a1×(86400+N×3600)+a2×(86400+N×3600)2;
每一调整周期终止时刻的预测时差值记为y1,
y1=y(86400+(N+1)×3600)=a0+a1×(86400+(N+1)×3600)+a2×(86400+(N+1)×3600)2;
其中N=0,1,2,3……是时间统一系统进入自主同步工作模式后经历的小时数;
步骤8,若|y0|>|y1|,则不调整频率标准的频率,维持当前时刻系统状态;若|y0|≤|y1|,转入步骤9;
步骤9,以频率标准的万秒稳定度σy(10k)作为频率调整增量,在系统异常后的每小时起始时刻调整频率标准输出信号的频率,并修正方程中的a1项,修正后的a1记为a1′,预测时差的拟合函数修正为y(t)=a0+a1′t+a2t2,当y1>0时,fadj的增量为-σy(10k),a1′=a1-σy(10k);当y1≤0时,fadj的增量为+σy(10k),a1′=a1+σy(10k);重复本步骤,直至|y0|>|y1|时,转入步骤10;
步骤10,若y0×y1>0,系统维持当前状态,若y0×y1≤0,则将频率控制量设置为fadj=-fd,直至系统恢复正常,转入步骤2。
所述的步骤1中,没有送检条件时,参照频率标准出厂说明书给出的频率稳定度及频偏。
所述的步骤2中,保证时间统一系统正常运行超过24小时,然后由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差数据。
所述的步骤2中,时间统一系统连续输出数据,则视为系统正常运行,数据缺失2小时以上,则视为系统异常。
所述的步骤3中,以3σ准则判别步骤2得到的时间差,剔除奇异值,并将奇异值以及发生的时刻记录到系统日志文件中。
本发明的有益效果是:在无法获取与时间溯源参考时差的情况下,使时间统一系统输出的时频信号依然可以与时间溯源参考保持高精度的同步。同时,本发明提出的小量步进调整的控制策略,保证了系统可以平稳、连续的过渡到自主同步状态,避免了系统输出的时频信号因溯源链路中断,而引起的信号频率或相位的突变,致使系统失效的后果。
附图说明
图1是时间统一系统原理图。
具体实施方式
下面结合附图和实施例对本发明进一步说明,本发明包括但不仅限于下述实施例。
本发明的技术特征在于通过准确计量时间统一系统中频率标准的频率稳定度与漂移特性,在保持频率标准自身稳定度特性的同时,保证时间统一系统在溯源链路异常的条件下,维持与溯源参考时间的高精度同步。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:将搭建时间统一系统所需要使用的频率标准(原子钟)送到保持有UTC(k)的时间实验室或国家一级时间计量单位,由实验室或计量单位准确计量频率标准的频率稳定度σy(τ)、频偏fd。没有送检条件时,参照频率标准出厂说明书给出的频率稳定度及频偏fd。
步骤2:当时间统一系统软硬件搭建完成后,保证系统正常运行超过24小时。由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差。数据连续输出,视为系统正常运行,转入步骤3,数据缺失2小时以上,视为系统异常,转入步骤7,直到数据正常后转入步骤3。
步骤3:对步骤2中得到的时间差数据进行预处理。以3σ准则进行判别,剔除奇异值,并将奇异值以及发生的时刻记录到系统日志文件中。
步骤4:对时间差数据中已经剔除的数据进行内插,保证时间差数据的连续性,并将此时间差数据存储在控制计算机内。
步骤5:根据步骤4中得到的这组时间差数据作最小二乘拟合,求出拟合二次方程,并以此方程预测下一小时整点时刻时间统一系统输出时间与溯源参考时间之差。如果预测值超过系统设定的阈值(由系统频率标准性能决定),需要对频率标准的频率进行调整,转入步骤6。反之,不调整频率标准的频率,维持当前时刻系统状态。转入步骤2。
步骤6:根据步骤5中时差数据的拟合函数y(t)=a0+a1t+a2t2可计算得到频率控制量,记为fadj。系统初始时刻频率控制量为频率标准的反向,即f’adj=—fd,则当前时刻的频率控制量为f’adj+fadj。通过控制计算机软件将频率控制量置入调整设备中。转入步骤2,继续测量时间差数据。
步骤7:当无法获取时间统一系统输出时间与溯源参考时间的时间差,时间统一系统进入自主同步的工作模式。根据系统异常时刻之前24小时的时差数据,以最小二乘法求出以上时差数据序列的拟合函数,拟合函数以二次多项式的形式表达:y(t)=a0+a1t+a2t2,并且设定频率调整周期为一小时,以此函数预测下一调整周期内,时间统一系统输出时间与溯源参考时间的时差,每一个调整周期起始时刻的预测时差值记为y0,y0=y(86400+N×3600)=a0+a1×(86400+N×3600)+a2×(86400+N×3600)2;每一调整周期终止时刻的预测时差值记为y1,y1=y(86400+(N+1)×3600)=a0+a1×(86400+(N+1)×3600)+a2×(86400+(N+1)×3600)2,其中N=0,1,2,3……是时间统一系统进入自主同步工作模式后,经历的小时数。
步骤8:若|y0|>|y1|,则不调整频率标准的频率,维持当前时刻系统状态,若|y0|≤|y1|,转入步骤9。
步骤9:以频率标准的万秒稳定度σy(10k)作为频率调整增量,在系统异常后的每小时起始时刻调整频率标准输出信号的频率,并修正方程中的a1项,修正后的a1记为a1′,预测时差的拟合函数修正为y(t)=a0+a1′t+a2t2,当y1>0时,fadj的增量为-σy(10k),a1′=a1-σy(10k);当y1≤0时,fadj的增量为+σy(10k),a1′=a1+σy(10k)。重复步骤9的计算过程,直至|y0|>|y1|时,转入步骤10。
步骤10:若y0×y1>0,系统维持当前状态,若y0×y1≤0,则将频率控制量设置为fadj=-fd,直至系统恢复正常。转入步骤2。
Claims (5)
1.一种无溯源参考条件下时间统一系统时间自主同步方法,其特征在于包括下述步骤:
步骤1,将搭建时间统一系统所需要使用的频率标准送到保持有UTC(k)的时间实验室或国家一级时间计量单位,由实验室或计量单位准确计量频率标准的频率稳定度σy(τ)和频偏fd;
步骤2,由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差数据;判断当前系统运行状态,若系统正常运行,则转入步骤3;若系统异常,则转入步骤7;
步骤3,对步骤2得到的时间差数据进行预处理,剔除奇异值;
步骤4,对时间差数据中已经剔除的数据位进行内插,保证时间差数据的连续性;
步骤5,对步骤4得到的时间差数据作最小二乘拟合,预测下一小时整点时刻时间统一系统输出时间与溯源参考时间之差;如果预测值超过设定的阈值,则进入步骤6;否则,返回步骤2;
步骤6,根据时间差数据的拟合函数y(t)=a0+a1t+a2t2计算得到频率控制量fadj,即根据步骤5中得到整点时刻时间统一系统输出时间与溯源参考时间之差的预报值y(t),计算得到频率控制量fadj,a0为常数项,a1为一次项系数,a2为二次项系数;系统初始时刻频率控制量f’adj=—fd,则当前时刻的频率控制量为f’adj+fadj;通过频率控制量对频率标准的频率进行调整;然后返回步骤2;
步骤7,以最小二乘法求出时间差数据序列的拟合函数,拟合函数y(t)=a0+a1t+a2t2,设定频率调整周期为一小时,以此函数预测下一调整周期内,时间统一系统输出时间与溯源参考时间的时差,每一个调整周期起始时刻的预测时差值记为y0,
y0=y(86400+N×3600)=a0+a1×(86400+N×3600)+a2×(86400+N×3600)2;
每一调整周期终止时刻的预测时差值记为y1,
y1=y(86400+(N+1)×3600)=a0+a1×(86400+(N+1)×3600)+a2×(86400+(N+1)×3600)2;
其中N=0,1,2,3……是时间统一系统进入自主同步工作模式后经历的小时数;
步骤8,若|y0|>|y1|,则不调整频率标准的频率,维持当前时刻系统状态;若|y0|≤|y1|,转入步骤9;
步骤9,以频率标准的万秒稳定度σy(10k)作为频率调整增量,在系统异常后的每小时起始时刻调整频率标准输出信号的频率,并修正方程中的a1项,修正后的a1记为a1′,预测时差的拟合函数修正为y(t)=a0+a1′t+a2t2,当y1>0时,fadj的增量为-σy(10k),a1′=a1-σy(10k);当y1≤0时,fadj的增量为+σy(10k),a1′=a1+σy(10k);重复本步骤,直至|y0|>|y1|时,转入步骤10;
步骤10,若y0×y1>0,系统维持当前状态,若y0×y1≤0,则将频率控制量设置为fadj=-fd,直至系统恢复正常,转入步骤2。
2.根据权利要求1所述的无溯源参考条件下时间统一系统时间自主同步方法,其特征在于:所述的步骤1中,没有送检条件时,参照频率标准出厂说明书给出的频率稳定度及频偏。
3.根据权利要求1所述的无溯源参考条件下时间统一系统时间自主同步方法,其特征在于:所述的步骤2中,保证时间统一系统正常运行超过24小时,然后由时间溯源设备测量得到时间统一系统输出时间与溯源参考之间的时间差数据。
4.根据权利要求1所述的无溯源参考条件下时间统一系统时间自主同步方法,其特征在于:所述的步骤2中,时间统一系统连续输出数据,则视为系统正常运行,数据缺失2小时以上,则视为系统异常。
5.根据权利要求1所述的无溯源参考条件下时间统一系统时间自主同步方法,其特征在于:所述的步骤3中,以3倍标准差准则判别步骤2得到的时间差,剔除奇异值,并将奇异值以及发生的时刻记录到系统日志文件中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810825290.1A CN109100927B (zh) | 2018-07-25 | 2018-07-25 | 无溯源参考条件下时间统一系统时间自主同步方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810825290.1A CN109100927B (zh) | 2018-07-25 | 2018-07-25 | 无溯源参考条件下时间统一系统时间自主同步方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100927A CN109100927A (zh) | 2018-12-28 |
CN109100927B true CN109100927B (zh) | 2021-02-26 |
Family
ID=64847419
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810825290.1A Expired - Fee Related CN109100927B (zh) | 2018-07-25 | 2018-07-25 | 无溯源参考条件下时间统一系统时间自主同步方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100927B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113078901B (zh) * | 2021-03-30 | 2024-02-13 | 中国科学院国家授时中心 | 基于脉冲星驾驭的原子钟准确度提高方法 |
CN113377003B (zh) * | 2021-06-07 | 2022-08-30 | 北京酷鲨科技有限公司 | 一种铯和氢原子频标的相位调整方法及装置 |
CN114114892B (zh) * | 2021-11-29 | 2023-07-25 | 成都天奥电子股份有限公司 | 一种产生可信时间的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7289820B2 (en) * | 2002-11-28 | 2007-10-30 | Nec Infrontia Corporation | Time-matching system and time-matching method |
CN101231337A (zh) * | 2008-02-15 | 2008-07-30 | 哈尔滨工程大学 | 高精度时间同步装置 |
CN104238352A (zh) * | 2014-09-24 | 2014-12-24 | 北京市计量检测科学研究院 | 一种国家基准驾驭的地方时间标准生成系统及方法 |
CN105676627A (zh) * | 2015-12-25 | 2016-06-15 | 中国科学院国家授时中心 | 守时系统主备主钟无缝切换系统与方法 |
CN105867108A (zh) * | 2016-06-20 | 2016-08-17 | 中国科学院国家授时中心 | 一种国家标准时间远程复现方法 |
-
2018
- 2018-07-25 CN CN201810825290.1A patent/CN109100927B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7289820B2 (en) * | 2002-11-28 | 2007-10-30 | Nec Infrontia Corporation | Time-matching system and time-matching method |
CN101231337A (zh) * | 2008-02-15 | 2008-07-30 | 哈尔滨工程大学 | 高精度时间同步装置 |
CN104238352A (zh) * | 2014-09-24 | 2014-12-24 | 北京市计量检测科学研究院 | 一种国家基准驾驭的地方时间标准生成系统及方法 |
CN105676627A (zh) * | 2015-12-25 | 2016-06-15 | 中国科学院国家授时中心 | 守时系统主备主钟无缝切换系统与方法 |
CN105867108A (zh) * | 2016-06-20 | 2016-08-17 | 中国科学院国家授时中心 | 一种国家标准时间远程复现方法 |
Non-Patent Citations (1)
Title |
---|
基于UTCr的主钟频率驾驭方法研究;闫敏;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20151215;C042-605页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109100927A (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100927B (zh) | 无溯源参考条件下时间统一系统时间自主同步方法 | |
CN101986555B (zh) | 用于定时模块故障保持的嵌入式自测试的系统和方法 | |
CN109299496B (zh) | 一种高精度同步时钟产生方法 | |
CN1918940B (zh) | 在网络分段上保持公共时间感测的系统和方法 | |
CN109581856B (zh) | 一种基于高性能晶振频率校准的对时守时方法 | |
US10505652B2 (en) | Methods and systems for estimating offset skew and drift | |
CN103141039A (zh) | 用于在通讯网络中的时间同步的方法 | |
CN113992296B (zh) | 时钟驯服方法、时间码监测装置及时间同步系统 | |
CN107086901B (zh) | 一种bdt建立方法及utc(ntsc)建立方法 | |
CN109302255B (zh) | 时间同步控制方法、装置、系统及计算机可读存储介质 | |
CN114205045B (zh) | 一种tte网络时钟校准方法及系统 | |
CN116125512B (zh) | 顾及时钟频率时变特性的ppp自适应钟差模型估算方法 | |
CN112636857A (zh) | 一种基于残差修正的高精度时频信号动态驾驭方法 | |
CN112433235B (zh) | 一种用于确定时间基准的方法、系统和介质 | |
CN104702214A (zh) | 一种晶体振荡器频率补偿的方法 | |
CN109557371A (zh) | 一种用于配电网相量测量的同步授时守时方法 | |
CN110752890A (zh) | 一种时间同步方法、时间同步系统及车辆 | |
CN110138491B (zh) | 一种基于广义加权最小二乘法的gps时钟补偿方法 | |
Kaczmarek et al. | Integrated system for monitoring and control of the national time and frequency standard | |
CN112969229B (zh) | 一种时钟校正方法、装置以及网络设备 | |
CN111865465B (zh) | 一种物联网设备精准校时方法 | |
CN110611544A (zh) | 一种利用rs232串口精确对时的方法及装置 | |
CN115639743B (zh) | 一种基于整网时间比对的天基时间基准建立方法及系统 | |
CN114371489B (zh) | 基于星钟参数线性预测的北斗三号星间链路快速恢复方法 | |
CN110618604B (zh) | 一种利用ntp辅助作源提高守时精度的方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210226 |
|
CF01 | Termination of patent right due to non-payment of annual fee |