A three-dimensional weighted comprehensive atomic time measurement method and device

Abstract

The present invention discloses a three-dimensional weighted comprehensive atomic time measurement method and device, including: obtaining the clock error data of each atomic clock in a clock group, decomposing the clock error data into three dimensions: frequency error, frequency drift, and initial clock error superimposed with atomic clock noise. Construct an initial matrix according to the decomposition result, and weight the decomposition result according to the corresponding weighting rule. Construct a weight matrix according to the weighted result, and perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to output a comprehensive atomic time. The present invention generates a comprehensive atomic time after three-dimensional decomposition, three-dimensional weight calculation, and three-dimensional weighting of the clock error data. The comprehensive atomic time takes into account the stability within multiple smoothing time periods, takes into account the long, medium and short-term stability performance, and improves the long, medium and short-term stability of the generated time scale, which can be used as the basis for adjusting the atomic clock to adjust the atomic clock to obtain the time and frequency actually used.

Description

Technical Field

[0001] The invention relates to the field of time frequency, and in particular to a three-dimensional weighted comprehensive atomic time measurement method. Background Art

[0002] With the advancement of science and technology, countries now generally establish timekeeping laboratories and widely use cesium atomic clocks and hydrogen atomic clocks to establish comprehensive atomic time. However, due to the risk of failure of actual physical clocks such as single cesium atomic clocks and single hydrogen atomic clocks, such as each individual clock may suddenly fail. Therefore, a clock group is usually used to give a comprehensive time scale through calculation. At present, the comprehensive atomic time algorithms used at home and abroad are mainly: weighted average algorithm, Kalman filter algorithm, wavelet decomposition algorithm. The weighted average algorithm is a classic combined clock time scale algorithm.

[0003] Typical weighted average algorithms include the ALGOS algorithm and the exponential filtering algorithm (AT1 algorithm). The ALGOS algorithm calculates the weight of each clock based on the frequency stability of each atomic clock; the AT1 algorithm calculates the weight of each clock based on the predicted error of the atomic clock. The predicted error refers to the absolute deviation between the predicted clock error of each clock and the reference clock and the actual clock error.

[0004] Since the weighted average algorithm takes a single weight for the atomic clock in each calculation cycle, it can only suppress one major noise or optimize the output time scale stability at a specified smoothing time point. It cannot consider the stability in multiple smoothing time periods, and therefore cannot take into account both long-term and short-term stability performance. In situations where short-term stability is important, the weighted average algorithm takes weights based on short-term stability, resulting in a large weight for the hydrogen clock and a small weight for the cesium clock, which cannot fully utilize the excellent long-term stability of the cesium clock, resulting in good short-term stability for the generated time scale, but poor long-term stability; in situations where long-term stability is important, the weighted average algorithm takes weights based on long-term stability, resulting in a large weight for the cesium clock and a small weight for the hydrogen clock, which cannot fully utilize the excellent short-term stability of the hydrogen clock, resulting in good long-term stability for the generated time scale, but poor short-term stability. Summary of the invention

[0005] The technical problem to be solved by the present invention is to provide a three-dimensional weighted integrated atomic time measurement method and device, which can take into account both long-term and short-term stability and improve the stability performance of the integrated atomic time.

[0006] In order to solve the above technical problems, the present invention provides a three-dimensional weighted comprehensive atomic time measurement method and device, comprising:

[0007] Obtain the clock error data of each atomic clock in the clock group;

[0008] Decomposing the clock error data into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise, and constructing an initial matrix according to the decomposition results;

[0009] According to a preset weighting rule, the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise are weighted respectively, and a weight matrix is ​​constructed according to the weighting results;

[0010] According to the weight matrix and the initial matrix, each atomic clock in the clock group is weighted averaged to obtain the comprehensive atomic time.

[0011] Furthermore, the clock error data is decomposed into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise, specifically:

[0012] Calculating the clock error data by sliding according to the least square method to obtain the frequency drift;

[0013] Calculating the clock difference data according to the least squares fitting method to obtain the frequency difference;

[0014] Subtraction calculation is performed based on the clock error data, the frequency drift and the frequency error to obtain the initial clock error of the superimposed atomic clock noise.

[0015] Furthermore, the initial matrix is ​​constructed according to the weighting results, specifically:

[0016] The initial matrix includes a decomposition matrix and a coefficient matrix;

