CN115128937B - 防欺骗式干扰铷原子钟驯服方法及系统 - Google Patents

防欺骗式干扰铷原子钟驯服方法及系统 Download PDF

Info

Publication number
CN115128937B
CN115128937B CN202211025733.1A CN202211025733A CN115128937B CN 115128937 B CN115128937 B CN 115128937B CN 202211025733 A CN202211025733 A CN 202211025733A CN 115128937 B CN115128937 B CN 115128937B
Authority
CN
China
Prior art keywords
signal
atomic clock
time difference
rubidium atomic
taming
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
Application number
CN202211025733.1A
Other languages
English (en)
Other versions
CN115128937A (zh
Inventor
吴宏硕
杨国文
林杰
王岭
许晨
刘敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN202211025733.1A priority Critical patent/CN115128937B/zh
Publication of CN115128937A publication Critical patent/CN115128937A/zh
Application granted granted Critical
Publication of CN115128937B publication Critical patent/CN115128937B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及铷原子钟驯服技术领域，提供一种防欺骗式干扰铷原子钟驯服方法及系统。该防欺骗式干扰铷原子钟驯服方法，包括如下步骤：对第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号进行分路处理；对分路信号进行频率合成；将分路信合成信号对应转化并输出1PPS信号；监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰，计算1PPS信号与外部输入1PPS信号的时差；对时差进行驯服校准，获得控制数据。本发明实现对双铷原子钟驯服过程中数据的综合判定，能够有效识别出缓慢变化的欺骗式干扰，可有效避免铷原子钟的频率输出被欺骗式干扰影响，保证铷原子钟在驯服期间频率输出的准确性。

Description

防欺骗式干扰铷原子钟驯服方法及系统
技术领域
本发明涉及铷原子钟驯服技术领域，尤其涉及一种防欺骗式干扰铷原子钟驯服方法及系统。
背景技术
目前，铷原子钟驯服方法，通常是基于外部参考时间源信号的输出特性，用来调整本地铷原子钟的频率输出准确度，以达到提高铷原子钟输出频率准确度的目的。
然而，当人为缓慢改变外部参考时间源的输出特性，对铷原子钟的驯服过程进行欺骗式干扰时，如果无法及时识别出干扰信号，铷原子钟的驯服过程会被外部信号欺骗，导致铷原子钟的频率信号输出也将被缓慢拉偏，其输出准确度变差。
发明内容
本发明提供一种防欺骗式干扰铷原子钟驯服方法及系统，用以解决现有技术中在铷原子钟的驯服过程，无法及时识别出干扰信号的缺陷，实现对双铷原子钟驯服过程中数据的综合判定，能够有效识别出缓慢变化的欺骗式干扰，可有效避免铷原子钟的频率输出被欺骗式干扰影响，保证铷原子钟在驯服期间频率输出的准确性。
本发明提供一种防欺骗式干扰铷原子钟驯服方法，包括如下步骤：
S10、控制第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号，对第一信号进行分路处理，输出第一分路信号、第二分路信号和第三分路信号，对第二信号进行分路处理，输出第四分路信号、第五分路信号和第六分路信号；
S20、分别对第一分路信号和第四分路信号进行频率合成，并输出第一合成信号和第二合成信号；
S30、将第二分路信号、第五分路信号、第一合成信号和第二合成信号对应转化并输出第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
S40、以第三分路信号和第六分路信号作为频率输入基准，分别计算第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号与外部输入1PPS信号的时差，对应获得第一时差、第二时差、第三时差和第四时差；
S50、利用第一时差和第二时差监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰，对第三时差和第四时差进行驯服校准，获得控制数据，并将控制数据输入S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S50步骤中，还包括如下步骤：
S51、在第一运行阶段内，确认外部输入1PPS信号无欺骗式干扰；
S52、在第一运行阶段的运行时间内，利用第一时差和第二时差通过时差法获得第一铷原子钟的准确度和第二铷原子钟的准确度；利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S51步骤中，第一运行阶段的运行时间大于或等于2小时。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S50步骤中，还包括如下步骤：
S53、在第一运行阶段结束后，进入第二运行阶段，利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业；
S54、利用第一铷原子钟的准确度、第二运行阶段开始时刻的初始的第一时差、第二运行阶段累计运行时长，获得第一1PPS信号和外部输入1PPS信号的第一预计时差；
S55、利用第二铷原子钟的准确度、第二运行阶段开始时刻的初始的第二时差、第二运行阶段累计运行时长，获得第二1PPS信号和外部输入1PPS信号的第二预计时差；
S56、计算第一时差与第一预计时差的第一相对误差，计算第二时差与第二预计时差的第二相对误差；
S57、通过最小二乘法对第一相对误差和第二相对误差进行曲线拟合，分别获得第一铷原子钟的第一误差曲线斜率和第二铷原子钟的第二误差曲线斜率。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S50步骤中，还包括如下步骤：
S58、将第一误差曲线斜率和第二误差曲线斜率分别与阈值判定，若第一误差曲线斜率和第二误差曲线斜率均大于阈值，则判定外部输入1PPS信号被注入欺骗式干扰；
S59、获得判定外部输入1PPS信号被注入欺骗式干扰的注入时刻，计算注入时刻以后的控制数据的输出总量，完成回调校准并停止驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S50步骤中，还包括如下步骤：
S591、当判断第一误差曲线斜率和第二误差曲线斜率均小于或等于阈值后，判定外部输入1PPS信号欺骗式干扰消失，再次启动驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S54步骤中，第一预计时差的计算式如下：
Figure 976050DEST_PATH_IMAGE002
其中，
Figure 100002_DEST_PATH_IMAGE003
为第一预计时差；
Figure 844780DEST_PATH_IMAGE004
为第二运行阶段开始时刻的初始的第一时差；
Figure 100002_DEST_PATH_IMAGE005
为第一铷原子钟的准确度；
Figure 130268DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S55步骤中，第二预计时差的计算式如下：
Figure 208558DEST_PATH_IMAGE008
其中，
Figure 100002_DEST_PATH_IMAGE009
为第二预计时差；
Figure 263102DEST_PATH_IMAGE010
为第二运行阶段开始时刻的初始的第二时差；
Figure 100002_DEST_PATH_IMAGE011
为第二铷原子钟的准确度；
Figure 365050DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，所述S56步骤中，第一相对误差的计算式如下：
Figure 100002_DEST_PATH_IMAGE013
其中，
Figure 888566DEST_PATH_IMAGE014
为第一相对误差；
Figure 100002_DEST_PATH_IMAGE015
为S40步骤中实际测量获得的第一时差；
第二相对误差的计算式如下：
Figure 100002_DEST_PATH_IMAGE017
其中，
Figure 39056DEST_PATH_IMAGE018
为第二相对误差；
Figure 100002_DEST_PATH_IMAGE019
为S40步骤中实际测量获得的第二时差。
本发明还提供一种防欺骗式干扰铷原子钟驯服系统，用以执行如上所述的防欺骗式干扰铷原子钟驯服方法，包括第一铷原子钟、第二铷原子钟、分路单元、DDS单元、信号转化单元、时差测量处理单元和驯服单元；
其中，所述第一铷原子钟和所述第二铷原子钟的输出端分别与所述分路单元的输入端连接，所述分路单元的第一输出端和第四输出端均与所述DDS单元的输入端连接，所述分路单元的第二输出端和第五输出端以及所述DDS单元的输出端均与所述信号转化单元的输入端连接，所述信号转化单元还用以接收外部输入1PPS信号；
所述分路单元的第三输出端和第六输出端以及所述信号转化单元的输出端均与所述时差测量处理单元的输入端连接，所述时差测量处理单元的输出端与所述驯服单元的输入端连接；
所述驯服单元的输出端与所述DDS单元的输入端连接。
本发明的有益效果如下：
本发明提供的防欺骗式干扰铷原子钟驯服方法及系统，包括如下步骤：
S10、控制第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号，对第一信号进行分路处理，输出第一分路信号、第二分路信号和第三分路信号，对第二信号进行分路处理，输出第四分路信号、第五分路信号和第六分路信号；
S20、分别对第一分路信号和第四分路信号进行频率合成，并输出第一合成信号和第二合成信号；
S30、将第二分路信号、第五分路信号、第一合成信号和第二合成信号对应转化并输出第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
S40、以第三分路信号和第六分路信号作为频率输入基准，分别计算第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号与外部输入1PPS信号的时差，对应获得第一时差、第二时差、第三时差和第四时差；
S50、利用第一时差和第二时差监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰；对第三时差和第四时差进行驯服校准，获得控制数据，并将控制数据输入S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业，通过利用双铷原子钟的频率特性对外部输入1PPS信号进行监测识别，当外部输入1PPS信号出现欺骗式干扰时，能够有效识别出干扰，同时能够在一定程度上消除欺骗式干扰识别过程中驯服被带偏的影响，可有效避免铷原子钟的频率输出被欺骗式干扰影响，保证铷原子钟在驯服期间频率输出的准确性，具有较高的应用价值。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的防欺骗式干扰铷原子钟驯服系统的结构框式示意图；
图2是本发明提供的防欺骗式干扰铷原子钟驯服方法中监测识别策略的流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明中的附图，对本发明中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
下面结合图1和图2描述本发明的一种防欺骗式干扰铷原子钟驯服方法，包括如下步骤：
S10、控制第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号，对第一信号进行分路处理，输出第一分路信号、第二分路信号和第三分路信号，对第二信号进行分路处理，输出第四分路信号、第五分路信号和第六分路信号；
S20、分别对第一分路信号和第四分路信号进行频率合成，并输出第一合成信号和第二合成信号；
S30、将第二分路信号、第五分路信号、第一合成信号和第二合成信号对应转化并输出第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
S40、以第三分路信号和第六分路信号作为频率输入基准，分别计算第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号与外部输入1PPS信号的时差，对应获得第一时差、第二时差、第三时差和第四时差；
S50、利用第一时差和第二时差监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰；对第三时差和第四时差进行驯服校准，获得控制数据，并将控制数据输入S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S50步骤中，还包括如下步骤：
S51、在第一运行阶段内，确认外部输入1PPS信号无欺骗式干扰；
S52、在第一设定时间内，利用第一时差和第二时差通过时差法获得第一铷原子钟的准确度和第二铷原子钟的准确度；利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S51步骤中，第一运行阶段的时间大于或等于2小时。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S50步骤中，还包括如下步骤：
S53、在第一运行阶段结束后，进入第二运行阶段，利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业；
S54、利用第一铷原子钟的准确度、第二运行阶段开始时刻的初始的第一时差、第二运行阶段累计运行时长，获得第一1PPS信号和外部输入1PPS信号的第一预计时差；
S55、利用第二铷原子钟的准确度、第二运行阶段开始时刻的初始的第二时差、第二运行阶段累计运行时长，获得第二1PPS信号和外部输入1PPS信号的第二预计时差；
S56、计算第一时差与第一预计时差的第一相对误差，计算第二时差与第二预计时差的第二相对误差；
S57、通过最小二乘法对第一相对误差和第二相对误差进行曲线拟合，分别获得第一铷原子钟的第一误差曲线斜率和第二铷原子钟的第二误差曲线斜率。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S50步骤中，还包括如下步骤：
S58、将第一误差曲线斜率和第二误差曲线斜率分别与阈值判定，若第一误差曲线斜率和第二误差曲线斜率均大于阈值，则判定外部输入1PPS信号被注入欺骗式干扰；
S59、获得判定外部输入1PPS信号被注入欺骗式干扰的注入时刻，计算注入时刻以后的控制数据的输出总量，完成回调校准并停止驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S50步骤中，还包括如下步骤：
S591、当判断第一误差曲线斜率和第二误差曲线斜率均小于或等于阈值后，判定外部输入1PPS信号欺骗式干扰消失，再次启动驯服作业。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S54步骤中，第一预计时差的计算式如下：
Figure 479264DEST_PATH_IMAGE002
其中，
Figure 890130DEST_PATH_IMAGE003
为第一预计时差；
Figure 212527DEST_PATH_IMAGE004
为第二运行阶段开始时刻的初始的第一时差；
Figure 104390DEST_PATH_IMAGE005
为第一铷原子钟的准确度；
Figure 71209DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S55步骤中，第二预计时差的计算式如下：
Figure 436332DEST_PATH_IMAGE008
其中，
Figure 590232DEST_PATH_IMAGE009
为第二预计时差；
Figure 285787DEST_PATH_IMAGE010
为第二运行阶段开始时刻的初始的第二时差；
Figure 231746DEST_PATH_IMAGE011
为第二铷原子钟的准确度；
Figure 252923DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
根据本发明提供的一种防欺骗式干扰铷原子钟驯服方法，S56步骤中，第一相对误差的计算式如下：
Figure 894120DEST_PATH_IMAGE013
其中，
Figure 642633DEST_PATH_IMAGE014
为第一相对误差；
Figure 456481DEST_PATH_IMAGE015
为S40步骤中实际测量获得的第一时差；
第二相对误差的计算式如下：
Figure 773193DEST_PATH_IMAGE017
其中，
Figure 291899DEST_PATH_IMAGE018
为第二相对误差；
Figure 594835DEST_PATH_IMAGE019
为S40步骤中实际测量获得的第二时差。
如图1所示，本发明还提供一种防欺骗式干扰铷原子钟驯服系统，用以执行如上的防欺骗式干扰铷原子钟驯服方法，包括第一铷原子钟、第二铷原子钟和时频驯服模块，时频驯服模块包括分路单元、DDS单元、信号转化单元、时差测量处理单元和驯服单元；
其中，第一铷原子钟和第二铷原子钟的输出端分别与分路单元的输入端连接，分路单元的第一输出端和第四输出端均与DDS单元的输入端连接，分路单元的第二输出端和第五输出端以及DDS单元的输出端均与信号转化单元的输入端连接，信号转化单元还用以接收外部输入1PPS信号；
分路单元的第三输出端和第六输出端以及信号转化单元的输出端均与时差测量处理单元的输入端连接，时差测量处理单元的输出端与驯服单元的输入端连接；
驯服单元的输出端与DDS单元的输入端连接。
下面对本发明提供的防欺骗式干扰铷原子钟驯服方法进行详细描述，下文描述的防欺骗式干扰铷原子钟驯服方法与上文描述的防欺骗式干扰铷原子钟驯服系统可相互对应参照。
本发明体用的防欺骗式干扰铷原子钟驯服方法，如图1和图2所示，具体包括如下步骤：
S10、第一铷原子钟输出一路10MHz信号即第一信号至分路单元，记为10MHz_Rb1；第二铷原子钟输出一路10MHz信号即第二信号至分路单元，记为10MHz_Rb2；
其中，10MHz_Rb1输出连接至第一分路子单元，第一分路子单元将10MHz_Rb1分路成3路10MHz信号，分别连接至DDS单元、信号转化单元、时差测量处理单元；10MHz_Rb2输出连接至第二分路子单元，第二分路子单元将10MHz_Rb2分路成3路10MHz信号，也分别连接至DDS单元、信号转化单元、时差测量处理单元；
需要说明的是，本实施例中，分路单元包括第一分路子单元和第二分路子单元。其中，第一分路子单元的输入端与第一铷原子钟的输出端连接，第一分路子单元设置有第一输出端、第二输出端和第三输出端，即上述的分路单元的第一输出端、第二输出端和第三输出端；第二分路子单元的输入端与第二铷原子钟的输出端连接，第二分路子单元设置有第四输出端、第五输出端和第六输出端，即上述的分路单元的第四输出端、第五输出端和第六输出端。
S20、DDS单元的输入端分别与第一分路子单元额第一输出端和第二分路子单元的第四输出端连接，用以接收第一分路信号10MHz_Rb11和第四分路信号10MHz_Rb21，作为频率输入基准；
同时，接收驯服单元输出的控制数据作为DDS单元控制字，实现对第一分路信号10MHz_Rb11和第四分路信号10MHz_Rb21的高精度频率合成，并输出4路合成信号，其中2路合成信号，第一合成信号和第二合成信号输出至信号转化单元，另外2路合成信号连接至对外输出接口；
S30、信号转化单元的输入端均与DDS单元的输出端、第一分路子单元的第二输出端和第二分路子单元的第五输出端连接，引入外部输入源1PPS信号、分路单元输出的第二分路信号10MHz_Rb12和第五分路信号10MHz_Rb22、DDS单元输出的第一合成信号10MHz_DDS1和第二合成信号10MHz_DDS2；
信号转化单元利用第二分路信号10MHz_Rb12和第五分路信号10MHz_Rb22作为频率基准，转化并输出2路与外部基准信号粗同步后的1PPS信号，分别记为第一1PPS信号1PPS_Rb1和第二1PPS信号1PPS_Rb2；
信号转化单元利用第一合成信号10MHz_DDS1和第二合成信号10MHz_DDS2，转化输出2路与外部基准信号粗同步后的1PPS信号，记为第三1PPS信号1PPS_DDS1和第四1PPS信号1PPS_DDS2，将4路1PPS输出信号第一1PPS信号1PPS_Rb1、第二1PPS信号1PPS_Rb2、第三1PPS信号1PPS_DDS1和第四1PPS信号1PPS_DDS2输出至时差测量处理单元；
S40、时差测量处理单元的输入端均与信号转化单元的输出端、第一分路子单元的第三输出端、第二分路子单元的第六输出端连接，接收第三分路信号10MHz_Rb13和第六分路信号10MHz_Rb23，同时接入外部输入1PPS信号以及第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
利用时差法分别测量外部输入1PPS信号和1PPS_Rb1、1PPS_Rb2、1PPS_DDS1和1PPS_DDS2之间的时差，并对时差数据进行卡尔曼滤波处理，得到4组时差数据，分别记为第一时差TERb1、第二时差TERb2、第三时差TEDDS1和第四时差TEDDS2，并输出至驯服单元；
S50、驯服单元的输入端与时差测量处理单元的输出端连接，接收时差测量处理单元输出的时差数据TERb1、TERb2、TEDDS1和TEDDS2，利用TEDDS1，TEDDS2数据进行驯服校准，并输出控制数据至DDS单元开展驯服工作；
利用TERb1，TERb2时差数据对外部输入信号的准确性进行监测，当外部输入源信号被注入欺骗式干扰时，驯服单元采取监测识别策略进行有效识别，完成防欺骗式干扰驯服控制；
其中，欺骗式干扰具体为外部输入1PPS信号的上升沿朝着单一方向缓慢渐变式变化。
其中，如图2所示，S50步骤中监测识别策略的具体流程为：
S51、系统工作开机后进入第一运行阶段，在第一运行阶段内判断外部输入1PPS信号不存在欺骗式干扰时，其中，第一运行阶段时长大于或等于2小时，本实施例中设定为2小时；
S52、在第一运行阶段内，驯服单元利用TERb1、TERb2数据通过时差法，获得第一铷原子钟和第二铷原子钟在第一运行阶段内的频率准确度A1和频率准确度A2；
驯服单元利用TEDDS1、TEDDS2数据通过驯服算法完成DDS单元控制数据输出，并将控制数据输入DDS单元，对第一分路信号和第四分路信号进行校准，开展对铷原子钟的驯服工作；需要说明的是，第一运行阶段的目的是对第一铷原子钟的准确度和第二铷原子钟的准确度重复计算，并与外部参考源进行判断对比，以提高第一铷原子钟的准确度和第二铷原子钟的准确度的准确性；
S53、在第一运行阶段结束后，系统进入第二运行阶段。在第二运行阶段期间，驯服单元利用TEDDS1、TEDDS2数据通过驯服算法完成DDS控制数据输出，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展铷原子钟驯服工作；
需要说明的是，第二运行阶段的目的是对外部输入1PPS信号的检测，并确认外部输入1PPS信号有无欺骗式干扰；
还需要说明的是，在第一运行阶段和第二运行阶段始终对第一分路信号和第四分路信号进行驯服作业，若第二运行阶段检测外部输入1PPS信号存在欺骗式干扰，则停止驯服作业；
S54、驯服单元综合利用第一铷原子钟的准确度A1、第二运行阶段开始时刻的初始 的第以时差
Figure 125174DEST_PATH_IMAGE004
、第二运行阶段累计运行时长
Figure 737421DEST_PATH_IMAGE006
计算获得第二1PPS信号1PPS_Rb1和外部 输入1PPS信号之间的第一预计时差
Figure 228576DEST_PATH_IMAGE003
，具体计算式为：
Figure 459837DEST_PATH_IMAGE020
；
S55、综合利用第二铷原子钟的准确度A2、第二运行阶段开始时刻的初始的第二时 差
Figure 969316DEST_PATH_IMAGE010
、第二运行阶段累计运行时长
Figure 503197DEST_PATH_IMAGE006
，计算获得第二1PPS信号和外部输入1PPS信号的第 二预计时差
Figure 606282DEST_PATH_IMAGE009
，具体计算式为：
Figure DEST_PATH_IMAGE021
；
S56、驯服单元利用时差测量处理单元实际测量值TERb1和第一预计时差
Figure 841567DEST_PATH_IMAGE003
，获得 两者之间的第一相对误差，计算公式如下：
Figure 346497DEST_PATH_IMAGE022
；
利用时差测量处理单元实际测量值TERb2和第二预计时差
Figure 769388DEST_PATH_IMAGE009
获得两者之间的第 二相对误差，计算公式如下：
Figure DEST_PATH_IMAGE023
；
S57、驯服单元利用最小二乘法对第一相对误差和第二相对误差进行曲线拟合，分 别获得第一铷原子钟的第一误差曲线斜率
Figure 969557DEST_PATH_IMAGE024
和第二铷原子钟的第二误差曲线斜率
Figure DEST_PATH_IMAGE025
；
S58、驯服单元将第一误差曲线斜率和第二误差曲线斜率分别与阈值判定，若第一误差曲线斜率和第二误差曲线斜率均大于阈值，则判定外部输入1PPS信号被注入欺骗式干扰；
S59、驯服单元判定外部输入1PPS信号出现欺骗式干扰后，给出外部输入1PPS信号出现欺骗式干扰的注入时刻，计算注入时刻后驯服单元输出的控制数据总量，并输出至DDS单元，完成回调校准并停止驯服工作；
S591、当判断第一误差曲线斜率和第二误差曲线斜率均小于或等于阈值后，判定外部输入1PPS信号欺骗式干扰消失，再次重新启动驯服作业。
需要说明的是，本实施例中，如图1所示，时频驯服模块利用双铷原子钟同时交叉进行时差数据测量和处理，其目的是提高防欺骗干扰监测识别策略可靠性。
进一步地，时频驯服模块利用时差数据TERb1 和TERb2作为监测依据，其目的是利用铷原子钟自身的频率特性来消除驯服结果对监测识别过程的影响，提高防欺骗干扰监测识别策略的准确性。
本发明提供的防欺骗式干扰铷原子钟驯服方法及系统，包括如下步骤：
S10、控制第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号，对第一信号进行分路处理，输出第一分路信号、第二分路信号和第三分路信号，对第二信号进行分路处理，输出第四分路信号、第五分路信号和第六分路信号；
S20、分别对第一分路信号和第四分路信号进行频率合成，并输出第一合成信号和第二合成信号；
S30、将第二分路信号、第五分路信号、第一合成信号和第二合成信号对应转化并输出第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
S40、以第三分路信号和第六分路信号作为频率输入基准，分别计算第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号与外部输入1PPS信号的时差，对应获得第一时差、第二时差、第三时差和第四时差；
S50、利用第一时差和第二时差监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰；对第三时差和第四时差进行驯服校准，获得控制数据，并将控制数据输入S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业，通过利用双铷原子钟的频率特性对外部输入1PPS信号进行监测识别，当外部输入1PPS信号出现欺骗式干扰时，能够有效识别出干扰，同时能够在一定程度上消除欺骗式干扰识别过程中驯服被带偏的影响，可有效避免铷原子钟的频率输出被欺骗式干扰影响，保证铷原子钟在驯服期间频率输出的准确性，具有较高的应用价值。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种防欺骗式干扰铷原子钟驯服方法，其特征在于，包括如下步骤：
S10、控制第一铷原子钟和第二铷原子钟分别输出第一信号和第二信号，对第一信号进行分路处理，输出第一分路信号、第二分路信号和第三分路信号，对第二信号进行分路处理，输出第四分路信号、第五分路信号和第六分路信号；
S20、分别对第一分路信号和第四分路信号进行频率合成，并输出第一合成信号和第二合成信号；
S30、将第二分路信号、第五分路信号、第一合成信号和第二合成信号对应转化并输出第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号；
S40、以第三分路信号和第六分路信号作为频率输入基准，分别计算第一1PPS信号、第二1PPS信号、第三1PPS信号和第四1PPS信号与外部输入1PPS信号的时差，对应获得第一时差、第二时差、第三时差和第四时差；
S50、利用第一时差和第二时差监测第一基准1PPS信号的准确性，判定外部输入1PPS信号被注入欺骗式干扰，对第三时差和第四时差进行驯服校准，获得控制数据，并将控制数据输入S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
2.根据权利要求1所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S50步骤中，还包括如下步骤：
S51、在第一运行阶段内，确认外部输入1PPS信号无欺骗式干扰；
S52、在第一运行阶段的运行时间内，利用第一时差和第二时差通过时差法获得第一铷原子钟的准确度和第二铷原子钟的准确度；利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业。
3.根据权利要求2所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S51步骤中，第一运行阶段的运行时间大于或等于2小时。
4.根据权利要求2所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S50步骤中，还包括如下步骤：
S53、在第一运行阶段结束后，进入第二运行阶段，利用第三时差和第四时差通过驯服算法获得控制数据，并将控制数据输入至S20步骤中，对第一分路信号和第四分路信号进行校准，开展对第一铷原子钟和第二铷原子钟的驯服作业；
S54、利用第一铷原子钟的准确度、第二运行阶段开始时刻的初始的第一时差、第二运行阶段累计运行时长，获得第一1PPS信号和外部输入1PPS信号的第一预计时差；
S55、利用第二铷原子钟的准确度、第二运行阶段开始时刻的初始的第二时差、第二运行阶段累计运行时长，获得第二1PPS信号和外部输入1PPS信号的第二预计时差；
S56、计算第一时差与第一预计时差的第一相对误差，计算第二时差与第二预计时差的第二相对误差；
S57、通过最小二乘法对第一相对误差和第二相对误差进行曲线拟合，分别获得第一铷原子钟的第一误差曲线斜率和第二铷原子钟的第二误差曲线斜率。
5.根据权利要求4所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S50步骤中，还包括如下步骤：
S58、将第一误差曲线斜率和第二误差曲线斜率分别与阈值判定，若第一误差曲线斜率和第二误差曲线斜率均大于阈值，则判定外部输入1PPS信号被注入欺骗式干扰；
S59、获得判定外部输入1PPS信号被注入欺骗式干扰的注入时刻，计算注入时刻以后的控制数据的输出总量，完成回调校准并停止驯服作业。
6.根据权利要求5所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S50步骤中，还包括如下步骤：
S591、当判断第一误差曲线斜率和第二误差曲线斜率均小于或等于阈值后，判定外部输入1PPS信号欺骗式干扰消失，再次启动驯服作业。
7.根据权利要求4所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S54步骤中，第一预计时差的计算式如下：
Figure 261257DEST_PATH_IMAGE002
其中，
Figure DEST_PATH_IMAGE003
为第一预计时差；
Figure 725868DEST_PATH_IMAGE004
为第二运行阶段开始时刻的初始的第一时差；
Figure DEST_PATH_IMAGE005
为 第一铷原子钟的准确度；
Figure 302343DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
8.根据权利要求7所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S55步骤中，第二预计时差的计算式如下：
Figure 776181DEST_PATH_IMAGE008
其中，
Figure DEST_PATH_IMAGE009
为第二预计时差；
Figure 387290DEST_PATH_IMAGE010
为第二运行阶段开始时刻的初始的第二时差；
Figure DEST_PATH_IMAGE011
为 第二铷原子钟的准确度；
Figure 22802DEST_PATH_IMAGE006
为第二运行阶段的累计运行时长。
9.根据权利要求8所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，所述S56步骤中，第一相对误差的计算式如下：
Figure DEST_PATH_IMAGE013
其中，
Figure 86573DEST_PATH_IMAGE014
为第一相对误差；
Figure DEST_PATH_IMAGE015
为S40步骤中实际测量获得的第一时差；
第二相对误差的计算式如下：
Figure DEST_PATH_IMAGE017
其中，
Figure 423489DEST_PATH_IMAGE018
为第二相对误差；
Figure DEST_PATH_IMAGE019
为S40步骤中实际测量获得的第二时差。
10.一种防欺骗式干扰铷原子钟驯服系统，用以执行如权利要求1至9任一项所述的防欺骗式干扰铷原子钟驯服方法，其特征在于，包括第一铷原子钟、第二铷原子钟、分路单元、DDS单元、信号转化单元、时差测量处理单元和驯服单元；
其中，所述第一铷原子钟和所述第二铷原子钟的输出端分别与所述分路单元的输入端连接，所述分路单元的第一输出端和第四输出端均与所述DDS单元的输入端连接，所述分路单元的第二输出端和第五输出端以及所述DDS单元的输出端均与所述信号转化单元的输入端连接，所述信号转化单元还用以接收外部输入1PPS信号；
所述分路单元的第三输出端和第六输出端以及所述信号转化单元的输出端均与所述时差测量处理单元的输入端连接，所述时差测量处理单元的输出端与所述驯服单元的输入端连接；
所述驯服单元的输出端与所述DDS单元的输入端连接。
CN202211025733.1A 2022-08-25 2022-08-25 防欺骗式干扰铷原子钟驯服方法及系统 Active CN115128937B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211025733.1A CN115128937B (zh) 2022-08-25 2022-08-25 防欺骗式干扰铷原子钟驯服方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211025733.1A CN115128937B (zh) 2022-08-25 2022-08-25 防欺骗式干扰铷原子钟驯服方法及系统

Publications (2)

Publication Number Publication Date
CN115128937A CN115128937A (zh) 2022-09-30
CN115128937B true CN115128937B (zh) 2022-11-04

Family

ID=83387905

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211025733.1A Active CN115128937B (zh) 2022-08-25 2022-08-25 防欺骗式干扰铷原子钟驯服方法及系统

Country Status (1)

Country Link
CN (1) CN115128937B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205375018U (zh) * 2016-02-03 2016-07-06 石家庄市经纬度科技有限公司 基于卫星授时的可驯服铷原子钟同步时钟设备
CN112147874A (zh) * 2020-11-05 2020-12-29 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
CN114884507A (zh) * 2022-03-24 2022-08-09 北京飞秒留声科技有限公司 一种铷原子钟驯服系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN205375018U (zh) * 2016-02-03 2016-07-06 石家庄市经纬度科技有限公司 基于卫星授时的可驯服铷原子钟同步时钟设备
CN112147874A (zh) * 2020-11-05 2020-12-29 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
WO2022095552A1 (zh) * 2020-11-05 2022-05-12 北京航天发射技术研究所 基于卫星授时和cpt原子钟守时的时频基准生成装置和方法
CN114884507A (zh) * 2022-03-24 2022-08-09 北京飞秒留声科技有限公司 一种铷原子钟驯服系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
原子钟与相关物理学的研究;王义遒;《物理》;20090512(第05期);全文 *
基于ARM的频率源驯服设计与实现;杨剑青等;《时间频率学报》;20141015(第04期);全文 *
频率源对雷达测速精度影响分析;杨博等;《宇航学报》;20180330(第03期);全文 *

Also Published As

Publication number Publication date
CN115128937A (zh) 2022-09-30

Similar Documents

Publication Publication Date Title
CN109299496B (zh) 一种高精度同步时钟产生方法
CN102998681B (zh) 一种卫星导航系统的高频钟差估计方法
CN110133997B (zh) 一种检测卫星时钟异常的方法
US20100141308A1 (en) Method and device for clock-data recovery
GB2462926A (en) GPS multipath reduction
CN115128937B (zh) 防欺骗式干扰铷原子钟驯服方法及系统
US7742884B2 (en) Sampling frequency control method and protective relay
EP2381266B1 (en) Self-diagnosis system and test circuit determination method
CN103675854A (zh) 基于北斗卫星的多通道授时拟合方法
US9871609B2 (en) Network system, time master station, and time slave station
CA2741758C (en) Method for correcting prediction errors of signal values with time variation subjected to interference by various uncontrollable systematic effects
WO2014161360A1 (zh) 一种同步时钟信号补偿方法、装置及计算机存储介质
CN108873023B (zh) 一种提高定位精度的观测量的处理方法
CN113992296A (zh) 时钟驯服方法、时间码监测装置及时间同步系统
US7461317B2 (en) System and method for aligning a quadrature encoder and establishing a decoder processing speed
US8482444B1 (en) Calibration of analog-to-digital converter data capture
CN109631892B (zh) 一种imu数据中断的组合导航数据处理方法
CN113114233B (zh) 星载时钟信号频率和相位跳变监测方法及系统
US8659475B2 (en) Method for detecting the distortion of a GNSS signal
US12003087B2 (en) Time synchronization between IEDs of different substations
CN113866808B (zh) 一种gnss振荡器级实时周跳探测与修复方法
KR20090024055A (ko) 산업플랜트의 트립원인 추적 방법 및 그 시스템
JP5654908B2 (ja) クロック信号生成装置およびクロック信号生成方法
KR20160134923A (ko) 아날로그-디지털 변환 장치 및 아날로그-디지털 변환 장치의 동작 방법
CN104253612B (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