A timing system interference detection method and a timing system

Abstract

The embodiment of the application discloses a timing system interference detection method and a timing system. The method comprises the following steps: determining an actual clock difference curve of the timing system according to clock difference data of the timing system and a standard clock system; calculating a timing deviation factor, a timing deviation speed factor and a timing deviation offset acceleration factor of the timing system respectively based on the actual clock difference curve; and performing interference detection on the timing system according to the timing deviation factor, the timing deviation speed factor and the timing deviation offset acceleration factor to obtain an interference situation of the timing system. Through the application, the technical problem that the interference situation of the timing system cannot be effectively detected in the related art is solved, the interference situation of the timing system can be detected in real time and accurately, and the technical effect that the time information output by the timing system can be effectively prevented from being cheated is achieved.

Description

Technical Field

[0001] This invention relates to the field of time signal processing technology, and in particular to a method for detecting interference in a timing system and a timing system. Background Technology

[0002] Time is a fundamental element in all stages of information acquisition, transmission, fusion, and application. Accurate and reliable acquisition and determination of time information by timing systems is a prerequisite for the stable operation of critical national sectors such as power, transportation, finance, communications, and national defense. Time-using devices rely on their own timing systems to synchronize time information through external time sources or to maintain time autonomously using internal clocks. If external time sources are susceptible to interference, it can lead to time rejection or time spoofing; furthermore, internal clocks, due to their inherent clock drift characteristics, cannot maintain accurate and stable time for extended periods. However, current technologies cannot effectively detect interference with timing systems, resulting in the output time information being spoofed, and the timing system failing to acquire accurate and reliable time information.

[0003] No effective solution has yet been proposed to address the above problems. Summary of the Invention

[0004] This invention provides a timing system interference detection method and a timing system, to at least solve the technical problem in related technologies that cannot effectively detect interference in timing systems.

[0005] According to one aspect of the present invention, a timing system interference detection method is provided, comprising: determining the actual clock difference curve of the timing system based on clock difference data between the timing system and a standard clock system; calculating a timing deviation factor, a timing skew factor, and a timing offset acceleration factor of the timing system based on the actual clock difference curve; and performing interference detection on the timing system based on the timing deviation factor, the timing skew factor, and the timing offset acceleration factor to obtain the interference status of the timing system.

[0006] Optionally, calculating the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system based on the actual clock difference curve includes: processing the actual clock difference curve to obtain the timing deviation fluctuation, timing skew difference, and timing offset acceleration of the timing system; calculating the timing deviation factor of the timing system based on the timing deviation fluctuation; calculating the timing skew factor of the timing system based on the timing skew difference; and calculating the timing offset acceleration factor of the timing system based on the timing offset acceleration.

[0007] Optionally, the actual clock error curve is processed to obtain the timing deviation fluctuation of the timing system, including: using a singular value removal function to perform proximity smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; obtaining a calibration clock error model of the timing system under interference-free conditions; calculating the difference between the calibration clock error model and the processing result, and using an absolute value function to calculate the absolute value of the difference; and using a mean function to process the absolute value of the difference to obtain the timing deviation fluctuation of the timing system.

[0008] Optionally, obtaining the calibration clock difference model of the timing system under interference-free conditions includes: obtaining first clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the first clock difference data is the clock difference data between the timing system and the standard clock system collected when calibrating with a stable external timing signal; determining the clock difference calibration curve of the timing system based on the first clock difference data; and using a distortion smoothing template to... The clock error calibration curve is subjected to distortion smoothing processing to obtain the calibration time clock error model.

[0009] Optionally, the actual clock error curve is processed to obtain the timing skew difference of the timing system, including: using a singular value removal function to perform proximity smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; obtaining an autonomous timing clock error model of the timing system under interference-free conditions; using a difference function to process the processing result and the autonomous timing clock error model respectively to obtain the clock error skew of the timing system; and using a mean function to process the clock error skew to obtain the timing skew difference of the timing system.

[0010] Optionally, obtaining the autonomous time difference model of the timing system under interference-free conditions includes: obtaining second clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the second clock difference data is the clock difference data between the timing system and the standard clock system collected when the external time signal is disconnected; determining the clock difference drift curve of the timing system based on the second clock difference data; and performing distortion smoothing processing on the clock difference drift curve using a distortion smoothing template to obtain the autonomous time difference model.

