CN114818247A - Atomic time calculation method and device based on hydrogen atomic clock drift prediction - Google Patents
Atomic time calculation method and device based on hydrogen atomic clock drift prediction Download PDFInfo
- Publication number
- CN114818247A CN114818247A CN202111669122.6A CN202111669122A CN114818247A CN 114818247 A CN114818247 A CN 114818247A CN 202111669122 A CN202111669122 A CN 202111669122A CN 114818247 A CN114818247 A CN 114818247A
- Authority
- CN
- China
- Prior art keywords
- clock
- difference
- drift
- atomic
- time
- 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
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 title claims abstract description 72
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 72
- 239000001257 hydrogen Substances 0.000 title claims abstract description 72
- 238000004364 calculation method Methods 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000001914 filtration Methods 0.000 claims description 17
- 230000008569 process Effects 0.000 abstract description 12
- 230000007774 longterm Effects 0.000 abstract description 8
- 230000006870 function Effects 0.000 description 17
- 238000004422 calculation algorithm Methods 0.000 description 9
- 229910052792 caesium Inorganic materials 0.000 description 5
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 description 5
- 238000004590 computer program Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 125000004429 atom Chemical group 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03L—AUTOMATIC CONTROL, STARTING, SYNCHRONISATION OR STABILISATION OF GENERATORS OF ELECTRONIC OSCILLATIONS OR PULSES
- H03L7/00—Automatic control of frequency or phase; Synchronisation
- H03L7/26—Automatic control of frequency or phase; Synchronisation using energy levels of molecules, atoms, or subatomic particles as a frequency reference
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application discloses an atomic time calculation method and device based on hydrogen atomic clock drift prediction, wherein the method comprises the following steps: acquiring historical clock error data from a database; obtaining a clock difference initial value between the combined clock and the reference clock and a clock difference initial value between each clock and the combined clock according to the clock difference data of each atomic clock; performing time function fitting according to the clock difference data and the drift values between the clocks and the combined clock to obtain a frequency difference drift predicted value changing along with time; and deducting a frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale. The method and the device solve the problem of long-time frequency drift in atomic time scale calculation based on the hydrogen atomic clock, thereby realizing real-time frequency drift updating of the atomic clock in real time in the atomic time scale calculation process, accurately deducting a drift item, and finally obtaining the atomic time scale with both short-term stability and long-term stability.
Description
Technical Field
The application relates to the field of atomic clocks, in particular to an atomic time calculation method and device based on hydrogen atomic clock drift prediction.
Background
Atomic time computing is a time keeping method used in time keeping systems. A time keeping system typically consists of a plurality of atomic clocks, each of which generates an atomic time stamp. Because each atomic clock has the possibility of abnormal conditions, the atomic time scale of the combined clock is calculated by collecting the clock difference data of each atomic clock and utilizing a statistical method, so that the atomic time scale has higher stability and plays a key role in establishing the atomic time scale of the time keeping system.
Atomic clocks in a timekeeping system generally comprise a hydrogen atomic clock and a cesium atomic clock, wherein the cesium atomic clock and the hydrogen atomic clock are different in that the short-term stability of the hydrogen atomic clock is better than that of the cesium atomic clock, but the long-term stability of the hydrogen atomic clock is worse than that of the cesium atomic clock because the hydrogen atomic clock has severe frequency drift in the long-term operation process. The atomic time algorithms that are currently used mainly include ALGOS, AT1, and Kalman algorithm. The Kalman algorithm is based on a Kalman filter, and a classical weight calculation mode is not used. The ALGOS algorithm only considers the change influence of the frequency difference of the atomic clock in the atomic time calculation process, and is more suitable for atomic time calculation based on the cesium atomic clock. The AT1 algorithm considers the problem of frequency drift of the hydrogen atomic clock in atomic time calculation, but only takes the frequency drift as a constant obtained by fitting historical data in the calculation process, and does not consider that the frequency drift of the hydrogen atomic clock changes over time in practice, so that the frequency drift cannot be accurately subtracted from the calculated time scale, and the final result still has frequency drift.
Disclosure of Invention
The embodiment of the application provides an atomic time calculation method and device based on hydrogen atomic clock drift prediction, and aims to at least solve the problem of long-time frequency drift in atomic time scale calculation based on a hydrogen atomic clock.
According to one aspect of the application, an atomic time calculation method based on hydrogen atomic clock drift prediction is provided, and comprises the following steps: acquiring historical clock difference data from a database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock; obtaining a clock difference initial value between each clock and a reference clock and a clock difference initial value between each clock and a combined clock according to the clock difference data of each atomic clock; performing time function fitting according to the clock difference data and clock difference values between each clock and the combined clock to obtain a frequency difference drift predicted value changing along with time; and obtaining the clock difference between the combined clock with the atomic clock speed and the clock drift deducted and the reference clock according to the clock difference data, the frequency difference and the frequency difference drift predicted value which changes along with time, and obtaining the atomic time scale.
Further, a clock difference initial value between each clock and the combined clock is obtained according to the clock difference data of each atomic clock, and a clock difference initial value between the combined clock and the reference clock is obtained: and carrying out average weighting on the clock difference data of each atomic clock to obtain an initial value of the clock difference between the combined clock and the reference clock.
Further, performing time function fitting according to the clock difference data and the clock difference value between the combined clock and the reference clock to obtain a frequency difference drift predicted value changing along with time: estimating the clock difference of the next time according to the clock difference, the frequency difference and the frequency drift of the previous time; calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock; calculating clock difference data between the combined clock and the reference clock according to the weight value; and performing time function fitting on the drift values between the clocks and the reference clock to obtain a frequency difference drift predicted value changing along with time.
Further, deducting the frequency offset drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency offset drift predicted value changing along with the time comprises: obtaining a frequency difference drift predicted value of the hydrogen atomic clock at the next moment according to the frequency difference drift predicted value changing along with the time; and deducting the predicted frequency difference drift value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock.
Further, obtaining the predicted value of the frequency difference drift of the hydrogen atomic clock at the next moment according to the predicted value of the frequency difference drift changing along with the time comprises: the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time; and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
According to another aspect of the present application, there is also provided an atomic time computing device based on hydrogen atomic clock drift prediction, including: the acquisition module is used for acquiring historical clock difference data from a database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock; the obtaining module is used for obtaining clock difference initial values between each clock and the combined clock and clock difference initial values between the combined clock and the reference clock according to the clock difference data of each atomic clock; the fitting module is used for performing time function fitting according to the clock difference data and clock difference values between the clocks and the combined clock to obtain a frequency difference drift predicted value changing along with time; and the deduction module is used for deducting the frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale.
Further, the obtaining module is configured to: obtaining clock difference initial values between each clock and the combined clock according to the clock difference data; and carrying out average weighting on the clock difference data of each atomic clock to obtain an initial value of the clock difference between the combined clock and the reference clock.
Further, the fitting module is to: estimating clock difference between each clock and the combined clock at the next moment according to the clock difference, the frequency difference and the frequency drift at the previous moment; calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock; and performing time function fitting on the drift value between each clock and the combined clock according to the clock difference of the next moment of each atomic clock to obtain a frequency difference drift predicted value changing along with time.
Further, the deduction module is to: obtaining a frequency difference drift predicted value of the hydrogen atomic clock at the next moment according to the frequency difference drift predicted value changing along with the time; and deducting the predicted frequency difference drift value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock.
Further, the deduction module is to: the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time; and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
In the embodiment of the application, historical clock difference data are obtained from a database, wherein the clock difference data comprise date and time data and clock difference data between each atomic clock and a reference clock; obtaining a clock difference initial value between each clock and the combined clock and a clock difference initial value between the combined clock and the reference clock according to the clock difference data of each atomic clock; performing time function fitting according to the clock difference data and clock difference values between each clock and the combined clock to obtain a frequency difference drift predicted value changing along with time; and deducting a frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale. The method and the device solve the problem of long-time frequency drift in atomic time scale calculation based on the hydrogen atomic clock, thereby realizing real-time frequency drift updating of the atomic clock in real time in the atomic time scale calculation process, accurately deducting a drift item, and finally obtaining the atomic time scale with both short-term stability and long-term stability.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this application, illustrate embodiments of the application and, together with the description, serve to explain the application and are not intended to limit the application. In the drawings:
FIG. 1 is a flow chart of atomic time calculations based on hydrogen atomic clock drift prediction according to an embodiment of the present application;
fig. 2 is a flowchart of an atomic time calculation method based on hydrogen atomic clock drift prediction according to an embodiment of the application.
Detailed Description
It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict. The present application will be described in detail below with reference to the embodiments with reference to the attached drawings.
It should be noted that the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions and that, although a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than presented herein.
In the present embodiment, a method for calculating atomic time based on hydrogen atomic clock drift prediction is provided, fig. 1 is a flowchart of atomic time calculation based on hydrogen atomic clock drift prediction according to an embodiment of the present application, and steps involved in the flowchart in fig. 1 are described below.
Step S102, historical clock error data are obtained from a database, wherein the clock error data comprise date and time data and clock error data between each atomic clock and a reference clock;
step S104, obtaining a clock difference initial value between each clock and the combined clock and a clock difference initial value between the combined clock and the reference clock according to the clock difference data of each atomic clock; for example, the clock difference data of each atomic clock may be weighted averagely to obtain an initial value of the clock difference between the combined clock and the reference clock.
Step S106, performing time function fitting according to the clock difference data and the clock difference values between each clock and the combined clock to obtain a frequency difference drift predicted value changing along with time; for example, the next clock offset may be estimated from the previous clock offset, the frequency offset, and the frequency drift; calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock; and performing time function fitting on the drift value between each clock and the combined clock according to the clock difference of the next moment of each atomic clock to obtain a frequency difference drift predicted value changing along with time.
And S108, deducting the frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale.
In this step, the predicted value of the frequency difference drift of the hydrogen atomic clock at the next moment can be obtained according to the predicted value of the frequency difference drift changing along with the time; and deducting the predicted frequency difference drift value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock. Optionally, obtaining the predicted frequency difference drift value of the hydrogen atomic clock at the next time according to the predicted frequency difference drift value varying with time includes: the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time; and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
The method solves the problem of long-time frequency drift in atomic time scale calculation based on the hydrogen atomic clock, thereby realizing real-time frequency drift updating of the real-time atomic clock in the atomic time scale calculation process, accurately deducting a drift term, and finally obtaining the atomic time scale with both short-term stability and long-term stability.
This is described below in connection with an alternative embodiment. In this alternative embodiment, an atomic time calculation method based on hydrogen atomic clock drift prediction is provided, and the method is used to record clock difference data of each hydrogen atomic clock in real time, and update a combined clock difference value over time and update a frequency difference value and a clock drift value in real time. The steps in the method are described below with reference to fig. 2.
And S1, extracting clock difference data collected by the time-keeping laboratory from the database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock.
S2: and performing data preprocessing on the extracted data, including abnormal point deletion and data averaging.
S3: obtaining a clock difference initial value between the combined clock and a reference clock and a clock difference initial value between each clock and the reference clock by an average weighting method; respectively obtaining an initial frequency difference value and an initial frequency drift term through a historical data fitting method; the initial weights of the clocks are set to equal weights.
S4: based on an AT1 algorithm, an atom time scale calculation method is improved, and an exponential filtering method based on time function limitation is adopted to obtain an accurate predicted value of frequency drift introduced by a hydrogen atom clock.
S5: and in the atomic time scale calculation process, updating the frequency drift value introduced by the hydrogen atomic clock in real time, and finally obtaining the atomic time scale after the frequency drift is accurately deducted.
The specific implementation process in this optional embodiment includes estimating the next clock offset according to the previous clock offset, the frequency offset, and the clock drift; calculating the weight value of each clock at the next moment according to the expected error of each clock; and updating the frequency difference value at the next moment by an exponential filtering method. Specifically, the exponential filtering method based on the time function limit is used for updating the hydrogen clock drift value at the next moment. The method specifically comprises the following steps: the first difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value; and calculating the obtained frequency drift predicted value according to a time function fitted by historical frequency difference data. And adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the accurate predicted value of the clock drift at the next moment.
The first order difference formula of the frequency difference is as follows:
the exponential filtering formula for calculating the frequency drift value is as follows:
the calculation formula of the accurate frequency difference prediction value is as follows:
wherein, y i (t) represents the clock speed of the ith clock at time t,
and
respectively representing the frequency drift calculation value and the frequency drift prediction value of the ith clock at the time of (t + tau). By minimizing the Allan variance, the corresponding time interval τ is obtained min Is used forExponential frequency average time constant m i And (4) calculating.
M is said i The calculation formula of (2) is as follows:
the atomic time calculation method based on hydrogen atomic clock drift prediction in the optional embodiment can realize real-time frequency drift real-time updating of the atomic clock in the atomic time scale calculation process, so that the drift term is accurately deducted, and finally, the atomic time scale with both short-term stability and long-term stability is obtained.
The atomic time calculation method based on hydrogen atomic clock drift prediction provided by the optional embodiment may include original clock difference data extraction, clock difference data preprocessing, and algorithm calculation and atomic time scale update. Specifically, the atomic time calculation method described in this example is improved on the basis of the existing AT1 algorithm, and the frequency drift of the hydrogen atomic clock is calculated in real time by using an exponential filtering method based on a time function, so that the obtained atomic time scale accurately subtracts the frequency drift error introduced by the hydrogen atomic clock, and higher long-term stability is obtained.
First, clock difference data collected by a time keeping laboratory, including date and time data and clock difference data between each atomic clock and a reference clock, is extracted from a database.
Then, the extracted data is subjected to data preprocessing including deletion of outliers and data averaging.
Obtaining a clock difference initial value between the combined clock and the reference clock by an average weighting method; respectively obtaining an initial frequency difference value and an initial frequency drift term through a historical data fitting method; the initial weights of the clocks are set to equal weights.
Suppose that the clock difference value of the ith clock relative to the combined clock is x i Frequency difference of y i The clock drift is z i The clock difference of the ith clock relative to the reference clock is x ij . And predicting the clock difference value of the next moment (t +1 moment) according to the clock difference and the frequency difference of the previous moment (t moment) and the frequency drift, and obtaining a formula (1).
Further, a clock difference value x of the combined clock relative to the reference clock is obtained r And the calculated value x of the clock difference between each clock and the combined clock i See equations (2) and (3).
x i (t+τ)＝x ij (t+τ)-x r (t+τ) (3)
Wherein, w i The (t + τ) is a weight value corresponding to each atomic clock at the time (t + τ). In the AT1 algorithm, the weight update method relies on the error between the predicted clock difference value and the calculated clock difference value of each clock and the combined clock. Under the condition of the maximum weight value, the weight value of each clock is inversely related to the corresponding expected error. The method specifically comprises the following steps:
calculating the frequency difference value at the (t + tau) moment by an exponential filtering mode
Using clock error data at time (t + tau)And (4) carrying out primary difference and then carrying out exponential filtering to obtain the difference, which is shown in formulas (8) - (9).
Wherein m is i For exponential frequency-averaged time constants, the corresponding time interval τ is obtained by minimizing the Allan variance min For exponential frequency average time constant m i And (4) calculating.
M is i The calculation formula of (2) is as follows:
specifically, the frequency drift calculation value at the time (t + τ) is calculated by an exponential filtering method, and a first difference of the frequency difference data at the time (t + τ) is adopted and then obtained by exponential filtering, as shown in equations (11) - (12). Obtaining a frequency drift predicted value at the (t + tau) moment through least square fitting of historical time period frequency difference data
And the final clock drift predicted value at the (t + tau) moment can be obtained by carrying out equal weight combination on the calculated value of the clock drift and the predicted value.
Wherein, y i (t) represents the clock speed of the ith clock at time t,
and
respectively representing the frequency drift calculation value and the frequency drift prediction value of the ith clock at the moment of (t + tau).
According to the atomic time calculation method based on the hydrogen atomic clock drift prediction, the frequency drift term introduced by the hydrogen atomic clock can be ensured to be updated in real time along with time in the atomic time calculation process, and further the calculated atomic time can be accurately deducted by the frequency drift error.
In this embodiment, an electronic device is provided, comprising a memory in which a computer program is stored and a processor arranged to run the computer program to perform the method in the above embodiments.
The programs described above may be run on a processor or may also be stored in memory (or referred to as computer-readable media), which includes both non-transitory and non-transitory, removable and non-removable media, that implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), Static Random Access Memory (SRAM), Dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), Read Only Memory (ROM), Electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), Digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
These computer programs may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks, and corresponding steps may be implemented by different modules.
Such an apparatus or system is provided in this embodiment. The device is called an atomic time computing device based on hydrogen atomic clock drift prediction, and comprises: the acquisition module is used for acquiring historical clock difference data from a database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock; the obtaining module is used for obtaining a frequency offset drift value between the combined clock and the reference clock according to the clock offset data of each atomic clock; the fitting module is used for performing time function fitting according to the clock difference data and the drift value between the combined clock and the reference clock to obtain a frequency difference drift predicted value changing along with time; and the deduction module is used for deducting the frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale.
The system or the apparatus is used for implementing the functions of the method in the foregoing embodiments, and each module in the system or the apparatus corresponds to each step in the method, which has been described in the method and is not described herein again.
For example, the obtaining module is configured to: and carrying out average weighting on the clock difference data of each atomic clock to obtain an initial value of the clock difference between the combined clock and the reference clock. Optionally, the fitting module is configured to: estimating the next clock difference according to the previous clock difference, the frequency difference and the clock drift; calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock; and performing time function fitting on the drift value between each clock and the combined clock according to the clock difference of the next moment of each atomic clock to obtain a frequency difference drift predicted value changing along with time.
For another example, the deduction module is configured to: obtaining a frequency difference drift predicted value of the hydrogen atomic clock at the next moment according to the frequency difference drift predicted value changing along with the time; and deducting the frequency difference drift predicted value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock. Optionally, the deduction module is configured to: the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time; and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
The atomic time calculation method based on the hydrogen atomic clock drift prediction provided by the embodiment has simple process and easy realization, and can ensure that the finally obtained atomic time scale reaches higher short-term and long-term stability; the clock drift predicted value changing along with the time can be obtained only by adding the frequency drift calculated value obtained by exponential filtering and the frequency drift predicted value obtained by calculating according to the time function fitted by the historical frequency difference data in an equal weight way.
The above are merely examples of the present application and are not intended to limit the present application. Various modifications and changes may occur to those skilled in the art. Any modification, equivalent replacement, improvement or the like made within the spirit and principle of the present application shall be included in the scope of the claims of the present application.
Claims (10)
1. An atomic time calculation method based on hydrogen atomic clock drift prediction is characterized by comprising the following steps:
acquiring historical clock difference data from a database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock;
obtaining a clock difference initial value between the combined clock and the reference clock and a clock difference initial value between each clock and the combined clock according to the clock difference data of each atomic clock;
performing time function fitting according to the clock difference data and clock difference values between each clock and the combined clock to obtain a frequency difference drift predicted value changing along with time;
and deducting a frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale.
2. The method of claim 1, wherein obtaining an initial value of the clock difference between each atomic clock and the combined clock from the clock difference data of each atomic clock and obtaining a value of the clock difference between the combined clock and the reference clock from the clock difference data of each atomic clock comprises:
and carrying out average weighting on the clock difference data of each atomic clock to obtain an initial value of the clock difference between the combined clock and the reference clock.
3. The method of claim 2, wherein a time function fit is performed based on the clock difference data and clock differences between the respective clocks and the combination clock to obtain predicted values of frequency difference drift over time:
estimating the next clock difference according to the previous clock difference, the frequency difference and the frequency difference drift;
calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock;
calculating clock difference data between the combined clock and the reference clock according to the weight value; and performing time function fitting on the drift values between the clocks and the reference clock to obtain a frequency difference drift predicted value changing along with time.
4. The method of claim 3, wherein subtracting a frequency offset drift value introduced by a hydrogen atomic clock of the combined clock from the time-varying frequency offset drift prediction value comprises:
obtaining a frequency difference drift predicted value of the hydrogen atomic clock at the next moment according to the frequency difference drift predicted value changing along with the time;
and deducting the frequency difference drift predicted value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock.
5. The method of claim 4, wherein obtaining the predicted frequency offset drift value of the hydrogen atomic clock at the next moment according to the predicted frequency offset drift value over time comprises:
the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time;
and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
6. An atomic time computing device based on hydrogen atomic clock drift prediction, comprising:
the acquisition module is used for acquiring historical clock difference data from a database, wherein the clock difference data comprises date and time data and clock difference data between each atomic clock and a reference clock;
the obtaining module is used for obtaining clock difference initial values between each clock and the combined clock and clock difference initial values between the combined clock and the reference clock according to the clock difference data of each atomic clock;
the fitting module is used for performing time function fitting according to the clock difference data and clock difference values between the clocks and the combined clock to obtain a frequency difference drift predicted value changing along with time;
and the deduction module is used for deducting the frequency difference drift value introduced by the hydrogen atomic clock of the combined clock according to the frequency difference drift predicted value changing along with the time to obtain an atomic time scale.
7. The apparatus of claim 6, wherein the means for obtaining is configured to:
obtaining clock difference initial values between each clock and the combined clock according to the clock difference data; and carrying out average weighting on the clock difference data of each atomic clock to obtain an initial value of the clock difference between the combined clock and the reference clock.
8. The apparatus of claim 7, wherein the fitting module is configured to:
estimating clock difference between each clock and the combined clock at the next moment according to the clock difference, the frequency difference and the frequency drift at the previous moment; calculating the corresponding weight value of each atomic clock according to the clock difference of the next time of each atomic clock; and performing time function fitting on the drift value between each clock and the combined clock according to the clock difference of the next moment of each atomic clock to obtain a frequency difference drift predicted value changing along with time.
9. The apparatus of claim 8, wherein the subtraction module is configured to:
obtaining a frequency difference drift predicted value of the hydrogen atomic clock at the next moment according to the frequency difference drift predicted value changing along with the time;
and deducting the predicted frequency difference drift value of the hydrogen atomic clock at the next moment from the hydrogen atomic clock of the combined clock.
10. The apparatus of claim 9, wherein the deduction module is configured to:
the primary difference of the frequency difference is used for exponential filtering calculation to obtain a frequency drift calculation value, and an obtained frequency drift prediction value is calculated according to the frequency difference drift prediction value which changes along with time;
and adding the frequency drift calculation value and the frequency drift predicted value in equal weight to obtain the frequency difference drift predicted value of the hydrogen atomic clock at the next moment.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111669122.6A CN114818247A (en) | 2021-12-31 | 2021-12-31 | Atomic time calculation method and device based on hydrogen atomic clock drift prediction |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111669122.6A CN114818247A (en) | 2021-12-31 | 2021-12-31 | Atomic time calculation method and device based on hydrogen atomic clock drift prediction |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114818247A true CN114818247A (en) | 2022-07-29 |
Family
ID=82527612
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111669122.6A Pending CN114818247A (en) | 2021-12-31 | 2021-12-31 | Atomic time calculation method and device based on hydrogen atomic clock drift prediction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114818247A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116819932A (en) * | 2023-08-31 | 2023-09-29 | 湖南时空信安科技有限公司 | Atomic clock steering parameter determining method, device, terminal equipment and storage medium |
-
2021
- 2021-12-31 CN CN202111669122.6A patent/CN114818247A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116819932A (en) * | 2023-08-31 | 2023-09-29 | 湖南时空信安科技有限公司 | Atomic clock steering parameter determining method, device, terminal equipment and storage medium |
| CN116819932B (en) * | 2023-08-31 | 2023-12-01 | 湖南时空信安科技有限公司 | Atomic clock steering parameter determining method, device, terminal equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108563548B (en) | Abnormality detection method and apparatus | |
| CN107480028B (en) | Method and device for acquiring usable residual time of disk | |
| CN110209560B (en) | Data anomaly detection method and detection device | |
| US20120239627A1 (en) | Data storage apparatus and data storage method | |
| EP2613263B1 (en) | Operations management device, operations management method, and program | |
| CN114818247A (en) | Atomic time calculation method and device based on hydrogen atomic clock drift prediction | |
| CN110909306A (en) | Service abnormity detection method and device, electronic equipment and storage equipment | |
| JP7393244B2 (en) | Time series data prediction device and time series data prediction method | |
| CN114708007A (en) | Intelligent decomposition method and system for store sales plan | |
| CN109922212B (en) | Method and device for predicting time-interval telephone traffic ratio | |
| CN111475772B (en) | Capacity assessment method and device | |
| CN114842648B (en) | Tunnel operation state early warning method, device and medium based on traffic flow | |
| CN110020744A (en) | Dynamic prediction method and its system | |
| CN111258854A (en) | Model training method, alarm method based on prediction model and related device | |
| CN113962741B (en) | Coal sales data prediction method, equipment and medium | |
| CN112329197A (en) | Comprehensive atomic time establishing method based on gray model | |
| CN112988536A (en) | Data anomaly detection method, device, equipment and storage medium | |
| CN112182864A (en) | Method for selecting clock error prediction based on drift condition of hydrogen atomic clock | |
| CN117389472A (en) | Cold and hot data migration method and system for transaction data | |
| CN111796995B (en) | Integrated learning-based cyclic serial number usage early warning method and system | |
| CN115396502B (en) | Multi-system message time conversion method and device | |
| CN115019551B (en) | Correction method of redundant bits of field library, training method and device of redundant bit prediction model | |
| CN116483622A (en) | Method, device and equipment for generating abnormal recovery prediction model | |
| CN116959259B (en) | Traffic prediction method based on school student events, electronic equipment and storage medium | |
| CN111143746B (en) | Gray task publishing method, device and equipment |
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 |