Analysis method and system for the integrity of multiple reference source inputs in time and frequency synchronization network

Abstract

The present invention provides a method and system for analyzing the integrity of multiple reference source inputs of a time and frequency synchronization network, which relates to the technical field of time and frequency synchronization and is used for evaluating and selecting the best clock reference source among at least two clock reference sources in a time and frequency synchronization network. The method comprises obtaining a data set of at least two clock reference sources to be compared, wherein each data set contains a plurality of attribute values ​​of the clock reference sources to which it belongs; comparing the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order, and performing long-term measurement on the phase center deviations of multiple clock reference sources by utilizing the good long-term stability of a local rubidium atomic clock; calculating a weighted clock quality evaluation value by comparing the attribute values ​​and phase center deviations of the clock reference sources, and selecting the best clock reference source for taming a local rubidium atomic clock; and generating and outputting a high-precision time signal and a frequency signal based on the tamed local rubidium atomic clock.

Description

Technical Field

[0001] The present invention relates to the technical field of time-frequency synchronization, and in particular to a method and system for analyzing the integrity of multiple reference source inputs in a time-frequency synchronization network. Background Art

[0002] my country's time and frequency synchronization network is divided into inter-provincial backbone transmission network layer, intra-provincial backbone transmission network layer and local transmission network layer. There are various high-precision time synchronization transmission methods, mainly including Beidou / GPS satellites, time and frequency synchronization network, etc. Different timing methods, due to factors such as transmission paths and timing systems, not only have fixed time deviations, but also different time stability. At present, most timing equipment has multi-source reception function, which can receive timing signals from Beidou satellites, GPS satellites, time and frequency synchronization networks, etc., to tame local clocks and calibrate local time.

[0003] The standard BMC algorithm published in the YD / T2375-2019 Technical Requirements for High-Precision Time Synchronization standard includes a data set comparison algorithm and a state decision algorithm. Among them, the data set comparison algorithm is widely used in the reference source selection of the time and frequency synchronization network.

[0004] The YD / T2375-2019 Technical Requirements for High-Precision Time Synchronization standard is applicable to high-precision time synchronization networks that use PTP technology and are mainly based on ground transmission. This standard publishes an optimal clock algorithm (BMC algorithm). The data set comparison algorithm in the standard BMC algorithm is currently widely used in time and frequency synchronization network equipment. The algorithm compares the quality of the data group and the information contained in the synchronization message, and calculates the recommended state of each port based on the data group comparison results to select the best clock.

[0005] However, the existing data set comparison algorithm only compares the values ​​of GM priority 1, GM priority II, GM level, GM accuracy, GM deviation change, GM identifier, and sender ID in the timing message. The long-term timing quality of the superior source cannot be accurately evaluated using the multiple attribute values. In particular, the impact of slow changes in the accuracy and deviation values ​​of the superior source on the timing system cannot be accurately evaluated. This problem will not only affect the timing system's choice of reference source, but will also affect the timing accuracy of the timing system.

[0006] Therefore, it is necessary to provide a method and system for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network to solve the above technical problems. Summary of the invention

[0007] In order to solve the above technical problems, the present invention provides a method and system for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network. By adding a weighted clock quality assessment method to the original algorithm, the long-term clock accuracy of the reference source is assessed. At the same time, it can accurately reflect the long-term clock quality of the reference source, ensure the reliability and accuracy of the reference source selection, and significantly improve the local clock accuracy level.

[0008] The present invention provides a method for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network, which is used to evaluate and select the optimal clock reference source among at least two clock reference sources in the time and frequency synchronization network. The method comprises the following steps:

[0009] In a time and frequency synchronization network, data sets of at least two clock reference sources to be compared are obtained, each data set including multiple attribute values ​​of the clock reference source;

[0010] Compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order;

[0011] The phase center deviation of multiple clock reference sources is measured using a local rubidium atomic clock.

[0012] A weighted clock quality assessment value is obtained based on the comparison results of the attribute values ​​of the clock reference source and the phase center deviation calculation;

[0013] Screening out an optimal clock reference source based on the weighted clock quality evaluation value, and using the optimal clock reference source to tame a local rubidium atomic clock;

[0014] According to the disciplined local rubidium atomic clock, the time signal and frequency signal are generated and output.

[0015] Preferably, the attribute values ​​include GM identifier, GM priority I, GM level, GM accuracy, GM deviation change value, GM priority II and sender ID.

[0016] Preferably, when comparing attribute values, first compare whether the GM identifiers of any two clock reference sources are the same. If they are the same, enter the clock reference source hop count comparison processing flow; if they are not the same, compare the remaining attribute values ​​in sequence according to a predefined comparison order.

[0017] Preferably, when the GM identifiers are different, the predefined comparison order is:

[0018] a. Compare whether the values ​​of GM priority I of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority I, and take the clock reference source with a smaller GM priority I value as the optimal clock reference source, and output the result;