[0011] Optionally, the actual clock error curve is processed to obtain the timing offset acceleration of the timing system, including: using a singular value removal function to perform neighbor smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; using a quadratic difference function to process the processing result to obtain the clock error offset acceleration of the timing system; and using a mean function to process the clock error offset acceleration to obtain the timing offset acceleration of the timing system.

[0012] Optionally, interference detection is performed on the timing system based on the timing deviation factor, the timing speed deviation factor, and the timing offset acceleration factor to obtain the interference status of the timing system, including: establishing a quality degradation model of the timing system, wherein the quality degradation model is used to identify the interference status of the timing system; inputting the timing deviation factor, the timing speed deviation factor, and the timing offset acceleration factor into the quality degradation model, and outputting detection values ​​from the quality degradation model; and determining the interference status of the timing system based on the detection values.

[0013] Optionally, determining the interference status of the timing system based on the detection value includes: when the detection value is 0, it indicates that the timing system is not interfered with; when the detection value is 1, it indicates that the timing system is subjected to time synchronization rejection interference; when the detection value is 2, it indicates that the timing system is subjected to time synchronization spoofing interference.

[0014] According to another aspect of the present invention, a timing system is also provided, comprising: a processor and a memory for storing processor-executable instructions, wherein the processor is configured to perform the timing system interference detection method described in any one of the above embodiments.

[0015] In this embodiment of the invention, the actual clock difference curve of the timing system is determined based on the clock difference data between the timing system and the standard clock system. Based on the actual clock difference curve, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated. Interference detection is then performed on the timing system based on these factors to determine the interference status of the timing system. In other words, this embodiment of the invention uses the clock difference data between the timing system and the standard clock system to generate the actual clock difference curve of the timing system. Then, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated using the actual clock difference curve. Finally, the interference status of the timing system is detected based on these factors, thereby solving the technical problem in related technologies where interference of the timing system cannot be effectively detected. This achieves the technical effect of being able to detect the interference status of the timing system in real time and accurately, thus effectively preventing the time information output by the timing system from being deceived. Attached Figure Description

[0016] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 A flowchart of a timing system interference detection method provided in an embodiment of the present invention; Figure 2 This is a schematic diagram of a timing system interference detection device provided in an embodiment of the present invention. Detailed Implementation

[0017] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0018] It should be noted that the terms "first," "second," etc., in the specification, claims, and drawings of this invention are used to distinguish different objects, rather than to limit a specific order.

[0019] According to one aspect of the present invention, a timing system interference detection method is provided. It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions. Furthermore, although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.

[0020] Figure 1 The flowchart of the timing system interference detection method provided in the embodiments of the present invention is as follows: Figure 1 As shown, the method includes the following steps: Step S102: Determine the actual clock difference curve of the timing system based on the clock difference data between the timing system and the standard clock system; Step S104: Calculate the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system based on the actual clock difference curve. Step S106: Perform interference detection on the timing system based on the timing deviation factor, timing speed deviation factor, and timing offset acceleration factor to obtain the interference status of the timing system.

[0021] The interference situations of the aforementioned timing system include, but are not limited to, no interference, interference caused by time synchronization rejection, and interference caused by time synchronization deception.

[0022] Optionally, the timing system's timing deviation, timing skew, and timing offset acceleration are analyzed using actual clock difference curves during operation. The timing deviation factor, timing skew factor, and timing offset acceleration factor are calculated. These factors are then used to comprehensively analyze the interference situation of the timing system, and a warning is issued when interference occurs. In other words, this invention combines the operating characteristics of the timing system with its own clock characteristics to analyze in real time whether the timing system is subject to timekeeping rejection interference or timekeeping deception interference, and provides warnings. This prevents the time information of time-using devices from being maliciously deceived. It can be applied to timing interference detection in complex scenarios and can also provide support for research on timing anti-interference technology.

[0023] It should be noted that when the interference to the timing system is analyzed in real time, it prompts the use of backup methods to enhance timing and resist interference according to actual timing needs, thereby preventing the time information output by the timing system from being deceived, and providing necessary support for the timing system to obtain time information accurately and reliably and maintain time synchronization for a long time in complex environments.