[0017] Constructing the decomposition matrix according to the frequency difference, frequency drift and initial clock difference of superimposed atomic clock noise corresponding to each atomic clock;

[0018] The coefficient matrix is ​​constructed according to the clock error data of each atomic clock.

[0019] Furthermore, according to a preset weighting rule, the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise are weighted respectively, and a weight matrix is ​​constructed according to the weighting result, specifically:

[0020] According to a short-term stability rule, a first weight is assigned to the initial clock error of the superimposed atomic clock noise;

[0021] According to the mid-term stability rule, a second weight is assigned to the frequency difference;

[0022] According to the long-term stability rule, a third weight is applied to the frequency drift;

[0023] The weight matrix is ​​constructed according to the first weight, the second weight and the third weight.

[0024] Furthermore, according to the weight matrix and the initial matrix, each atomic clock in the clock group is weighted averaged to obtain the comprehensive atomic time, which is specifically:

[0025] According to the following calculation formula, the weighted average of each atomic clock is used to obtain the comprehensive atomic time:

[0026] W = (A.*WT)′[1 1 1]′;

[0027] In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the atomic time after weight adjustment; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix;

[0028] The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

[0029] Furthermore, the present invention also provides a three-dimensional weighted comprehensive atomic time measurement device, comprising: an acquisition module, a decomposition module, a weight module and a calculation module;

[0030] The acquisition module is used to acquire the clock error data of each atomic clock in the clock group;

[0031] The decomposition module is used to decompose the clock error data into three dimensions: frequency error, frequency drift and initial clock error superimposed with atomic clock noise, and construct an initial matrix according to the decomposition result;

[0032] The weight module is used to weight the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise according to a preset weighting rule, and construct a weight matrix according to the weighting result;

[0033] The comprehensive atomic time calculation module is used to perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the comprehensive atomic time.

[0034] Furthermore, the decomposition module includes a decomposition unit and an initial matrix construction unit, specifically:

[0035] The decomposition unit includes a first decomposition subunit, a second decomposition subunit and a third decomposition subunit;

[0036] The first decomposition subunit is used to calculate the clock error data according to the least squares sliding method to obtain the frequency drift;

[0037] The second decomposition subunit is used to calculate the clock error data according to the least squares fitting method to obtain the frequency error;

[0038] The third decomposition subunit is used to perform subtraction calculation based on the clock error data, the frequency drift and the frequency difference to obtain the initial clock error of the superimposed atomic clock noise.

[0039] Further, the initial matrix construction unit includes a first matrix construction sub-unit and a second matrix sub-construction unit, specifically:

[0040] The first matrix construction subunit is used to construct a decomposition matrix, and the decomposition matrix is ​​constructed according to the frequency difference, frequency drift and initial clock difference superimposed with atomic clock noise corresponding to each atomic clock;

[0041] The second matrix construction subunit is used to construct a coefficient matrix according to the clock error data of each atomic clock.

[0042] Furthermore, the weight module includes a weighting unit and a matrix construction unit, specifically:

[0043] The weighting unit takes a first weight for the initial clock error of the superimposed atomic clock noise according to the short-term stability rule; takes a second weight for the frequency error according to the medium-term stability rule; and takes a third weight for the frequency drift according to the long-term stability rule;

[0044] The matrix construction unit constructs the weight matrix according to the first weight, the second weight and the third weight.

[0045] Furthermore, the comprehensive atomic time calculation module is used to perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the comprehensive atomic time as follows:

[0046] According to the following calculation formula, the comprehensive atomic time obtained by each atomic clock after weighted averaging is calculated:

[0047] Y = (A.*WT)′[1 1 1]′;

[0048] In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the comprehensive atomic time obtained after weighted averaging; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix;

[0049] The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

[0050] Compared with the prior art, the three-dimensional weighted comprehensive atomic time measurement method and device of the embodiment of the present invention has the beneficial effect of: decomposing the clock error data in three dimensions into data of three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise. Then, the generated frequency error, frequency drift, and initial clock error with superimposed atomic clock noise are weighted according to medium-term stability, long-term stability, and short-term stability, respectively, to obtain corresponding weights. After three-dimensional decomposition and three-dimensional weight calculation, three-dimensional weighted averaging is performed to obtain comprehensive atomic time. The comprehensive atomic time can take into account the stability within multiple smoothing time periods, taking into account the long, medium and short-term stability performance, while improving the long, medium and short-term stability of the generated time scale, which can be used as the basis for adjusting the atomic clock in the clock group, and adjusting the atomic clock to obtain the actual time and frequency used. BRIEF DESCRIPTION OF THE DRAWINGS

