CN104199274A - 一种铷钟频率修正值的预估方法 - Google Patents
一种铷钟频率修正值的预估方法 Download PDFInfo
- Publication number
- CN104199274A CN104199274A CN201410494230.8A CN201410494230A CN104199274A CN 104199274 A CN104199274 A CN 104199274A CN 201410494230 A CN201410494230 A CN 201410494230A CN 104199274 A CN104199274 A CN 104199274A
- Authority
- CN
- China
- Prior art keywords
- rubidium clock
- frequency
- nim
- utc
- relative
- 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
Landscapes
- Electric Clocks (AREA)
Abstract
本发明提供了一种铷钟频率修正值的预估方法,属于时间频率校准领域。本方法包括:S1,建立铷钟相对UTC(NIM)的频差的预估模型:y=D*(t-t0)+y0(1)其中,式中,y表示在未来时刻t,铷钟相对UTC(NIM)的预估频差;D是铷钟的频率漂移;y0表示在时刻t0铷钟相对UTC(NIM)的频差;S2,分别计算铷钟的频率漂移D和在时刻t0铷钟相对UTC(NIM)的频差y0;S3,用预估方程预测未来m天数据;S4,用补偿量对S3得到的预测值进行补偿。
Description
技术领域
本发明属于时间频率校准领域,具体涉及一种铷钟频率修正值的预估方法。
背景技术
建立准确、稳定、可靠的地方原子频标是建立北斗卫星导航检测平台的基础,是实现对平台关键参数进行量值传递的必要条件。
以北京计量院时间频率计量标准铷原子钟作为参考钟,通过研制的时间比对接收系统,每隔16分钟,测量GPS或北斗时间与这台铷参考钟之间的时差,并按照CGGTTS格式文件保存测量数据。在中国计量院网站上下载UTC(NIM)与GPS或北斗之差(CGGTTS格式文件),相减得到UTC(NIM)与铷钟之间的修正值(每16分钟/个)。
通过这种方法使铷钟的频率受到UTC(NIM)的驾驭。铷钟的频率准确度一般为5×10-12,修正后铷钟的日平均准确度达到5×10-14。在UTC(NIM)的驾驭下,由于每16分钟与UTC(NIM)进行一次比对,一旦突然发生时刻和速率的突变都会造成时标的不连续和频率的偏差,因此需要利用UTC(NIM)的频率进行修正,在累积的长期测试数据基础上,建立铷钟相对UTC(NIM)的频差修正值预估计算模型,从而实现对铷钟频率的预估,提高铷钟的频率准确度。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种铷钟频率修正值的预估方法,提高铷钟频率准确度。在累积的铷钟日常测试数据的基础上,分析铷钟的频率漂移规律,拟合频率漂移曲线,建立铷钟相对UTC(NIM)的频差修正值预估计算模型,同时利用相位补偿技术,建立相位补偿量方程,从而实现对本地铷钟频率的预估,提高铷钟频率准确度,实现与国家时间基准UTC(NIM)同步在±20ns以内。
本发明是通过以下技术方案实现的:
一种铷钟频率修正值的预估方法,包括:
S1,建立铷钟相对UTC(NIM)的频差的预估模型:
y=D*(t-t0)+y0 (1)
其中,y表示在未来时刻t,铷钟相对UTC(NIM)的预估频差;D是铷钟的频率漂移;y0表示在时刻t0铷钟相对UTC(NIM)的频差;
S2,分别计算铷钟的频率漂移D和在时刻t0铷钟相对UTC(NIM)的频差y0;
S3,用所述预估模型预测未来m天数据;
S4,用补偿量对S3得到的预测值进行补偿。
所述S2是这样实现的:
使用公式(2)获得铷钟的频率漂移D:
Dn=a×(sn-Dn-1)+Dn-1 (2)
其中,Dn表示在第n次计算出的铷钟频率漂移;Dn-1表示在第n-1次计算出的铷钟频率漂移;sn表示最近60天铷钟相对UTC(NIM)频差数据作一次线性拟合计算出的斜率;系数a取值范围是0.5~1.5;
使用公式(3)获得在时刻t0铷钟相对UTC(NIM)的频差y0:
其中,y0n表示在第n次计算出的铷钟相对UTC(NIM)的频差;y0n-1表示在第n-1次计算出的铷钟相对UTC(NIM)的频差;f表示实际测量的铷钟相对UTC(NIM)的频差;y表示预估的铷钟相对UTC(NIM)的频差(即前一次预估的频差);m表示使用第n次计算前多少天来计算本次y0n;系数b取值范围是0.1~2.0。
所述a取1。
所述b取0.8。
所述S3是这样实现的:
将S2得到的铷钟的频率漂移D和在时刻t0铷钟相对UTC(NIM)的频差y0代入到公式(1)中得到直线方程,将直线外推即得到未来m天数据。
所述S4是这样实现的:
补偿量由公式(4)计算:
式中,Yn表示在第n次计算的补偿量;m表示第n-1次计算到第n次计算间的天数;
将补偿量取均值加入到预测的未来m天数据中。
与现有技术相比,本发明的有益效果是:通过建立铷钟相对UTC(NIM)的频差修正值预估计算模型,和相位补偿量方程,使铷钟的频率与UTC(NIM)保持一致,从而实现对本地铷钟频率的预估,提高铷钟频率准确度,实现与国家时间基准UTC(NIM)同步在±20ns以内。
附图说明
图1铷钟频差数据作一次线性拟合,拟合长度为60天和75天,采样周期为10天
图2铷钟发生频率跳变时的规律
图3实测修正值预估曲线。
图4本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
(1)建立频率漂移拟合曲线方程
Dn=a×(sn-Dn-1)+Dn-1
式中,Dn表示在第n次计算出的铷钟频率漂移,Dn-1表示在第n-1次计算出的铷钟频率漂移,sn表示最近60天铷钟相对UTC(NIM)频差数据作一次线性拟合计算出的斜率。系数a取值范围是0.5~1.5,由操作人员根据经验判断给出,一般情况下a取1。
(2)建立频差修正值评估方程
式中,y0n表示在第n次计算出的铷钟相对UTC(NIM)的频差。y0n-1表示在第n-1次计算出的相对频差,可用铷钟实测值减去UTC(NIM)的频差赋初值。f表示实际测量的铷钟相对UTC(NIM)的频差。y表示预估的铷钟相对UTC(NIM)的频差。m表示使用第n次计算前多少天计算本次y0n,系数b取值范围是0.1~2.0,由操作人员根据跳变幅度确定取值大小,通常情况下b为0.8。
(3)建立相位补偿量方程
式中,Yn表示在第n次计算的补偿量。f表示实际测量的铷钟相对UTC(NIM)的频差。y表示预估的铷钟相对UTC(NIM)的频差。m表示第n-1次计算到第n次计算间的天数。一般会用5天时间完成补偿过程,用5天平均可以减少因特殊原因造成铷钟频差突变造成的影响。m取值可从1~60。
如图4所示,本发明的一个实施例如下:
通过下面的公式(1)、(2)和(3),建立起频差修正值预估计算模型。预估的铷钟相对UTC(NIM)的频差可表示为公式(1),
y=D*(t-t0)+y0 (1)
式中,y表示在未来时刻t,铷钟相对UTC(NIM)的预估频差。D是铷钟的频率漂移。y0表示在时刻t0铷钟相对UTC(NIM)的频差。
通过研究(见图1)发现,铷钟的频率漂移具有以下规律:
1、频率漂移随时间发生变化,变化曲线为类似正弦波的不规则曲线。
2、使用一次线性拟合评估频率漂移时,拟合长度为60天较合适。
根据上述规律,使用公式(2)评估铷钟相对UTC(NIM)的频率漂移
Dn=a×(sn-Dn-1)+Dn-1 (2)
式中,Dn表示在第n次计算出的铷钟频率漂移。Dn-1表示在第n-1次计算出的铷钟频率漂移。sn表示最近60天铷钟相对UTC(NIM)频差数据作一次线性拟合计算出的斜率。系数a取值范围是0.5~1.5,由操作人员根据经验判断给出,一般情况下a取1。
图2展示了y0的变化规律,图中实线是所有数据作一次线性拟合得出的方程直线。三段虚线是发生频率跳变前后三个时间段内各自拟合得出的方程直线。其规律为:
1、铷钟会发生频率跳变,发生跳变的周期可长可短。
2、铷钟发生频率跳变前后的频率漂移可以认为保持不变。
根据上述规律,使用公式(3)评估y0
式中,y0n表示在第n次计算出的铷钟相对UTC(NIM)的频差。y0n-1表示在第n-1次计算出的相对频差。f表示实际测量的铷钟相对UTC(NIM)的频差。y表示预估的铷钟相对UTC(NIM)的频差。m表示使用第n次计算前多少天计算本次y0n。系数b取值范围是0.1~2.0,通常情况下b为0.8。
图3为实测的修正值预估曲线,包括:外部数据源(UTC(NIM)数据),ax+b(预估频率),SDI铷钟实际运行频率,7天预测频率(拟合曲线)。界面显示最近60天的相关数据,横坐标是日期,纵坐标是相对频差值。从图中可以看出修正值预估可以大大提升铷钟的准确度。
预估的铷钟相对UTC(NIM)间的频率差与实际测量的频率差之间一定存在着偏差,把这些偏差补偿上,铷钟的频率才能与UTC(NIM)的频率保持一致。补偿量由公式(4)计算
式中,Yn表示在第n次计算的补偿量。f表示实际测量的铷钟相对UTC(NIM)的频差。y表示预估的铷钟相对UTC(NIM)的频差。m表示第n-1次计算到第n次计算间的天数。
为保持时标的频率稳定度,一般会用5天时间完成补偿过程。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (6)
1.一种铷钟频率修正值的预估方法,其特征在于:所述方法包括:
S1,建立铷钟相对UTC(NIM)的频差的预估模型:
y=D*(t-t0)+y0 (1)
其中,y表示在未来时刻t,铷钟相对UTC(NIM)的预估频差;D是铷钟的频率漂移;y0表示在时刻t0铷钟相对UTC(NIM)的频差;
S2,分别计算铷钟的频率漂移D和在时刻t0铷钟相对UTC(NIM)的频差y0;
S3,用预估方程预测未来m天数据;
S4,用补偿量对S3得到的预测值进行补偿。
2.根据权利要求1所述的铷钟频率修正值的预估方法,其特征在于:所述S2是这样实现的:
使用公式(2)获得铷钟的频率漂移D:
Dn=a×(sn-Dn-1)+Dn-1 (2)
其中,Dn表示在第n次计算出的铷钟频率漂移;Dn-1表示在第n-1次计算出的铷钟频率漂移;sn表示最近60天铷钟相对UTC(NIM)频差数据作一次线性拟合计算出的斜率;系数a取值范围是0.5~1.5;
使用公式(3)获得在时刻t0铷钟相对UTC(NIM)的频差y0:
其中,y0n表示在第n次计算出的铷钟相对UTC(NIM)的频差;y0n-1表示在第n-1次计算出的铷钟相对UTC(NIM)的频差;f表示实际测量的铷钟相对UTC(NIM)的频差;y表示预估的铷钟相对UTC(NIM)的频差;m表示使用第n次计算前多少天来计算本次y0n;系数b取值范围是0.1~2.0。
3.根据权利要求2所述的铷钟频率修正值的预估方法,其特征在于:所述a取1。
4.根据权利要求2所述的铷钟频率修正值的预估方法,其特征在于:所述b取0.8。
5.根据权利要求2所述的铷钟频率修正值的预估方法,其特征在于:所述S3是这样实现的:将S2得到的铷钟的频率漂移D和在时刻t0铷钟相对UTC(NIM)的频差y0代入到公式(1)中得到直线方程,将直线外推即得到未来m天数据。
6.根据权利要求5所述的铷钟频率修正值的预估方法,其特征在于:所述S4是这样实现的:
补偿量由公式(4)计算:
式中,Yn表示在第n次计算的补偿量;m表示第n-1次计算到第n次计算间的天数;
将补偿量取均值加入到预测的未来m天数据中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410494230.8A CN104199274B (zh) | 2014-09-24 | 2014-09-24 | 一种铷钟频率修正值的预估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410494230.8A CN104199274B (zh) | 2014-09-24 | 2014-09-24 | 一种铷钟频率修正值的预估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104199274A true CN104199274A (zh) | 2014-12-10 |
CN104199274B CN104199274B (zh) | 2017-03-22 |
Family
ID=52084577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410494230.8A Expired - Fee Related CN104199274B (zh) | 2014-09-24 | 2014-09-24 | 一种铷钟频率修正值的预估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199274B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898404A (zh) * | 2015-06-25 | 2015-09-09 | 中国人民解放军国防科学技术大学 | 等价于二状态变量Kalman滤波器的数字锁相环原子钟驾驭方法 |
CN106154814A (zh) * | 2016-07-20 | 2016-11-23 | 中国计量科学研究院 | 时间信号产生系统、校准控制装置及校准控制方法 |
CN106773610A (zh) * | 2016-11-28 | 2017-05-31 | 北京工业大学 | 一种铯原子钟和氢钟频差预估方法 |
CN108229747A (zh) * | 2018-01-12 | 2018-06-29 | 中国计量科学研究院 | 校准控制方法、装置及时间信号产生系统 |
CN113641087A (zh) * | 2021-07-27 | 2021-11-12 | 中国科学院国家授时中心 | 光钟间歇运行驾驭微波钟产生高精度时间的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1973248A (zh) * | 2004-03-30 | 2007-05-30 | 国家科研中心 | 利用相干群体俘获调制原子钟信号的方法以及相应的原子钟 |
FR2924826A1 (fr) * | 2007-12-11 | 2009-06-12 | Commissariat Energie Atomique | Horloge atomique a correction du champ magnetique ambiant |
CN202102264U (zh) * | 2011-06-24 | 2012-01-04 | 成都可为科技发展有限公司 | 铷钟驯服系统 |
US20130335154A1 (en) * | 2011-03-01 | 2013-12-19 | National Research Council Of Canada | Frequency Stabilization of An Atomic Clock Against Variations of the C-Field |
CN203366003U (zh) * | 2013-05-31 | 2013-12-25 | 江汉大学 | 一种漂移优化的原子钟 |
-
2014
- 2014-09-24 CN CN201410494230.8A patent/CN104199274B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1973248A (zh) * | 2004-03-30 | 2007-05-30 | 国家科研中心 | 利用相干群体俘获调制原子钟信号的方法以及相应的原子钟 |
JP2007530965A (ja) * | 2004-03-30 | 2007-11-01 | サントル ナシオナル ドゥ ラ ルシェルシェサイアンティフィク(セエヌエールエス) | コヒーレント・ポピュレーション・トラッピングにより原子時計信号を変調する方法および対応する原子時計 |
FR2924826A1 (fr) * | 2007-12-11 | 2009-06-12 | Commissariat Energie Atomique | Horloge atomique a correction du champ magnetique ambiant |
US20130335154A1 (en) * | 2011-03-01 | 2013-12-19 | National Research Council Of Canada | Frequency Stabilization of An Atomic Clock Against Variations of the C-Field |
CN202102264U (zh) * | 2011-06-24 | 2012-01-04 | 成都可为科技发展有限公司 | 铷钟驯服系统 |
CN203366003U (zh) * | 2013-05-31 | 2013-12-25 | 江汉大学 | 一种漂移优化的原子钟 |
Non-Patent Citations (1)
Title |
---|
王亚莉: "精密时频系统中漂移率补偿处理技术研究", 《中国优秀硕士论文全文库 工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104898404A (zh) * | 2015-06-25 | 2015-09-09 | 中国人民解放军国防科学技术大学 | 等价于二状态变量Kalman滤波器的数字锁相环原子钟驾驭方法 |
CN104898404B (zh) * | 2015-06-25 | 2017-03-29 | 中国人民解放军国防科学技术大学 | 等价于二状态变量Kalman滤波器的数字锁相环原子钟驾驭方法 |
CN106154814A (zh) * | 2016-07-20 | 2016-11-23 | 中国计量科学研究院 | 时间信号产生系统、校准控制装置及校准控制方法 |
CN106154814B (zh) * | 2016-07-20 | 2018-12-07 | 中国计量科学研究院 | 时间信号产生系统、校准控制装置及校准控制方法 |
CN106773610A (zh) * | 2016-11-28 | 2017-05-31 | 北京工业大学 | 一种铯原子钟和氢钟频差预估方法 |
CN108229747A (zh) * | 2018-01-12 | 2018-06-29 | 中国计量科学研究院 | 校准控制方法、装置及时间信号产生系统 |
CN108229747B (zh) * | 2018-01-12 | 2020-12-15 | 中国计量科学研究院 | 校准控制方法、装置及时间信号产生系统 |
CN113641087A (zh) * | 2021-07-27 | 2021-11-12 | 中国科学院国家授时中心 | 光钟间歇运行驾驭微波钟产生高精度时间的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104199274B (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104199274A (zh) | 一种铷钟频率修正值的预估方法 | |
CN109425844B (zh) | 一种数据采样的校准方法和系统 | |
CN109428598B (zh) | 数据采样的校准方法和系统 | |
Li et al. | Modeling of the GPS satellite clock error and its performance evaluation in precise point positioning | |
US20140136093A1 (en) | Systems and methods of reducing timing measurement error due to clock offset | |
CN103336891A (zh) | 一种用于配电网状态估计的伪量测生成方法 | |
CN109283864B (zh) | 一种数据采样的时间同步、校准方法和系统 | |
Wöppelmann et al. | Rescue of the historical sea level record of Marseille (France) from 1885 to 1988 and its extension back to 1849–1851 | |
CN104238352A (zh) | 一种国家基准驾驭的地方时间标准生成系统及方法 | |
CN107086901B (zh) | 一种bdt建立方法及utc(ntsc)建立方法 | |
CN104040378B (zh) | 气象预测装置以及气象预测方法 | |
CN103119453A (zh) | 基于二次型的数字频率估计 | |
US8224606B2 (en) | Measuring clock jitter | |
CN116125512A (zh) | 顾及时钟频率时变特性的ppp自适应钟差模型估算方法 | |
CN104537233B (zh) | 一种基于核密度估计的配电网伪量测生成方法 | |
CN101483432A (zh) | 用于雷电时差探测站中的高精度gps时钟及其同步方法 | |
CN107678048B (zh) | 卫星钟的超快钟差预报模型的预先选择方法 | |
CN103605904B (zh) | 基于误差估算的自补偿电力系统幅值算法 | |
CN201323567Y (zh) | 用于雷电时差探测站中的高精度gps时钟 | |
CN102937710B (zh) | 一种确定原子钟非理想特性引起用户测距误差的方法 | |
CN110933652B (zh) | Nbiot传感设备及其周期同步采集方法、装置和介质 | |
CN201707439U (zh) | 北斗一号双向授时型用户机 | |
CN204065237U (zh) | 授时用时基参数测量装置 | |
CN205506934U (zh) | 一种冲击负荷电能计量系统 | |
Yang et al. | Numerical Simulation Research of Advanced Multiple Aperture Seeing Profiler |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170322 Termination date: 20190924 |