[0024] In this embodiment of the invention, the actual clock difference curve of the timing system is determined based on the clock difference data between the timing system and the standard clock system. Based on the actual clock difference curve, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated. Interference detection is then performed on the timing system based on these factors to determine the interference status of the timing system. In other words, this embodiment of the invention uses the clock difference data between the timing system and the standard clock system to generate the actual clock difference curve of the timing system. Then, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated using the actual clock difference curve. Finally, the interference status of the timing system is detected based on these factors, thereby solving the technical problem in related technologies where interference of the timing system cannot be effectively detected. This achieves the technical effect of being able to detect the interference status of the timing system in real time and accurately, thus effectively preventing the time information output by the timing system from being deceived.

[0025] In one optional embodiment, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated based on the actual clock difference curve, including: processing the actual clock difference curve to obtain the timing deviation fluctuation, timing skew difference, and timing offset acceleration of the timing system; calculating the timing deviation factor of the timing system based on the timing deviation fluctuation; calculating the timing skew factor of the timing system based on the timing skew difference; and calculating the timing offset acceleration factor of the timing system based on the timing offset acceleration.

[0026] Optionally, based on the timing deviation fluctuation W of the timing system d Calculate the timing deviation factor T of the timing system. d as follows: Among them, TR d Indicates the timing bias differential threshold; optionally, TR d =50ns. When At that time, T d =0 indicates that the timing deviation of the timing system is within the allowable range; when At that time, T d =1 indicates that the timing deviation of the timing system exceeds the allowable range.

[0027] Optionally, based on the timing skew difference W of the timing system s Calculate the timing skew factor T of the timing system. s as follows:

[0028] Among them, when W s When ≤0, T s =0, indicating that the timing skew rate of the timing system does not exceed its own clock drift rate; when W s When T > 0, s =1 indicates that the timing skew rate of the timing system exceeds its own clock drift rate.

[0029] Optionally, based on the timing offset acceleration W of the timing system a Calculate the timing offset acceleration factor T of the timing system. a as follows:

[0030] Among them, when W a When ≥0, T a =0 indicates that the timing eccentricity of the timing system remains constant or increases; when W a When <0, T a =1 indicates that the timing skew speed of the timing system is getting slower and slower.

[0031] In one optional embodiment, the actual clock error curve is processed to obtain the timing deviation fluctuation of the timing system, including: using a singular value removal function to perform proximity smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; obtaining a calibration clock error model of the timing system under interference-free conditions; calculating the difference between the calibration clock error model and the processing result, and using an absolute value function to calculate the absolute value of the difference; and using a mean function to process the absolute value of the difference to obtain the timing deviation fluctuation of the timing system.

[0032] In actual operation of the timing system, clock difference data between the timing system and the standard clock system within a time interval of T seconds is collected in real time at 1-second intervals to obtain the actual clock difference curve R of the timing system. t Alternatively, T = 600.

[0033] Furthermore, the timing deviation fluctuation W of the timing system is calculated. d as follows: W d =mean(abs(M) c -clip(R t ,5%))) Among them, M c This represents the calibration clock difference model of the timing system when it operates without interference. `clip` is the singular value removal function, which represents the expression for the actual clock difference curve R of the timing system. t The largest and smallest 5% of the singular values ​​are processed by nearest-neighbor smoothing, where abs represents the absolute value function and mean represents the mean function.

[0034] In one optional embodiment, obtaining the calibration time difference model of the timing system under interference-free conditions includes: obtaining first clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the first clock difference data is the clock difference data between the timing system and the standard clock system collected when a stable external timing signal is connected for calibration; determining the clock difference calibration curve of the timing system based on the first clock difference data; and using a distortion smoothing template to perform distortion smoothing processing on the clock difference calibration curve to obtain the calibration time difference model.

[0035] Under interference-free conditions, the 24-hour clock difference between the timing system and the standard clock system is sampled in 1-second increments, using the standard clock system as a reference.

[0036] Under the condition of real-time calibration with access to a stable external time signal, the 24-hour clock difference data between the timing system and the standard clock system (corresponding to the first clock difference data mentioned above) is sampled in 1-second increments to obtain the clock difference calibration curve C of the timing system. t Furthermore, through distortion smoothing, a calibration time difference model M for the timing system under interference-free operation is established. c as follows: M c =convsld(C t ,fs) Where fs = (0.03, 0.07, 0.15, 0.25, 0, 0.25, 0.15, 0.07, 0.03) represents the distortion smoothing template, convsld represents the sliding convolution function, and the clock error calibration curve C is calibrated using the distortion smoothing template fs. t Perform distortion smoothing.