[0051] Figure 1 It is a flow chart of an embodiment of a three-dimensional weighted integrated atomic time measurement method provided by an embodiment of the present invention;

[0052] Figure 2 This is a frequency stability comparison diagram of the integrated atomic time obtained by a three-dimensional weighted integrated atomic time measurement method provided by an embodiment of the present invention and the integrated atomic time obtained by the ALGOS and AT1 algorithms;

[0053] Figure 3 It is a structural schematic diagram of an embodiment of a three-dimensional weighted integrated atomic time generation device provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0054] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0055] Example 1

[0056] See also Figure 1 , Figure 1 1 is a flow chart of an embodiment of a three-dimensional weighted integrated atomic time measurement method provided by an embodiment of the present invention. The method includes steps 101 to 105, and each step is specifically as follows:

[0057] Step 101: Obtain clock error data of each atomic clock in the clock group.

[0058] The clock group is composed of at least three atomic clocks, and can be applied to various clock group configurations such as two hydrogen and one cesium, two cesium and one hydrogen, and multiple cesium and multiple hydrogen. In the embodiment of the present invention, the clock group is composed of two hydrogen clocks H1 and H2 and one cesium clock Cs. The clock error data of each atomic clock in the clock group is obtained, and the clock error data of the hydrogen clock H1 is The clock error data of hydrogen clock H2 is The clock error data of the cesium clock Cs is

[0059] Step 102: Decompose the clock error data into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise.

[0060] In an embodiment of the present invention, and It represents the initial clock error of the three atomic clocks in the clock group, which is the sum of the atomic clock noise; and represents the frequency difference of the three atomic clocks in the clock group; and The relationship between the clock error data X and the initial clock error a0, frequency error a1 and frequency drift a2 superimposed with atomic clock noise can be expressed by the following formula:

[0061]

[0062] Where t is time.

[0063] In the embodiment of the present invention, the relationship between the clock error of each atomic clock and the initial clock error, frequency error and frequency drift of the corresponding superimposed atomic clock noise can be expressed by the following formula:

[0064] The first relation:

[0065] The second relation:

[0066] The third relation:

[0067] The frequency drift a2 is calculated by sliding according to the least square method; the frequency difference a1 is calculated by fitting according to the least square method; the initial clock difference a0 of the superimposed atomic clock noise is calculated according to the relationship between the clock difference data X and the initial clock difference a0 of the superimposed atomic clock noise, the frequency difference a1 and the frequency drift a2.

[0068] Step 103: construct an initial matrix according to the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise.

[0069] The initial matrix includes a decomposition matrix and a coefficient matrix.

[0070] The decomposition matrix is ​​constructed according to the frequency difference a1, frequency drift a2 and initial clock difference a0 superimposed with atomic clock noise corresponding to each atomic clock.

[0071]

[0072] In the embodiment of the present invention, A represents a decomposition matrix, which is obtained by performing a three-dimensional decomposition of the clock error of each atomic clock in the clock group, and is composed of an initial clock error a0, a frequency error a1, and a frequency drift a2 superimposed with atomic clock noise. Each row of parameters represents the result of the three-dimensional decomposition of an atomic clock. There are three columns of parameters, which are the initial clock error a0, frequency error a1, and frequency drift a2 of the superimposed atomic clock noise corresponding to each atom.

[0073] According to the relationship formula between the clock difference data X and the initial clock difference a0, frequency difference a1 and frequency drift a2 superimposed with atomic clock noise, the coefficients of the initial clock difference a0, frequency difference a1 and frequency drift a2 superimposed with atomic clock noise in the clock difference data X are extracted to construct a coefficient matrix.

[0074]

[0075] T represents the coefficient matrix. The first row of parameters indicates that the initial clock difference a0 superimposed with atomic clock noise appears as a noise superposition constant term in the clock difference data; the second row of parameters indicates that the frequency difference a1 appears as a linear trend term in the clock difference data; the third row of parameters indicates that the frequency drift a2 appears as a quadratic term in the clock difference data.