[0019] b. Compare whether the GM levels of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of the GM levels, and use the clock reference source with a smaller GM level as the optimal clock reference source, and output the result;

[0020] c. Compare whether the GM precision of any two clock reference sources is the same. If they are the same, proceed to the next step. If they are not the same, compare the numerical values ​​of the GM precision, and use the clock reference source with the smaller GM precision value as the optimal clock reference source, and output the result;

[0021] d. Compare the GM deviation change values ​​of any two clock reference sources to see if they are the same. If they are the same, proceed to the next step. If they are not the same, compare the GM deviation change values, and use the clock reference source with the smaller GM deviation change value as the optimal clock reference source, and output the result;

[0022] e. Compare whether the phase center offsets of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, calculate the weighted clock quality assessment value based on the attribute value of the clock reference source and the phase center offset, and output the result;

[0023] f. Compare whether the values ​​of GM priority II of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority II, and use the clock reference source with a smaller GM priority II value as the optimal clock reference source, and output the result;

[0024] g. Compare the values ​​of the transmitter IDs of any two clock reference sources, and use the clock reference source with the smaller transmitter ID value as the optimal clock reference source, and output the result.

[0025] Preferably, when the phase center offsets of any two clock reference sources are different, the method further includes the following steps:

[0026] Dynamically assign weighted values ​​of the arbitrary two clock reference sources according to the magnitude of the phase center offset of the clock reference source, and calculate a weighted clock quality assessment value according to the attribute values ​​and the phase center offset of the arbitrary two clock reference sources;

[0027] If the weighted clock quality evaluation values ​​are the same, the values ​​of GM priority II and sender ID are compared in turn, and the clock reference source with a smaller attribute value is used as the optimal clock reference source. If they are not the same, the weighted clock quality evaluation values ​​are compared, and the clock reference source with a smaller weighted clock quality evaluation value is used as the optimal clock reference source, and the result is output.

[0028] Preferably, the calculation formula of the phase center offset is:

[0029]

[0030] in, Indicates the phase center deviation, Indicates The phase deviation between the clock reference source and the local oscillator in seconds. Indicates the recording duration.

[0031] Preferably, the calculation formulas for the weighted values ​​of any two clock reference sources are:

[0032]

[0033]

[0034] in, Indicates the clock reference source The weighted value of Indicates the clock reference source The weighted value of and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0035] Preferably, the calculation formula of the weighted clock quality evaluation value is:

[0036]

[0037]

[0038] in, and Respectively represent the clock reference source and clock reference source The weighted clock quality assessment value of Indicates GM accuracy, Indicates the GM deviation change value.

[0039] The present invention also provides an analysis system for the integrity of multiple reference source inputs in a time and frequency synchronization network, which is used to execute an analysis method for the integrity of multiple reference source inputs in a time and frequency synchronization network, so as to evaluate and select the best clock reference source among at least two clock reference sources in the time and frequency synchronization network, and the analysis system comprises:

[0040] A data acquisition module is used to acquire data sets of at least two clock reference sources to be compared in a time and frequency synchronization network, each data set including multiple attribute values ​​of the clock reference source to which it belongs;

[0041] An attribute comparison module, used to compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order;

[0042] A phase center deviation measurement module is used to measure the phase center deviation of multiple clock reference sources using a local rubidium atomic clock;

[0043] A weighted clock quality assessment module, used to calculate a weighted clock quality assessment value based on the comparison results of various attribute values ​​of the clock reference source and the phase center offset;

[0044] A local rubidium atomic clock taming module, used for selecting an optimal clock reference source based on the weighted clock quality evaluation value, and taming the local rubidium atomic clock using the optimal clock reference source;

[0045] The time-frequency output module is used to generate and output time signals and frequency signals according to the disciplined local rubidium atomic clock.

[0046] Compared with the related art, the method and system for analyzing the integrity of multiple reference source inputs of a time and frequency synchronization network provided by the present invention have the following beneficial effects:

[0047] The present invention can more accurately evaluate the timing accuracy of the reference source, ensure the accuracy of the reference source selection, and significantly improve the local timing level.

[0048] At the same time, the improved data set comparison algorithm enables the equipment that tracks the time and frequency synchronization network clock to accurately select the best master clock after the master clock of the time and frequency synchronization network has been transmitted through multiple stages and is subject to large background noise, thereby ensuring the precision and accuracy of the equipment's timing signal. The algorithm is also compatible with the original data set comparison algorithm, making it easy to transplant and quickly implement. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] Figure 1 A flowchart of a method for analyzing the integrity of multiple reference source inputs in a time-frequency synchronization network provided in Embodiment 1 of the present invention;

[0050] Figure 2 A module structure diagram of a system for analyzing the integrity of multiple reference source inputs of a time-frequency synchronization network provided in Embodiment 2 of the present invention;

[0051] Figure 3 A flowchart of a method for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network provided in Embodiment 3 of the present invention. DETAILED DESCRIPTION

