CN113078901B - 基于脉冲星驾驭的原子钟准确度提高方法 - Google Patents
基于脉冲星驾驭的原子钟准确度提高方法 Download PDFInfo
- Publication number
- CN113078901B CN113078901B CN202110337703.3A CN202110337703A CN113078901B CN 113078901 B CN113078901 B CN 113078901B CN 202110337703 A CN202110337703 A CN 202110337703A CN 113078901 B CN113078901 B CN 113078901B
- Authority
- CN
- China
- Prior art keywords
- pulsar
- timing
- time
- atomic clock
- toa
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 8
- 239000006185 dispersion Substances 0.000 claims description 3
- MRMBZHPJVKCOMA-YJFSRANCSA-N biapenem Chemical compound C1N2C=NC=[N+]2CC1SC([C@@H]1C)=C(C([O-])=O)N2[C@H]1[C@@H]([C@H](O)C)C2=O MRMBZHPJVKCOMA-YJFSRANCSA-N 0.000 abstract 1
- 229910052792 caesium Inorganic materials 0.000 description 4
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 230000007123 defense Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/26—Automatic control of frequency or phase; Synchronisation using energy levels of molecules, atoms, or subatomic particles as a frequency reference
-
- 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)
- Spectroscopy & Molecular Physics (AREA)
- Electric Clocks (AREA)
Abstract
本发明提供了一种基于脉冲星驾驭的原子钟准确度提高方法，获得一段时间的脉冲星TOA观测数据，拟合计时模型参数，获得拟合后计时残差和更新后的计时模型参数；利用要被驾驭的原子钟作为参考时继续开展脉冲星计时观测，得到由该参考钟记录的实测TOA数据；在该观测时段内利用计时模型参数在SSB处预报TOA，结合原子钟实测的TOA序列，得到拟合前计时残差；将拟合前计时残差数据做线性拟合，得到斜率k；将被驾驭的原子钟与TT(BIPM)的钟差补偿一个(1‑k)因子，完成频率驾驭。本发明能够实时地对原子钟进行频率驾驭，提高原子钟的准确度，并在UTC中断时保持我国标准时间持续高性能。
Description
技术领域
本发明属于脉冲星计时应用的技术领域,涉及脉冲星钟模型参数建立，以及如何利用自转高度稳定的脉冲星对原子钟实现频率驾驭的方法。
背景技术
一台守时型铯原子钟的输出频率通常会偏离其标称频率，这会导致它与标准时间的钟差越来越大，因此需要外部的频率基准去校准其频率，提高其准确度。一个守时实验室中的多台原子钟建立的综合原子时尽管缩小了相对于标准时间的频率偏差，但偏差仍然是不可忽略的，尤其是在各行各业对时间精度需求日益提高的背景下。而标准时间通常由全世界的守时实验室的500多台原子钟综合得到的，即国际原子时(TAI)，或者是与其只有闰秒之差的协调世界时(UTC)，每月发布一次。当前我国标准时间保持体系基于国际权度局(BIPM)每月公布的UTC(NTSC)相对于UTC的偏差和原子钟的相关参数，其可靠性和准确性过度依赖国际组织，并且实时性较差。尤其当比对链路中断，地方原子时UTC(k)的准确度将不受保障。此外，依赖于地面原子钟建立的原子时，易受打击，抗摧毁能力差，特别是非常时期，一旦中断，便无法恢复，难以满足社会稳定、国防安全的需求。
脉冲星是宇宙中存在的一类致密天体，沿磁极辐射电磁波束，其自转非常稳定。当电磁波束扫过地球时，人们可以观测到稳定的周期性的脉冲信号，这种稳定的周期脉冲信号形成一种天然的频率源，就像是处于宇宙深处的原子钟，称为脉冲星钟。尽管单脉冲的变化比较大，但经几百到几千个单脉冲叠加后可获得稳定的脉冲轮廓，脉冲轮廓尖峰对应的到达时间被定义为脉冲达到时间(TOA)，可以被原子钟记录。对于毫秒脉冲星而言，TOA的测量精度很高，可以达到百纳秒量级甚至更小。经过长期的计时观测后，人们可以通过TOA数据拟合得到毫秒脉冲星的自转频率、自转频率变化率、位置、自行、视差等一系列计时模型参数。其中脉冲星的自转频率及其变化率实际上就是脉冲星钟的钟模型参数。由于毫秒脉冲星的自转频率变化率很小，因此短时间内脉冲星的自转频率非常稳定，构成一种频率基准，而这种频率基准可以替代TAI去驾驭原子钟，提高其准确度。与原子钟相比，脉冲星是自然天体，具有寿命长、可靠性高、服务范围广，以及不受环境和人为等因素的影响等优势。如果能密集地进行脉冲星计时观测，则可以实现比TAI实时性更强的频率驾驭，以此可以提供更高精度的实时性时间服务。作为独立于UTC的时间保持体系，能够解决非常时期UTC中断时我国标准时间持续高性能保持的问题，对我国国防安全和社会稳定具有重要的战略意义。
在长期的脉冲星计时观测中，通常将参考原子时溯源到由BIPM发布的地球时(TT(BIPM))。TT(BIPM)是在TAI实现的地球时TT(TAI)的基础上由BIPM再次数据处理得到，每年发布一次，是目前通过原子钟得到的精度最高的时间尺度。脉冲星的计时模型参数是通过多次迭代拟合TOA数据得到的。首先根据初始的计时模型参数在太阳系质心(SSB)处预报TOA，利用参考原子时实际测量的TOA减去模型预报的TOA便得到拟合前脉冲星计时残差。然后通过最小二乘法多次迭代拟合直至收敛，得到拟合后计时残差和更新后的脉冲星计时模型参数。因此当有了新的观测数据，又会重复上述过程再次更新脉冲星的计时模型参数。脉冲星计时过程实际上是参考原子钟和脉冲星钟时间比对的过程，因此计时残差实际上代表了参考原子钟和脉冲星钟的钟差。对毫秒脉冲星来说，以TT(BIPM)作为参考原子时的拟合后计时残差往往呈现出均值接近零的白噪声特征。噪声的均方根通常在微秒或百纳秒量级，主要是由脉冲星TOA的测量噪声引起的。随着FAST大型射电望远镜对毫秒脉冲星的计时观测，白噪声水平有望达到10纳秒量级。如果新观测的TOA数据时间跨度不是很长，例如1年，那么这1年的拟合前计时残差依然呈现均值为零的白噪声特征。这主要是由于经多年计时观测建立的脉冲星钟的稳定性很高，计时模型参数拟合前后变化很小，即使是影响最大的自转频率变化也很小。但是如果计时观测中的参考原子时是一个普通的原子钟，由于原子钟相对于TT(BIPM)有频率偏差，那么拟合前计时残差会有相同的频率偏差。
发明内容
为了克服现有技术的不足，本发明提供一种基于脉冲星驾驭原子钟频率的方法，利用脉冲星稳定的自转频率来驾驭原子钟频率，能够实时地对原子钟进行频率驾驭，提高原子钟的准确度，并在UTC中断时保持我国标准时间持续高性能。
本发明解决其技术问题所采用的技术方案包括以下步骤：
1)获得一段时间的脉冲星TOA观测数据，拟合计时模型参数，获得拟合后计时残差和更新后的计时模型参数；
2)利用要被驾驭的原子钟作为参考时继续开展脉冲星计时观测，得到由该参考钟记录的实测TOA数据；在该观测时段内，利用步骤1)得到的计时模型参数在SSB处预报TOA，结合原子钟实测的TOA序列，得到拟合前计时残差；
3)将拟合前计时残差数据做线性拟合，得到斜率k；
4)将被驾驭的原子钟与TT(BIPM)的钟差补偿一个(1-k)因子，完成频率驾驭。
所述的步骤1)采用BIPM发布的地球时TT(BIPM)作为参考时获得一段时间的脉冲星TOA观测数据。
所述的计时模型参数包括脉冲星的自转频率及其一阶导数、位置、自行、视差、星际介质色散和双星轨道模型参数。
本发明的有益效果是：如果计时观测中的参考原子时是一个普通的原子钟，由于原子钟相对于TT(BIPM)有频率偏差，那么拟合前计时残差会有相同的频率偏差。本发明基于这一特点，利用脉冲星稳定的自转频率对参考原子钟实施频率驾驭，可以改正原子钟的输出频率，提高其准确度，从而提高其时间信号的准确性。
附图说明
图1是脉冲星频率驾驭原子钟流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明主要利用被驾驭的原子钟作为参考时时间来进行脉冲星计时，根据拟合前计时残差来驾驭原子钟频率，具体包括以下步骤：
1)首先利用高精度的原子时作为参考时获得一段时间的脉冲星TOA观测数据，这里的高精度原子时采用BIPM发布的地球时TT(BIPM)，因此需要通过高精度时间传递将计时的本地钟溯源到TT(BIPM)。然后拟合计时模型参数，包括脉冲星的自转频率及其一阶导数、位置、自行、视差、星际介质色散(DM)、双星轨道模型参数等，获得拟合后计时残差和更新后的计时模型参数。
2)利用要被驾驭的原子钟作为参考时继续开展脉冲星计时观测，得到由该参考钟记录的实测TOA数据。在这段观测的时段内，利用步骤1)得到的计时模型参数在SSB处预报TOA，结合原子钟实测的TOA序列，得到拟合前计时残差，不再做计时模型参数拟合，这实际上代表了参考原子钟和脉冲星钟的钟差。
3)将拟合前计时残差数据做线性拟合，得到斜率k。斜率k代表了被驾驭的原子钟相对于脉冲星钟的相对频率偏差。
4)将被驾驭的原子钟与TT(BIPM)的钟差补偿一个(1-k)因子，则原子钟与TT(BIPM)之差的斜率趋近于零，完成频率驾驭。
下面利用两颗毫秒脉冲星的仿真数据来驾驭一台原子钟(中国科学院国家授时中心的Cs2098)为例来对本发明进一步说明，本发明包括但不仅限于下述实例。
(1)首先以TT(BIPM2018)为时间基准并将其作为脉冲星计时的参考时间，基于毫秒脉冲星J0437-4715的计时模型参数分别仿真其TOA数据，数据时间跨度均为3年(MJD：54834-55924)，TOA数据均只含有白噪声，均方根(RMS)分别100ns和10ns，采样频率均为每5天一个数据点。实际情况下，毫秒脉冲星的计时模型参数可通过长时间的TOA实测数据拟合而获得，TOA观测数据越长，计时模型参数越精确。
(2)截取最后一年的模拟数据作为检验数据，即前两年的数据参考时为TT(BIPM2018)，利用脉冲星计时数据处理软件将这两年的数据重新做计时模型参数拟合，得到一组更新后的计时模型参数，并以此在SSB处往后预报TOA。
而最后一年的数据(MJD：55559-55924)参考时为铯钟Cs2098，得到Cs2098记录的TOA。由于脉冲星时是在太阳系质心(SSB)处预报TOA，需要将Cs2098在地球记录的TOA转换到SSB处。
(3)将步骤(2)中SSB处Cs2098记录的TOA与更新后计时模型预报的TOA做差，即得到铯原子钟与脉冲星钟的钟差。分别讨论三种时间长度的钟差情况，即一年、半年和一个月。将未经驾驭的Cs2098的输出时间记为AT，将基于PSR J0437-4715的脉冲星时记为PT，二者之差为：
AT-PT
这个钟差序列在数据处理中可由脉冲星拟合前计时残差获得，只需将参考时设为Cs2098。在把TT(BIPM2018)替换为AT的时刻(MJD＝55559)，利用AT-TT(BIPM2018)的历史数据将AT补偿一个钟差，使其校准至TT(BIPM2018)。
(4)将步骤(3)得到的AT-PT做线性拟合得到斜率k，k代表了AT偏离于标准时间信号的相对频率偏差，k的误差决定了频率驾驭精度。将同一时间段AT-TT(BIPM2018)的斜率记为k0。二者的相对偏差代表了频率驾驭的***误差。
(5)将k作为输入量利用相位微调仪将铯原子钟的输出频率由原来的f0调至f0(1-k)，达到驾驭的目的。将脉冲星驾驭后铯钟输出的时间信号记为ATsteer，则
ATsteer(j)＝AT(1)+[AT(j)-AT(1)](1-k)
其中j表示从MJD＝55560起算的第j个时间序列，因此经脉冲星频率驾驭后的Cs2098与TT(BIPM2018)的钟差为
[ATsteer-TT(BIPM2018)](j)＝[AT-TT(BIPM2018)](j)+[AT(j)-AT(1)](-k)
完成频率驾驭过程。表1列出了利用脉冲星时的驾驭量与目标驾驭量(AT与TT(BIPM2018)本身的相对频率偏差)之间的比较结果。
表1驾驭量与目标驾驭量之间的比较结果
为了增加显示度，利用同一时间段的TAI数据计算的驾驭量也列在了表中。可以看出，对于1年和半年而言，PT的驾驭量比TAI的驾驭量更接近目标驾驭量；而1个月的情况，TAI的驾驭量比PT的驾驭量更接近目标驾驭量。另一方面，TOA的测量噪声水平严重影响PT驾驭量的计算，如果能进一步减小TOA的测量误差，那么在1个月的尺度上PT的频率驾驭水平有望会超过TAI。
Claims (2)
1.一种基于脉冲星驾驭的原子钟准确度提高方法，其特征在于，包括以下步骤：
1)获得一段时间的脉冲星脉冲达到时间TOA观测数据，拟合计时模型参数，获得拟合后计时残差和更新后的计时模型参数；
所述的计时模型参数包括脉冲星的自转频率及其一阶导数、位置、自行、视差、星际介质色散和双星轨道模型参数；
2)利用要被驾驭的原子钟作为参考时继续开展脉冲星计时观测，得到由该参考钟记录的实测脉冲达到时间TOA数据；在该观测时段内，利用步骤1)得到的计时模型参数在太阳系质心SSB处预报脉冲达到时间TOA，结合原子钟实测的脉冲达到时间TOA序列，得到拟合前计时残差；
3)将拟合前计时残差数据做线性拟合，得到斜率k；
4)将被驾驭的原子钟与BIPM发布的地球时TT(BIPM)的钟差补偿一个(1-k)因子，完成频率驾驭。
2.根据权利要求1所述的基于脉冲星驾驭的原子钟准确度提高方法，其特征在于，所述的步骤1)采用BIPM发布的地球时TT(BIPM)作为参考时获得一段时间的脉冲星脉冲达到时间TOA观测数据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110337703.3A CN113078901B (zh) | 2021-03-30 | 2021-03-30 | 基于脉冲星驾驭的原子钟准确度提高方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110337703.3A CN113078901B (zh) | 2021-03-30 | 2021-03-30 | 基于脉冲星驾驭的原子钟准确度提高方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113078901A CN113078901A (zh) | 2021-07-06 |
| CN113078901B true CN113078901B (zh) | 2024-02-13 |
Family
ID=76611414
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110337703.3A Active CN113078901B (zh) | 2021-03-30 | 2021-03-30 | 基于脉冲星驾驭的原子钟准确度提高方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113078901B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113987803B (zh) * | 2021-10-29 | 2022-05-13 | 星汉时空科技(北京)有限公司 | 基于模型修正的频率驾驭方法、装置、设备及存储介质 |
| CN114035418B (zh) * | 2021-11-01 | 2022-11-15 | 中国科学院国家授时中心 | 基于脉冲星的共视时间比对方法 |
| CN114114883B (zh) * | 2021-11-01 | 2022-12-23 | 中国科学院国家授时中心 | 利用观测多颗脉冲星检测原子时钟跳的方法 |
| CN114002933B (zh) * | 2021-11-01 | 2022-10-21 | 中国科学院国家授时中心 | 基于脉冲星搜寻技术测量原子钟频率漂移的方法 |
| CN117014007B (zh) * | 2023-09-07 | 2024-01-05 | 湖南时空信安科技有限公司 | 钟差驾驭方法、装置及终端设备 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759885A (zh) * | 2012-07-13 | 2012-10-31 | 西安交通大学 | 一种高可靠性的脉冲星时间同步装置 |
| CN103017774A (zh) * | 2012-12-24 | 2013-04-03 | 中国科学院国家授时中心 | 单探测器脉冲星导航方法 |
| CN103674032A (zh) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | 融合脉冲星辐射矢量和计时观测的卫星自主导航***及方法 |
| CN109100927A (zh) * | 2018-07-25 | 2018-12-28 | 中国科学院国家授时中心 | 无溯源参考条件下时间统一***时间自主同步方法 |
| CN110780583A (zh) * | 2019-10-29 | 2020-02-11 | 中国科学院国家天文台 | 一种月基脉冲星时间基准生成*** |
| CN111142358A (zh) * | 2020-01-20 | 2020-05-12 | 山东科技大学 | 一种基于自然x射线源观测的深空自主守时方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7034695B2 (en) * | 2000-12-26 | 2006-04-25 | Robert Ernest Troxler | Large area position/proximity correction device with alarms using (D)GPS technology |
-
2021
- 2021-03-30 CN CN202110337703.3A patent/CN113078901B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759885A (zh) * | 2012-07-13 | 2012-10-31 | 西安交通大学 | 一种高可靠性的脉冲星时间同步装置 |
| CN103674032A (zh) * | 2012-09-04 | 2014-03-26 | 西安电子科技大学 | 融合脉冲星辐射矢量和计时观测的卫星自主导航***及方法 |
| CN103017774A (zh) * | 2012-12-24 | 2013-04-03 | 中国科学院国家授时中心 | 单探测器脉冲星导航方法 |
| CN109100927A (zh) * | 2018-07-25 | 2018-12-28 | 中国科学院国家授时中心 | 无溯源参考条件下时间统一***时间自主同步方法 |
| CN110780583A (zh) * | 2019-10-29 | 2020-02-11 | 中国科学院国家天文台 | 一种月基脉冲星时间基准生成*** |
| CN111142358A (zh) * | 2020-01-20 | 2020-05-12 | 山东科技大学 | 一种基于自然x射线源观测的深空自主守时方法 |
Non-Patent Citations (1)
| Title |
|---|
| 脉冲星时间开环驾驭原子钟的计量技术标准化探讨;董江;;计量学报(第02期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113078901A (zh) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113078901B (zh) | 基于脉冲星驾驭的原子钟准确度提高方法 | |
| CN109525351A (zh) | 一种用于与时间基准站实现时间同步的设备 | |
| Guinot et al. | Atomic time-keeping from 1955 to the present | |
| CN108387915B (zh) | 一种融合钟差预报方法 | |
| CN109633701A (zh) | 基于守时实验室资源的gnss定时接收机***延迟校准方法 | |
| CN113219815B (zh) | 一种基于x射线脉冲星的深空授时方法 | |
| CN107086901A (zh) | 一种bdt建立方法及utc(ntsc)建立方法 | |
| CN113110016A (zh) | 共视数据生成方法、接收机和时间校准*** | |
| McCarthy | Evolution of timescales from astronomy to physical metrology | |
| Norton et al. | Solar cycle observations | |
| CN113985719B (zh) | 基于滑动窗口的脉冲星时驯服铯原子钟的方法 | |
| CN114911153B (zh) | 基于星-地超长基线干涉测量的世界时ut1测定方法 | |
| CN104635193B (zh) | 一种电表用的soc芯片实时时钟高精度补偿方法 | |
| CN103278128A (zh) | 一种快速的kdp晶体最佳匹配角精确测量方法 | |
| CN113341679A (zh) | 基于北斗空间铷原子钟的高精度信号频率控制方法及*** | |
| Essen | Time scales | |
| CN112003583A (zh) | 一种基于温度传感器的rtc自补偿电路 | |
| CN105137753A (zh) | 北斗多功能时统*** | |
| Guinot | Comments on the Terrestrial Pole of Reference, the Origin of the Longitudes, and on the definition of UT1 | |
| Zhao et al. | A new steering strategy for UTC (NTSC) based on hydrogen maser | |
| Lewandowski et al. | A Brief History of UTC Leap Second | |
| CN112462394B (zh) | 一种基于gnss全视比对的远程时间频率配送方法 | |
| CN111025333B (zh) | 基于导航卫星信号短基线定姿的安装偏差标定与修正方法 | |
| CN111123311B (zh) | 星载原子钟调频调相方法 | |
| Hudson | Some characteristics of commonly used time scales |
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 |