[0076] Step 104: According to a preset weighting rule, the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise are weighted respectively, and a weight matrix is ​​constructed according to the weighting results.

[0077] In order to ensure short-term stability, the initial clock error a0 superimposed with atomic clock noise is weighted according to the short-term stability, which can be calculated using the following formula:

[0078] First right formula:

[0079] In the formula, represents the weight of the initial clock error a0 term superimposed with atomic clock noise, σ 2 (1000s) indicates short-term stability with a smoothing time of 1000s.

[0080] In order to ensure medium-term stability, the frequency difference a1 is weighted according to the medium-term stability. Since the intersection of the stability curves of the hydrogen clock and the cesium clock is near the smoothing time of 100,000 seconds, it indicates that the stability performance of the hydrogen clock and the cesium clock is close near the smoothing time of 100,000 seconds. Therefore, in the embodiment of the present invention, the frequency difference a1 is weighted according to the stability with a smoothing time of 100,000 seconds, which can be expressed by the following formula:

[0081] The second weighting formula:

[0082] In the formula, represents the weight of the frequency difference a1 term, σ 2 (100000s) means the smoothing time is 100000s of stability.

[0083] In order to ensure long-term stability, the frequency drift a2 is weighted according to the long-term stability. The absolute value of the frequency drift a2 reflects the magnitude of the frequency drift. The larger the absolute value of the frequency drift a2, the larger the frequency drift and the worse the long-term stability. In order to improve long-term stability, the larger the absolute value of the frequency drift a2, the smaller the weight should be assigned. In the embodiment of the present invention, the inverse of the absolute value of the frequency drift a2 is used as the weight of the frequency drift a2, which can be expressed by the following formula:

[0084] The third weighting formula:

[0085] In the formula, Represents the weight of the frequency drift a2 term.

[0086] The weights corresponding to the initial clock error a0, frequency error a1 and frequency drift a2 of each atomic clock in the clock group superimposed with the atomic clock noise are calculated according to the first weighting formula, the second weighting formula and the third weighting formula, and a weight matrix is ​​constructed according to the weighting results.

[0087] In the embodiment of the present invention, a weight matrix is ​​constructed according to the weighting result as follows:

[0088]

[0089] W represents the weight matrix. Each row of parameters represents the corresponding weights of the initial clock error a0, frequency error a1 and frequency drift a2 of an atomic clock superimposed with the atomic clock noise. Parameters in different rows represent the weights of different atomic clocks.

[0090] Step 105: Perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the integrated atomic time, specifically:

[0091] According to the following calculation formula, the weighted average of each atomic clock is calculated to obtain the comprehensive atomic time:

[0092] Y = (A.*WT)′[1 1 1]′;

[0093] In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the atomic time after weight adjustment; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix;

[0094] In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the comprehensive atomic time obtained after weighted averaging; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix;

[0095] The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

[0096] Stability represents the ability to maintain a constant scale interval. Stability is represented by Allan Variance in the time domain. The calculation formula of Allan Variance is:

[0097]

[0098] Where τ represents the smoothing time, which is usually an integer multiple of the measurement time interval; x(i+1) and x(i) are time differences, and their time interval is τ; represents the Allan variance; N represents the number of time difference data points. The larger the Allan variance value, the larger the time scale fluctuation of the corresponding smoothing time and the worse the stability; the smaller the Allan variance value, the smaller the time scale fluctuation of the corresponding smoothing time and the better the stability.

[0099] In the embodiment of the present invention, the Allan variance of the comprehensive atomic time Y obtained by calculating the three-dimensional weighted average is Compare the Allan variances corresponding to different smoothing times and evaluate the stability of the comprehensive atomic time Y obtained by three-dimensional weighted averaging. Figure 2 , Figure 2 This is a comparison chart of the frequency stability of the comprehensive atomic time obtained by a three-dimensional weighted comprehensive atomic time calculation method provided by an embodiment of the present invention and the comprehensive atomic time calculated by the ALGOS and AT1 algorithms. Curve 1 corresponds to the ALGOS algorithm, which calculates the weight of each clock based on the frequency stability of each atomic clock; Curve 2 corresponds to the AT1 algorithm, which calculates the weight of each clock based on the predicted error of the atomic clock, and the predicted error refers to the absolute deviation between the predicted clock difference of each clock and the reference clock and the actual clock difference; Curve 3 corresponds to a three-dimensional weighted comprehensive atomic time measurement method according to an embodiment of the present invention. According to the stability of each curve at different smoothing times, it can be seen that the comprehensive atomic time Y output by three-dimensional weighted average takes into account the long-, medium- and short-term stability performance, and the long-, medium- and short-term stability at different smoothing times are all improved.