[0052] The present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It is to be understood that the specific embodiments described herein are only used to explain the present invention, rather than to limit the present invention. It should also be noted that, for ease of description, only the parts related to the present invention, rather than all structures, are shown in the accompanying drawings. In addition, the embodiments of the present invention and the features in the embodiments may be combined with each other without conflict.

[0053] It should also be noted that, for ease of description, only the parts related to the present invention, but not all of the contents, are shown in the accompanying drawings. Before discussing the exemplary embodiments in more detail, it should be mentioned that some exemplary embodiments are described as processes or methods depicted as flow charts. Although the flow charts describe the operations (or steps) as sequential processes, many of the operations therein can be implemented in parallel, concurrently or simultaneously. In addition, the order of the operations can be rearranged. The process can be terminated when its operation is completed, but it can also have additional steps not included in the accompanying drawings. The process can correspond to methods, functions, procedures, subroutines, subprograms, etc.

[0054] Embodiment 1

[0055] The present invention provides a method for analyzing the integrity of multiple reference source inputs in a time-frequency synchronization network, which is used to evaluate and select the optimal clock reference source from at least two clock reference sources in the time-frequency synchronization network. Figure 1 As shown, the analysis method includes the following steps:

[0056] S1: In a time and frequency synchronization network, data sets of at least two clock reference sources to be compared are obtained, each data set including multiple attribute values ​​of the corresponding clock reference source.

[0057] S2: Compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order.

[0058] S3: Use the local rubidium atomic clock to measure the phase center deviation of multiple clock reference sources.

[0059] S4: A weighted clock quality assessment value is obtained based on the comparison results of the attribute values ​​of the clock reference source and the phase center deviation calculation.

[0060] When calculating the weighted clock quality evaluation value, the phase center offset is taken as an absolute value.

[0061] S5: Filter out an optimal clock reference source based on the weighted clock quality evaluation value, and use the optimal clock reference source to tame the local rubidium atomic clock.

[0062] S6: Generate and output time signal and frequency signal according to the disciplined local rubidium atomic clock.

[0063] Specifically, the attribute values ​​include GM identifier, GM priority I, GM level, GM accuracy, GM deviation change value, GM priority II and sender ID. This information is usually included in the message sent by the clock reference source. The required data set can be obtained by parsing the message.

[0064] When comparing attribute values, first compare whether the GM identifiers of any two clock reference sources are the same. If they are the same, it means that the received reference sources come from the same source, then enter the clock reference source hop count comparison processing flow; if they are not the same, it means that the received reference sources come from different sources, then compare the remaining attribute values ​​in sequence according to the predefined comparison order.

[0065] Specifically, when the GM identifiers are different, the predefined comparison order is:

[0066] a. Compare whether the values ​​of GM priority I of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority I, and take the clock reference source with a smaller GM priority I value as the optimal clock reference source, and output the result;

[0067] b. Compare whether the GM levels of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of the GM levels, and use the clock reference source with a smaller GM level as the optimal clock reference source, and output the result;

[0068] c. Compare whether the GM precision of any two clock reference sources is the same. If they are the same, proceed to the next step. If they are not the same, compare the numerical values ​​of the GM precision, and use the clock reference source with the smaller GM precision value as the optimal clock reference source, and output the result;

[0069] d. Compare the GM deviation change values ​​of any two clock reference sources to see if they are the same. If they are the same, proceed to the next step. If they are not the same, compare the GM deviation change values, and use the clock reference source with the smaller GM deviation change value as the optimal clock reference source, and output the result;

[0070] e. Compare whether the phase center offsets (absolute values) of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, calculate the weighted clock quality assessment value based on the attribute value of the clock reference source and the phase center offset, and output the result;

[0071] f. Compare whether the values ​​of GM priority II of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority II, and use the clock reference source with a smaller GM priority II value as the optimal clock reference source, and output the result;

[0072] g. Compare the values ​​of the transmitter IDs of any two clock reference sources, and use the clock reference source with the smaller transmitter ID value as the optimal clock reference source, and output the result.

[0073] In the time and frequency synchronization network, the attribute value of the clock reference source represents the quality of the original clock of the reference source. For example, if the value of GM priority I of one clock reference source is greater than the value of GM priority I of another clock reference source, it means that the original clock quality of the former clock reference source is worse than that of the latter clock reference source. Conversely, the original clock quality of the former clock reference source is better than that of the latter clock reference source.

[0074] Exemplary, predefined comparison order is:

[0075] GM identifier comparison: First, compare the GM identifiers of the two reference sources. If the GM identifiers are the same, continue to compare the subsequent attribute values; if the GM identifiers are different, enter the hop count comparison process, which is not within the scope of the present invention.