[0037] In one optional embodiment, the actual clock error curve is processed to obtain the timing skew difference of the timing system, including: using a singular value removal function to perform proximity smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; obtaining an autonomous timing clock error model of the timing system under interference-free conditions; using a difference function to process the processing result and the autonomous timing clock error model respectively to obtain the clock error skew of the timing system; and using a mean function to process the clock error skew to obtain the timing skew difference of the timing system.

[0038] In actual operation of the timing system, clock difference data between the timing system and the standard clock system within a time interval of T seconds is collected in real time at 1-second intervals to obtain the actual clock difference curve R of the timing system. t Alternatively, T = 600.

[0039] Furthermore, the timing skew difference W of the timing system is calculated. s as follows: W s =mean(diff(clip(R) t ,5%))-diff(M s )) Among them, M s This represents the autonomous clock bias model of the timing system when it operates without interference. `clip` is the singular value removal function, representing the transformation of the actual clock bias curve R of the timing system. t The largest and smallest 5% of the singular values ​​are smoothed by proximity processing. diff is the difference function, which is used to calculate the clock bias, and mean is used to calculate the mean function.

[0040] In one optional embodiment, obtaining the autonomous time difference model of the timing system under interference-free conditions includes: obtaining second clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the second clock difference data is the clock difference data between the timing system and the standard clock system collected when the external time signal is disconnected; determining the clock difference drift curve of the timing system based on the second clock difference data; and using a distortion smoothing template to perform distortion smoothing processing on the clock difference drift curve to obtain the autonomous time difference model.

[0041] Under interference-free conditions, the 24-hour clock difference between the timing system and the standard clock system is sampled in 1-second increments, using the standard clock system as a reference.

[0042] With the external time signal disconnected, the 24-hour clock difference data between the timing system and the standard clock system (corresponding to the second clock difference data mentioned above) is sampled in 1-second increments to obtain the clock difference drift curve S of the timing system. t And through distortion Smoothing processing is used to establish an autonomous time-difference model M for the timing system to operate without interference. s as follows: M s =convsld(S t ,fs) Where fs = (0.03, 0.07, 0.15, 0.25, 0, 0.25, 0.15, 0.07, 0.03) represents the distortion smoothing template, convsld represents the sliding convolution function, and the distortion smoothing template fs is used to smooth the clock drift curve S. t Perform distortion smoothing.

[0043] In one optional embodiment, the actual clock error curve is processed to obtain the timing offset acceleration of the timing system, including: using a singular value removal function to perform neighbor smoothing on a predetermined maximum singular value and a predetermined minimum singular value in the actual clock error curve to obtain a processing result; using a quadratic difference function to process the processing result to obtain the clock error offset acceleration of the timing system; and using a mean function to process the clock error offset acceleration to obtain the timing offset acceleration of the timing system.

[0044] In actual operation of the timing system, clock difference data between the timing system and the standard clock system within a time interval of T seconds is collected in real time at 1-second intervals to obtain the actual clock difference curve R of the timing system. t Alternatively, T = 600.

[0045] Furthermore, the timing offset acceleration Wa of the timing system is calculated as follows: W a =mean(diff2(clip(R) t ,5%))) Here, clip is the singular value removal function, which represents the function that removes singular values ​​from the actual clock error curve R of the timing system. t The largest and smallest 5% of the singular values ​​are smoothed by proximity processing. diff2 is a quadratic difference function that calculates the clock offset acceleration, and mean is a function that calculates the mean.

[0046] In one optional embodiment, interference detection is performed on the timing system based on the timing deviation factor, timing speed deviation factor, and timing offset acceleration factor to obtain the interference status of the timing system. This includes: establishing a quality degradation model of the timing system, wherein the quality degradation model is used to identify the interference status of the timing system; inputting the timing deviation factor, timing speed deviation factor, and timing offset acceleration factor into the quality degradation model, and outputting detection values ​​from the quality degradation model; and determining the interference status of the timing system based on the detection values.

[0047] It should be noted that when the timing system operates normally under interference-free conditions, its internal clock performs real-time calibration by receiving external time signals to achieve precise time synchronization with the external time source, thus outputting an accurate time. When the external time signal is interrupted, the timing system uses its internal clock for autonomous time maintenance, and the output time shifts due to the inherent clock drift characteristic of the internal clock. Under prolonged interruption of the external time signal, the cumulative deviation of the timing system's output time will exceed the allowable range, which is considered to be time signal rejection interference. When the received external time signal is a spurious signal, the internal clock will be gradually pulled off course according to the deception paradigm of the spurious signal, causing the timing system's output time to be induced to deviate, which is considered to be time signal deception interference.