[0100] Example 2

[0101] See also Figure 3 , Figure 3The schematic diagram of the structure of an embodiment of a three-dimensional weighted integrated atomic time generation device provided by the embodiment of the present invention comprises: an acquisition module 301 , a decomposition module 302 , a weight module 303 and a calculation module 304 .

[0102] In the embodiment of the present invention, the acquisition module 301 is used to acquire the clock error data of each atomic clock in the clock group. The clock group is composed of two hydrogen clocks H1 and H2 and a cesium clock Cs. The clock error data of each atomic clock in the clock group is acquired. The clock error data of the hydrogen clock H1 is The clock error data of hydrogen clock H2 is The clock error data of the cesium clock Cs is

[0103] In the embodiment of the present invention, the decomposition module 302 is used to decompose the clock error data into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise, and construct an initial matrix according to the decomposition results. The decomposition module 302 includes a decomposition unit and an initial matrix construction unit. The decomposition unit includes a first decomposition subunit, a second decomposition subunit, and a third decomposition subunit. In the embodiment of the present invention, and It represents the initial clock error of the three atomic clocks in the clock group, which is the sum of the atomic clock noise; and represents the frequency difference of the three atomic clocks in the clock group; and The relationship between the clock error data X and the initial clock error a0, frequency error a1 and frequency drift a2 superimposed with atomic clock noise can be expressed by the following formula:

[0104]

[0105] Where t is time.

[0106] In the embodiment of the present invention, the relationship between the initial clock error, frequency error and frequency drift of each atomic clock and its corresponding superimposed atomic clock noise can be expressed by the following formula:

[0107] The first relation:

[0108] The second relation:

[0109] The third relation:

[0110] The frequency drift a2 is calculated by sliding according to the least square method; the frequency difference a1 is calculated by fitting according to the least square method; the initial clock difference a0 of the superimposed atomic clock noise is calculated according to the relationship between the clock difference data X and the initial clock difference a0 of the superimposed atomic clock noise, the frequency difference a1 and the frequency drift a2.

[0111] In the embodiment of the present invention, the decomposition module 302 further includes an initial matrix construction unit, which is used to construct an initial matrix according to the decomposition result. The initial matrix construction unit includes a first matrix construction subunit and a second matrix construction subunit.

[0112] The first matrix construction subunit is used to construct a decomposition matrix according to the frequency difference a1, frequency drift a2 and initial clock difference a0 superimposed with atomic clock noise corresponding to each atomic clock.

[0113]

[0114] In the embodiment of the present invention, A represents a decomposition matrix, which is obtained by performing a three-dimensional decomposition of the clock error of each atomic clock in the clock group, and is composed of an initial clock error a0, a frequency error a1, and a frequency drift a2 superimposed with atomic clock noise. Each row of parameters represents the result of the three-dimensional decomposition of an atomic clock. There are three columns of parameters, which are the initial clock error a0, frequency error a1, and frequency drift a2 of the superimposed atomic clock noise corresponding to each atom.

[0115] The second matrix construction subunit is used to construct a coefficient matrix. According to the relationship formula between the clock difference data X and the initial clock difference a0, frequency difference a1 and frequency drift a2 of the superimposed atomic clock noise, the coefficients of the initial clock difference a0, frequency difference a1 and frequency drift a2 of the superimposed atomic clock noise in the clock difference data X are extracted to construct the coefficient matrix.

[0116]

[0117] T represents the coefficient matrix. The first row of parameters indicates that the initial clock difference a0 superimposed with atomic clock noise appears as a noise superposition constant term in the clock difference data; the second row of parameters indicates that the frequency difference a1 appears as a linear trend term in the clock difference data; the third row of parameters indicates that the frequency drift a2 appears as a quadratic term in the clock difference data.

[0118] In the embodiment of the present invention, the weight module 303 is used to weight the frequency difference a1, the frequency drift a2 and the initial clock difference a0 of the superimposed atomic clock noise according to a preset weighting rule, and construct a weight matrix according to the weighting result. The weight module 303 includes a weighting unit and a matrix construction unit.