[0076] GM priority I comparison: If the GM identifiers are the same, compare the GM priority I values. If the GM priority I value of clock reference source C is greater than the GM priority I value of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0077] GM level comparison: If the GM priority levels I are the same, compare the GM level values. If the GM level value of clock reference source C is greater than the GM level value of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0078] GM accuracy comparison: If the GM levels are the same, compare the GM accuracy values. If the GM accuracy value of clock reference source C is greater than that of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0079] GM deviation change value comparison: If the GM accuracies are the same, compare the GM deviation change values. If the GM deviation change value of clock reference source C is greater than that of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0080] Phase center deviation comparison: If the GM deviation change values ​​are the same, the absolute value of the phase center deviation is taken for comparison, and the phase center deviation is taken as an absolute value because there is a distinction between positive and negative when measuring. If the phase center deviation of clock reference source C is greater than the phase center deviation of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0081] GM priority II comparison: If the phase center offsets are the same, compare the GM priority II values. If the GM priority II value of clock reference source C is greater than the GM priority II value of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0082] Transmitter ID comparison: If the GM priority IIs are the same, compare the transmitter ID values. If the transmitter ID value of clock reference source C is greater than the transmitter ID value of clock reference source D, it means that the original clock quality of clock reference source C is worse than that of clock reference source D; otherwise, the original clock quality of clock reference source C is better than that of clock reference source D.

[0083] Wherein, when the phase center offsets of any two clock reference sources are different, the method further includes the following steps:

[0084] Dynamically assign weighted values ​​of any two clock reference sources according to the magnitude of the phase center offset of the clock reference sources, and calculate a weighted clock quality assessment value according to the attribute values ​​and the phase center offset of the any two clock reference sources;

[0085] If the weighted clock quality evaluation values ​​are the same, the values ​​of GM priority II and sender ID are compared in turn, and the clock reference source with a smaller attribute value is used as the optimal clock reference source. If they are not the same, the weighted clock quality evaluation values ​​are compared, and the clock reference source with a smaller weighted clock quality evaluation value is used as the optimal clock reference source, and the result is output.

[0086] Specifically, the calculation formula of the phase center offset is:

[0087]

[0088] in, Indicates the phase center deviation, Indicates The phase deviation between the clock reference source and the local oscillator in seconds. Indicates the recording duration.

[0089] Specifically, the calculation formulas for the weighted values ​​of any two clock reference sources are:

[0090]

[0091]

[0092] in, Indicates the clock reference source The weighted value of Indicates the clock reference source The weighted value of and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0093] Specifically, the calculation formula of the weighted clock quality evaluation value is:

[0094]

[0095]

[0096] in, and Respectively represent the clock reference source and clock reference source The weighted clock quality assessment value of Indicates GM accuracy, Indicates the GM deviation change value.

[0097] Specifically, in steps S5 and S6, the optimal clock reference source is first determined, and the selected optimal reference source will be used to discipline the rubidium atomic clock. Based on the disciplined rubidium atomic clock, a high-precision time signal and frequency signal are generated and output.

[0098] The more specific steps are as follows:

[0099] According to the frequency information of the optimal clock reference source and the local rubidium atomic clock, the frequency of the local rubidium atomic clock is adjusted to keep consistent with the frequency of the optimal clock reference source.

[0100] According to the phase information of the optimal clock reference source and the local rubidium atomic clock, the phase of the local rubidium atomic clock is adjusted to keep it consistent with the phase of the optimal clock reference source. This is usually achieved through a phase-locked loop (PLL) or other phase synchronization technology.

[0101] After the local rubidium atomic clock is adjusted, a high-precision time signal is generated. The time signal is usually output in the form of pulses, such as a PPS (Pulse Per Second) signal once per second.

[0102] At the same time, a high-precision frequency signal is generated. The frequency signal is usually output in the form of a continuous sine wave or square wave, and the frequency value can be configured according to application requirements, such as 10 MHz, 5 MHz, etc.

[0103] Embodiment 2

[0104] The present invention also provides an analysis system for the integrity of multiple reference source inputs in a time-frequency synchronization network, which is used to evaluate and select the best clock reference source from at least two clock reference sources in the time-frequency synchronization network. Figure 2 As shown, the analysis system includes the following modules:

[0105] The data acquisition module 100 is used to acquire data sets of at least two clock reference sources to be compared in the time and frequency synchronization network, each data set including multiple attribute values ​​of the clock reference source to which it belongs;

[0106] The attribute comparison module 200 is used to compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order;

[0107] The phase center deviation measurement module 300 is used to measure the phase center deviation of multiple clock reference sources using a local rubidium atomic clock;

[0108] A weighted clock quality evaluation module 400, configured to calculate a weighted clock quality evaluation value based on the comparison results of the attribute values ​​of the clock reference source and the phase center offset;

[0109] A local rubidium atomic clock taming module 500, configured to select an optimal clock reference source based on the weighted clock quality evaluation value, and tame the local rubidium atomic clock using the optimal clock reference source;