[0048] Furthermore, a quality degradation model M for the timing system is constructed. t as follows:

[0049] Among them, T d T represents the timing deviation factor of the timing system. s T represents the timing skew factor of the timing system. a This represents the timing offset acceleration factor of the timing system.

[0050] Furthermore, when M t When M = 0, it indicates that the timing system is not disturbed; when M = 0, it indicates that the timing system is not disturbed. t When M = 1, it indicates that the timing system is experiencing timing rejection interference; when M = 1, it indicates that the timing system is experiencing timing rejection interference. t When the value is 2, it indicates that the timing system is being interfered with by time synchronization fraud.

[0051] In one optional embodiment, determining the interference status of the timing system based on the detection value includes: when the detection value is 0, it indicates that the timing system is not interfered with; when the detection value is 1, it indicates that the timing system is subjected to time synchronization rejection interference; and when the detection value is 2, it indicates that the timing system is subjected to time synchronization spoofing interference.

[0052] In actual operation of the timing system, clock difference data between the timing system and the standard clock system is collected in real time within a time interval of T seconds, with an interval of X seconds and a unit of 1 second, to obtain the actual clock difference curve R of the timing system. t Alternatively, X = 1800, T = 600.

[0053] Furthermore, the timing deviation factor T of the timing system is calculated separately. d Timing eccentricity factor T s and timing offset acceleration factor T a Then calculate the timing system quality degradation model M. t Then according to M t The value is analyzed in real time to assess the interference situation of the timing system, as follows: When T d When = 0, M t =0 indicates that the timing deviation of the timing system is within the allowable range, meaning that the timing system has not been disturbed.

[0054] When T d =1, and T s When = 0, M t =1 indicates that the timing deviation of the timing system exceeds the allowable range, and the timing deviation rate of the timing system does not exceed its own clock drift rate. That is, the timing system is subjected to time synchronization rejection interference and immediately issues a prompt.

[0055] When T d =1, and T s =1, and T a When = 0, M t =2 indicates that the timing deviation of the timing system exceeds the allowable range, and the timing deviation speed of the timing system exceeds its own clock drift rate. Furthermore, the timing deviation speed of the timing system remains constant or increases, meaning that the timing system is being subjected to timekeeping deception interference, and an immediate warning is issued.

[0056] According to one aspect of the present invention, a timing system interference detection device is provided. Figure 2 This is a schematic diagram of a timing system interference detection device provided in an embodiment of the present invention, as shown below. Figure 2 As shown, the timing system interference detection device includes: a determination module 22, a calculation module 24, and a detection module 26. The timing system interference detection device will now be described in detail.

[0057] The determination module 22 is used to determine the actual clock difference curve of the timing system based on the clock difference data between the timing system and the standard clock system; The calculation module 24, connected to the determination module 22, is used to calculate the timing deviation factor, timing skew factor and timing offset acceleration factor of the timing system based on the actual clock difference curve. The detection module 26 is connected to the calculation module 24 and is used to perform interference detection on the timing system based on the timing deviation factor, timing speed factor and timing offset acceleration factor to obtain the interference status of the timing system.

[0058] It should be noted that the above-mentioned determining module 22, calculation module 24 and detection module 26 correspond to steps S102 to S106 in the method embodiment. The examples and application scenarios implemented by the above modules and corresponding steps are the same, but are not limited to the content disclosed in the above method embodiment.

[0059] In the above embodiments of the present invention, the timing system interference detection device can generate an actual clock difference curve of the timing system using clock difference data between the timing system and the standard clock system. Then, it can calculate the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system using the actual clock difference curve. Based on the timing deviation factor, timing skew factor, and timing offset acceleration factor, it can detect the interference of the timing system. This solves the technical problem in related technologies that it is impossible to effectively detect the interference of the timing system, and achieves the technical effect of being able to detect the interference of the timing system in real time and accurately, thereby effectively avoiding the time information output by the timing system being deceived.