[0119] In order to ensure short-term stability, the weighting unit takes the first weight for the initial clock error a0 superimposed with atomic clock noise according to the short-term stability rule, which can be calculated using the following formula:

[0120] First right formula:

[0121] In the formula, represents the weight of the initial clock error a0 term superimposed with atomic clock noise, and σ2(1000s) represents the short-term stability with a smoothing time of 1000s.

[0122] In order to ensure the medium-term stability, the weighting unit takes the second weight for the frequency difference a1 according to the medium-term stability taking rule. Since the intersection of the stability curves of the hydrogen clock and the cesium clock is near the smoothing time of 100,000 seconds, it indicates that the stability performance of the hydrogen clock and the cesium clock is close near the smoothing time of 100,000 seconds. Therefore, in the embodiment of the present invention, the frequency difference a1 is weighted by the stability with a smoothing time of 100,000 seconds, which can be expressed by the following formula:

[0123] The second weighting formula:

[0124] In the formula, represents the weight of the frequency difference a1 term, σ 2 (100000s) means the smoothing time is 100000s of stability.

[0125] In order to ensure long-term stability, the frequency drift a2 is given a third weight according to the long-term stability. The absolute value of the frequency drift a2 reflects the magnitude of the frequency drift. The larger the absolute value of the frequency drift a2, the larger the frequency drift and the worse the long-term stability. In order to improve long-term stability, the larger the absolute value of the frequency drift a2, the smaller the weight should be assigned. In the embodiment of the present invention, the inverse of the absolute value of the frequency drift a2 is used as the weight of the frequency drift a2, which can be expressed by the following formula:

[0126] The third weighting formula:

[0127] In the formula, Represents the weight of the frequency drift a2 term.

[0128] The construction unit constructs the weight matrix according to the first weight, the second weight and the third weight, and the result is as follows:

[0129]

[0130] W represents the weight matrix. Each row of parameters represents the corresponding weights of the initial clock error a0, frequency error a1 and frequency drift a2 of an atomic clock superimposed with the atomic clock noise. Parameters in different rows represent the weights of different atomic clocks.

[0131] In the embodiment of the present invention, the operation module 304 is used to perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the comprehensive atomic time, specifically:

[0132] According to the following calculation formula, the comprehensive atomic time obtained after weighted averaging for each atomic clock is calculated:

[0133] Y = (A.*WT)′[1 1 1]′;

[0134] In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the comprehensive atomic time obtained after weighted averaging; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix;

[0135] The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

[0136] In summary, the embodiment of the present invention provides a three-dimensional weighted atomic time generation device, which decomposes the atomic time clock error data according to the three dimensions of initial clock error, frequency error and frequency drift superimposed with atomic time noise, and then calculates weights for the three dimensions respectively. It overcomes the defect that the weighted average algorithm takes a single weight for the atomic clock in each calculation cycle, can only suppress one major noise or can only make the output time scale stability of a certain smooth time point reach the optimal level. The generated comprehensive atomic time can take into account the stability performance of the long, medium and short term.

[0137] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and substitutions can be made without departing from the technical principles of the present invention. These improvements and substitutions should also be regarded as the scope of protection of the present invention.

[0138] Those skilled in the art can understand that all or part of the processes in the above embodiments can be implemented by instructing related hardware through a computer program, and the program can be stored in a computer-readable storage medium, and when the program is executed, it can include the processes in the above embodiments. The storage medium can be a disk, an optical disk, a read-only memory (ROM) or a random access memory (RAM), etc.

Claims

1. A three-dimensional weighted comprehensive atomic time measurement method, characterized in that: include: Obtain the clock error data of each atomic clock in the clock group; Decomposing the clock error data into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise, and constructing an initial matrix according to the decomposition results; Specifically: the initial matrix includes a decomposition matrix and a coefficient matrix; the decomposition matrix is ​​constructed according to the frequency difference, frequency drift and initial clock difference of superimposed atomic clock noise corresponding to each atomic clock; Constructing the coefficient matrix according to the clock error data of each atomic clock; According to a preset weighting rule, the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise are weighted respectively, and a weight matrix is ​​constructed according to the weighting results; specifically, according to a short-term stability rule, a first weight is taken for the initial clock difference of the superimposed atomic clock noise; according to a medium-term stability rule, a second weight is taken for the frequency difference; according to a long-term stability rule, a third weight is taken for the frequency drift; and the weight matrix is ​​constructed according to the first weight, the second weight and the third weight; According to the weight matrix and the initial matrix, each atomic clock in the clock group is weighted averaged to obtain the comprehensive atomic time.