[0110] The time-frequency output module 600 is used to generate and output a time signal and a frequency signal according to the disciplined local rubidium atomic clock.

[0111] Embodiment 3

[0112] In this embodiment, taking clock reference source A and clock reference source B as an example, the working process of the analysis method is described in detail below:

[0113] Step 1: Obtain data set A and data set B of reference source A and reference source B;

[0114] Step 2: Compare the GM identifiers in the messages of reference source A and reference source B. If GM identifier_A is equal to GM identifier_B, then enter the hop count comparison process of reference source A and reference source B. The hop count comparison process of reference source A and reference source B is not within the scope of the present invention. If GM identifier_A is not equal to GM identifier_B, then enter step 3. At the same time, the GM identifiers of the reference sources are recorded in GM identifier_A and GM identifier_B respectively;

[0115] Step 3: Compare the GM priority 1 values ​​in the messages of reference source A and reference source B (if GM identifier_A is not equal to GM identifier_B)

[0116] If the value of GM identifier_A is equal to the value of GM identifier_B, compare the GM priority I values ​​in the messages of reference source A and reference source B;

[0117] If the value of GM priority Ⅰ_A is equal to the value of GM priority Ⅰ_B, proceed to step 4;

[0118] If the value of GM priority Ⅰ_A is greater than the value of GM priority Ⅰ_B, then the result of B being better than A is returned;

[0119] If the value of GM priority Ⅰ_A is smaller than the value of GM priority Ⅰ_B, then the result that A is better than B is returned, and the value of GM priority Ⅰ of the reference source is recorded in GM priority Ⅰ_A and GM priority Ⅰ_B respectively.

[0120] Step 4: Compare the GM level values ​​in the messages of reference source A and reference source B (if GM priority I_A is equal to GM priority I_B).

[0121] If the value of GM priority I_A is equal to the value of GM priority I_B, compare the values ​​of the GM levels in the messages of reference source A and reference source B;

[0122] If the value of GM Level_A is equal to the value of GM Level_B, proceed to step 5;

[0123] If the value of GM Level_A is greater than the value of GM Level_B, then the result that B is better than A is returned;

[0124] If the value of GM Level_A is less than the value of GM Level_B, then the result that A is better than B is returned, and the values ​​of the GM Level of the reference source are recorded in GM Level_A and GM Level_B respectively.

[0125] Step 5: Compare the GM accuracy values ​​in the messages of reference source A and reference source B (if the value of GM Level_A is equal to the value of GM Level_B);

[0126] If the value of GM Level_A is equal to the value of GM Level_B, compare the GM accuracy values ​​in the messages of reference source A and reference source B;

[0127] If the value of GM precision_A is equal to the value of GM precision_B, proceed to step 6;

[0128] If the value of GM precision_A is greater than GM precision_B, then the result of B being better than A is returned;

[0129] If the value of GM Precision_A is smaller than the value of GM Precision_B, then the result that A is better than B is returned, and the values ​​of the GM Precision of the reference source are recorded in GM Precision_A and GM Precision_B respectively.

[0130] Step 6: Compare the GM deviation change values ​​in the messages of reference source A and reference source B (if the value of GM accuracy_A is equal to the value of GM accuracy_B);

[0131] If the value of GM Precision_A is equal to the value of GM Precision_B, compare the values ​​of the GM Deviation Change values ​​in the messages of Reference Source A and Reference Source B;

[0132] If the value of GM deviation change value_A is equal to the value of GM deviation change value_B, proceed to step 7;

[0133] If the value of GM Deviation Change Value_A is greater than the value of GM Deviation Change Value_B, then the result that B is better than A is returned;

[0134] If the value of GM deviation change value_A is less than the value of GM deviation change value_B, then the result of A being better than B is returned. At the same time, the values ​​of the GM accuracy of the reference source are recorded in GM deviation change value_A and GM deviation change value_B respectively.

[0135] Step 7: Calculate the phase center deviation of reference source A and reference source B respectively (if the value of GM deviation change value_A is equal to the value of GM deviation change value_B)

[0136] Since the time-frequency synchronization equipment is equipped with a local oscillator, which is generally a rubidium atomic clock, the local oscillator can be used as a local reference by recording the phase deviation between reference source A and the local oscillator, and the phase deviation between reference source B and the local oscillator in real time. Record for a period of time (n is the recording time, generally 3600 seconds or 7200 seconds), and then calculate the phase center deviation of reference source A. , calculate the phase center deviation of reference source B .

[0137] The calculation formula is as follows:

[0138]

[0139] in, Indicates The phase deviation between the clock reference source A and the local oscillator in seconds, Indicates the recording duration.

[0140] The calculation formula is as follows:

[0141]

[0142] in, Indicates The phase deviation between the clock reference source B and the local oscillator in seconds, Indicates the recording duration.

[0143] If the value of GM deviation change value_A is equal to the value of GM deviation change value_B, compare the absolute values ​​of the phase center deviations of reference source A and reference source B respectively;

