CN110531380A - The device and method of satellite clock source low amplitude persistent anomaly for identification - Google Patents
The device and method of satellite clock source low amplitude persistent anomaly for identification Download PDFInfo
- Publication number
- CN110531380A CN110531380A CN201910813061.2A CN201910813061A CN110531380A CN 110531380 A CN110531380 A CN 110531380A CN 201910813061 A CN201910813061 A CN 201910813061A CN 110531380 A CN110531380 A CN 110531380A
- Authority
- CN
- China
- Prior art keywords
- frequency
- time interval
- satellite clock
- pps
- crystal oscillator
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000002085 persistent effect Effects 0.000 title claims abstract description 17
- 239000013078 crystal Substances 0.000 claims abstract description 61
- 238000005259 measurement Methods 0.000 claims abstract description 27
- 238000009825 accumulation Methods 0.000 claims abstract description 18
- 230000002159 abnormal effect Effects 0.000 claims description 9
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000010200 validation analysis Methods 0.000 abstract description 2
- 230000032683 aging Effects 0.000 description 5
- 230000002045 lasting effect Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 241001123248 Arma Species 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/02—Details of the space or ground control segments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/25—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
- G01S19/256—Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Electric Clocks (AREA)
Abstract
A kind of device and method of source of satellite clock for identification low amplitude persistent anomaly is to be based on the scale frequency drift discreet value of constant-temperature crystal oscillator second by CPU to issue control signal control frequency division coefficient control module generation frequency dividing control parameter；Frequency divider adjusts the frequency division coefficient corresponding with nominal count value of its setting according to frequency dividing control parameter and generates PPS；The calculating of time interval difference is carried out using satellite clock signal 1PPS and above-mentioned PPS of the time interval measurement module to input；Time interval difference is compared with threshold value 1us again, satellite clock exception is then judged if more than A；Otherwise, continuous accumulation operation is carried out to time interval difference until judging satellite clock source exception when accumulating operation values >=preset threshold B.For the present invention by carrying out outside this second threshold decision to satellite clock source 1PPS, the validation checking under another rise time scale can effectively identify duration synchronization attack, it is ensured that each real-time control system time synchronization safety.
Description
Technical field
The present invention relates to a kind of methods for detecting satellite clock source state, and in particular to a kind of satellite clock source for identification
The device and method of low amplitude persistent anomaly.
Background technique
Domestic Beidou satellite navigation system (BDS) and american global positioning system (GPS) can provide round-the-clock satellite time transfer
Signal is widely used in the infrastructure such as electric power, broadcast, traffic, real-time industry control.The reliability in satellite clock source is always to close
The Important Problems of note.As national network space militarization trend is aggravated, Global Satellite Navigation System (including DBS and GPS) is latent
Security risk come true threat.Civilian satellite channel is used except the military outer synchronization system plenty of time, shortage is recognized
The guarantee for demonstrate,proving encryption, provides opportunity to malicious attacker --- lure that time synchronization fills by sending false satellite-signal
The time signal for generating and having certain deviation is set, and then influences goal systems normal function.
To enhance each system time source reliability, China mainly uses the mode based on space-based time service, supplemented by ground time service
Accurate clock synchronization is carried out, ensures to select reliable clock reference by the clock deviation comparison to different clocks source signal and local clock.
In addition, the validity for having correlative study to judge satellite clock source by predetermined threshold value, avoids because of satellite clock source anomalous effects
Time service, but duration synchronization attack is not fully taken into account, it can control exceptional value in threshold value, by constantly tired
Long-pending mode causes clock deviation to be more than system permission, and therefore, how to detect identification to continue satellite clock source is to study herein extremely
Emphasis.
Constant-temperature crystal oscillator has the characteristics that short-term stability is high, due to crystal oscillator aging drift phenomenon, crystal oscillator under long time scale
Frequency will appear the error of lasting accumulation.The a large amount of history numbers accumulated in stage test are being tamed and kept based on constant-temperature crystal oscillator
According to establishing age compensation model and predict aging characteristics trend, can be calculated according to different time sections deviation corresponding
Compensation rate compensate, so that constant-temperature crystal oscillator is provided relatively reliable stable precise frequency.
Summary of the invention
For the reliability for ensuring the real-time industrials control system time synchronization process Satellite clock source such as electric power, realize to defending
The lasting effectiveness of star clock source detects, and the present invention provides a kind of dresses of satellite clock source low amplitude persistent anomaly for identification
It sets and method, constructs the standard time reference hierarchy based on High Accuracy Constant Temperature crystal oscillator by mending crystal oscillator frequency drift (frequency drift) in advance,
Validation checking is carried out to satellite clock source under long time scale, can effectively identify and be defended caused by duration synchronization attack
Star clock source is abnormal.
In order to solve the above technical problems, the technical solution adopted by the present invention is, a kind of satellite clock source for identification is by a narrow margin
Spend the device of persistent anomaly comprising CPU, frequency division coefficient control module, frequency divider, time interval measurement module and register,
Wherein: the CPU issues control signal to frequency division coefficient control module based on the discreet value of constant-temperature crystal oscillator second scale frequency drift；The frequency dividing
Coefficient control module is connect with CPU, for receiving the control signal of CPU sending, and generates frequency dividing control parameter, while will be generated
Frequency dividing control parameter be sent to frequency divider；The frequency divider is connect with frequency division coefficient control module and constant-temperature crystal oscillator respectively, is used for
It receives the practical high-frequency count signal generated of constant-temperature crystal oscillator and receives the frequency dividing control parameter that frequency division coefficient control module generates, and
The constant-temperature crystal oscillator set the corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances is corrected according to the frequency dividing control coefficient；Basis is repaired simultaneously
The practical high-frequency count signal generated of frequency division coefficient and constant-temperature crystal oscillator after just generates PPS, then the PPS of generation is sent to the time
Interval measurement module；The time interval measurement module is connect with the frequency divider and satellite receiver respectively, is connect for receiving satellite
The PPS that the satellite clock signal 1PPS and frequency divider that receipts machine is sent are sent, and the survey of time interval difference is carried out to 1PPS and PPS
Amount, while the time interval difference obtained after measurement is sent to register；The register respectively with time interval measurement module
It is connected with CPU, for receiving the time interval difference of the 1PPS and PPS, and it is default to judge whether the time interval difference is greater than
Otherwise threshold value A, carries out time interval difference from the time interval difference continuous if more than satellite clock exception is then judged
Operation is accumulated, until accumulation operation values >=preset threshold B, judges satellite clock exception, otherwise satellite clock is normal at this time.
The register is also used to the normal or abnormal judging result in satellite clock source being sent to CPU to carry out in next step
Processing.
Present invention simultaneously provides a kind of method using above-mentioned apparatus identification satellite clock source low amplitude persistent anomaly, the party
Steps are as follows for method:
Step 1: estimating the frequency drift of constant-temperature crystal oscillator, crystal oscillator frequency drift is reflected on second scale；
Step 2:CPU issues control signal to frequency division coefficient control module according to the crystal oscillator frequency drift on the second scale estimated,
Frequency division coefficient control module generates frequency dividing control parameter after receiving the control signal that CPU is issued, while by the frequency dividing control of generation
Parameter is sent to frequency divider；
Step 3: frequency divider receives the frequency dividing control parameter that frequency division coefficient control module generates, according to the frequency dividing control coefficient
The constant-temperature crystal oscillator set the corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances is corrected, and according to revised frequency division coefficient and constant temperature
The practical high-frequency count signal generated of crystal oscillator generates PPS, while the PPS of generation is sent to time interval measurement module；
Step 4: it is poor that time interval measurement module carries out time interval to the satellite clock signal 1PPS and above-mentioned PPS of input
Value measurement, and the time interval difference is sent to register；
Step 5: register receives the time interval difference of the 1PPS and PPS, and judges whether the time interval difference is big
In preset threshold 1us, if more than satellite clock exception is then judged, otherwise, further between the time from the time interval difference
Continuous accumulation operation is carried out every difference until accumulation operation values >=preset threshold B, judges satellite clock exception, otherwise satellite at this time
Clock is normal.
The process of the frequency division coefficient according to frequency dividing control coefficient amendment setting referred in above-mentioned steps 3 are as follows: setting point
Frequency control parameter Nerr, then the aborning frequency division coefficient of frequency divider is adjusted to Nrev=Nstd+Nerr, wherein NstdJust for constant-temperature crystal oscillator
The corresponding frequency division coefficient of pulse per second (PPS) in normal situation.
The above-mentioned preset threshold B referred to is to set according to different real-time device systems to time precision prescribed.
The solution have the advantages that: for satellite clock source persistent anomaly caused by duration synchronization attack, we
Method eliminates the accumulated error of constant-temperature crystal oscillator by establishing constant-temperature crystal oscillator frequency-drift compensation model, and then by adjusting in frequency divider
Frequency division coefficient generates PPS after compensation, carries out the inspection of the validity under this second threshold decision and long time scale to 1PPS by PPS
It surveys.It is low that the simulation experiment result shows that this method can effectively identify that satellite clock source caused by duration synchronization attack continues
Amplitude is abnormal.
Detailed description of the invention
Fig. 1 is the functional block diagram of satellite clock source persistent anomaly for identification of the invention；
Fig. 2 is recognition methods flow chart of the invention；
Fig. 3 is the actual measurement of constant-temperature crystal oscillator frequency drift and prognostic chart；
Fig. 4 is 1PPS random error distribution map when satellite clock source is normal and abnormal；
Fig. 5 is the accumulated error effect picture that satellite clock source is normally identified with present invention when exception.
Specific embodiment
For the reliability for ensuring the real-time industrials control system time synchronization process Satellite clock source such as electric power, realize to defending
The lasting effectiveness of star clock source detects, and the invention proposes a kind of dresses of satellite clock source low amplitude persistent anomaly for identification
It sets, in conjunction with referring to Fig. 1 comprising CPU, frequency division coefficient control module, frequency divider, time interval measurement module and register,
In:
The CPU is for issuing control according to the crystal oscillator frequency drift on the second scale estimated after the crystal oscillator frequency drift on second scale is estimated
Signal processed is to frequency division coefficient control module；
The frequency division coefficient control module is connect with CPU, for receiving the control signal of CPU sending, and generates frequency dividing control
Parameter, while the frequency dividing control parameter of generation is sent to frequency divider；
The frequency divider is connect with frequency division coefficient control module and constant-temperature crystal oscillator respectively, and for receiving, constant-temperature crystal oscillator is practical to be generated
High-frequency count signal (f) and receive the frequency dividing control parameter that frequency division coefficient control module generates, and according to the frequency dividing control system
Number corrects the constant-temperature crystal oscillator set the corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances；Simultaneously according to revised frequency division coefficient and
The practical high-frequency count signal generated of constant-temperature crystal oscillator generates PPS, then the PPS of generation is sent to time interval measurement module；
The time interval measurement module is connect with the frequency divider and satellite receiver respectively, for receiving satellite receiver hair
The PPS that the satellite clock signal 1PPS and frequency divider sent is sent, and the measurement of time interval difference is carried out to 1PPS and PPS, simultaneously
The time interval difference obtained after measurement is sent to register；
The register is connect with time interval measurement module and CPU respectively, for receiving the time interval of the 1PPS and PPS
Difference, and judge whether the time interval difference is greater than preset threshold A i.e. 1us, it is no if more than satellite clock exception is then judged
Then, continuous accumulation operation further is carried out to time interval difference from the time interval difference, until accumulation operation values >=pre-
If threshold value B, satellite clock exception is judged at this time, otherwise satellite clock is normal.The register is also used to judging result (when satellite
The normal or abnormal result in clock source) it is sent to CPU and is further processed.
The present invention using above-mentioned apparatus carry out satellite clock source low amplitude persistent anomaly recognition methods specific steps such as
Under, it ties referring to fig. 2:
Step 1: estimating the frequency drift of constant-temperature crystal oscillator, crystal oscillator frequency drift is reflected on second scale；
Step 2:CPU issues control signal to frequency division coefficient control module according to the crystal oscillator frequency drift on the second scale estimated,
Frequency division coefficient control module generates frequency dividing control parameter after receiving the control signal that CPU is issued, while by the frequency dividing control of generation
Parameter is sent to frequency divider；
Step 3: frequency divider receives the frequency dividing control parameter that frequency division coefficient control module generates, according to the frequency dividing control coefficient
The constant-temperature crystal oscillator set the corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances is corrected, and according to revised frequency division coefficient and constant temperature
The practical high-frequency count signal generated of crystal oscillator generates PPS, while the PPS of generation is sent to time interval measurement module；
Step 4: it is poor that time interval measurement module carries out time interval to the satellite clock signal 1PPS and above-mentioned PPS of input
Value measurement, and the time interval difference is sent to register；
Step 5: register receives the time interval difference of the 1PPS and PPS, and judges whether the time interval difference is big
In preset threshold A, if more than satellite clock exception is then judged, otherwise, further to time interval from the time interval difference
Difference carries out continuous accumulation operation, until accumulation operation values >=preset threshold B, judges satellite clock exception, otherwise satellite at this time
Clock is normal.
It is next that the register is also used to for judging result (the normal or abnormal result in satellite clock source) being sent to CPU progress
Step processing.
The predictor method of the frequency drift of the above-mentioned constant-temperature crystal oscillator referred to is being tamed and is being kept for the stage based on constant-temperature crystal oscillator
The a large amount of historical datas accumulated in test, using the methods of empirical model (logarithmic model, arma modeling etc.), neural network to accidentally
Difference measures.By taking logarithmic model as an example, specifically:
F (t)=A+Bln (t+C)
In formula, A, B, C are ageing parameter.A, the value of B, C can be estimated by historical data.F (t) is that t moment constant temperature is brilliant
Vibration frequency.This method is the prior art, specifically reference can be made to " the aging modeling and its intelligence benefit based on crystal oscillator of Xue Chong
Repay " [D] Xian Electronics Science and Technology University, 2014.
In addition, crystal oscillator frequency drift is assigned to if being estimated with day scale or hour scale to constant-temperature crystal oscillator frequency drift
On second scale.Specifically: constant-temperature crystal oscillator frequency is set as f0(unit HZ), crystal oscillator frequency drift is f under day scaleday, then compensation per second
fday/86400。
The above-mentioned frequency division coefficient process according to frequency dividing control coefficient amendment setting referred to are as follows: setting frequency dividing control parameter
Nerr, then the aborning frequency division coefficient of frequency divider is adjusted to Nrev=Nstd+Nerr, wherein NstdFor constant-temperature crystal oscillator second under normal circumstances
The corresponding frequency division coefficient of pulse, numerical value are equal to crystal oscillator frequency；The constant-temperature crystal oscillator refers to that constant-temperature crystal oscillator is unaged under normal circumstances
In the case where.
The generation process of the above-mentioned PPS referred to is the ordinary skill in the art.
The above-mentioned preset threshold A referred to is 1us；Preset threshold B then wants the time according to different real-time device control systems
Refinement degree is set, and such as us, ms and s grades.
Embodiment 1
Continue the feasible of the method and apparatus of low amplitude exception for the verifying source of satellite clock for identification of the present invention
Property, emulation experiment verifying is carried out using matlab.Mean value is used to generate satellite clock for zero random function simulation in experiment
Random error, σ=50ns.The random function simulation for being 20 by using mean value generates lasting low amplitude satellite clock and believes extremely
Number.
As shown in figure 3, to tame and keeping a large amount of historical datas accumulated in stage test according to practical constant-temperature crystal oscillator
Curve graph with prediction clock deviation is surveyed to constant-temperature crystal oscillator frequency drift；The present invention carries out the pre- of constant-temperature crystal oscillator frequency drift using logarithmic model
It surveys.Wherein, curve 1 is actual measurement crystal oscillator frequency drift curve, and curve 2 is prediction aging curve, and curve 3 is curve after compensation, Cong Tuke
Know, the foundation for realizing standard time reference source can be compensated to crystal oscillator by the method for estimating crystal oscillator frequency drift.In addition, with people
The development of work intelligence, more more adaptable algorithm estimation effect may be more preferable, herein only with logarithmic model example.
Fig. 4 is satellite clock source normally with 1PPS random error distribution map when exception, and (a) is satellite clock source normal condition
The distribution of lower random error, (b) distribution of random error is descended in satellite clock source extremely.Fig. 5 is that satellite clock source is normal and abnormal
When the accumulated error effect picture that identifies of the present invention, curve 1 is that PPS and 1PPS long-time clock deviation are missed under normal circumstances in satellite clock source
Poor accumulated value (clock deviation error accumulation value is indicated in Fig. 5 with accumulating clock deviation), curve 2 is PPS under the abnormal conditions of satellite clock source
With 1PPS long-time clock deviation error accumulation value.It as seen from the figure, is 0 just since 1PPS obeys mean value when satellite clock source is normal
State distribution, accumulation clock deviation value (the accumulation operation values carried out to time interval difference referred in step 5) is in long time scale
Lower stabilization is intended to a steady state value；When the exception of satellite clock source, accumulation clock deviation value will be more and more big, as the accumulation of time will
It can be more than threshold value.By simulation result it is found that this method can effectively identify that satellite clock source continues the situation of low amplitude exception.
Claims (6)
1. a kind of device of satellite clock source low amplitude persistent anomaly for identification, which is characterized in that it includes CPU, frequency dividing system
Number control module, frequency divider, time interval measurement module and register, in which:
The CPU is based on the scale frequency drift discreet value of constant-temperature crystal oscillator second and issues control signal to frequency division coefficient control module；
The frequency division coefficient control module is connect with CPU, for receiving the control signal of CPU sending, and generates frequency dividing control parameter,
The frequency dividing control parameter of generation is sent to frequency divider simultaneously；
The frequency divider is connect with frequency division coefficient control module and constant-temperature crystal oscillator respectively, for receiving the practical height generated of constant-temperature crystal oscillator
Frequency count signal and the frequency dividing control parameter for receiving the generation of frequency division coefficient control module, and set according to frequency dividing control coefficient amendment
The fixed constant-temperature crystal oscillator corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances；Simultaneously according to revised frequency division coefficient and constant-temperature crystal oscillator
The high-frequency count signal actually generated generates PPS, then the PPS of generation is sent to time interval measurement module；
The time interval measurement module is connect with the frequency divider and satellite receiver respectively, for receiving satellite receiver transmission
The PPS that satellite clock signal 1PPS and frequency divider are sent, and the measurement of time interval difference is carried out to 1PPS and PPS, while will survey
The time interval difference obtained after amount is sent to register；
The register is connect with time interval measurement module and CPU respectively, and the time interval for receiving the 1PPS and PPS is poor
Value, and judge whether the time interval difference is greater than preset threshold A i.e. 1us, if more than satellite clock exception is then judged, otherwise,
To the continuous accumulation operation of time interval difference progress up to accumulating operation values >=preset threshold B from the time interval difference, this
When judge satellite clock exception, otherwise satellite clock is normal.
2. a kind of method of satellite clock source low amplitude persistent anomaly for identification as described in claim 1, which is characterized in that
The preset threshold B is to set according to different real-time device systems to time precision prescribed.
3. the device of satellite clock source low amplitude persistent anomaly for identification as described in claim 1, which is characterized in that described
Register is also used to the normal or abnormal judging result in satellite clock source being sent to CPU to be further processed.
4. a kind of method using device identification satellite clock source low amplitude persistent anomaly described in claims 1 or 22 or 3,
It is characterized in that, the method steps are as follows:
Step 1: estimating the frequency drift of constant-temperature crystal oscillator, crystal oscillator frequency drift is reflected on second scale；
Step 2:CPU issues control signal to frequency division coefficient control module, frequency dividing according to the crystal oscillator frequency drift on the second scale estimated
Coefficient control module generates frequency dividing control parameter after receiving the control signal that CPU is issued, while by the frequency dividing control parameter of generation
It is sent to frequency divider；
Step 3: frequency divider receives the frequency dividing control parameter that frequency division coefficient control module generates, and is corrected according to the frequency dividing control coefficient
The constant-temperature crystal oscillator of the setting corresponding frequency division coefficient of pulse per second (PPS) under normal circumstances, and according to revised frequency division coefficient and constant-temperature crystal oscillator
The high-frequency count signal actually generated generates PPS, while the PPS of generation is sent to time interval measurement module；
Step 4: time interval measurement module carries out the survey of time interval difference to the satellite clock signal 1PPS and above-mentioned PPS of input
Amount, and the time interval difference is sent to register；
Step 5: register receives the time interval difference of the 1PPS and PPS, and it is pre- to judge whether the time interval difference is greater than
If threshold value 1us, if more than satellite clock exception is then judged, otherwise, time interval difference is carried out from the time interval difference
Continuous accumulation operation judges satellite clock exception, otherwise satellite clock is normal until accumulation operation values >=preset threshold B at this time.
5. a kind of method of satellite clock source low amplitude persistent anomaly for identification as claimed in claim 4, which is characterized in that
The process of the frequency division coefficient according to frequency dividing control coefficient amendment setting referred in the step 3 are as follows: setting frequency dividing control ginseng
Number Nerr, then the aborning frequency division coefficient of frequency divider is adjusted to Nrev=Nstd+Nerr, wherein NstdUnder normal circumstances for constant-temperature crystal oscillator
The corresponding frequency division coefficient of pulse per second (PPS).
6. a kind of method of satellite clock source low amplitude persistent anomaly for identification as claimed in claim 4, which is characterized in that
Preset threshold B in the step 5 is to set according to different real-time device systems to time precision prescribed.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910813061.2A CN110531380A (en) | 2019-08-30 | 2019-08-30 | The device and method of satellite clock source low amplitude persistent anomaly for identification |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910813061.2A CN110531380A (en) | 2019-08-30 | 2019-08-30 | The device and method of satellite clock source low amplitude persistent anomaly for identification |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110531380A true CN110531380A (en) | 2019-12-03 |
Family
ID=68665509
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910813061.2A Pending CN110531380A (en) | 2019-08-30 | 2019-08-30 | The device and method of satellite clock source low amplitude persistent anomaly for identification |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110531380A (en) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111600649A (en) * | 2020-06-02 | 2020-08-28 | 中国科学院国家授时中心 | Frequency signal on-line monitoring switching equipment |
| CN111711448A (en) * | 2020-08-07 | 2020-09-25 | 电信科学技术第五研究所有限公司 | Rubidium atomic clock taming system and method |
| CN111787013A (en) * | 2020-07-03 | 2020-10-16 | 中国电子科技集团公司第三十研究所 | Attack monitoring device and method of optical fiber time synchronization system |
| CN112637137A (en) * | 2020-12-08 | 2021-04-09 | 中国电子科技集团公司第三十研究所 | Optical fiber time synchronization monitoring method and system based on clock error dynamic model |
| CN113419286A (en) * | 2021-06-21 | 2021-09-21 | 中国科学院地质与地球物理研究所 | High-precision frequency synthesis compensation method for complementation of satellite and constant-temperature crystal oscillator |
| CN113848568A (en) * | 2021-09-24 | 2021-12-28 | 上海精密计量测试研究所 | Time correction system and method |
| US20220174629A1 (en) * | 2020-11-30 | 2022-06-02 | Viettel Group | Method and apparatus for data frame synchronization of 5g base station |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008019605A1 (en) * | 2006-08-11 | 2008-02-21 | Huawei Technologies Co., Ltd. | A method and device for improving satellite time synchronization pulse hold performance |
| CN202759413U (en) * | 2012-09-17 | 2013-02-27 | 南京澳德思电气有限公司 | An on-line time error correcting system based on a satellite and a constant temperature crystal oscillator |
| CN104300969A (en) * | 2014-05-12 | 2015-01-21 | 长沙理工大学 | High-precision synchronization clock realization method based on full digital phase-locked loop |
| CN105259561A (en) * | 2015-11-24 | 2016-01-20 | 郑州威科姆科技股份有限公司 | Device and method for testing exception handling ability of satellite time synchronization device |
| CN105446128A (en) * | 2015-12-24 | 2016-03-30 | 张宏伟 | Navigation satellite high-precision time service system with constant temperature circuit and method |
| CN108809463A (en) * | 2017-04-28 | 2018-11-13 | 工业和信息化部电信研究院 | A kind of multiple reference clock weighting synthetic method |
| CN109557577A (en) * | 2019-01-16 | 2019-04-02 | 合肥国为电子有限公司 | Suitable for the seismographic clock-generating device of node and method |
| CN110133997A (en) * | 2019-05-17 | 2019-08-16 | 长沙理工大学 | A method of detection satellite clock is abnormal |
| CN110133998A (en) * | 2019-05-17 | 2019-08-16 | 长沙理工大学 | A kind of method of anti-duration synchronization attack |
-
2019
- 2019-08-30 CN CN201910813061.2A patent/CN110531380A/en active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008019605A1 (en) * | 2006-08-11 | 2008-02-21 | Huawei Technologies Co., Ltd. | A method and device for improving satellite time synchronization pulse hold performance |
| CN202759413U (en) * | 2012-09-17 | 2013-02-27 | 南京澳德思电气有限公司 | An on-line time error correcting system based on a satellite and a constant temperature crystal oscillator |
| CN104300969A (en) * | 2014-05-12 | 2015-01-21 | 长沙理工大学 | High-precision synchronization clock realization method based on full digital phase-locked loop |
| CN105259561A (en) * | 2015-11-24 | 2016-01-20 | 郑州威科姆科技股份有限公司 | Device and method for testing exception handling ability of satellite time synchronization device |
| CN105446128A (en) * | 2015-12-24 | 2016-03-30 | 张宏伟 | Navigation satellite high-precision time service system with constant temperature circuit and method |
| CN108809463A (en) * | 2017-04-28 | 2018-11-13 | 工业和信息化部电信研究院 | A kind of multiple reference clock weighting synthetic method |
| CN109557577A (en) * | 2019-01-16 | 2019-04-02 | 合肥国为电子有限公司 | Suitable for the seismographic clock-generating device of node and method |
| CN110133997A (en) * | 2019-05-17 | 2019-08-16 | 长沙理工大学 | A method of detection satellite clock is abnormal |
| CN110133998A (en) * | 2019-05-17 | 2019-08-16 | 长沙理工大学 | A kind of method of anti-duration synchronization attack |
Non-Patent Citations (7)
| Title |
|---|
| 李泽文 等: "基于全数字锁相环的电力系统高精度同步时钟" * |
| 李泽文 等: "基于全数字锁相环的电力系统高精度同步时钟", 《电力自动化设备》 * |
| 李泽文 等: "基于数字锁相原理的GPS 高精度同步时钟产生新方法" * |
| 田禹泽 等: "基于数字锁相环的星载光谱仪本地时钟源设计", 《电子与信息学报》 * |
| 程骏: "基于北斗卫星的故障行波定位系统的设计研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
| 蒋陆萍 等: "基于GPS 实现电力系统高精度同步时钟" * |
| 蒋陆萍 等: "基于GPS 实现电力系统高精度同步时钟", 《电网技术》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111600649A (en) * | 2020-06-02 | 2020-08-28 | 中国科学院国家授时中心 | Frequency signal on-line monitoring switching equipment |
| CN111787013A (en) * | 2020-07-03 | 2020-10-16 | 中国电子科技集团公司第三十研究所 | Attack monitoring device and method of optical fiber time synchronization system |
| CN111787013B (en) * | 2020-07-03 | 2022-02-25 | 中国电子科技集团公司第三十研究所 | Attack monitoring device and method of optical fiber time synchronization system |
| CN111711448A (en) * | 2020-08-07 | 2020-09-25 | 电信科学技术第五研究所有限公司 | Rubidium atomic clock taming system and method |
| US20220174629A1 (en) * | 2020-11-30 | 2022-06-02 | Viettel Group | Method and apparatus for data frame synchronization of 5g base station |
| US11683771B2 (en) * | 2020-11-30 | 2023-06-20 | Viettel Group | Method and apparatus for data frame synchronization of 5G base station |
| CN112637137A (en) * | 2020-12-08 | 2021-04-09 | 中国电子科技集团公司第三十研究所 | Optical fiber time synchronization monitoring method and system based on clock error dynamic model |
| CN112637137B (en) * | 2020-12-08 | 2022-02-25 | 中国电子科技集团公司第三十研究所 | Optical fiber time synchronization monitoring method and system based on clock error dynamic model |
| CN113419286A (en) * | 2021-06-21 | 2021-09-21 | 中国科学院地质与地球物理研究所 | High-precision frequency synthesis compensation method for complementation of satellite and constant-temperature crystal oscillator |
| CN113419286B (en) * | 2021-06-21 | 2022-01-14 | 中国科学院地质与地球物理研究所 | High-precision frequency synthesis compensation method for complementation of satellite and constant-temperature crystal oscillator |
| CN113848568A (en) * | 2021-09-24 | 2021-12-28 | 上海精密计量测试研究所 | Time correction system and method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110531380A (en) | The device and method of satellite clock source low amplitude persistent anomaly for identification | |
| CN110133997B (en) | Method for detecting satellite clock abnormity | |
| CN110133998B (en) | Method for preventing duration time synchronization attack | |
| KR101930354B1 (en) | Apparatus and method for detecting deception signal in global navigation satellite receiver | |
| CN203164620U (en) | High precision time synchronization device | |
| CN110058287A (en) | A kind of Orbit determination for LEOs method, apparatus and system | |
| US10310091B2 (en) | GPS-based time stamp system | |
| US10274601B2 (en) | Communications system | |
| CN103616814A (en) | Synchronous sampling clock closed loop correcting method and system based on FPGA | |
| CN107247405A (en) | A kind of satellite navigation time service is kept time the solution in clock system | |
| EP3043493B1 (en) | Frequency shift detection method and device | |
| CN103913987B (en) | GPS time dissemination systems and its method for obtaining accurate time reference | |
| CN109085616B (en) | Satellite time service method, device and storage medium | |
| CN104702214B (en) | A kind of method of crystal oscillator frequency compensation | |
| WO2022001368A1 (en) | Radar and radar system | |
| US12061223B2 (en) | Method and system for traveling wave-based fault location of AC and DC transmission lines based on 5G communication | |
| AU2009309734B2 (en) | Method for correcting prediction errors of signal values with time variation subjected to interference by various uncontrollable systematic effects | |
| US11456808B2 (en) | Data receiver and method of identifying an operational state of a data transmitter by monitoring environment parameters | |
| CN102004441A (en) | Adaptive crystal oscillator frequency timekeeping method | |
| WO2020149103A1 (en) | Monitoring system and synchronization method | |
| CN110138491A (en) | GPS clock compensation method based on generalized weighted least square method | |
| CN104808480A (en) | Pulse per second (PPS) generating method and device | |
| US11209469B2 (en) | Timestamp inconsistency and shift detection for synchrophasor data based on correlation between relative phase angle and frequency | |
| CN103116172B (en) | False detection method of satellite navigation signal and satellite navigation positioning receiver | |
| CN109100927A (en) | Without timing system time autonomous synchronous method under reference conditions of tracing to the source |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20191203 |