CN109100927A - 无溯源参考条件下时间统一***时间自主同步方法 - Google Patents
无溯源参考条件下时间统一***时间自主同步方法 Download PDFInfo
- Publication number
- CN109100927A CN109100927A CN201810825290.1A CN201810825290A CN109100927A CN 109100927 A CN109100927 A CN 109100927A CN 201810825290 A CN201810825290 A CN 201810825290A CN 109100927 A CN109100927 A CN 109100927A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- source
- time difference
- timing system
- 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.)
- Granted
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 18
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005259 measurement Methods 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 4
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 238000007689 inspection Methods 0.000 claims description 3
- 230000000630 rising effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000013632 homeostatic process Effects 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
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;***初始时刻频率控制量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 true CN109100927A (zh) | 2018-12-28 |
CN109100927B 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) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113078901A (zh) * | 2021-03-30 | 2021-07-06 | 中国科学院国家授时中心 | 基于脉冲星驾驭的原子钟准确度提高方法 |
CN113377003A (zh) * | 2021-06-07 | 2021-09-10 | 北京酷鲨科技有限公司 | 一种铯和氢原子频标的相位调整方法及装置 |
CN114114892A (zh) * | 2021-11-29 | 2022-03-01 | 成都天奥电子股份有限公司 | 一种产生可信时间的方法 |
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的主钟频率驾驭方法研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113078901A (zh) * | 2021-03-30 | 2021-07-06 | 中国科学院国家授时中心 | 基于脉冲星驾驭的原子钟准确度提高方法 |
CN113078901B (zh) * | 2021-03-30 | 2024-02-13 | 中国科学院国家授时中心 | 基于脉冲星驾驭的原子钟准确度提高方法 |
CN113377003A (zh) * | 2021-06-07 | 2021-09-10 | 北京酷鲨科技有限公司 | 一种铯和氢原子频标的相位调整方法及装置 |
CN114114892A (zh) * | 2021-11-29 | 2022-03-01 | 成都天奥电子股份有限公司 | 一种产生可信时间的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109100927B (zh) | 2021-02-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100927A (zh) | 无溯源参考条件下时间统一***时间自主同步方法 | |
US8797827B2 (en) | Seismic acquisition system | |
CN103913753B (zh) | 一种利用导航卫星的高精度授时***及方法 | |
CN109302255B (zh) | 时间同步控制方法、装置、***及计算机可读存储介质 | |
Kebkal et al. | Underwater acoustic modems with integrated atomic clocks for one-way travel-time underwater vehicle positioning | |
CN107247405A (zh) | 一种卫星导航授时守时在时钟同步***中的解决方法 | |
EP2811672A1 (en) | Communication apparatus, communication system, and communication method for time synchronization | |
CN104238352A (zh) | 一种国家基准驾驭的地方时间标准生成***及方法 | |
CN101022373A (zh) | 一种测试远程帧时延的***及方法及装置 | |
CN105049040A (zh) | 一种利用gnss校正相干布居囚禁原子钟输出频率的方法 | |
CN108803300B (zh) | 基于恒温晶振的时间同步装置守时方法和时间同步装置 | |
CN103529457A (zh) | 一种用于卫星导航信号模拟器自身时延校准的***及方法 | |
CN114205045B (zh) | 一种tte网络时钟校准方法及*** | |
CN103117742A (zh) | Gps/北斗双模卫星时钟晶体振荡器驯服*** | |
CN112636857A (zh) | 一种基于残差修正的高精度时频信号动态驾驭方法 | |
US8224606B2 (en) | Measuring clock jitter | |
CN104716904A (zh) | 一种晶体振荡器频率补偿的方法 | |
CN106383438B (zh) | 一种基于滑动窗口时间扩展高精度驯钟方法 | |
CN107145061A (zh) | 卫星授时信号的智能抗干扰技术在时间同步中的防御方法 | |
WO2022001368A1 (zh) | 雷达及雷达*** | |
CN110611544A (zh) | 一种利用rs232串口精确对时的方法及装置 | |
CN110471087A (zh) | 一种空间飞行器的时间漂移计算方法和*** | |
CN110618604B (zh) | 一种利用ntp辅助作源提高守时精度的方法及装置 | |
CN105589328B (zh) | 时间同步测试方法、测试精度确定方法及装置 | |
CN110687773B (zh) | 时间统一***授时精度的测量方法、装置和*** |
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 |