[0144] If the value of the phase center offset_A is equal to the value of the phase center offset_B, the process proceeds to step eight;

[0145] If the value of phase center offset_A is greater than the value of phase center offset_B, then go to step ten;

[0146] If the value of phase center offset_A is smaller than the value of phase center offset_B, then enter step thirteen and record the value of the phase center offset change of the reference source into phase center offset_A and phase center offset_B respectively.

[0147] Step 8: Compare the GM priority II values ​​in the messages of reference source A and reference source B (if the value of phase center offset_A is equal to the value of phase center offset_B)

[0148] If the value of phase center offset_A is equal to the value of phase center offset_B, compare the GM priority II values ​​in the messages of reference source A and reference source B;

[0149] If the value of GM priority II_A is equal to the value of GM priority II_B, proceed to step nine;

[0150] If the value of GM priority Ⅱ_A is greater than the value of GM priority Ⅱ_B, then the result of B being better than A is returned;

[0151] If the value of GM priority II_A is less than that of GM priority II_B, then the result that A is better than B is returned, and the value of GM priority II of the reference source is recorded in GM priority II_A and GM priority II_B respectively;

[0152] Step 9: Compare the sender IDs in the messages of reference source A and reference source B (if the value of GM priority II_A is equal to the value of GM priority II_B)

[0153] If the value of GM priority II_A is equal to the value of GM priority II_B, compare the values ​​of the sender ID in the messages of reference source A and reference source B;

[0154] If the value of sender ID_A is greater than the value of sender ID_B, then the result that B is better than A is returned;

[0155] If the value of sender ID_A is less than the value of sender ID_B, then the result that A is better than B is returned, and the value of the sender ID of the reference source is recorded in sender ID_A and sender ID_B respectively;

[0156] Step 10: Dynamically assign weighted values ​​for clock reference source A and clock reference source B (if phase center offset_A is greater than phase center offset_B)

[0157] Weighted value of clock reference source A The calculation formula is as follows:

[0158]

[0159] in, Indicates the weighted value of clock reference source A. and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0160] Weighted value of clock reference source B The calculation formula is as follows:

[0161]

[0162] in, Indicates the weighted value of clock reference source B. and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0163] If the phase center offset_A is greater than the phase center offset_B, the dynamic weighted value of the clock reference source is recorded in the weighted value , weighted value , determine whether the result calculated in this step is a weighted value Greater than weighted value If the weighted values Greater than weighted value , proceed to step 11 to calculate and record the weighted clock quality assessment values ​​of clock reference source A and clock reference source B;

[0164] Step 11: Calculate and compare the weighted clock quality evaluation values ​​of clock reference source A and clock reference source B (if the weighted value Greater than weight ).

[0165] The weighted clock quality evaluation value of clock reference source A is calculated as follows:

[0166]

[0167] in, Indicates the weighted clock quality assessment value of clock reference source A. Indicates GM accuracy, Indicates the weight of GM deviation change value.

[0168] The weighted quality number of clock reference source B is calculated as follows:

[0169]

[0170] in, Indicates the weighted clock quality assessment value of clock reference source B. Indicates GM accuracy, Indicates the weight of GM deviation change value.

[0171] If the weighted value Greater than weight , compare the weighted clock quality evaluation values ​​of clock reference source A and clock reference source B. If the weighted clock quality evaluation value of clock reference source A is Weighted clock quality assessment value with clock reference source B If the values ​​are equal, proceed to step 12. Greater than the weighted clock quality assessment value , then return the result that B is better than A; if the weighted clock quality evaluation value Less than the weighted clock quality assessment value , then the result that A is better than B is returned. At the same time, the reference source weighted clock quality evaluation value is recorded in the weighted clock quality evaluation value , weighted clock quality assessment value .

[0172] Step 12: Compare the sender IDs in the messages of clock reference source A and clock reference source B (if the weighted clock quality assessment value Equal to the weighted clock quality assessment value ).

[0173] If the weighted clock quality assessment value of clock reference source A is Weighted clock quality assessment value with clock reference source B If the values ​​are equal, compare the values ​​of the sender ID in the messages of clock reference source A and clock reference source B. If the value of sender ID_A is greater than the value of sender ID_B, then return the result that B is better than A. If the value of sender ID_A is less than the value of sender ID_B, then return the result that A is better than B. At the same time, record the values ​​of the sender ID of the reference source to sender ID_A and sender ID_B respectively.

[0174] Step 13: Dynamically assign weighted values ​​for clock reference source A and clock reference source B (if phase center offset_A is less than phase center offset_B)

[0175] Weighted value of clock reference source A The calculation formula is as follows:

[0176]

[0177] in, Indicates the weighted value of clock reference source A. and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0178] Weighted value of clock reference source B The calculation formula is as follows:

[0179]

[0180] in, shows the weighted value of clock reference source B, and Respectively represent the clock reference source and clock reference source The phase center deviation.

