CN117014007B - Clock difference driving method and device and terminal equipment - Google Patents
Clock difference driving method and device and terminal equipment Download PDFInfo
- Publication number
- CN117014007B CN117014007B CN202311151079.3A CN202311151079A CN117014007B CN 117014007 B CN117014007 B CN 117014007B CN 202311151079 A CN202311151079 A CN 202311151079A CN 117014007 B CN117014007 B CN 117014007B
- Authority
- CN
- China
- Prior art keywords
- pulse
- clock
- clock difference
- simplified
- difference signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000004590 computer program Methods 0.000 claims description 18
- 229910052792 caesium Inorganic materials 0.000 claims description 10
- TVFDJXOCXUVLDH-UHFFFAOYSA-N caesium atom Chemical compound [Cs] TVFDJXOCXUVLDH-UHFFFAOYSA-N 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 208000002569 Machado-Joseph Disease Diseases 0.000 description 16
- 208000036834 Spinocerebellar ataxia type 3 Diseases 0.000 description 16
- MRMBZHPJVKCOMA-YJFSRANCSA-N biapenem Chemical compound C1N2C=NC=[N+]2CC1SC([C@@H]1C)=C(C([O-])=O)N2[C@H]1[C@@H]([C@H](O)C)C2=O MRMBZHPJVKCOMA-YJFSRANCSA-N 0.000 description 12
- 238000005295 random walk Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 229910052739 hydrogen Inorganic materials 0.000 description 4
- 239000001257 hydrogen Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 108091008152 Machado-Josephin domain proteases Proteins 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G04—HOROLOGY
- G04F—TIME-INTERVAL MEASURING
- G04F5/00—Apparatus for producing preselected time intervals for use as timing standards
- G04F5/14—Apparatus for producing preselected time intervals for use as timing standards using atomic clocks
-
- 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/06—Automatic control of frequency or phase; Synchronisation using a reference signal applied to a frequency- or phase-locked loop
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Manipulation Of Pulses (AREA)
Abstract
The application is applicable to the technical field of pulsar time, and particularly relates to a clock error control method, a clock error control device and terminal equipment. The clock difference driving method comprises the following steps: acquiring observation points of a plurality of pulsars, simplified julian days and timing residual errors, clock difference data of a plurality of atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days; determining a first second pulse clock difference signal based on the observation point number of each pulsar, the simplified julian day, the timing residual error and the first tracing data, wherein the first second pulse clock difference signal is a second pulse clock difference signal when the pulsar is synthesized; determining a second pulse clock difference signal based on clock difference data of each atomic clock, wherein the second pulse clock difference signal is a second pulse clock difference signal when atoms are synthesized; based on the first second pulse clock difference signal, the second pulse clock difference signal is controlled by a phase-locked loop model, and the accuracy of pulsar control atomic clock is improved.
Description
Technical Field
The application belongs to the technical field of pulsar time, and particularly relates to a clock error control method, a clock error control device and terminal equipment.
Background
Pulsar is a compact celestial body, has the characteristics of strong magnetic field and strong electric field, can radiate stable periodic pulse signals, especially the rotation of millisecond pulsar is very stable, and stable and accurate rotation frequency signals can be used as a natural reference frequency source to calibrate an atomic clock.
Pulsar has two natural advantages: firstly, the long-term stability is high, and the defect of atomic time in the aspect can be overcome; secondly, the anti-deception performance is strong, and due to extremely complex pulsar signals and various noise such as cosmic background radiation, deceptions are difficult to generate or tamper the pulsar signals without being detected. Although the atomic time scale has smaller high-frequency noise and higher short-term stability than the pulsar time scale, the pulsar has excellent long-term stability, can be used as a standby time service source of a time service system, and can be used for replacing the atomic clock time service system when the atomic clock time service system is unavailable or can not be calibrated for a long time or needs autonomous operation, and the time service system can maintain reliability by virtue of excellent long-term stability.
Currently, pulsar control atomic clocks can adopt a dual control method, wherein a month control utilizes a predicted pulsar signal to control the atomic clock, the control period is overlong and the reliability is not high, and the accuracy of the pulsar control atomic clock is low.
Therefore, how to improve the accuracy of pulsar-driven atomic clocks is a challenge.
Disclosure of Invention
The embodiment of the application provides a clock difference driving method, a clock difference driving device and terminal equipment, which improve the accuracy of pulsar atomic clock driving.
In a first aspect, embodiments of the present application provide a method of clock-level steering, the method comprising: acquiring observation points, simplified julian days and timing residual errors of a plurality of pulsars, clock difference data of a plurality of atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days; determining a first second pulse clock difference signal based on the observation point number of each pulsar, the simplified julian day, the timing residual error and the first tracing data, wherein the first second pulse clock difference signal is a second pulse clock difference signal when the pulsar is synthesized; determining a second pulse clock difference signal based on clock difference data of each atomic clock, wherein the second pulse clock difference signal is a second pulse clock difference signal when atoms are synthesized; the second pulse clock signal is steered by a phase-locked loop model based on the first second pulse clock signal.
In a possible implementation manner, after the step of steering the second pulse clock difference signal through a phase-locked loop model based on the first pulse second clock difference signal, the method further includes: calibrating a second pulse clock signal of a cesium clock based on a second pulse clock signal of a last second of the first second pulse clock signal and a second pulse clock signal of a last second of the second pulse clock signal after steering.
In one possible implementation manner, the determining the first second pulse clock difference signal based on the number of observation points of each pulsar, the simplified julian day, the timing residual and the first trace data includes: determining second tracing data based on simplified julian days of each pulsar and the first tracing data, wherein the second tracing data comprises timing residual error data from tracing of each pulsar to coordination of universal time; determining a third second pulse clock difference signal based on the simplified julian date of each pulsar and the second tracing data, wherein the third second pulse clock difference signal is a second pulse clock difference signal corresponding to any pulsar; determining a first weight value corresponding to any pulsar based on the observation point number and the timing residual error of each pulsar; the first second pulse clock signal is determined based on the first weight value and the third second pulse clock signal.
In one possible implementation manner, the determining the second tracing data based on the simplified julian day of each pulsar and the first tracing data includes: determining a first simplified julian day and a second simplified julian day based on a target simplified julian day and the first traceable data, the first simplified julian day being a simplified julian day adjacent to the target simplified julian day and less than the target simplified julian day in the first traceable data, the second simplified julian day being a simplified julian day adjacent to the target simplified julian day and greater than the target simplified julian day in the first traceable data, the target simplified julian day being any one of the simplified julian days of each pulsar; determining a first tracing deviation value and a second tracing deviation value based on the first simplified julian day, the second simplified julian day and the first tracing data, wherein the first tracing deviation value is a tracing deviation value corresponding to the first simplified julian day in the first tracing data, and the second tracing deviation value is a tracing deviation value corresponding to the second simplified julian day in the first tracing data; determining a target tracing bias value based on the target simplified julian day, the first simplified julian day, the second simplified julian day, the first tracing bias value and the second tracing bias value, wherein the target tracing bias value is the tracing bias value corresponding to the target simplified julian day; and determining the second tracing data through tracing calculation based on the target tracing deviation value and the first tracing data.
In one possible implementation of the present invention,wherein->For the target trace-source bias value, +.>For the first trace data, +.>In the case of international atomic->And the second tracing data.
In one possible implementation manner, the determining the second pulse clock difference signal based on the clock difference data of each atomic clock includes: based on clock difference data of each atomic clock, determining a fourth second pulse clock difference signal by adopting an atomic clock error model, wherein the fourth second pulse clock difference signal is a second pulse clock difference signal of any atomic clock; determining a second weight value corresponding to any atomic clock based on the fourth second pulse clock difference signal; the second pulse-with-clock-difference signal is determined based on the fourth pulse-with-second-clock-difference signal and the second weight value.
In one possible implementation manner, the determining, based on clock difference data of the plurality of atomic clocks and using an atomic clock error model, a fourth second pulse clock difference signal includes: determining a clock difference, a clock speed, and a parameter value of Zhong Piao of the atomic clock error model based on clock difference data of the plurality of atomic clocks; the fourth second pulse clock difference signal is determined based on clock difference data and clock differences of the plurality of atomic clocks, clock speeds, and parameter values of Zhong Piao.
In one possible implementation, the phase-locked loop model includes loop parameters, the steering the second pulse clock signal by the phase-locked loop model based on the first second pulse clock signal, comprising: determining a pulse per second clock difference frequency signal of the integrated pulsar based on the first pulse per second clock difference signal; determining a target loop parameter based on the loop parameter; the second pulse-second clock signal is steered based on the target loop parameter and the integrated pulsar-time pulse-second clock frequency signal.
In a second aspect, embodiments of the present application provide a clock-difference handling device, the device comprising: the acquisition module is used for acquiring the observation points of the pulsar, the simplified julian days and timing residual errors, the clock difference data of the atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values corresponding to the simplified julian days respectively; the first determining module is used for determining a first second pulse clock difference signal based on the observation point number, the simplified julian day, the timing residual error and the first tracing data of each pulsar, wherein the first second pulse clock difference signal is a second pulse clock difference signal when the pulsar is synthesized; the second determining module is used for determining a second pulse clock difference signal based on clock difference data of each atomic clock, wherein the second pulse clock difference signal is a second pulse clock difference signal when atomic atoms are synthesized; a steering module for steering the second pulse clock signal through a phase-locked loop model based on the first second pulse clock signal.
In a third aspect, an embodiment of the present application provides a terminal device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the method according to the first aspect or any implementation manner of the first aspect when executing the computer program.
In a fourth aspect, embodiments of the present application provide a computer-readable storage medium storing a computer program that, when executed by a processor, implements a method according to the first aspect or any one of the implementations.
In a fifth aspect, embodiments of the present application provide a computer program product, which when run on a terminal device, causes the terminal device to perform the method according to the first aspect or any implementation manner of the first aspect.
Compared with the prior art, the embodiment of the application has the beneficial effects that: acquiring observation points of a plurality of pulsars, simplified julian days and timing residual errors, clock difference data of a plurality of atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days; determining a second pulse clock difference signal when the pulsar is synthesized based on the observation point number of each pulsar, the simplified julian day, the timing residual error and the first tracing data; determining a second pulse clock difference signal when the atoms are synthesized based on clock difference data of each atomic clock; based on the second pulse clock difference signal of the comprehensive pulsar time, the second pulse clock difference signal of the comprehensive atomic time is driven by a phase-locked loop model, compared with a method (daily driving and month driving) for realizing the driving of the pulsar time to the atomic clock by using double driving, the accuracy of the pulsar driving atomic clock is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following description will briefly introduce the drawings that are needed in the embodiments or the description of the prior art, it is obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a method of clock steering according to one embodiment of the present application;
FIG. 2 is a flow chart of one implementation of S120 in a method of clock steering, according to one embodiment of the present application;
FIG. 3 is a flow chart of one implementation of S121 in a method of clock-level steering, as provided by an embodiment of the present application;
FIG. 4 is a flow chart of one implementation of S130 in a method of clock steering, provided in an embodiment of the present application;
FIG. 5 is a flow chart of one implementation of S131 in a clock-level steering method, according to one embodiment of the present application;
FIG. 6 is a flow chart of one implementation of S140 in a method of clock steering, provided in an embodiment of the present application;
FIG. 7 is a block diagram of a clock steering apparatus according to one embodiment of the present application;
Fig. 8 is a schematic structural diagram of a terminal device according to an embodiment of the present application.
Detailed Description
In the following description, for purposes of explanation and not limitation, specific details are set forth, such as particular system configurations, techniques, etc. in order to provide a thorough understanding of the embodiments of the present application. It will be apparent, however, to one skilled in the art that the present application may be practiced in other embodiments that depart from these specific details. In other instances, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present application with unnecessary detail.
It should be understood that the terms "comprises" and/or "comprising," when used in this specification and the appended claims, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It should also be understood that the term "and/or" as used in this specification and the appended claims refers to any and all possible combinations of one or more of the associated listed items, and includes such combinations.
As used in this specification and the appended claims, the term "if" may be interpreted as "when..once" or "in response to positioning" or "in response to detection" depending on the context. Similarly, the phrase "if located" or "if [ a described condition or event ] is detected" may be interpreted in the context of meaning "upon locating" or "in response to locating" or "upon detecting [ a described condition or event ]" or "in response to detecting [ a described condition or event ]".
In addition, in the description of the present application and the appended claims, the terms "first," "second," "third," and the like are used merely to distinguish between descriptions and are not to be construed as indicating or implying relative importance.
Reference in the specification to "one embodiment" or "some embodiments" or the like means that a particular feature, structure, or characteristic described in connection with the embodiment is included in one or more embodiments of the application. Thus, appearances of the phrases "in one embodiment," "in some embodiments," "in other embodiments," and the like in the specification are not necessarily all referring to the same embodiment, but mean "one or more but not all embodiments" unless expressly specified otherwise. The terms "comprising," "including," "having," and variations thereof mean "including but not limited to," unless expressly specified otherwise.
For easy understanding, the technical solutions of the present application will be described in detail below with reference to the accompanying drawings.
Fig. 1 is a flow chart of a method of clock steering according to an embodiment of the present application. As shown in fig. 1, the method at least includes S110 to S140.
S110, obtaining observation points of a plurality of pulsars, simplified julian days and timing residual errors, clock difference data of a plurality of atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days.
In one possible implementation manner, based on the number of observation points and the observation time of each of the 37 pulse stars recorded in the pulse star data set, and a preset number of observation points threshold value and a preset observation time threshold value, a plurality of pulse stars are screened out as frequency sources, and according to ephemeris files and model parameters in the pulse star data set, data such as simplified julian days (Modified Julian Data, MJD), timing residuals, errors of timing residuals and the like of each of the screened plurality of pulse stars are fitted and derived by using Tempo2 software. The number of observation points of each screened pulsar is larger than a preset observation point threshold value, and the observation time is larger than a preset observation time threshold value.
Illustratively, the preset number of observation points threshold is 5000 and the preset observation time threshold is 10 years.
It should be noted that, the preset observation point number threshold value, the preset observation time threshold value and the number of the screened multiple pulsars can be set in a self-defined manner according to actual conditions, which is not limited in this application.
In one possible implementation, clock difference data for a plurality of atomic clocks is obtained. The clock difference data of each atomic clock comprises actual measurement clock difference data of each atomic clock for 5 days/point, the MJD range is 53739-55584, and the data time span is 5 years.
As an example, for the case that data of a few observation points of a part of the acquired clock difference data of each atomic clock are missing, linear interpolation is performed first, noise with the average value of the secondary difference of the data obtained after the linear interpolation as the average value of the amplitude is added, clock difference data corresponding to the observation points with the data missing are determined, and individual missing signals in the clock difference data of each atomic clock are supplemented.
In one possible implementation, earth Time (TT) data, such as TT (BIPM 19), is acquired, and data between MJD53739 and MJD54104 is truncated as first trace data, and since the trace data is 10 days apart, linear interpolation is required for timing residual points of each screened pulsar.
As an example, the first trace data includes, but is not limited to, a plurality of simplified julian days and trace offset values corresponding to each simplified julian day, an adjustment amount when TT (BIPM) traces to a free atom, an adjustment amount when TT (BIPM) traces to an international atom (International Atomic Time, TAI). The tracing deviation value is a difference between a value of pulsar tracing to the earth and a value of tracing to Coordinated world time (Coordinated UniversalTime, UTC), for example, the tracing deviation value includes a difference between a timing residual of pulsar tracing to the earth and a timing residual of tracing to Coordinated world time.
S120, determining a first second pulse clock difference signal based on the observation point number of each pulsar, the simplified julian day, the timing residual error and the first tracing data, wherein the first second pulse clock difference signal is a second pulse clock difference signal when the pulsar is synthesized.
As an example, fig. 2 is a flow chart of one implementation of S120 in a method of clock-level driving, according to an embodiment of the present application. As shown in fig. 2, the method includes at least S121 to S124.
S121, determining second tracing data based on the simplified julian day and the first tracing data of each screened pulsar, wherein the second tracing data comprises timing residual data from tracing each pulsar to coordinating universal time.
As an example, fig. 3 is a flow chart illustrating an implementation of S121 in a method of clock-level driving according to an embodiment of the present application. Referring to fig. 3, the method includes at least S1211 to S1214.
S1211, determining a first simplified julian day and a second simplified julian day based on the target simplified julian day and the first tracing data, the first simplified julian day being a simplified julian day adjacent to and less than the target simplified julian day in the first tracing data, the second simplified julian day being a simplified julian day adjacent to and greater than the target simplified julian day in the first tracing data, the target simplified julian day being any one of the simplified julian days of each pulsar.
As one example, the target simplified julian day is 51910, the first simplified julian day is 51909, and the second simplified julian day is 51919, by looking up a plurality of simplified julian days in the first traceable data.
S1212, determining a first tracing bias value and a second tracing bias value based on the first simplified julian day, the second simplified julian day and the first tracing data, wherein the first tracing bias value is a tracing bias value corresponding to the first simplified julian day in the first tracing data, and the second tracing bias value is a tracing bias value corresponding to the second simplified julian day in the first tracing data.
As an example, based on the first simplified julian day 51909 and the second simplified julian day 51919, the first tracing deviation value and the second tracing deviation value corresponding to the first simplified julian day and the second simplified julian day are searched in the first tracing data, and finally the first tracing deviation value is 25.677 and the second tracing deviation value is 25.682.
For example, the first tracing bias value may be a difference between a timing residual of pulsar tracing to TT (BIPM 19) and a timing residual of tracing to UTC corresponding to the first simplified julian day, and the second tracing bias value may be a difference between a timing residual of pulsar tracing to TT (BIPM 19) and a timing residual of tracing to UTC corresponding to the second simplified julian day.
S1213, determining a target tracing bias value based on the target simplified julian day, the first simplified julian day, the second simplified julian day, the first tracing bias value and the second tracing bias value, wherein the target tracing bias value is the tracing bias value corresponding to the target simplified julian day.
As an example, a linear interpolation process is performed on the target simplified julian day 51910 based on the first simplified julian day 51909, the second simplified julian day 51919, the first tracing bias value 25.677, and the second tracing bias value 25.682 to obtain a target tracing bias value.
Illustratively, the target traceability offset value is obtained by calculating 25.677+ (25.682-25.677) × (51910-51909)/(51919-51909).
For example, when the first tracing bias value is a difference between a timing residual of pulsar tracing to TT (BIPM 19) corresponding to the first simplified julian day and a timing residual of tracing to UTC, the second tracing bias value is a difference between a timing residual of pulsar tracing to TT (BIPM 19) corresponding to the second simplified julian day and a timing residual of tracing to UTC, the target tracing bias value is a difference between a timing residual of pulsar tracing to TT (BIPM 19) corresponding to the target simplified julian day and a timing residual of tracing to UTC.
S1214, determining second tracing data through a tracing calculation formula based on the target tracing deviation value and the first tracing data.
As an example, the traceability calculation is。
Wherein,for the target tracing deviation value, +.>For the first trace data,/o>In the case of an international atom, the atomic number,is the second trace data. And (3) performing simultaneous computation on two computing formulas in the tracing computing formulas to obtain second tracing data.
Illustratively, whenThe represented target tracing bias value is the difference between the timing residual error from pulsar tracing corresponding to the simplified julian day to TT (BIPM 19) and the timing residual error from tracing to UTC, and when TT (BIPM 19) is the timing residual error from pulsar tracing corresponding to the simplified julian day to TT (BIPM 19), the second tracing data obtained by the tracing formula is the timing residual error from pulsar tracing corresponding to the simplified julian day to UTC.
Based on the simplified julian day and the first tracing data of each screened pulsar, calculating by using the tracing calculation formula, and obtaining the timing residual error from tracing each screened pulsar to UTC.
S122, determining a third second pulse clock difference signal based on the simplified julian date of each pulsar and the second tracing data, wherein the third second pulse clock difference signal is a second pulse clock difference signal corresponding to any pulsar.
In one possible implementation, an MJD index table is built with minimum value of 0 and maximum value of the number of points of the timing residual between MJD53739 and MJD 54104. Since 365 days exist a year and 31536000 seconds are left for 365 days, 31536000 values corresponding to 31536000 seconds are required to be linearly interpolated in the second traceable data and uniformly distributed into gaps of the number of timing residual points, and the table is built, namely, the third second pulse clock difference signal is determined.
As an example, three cases need to be considered when performing linear interpolation, which correspond to the head, middle and tail of the existing data, respectively. The head is the interpolation point before all the existing timing residual points, the middle is the interpolation point between all the existing residual points, and the tail is the interpolation point after all the existing residual points.
For the middle part of the data, directly carrying out linear interpolation by using the front timing residual value and the rear timing residual value of the MJD index; and for the head and tail, the first two points and the last two points of the prior timing residual data are respectively used for interpolation. After the timing residual error interpolation of each pulsar is finished, the original timing residual error point can be abandoned, so that a third second pulse signal is obtained.
As an example, for header data, assume that the abscissa of the first two points of the existing timing residual data before interpolation are respectivelyAnd->The value of the respective corresponding timing residuals is +.>And->. The abscissa of the point to be interpolated is +.>Interpolation is +.>. Then->Time->. And similarly, performing linear interpolation calculation by using two points at the tail of the original sequence.
In one possible implementation, before performing linear interpolation, additive noise is generated according to the average value of the second order difference of the timing residual error obtained in S110, and after the linear interpolation is completed, the generated additive noise is added to the second pulse clock difference signal of each pulsar, so as to approximate to the real signal noise characteristic.
S123, determining a first weight value corresponding to any pulsar based on the observation point number and the timing residual error of each pulsar.
In one possible implementation, the screened pulsars are arranged in an ascending order of MJD. The time scale of the second pulse clock difference is set to be 1 year, because the acquired atomic clock difference data starts from MJD53739, the data of each pulsar after MJD54104 needs to be removed first, but the data before MJD53739 is temporarily reserved, after the data after MJD54104 is deleted, among the screened pulsars, the pulsar with the least observation point number for timing residual is selected as a reference pulsar, the MJDs of other pulsars are aligned to the MJD of the reference pulsar, and the MJD of the reference pulsar is recorded as the reference MJD. Based on the reference MJD and the number of observation points of each screened pulsar, the timing residual errors and errors of the rest pulsars except the reference pulsar are downsampled, the timing residual errors and downsampled errors of the rest pulsars are obtained, and then the stability of each pulsar is calculated based on the downsampled timing residual errors and downsampled errors of each pulsar.
As an example, based on the post-pulsar downsampled timing residuals, the downsampled errors, andvariance (also known as Hadamard variance) was used to calculate the stability of each pulsar.
Wherein the stability of each pulsarThe calculation formula of (2) is specifically +.>。
Wherein the timing residual error after each pulsar downsampling and the error after downsampling are equally divided into equal intervalsIs fitted with a cubic polynomial in each subsequence, note +.>For the coefficients of the highest order term, +.>Represented on all subsequences to interact withThe square of the error is inversely proportional to the weighted average.
It should be noted that, based on the reference MJD and the number of observation points of each pulsar, the timing residual error and error of the rest of the pulsars except the reference pulsar are downsampled, and specific implementation steps of the timing residual error after downsampling of the rest of the pulsars and the error after downsampling may refer to implementation methods in the prior art, which are not described in detail in the present application.
Specifically, based on calculated pulsarsThe stability is calculated by the following methodIs a first weight value of (a): />。
Wherein,is->Permallet->Stability. />For the number of pulsars selected, +.>And the first weight value corresponds to each pulsar.
S124, determining a first second pulse clock difference signal based on the first weight value and the third second pulse clock difference signal.
In one possible implementation, the third second pulse clock difference signal corresponding to each first weight value is weighted based on the first weight value of each screened pulse star, so as to obtain a total time residual, namely the first second pulse clock difference signal.
S130, determining a second pulse clock difference signal based on clock difference data of each atomic clock, wherein the second pulse clock difference signal is a second pulse clock difference signal when atomic atoms are synthesized.
As an example, fig. 4 is a flow chart illustrating an implementation of S130 in a method of clock-level driving according to an embodiment of the present application. As shown in fig. 4, the method includes at least S131 to S133.
S131, determining a fourth second pulse clock difference signal based on clock difference data of a plurality of atomic clocks by adopting an atomic clock error model, wherein the fourth second pulse clock difference signal is a second pulse clock difference signal of any atomic clock.
As an example, fig. 5 is a flow chart illustrating an implementation of S131 in a method of clock-level driving according to an embodiment of the present application. As shown in fig. 5, the method includes at least S1311 to S1312.
S1311, determining the clock difference, clock speed, and parameter values of Zhong Piao of the atomic clock error model based on the clock difference data of the plurality of atomic clocks.
As an example, the calculation formula of the atomic clock error model is specifically as follows。
Wherein,，/>and->The initial clock difference, initial clock speed, and initial Zhong Piao, respectively, are deterministic components, with the remainder being noise. />Is the atomic clock error obtained by fitting.
Based on the atomic clock error model, performing quadratic polynomial fitting on the acquired clock difference data of the atomic clocks to complete approximation of clock difference, clock speed and Zhong Piao parameters, and determining the clock difference, clock speed and Zhong Piao parameter values of the atomic clock error model.
S1312 determines a fourth second pulse clock difference signal based on clock difference data and clock differences of the plurality of atomic clocks, clock speeds, and parameter values of Zhong Piao.
In one possible implementation, the MJD of the obtained plurality of atomic clocks is aligned with the MJD of the first second pulse clock difference signal, that is, clock difference data between MJD53739-MJD54104 in clock difference data of the plurality of atomic clocks is intercepted, and linear interpolation is performed on the intercepted clock difference data of the plurality of atomic clocks based on the clock difference, clock speed and parameter values of Zhong Piao of the atomic clock error model, so as to obtain second pulse clock difference signals after interpolation of each atomic clock.
In one possible implementation, noise is simulated and added to the interpolated second pulse clock difference signal for each atomic clock to obtain a fourth second pulse clock difference signal.
Specifically, for a hydrogen clock and a cesium clock, the main noise is frequency white noise and frequency random walk noise, so that the noise can be modeled as a superposition of the frequency white noise and the frequency random walk noise, but the noise and the frequency random walk noise have different proportionality coefficients, the proportion of the frequency white noise of the hydrogen clock is higher, and the proportion of the frequency random walk noise of the cesium clock is higher. And respectively adding proper simulation noise into clock differences of corresponding atomic clocks to generate fourth second pulse clock difference signals of the atomic clocks.
Illustratively, the coefficients of the frequency white noise and the frequency random walk noise of the hydrogen clock are 1×10-25 and 1.1×10-35, respectively, and the coefficients of the frequency white noise and the frequency random walk noise of the cesium clock are 4.8×10-23 and 1.9×10-36, respectively.
The above-mentioned coefficients of the frequency white noise and the frequency random walk noise of the hydrogen clock and the coefficients of the frequency white noise and the frequency random walk noise of the cesiated clock are only examples, and the present application is not limited thereto.
S132, determining a second weight value corresponding to any atomic clock based on the fourth second pulse clock difference signal.
In one possible implementation, after the fourth second pulse clock difference signal is generated, an alan deviation of the fourth second pulse clock difference signal is evaluated, and based on the alan deviation, a second weight value of each atomic clock is calculated by:
。
wherein,is->+.>Stability (i.e. Allan deviation),>for the second weight value corresponding to each atomic clock, < ->Indicating the number of atomic clocks.
And S133, determining a second pulse clock difference signal based on the fourth pulse clock difference signal and the second weight value.
In one possible implementation, because the calculated second weight values are normalized, each corresponding atomic clock may be weighted directly based on the second weight values, generating a second pulse-second clock difference signal.
S140, based on the first second pulse clock difference signal, the second pulse clock difference signal is driven through a phase-locked loop model.
As an example, fig. 6 is a flow chart of one implementation of S140 in a method of clock-level driving according to an embodiment of the present application. As shown in fig. 6, the method includes at least S141 to S143.
S141, determining a second pulse clock difference frequency signal of the integrated pulsar based on the first second pulse clock difference signal.
In one possible implementation, the pulsar signal is used as a steering frequency source, and the steering physical quantity to be provided is frequency, not phase, so that the first second pulse clock difference signal clock difference needs to be subjected to first-order difference to obtain a second pulse clock difference frequency signal when the pulsar is synthesized, and the second pulse clock difference frequency signal is changed from a phase signal to a frequency signal, thereby completing steering of the second pulse clock difference signal in a subsequent phase-locked loop steering algorithm.
S142, determining a target loop parameter based on the loop parameter.
In one possible implementation, loop parameters are included in the phase-locked loop model, loop divergence refers to loss of steering function by the loop, and explosive growth of the output signal, such that the local frequency source cannot be tracked during the steering period, and the steering frequency source cannot be tracked during the steering period. Therefore, the loop parameters of the phase locked loop must be correctly adjusted to converge the loop. Considering that the third-order phase-locked loop can stably track the frequency ramp signal, the performance is superior to that of the second-order phase-locked loop, the control stability of the first second pulse clock difference signal to the second pulse clock difference signal can be ensured, and the system transfer function of the digital phase-locked loop (digital phase locked loop, DPLL) is selected as follows
。
Wherein,is the sampling period +.>Is the phase detector gain,/>Is the voltage controlled oscillator gain, and +.>，/>And->Is a time constant.
The intersection bandwidth of the system transfer function and the error transfer function is。
If and only if the single sideband phase noise crossing of the first second pulse clock signal and the second pulse clock signal is equal toAnd when the control device is used, the optimal loop parameters are obtained, the loop convergence is ensured, and the control effect is optimal.
By way of example only, and not by way of limitation,the preferred value range of (2) is [1 x 10-7, 1 x 10-5 ]]。
S143, based on the target loop parameter and the second pulse clock difference frequency signal when the pulsar is synthesized, the second pulse clock difference signal is driven.
In one possible implementation, when the output signal of the phase-locked loop model is in a first steering period, the steered source, i.e., the second pulse clock difference signal, is output; in the subsequent control period, the frequency of the second pulse clock difference signal is controlled in a feedback way by calculating the difference value between the frequency of the first second pulse clock difference signal and the frequency of the second pulse clock difference signal, so that the uncontrolled phase drift of the second pulse clock difference signal is avoided; meanwhile, the excellent short-term stability of the second pulse clock difference signal is utilized, so that a large amount of high-frequency noise of the first second pulse clock difference signal is avoided, and the second-order difference of the output signal approximates to the second-order difference of the second pulse clock difference signal.
After the phase-locked loop is used for controlling, the finally generated first second pulse clock difference signal and the second pulse clock difference signal inherit the good short stability of the second pulse clock difference signal and the good long stability of the first second pulse clock difference signal.
In one possible implementation, after performing S140, a clock-skew steering method provided herein further includes: the second pulse clock signal of the cesium clock is calibrated based on the second pulse clock signal of the last second of the first second pulse clock signal and the second pulse clock signal of the last second of the second pulse clock signal after steering.
As an example, of the first second pulse clock signal and the second pulse clock signal after the steering, which are obtained with high accuracy and high stability within 1 year based on the above steps, only the last second pulse clock signal can be called time service because of the two input signals of the first second pulse clock signal and the second pulse clock signal after the steering, only the last second pulse clock signal is real-time, and all the rest of the data is historic. Therefore, the second pulse clock difference of the last second in the first second pulse clock difference signal and the second pulse clock difference signal after being steered are selected as the calibration sources of the local cesium clock.
When the phase-locked loop is operated once a day, the method is used for introducing new data of one day, deleting the earliest data of one day in the last group of input data, and then the second pulse clock difference of the last second in the first second pulse clock difference signal and the second pulse clock difference of the last second in the second pulse clock difference signal after being driven are calculated by driving in the same method. At this point, the phase of the local Cesium clock is calibrated with this calculated phase without changing the speed of the Cesium clock and Zhong Piao. The actual time service day error of the time service system can be calculated by adding the clock difference introduced by the local cesium in one day with the errors of the first second pulse clock difference signal and the second pulse clock difference signal after being controlled.
According to the technical scheme, the true timing residual error of each screened pulsar is downsampled according to the reference MJD, is uniformly distributed, the data of TT (BIPM 19) is traced to UTC, the weight value of each pulsar in the comprehensive pulsar is calculated, the MJD of each pulsar is divided according to the second pulse clock difference, and the true single-pulsar second pulse clock difference signal is generated by linearly interpolating according to the true timing residual error at each point, the aligned pulsars generate the second pulse signal of the comprehensive pulsar according to the calculated weight, future data do not need to be predicted, and the accuracy of the second pulse clock difference signal when the comprehensive pulsar is generated is improved; in addition, the clock difference data of each atomic clock are fitted by a quadratic model, the weight of each atomic clock in the comprehensive atomic time is calculated, the noise proportion coefficient of each atomic clock is determined, timing noise is added to each atomic clock according to the noise proportion coefficient, a corresponding second pulse clock difference signal is generated, the second pulse clock difference signal of the comprehensive atomic time is generated according to the calculated weight of each atomic clock and the second pulse clock difference signal of each atomic clock, and the accuracy of the second pulse clock difference signal of the comprehensive atomic time is improved; and finally, the second pulse clock difference signal of the comprehensive atomic clock is controlled according to the generated second pulse clock difference signal of the comprehensive pulsar, so that the accuracy of the second pulse clock difference signal of the comprehensive pulsar is improved.
Fig. 7 is a block diagram of the structure of the clock steering apparatus provided in an embodiment of the present application, showing only those parts relevant to the embodiments of the present application for ease of illustration. Referring to fig. 7, the clock steering device 700 may include an acquisition module 701, a first determination module 702, a second determination module 703, and a steering module 704.
In one implementation, the apparatus 700 may be used to implement the method illustrated in FIG. 1 described above. For example, the acquisition module 701 is used to implement S110, the first determination module 702 is used to implement S120, the second determination module 703 is used to implement S130, and the steering module 704 is used to implement S140.
In another implementation, the apparatus 700 may be used to implement the method illustrated in fig. 2 described above. For example, the first determining module 702 is used to implement S121 to S124.
In yet another implementation, the apparatus 700 may be used to implement the method illustrated in fig. 3 described above. For example, the first determination module 702 is used to implement S1211 to S1214.
In yet another implementation, the apparatus 700 may be used to implement the method illustrated in fig. 4 described above. For example, the second determining module 703 is used to implement S131 to S133.
In yet another implementation, the apparatus 700 may be used to implement the method illustrated in fig. 5 described above. For example, the second determining module 703 is used to implement S1311 to S1312.
In yet another implementation, the apparatus 700 may be used to implement the method illustrated in fig. 6 described above. For example, the steering module 704 is used to implement S141 to S143.
In the embodiment, the number of observation points of a plurality of pulsars, simplified julian days and timing residual errors, clock difference data of a plurality of atomic clocks and first tracing data are obtained, wherein the first tracing data comprise a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days; determining a first second pulse clock difference signal based on the observation point number of each pulsar, the simplified julian day, the timing residual error and the first tracing data, wherein the first second pulse clock difference signal is a second pulse clock difference signal when the pulsar is synthesized; determining a second pulse clock difference signal based on clock difference data of each atomic clock, wherein the second pulse clock difference signal is a second pulse clock difference signal when atoms are synthesized; based on the first second pulse clock difference signal, the second pulse clock difference signal is controlled by a phase-locked loop model, and the accuracy of pulsar control atomic clock is improved.
It should be noted that, because the content of information interaction and execution process between the above devices/units is based on the same concept as the method embodiment of the present application, specific functions and technical effects thereof may be referred to in the method embodiment section, and will not be described herein again.
Fig. 8 is a schematic structural diagram of a terminal device according to an embodiment of the present application. As shown in fig. 8, the terminal device 8 of this embodiment includes: at least one processor 80 (only one shown in fig. 8), a memory 81, and a computer program 82 stored in the memory 81 and executable on the at least one processor 80, the processor 80 implementing the steps in any embodiment of the above-described clock-level-ride-through method when executing the computer program 82.
The terminal device 8 may be a computing device such as a desktop computer, a notebook computer, a palm computer, a cloud server, etc. The terminal device may include, but is not limited to, a processor 80, a memory 81. It will be appreciated by those skilled in the art that fig. 8 is merely an example of the terminal device 8 and is not limiting of the terminal device 8, and may include more or fewer components than shown, or may combine certain components, or different components, such as may also include input-output devices, network access devices, etc.
The processor 80 may be a central processing unit (Central Processing Unit, CPU), the processor 80 may also be other general purpose processors, digital signal processors (Digital Signal Processor, DSP), application specific integrated circuits (Application Specific Integrated Circuit, ASIC), off-the-shelf programmable gate arrays (Field-Programmable Gate Array, FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, or the like. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The memory 81 may in some embodiments be an internal storage unit of the terminal device 8, such as a hard disk or a memory of the terminal device 8. The memory 81 may in other embodiments also be an external storage device of the terminal device 8, such as a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card) or the like, which are provided on the terminal device 8. Further, the memory 81 may also include both an internal storage unit and an external storage device of the terminal device 8. The memory 81 is used for storing an operating system, application programs, boot loader (BootLoader), data, other programs etc., such as program codes of the computer program etc. The memory 81 may also be used to temporarily store data that has been output or is to be output.
It will be apparent to those skilled in the art that, for convenience and brevity of description, only the above-described division of the functional units and modules is illustrated, and in practical application, the above-described functional distribution may be performed by different functional units and modules according to needs, i.e. the internal structure of the apparatus is divided into different functional units or modules to perform all or part of the above-described functions. The functional units and modules in the embodiment may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit, where the integrated units may be implemented in a form of hardware or a form of a software functional unit. In addition, specific names of the functional units and modules are only for convenience of distinguishing from each other, and are not used for limiting the protection scope of the present application. The specific working process of the units and modules in the above system may refer to the corresponding process in the foregoing method embodiment, which is not described herein again.
The embodiment of the application also provides a network device, which comprises: at least one processor, a memory, and a computer program stored in the memory and executable on the at least one processor, which when executed by the processor performs the steps of any of the various method embodiments described above.
Embodiments of the present application also provide a computer readable storage medium storing a computer program which, when executed by a processor, implements steps that may implement the various method embodiments described above.
Embodiments of the present application provide a computer program product which, when run on a mobile terminal, causes the mobile terminal to perform steps that may be performed in the various method embodiments described above.
The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the present application implements all or part of the flow of the method of the above embodiments, and may be implemented by a computer program to instruct related hardware, where the computer program may be stored in a computer readable storage medium, where the computer program, when executed by a processor, may implement the steps of each of the method embodiments described above. Wherein the computer program comprises computer program code which may be in source code form, object code form, executable file or some intermediate form etc. The computer readable medium may include at least: any entity or device capable of carrying computer program code to a photographing device/terminal apparatus, recording medium, computer Memory, read-Only Memory (ROM), random access Memory (Random Access Memory, RAM), electrical carrier signals, telecommunications signals, and software distribution media. Such as a U-disk, removable hard disk, magnetic or optical disk, etc.
In the foregoing embodiments, the descriptions of the embodiments are emphasized, and in part, not described or illustrated in any particular embodiment, reference is made to the related descriptions of other embodiments.
Those of ordinary skill in the art will appreciate that the various illustrative elements and algorithm steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, or combinations of computer software and electronic hardware. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the solution. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present application.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus/network device and method may be implemented in other manners. For example, the apparatus/network device embodiments described above are merely illustrative, e.g., the division of the modules or units is merely a logical functional division, and there may be additional divisions in actual implementation, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed may be an indirect coupling or communication connection via interfaces, devices or units, which may be in electrical, mechanical or other forms.
The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
The above embodiments are only for illustrating the technical solution of the present application, and are not limiting; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present application, and are intended to be included in the scope of the present application.
Claims (7)
1. A method of clock-level ride, the method comprising:
acquiring observation points, simplified julian days and timing residual errors of a plurality of pulsars, clock difference data of a plurality of atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values respectively corresponding to the simplified julian days;
Determining second tracing data based on simplified julian days of each pulsar and the first tracing data, wherein the second tracing data comprises timing residual error data from tracing of each pulsar to coordination of universal time;
determining a third second pulse clock difference signal based on simplified julian days of each pulsar and the second tracing data, wherein the third second pulse clock difference signal is a second pulse clock difference signal corresponding to any pulsar;
determining a first weight value corresponding to any pulsar based on the observation point number and the timing residual error of each pulsar;
determining a first second pulse clock difference signal based on the first weight value and the third second pulse clock difference signal, wherein the first second pulse clock difference signal is a second pulse clock difference signal when pulsars are synthesized;
based on clock difference data of the atomic clocks, determining a fourth second pulse clock difference signal by adopting an atomic clock error model, wherein the fourth second pulse clock difference signal is a second pulse clock difference signal of any atomic clock;
determining a second weight value corresponding to any atomic clock based on the fourth second pulse clock difference signal;
determining a second pulse clock difference signal based on the fourth pulse clock difference signal and the second weight value, wherein the second pulse clock difference signal is a pulse second clock difference signal when atoms are synthesized;
Determining a pulse per second clock difference frequency signal of the integrated pulsar based on the first pulse per second clock difference signal;
determining a target loop parameter based on the loop parameter in the phase-locked loop model;
the second pulse-second clock signal is steered based on the target loop parameter and the integrated pulsar-time pulse-second clock frequency signal.
2. The method of claim 1, wherein said step of steering said second pulse-second clock signal based on said target loop parameter and said integrated pulsar-time pulse-second clock signal further comprises, after said step of steering:
calibrating a second pulse clock signal of a cesium clock based on a second pulse clock signal of a last second of the first second pulse clock signal and a second pulse clock signal of a last second of the second pulse clock signal after steering.
3. The method of claim 1, wherein said determining second trace-out data based on the simplified julian day of each of said pulsars and said first trace-out data comprises:
determining a first simplified julian day and a second simplified julian day based on a target simplified julian day and the first traceable data, the first simplified julian day being a simplified julian day adjacent to the target simplified julian day and less than the target simplified julian day in the first traceable data, the second simplified julian day being a simplified julian day adjacent to the target simplified julian day and greater than the target simplified julian day in the first traceable data, the target simplified julian day being any one of the simplified julian days of each pulsar;
Determining a first tracing deviation value and a second tracing deviation value based on the first simplified julian day, the second simplified julian day and the first tracing data, wherein the first tracing deviation value is a tracing deviation value corresponding to the first simplified julian day in the first tracing data, and the second tracing deviation value is a tracing deviation value corresponding to the second simplified julian day in the first tracing data;
determining a target tracing bias value based on the target simplified julian day, the first simplified julian day, the second simplified julian day, the first tracing bias value and the second tracing bias value, wherein the target tracing bias value is the tracing bias value corresponding to the target simplified julian day;
and determining the second tracing data through tracing calculation based on the target tracing deviation value and the first tracing data.
4. The method of claim 3, wherein the traceability calculation is:
，
wherein,for the target trace-source bias value, +.>For the first trace data, +.>In the case of international atomic->And the second tracing data.
5. The method of claim 1, wherein determining a fourth pulse-per-second clock difference signal using an atomic clock error model based on clock difference data of the plurality of atomic clocks comprises:
Determining a clock difference, a clock speed, and a parameter value of Zhong Piao of the atomic clock error model based on clock difference data of the plurality of atomic clocks;
the fourth second pulse clock difference signal is determined based on clock difference data and clock differences of the plurality of atomic clocks, clock speeds, and parameter values of Zhong Piao.
6. A clock-difference handling device, the device comprising:
the acquisition module is used for acquiring the observation points of the pulsar, the simplified julian days and timing residual errors, the clock difference data of the atomic clocks and first tracing data, wherein the first tracing data comprises a plurality of simplified julian days and tracing deviation values corresponding to the simplified julian days respectively;
the first determining module is used for determining second tracing data based on the simplified julian day of each pulsar and the first tracing data, wherein the second tracing data comprises timing residual error data from the tracing of each pulsar to the coordination of the universal time;
the first determining module is further configured to determine a third second pulse clock difference signal based on the simplified julian day of each pulsar and the second tracing data, where the third second pulse clock difference signal is a second pulse clock difference signal corresponding to any pulsar;
The first determining module is further configured to determine a first weight value corresponding to any pulsar based on the number of observation points and the timing residual error of each pulsar;
the first determining module is further configured to determine a first second pulse clock difference signal based on the first weight value and the third second pulse clock difference signal, where the first second pulse clock difference signal is a second pulse clock difference signal when pulsars are synthesized;
the second determining module is used for determining a fourth second pulse clock difference signal based on clock difference data of the atomic clocks by adopting an atomic clock error model, wherein the fourth second pulse clock difference signal is a second pulse clock difference signal of any atomic clock;
the second determining module is further configured to determine a second weight value corresponding to any one of the atomic clocks based on the fourth second pulse clock difference signal;
the second determining module is further configured to determine a second pulse clock difference signal based on the fourth pulse second clock difference signal and the second weight value, where the second pulse second clock difference signal is a pulse second clock difference signal when the atoms are synthesized;
a steering module for determining a pulse-per-second clock difference frequency signal at the integrated pulsar based on the first pulse-per-second clock difference signal;
The steering module is further used for determining a target loop parameter based on the loop parameter in the phase-locked loop model;
the control module is further configured to control the second pulse-second clock signal based on the target loop parameter and the integrated pulsar-time pulse-second clock signal.
7. A terminal device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, characterized in that the processor implements the method according to any of claims 1 to 5 when executing the computer program.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311151079.3A CN117014007B (en) | 2023-09-07 | 2023-09-07 | Clock difference driving method and device and terminal equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311151079.3A CN117014007B (en) | 2023-09-07 | 2023-09-07 | Clock difference driving method and device and terminal equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117014007A CN117014007A (en) | 2023-11-07 |
| CN117014007B true CN117014007B (en) | 2024-01-05 |
Family
ID=88571129
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311151079.3A Active CN117014007B (en) | 2023-09-07 | 2023-09-07 | Clock difference driving method and device and terminal equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117014007B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2378676C1 (en) * | 2008-05-29 | 2010-01-10 | Общество с ограниченной ответственностью Научно-исследовательский центр "НИКА" | Method of determining world (universal) time from pulsar pulses |
| CN112147874A (en) * | 2020-11-05 | 2020-12-29 | 北京航天发射技术研究所 | Time-frequency reference generation device and method based on satellite time service and CPT atomic clock timekeeping |
| CN112433235A (en) * | 2020-11-19 | 2021-03-02 | 北京卫星导航中心 | Method, system, and medium for determining a time reference |
| CN112636857A (en) * | 2020-11-29 | 2021-04-09 | 中国科学院国家授时中心 | High-precision time-frequency signal dynamic driving method based on residual error correction |
| CN113078901A (en) * | 2021-03-30 | 2021-07-06 | 中国科学院国家授时中心 | Method for improving accuracy of atomic clock based on pulsar control |
| CN113985719A (en) * | 2021-10-25 | 2022-01-28 | 中国科学院国家授时中心 | Sliding window-based pulsar time taming cesium atomic clock method |
| CN116627015A (en) * | 2023-05-06 | 2023-08-22 | 中国科学院国家授时中心 | Pulsar time and atomic time fusion method based on wavelet analysis method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US10802036B2 (en) * | 2018-09-14 | 2020-10-13 | Booz Allen Hamilton Inc. | System and methods for updating a reference time from a decaying rotational period of a pulsar |
| EP4053648A4 (en) * | 2019-10-29 | 2023-01-25 | Innovation Academy for Microsatellites of CAS | Navigation satellite time system and autonomous recovery method therefor |
-
2023
- 2023-09-07 CN CN202311151079.3A patent/CN117014007B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2378676C1 (en) * | 2008-05-29 | 2010-01-10 | Общество с ограниченной ответственностью Научно-исследовательский центр "НИКА" | Method of determining world (universal) time from pulsar pulses |
| CN112147874A (en) * | 2020-11-05 | 2020-12-29 | 北京航天发射技术研究所 | Time-frequency reference generation device and method based on satellite time service and CPT atomic clock timekeeping |
| CN112433235A (en) * | 2020-11-19 | 2021-03-02 | 北京卫星导航中心 | Method, system, and medium for determining a time reference |
| CN112636857A (en) * | 2020-11-29 | 2021-04-09 | 中国科学院国家授时中心 | High-precision time-frequency signal dynamic driving method based on residual error correction |
| CN113078901A (en) * | 2021-03-30 | 2021-07-06 | 中国科学院国家授时中心 | Method for improving accuracy of atomic clock based on pulsar control |
| CN113985719A (en) * | 2021-10-25 | 2022-01-28 | 中国科学院国家授时中心 | Sliding window-based pulsar time taming cesium atomic clock method |
| CN116627015A (en) * | 2023-05-06 | 2023-08-22 | 中国科学院国家授时中心 | Pulsar time and atomic time fusion method based on wavelet analysis method |
Non-Patent Citations (4)
| Title |
|---|
| Analysis of the Precision of Pulsar Time Clock Modeltwo;Cheng-shi Zhao;《Chinese Astronomy and Astrophysics》;第291-302页 * |
| 利用Crab脉冲星X射线观测校准星载原子钟频率;童明雷;《航空学报》;第1-11页 * |
| 脉冲星时间开环驾驭原子钟的计量技术标准化探讨;董江;;计量学报(第02期);第243-249页 * |
| 董江 ; .脉冲星时间开环驾驭原子钟的计量技术标准化探讨.计量学报.2020,(第02期),第243-249页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117014007A (en) | 2023-11-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103226324B (en) | A kind of high-accuracy temporal frequency source taming time and frequency standard in real time | |
| ES2713148T3 (en) | Method and apparatus for processing radionavigation signals for atmospheric monitoring | |
| CN112540388B (en) | Satellite communication module and uplink signal Doppler compensation method thereof | |
| Weiss et al. | AT2, a new time scale algorithm: AT1 plus frequency variance | |
| CN112198557B (en) | Data correction method, device, terminal equipment and storage medium | |
| CN113311694A (en) | Method and device for jointly taming rubidium clock by Beidou satellite common vision and unidirectional time service | |
| Cao et al. | Assessment and comparison of satellite clock offset between BeiDou-3 and other GNSSs | |
| CN111948686A (en) | Time synchronization method and device | |
| CN114002712B (en) | GNSS real-time satellite clock error determination system and method based on data adaptive selection | |
| CN110031873B (en) | GNSS multipath signal simulation method and GNSS multipath signal simulator | |
| CN117014007B (en) | Clock difference driving method and device and terminal equipment | |
| Cano et al. | Improving orbital uncertainty realism through covariance determination in geo | |
| Píriz et al. | PulChron: A pulsar time scale demonstration for PNT systems | |
| CN115903444B (en) | Local time generation method and device for establishing distributed constellation time reference | |
| WO2020231641A1 (en) | Methods and enhancing non-global navigation satellite system location and timing pseudorange positioning calculations and systems thereof | |
| US11824635B1 (en) | Hardware module for determining a clock value based on multiple timing references | |
| Lawrence et al. | Navigation sensor, filter, and failure mode simulation results using the distributed Kalman filter simulator (DKFSIM) | |
| US11616587B1 (en) | Remote clock synchronization using network communication and satellite signals | |
| Innerkofler et al. | GNSS radio occultation excess phase processing for climate applications including uncertainty estimation | |
| CN113703017B (en) | Satellite antenna phase center deviation calculation method and device | |
| CN112286039A (en) | Clock calibration method and device and readable storage medium | |
| CN112363184A (en) | Satellite signal simulation system and method | |
| Doberstein et al. | Introduction to the global positioning system | |
| CA2704385A1 (en) | Method and device for integrating a rf gps emulator in a simulation environment | |
| Herring et al. | Gamit |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |