Digital phase-locked loop atomic clock control method equivalent to two-state variable Kalman filter
Technical Field
The invention relates to the field of time frequency and signal processing, in particular to a digital phase-locked loop atomic clock control method which is equivalent to a two-state variable Kalman filter.
Background
Atomic clock steering technology plays an important role in time-keeping laboratories and satellite navigation systems. There are two main purposes of handling atomic clocks: firstly, the time of the driven atomic clock is synchronized with the time of the reference atomic clock, and the time deviation between the driven atomic clock and the reference atomic clock is reduced as much as possible; and secondly, the long-term stability of the driven atomic clock is improved.
Typical atomic clock manipulation methods include two broad categories: an open loop steering method and a closed loop steering method. The core of the open loop driving method is to design a reasonable prediction algorithm. The closed loop control method mainly comprises the following steps: a linear quadratic (LGQ) control method, a switch (Bang-Bang) control method, and a Digital Phase Locked Loop (DPLL) method. However, the parameters of these closed-loop control methods are not easily selected, and generally, for each specific application, a large number of parameters are compared through a large number of simulation experiments to select an optimal set of parameters. In addition, for the DPLL method, if the parameters are not properly selected, not only the handling performance cannot be guaranteed, but also the control system is unstable.
Disclosure of Invention
In order to solve the problems in the prior art, the invention aims to provide a digital phase-locked loop atomic clock control method which is equivalent to a two-state variable Kalman filter.
The principle of the invention is as follows: the hydrogen clock is steered with the cesium clock through a DPLL equivalent to a two-state variable Kalman filter. The invention derives the closed-loop system transfer function of the DPLL, and proves that the DPLL is a second-order class-2 DPLL. The present invention derives a Kalman gain determined approximately by the ratio of the process noise variance and the observed noise variance. Thus, as can be seen from the closed loop transfer function of the DPLL, the 3 parameters of the ratio, as well as the loop gain and the sampling interval, determine the performance of the DPLL. Furthermore, the invention provides a simple and effective method for selecting parameters. Through the equivalent transformation, an Analog Phase Locked Loop (APLL) equivalent to the DPLL is obtained. In the frequency domain, by adjusting parameters, the intersection frequency of the system transfer function and the error transfer function of the equivalent APLL is equal to the intersection frequency of SSB phase noise curves of a hydrogen clock and a cesium clock, and finally the optimal parameters are obtained. Briefly, a reference atomic clock (or timescale) is used to steer the atomic clock (or timescale) such that the steered atomic clock (or timescale) remains synchronized with the reference atomic clock (or timescale).
The technical scheme adopted by the invention is as follows:
a digital phase-locked loop atomic clock control method equivalent to a two-state variable Kalman filter adopts the following steps:
s.1, deriving a transfer function of a DPLL equivalent to a two-state variable Kalman filter;
s.1.1 use DPLL and cesium bell to drive hydrogen bell. Fig. 1 depicts the principle structure of the DPLL, summarized as: acquiring the time deviation of the cesium clock and the hydrogen clock after control; filtering the time deviation by using a loop filter to filter high-frequency noise; and then adjusting the time deviation of the hydrogen clock by using the filtered time deviation to obtain the time deviation of the hydrogen clock after driving. In fig. 1: cs represents the output signal of the cesium clock, symbol Hm represents the output signal of the hydrogen clock, symbol Steered Hm represents the output signal of the hydrogen clock being driven, and LF represents the loop filter.
In the Z domain, the transfer function of the loop filter LF of the DPLL is denoted g (Z); the output of the DPLL is then:
it is clear that,anda system transfer function and an error transfer function of the DPLL respectively; cs represents the output signal of the cesium clock, symbol Hm represents the output signal of the hydrogen clock, and symbol Steered Hm represents the output signal of the hydrogen clock being driven.
By z ═ ej2πf·T(where f represents frequency and T represents sampling period) to equate the DPLL to a corresponding APLL; thus, equation (2) is expressed in terms of SSB phase noise, written as:
wherein, symbol L represents SSB phase noise, symbol G represents transfer function; subscript Cs represents the output signal of the cesium clock, symbol Hm represents the output signal of the hydrogen clock, and symbol Steered Hm represents the output signal of the hydrogen clock being driven;
s.1.2, deriving expressions of a system transfer function and an error transfer function of the DPLL equivalent to the two-state variable Kalman filter according to the system transfer function and the error transfer function of the step S.1.1;
finally, the system transfer function of the DPLL is expressed as:
wherein, K0Is the loop gain. Ks11And Ks21Respectively, the steady state gains of the Kalman filters; t is the sampling interval.
The error transfer function of the DPLL is expressed as:
to this endThe system transfer function and the error transfer function of the DPLL equivalent to a two-state variable Kalman filter are obtained. Represented by formula (22) and formula (23): parameter (Ks)11,Ks21,T,K0) Determines the performance of the DPLL.
S.2, deducing Ks in the system transfer function and the error transfer function according to the system transfer function and the error transfer function of the DPLL obtained in the step S.111And Ks21By process noise variance Q22Observation noise variance R, and sampling interval T.
When the Kalman filter enters a steady state, obtaining
And
from equations (30) and (35), the Kalman gain is given by Q22And R.
Let Q22=1s2Equation (30) and equation (35) are substituted for equation (22) and equation (23) to obtain a system transfer function and an error transfer function, which are respectively expressed as:
and
as can be seen from equations (36) and (37): when Q is22=1s2When, the parameters (R, T, K)0) Determines the performance of the DPLL. Substituting equations (36) and (37) for equation (1) yields the relation between DPLL input and output:
s.3, determining parameters (R, T, K) according to the system transfer function and the error transfer function obtained in the step S.20) Optimizing the frequency stability of the DPLL output specifically comprises the following steps:
first, the phase noise of the cesium clock and the hydrogen clock is expressed as:
and
wherein the symbol L represents SSB phase noise, f0Is the carrier frequency, hiFor noise figure, f is the sideband frequency, i is used to indicate the noise type.
By solving equation LCs(f)-LHm(f) The frequency at the intersection of the cesium clock phase noise curve and the hydrogen clock phase noise curve can be found as f' 0.
Secondly, changing z to ej2πf·TThe system transfer function and the error transfer function of the equivalent APLL are obtained by substituting equations (36) and (37), which are respectively expressed as:
and
from equations (41) and (42), the amplitude-frequency responses 20log | h (f) and 20log | he (f) of the system transfer function and the error transfer function are obtained, and the frequency at the intersection of these two curves is obtained and is denoted as f ".
Third, for a given T, K is given first0Set a value (for convenience of analysis, it can be set as an integer of 1, 2, etc. first, and then gradually set as other values) where (T, K)0) Under the condition, the value of f ' is changed by changing the value of R, so that f ' is equal to f '. At this time, the frequencies of the two intersections are equal. In this case, given T, the corresponding R can be determined. Thus, a set of parameters, denoted (R', K), is obtained that makes the frequencies of the two intersections equal0'). When the frequencies of the two intersection points are equal, the DPLL filters short-term noise of the cesium clock and long-term noise of the hydrogen clock to the maximum extent, so that the phase noise of the output of the DPLL is optimal, and the frequency stability of the output of the DPLL is optimal at the moment.
Fourthly, for a given T, K is changed0Using the same method, a number of sets of optimum parameters (R', K) can be obtained0'). These parameters (R', K)0') are all optimal parameters, all ensure that the frequency stability of the DPLL output is optimal.
Specifically, in step s.1.2 of the present invention, the derivation process of equations (22) and (23) is as follows:
for a two-state variable system, the state equation is expressed as:
wherein x iskAnd ykIs two state variables, T is the sampling interval, ukIs process noise;
the observation equation is expressed as:
zk=xk+wk(4)
wherein z iskAs an observed quantity, wkTo observe noise;
the two equations, the state equation and the observation equation, are expressed in the form of a matrix as:
wherein s isk=[xkyk]T;Jk=[0uk]T;H=[10];
The variances of the process noise and the observation noise are respectively: r ═ E [ wk 2],Wherein symbol E represents a mathematical expectation;
the Kalman filter is described by the following 5 equations:
Pk,k-1=φPk-1,k-1φT+Q (7)
Kk=Pk,k-1·HT(H·Pk,k-1·HT+R)-1(8)
Pk,k=(I-Kk·H)·Pk,k-1(10)
wherein,is an estimate of the state variable at time k, Pk,kIs to estimate the error matrix of the received signal,is a single step prediction value, P, of the state variable at time k-1 to time kk,k-1Is a single step prediction error matrix, KkIs a Kalman gain matrix;
it can be demonstrated that the system defined by equation (5) is fully observable, so Pk,k-1、Pk,kAnd KkAll converge to Pk,k、Pk,k-1And KkRespectively, are recorded as: ps, Ps-And Ks;
from equations (6) and (9), when the Kalman filter enters the steady state, there are:
wherein,is an estimate of the state variable x at time k,is an estimate of the state variable y at time k, KsijAn element representing the ith row and the jth column in the Ks matrix;
defining:
by substituting formula (12) for formula (11), we obtain:
equation (13) is expressed in the Z domain as:
wherein X, Y and V are respectivelyvkZ-transformation of (a);
obtained from formula (14):
from equations (14) and (15), equation (12) is represented in the Z domain as:
wherein Z represents ZkZ-transformation of (a);
from formula (16), we obtain:
defining:
from the formulae (15), (17) and (18) to give
Obviously, z iskIs an input to a Kalman filter, andis the output of the Kalman filter; equation (19) thus shows the relationship between the steady state Kalman filter input and output.
The system transfer function is then expressed as:
from equation (20), the transfer function of the Kalman filter is the same as the closed loop system transfer function of the second-order class 2 DPLL. Thus, the Kalman filter is equivalent to a second order class 2 DPLL.
Due to Ks111, then equation (20) is written approximately as:
for the DPLL to work properly, a delay z is introduced-1The system transfer function of the DPLL is then expressed as:
wherein, K0Is the loop gain;
the error transfer function of the DPLL is expressed as:
thus, a system transfer function and an error transfer function equivalent to the DPLL of the two-state variable Kalman filter are obtained.
In s.2 of the present invention, the derivation process of equations (30) and (35) is as follows:
when the Kalman filter enters steady state, equation (7) is written as:
equation (8) is written as:
equation (10) is written as:
in equations (25) to (27), subscript ij denotes an element corresponding to ith row and jth column of the matrix;
from formula (25), formula (26), formula (2), formula (27), formula 3, we obtain:
from formula (28) and formula (26), formula 2, we obtain:
due to the fact thatIs obtained from formula (29)
From formulae (25) to (27) to give
From formula (31) to yield
Due to Ks21·T<<Ks11From formula (32), we obtain:
by substituting formula (33) for formula (26), the following can be obtained:
from equation (34), we obtain:
the present invention can be applied to: using one atomic clock (such as cesium clock) to drive another atomic clock (such as hydrogen clock); using an atomic clock to drive a Numerical Control Oscillator (NCO); using a paper time scale (for example: TA (k)) to drive an atomic clock; steering one paper time scale (e.g., UTC) over another paper time scale (e.g., TA (k)); an atomic clock (e.g., UTC (NTSC)) is used to drive a paper time scale (e.g., BDT).
In each of the above applications, there is one reference atomic clock (or time scale), and one driven atomic clock (or time scale). In these applications, the reference atomic clocks (or time scales) are: an atomic clock (such as cesium clock), an atomic clock, a paper time scale, and an atomic clock; the atomic clocks (or time scales) driven are respectively: another atomic clock (e.g. hydrogen clock), a numerically controlled oscillator, an atomic clock, another paper time scale, and a paper time scale.
The principles and procedures of the present invention are the same in all of these applications. Technical solution the present invention is described in the context of using cesium bell to drive hydrogen bell. Indeed, the present invention can be considered to be applicable in the context of other applications as described above, provided that the cesium clock is replaced with other reference atomic clocks (or time scales), and the hydrogen clock is replaced with other driven atomic clocks (or time scales).
The present invention uses a cesium clock to steer the hydrogen clock through a DPLL equivalent to a two-state variable Kalman filter. The invention derives the closed-loop system transfer function of the DPLL, and proves that the DPLL is a second-order class-2 DPLL. The present invention derives a Kalman gain determined approximately by the ratio of the process noise variance and the observed noise variance. Thus, as can be seen from the closed loop transfer function of the DPLL, the 3 parameters of the ratio, as well as the loop gain and the sampling interval, determine the performance of the DPLL. Furthermore, the invention provides a simple and effective method for selecting parameters. Through the equivalent transformation, an Analog Phase Locked Loop (APLL) equivalent to the DPLL is obtained. In the frequency domain, by adjusting parameters, the intersection frequency of the system transfer function and the error transfer function of the equivalent APLL is equal to the intersection frequency of SSB phase noise curves of a hydrogen clock and a cesium clock, and finally the optimal parameters are obtained.
The invention has the beneficial effects that:
1. the cesium clock can be effectively used for controlling the hydrogen clock, and the control performance (control error and frequency stability) is very good;
2. the parameter selection method of the invention is relatively simple.
Drawings
FIG. 1 illustrates the implementation principle of the present invention;
FIG. 2 illustrates the algorithm principle of the present invention;
FIG. 3 illustrates SSB phase noise curves for a hydrogen clock and a cesium clock
FIG. 4 illustrates the magnitude-frequency response of the system transfer function and the error transfer function of the DPLL;
FIG. 5 illustrates SSB phase noise curves for a hydrogen clock, a cesium clock, and a steered hydrogen clock;
FIG. 6 illustrates the time difference between a hydrogen clock, a cesium clock, and a hydrogen clock being driven when there is no drift in the hydrogen clock;
FIG. 7 illustrates the Allan bias of the hydrogen clock, cesium clock, and hydrogen clock being driven when there is no drift in the hydrogen clock;
FIG. 8 illustrates the Allan deviation of hydrogen bell, cesium bell, and hydrogen bell under control when there is a drift in hydrogen bell;
FIG. 9 illustrates the ride error when the hydrogen clock has no drift;
Detailed Description
The following description of the embodiments of the present invention will be made with reference to the accompanying drawings.
The invention relates to a digital phase-locked loop atomic clock control method equivalent to a two-state variable Kalman filter, which comprises the following steps.
S.1, deriving a transfer function of a DPLL equivalent to a two-state variable Kalman filter;
s.1.1 a schematic diagram of the implementation of the present invention is shown in fig. 1. The principle of using DPLL with cesium clock to drive the hydrogen clock is depicted in fig. 1. In the figure, symbol Cs represents the output signal of the cesium clock, symbol Hm represents the output signal of the hydrogen clock, symbol steeedhm represents the output signal of the hydrogen clock, and LF represents the loop filter.
The performance of the DPLL is analyzed in the Z-domain, and the transfer function of the loop filter LF of the DPLL is denoted as g (Z). The output of the DPLL is then:
it is clear that,andthe system transfer function and the error transfer function of the DPLL, respectively.
By z ═ ej2πf·T(where f represents frequency and T represents sampling period) may be equated to a corresponding APLL. Thus, equation (2) is expressed in terms of SSB phase noise, written as:
wherein, symbol L represents SSB phase noise, symbol G represents transfer function; the symbol Cs in the subscript represents the output signal of the cesium clock, the symbol Hm represents the output signal of the hydrogen clock, and the symbol steeedhm represents the output signal of the hydrogen clock;
s.1.2 derives expressions of the system transfer function and the error transfer function of the DPLL equivalent to the two-state variable Kalman filter from the system transfer function and the error transfer function of step s.1.1.
For a two-state variable system, the state equation is expressed as:
wherein x iskAnd ykIs two state variables, T is the sampling interval, ukIs process noise.
The observation equation is expressed as:
zk=xk+wk(4)
wherein z iskAs an observed quantity, wkTo observe the noise.
These two equations are expressed in matrix form as:
wherein s isk=[xkyk]T;Jk=[0uk]T;H=[10]。
The variances of the process noise and the observation noise are respectively: r ═ E [ wk 2],Where the symbol E represents the mathematical expectation.
The Kalman filter can be described with the following 5 steps:
Pk,k-1=φPk-1,k-1φT+Q (7)
Kk=Pk,k-1·HT(H·Pk,k-1·HT+R)-1(8)
Pk,k=(I-Kk·H)·Pk,k-1(10)
wherein,is an estimate of the state variable at time k, Pk,kIs to estimate the error matrix of the received signal,is a single step prediction value, P, of the state variable at time k-1 to time kk,k-1Is a single step prediction error matrix, KkIs a Kalman gain matrix.
It can be demonstrated that the system defined by equation (5) is fully observable, hence Pk,k-1、Pk,kAnd KkAre all converged. Handle Pk,k、Pk,k-1And KkRespectively, are recorded as: ps, Ps-And Ks.
From equations (6) and (9), when the Kalman filter enters the steady state, there are:
wherein,is an estimate of the state variable x at time k,is an estimate of the state variable y at time k, KsijRepresenting the element in the ith row and the jth column of the Ks matrix.
Defining:
by substituting formula (12) for formula (11), we obtain:
equation (13) is expressed in the Z domain as:
wherein X, Y and V are respectivelyvkZ transformation of (1).
This can be obtained from equation (14):
from equations (14) and (15), equation (12) is represented in the Z domain as:
wherein Z represents ZkZ transformation of (1).
From formula (16), we obtain:
defining:
from the formulae (15), (17) and (18) to give
Obviously, z iskIs an input to a Kalman filter, andis the output of the Kalman filter. Equation (19) thus shows the relationship between the steady state Kalman filter input and output.
The system transfer function is then expressed as:
from equation (20), the transfer function of the Kalman filter is the same as the closed loop system transfer function of the second-order class 2 DPLL. Thus, the Kalman filter may be equivalent to a second order class 2 DPLL.
Through experiments, Ks is found111, then equation (20) can be written approximately as:
for the DPLL to work properly, a delay z is introduced-1The system transfer function of the DPLL is then expressed as:
wherein, K0Is the loop gain.
The error transfer function of the DPLL is expressed as:
thus, a system transfer function and an error transfer function equivalent to the DPLL of the two-state variable Kalman filter are obtained. As can be seen from equations (22) and (23): parameter (Ks)11,Ks21,T,K0) Determines the performance of the DPLL.
Fig. 2 depicts the algorithmic principles of the DPLL. As can be seen from FIG. 2, the amount of control over the hydrogen clock during each sampling interval is represented in the Z domain asThus, in the time domain, the relationship between the time difference of the hydrogen clock and the time difference of the hydrogen clock is expressed as:
in the formula (25), err (j) is the time difference cs (j) -Hm between the cesium clock and the hydrogen clock at the jth timesteered(j)。
S.2, deducing Ks in the system transfer function and the error transfer function according to the system transfer function and the error transfer function of the DPLL obtained in the step S.111And Ks21Is composed of Q22R, T.
In fact, when the Kalman filter is operated and enters a steady state, the Ks can be obtained11And Ks21The true value of (d). The idea of this stepMeaning that the Kalman filter may not be run, pass Q22R, T to calculate Ks11And Ks21An approximation of (d). Simulation experiments show that: the error of the approximate value from the true value is small. Thus, the Ks may be substituted for the true value by an approximate value11And Ks21The expression of the approximate value of (A) is substituted for the expressions (23) and (24) to obtain the formula (Q)22,R,T,K0) A system transfer function and an error transfer function for the expression. Thus, the parameter (Q)22,R,T,K0) Determines the performance of the DPLL.
When the Kalman filter enters steady state, equation (7) can be written as:
equation (8) can be written as:
equation (10) can be written as:
in equations (25) to (27), subscript ij denotes an element corresponding to ith row and jth column of the matrix.
From formula (25), formula (26), formula (2), formula (27), formula 3, we obtain:
from formula (28) and formula (26), formula 2, we obtain:
due to the fact thatIs obtained from formula (29)
From formulae (25) to (27) to give
From formula (31) to yield
Due to Ks21·T<<Ks11From formula (32), we obtain:
by substituting formula (33) for formula (26), the following can be obtained:
from equation (34), we obtain:
from equations (30) and (35), it can be seen that: the Kalman gain is approximated by Q22And the proportional relation between RAnd (6) determining.
Simulation experiments show that the error of the approximation is small. For example: when Q is22=1s2,R=6×1020s2Then, the approximate value Ks is calculated from the formula (30) and the formula (35)11≈5.42×10-4And Ks21≈4.08×10-11(ii) a And the Kalman filter is operated to obtain the true value Ks11=5.32×10-4And Ks21=3.95×10-11(ii) a It is clear that the error between the approximated and true values is less than 3%. Then, the calculated approximate values of the equations (30) and (35) may be used to approximate instead of the true values.
Let Q22=1s2Equation (30) and equation (35) are substituted for equation (22) and equation (23) to obtain a system transfer function and an error transfer function, which are respectively expressed as:
and
as can be seen from equations (36) and (37): when Q is22=1s2When, the parameters (R, T, K)0) Determines the performance of the DPLL. Substituting equations (36) and (37) for equation (1) yields the relation between DPLL input and output:
s.3, determining parameters (R, T, K) according to the system transfer function and the error transfer function obtained in the step S.20) And optimizing the frequency stability of the DPLL output.
The frequency stability can be represented in the time domain by the Allan variance or bySSB phase noise is represented in the frequency domain. The parameters are easier to determine in the frequency domain than in the time domain. By transforming z to ej2πf·TThe DPLL is equivalent to an equivalent APLL. The performance of the equivalent APLL in the frequency domain can then be analyzed.
In the frequency domain, the parameters (R, T, K) are selected0) The purpose of the method is mainly to optimize the SSB phase noise performance of the equivalent APLL output signal, and the method comprises the following steps: by selecting the parameters (R, T, K)0) The frequency at the intersection of the system transfer function and the error transfer function is made equal to the frequency at the intersection of the phase noise curve of the cesium clock and the hydrogen clock SSB.
First, the phase noise of the cesium clock and the hydrogen clock is expressed as:
and
wherein f is0Is the carrier frequency, hiFor noise figure, f is the sideband frequency, i is used to indicate the noise type.
By solving equation LCs(f)-LHm(f) The frequency at the intersection of the cesium clock phase noise curve and the hydrogen clock phase noise curve can be found as f' 0.
Secondly, changing z to ej2πf·TThe system transfer function and the error transfer function of the equivalent APLL are obtained by substituting equations (36) and (37), which are respectively expressed as:
and
from equations (41) and (42), the amplitude-frequency responses 20log | h (f) and 20log | he (f) of the system transfer function and the error transfer function can be obtained, and the frequency at the intersection of these two curves can be obtained and is denoted as f ".
Third, for a given T, K is given first0Set a value at which (T, K)0) Under the condition, the value of f ' is changed by changing the value of R, so that f ' is equal to f '. At this time, the frequencies of the two intersections are equal. In this case, given T, the corresponding R can be determined. Thus, a set of parameters, denoted (R', K), is obtained that makes the frequencies of the two intersections equal0'). When the frequencies of the two intersection points are equal, the DPLL filters short-term noise of the cesium clock and long-term noise of the hydrogen clock to the maximum extent, so that the phase noise of the output of the DPLL is optimized, and therefore, the frequency stability of the output of the DPLL can be considered to be optimized at this time.
Fourthly, for a given T, K is changed0Using the same method, a number of sets of optimum parameters (R', K) can be obtained0'). These parameters (R', K)0') are all optimal parameters, all ensure that the frequency stability of the DPLL output is optimal, but there is still a subtle difference in the actual handling performance. Through simulation experiments, different groups (R', K) can be analyzed0') ride performance under conditions.
And generating a hydrogen clock and a cesium clock through simulation. The hydrogen clock and the cesium clock mainly contain frequency white noise and frequency random walk noise. The noise coefficient of the hydrogen clock is h(Hm)0=6×10-26And h(Hm)-2=2×10-36Noise figure of cesium clock is h(Cs)0=3×10-23And h(Cs)-2=3.33×10-38. The sequence of hydrogen clock and cesium clock has 10000 samples each, with a sample interval T of 3600 s.
FIG. 2 shows a hydrogen clock and a cesium clock at a sampling frequency f0When it is 10MHzThe theoretical SSB phase noise curve of (a). It can be seen from fig. 2 that the frequency at the intersection is about 10-6.6Hz. Firstly, selecting K012, it was found that when R6 × 1020At the intersection of the amplitude-frequency response 20log | H (f) and 20log | He (f) curves, the frequency is approximately equal to 10-6.6Hz. Thus, a set of parameters (R', K) is obtained0’)=(6×1020s212) by the same method, other parameters similar to (3 × 10) were obtained22s2,30),(5×1017s2And 2) and the like. FIG. 3 shows that (R', K)0’)=(6×1020s2And 12) the amplitude-frequency response of the system transfer function and the error transfer function.
From the formula (3), the formula (41) and the formula (42), the parameters (R', K) can be calculated0') and T3600 s, the theoretical SSB phase noise L of the equivalent APLL outputHm_steered(f) In that respect Returning to the Z domain and the system transfer function shown in equation (36). In case 1, when (R', K)0’)=(5×1017s22), the pole is 0.9968 +/-1 × 10-8i; in case 2, when (R', K)0’)=(6×1020s212), poles are 0.9997 and 0.9938; in case 3, when (R', K)0’)=(3×1022s230), the poles are 0.9999 and 0.9940. In these 3 cases, the system is stable since the poles are all within the unit circle. In case 1, the system is an under-damped control system. In cases 2 and 3, the system is an over-damped control system, but the system is more sensitive in case 2. To sum up, the parameters (R', K) of case 2 were selected0’)=(6×1020s2,12)。
FIG. 4 shows the current parameters (R', K)0’)=(6×1020s2And 12) hydrogen clock, cesium clock, and post-drive hydrogen clock theoretical SSB phase noise. As can be seen from fig. 4, the SSB phase noise performance of the post-drive hydrogen clock is very good.
By the formula (38), parameters (R', K)0’)=(6×1020s212), and the generated hydrogen clock and cesium clock sequences, the output sequence of the DPLL, i.e., the sequence of the hydrogen clock after mounting, can be derived.
Fig. 5 shows the time difference of the hydrogen clock, cesium clock, and post-handling hydrogen clock for the first 4800 data points. Fig. 6 shows the alan deviations of hydrogen chimes, cesium chimes, and post-drive hydrogen chimes. As can be seen in fig. 5, the DPLL can effectively use the cesium clock to drive the hydrogen clock with less driving error. As can be seen from fig. 6: the frequency stability of the hydrogen bell after driving represents the short-term stability of the hydrogen bell and the long-term stability of the cesium bell, and the short-term stability and the long-term stability of the hydrogen bell are only slightly deteriorated compared with the hydrogen bell and the cesium bell, respectively.
Adding amplitude 3 × 10 to original hydrogen clock sequence-21s-1The frequency drift of (2). Fig. 7 shows the alan bias of hydrogen bell, cesium bell, and post-handling hydrogen bell when there is a drift in hydrogen bell. Comparing fig. 6 and 7, it can be seen that the frequency drift does not affect the frequency stability of the post-driving hydrogen clock.
To more clearly demonstrate the handling error (i.e., the deviation of the hydrogen clock after handling relative to the cesium clock), hydrogen clock and cesium clock sequences were regenerated, with the parameters of each sequence being identical to the parameters of the original hydrogen clock and the original cesium clock, for a total of 200000 sampling points. By using the same DPLL method, a sequence of post-drive hydrogen bells was obtained. Fig. 8 shows the steering error. As can be seen in fig. 8: the maximum steering error is 30ns, and the standard deviation is 8.99 ns. In fact, the handling error is caused by phase noise of the untracked hydrogen clock and cesium clock, and the standard deviation of the handling error, i.e., jitter, can be accurately calculated by control theory. It can therefore be presumed that: in applications where paper time is steered by UTC, the steering error will be correspondingly smaller since both paper time and UTC are more stable, containing less equivalent SSB phase noise components.
Adding amplitude 3 × 10 to original hydrogen clock sequence-21s-1The frequency drift of (2). Fig. 9 illustrates the handling error when the hydrogen clock contains a frequency drift. Comparing fig. 8 and 9, it can be seen that: when the hydrogen clock has frequency drift, a steady state tracking is causedError, amplitude of about 20 ns.
In conclusion, the present invention can effectively use the reference atomic clock (or time scale) to drive the driven atomic clock (or time scale), and the performance of driving error and the frequency stability after driving is superior. The drift of the driven atomic clock (or time scale) causes a steady state tracking error whose value can be accurately calculated and compensated by other algorithms.
In summary, although the present invention has been disclosed in terms of the preferred embodiment (hydrogen clock driven by cesium clock), it should not be construed as limiting the present invention, but rather as providing various modifications and alterations to the present invention without departing from the spirit and scope of the present invention, as defined by the appended claims.