[0181] If the phase center offset_A is less than the phase center offset_B, the dynamic weight value of the clock reference source is recorded in the weight value , weighted value , determine whether the result calculated in this step is a weighted value Less than weighted value If the weighted values Less than weighted value , proceed to step 14 to calculate and record the weighted clock quality assessment values ​​of reference sources A and B;

[0182] Step 14: Calculate and compare the weighted clock quality assessment values ​​of reference sources A and B (if the weighted value Less than weighted value )

[0183] The weighted clock quality evaluation value of clock reference source A is calculated as follows:

[0184]

[0185] in, Indicates the weighted clock quality assessment value of clock reference source A. Indicates GM accuracy, Indicates the weight of GM deviation change value.

[0186] The weighted clock quality evaluation value of clock reference source B is calculated as follows:

[0187]

[0188] in, Indicates the weighted clock quality assessment value of clock reference source B. Indicates GM accuracy, Indicates the weight of GM deviation change value.

[0189] If the weighted value Less than weighted value , compare the weighted clock quality evaluation values ​​of clock reference source A and clock reference source B. If the weighted clock quality evaluation value and weighted clock quality assessment value If the values ​​are equal, proceed to step 15. Greater than the weighted clock quality assessment value , then return the result that B is better than A; if the weighted clock quality evaluation value Less than the weighted clock quality assessment value , then the result that A is better than B is returned. At the same time, the weighted clock quality evaluation value of the reference source is recorded in the weighted clock quality evaluation value , weighted clock quality assessment value .

[0190] Step 15: Compare the sender IDs in the messages of reference source A and reference source B (if the weighted clock quality assessment value Equal to the weighted clock quality assessment value )

[0191] If the weighted clock quality assessment value and weighted clock quality assessment value If the values ​​are equal, compare the values ​​of the sender ID in the messages of clock reference source A and clock reference source B. If the value of sender ID_A is greater than the value of sender ID_B, then return the result that B is better than A. If the value of sender ID_A is less than the value of sender ID_B, then return the result that A is better than B. At the same time, record the values ​​of the sender ID of the reference source to sender ID_A and sender ID_B respectively.

[0192] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0193] A person skilled in the art may understand that all or part of the steps in the various methods of the above embodiments may be completed by instructing related hardware through a program, and the program may be stored in a computer-readable storage medium, the storage medium including a read-only memory (ROM), a random access memory (RAM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), a one-time programmable read-only memory (OTPROM), an electronically-erasable programmable read-only memory (EEPROM), a compact disc (CD-ROM) or other optical disc storage, magnetic disk storage, magnetic tape storage, or any other computer-readable medium that can be used to carry or store data.

[0194] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

Claims

1. A method for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network, used to evaluate and select the best clock reference source from at least two clock reference sources in a time and frequency synchronization network, characterized in that: The analytical method comprises the following steps: In a time and frequency synchronization network, data sets of at least two clock reference sources to be compared are obtained, each data set including multiple attribute values ​​of the clock reference source; Compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order; The phase center deviation of multiple clock reference sources is measured using a local rubidium atomic clock. A weighted clock quality assessment value is obtained based on the comparison results of the attribute values ​​of the clock reference source and the phase center deviation calculation; Screening out an optimal clock reference source based on the weighted clock quality evaluation value, and using the optimal clock reference source to tame a local rubidium atomic clock; Generate and output time signal and frequency signal according to the disciplined local rubidium atomic clock; The attribute values ​​include GM identifier, GM priority I, GM level, GM accuracy, GM deviation change value, GM priority II and sender ID; When comparing attribute values, first compare whether the GM identifiers of any two clock reference sources are the same. If they are the same, then enter the clock reference source hop count comparison process; if they are not the same, then compare the remaining attribute values ​​in sequence according to the predefined comparison order; When the GM identifiers are different, the predefined comparison order is: a. Compare whether the values ​​of GM priority I of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority I, and take the clock reference source with a smaller GM priority I value as the optimal clock reference source, and output the result; b. Compare whether the GM levels of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of the GM levels, and use the clock reference source with a smaller GM level as the optimal clock reference source, and output the result; c. Compare whether the GM precision of any two clock reference sources is the same. If they are the same, proceed to the next step. If they are not the same, compare the numerical values ​​of the GM precision, and use the clock reference source with the smaller GM precision value as the optimal clock reference source, and output the result; d. Compare the GM deviation change values ​​of any two clock reference sources to see if they are the same. If they are the same, proceed to the next step. If they are not the same, compare the GM deviation change values, and use the clock reference source with the smaller GM deviation change value as the optimal clock reference source, and output the result; e. Compare whether the phase center offsets of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, calculate the weighted clock quality assessment value based on the attribute value of the clock reference source and the phase center offset, and output the result; f. Compare whether the values ​​of GM priority II of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority II, and use the clock reference source with a smaller GM priority II value as the optimal clock reference source, and output the result; g. Compare the values ​​of the transmitter IDs of any two clock reference sources, and use the clock reference source with the smaller transmitter ID value as the optimal clock reference source, and output the result.

