CN113985719A - 基于滑动窗口的脉冲星时驯服铯原子钟的方法 - Google Patents
基于滑动窗口的脉冲星时驯服铯原子钟的方法 Download PDFInfo
- Publication number
- CN113985719A CN113985719A CN202111242954.XA CN202111242954A CN113985719A CN 113985719 A CN113985719 A CN 113985719A CN 202111242954 A CN202111242954 A CN 202111242954A CN 113985719 A CN113985719 A CN 113985719A
- Authority
- CN
- China
- Prior art keywords
- atomic clock
- pulsar
- time
- cesium atomic
- window
- 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
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Unknown Time Intervals (AREA)
- Electric Clocks (AREA)
Abstract
一种基于滑动窗口的脉冲星时驯服铯原子钟方法,包括:脉冲星计时观测;脉冲星计时数据处理;利用活动窗口计算铯原子钟频率偏差;铯原子钟频率校准。本发明解决了脉冲星脉冲到达时间测量误差大,对原子钟频率漂移测量精度影响,在窗口滑动步长不变的情况下,增加滑动窗口大小,有效提高了频率偏差测量精度。本发明将为基于脉冲星时实现原子钟的高精度频率溯源提供支撑,具有预报时间短、测量误差小、识别精度高、速度快的优点。
Description
技术领域
本发明属于应用原子钟技术领域,具体涉及到一种基于滑动窗口的脉冲星时驯服铯原子钟的方法。
背景技术
传统原子钟频率校准是将本地原子钟与基准原子钟信号进行比对测量,根据比对测量钟差数据,计算本地原子钟相对参考基准钟的频率偏差,将当前的频率偏差值作为下一段时间内频率偏差估计值,通过频率驾驭控制***对原子钟频率进行校准,实现本地钟与基准钟频率一致。
脉冲星是致密天体,具有强磁场、强电场特性,辐射稳定的周期性脉冲信号,被誉为自然界最稳定的“天然钟”。利用天文观测技术可实现脉冲星自转参数的精确测量,目前脉冲星自转频率测量不确定度已达到10-15量级,可作为基准频率源。当没有原子时基准频率参考信号时,利用地面脉冲星观测装置以原子钟为时间参考观测脉冲星信号,通过计时数据处理技术可获得原子钟与脉冲星时的钟差序列,若参考原子钟频率存在误差,则测量脉冲峰值相位相对于本征预报脉冲峰值相位随时间具有误差累积效应,如图1所示。利用钟差序列拟合得到原子钟频率漂移值,进而可实现对原子钟频率校准。被驯服原子钟最重要的特点是短期稳定性好,且原子时走势可预测性强,这有利于建立准确的原子时预报模型。综合各类原子钟性能及脉冲星时特点,铯原子钟的中长期稳定特性(月稳)优于氢钟,且主要表现为频率漂移。最适合用于被脉冲星时信号驯服。
脉冲星时特点是长期稳定性(年稳)好,产生时间可长期自主保持。但目前脉冲信号到达时间(Time of Arrive,TOA)测量精度低,为百纳秒甚至微秒量级,原子钟测量精度已优于纳秒。利用测量噪声大的脉冲星时驾驭铯原子钟,至少需要3个月以上的脉冲星时与原子钟比对数据,才能计算得到原子钟频率漂移量。目前原子钟驾驭采用策略是原子钟建模数据与钟预报数据跨度一致,国家时间基准保持实验室的原子钟驾驭策略是建模数据跨度和预报时长均为30天。如果直接利用传统原子钟驾驭策略,基于30天跨度的脉冲星观测数据对原子钟进行建模时,由于脉冲星TOA测量误差大,脉冲星时与原子钟比对数据误差大,则利用钟差比对数据计算的原子钟频率偏差值测量误差大,从而导致驾驭后原子钟时频信号准确度低。若利用脉冲星时长期稳特性好特点,以6个月以上时间跨度数据对原子钟建模,基于钟模型预报时长也为6个月以上,则由于预报时间太长而影响最终的时频信号精准度。如何有效利用脉冲星时长期稳定度高特点,又降低TOA测量误差大对原子钟频率漂移量解算精度的影响是脉冲星时应用的关键。
发明内容
本发明所要解决的技术问题在于克服现有技术的缺点,提供一种预报时间短、测量误差小、识别精度高、速度快的基于滑动窗口的脉冲星时驯服铯原子钟的方法。
解决上述技术问题所采用的技术方案是:一种基于滑动窗口的脉冲星时驯服铯原子钟方法,由下步骤组成:
S1.脉冲星计时观测
以铯原子钟为参考对一颗毫秒脉冲星进行计时观测,并采集观测数据,每天观测1次,每次观测时长为30~60分钟;
S2.脉冲星计时数据处理
当观测数据时间跨度达到活动窗口大小时,利用脉冲星专业数据处理软件psrchive对观测数据进行处理,获得一系列脉冲到达时间数据序列,然后结合最新脉冲星星历表,利用计时专业数据处理软件Tempo2对脉冲到达时间数据序列进行计时分析,获得脉冲星时与原子钟的钟差序列;
S3.利用活动窗口计算铯原子钟频率偏差
利用滑动窗口内脉冲星时与铯原子钟时的钟差序列,按铯原子钟模型公式(1)进行最小二乘拟合,获得铯原子钟频率偏差k;
基于脉冲星时观测数据建立的铯原子钟模型为:
x(t)=x0+(t-t0)·k+δAT+δPT (1)
式中,x(t)是铯原子钟相对于脉冲星时钟差,x0为t0时刻铯原子钟相对于脉冲星时初始钟差,k为铯原子钟频率偏差,δAT为铯原子钟随机噪声,δPT为脉冲星时噪声;
S4.铯原子钟频率校准
利用相位微调仪对铯原子钟频率进行驾驭,驾驭频度为1次/窗口滑动步长,即当观测数据每增加一个窗口滑动步长,则以当前时刻为终点,按滑动窗口大小提取观测数据,计算更新铯原子钟频率偏差k,并利用相位微调仪对铯原子钟频率进行驾驭;所述窗口滑动步长是滑动窗口的1/4~1/3倍。
作为一种优选的技术方案,所述的步骤S2中活动窗口大小为90天~180天。
作为一种优选的技术方案,所述的步骤S3中脉冲星时噪声δPT选取白噪声的源。
本发明的有益效果如下:
本发明提出基于滑动窗口技术的铯原子钟相对于脉冲星钟频率偏差测量方法,解决了脉冲星脉冲到达时间测量误差大,对原子钟频率漂移测量精度影响,在窗口滑动步长不变的情况下,增加滑动窗口大小,有效提高了频率偏差测量精度。本发明将为基于脉冲星时实现原子钟的高精度频率溯源提供支撑,具有预报时间短、测量误差小、识别精度高、速度快的优点。
附图说明
图1是利用脉冲星测量原子钟误差原理图。
图2是本发明基于滑动窗口的脉冲星时驯服铯原子钟方法的流程图。
图3是滑动窗口的示意图。
具体实施方式
下面结合附图和实施例对本发明进一步详细说明,但本发明不限于下述的实施方式。
本实施例以基于滑动窗口利用脉冲星J0437-4715观测数据驯服铯原子钟的方法为例,对基于滑动窗口的脉冲星时驯服铯原子钟的方法进行说明,如图2和图3,具体步骤如下:
S1.脉冲星计时观测
利用地面射电望远镜脉冲星观测***,以铯原子钟为参考对脉冲星J0437-4715进行计时观测,并采集观测数据,每天观测1次,每次观测时长为60分钟;
S2.脉冲星计时数据处理
当观测数据时间跨度达到活动窗口大小即时,其中活动窗口大小为90天,也可以是180天,利用脉冲星专业数据处理软件psrchive对观测数据进行处理,获得一系列脉冲到达时间数据序列,然后结合最新脉冲星星历表,利用计时专业数据处理软件Tempo2对脉冲到达时间数据序列进行计时分析,获得脉冲星时与铯原子钟时的钟差序列;
上述计时分析时,不对脉冲星星历表参数进行拟合;
S3.利用活动窗口计算铯原子钟频率偏差
利用滑动窗口内脉冲星时与铯原子钟时的钟差序列,按铯原子钟模型公式(1)进行最小二乘拟合,获得铯原子钟频率偏差k;
基于脉冲星时观测数据建立的铯原子钟模型为:
x(t)=x0+(t-t0)·k+δAT+δPT (1)
式中,x(t)是铯原子钟相对于脉冲星时钟差,x0为t0时刻铯原子钟相对于脉冲星时初始钟差,k为铯原子钟频率偏差,δAT为铯原子钟随机噪声,δPT为脉冲星计时噪声,在候选观测脉冲星时优先δPT为白噪声的源,有助于降低脉冲星时噪声对铯原子钟频率偏差k拟合的影响;
S4.铯原子钟频率校准
利用相位微调仪对铯原子钟频率进行驾驭,驾驭频度为1次/窗口滑动步长,即当观测数据每增加一个窗口滑动步长,则以当前时刻为终点,按滑动窗口大小提取观测数据,计算更新铯原子钟频率偏差k,并利用相位微调仪对铯原子钟进行频率驾驭,所述窗口滑动步长为30天。
下面利用中科院国家授时中心自由运转的铯原子钟Cs2147与国家标准时间UTC(NTSC)实测钟差数据为参考,分析基于滑动窗口技术利用脉冲星时驾驭铯原子钟效果。首先分析利用UTC(NTSC)驾驭Cs2147钟的精度,原子钟建模数据跨度和预报时长均为30天。驾驭策略及结果见表1,驾驭后的Cs2147相对于UTC(NTSC)钟差RMS为18.86ns,极差值即钟差数据的极大值与极小值之差为107.58ns。
表1.利用UTC(NTSC)驾驭Cs2147结果
下面以目前计时精度最高的脉冲星J0437-4715为基准,开展基于滑动窗口的脉冲星时驾驭铯原子钟Cs2147精度模拟分析。利用专业脉冲计时软件Tempo2模拟昊平40米天线观测J0437-4715的计时数据,观测参考钟为Cs2147,模拟数据参数设置如下:观测频次为1次/天,脉冲到达时间TOA测量误差为0.1微秒。模拟观测数据跨度为铯原子钟Cs2147与UTC(NTSC)钟比对数据跨度。
为了比较基于滑动窗口技术驾驭铯原子钟效果,利用模拟脉冲星数据对Cs2147进行驾驭时,按以下两种方法分析:
1)与上面UTC(NTSC)驾驭Cs2147策略一致;
2)采用滑动窗口技术,设置滑动窗口大小为90天,窗口滑动步长为30天。下表为利用脉冲星J0437-4715驾驭铯原子钟Cs2147的结果。
表2.基于滑动窗口技术利用脉冲星时驾驭Cs2147结果
根据上表可得知,采用传统驾驭原子钟技术,当钟建模数据跨度与预报时长一致即1个月,利用脉冲星时驾驭后的Cs2147相对于UTC(NTSC)极差值为246.93ns,钟差RMS为44.72ns。利用滑动窗口技术,驾驭后的Cs2147相对于UTC(NTSC)极差值为103.5ns,RMS为17.58ns。根据结果得知,利用传统驾驭策略,基于脉冲星时驾驭后的原子钟精度,明显低于采用UTC(NTSC)为参考驾驭原子钟结果。当采用滑动窗口技术时,其驾驭结果略优于传统原子钟驾驭方法精度。
Claims (3)
1.一种基于滑动窗口的脉冲星时驯服铯原子钟方法,其特征在于,由下步骤组成:
S1.脉冲星计时观测
以铯原子钟为参考对一颗毫秒脉冲星进行计时观测,并采集观测数据,每天观测1次,每次观测时长为30~60分钟;
S2.脉冲星计时数据处理
当观测数据时间跨度达到活动窗口大小时,利用脉冲星专业数据处理软件psrchive对观测数据进行处理,获得一系列脉冲到达时间数据序列,然后结合最新脉冲星星历表,利用计时专业数据处理软件Tempo2对脉冲到达时间数据序列进行计时分析,获得脉冲星时与原子钟的钟差序列;
S3.利用活动窗口计算铯原子钟频率偏差
利用滑动窗口内脉冲星时与铯原子钟时的钟差序列,按铯原子钟模型公式(1)进行最小二乘拟合,获得铯原子钟频率偏差k;
基于脉冲星时观测数据建立的铯原子钟模型为:
x(t)=x0+(t-t0)·k+δAT+δPT (1)
式中,x(t)是铯原子钟相对于脉冲星时钟差,x0为t0时刻铯原子钟相对于脉冲星时初始钟差,k为铯原子钟频率偏差,δAT为铯原子钟随机噪声,δPT为脉冲星时噪声;
S4.铯原子钟频率校准
利用相位微调仪对铯原子钟频率进行驾驭,驾驭频度为1次/窗口滑动步长,即当观测数据每增加一个窗口滑动步长,则以当前时刻为终点,按滑动窗口大小提取观测数据,计算更新铯原子钟频率偏差k,并利用相位微调仪对铯原子钟频率进行驾驭;所述窗口滑动步长是滑动窗口的1/4~1/3倍。
2.根据权利要求1所述的基于滑动窗口的脉冲星时驯服铯原子钟方法,其特征在于:所述的步骤S2中活动窗口大小为90天~180天。
3.根据权利要求1所述的基于滑动窗口的脉冲星时驯服铯原子钟方法,其特征在于:所述的步骤S3中脉冲星时噪声δPT选取白噪声的源。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111242954.XA CN113985719B (zh) | 2021-10-25 | 2021-10-25 | 基于滑动窗口的脉冲星时驯服铯原子钟的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111242954.XA CN113985719B (zh) | 2021-10-25 | 2021-10-25 | 基于滑动窗口的脉冲星时驯服铯原子钟的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113985719A true CN113985719A (zh) | 2022-01-28 |
CN113985719B CN113985719B (zh) | 2022-09-16 |
Family
ID=79741185
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111242954.XA Active CN113985719B (zh) | 2021-10-25 | 2021-10-25 | 基于滑动窗口的脉冲星时驯服铯原子钟的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113985719B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117014007A (zh) * | 2023-09-07 | 2023-11-07 | 湖南时空信安科技有限公司 | 钟差驾驭方法、装置及终端设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050192719A1 (en) * | 2003-12-08 | 2005-09-01 | Suneel Ismail Sheikh | Navigational system and method utilizing sources of pulsed celestial radiation |
WO2012062207A1 (zh) * | 2010-11-11 | 2012-05-18 | 国网电力科学研究院 | 基于铷振荡器的标准频率与时间调整方法 |
CN106383438A (zh) * | 2016-11-14 | 2017-02-08 | 南京音视软件有限公司 | 一种基于滑动窗口时间扩展高精度驯钟方法 |
CN108196267A (zh) * | 2017-12-20 | 2018-06-22 | 中国科学院国家授时中心 | 一种基于gnss cp技术的不间断时间传递方法 |
CN110780588A (zh) * | 2019-10-16 | 2020-02-11 | 北京航空航天大学 | 一种广域精确授时wpt***和方法 |
CN113238474A (zh) * | 2021-05-06 | 2021-08-10 | 中国科学院国家授时中心 | 脉冲星时信号地面复现*** |
-
2021
- 2021-10-25 CN CN202111242954.XA patent/CN113985719B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050192719A1 (en) * | 2003-12-08 | 2005-09-01 | Suneel Ismail Sheikh | Navigational system and method utilizing sources of pulsed celestial radiation |
WO2012062207A1 (zh) * | 2010-11-11 | 2012-05-18 | 国网电力科学研究院 | 基于铷振荡器的标准频率与时间调整方法 |
CN106383438A (zh) * | 2016-11-14 | 2017-02-08 | 南京音视软件有限公司 | 一种基于滑动窗口时间扩展高精度驯钟方法 |
CN108196267A (zh) * | 2017-12-20 | 2018-06-22 | 中国科学院国家授时中心 | 一种基于gnss cp技术的不间断时间传递方法 |
CN110780588A (zh) * | 2019-10-16 | 2020-02-11 | 北京航空航天大学 | 一种广域精确授时wpt***和方法 |
CN113238474A (zh) * | 2021-05-06 | 2021-08-10 | 中国科学院国家授时中心 | 脉冲星时信号地面复现*** |
Non-Patent Citations (2)
Title |
---|
孙大双等: "北斗在轨卫星钟产品质量分析", 《测绘科学技术学报》 * |
赵齐乐等: "GNSS实时精密轨道快速计算方法及服务", 《武汉大学学报(信息科学版)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117014007A (zh) * | 2023-09-07 | 2023-11-07 | 湖南时空信安科技有限公司 | 钟差驾驭方法、装置及终端设备 |
CN117014007B (zh) * | 2023-09-07 | 2024-01-05 | 湖南时空信安科技有限公司 | 钟差驾驭方法、装置及终端设备 |
Also Published As
Publication number | Publication date |
---|---|
CN113985719B (zh) | 2022-09-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rivera et al. | A refined astronomically calibrated 40Ar/39Ar age for Fish Canyon sanidine | |
CN113078901B (zh) | 基于脉冲星驾驭的原子钟准确度提高方法 | |
CN109633701B (zh) | 基于守时实验室资源的gnss定时接收机***延迟校准方法 | |
Rodriguez-Trelles et al. | Time-series analysis of seasonal changes of the O inversion polymorphism of Drosophila subobscura | |
CN109581856B (zh) | 一种基于高性能晶振频率校准的对时守时方法 | |
CN105388445B (zh) | 一种单相电能表的自动校表方法 | |
CN113985719B (zh) | 基于滑动窗口的脉冲星时驯服铯原子钟的方法 | |
CN202008583U (zh) | 一种同步相量测量装置时钟源 | |
CN103176400A (zh) | 智能电表时钟校准方法 | |
CN108107233B (zh) | 加速度计标度因数的连续温度校正方法及*** | |
CN110837221B (zh) | 一种有效提高授时可靠性与连续性的方法 | |
CN110515294A (zh) | 一种嵌入式gps终端的快速高精度授时方法 | |
CN116125512A (zh) | 顾及时钟频率时变特性的ppp自适应钟差模型估算方法 | |
CN106383438B (zh) | 一种基于滑动窗口时间扩展高精度驯钟方法 | |
CN201540331U (zh) | 多通道高精度同步测频装置 | |
CN102928805A (zh) | 多表位数字化电能表校验仪 | |
CN104460313A (zh) | 用于提供高精度大步长时间同步信号的gps授时装置 | |
CN103278128B (zh) | 一种快速的kdp晶体最佳匹配角精确测量方法 | |
CN116627015A (zh) | 基于小波分析法的脉冲星时与原子时融合方法 | |
Gojcic et al. | Synchronization routine for real-time synchronization of robotic total stations | |
CN211180593U (zh) | 时频检测设备 | |
CN101464822B (zh) | 一种cpu定时器准确性的测试方法 | |
CN110187198B (zh) | 一种频率器件性能评估的方法及装置 | |
CN110908272B (zh) | 一种1pps脉冲信号校时方法 | |
CN110658536B (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 |