[0060] In an optional embodiment, the calculation module 24 includes: a first processing unit for processing the actual clock error curve to obtain the timing deviation fluctuation, timing skew speed difference, and timing offset acceleration of the timing system; a first calculation unit for calculating the timing deviation factor of the timing system based on the timing deviation fluctuation; a second calculation unit for calculating the timing skew speed factor of the timing system based on the timing skew speed difference; and a third calculation unit for calculating the timing offset acceleration factor of the timing system based on the timing offset acceleration.

Claims

1. A method for detecting interference in a timing system, characterized in that, include: The actual clock difference curve of the timing system is determined based on the clock difference data between the timing system and the standard clock system; Based on the actual clock error curve, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated respectively. A quality degradation model for the timing system is established, wherein the quality degradation model is used to identify interference conditions affecting the timing system; the quality degradation model of the timing system for: in, This represents the timing deviation factor of the timing system. This represents the timing skew factor of the timing system. This represents the timing offset acceleration factor of the timing system; The timing deviation factor, the timing speed deviation factor, and the timing offset acceleration factor are input into the quality degradation model, and the detection value is output from the quality degradation model. The interference status of the timing system is determined based on the detected values.

2. The method according to claim 1, characterized in that, Based on the actual clock error curve, the timing deviation factor, timing skew factor, and timing offset acceleration factor of the timing system are calculated, including: The actual clock error curve is processed to obtain the timing deviation fluctuation, timing skew speed difference, and timing offset acceleration of the timing system. The timing deviation factor of the timing system is calculated based on the timing deviation fluctuation. The timing skew factor of the timing system is calculated based on the timing skew difference. The timing offset acceleration factor of the timing system is calculated based on the timing offset acceleration.

3. The method according to claim 2, characterized in that, The actual clock error curve is processed to obtain the timing deviation fluctuation of the timing system, including: The singular value removal function is used to perform proximity smoothing on the predetermined maximum and minimum singular values ​​in the actual clock error curve to obtain the processing result; Obtain the calibration time difference model of the timing system under interference-free conditions; Calculate the difference between the calibration time difference model and the processing result, and use the absolute value function to calculate the absolute value of the difference; The absolute value of the difference is processed using a mean function to obtain the timing deviation fluctuation of the timing system.

4. The method according to claim 3, characterized in that, Obtaining the calibration time difference model of the timing system under interference-free conditions includes: Acquire first clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the first clock difference data is the clock difference data between the timing system and the standard clock system collected when calibrating by accessing a stable external timing signal; The clock error calibration curve of the timing system is determined based on the first clock error data. The clock error calibration curve is smoothed using a distortion smoothing template to obtain the calibration time-of-flight clock error model.

5. The method according to claim 2, characterized in that, The actual clock error curve is processed to obtain the timing skew difference value of the timing system, including: The singular value removal function is used to perform proximity smoothing on the predetermined maximum and minimum singular values ​​in the actual clock error curve to obtain the processing result; Obtain the autonomous time difference model of the timing system under interference-free conditions; The clock bias speed of the timing system is obtained by processing the processing results and the autonomous timer clock bias model using a difference function. The clock skew rate is processed using a mean function to obtain the timing skew rate difference of the timing system.

6. The method according to claim 5, characterized in that, Obtaining the autonomous time difference model of the timing system under interference-free conditions includes: Acquire second clock difference data between the timing system and the standard clock system under interference-free conditions, wherein the second clock difference data is the clock difference data between the timing system and the standard clock system collected when the external timing signal is disconnected; The clock drift curve of the timing system is determined based on the second clock drift data. The clock drift curve is smoothed using a distortion smoothing template to obtain the autonomous time-steering clock drift model.

7. The method according to claim 2, characterized in that, The actual clock error curve is processed to obtain the timing offset acceleration of the timing system, including: The singular value removal function is used to perform proximity smoothing on the predetermined maximum and minimum singular values ​​in the actual clock error curve to obtain the processing result; The clock offset acceleration of the timing system is obtained by processing the processing result using a quadratic difference function. The clock offset acceleration is processed using a mean function to obtain the timing offset acceleration of the timing system.

8. The method according to claim 1, characterized in that, Determining the interference status of the timing system based on the detected values ​​includes: When the detection value is 0, it indicates that the timing system has not been disturbed; When the detection value is 1, it indicates that the timing system is subject to timing rejection interference; When the detection value is 2, it indicates that the timing system is being interfered with by timing deception.

9. A timing system, characterized in that, include: A processor and a memory for storing processor-executable instructions, wherein the processor is configured to perform the timing system interference detection method according to any one of claims 1 to 8.

Images