2. According to claim 1, a method for analyzing the integrity of multiple reference source inputs of a time and frequency synchronization network is characterized in that: When the phase center offsets of any two clock reference sources are different, the method further includes the following steps: Dynamically assign weighted values ​​of the arbitrary two clock reference sources according to the magnitude of the phase center offset of the clock reference source, and calculate a weighted clock quality assessment value according to the attribute values ​​and the phase center offset of the arbitrary two clock reference sources; If the weighted clock quality evaluation values ​​are the same, the values ​​of GM priority II and sender ID are compared in turn, and the clock reference source with a smaller attribute value is used as the optimal clock reference source. If they are not the same, the weighted clock quality evaluation values ​​are compared, and the clock reference source with a smaller weighted clock quality evaluation value is used as the optimal clock reference source, and the result is output.

3. The method for analyzing the integrity of multiple reference source inputs of a time and frequency synchronization network according to claim 1, characterized in that: The calculation formula of the phase center deviation is: in, Indicates the phase center deviation, Indicates The phase deviation between the clock reference source and the local oscillator in seconds. Indicates the recording duration.

4. The method for analyzing the integrity of multiple reference source inputs of a time-frequency synchronization network according to claim 2, characterized in that: The calculation formulas for the weighted values ​​of any two clock reference sources are: in, Indicates the clock reference source The weighted value of Indicates the clock reference source The weighted value of and Respectively represent the clock reference source and clock reference source The phase center deviation.

5. The method for analyzing the integrity of multiple reference source inputs of a time and frequency synchronization network according to claim 2, characterized in that: The calculation formula of the weighted clock quality evaluation value is: in, and Respectively represent the clock reference source and clock reference source The weighted clock quality assessment value of Indicates GM accuracy, Indicates the GM deviation change value.

6. A system for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network, used to execute a method for analyzing the integrity of multiple reference source inputs in a time and frequency synchronization network as described in any one of claims 1 to 5, so as to evaluate and select the best clock reference source among at least two clock reference sources in the time and frequency synchronization network, characterized in that: The analysis system comprises: A data acquisition module is used to acquire data sets of at least two clock reference sources to be compared in a time and frequency synchronization network, each data set including multiple attribute values ​​of the clock reference source to which it belongs; The attribute comparison module is used to compare the attribute values ​​of the clock reference sources in sequence according to a predefined comparison order; A phase center deviation measurement module is used to measure the phase center deviation of multiple clock reference sources using a local rubidium atomic clock; A weighted clock quality assessment module, used to calculate a weighted clock quality assessment value based on the comparison results of various attribute values ​​of the clock reference source and the phase center offset; A local rubidium atomic clock taming module, used to screen out an optimal clock reference source based on the weighted clock quality evaluation value, and tame the local rubidium atomic clock using the optimal clock reference source; The time-frequency output module is used to generate and output time signals and frequency signals according to the disciplined local rubidium atomic clock; The attribute values ​​include GM identifier, GM priority I, GM level, GM accuracy, GM deviation change value, GM priority II and sender ID; The attribute comparison module is also used to compare whether the GM identifiers of any two clock reference sources are the same when comparing the attribute values. If they are the same, the hop count comparison processing flow of the clock reference source is entered; if they are not the same, the remaining attribute values ​​are compared in sequence according to a predefined comparison order; When the GM identifiers are different, the predefined comparison order is: a. Compare whether the values ​​of GM priority I of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority I, and take the clock reference source with a smaller GM priority I value as the optimal clock reference source, and output the result; b. Compare whether the GM levels of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of the GM levels, and use the clock reference source with a smaller GM level as the optimal clock reference source, and output the result; c. Compare whether the GM precision of any two clock reference sources is the same. If they are the same, proceed to the next step. If they are not the same, compare the numerical values ​​of the GM precision, and use the clock reference source with the smaller GM precision value as the optimal clock reference source, and output the result; d. Compare the GM deviation change values ​​of any two clock reference sources to see if they are the same. If they are the same, proceed to the next step. If they are not the same, compare the GM deviation change values, and use the clock reference source with the smaller GM deviation change value as the optimal clock reference source, and output the result; e. Compare whether the phase center offsets of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, calculate the weighted clock quality assessment value based on the attribute value of the clock reference source and the phase center offset, and output the result; f. Compare whether the values ​​of GM priority II of any two clock reference sources are the same. If they are the same, proceed to the next step. If they are not the same, compare the values ​​of GM priority II, and use the clock reference source with a smaller GM priority II value as the optimal clock reference source, and output the result; g. Compare the values ​​of the transmitter IDs of any two clock reference sources, and use the clock reference source with the smaller transmitter ID value as the optimal clock reference source, and output the result.

Images