2. A three-dimensional weighted integrated atomic time measurement method as claimed in claim 1, characterized in that: The clock error data is decomposed into three dimensions: frequency error, frequency drift, and initial clock error with superimposed atomic clock noise, specifically: Calculating the clock error data by sliding according to the least square method to obtain the frequency drift; Calculating the clock difference data according to the least squares fitting method to obtain the frequency difference; Subtraction calculation is performed based on the clock error data, the frequency drift and the frequency error to obtain the initial clock error of the superimposed atomic clock noise.

3. A three-dimensional weighted integrated atomic time measurement method as claimed in claim 1, characterized in that: According to the weight matrix and the initial matrix, each atomic clock in the clock group is weighted averaged to obtain the comprehensive atomic time, which is specifically: According to the following calculation formula, the comprehensive atomic time obtained by weighted averaging of each atomic clock is calculated: Y = (A.*WT)′[1 1 1]′; In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the comprehensive atomic time obtained after weighted averaging; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix; The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

4. A three-dimensional weighted integrated atomic time measurement device, characterized in that: include: Acquisition module, decomposition module, weight module and operation module; The acquisition module is used to acquire the clock error data of each atomic clock in the clock group; The decomposition module is used to decompose the clock difference data into three dimensions of frequency difference, frequency drift and initial clock difference of superimposed atomic clock noise, and construct an initial matrix according to the decomposition result; wherein the initial matrix construction unit includes a first matrix construction subunit and a second matrix construction subunit, specifically: the first matrix construction subunit is used to construct a decomposition matrix, and the decomposition matrix is ​​constructed according to the frequency difference, frequency drift and initial clock difference of superimposed atomic clock noise corresponding to each atomic clock; the second matrix construction subunit is used to construct a coefficient matrix, and the coefficient matrix is ​​constructed according to the clock difference data of each atomic clock; The weight module is used to weight the frequency difference, the frequency drift and the initial clock difference of the superimposed atomic clock noise according to a preset weighting rule, and construct a weight matrix according to the weighting results; wherein the weight module includes a weighting unit and a matrix construction unit, specifically: the weighting unit takes a first weight for the initial clock difference of the superimposed atomic clock noise according to a short-term stability rule; takes a second weight for the frequency difference according to a medium-term stability rule; takes a third weight for the frequency drift according to a long-term stability rule; the matrix construction unit constructs the weight matrix according to the first weight, the second weight and the third weight; The operation module is used to perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the comprehensive atomic time.

5. A three-dimensional weighted integrated atomic time measurement device as claimed in claim 4, characterized in that: The decomposition module includes a decomposition unit and an initial matrix construction unit, specifically: The decomposition unit includes a first decomposition subunit, a second decomposition subunit and a third decomposition subunit; The first decomposition subunit is used to calculate the clock error data according to the least squares sliding method to obtain the frequency drift; The second decomposition subunit is used to calculate the clock error data according to the least squares fitting method to obtain the frequency error; The third decomposition subunit is used to perform subtraction calculation based on the clock error data, the frequency drift and the frequency difference to obtain the initial clock error of the superimposed atomic clock noise.

6. A three-dimensional weighted integrated atomic time measurement device as claimed in claim 4, characterized in that: The operation module is used to perform weighted averaging on each atomic clock in the clock group according to the weight matrix and the initial matrix to obtain the comprehensive atomic time, specifically: According to the following calculation formula, the comprehensive atomic time obtained after weighted averaging for each atomic clock is calculated: Y = (A.*WT)′[1 1 1]′; In the formula, the operator .* is defined to represent the multiplication operation of the corresponding elements of two matrices of the same dimension; the operator ' is defined to represent the matrix transposition operation; Y is the comprehensive atomic time obtained after weighted averaging; A is the decomposition matrix; W is the weight matrix; T is the coefficient matrix; The comprehensive atomic time Y obtained after weighted averaging is used as the basis for adjusting the atomic clock driving, and the atomic clocks in the clock group are adjusted to obtain the time and frequency actually used.

Images