CN118092120A - 一种多主钟自动切换的守时装置及方法 - Google Patents
一种多主钟自动切换的守时装置及方法 Download PDFInfo
- Publication number
- CN118092120A CN118092120A CN202410180582.XA CN202410180582A CN118092120A CN 118092120 A CN118092120 A CN 118092120A CN 202410180582 A CN202410180582 A CN 202410180582A CN 118092120 A CN118092120 A CN 118092120A
- Authority
- CN
- China
- Prior art keywords
- clock
- atomic
- frequency
- time
- signal
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims description 73
- 238000004422 calculation algorithm Methods 0.000 claims description 18
- 238000005259 measurement Methods 0.000 claims description 18
- 238000005070 sampling Methods 0.000 claims description 11
- 230000004927 fusion Effects 0.000 claims description 10
- 238000005516 engineering process Methods 0.000 claims description 9
- 230000002159 abnormal effect Effects 0.000 claims description 8
- 230000005856 abnormality Effects 0.000 claims description 6
- 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 claims description 5
- 230000015572 biosynthetic process Effects 0.000 claims description 5
- 238000003786 synthesis reaction Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000004148 unit process Methods 0.000 claims description 4
- 206010000117 Abnormal behaviour Diseases 0.000 claims description 3
- 101100492780 Arabidopsis thaliana ATI1 gene Proteins 0.000 claims description 3
- 230000007274 generation of a signal involved in cell-cell signaling Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000003321 amplification Effects 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 1
- 230000002457 bidirectional effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 229910052792 caesium Inorganic materials 0.000 description 1
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009966 trimming Methods 0.000 description 1
Landscapes
- Electric Clocks (AREA)
Abstract
本发明提供了一种多主钟自动切换的守时装置及方法，依托于无缝切换设备实现守时系统输出物理信号的参考可以在多个原子钟之间选择，扩展主钟备选数量的同时节约相位微调器的使用，同时提供了一种主钟的切换策略，减少单台原子钟受随机噪声波动对系统时间信号的影响。本发明主钟的性能不再由单一原子钟决定，能够在多个原子钟频率信号中灵活选择；主钟切换过程与以往的切换方式不同，在主钟信号经过频率控制之前进行切换，是在守时系统频率源上的切换，能够减少相位微调仪的使用数量；同时该方法提供的主钟选择策略也保障了系统时间不会因为单台钟的噪声波动而出现较大的相位偏差。
Description
技术领域
本发明涉及守时领域，尤其是一种守时方法。
背景技术
时间作为现代社会运行的基础，人类的日常生活离不开时间。随着科技的日益发展，计时的精度也日益提高。为了产生并保持一个稳定的时间尺度，当前国际上大多数的守时实验室基本采用多台原子钟来搭建原子钟组，通过钟组的内部测量与分析以削弱单台钟噪声对时间尺度的影响。
守时简单来说就是利用原子钟对时间基准进行复现，输出实时的物理信号，为物理信号提供频率参考的原子钟便称作为主钟。守时系统中通常选择原子钟组中性能最稳定、可预测性最好的原子钟作为主钟，并使用相位微调器对主钟的频率和相位偏差进行调整。但无论如何，都会由于原子钟本身的问题而出现信号中断或性能变差的现象，无法保证输出物理信号的连续稳定。为解决这一问题，目前大多数守时实验室通常选择两台主钟与两台相位微调器构成输出主备两路信号，保持两路信号的相位与频率相一致，在主路信号发生故障或恢复时，通过切换开关手动控制实现系统时间信号在主备两路信号之间的快速切换，以保证系统时间的可靠性。但对于一些小型化守时系统来说，参与守时的原子钟数目较少，同时为节约成本，可能无法搭建两套相对独立的物理信号输出手段。另一方面，若参与守时的钟组中，各原子钟的性能相差不大，主钟的数量可以进一步扩展，以提高本地物理信号的可靠性。但目前的主备主钟切换都是对原子钟经过相位微调器驾驭后的信号进行切换，是在守时系统信号的输出端进行切换，因而每有一台主钟便需要一台相位微调仪，加大了对相位微调仪的需求，因此需要对现有的这种主钟频率控制方法进行改进。
发明内容
为了克服现有技术的不足，本发明提供一种多主钟自动切换的守时装置及方法，依照原子钟的频率驾驭原理，依托于无缝切换设备实现守时系统输出物理信号的参考可以在多个原子钟之间选择，扩展主钟备选数量的同时节约相位微调器的使用，同时提供了一种主钟的切换策略，减少单台原子钟受随机噪声波动对系统时间信号的影响。
本发明解决其技术问题所采用的技术方案是：
一种多主钟自动切换的守时装置，包括原子钟钟组单元、测量比对单元、主钟信号产生单元和频率控制单元；
所述的原子钟钟组单元提供多路主钟频率信号；
所述的测量比对单元通过计数器和比相仪完成各原子钟信号、远程溯源信号和守时系统时间信号之间的钟差和相差数据测量；
所述的主钟信号产生单元采用无缝切换设备，以输入的多个主钟频率信号为参考，实现其输出的信号在多个频率源参考之间的无缝切换，同时基于数字频率合成技术，实现输出信号的相位与频率调整；
所述的频率控制单元对测量比对单元所获的比对数据进行处理和分析，首先利用本地的比对数据基于原子时算法得到一个本地时间基准；之后利用远程溯源数据对本地时间基准进行修正；最后通过预报的方式计算每台主钟相对于修正后的本地时间基准的频率控制量，利用相位微调仪对当前无缝切换设备输出的信号进行频率控制，产生稳定的系统时间信号，并在主钟切换后更新频率控制量，保障切换后系统时间不发生频率跳变。
所述的原子钟钟组单元中，原子钟的数目大于三台，对各原子钟两两比对进行主钟的噪声波动判断。
所述的无损切换设备包括数据采样单元、频率融合单元和DDS频率模块；每一路主钟频率信号分别输入一个数据采样单元，采用下变频技术进行鉴相，测量原子钟相对本地时钟的频率差别与相位差别；数据采样单元的输出信号输入频率融合单元，以本地时钟为基准，通过频率融合单元与DDS频率模块，对所选择的原子钟10MHz频率进行输出；其输出信号的频率控制方式与相位微调器相同，在进行主钟切换时对输出的信号直接进行频率与相位调整。
所述的原子时算法采用Kalman算法、ALOGS算法或ATI1算法。
本发明还提供一种利用多主钟自动切换的守时装置的守时方法，包括以下步骤：
1)将相位微调仪输出的本地10MHz信号作为多通道时间间隔计数器的频率参考，输出的1PPS信号作为开门信号，将原子钟的1PPS信号和外部溯源的1PPS信号分别作为关门信号，由工控机中的上位机软件实现比对数据的采集，记为MC-clocki和MC-UTCk，MC指相位微调仪输出的1PPS秒信号，clocki和UTCk分别表示各原子钟i和远程溯源1PPS秒信号；外部溯源的1PPS信号由UTC(k)直接给出物理信号，若没有直连信号时，需通过卫星共视、双向卫星时间传递或其他方法获取与外部溯源的比对数据，并存储至工控机中；
获取原有比对数据后，为减少测量中的偶然误差影响，根据3-Sigma原则对MC-clocki数据中的粗差进行剔除，同时采用线性插值对剔除的部分进行插值；
2)综合原子时的计算；
综合原子时的计算参考BIPM在计算UTC时使用的ALGOS算法，基本核心思想均为以可预测性对参与计算的原子钟进行赋权，进而对原子钟的比对数据进行加权平均计算；计算流程图如图3所示，步骤如下：
2.1)首先依照设定的综合原子时计算频度与所使用的历史数据长度，从本地文件中读取比对数据；
综合原子时的计算是为了保障溯源中断后的本地时间保持能力，计算频度的设定与使用的溯源链路有关，若溯源采用直连、共视或双向时间比对，溯源的更新频率较快，原子时的计算频度便可以以较快的频率进行计算，若采用BIPM公报的UTC或UTC(r)来进行溯源，原子时的计算频度不宜太快；采用直连物理信号进行溯源，故设定综合原子时的计算频度为1小时；使用历史数据长度过长会导致原子时的计算无法准确反应原子钟近期的波动情况，过短则会导致原子时的稳定性下降，本发明选择36条历史数据；
2.2)原子时计算；
原子钟的权重取决于原子钟的可预测性，也就是原子钟对一个模型的符合程度，本发明采用一次多项式作为计算模型：
y＝a+b*t
在获取36条MC-clocki的小时比对数据后，对每个数据扣除该原子钟的在该时刻的模型值得到残差，初始计算时a和b为0，并上次计算获得的原子钟权重对残差进行加权平均得到原子时MC-TA，并将最后一条MC-TA保存在本地文件中，作为本次计算的原子时，初始计算时每台钟的权重为1/N，N为参与计算的原子钟数目，取N为3；
由MC-clocki与MC-TA得到每台钟的36条TA-clocki，采用最小二乘法更新每台钟的a和b值，并从TA-clocki减去新的模型值，依照得到残差的平方和大小计算每台钟的权重，权重计算公式如下：
权重计算公式需要进行迭代，以保障计算的稳定性，本发明迭代3次，同时原子钟的权重设定一个上限，以防止单台原子时完全却决于每台原子钟的运行，影响原子时的可靠性，本发明权重上限为2/N，对于超出上限的原子钟，超出的权重部分需要按照比例重新分配至其他原子钟，同时某台原子钟的权重过低时，说明原子钟的可预测性极差，需要在本次计算中设置权重为0，本发明权重的下限为1％；
2.3)原子时频率校准；
步骤2.2)中计算的MC-TA仅由本地原子钟比对数据进行加权计算，原子时需要与外部的溯源保持相对一致，需要对原子时进行相位和频率校准；
由MC-UTCk与MC-TA得到TA-UTCk，数据的条数选择24条，采用最小二乘法对24条TA-UTCk进行一次多项式建模y＝a+b*t，并从MC-TA中扣除拟合值即为较频后的MC-TA，每新产生的MC-TA均需要扣除拟合值y，模型中的a,b值的更新频度与外部溯源的更新频度和原子时计算频度有关，本发明中为1小时更新一次；
每次完成原子时相位和频率校准后，由MC-TA与MC-clocki得到TA-clocki，存储至本地文件中；
3)主钟频率驾驭量计算；
主钟的频率驾驭量与最新的MC-TA相位差值、当前主钟相对于TA的频率偏差和驾驭频度有关，主钟的驾驭频度应大于外部溯源和原子时的更新频度，本发明的主钟驾驭频度为1小时，与原子时的计算频度一致，确定好当前主钟后，驾驭量offset的计算公式为：
TA-clockm为最近获得的原子时与主钟的相位差值，k为调整系数，为经验常数，本发明取2，bm为主钟相对于原子时的频率偏差，对最近的24条TA-clockm采用最小二乘法进行一次多项式建模获得，数据条数可根据实际的需要进行调整；
4)主钟的切换控制；
主钟的切换判断需要根据当前主钟的性能状态来进行，当主钟发生相位跳频或频率跳变等异常行为时，需要及时对主钟信号进行切换，保障稳定的本地物理信号输出；本发明采用动态allan方差法来完成原子钟的异常判断，实现步骤如下：
4.1)对于每秒获取的MC-clocki，由两两比对获取clocki-clockj计算，与100s前的clocki-clockj之差，统计最近的50条差值，计算clocki-clockj的100s稳定度，当前时刻计算的100s稳定度超出使用原子钟的标称100s稳定度的2倍时，认为i和j中有一台原子钟发生异常，3组比对值判断完成后便可确定发生异常的原子钟；稳定度的平均时间和判断的放大倍数根据实际的原子钟性能进行确定；
4.2)主钟发生异常后，将主钟切换至未发生异常的原子钟即可，本发明中原子钟仅有3台，若同时有两台原子钟发生异常，则无法判断具体的原子钟，此时则不进行主钟切换；
4.3)主钟切换时需及时依照步骤3)中的驾驭量计算公式更正相位微调器的频率驾驭量。
本发明的有益效果在于主钟的性能不再由单一原子钟决定，能够在多个原子钟频率信号中灵活选择；主钟切换过程与以往的切换方式不同，在主钟信号经过频率控制之前进行切换，是在守时系统频率源上的切换，能够减少相位微调仪的使用数量；同时该方法提供的主钟选择策略也保障了系统时间不会因为单台钟的噪声波动而出现较大的相位偏差。
附图说明
图1是无损切换设备原理框图；
图2是多主钟守时方法设备连接实施例示意图；
图3是综合原子时计算流程图。
具体实施方式
下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。
本发明首先需要解决的是不同原子钟之间频率源参考的切换问题，该问题的解决方案如图1所示，依托于无损切换设备完成主钟信号在多个原子钟信号之间切换的功能。无损切换设备对多路外原子钟频率信号输入采用下变频技术进行高精度鉴相，可实时快速的测量多路原子钟相对本地时钟的频率差别与相位差别，进一步以本地为时钟基准，通过频率融合与DDS控制，对所选择的原子钟10MHz频率进行输出。其输出信号的频率控制方式与相位微调器相同，在使用该种方式进行主钟切换时也可以对输出的信号直接进行频率与相位调整，调整原理与常见的相位微调器相一致，从而起到替代相位微调器的作用，节省硬件成本。
本发明需要解决的第二个问题是当前主钟的驾驭量计算问题，由于主钟不再由单一的原子钟所确定，所以在主钟驾驭过程中为保持系统时间与时间基准的一致性，就需要对每一个参与主钟无损切换的原子钟速率进行预报，并根据当前所选择的主钟来实时更新对应的频率控制量，以保障输出信号与参考的频率始终一致。
本发明需要解决的第三个问题是何时进行主钟的切换。在一般的守时系统中，主钟的频率驾驭分为两个部分，一个是为了修正主钟相对于参考时间尺度的相对频率偏差，另一部分是为了修正由于主钟噪声波动所引起的相位偏差，保障系统时间与参考信号的绝对同步，因而主钟的选择上应当选取噪声波动小的原子钟来作为主钟，以避免在主钟的驾驭频度内出现系统时间与参考信号出现较大的相位偏差，因而需要实时监测当前主钟的噪声波动。通常，原子钟的噪声都是通过频率稳定度来表征，也就是通过allan方差来进行统计估计。在原子钟出厂时，不同型号的原子钟均会提供相应的频率稳定度标称值，在进行原子钟的噪声监测时可以依照该标称值或实际应用的频率稳定度指标需求来判断原子钟是否产生了较大的噪声波动，从而判断是否切换当前的主钟。
本发明首先保障守时系统能够输出一个稳定的系统时间，其次要在信号稳定的前提下保障多主钟的切换，就需要能够在切换时准确计算每台原子钟相对于参考信号的频率偏差，实时的更正主钟驾驭量，同时需要监测当前主钟的相位波动，当出现较大的相位波动时，对当前主钟进行切换。为达到这三个目的，本发明需包括原子钟钟组单元、测量比对单元、主钟无缝切换单元和频率控制单元。
原子钟钟组单元为守时系统提供多路稳定的主钟信号，保障守时系统的可靠性。主钟的噪声波动需要在原子钟的两两比对数据进行判断，所以对于守时系统来说，原子钟的数目最好大于三台。
测量比对单元通过计数器和比相仪完成各原子钟信号、远程溯源信号和守时系统时间信号之间的钟差和相差数据测量，为后续原子钟的主钟切换判断、原子时计算和主钟频率偏差控制量计算等提供数据保障。
主钟信号产生单元依托于现有的频率源无缝切换设备，以输入的多个主钟频率信号为参考，对压控振荡器进行锁相控制来实现其输出的信号在多个频率源参考之间的无缝切换，同时基于数字频率合成技术，实现输出信号的相位与频率调整。
频率控制单元对数据比对单元所获的比对数据进行处理和分析，首先利用本地的比对数据基于原子时算法(Kalman算法、ALOGS算法和ATI1算法等)得到一个本地时间基准，之后利用远程溯源数据对本地时间基准进行修正，计算本地时间基准的目的是在外部溯源链路中断后，提供自主守时能力。最后通过预报的方式计算每台主钟相对于修正后的本地时间基准的频率控制量，利用相位微调仪对当前无缝切换设备输出的信号进行频率控制，产生稳定的系统时间信号，并在主钟切换后及时更新频率控制量，保障切换后系统时间不发生频率跳变。
本实施例的硬件部分包括三台铯原子钟、无缝切换设备、相位微调器、多通道时间间隔计数器、多通道比相仪和工控机，在设计过程中原子钟的类型和数目，采集设备的个数需根据实际应用需求进行调整。具体设备连接方式如图2所示。
一种多主钟自动切换的守时装置，包括原子钟钟组单元、测量比对单元、主钟信号产生单元和频率控制单元；
所述的原子钟钟组单元提供多路主钟频率信号；
所述的测量比对单元通过计数器和比相仪完成各原子钟信号、远程溯源信号和守时系统时间信号之间的钟差和相差数据测量；
所述的主钟信号产生单元采用无缝切换设备，以输入的多个主钟频率信号为参考，实现其输出的信号在多个频率源参考之间的无缝切换，同时基于数字频率合成技术，实现输出信号的相位与频率调整；
所述的频率控制单元对测量比对单元所获的比对数据进行处理和分析，首先利用本地的比对数据基于原子时算法得到一个本地时间基准；之后利用远程溯源数据对本地时间基准进行修正；最后通过预报的方式计算每台主钟相对于修正后的本地时间基准的频率控制量，利用相位微调仪对当前无缝切换设备输出的信号进行频率控制，产生稳定的系统时间信号，并在主钟切换后更新频率控制量，保障切换后系统时间不发生频率跳变。
所述的原子钟钟组单元中，原子钟的数目大于三台，对各原子钟两两比对进行主钟的噪声波动判断。
所述的无损切换设备包括数据采样单元、频率融合单元和DDS频率模块；每一路主钟频率信号分别输入一个数据采样单元，采用下变频技术进行鉴相，测量原子钟相对本地时钟的频率差别与相位差别；数据采样单元的输出信号输入频率融合单元，以本地时钟为基准，通过频率融合单元与DDS频率模块，对所选择的原子钟10MHz频率进行输出；其输出信号的频率控制方式与相位微调器相同，在进行主钟切换时对输出的信号直接进行频率与相位调整。
本实施例的频率控制部分需要完成时差数据的采样和预处理，综合原子时的计算，频率驾驭量计算，主钟切换控制四个过程。
1)时差数据的采样和预处理
时差数据的采集通过工控机和多通道时间间隔计数器来实现，将相位微调器输出的本地10MHz信号作为多通道时间间隔计数器的频率参考，输出的1PPS信号作为开门信号，原子钟的1PPS信号和外部溯源的1PPS信号分别作为关门信号，由工控机中的上位机软件实现比对数据的采集。记为MC-clocki和MC-UTCk，MC指相位微调仪输出的1PPS秒信号，clocki和UTCk为各原子钟和远程溯源1PPS秒信号。本实施例的溯源1PPS信号由UTC(k)直接给出物理信号，若没有直连信号时，需通过卫星共视、双向卫星时间传递或其他方法获取与外部溯源的比对数据，并存储至工控机中。
获取原有比对数据后，为减少测量中的偶然误差影响，需要根据3-Sigma原则对MC-clocki数据中的粗差进行剔除，同时采用线性插值对剔除的部分进行插值。
2)综合原子时的计算
本实施例中综合原子时的计算参考BIPM在计算UTC时使用的ALGOS算法，基本核心思想均为以可预测性对参与计算的原子钟进行赋权，进而对原子钟的比对数据进行加权平均计算。计算流程图如图3所示，步骤如下：
2.1)首先依照设定的综合原子时计算频度与所使用的历史数据长度，从本地文件中读取比对数据。
综合原子时的计算是为了保障溯源中断后的本地时间保持能力，计算频度的设定与使用的溯源链路有关，若溯源采用直连、共视或双向时间比对，溯源的更新频率较快，原子时的计算频度便可以以较快的频率进行计算，若采用BIPM公报的UTC或UTC(r)来进行溯源，原子时的计算频度不宜太快。本实施例采用的直连物理信号进行溯源，故设定综合原子时的计算频度为1小时。使用历史数据长度过长会导致原子时的计算无法准确反应原子钟近期的波动情况，过短则会导致原子时的稳定性下降，本实施例选择36条历史数据。
2.2)原子时计算
原子钟的权重取决于原子钟的可预测性，也就是原子钟对一个模型的符合程度，在本示例中采用一次多项式作为计算模型：
y＝a+b*t
在获取36条MC-clocki的小时比对数据后，对每个数据扣除该原子钟的在该时刻的模型值得到残差，初始计算时a和b为0，并上次计算获得的原子钟权重对残差进行加权平均得到原子时MC-TA，并将最后一条MC-TA保存在本地文件中，作为本次计算的原子时，初始计算时每台钟的权重为1/N，N为参与计算的原子钟数目，本实施例中N为3。
由MC-clocki与MC-TA得到每台钟的36条TA-clocki。采用最小二乘法更新每台钟的a和b值，并从TA-clocki减去新的模型值，依照得到残差的平方和大小计算每台钟的权重，权重计算公式如下：
上述计算通常需要进行迭代，以保障计算的稳定性，本实施例迭代3次。同时原子钟的权重需要设定一个上限，以防止单台原子时完全却决于每台原子钟的运行，影响原子时的可靠性，本实施例采用权重上限为2/N，对于超出上限的原子钟，超出的权重部分需要按照比例重新分配至其他原子钟，同时某台原子钟的权重过低时，说明原子钟的可预测性极差，需要在本次计算中设置权重为0，本实施例权重的下限为1％。
2.3)原子时频率校准；
步骤2.2)中计算的MC-TA仅由本地原子钟比对数据进行加权计算，原子时需要与外部的溯源保持相对一致，需要对原子时进行相位和频率校准。
由MC-UTCk与MC-TA得到TA-UTCk，数据的条数本实施例中选择24条，采用最小二乘法对24条TA-UTCk进行一次多项式建模y＝a+b*t，并从MC-TA中扣除拟合值即为较频后的MC-TA。每新产生的MC-TA均需要扣除拟合值y，模型中的a,b值的更新频度与外部溯源的更新频度和原子时计算频度有关，本实施例中为1小时更新一次。
每次完成原子时相位和频率校准后，由MC-TA与MC-clocki得到TA-clocki，存储至本地文件中。
3)主钟频率驾驭量计算；
主钟的频率驾驭量与最新的MC-TA相位差值、当前主钟相对于TA的频率偏差和驾驭频度有关。主钟的驾驭频度应大于外部溯源和原子时的更新频度，本实施例的主钟驾驭频度为1小时,与原子时的计算频度一致。确定好当前主钟后，驾驭量offset的计算公式为：
TA-clockm为最近获得的原子时与主钟的相位差值，k为调整系数，为经验常数，本实施例为2，bm为主钟相对于原子时的频率偏差，对最近的24条TA-clockm采用最小二乘法进行一次多项式建模获得，数据条数可根据实际的需要进行调整。
4)主钟的切换控制；
主钟的切换判断需要根据当前主钟的性能状态来进行，当主钟发生相位跳频或频率跳变等异常行为时，需要及时对主钟信号进行切换，保障稳定的本地物理信号输出。本实施例采用动态allan方差法来完成原子钟的异常判断，实现步骤如下：
4.1)对于每秒获取的MC-clocki，由两两比对获取clocki-clockj计算，其与100s前的clocki-clockj之差，统计最近的50条差值，计算clocki-clockj的100s稳定度，当前时刻计算的100s稳定度超出使用原子钟的标称100s稳定度的2倍时，认为i和j中有一台原子钟发生异常，3组比对值判断完成后便可确定发生异常的原子钟。稳定度的平均时间和判断的放大倍数可根据实际的原子钟性能进行确定。
4.2)主钟发生异常后，将主钟切换至未发生异常的原子钟即可，本实施例中原子钟仅有3台，若同时有两台原子钟发生异常，则无法判断具体的原子钟，此时则不进行主钟切换。
4.3)主钟切换时需及时依照3)中的驾驭量计算公式更正相位微调器的频率驾驭量。
Claims (9)
1.一种多主钟自动切换的守时装置，包括原子钟钟组单元、测量比对单元、主钟信号产生单元和频率控制单元；其特征在于：
所述的原子钟钟组单元提供多路主钟频率信号；
所述的测量比对单元通过计数器和比相仪完成各原子钟信号、远程溯源信号和守时系统时间信号之间的钟差和相差数据测量；
所述的主钟信号产生单元采用无缝切换设备，以输入的多个主钟频率信号为参考，实现其输出的信号在多个频率源参考之间的无缝切换，同时基于数字频率合成技术，实现输出信号的相位与频率调整；
所述的频率控制单元对测量比对单元所获的比对数据进行处理和分析，首先利用本地的比对数据基于原子时算法得到一个本地时间基准；之后利用远程溯源数据对本地时间基准进行修正；最后通过预报的方式计算每台主钟相对于修正后的本地时间基准的频率控制量，利用相位微调仪对当前无缝切换设备输出的信号进行频率控制，产生稳定的系统时间信号，并在主钟切换后更新频率控制量，保障切换后系统时间不发生频率跳变。
2.根据权利要求1所述的多主钟自动切换的守时装置，其特征在于：
所述的原子钟钟组单元中，原子钟的数目大于三台，对各原子钟两两比对进行主钟的噪声波动判断。
3.根据权利要求1所述的多主钟自动切换的守时装置，其特征在于：
所述的无损切换设备包括数据采样单元、频率融合单元和DDS频率模块；每一路主钟频率信号分别输入一个数据采样单元，采用下变频技术进行鉴相，测量原子钟相对本地时钟的频率差别与相位差别；数据采样单元的输出信号输入频率融合单元，以本地时钟为基准，通过频率融合单元与DDS频率模块，对所选择的原子钟10MHz频率进行输出；其输出信号的频率控制方式与相位微调器相同，在进行主钟切换时对输出的信号直接进行频率与相位调整。
4.根据权利要求1所述的多主钟自动切换的守时装置，其特征在于：
所述的原子时算法采用Kalman算法、ALOGS算法或ATI1算法。
5.一种利用权利要求1所述的多主钟自动切换的守时装置的守时方法，其特征在于包括下述步骤：
1)将相位微调仪输出的本地10MHz信号作为多通道时间间隔计数器的频率参考，输出的1PPS信号作为开门信号，将原子钟的1PPS信号和外部溯源的1PPS信号分别作为关门信号，由工控机中的上位机软件实现比对数据的采集，记为MC-clocki和MC-UTCk，MC指相位微调仪输出的1PPS秒信号，clocki和UTCk分别表示各原子钟i和远程溯源1PPS秒信号；外部溯源的1PPS信号由UTC(k)直接给出物理信号，若没有直连信号时，需通过卫星共视、双向卫星时间传递或其他方法获取与外部溯源的比对数据，并存储至工控机中；
获取原有比对数据后，为减少测量中的偶然误差影响，根据3-Sigma原则对MC-clocki数据中的粗差进行剔除，同时采用线性插值对剔除的部分进行插值；
2)综合原子时的计算；
综合原子时的计算参考BIPM在计算UTC时使用的ALGOS算法，基本核心思想均为以可预测性对参与计算的原子钟进行赋权，进而对原子钟的比对数据进行加权平均计算；
3)主钟频率驾驭量计算；
主钟的频率驾驭量与最新的MC-TA相位差值、当前主钟相对于TA的频率偏差和驾驭频度有关，主钟的驾驭频度应大于外部溯源和原子时的更新频度，主钟驾驭频度为1小时，与原子时的计算频度一致，确定好当前主钟后，驾驭量offset的计算公式为：
TA-clockm为最近获得的原子时与主钟的相位差值，k为调整系数，bm为主钟相对于原子时的频率偏差，对最近的24条TA-clockm采用最小二乘法进行一次多项式建模获得，数据条数可根据实际的需要进行调整；
4)主钟的切换控制；
主钟的切换判断需要根据当前主钟的性能状态进行，当主钟发生相位跳频或频率跳变等异常行为时，需要及时对主钟信号进行切换，保障稳定的本地物理信号输出。
6.根据权利要求5所述的多主钟自动切换的守时装置的守时方法，其特征在于：所述步骤2)的具体步骤如下：
2.1)首先依照设定的综合原子时计算频度与所使用的历史数据长度，从本地文件中读取比对数据；
采用直连物理信号进行溯源，故设定综合原子时的计算频度为1小时；选择36条历史数据；
2.2)原子时计算；
原子钟的权重取决于原子钟的可预测性，也就是原子钟对一个模型的符合程度，采用一次多项式作为计算模型：
y＝a+b*t
在获取36条MC-clocki的小时比对数据后，对每个数据扣除该原子钟的在该时刻的模型值得到残差，初始计算时a和b为0，并上次计算获得的原子钟权重对残差进行加权平均得到原子时MC-TA，并将最后一条MC-TA保存在本地文件中，作为本次计算的原子时，初始计算时每台钟的权重为1/N，N为参与计算的原子钟数目，取N为3；
由MC-clocki与MC-TA得到每台钟的36条TA-clocki，采用最小二乘法更新每台钟的a和b值，并从TA-clocki减去新的模型值，依照得到残差的平方和大小计算每台钟的权重，权重计算公式如下：
权重计算公式需要进行迭代，以保障计算的稳定性，迭代3次，同时原子钟的权重设定一个上限，以防止单台原子时完全却决于每台原子钟的运行，影响原子时的可靠性，对于超出上限的原子钟，超出的权重部分需要按照比例重新分配至其他原子钟，同时某台原子钟的权重过低时，说明原子钟的可预测性极差，需要在本次计算中设置权重为0；
2.3)原子时频率校准；
步骤2.2)中计算的MC-TA仅由本地原子钟比对数据进行加权计算，原子时需要与外部的溯源保持相对一致，需要对原子时进行相位和频率校准；
由MC-UTCk与MC-TA得到TA-UTCk，数据的条数选择24条，采用最小二乘法对24条TA-UTCk进行一次多项式建模y＝a+b*t，并从MC-TA中扣除拟合值即为较频后的MC-TA，每新产生的MC-TA均需要扣除拟合值y，模型中的a,b值的更新频度与外部溯源的更新频度和原子时计算频度有关，1小时更新一次；
每次完成原子时相位和频率校准后，由MC-TA与MC-clocki得到TA-clocki，存储至本地文件中。
7.根据权利要求5所述的多主钟自动切换的守时装置的守时方法，其特征在于：所述步骤4)中采用动态allan方差法来完成原子钟的异常判断，实现步骤如下：
4.1)对于每秒获取的MC-clocki，由两两比对获取clocki-clockj计算，与100s前的clocki-clockj之差，统计最近的50条差值，计算clocki-clockj的100s稳定度，当前时刻计算的100s稳定度超出使用原子钟的标称100s稳定度的2倍时，认为i和j中有一台原子钟发生异常，3组比对值判断完成后便可确定发生异常的原子钟；稳定度的平均时间和判断的放大倍数根据实际的原子钟性能进行确定；
4.2)主钟发生异常后，将主钟切换至未发生异常的原子钟即可，原子钟仅有3台，若同时有两台原子钟发生异常，则无法判断具体的原子钟，此时则不进行主钟切换；
4.3)主钟切换时需及时依照步骤3)中的驾驭量计算公式更正相位微调器的频率驾驭量。
8.根据权利要求5所述的多主钟自动切换的守时装置的守时方法，其特征在于：所述的调整系数k取2。
9.根据权利要求5所述的多主钟自动切换的守时装置的守时方法，其特征在于：权重上限为2/N，权重的下限为1％。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410180582.XA CN118092120A (zh) | 2024-02-18 | 2024-02-18 | 一种多主钟自动切换的守时装置及方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410180582.XA CN118092120A (zh) | 2024-02-18 | 2024-02-18 | 一种多主钟自动切换的守时装置及方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118092120A true CN118092120A (zh) | 2024-05-28 |
Family
ID=91146923
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410180582.XA Pending CN118092120A (zh) | 2024-02-18 | 2024-02-18 | 一种多主钟自动切换的守时装置及方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118092120A (zh) |
-
2024
- 2024-02-18 CN CN202410180582.XA patent/CN118092120A/zh active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2553495C (en) | Apparatus and method for compensating the drift of a local clock used as sampling frequency | |
| CN109581856B (zh) | 一种基于高性能晶振频率校准的对时守时方法 | |
| US8954294B2 (en) | Measuring device and measuring method for continuous physical quantity | |
| CN105676627A (zh) | 守时系统主备主钟无缝切换系统与方法 | |
| CN112636857A (zh) | 一种基于残差修正的高精度时频信号动态驾驭方法 | |
| JP2007078405A (ja) | ソフトウェア時計の計時プログラム | |
| CN101090311B (zh) | Gps定时脉冲失效后产生维持时钟输出的方法和系统 | |
| CN104297716A (zh) | 一种智能电表rtc校准方法 | |
| CN109100927B (zh) | 无溯源参考条件下时间统一系统时间自主同步方法 | |
| JP5972450B2 (ja) | 信号処理装置 | |
| US20110309864A1 (en) | Apparatus and Method for Generating a Timing Signal | |
| CN116382058A (zh) | 一种基于北斗的卫星共视授时系统及方法 | |
| CN118092120A (zh) | 一种多主钟自动切换的守时装置及方法 | |
| CN105026938A (zh) | 信号处理装置 | |
| CN115733489B (zh) | 一种芯片原子钟驯服方法、装置和系统 | |
| CN110611544A (zh) | 一种利用rs232串口精确对时的方法及装置 | |
| CN102096440A (zh) | 一种保持实时时钟电路高精度的方法 | |
| JP4788260B2 (ja) | 多変量伝送器 | |
| CN115694703A (zh) | 一种时钟同步方法及装置 | |
| CN110618604B (zh) | 一种利用ntp辅助作源提高守时精度的方法及装置 | |
| CN210351192U (zh) | 一种利用rs232串口精确对时的装置 | |
| CN110824237B (zh) | 一种脉冲采样时间间隔调整方法及装置 | |
| CN116980065B (zh) | 时钟校准方法、装置、终端设备及存储介质 | |
| CN116909351B (zh) | 一种时钟芯片内部时钟精度校正方法 | |
| CN114978365B (zh) | 一种基于中心基准的多节点时频信号巡检及数据处理方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication |