A method, device and equipment for continuously calculating world time and a storage medium

Abstract

The application discloses a world time continuous solution method, comprising: obtaining and preprocessing the solution result files of each session in a predetermined time range; applying uniform constraints to each session after preprocessing, and performing parameter blocking and local parameter elimination; constructing a sliding time window in a target time range, parameterizing the world time by using a segmented linear function in the window, and establishing a session and time node mapping relationship; constructing a sliding time window on a global time axis, collecting sessions in the window, superimposing a session weight equation, adding a time node smoothing constraint, and jointly solving the world time residual increment; superimposing a priori on the joint solution result, recovering a tidal term, and generating a world time continuous sequence. The application solves the problems of the existing world time product, such as limited time resolution based on daily parameter level combination, insufficient utilization of session level observation geometry and equation information, sensitivity to modeling differences of different session types, and difficulty in adapting to the existing reference framework.

Description

Technical Field

[0001] This invention belongs to the field of space geodesy, and particularly relates to a method, apparatus, equipment and storage medium for continuous calculation of world time. Background Technology

[0002] Earth orientation parameters are key transformation parameters characterizing the spatial orientation of the International Earth Reference Frame relative to the International Celestial Reference Frame over time. These parameters mainly include polar motion, Universal Time (UTC), day length variation, and celestial pole offset. UTC, in particular, directly reflects the Earth's rotation angle and is a crucial parameter for high-precision applications such as establishing and maintaining time standards, satellite navigation and precise orbit determination, earthquake monitoring, and geodynamics research. Compared to other Earth orientation parameter components, UTC exhibits significant irregular variations on short timescales and has poor predictability. Currently, it can only be independently and accurately determined using Very Long Baseline Interferometry (VLBI) technology, utilizing a quasi-inertial reference frame constructed from distant quasars.

[0003] Existing UTC operational products are primarily released by the International Earth Rotation and Reference System Service (IERS), including the final series of products from the Earth Rotation and Reference Systems Centre (IERS C04) and rapid combination and forecast products from the Rapid Services and Forecasts Centre (IERS Bulletin A). These products typically employ a daily parameter-level combination strategy, which statistically weights daily UTC solutions provided by different space geodetic techniques or analysis centers, retaining only the daily UTC scale and its formal uncertainty. While this strategy is simple to implement and stable, it has the following limitations: First, session-level observational geometric information, normal equations, and parameter correlations are significantly weakened or discarded during the combination process. Differences in constraint strategies, parameterization methods, and stochastic models among different techniques or analysis centers are compressed into a single daily parameter, hindering optimal weighting and easily introducing systematic biases. Second, daily-scale combinations have inherent limitations in temporal resolution, making it difficult to effectively characterize intraday and higher-frequency UTC variations. High-frequency signals are easily smoothed or aliased, thus affecting the timeliness and short-period accuracy of the products.

[0004] At the observational level, the integrated UTC product provided by IERS still fundamentally relies primarily on very long baseline interferometry (VLBI) observations, especially the 1-hour enhanced observation sessions (INT1 / INT2 / INT3) and 24-hour regular observation sessions (R1 / R4) organized by the International VLBI Analysis Centre (IVS). The 24-hour regular sessions are typically conducted twice a week, employing a global multi-baseline network comprising approximately ten or more radio telescopes. This provides excellent all-sky coverage and high redundancy, enabling robust separation of UTC, tropospheric zenith delay, station clock errors, and station coordinate corrections within a single session, thus jointly solving the complete set of Earth orientation parameters. However, the operational delay from the end of the observation to the generation of SINEX equations or formal solutions for this type of session is typically around two weeks, and time sampling is only twice a week, making it difficult to meet the needs of quasi-continuous monitoring of UTC. To improve the frequency of UTC updates, IVS schedules almost daily 1-hour enhanced observation sessions, utilizing a few long east-west baselines to achieve rapid UTC calculation, offering advantages such as timely updates and lower operational delays. However, due to the limitations of short-arc observations and small network structures consisting of only 2–3 stations, the enhanced session observation geometry is weak and the redundancy is limited. It is generally only suitable for estimating UTC (or UTC under strong constraints), and its solution accuracy is usually about 2–3 times lower than that of a complete conventional observation multi-baseline network session.

[0005] Under the aforementioned business model, existing technologies generally adopt an operational scheme of "enhanced observation for rapid updates + high-precision supplementation by conventional observations," but there are still significant shortcomings: On the one hand, enhanced observations are short-arc, small-network observations, which are highly sensitive to system effects such as the status of individual station equipment, tropospheric modeling errors, and clock anomalies. The Universal Time (UTC) solution is easily affected by station network geometric degradation and residual modeling errors, resulting in limited stability and reliability. On the other hand, although conventional observations can provide more accurate UTC and complete Earth orientation parameter solutions, their observation intervals and data processing delays are relatively long, making it difficult to reflect rapid changes in UTC in a timely manner. In addition, when enhanced and conventional observations use simple parameter-level parameter combinations at the daily scale, the differences between different session types in terms of observation modeling, constraint methods, and stochastic models are often "folded" into the estimation of UTC on a single day. This is not conducive to fully leveraging the advantages of various sessions and can easily lead to smoothing or aliasing of high-frequency UTC signals, affecting the short-cycle accuracy and timeliness of the product. Summary of the Invention

[0006] The purpose of this invention is to overcome the problems of limited time resolution, insufficient utilization of session-level observation geometry and normal equation information, sensitivity to modeling differences between different session types, and difficulty in adapting to the current reference frame in existing Universal Time (UTC) products based on daily parameter-level combinations. This invention proposes a method to achieve high time resolution, continuity, and high accuracy in UTC solutions by utilizing the normal equations of different session types in VLBI under the current reference frame and constraint strategy, thereby obtaining a continuous UTC sequence that balances timeliness and accuracy.

[0007] According to one aspect of the present invention, a method for continuous calculation of Universal Time is provided, comprising:

[0008] Obtain the solution result files for each session within a predetermined time range and perform preprocessing;

[0009] Apply uniform constraints to each session after preprocessing, and perform parameter block division and local parameter elimination;

[0010] Construct a sliding time window within the target time range, parameterize world time using a piecewise linear function within the window, and establish a mapping relationship between sessions and time nodes;

[0011] Construct a sliding time window on the global time axis, collect sessions within the window, superimpose the equations according to session weights, add time node smoothing constraints, and jointly solve the world time margin increment.

[0012] Priors are superimposed on the joint solution results, and the tidal term is recovered to generate a continuous sequence of universal time.

[0013] As a further technical solution, the solution result files of each session within a predetermined time range are obtained and preprocessed, including:

[0014] Based on the solution result files of each VLBI session within a predetermined time range, the session normal equations are written as follows: ,in, Let be the normal equation coefficient matrix for the s-th session; Let be the column vector of parameters to be estimated for the s-th session; Let be the column vector of constant terms in the normal equation for the s-th session;

[0015] All station coordinates for all sessions are unified under the same Earth reference frame, and Universal Time tidal correction and incremental linearization are performed.

[0016] As a further technical solution, Universal Time tidal correction and incremental linearization processing are performed, including:

[0017] Tidal correction is applied to Universal Time to obtain the Universal Time margin at time t. for: ,in, Let t be the world time. Let t be the world time tidal correction at time t;

[0018] Based on the prior value and observed value of the world time margin, the world time margin increment is obtained, and the world time margin increment is substituted into the conversational method equation as one of the world time unknowns.

[0019] As a further technical solution, the column vector of parameters to be estimated includes at least: the world time surplus increment, the clock difference parameters of each VLBI station, and the tropospheric zenith delay parameters of each VLBI station.

[0020] As a further technical solution, a piecewise linear function is used within the window to parameterize world time, establishing a mapping relationship between sessions and time nodes, including:

[0021] At any point between adjacent time nodes within the window, a piecewise linear function is used to approximate the increment of the world time balance;

[0022] Introducing a representative epoch for the s-th session, where the representative epoch is located between the aforementioned adjacent time nodes, the world time residual increment is transformed using a piecewise linear function to obtain an approximate expression of the world time residual increment for the s-th session.

[0023] By introducing the row vector corresponding to the time node of the s-th session, the approximate expression of the world time margin increment of the s-th session is transformed to obtain the final approximate expression of the world time margin increment.

[0024] By substituting the final approximate expression of the world time margin increment into the reduction equation after local parameter elimination, a mapping relationship between sessions and time nodes is established.

[0025] As a further technical solution, the session weight superposition method equation is used, and time node smoothing constraints are added to jointly solve for the world time surplus increment, including:

[0026] Within the j-th time window, all items belonging to the set The normal equations of the conversation are weighted and superimposed to obtain the normal equation coefficient matrix synthesized within the j-th time window. And the column vector of normal equation constant terms synthesized within the j-th time window. ;

[0027] Time node parameters A second-order difference smoothing constraint is applied, where, Let T be the second-order difference operator matrix at time node k, and let T be the transpose. To smooth out the constraint weighting factors, Let be the column vector of the world time surplus increments for the time nodes to be estimated within the j-th time window. Then, the joint normal equations within the j-th time window are:

[0028] ,

[0029] Solving the joint method equations yields an estimate of the cosmic time surplus increment at the j-th time window:

[0030] ,

[0031] By sliding the time window along the time axis and repeating the solution process, continuous time node estimates are obtained over the entire target time range.

[0032] As a further technical solution, after obtaining the estimate of the world node's world time surplus increment, it also includes:

[0033] By using prior world time series to recover the world time margin and combining it with piecewise linear interpolation, a continuous world time margin estimate is obtained.

[0034] Based on the obtained Universal Time Residual Estimate, the tidal term is recovered, and a continuous Universal Time sequence is generated.

[0035] According to one aspect of the present invention, a continuous world time calculation apparatus is provided, comprising:

[0036] The first main module is used to obtain the solution result files of each session within a predetermined time range and perform preprocessing.

[0037] The second main module is used to apply uniform constraints to each session after preprocessing, and to perform parameter block division and local parameter elimination.

[0038] The third main module is used to construct a sliding time window within the target time range, parameterize world time using a piecewise linear function within the window, and establish a mapping relationship between sessions and time nodes.

[0039] The fourth main module is used to construct a sliding time window on the global time axis, collect sessions within the window, superimpose the normal equations according to session weights, add time node smoothing constraints, and jointly solve the world time margin increment.

[0040] The fifth main module is used to superimpose priors on the joint solution results and recover the tidal term to generate a continuous sequence of universal time.

[0041] According to one aspect of the present invention, a continuous commencing device for Universal Time (UTC) is provided, comprising a memory and a processor, wherein the memory stores program instructions that are executed by the processor, and the processor invokes the program instructions to execute the aforementioned continuous commencing method for UTC.

[0042] According to one aspect of the present invention, a non-transitory computer-readable storage medium is provided, the non-transitory computer-readable storage medium storing computer instructions that cause the computer to execute the aforementioned continuous world time calculation method.

[0043] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0044] This invention proposes a method for achieving high temporal resolution, continuity, and high precision solution of Universal Time (UTC) using VLBI normal equations for different session types under existing reference frames and constraint strategies. This method obtains a continuous UTC sequence that balances timeliness and accuracy, and solves the problems of limited temporal resolution, insufficient utilization of session-level observation geometry and normal equation information, sensitivity to modeling differences between different session types, and difficulty in adapting to existing reference frames in existing UTC products based on daily parameter-level combinations. Attached Figure Description

[0045] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the accompanying drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0046] Figure 1 This is a flowchart illustrating a continuous world time calculation method provided in an embodiment of the present invention.

[0047] Figure 2 This is a schematic diagram comparing the comprehensive solution of world time provided in an embodiment of the present invention with the official product. Detailed Implementation

[0048] The terms “comprising” and “having”, and any variations thereof, in the specification, claims, and accompanying drawings of this invention are intended to cover a non-exclusive inclusion, such as a process, method, system, product, or apparatus that includes a series of steps or units, not necessarily limited to those explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0049] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention. In addition, the technical features of the various embodiments or individual embodiments provided by the present invention can be arbitrarily combined to form new technical solutions. Such combinations are not bound by the order of steps and / or structural composition patterns, but must be based on the ability of those skilled in the art to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection claimed by the present invention.

[0050] This invention provides a method for continuous calculation of world time, referring to... Figure 1 As shown, the process includes five main steps: observation data acquisition and preprocessing, unified modeling of session-level normal equations, continuous parameterization of Universal Time and construction of time grids, fusion and joint solution of normal equations within a sliding time window, and restoration and generation of continuous Universal Time sequences. Under the existing reference framework and constraint strategy, the process utilizes normal equations of different session types of VLBI to achieve high temporal resolution, continuity, and high precision in Universal Time solution, thereby obtaining a continuous Universal Time sequence that balances timeliness and accuracy.

[0051] The term "session" in this invention refers to a single, independently scheduled and computed VLBI observation experiment (such as INT or RAP).

[0052] The specific implementation process of the continuous world time calculation method described in this embodiment of the invention is as follows:

[0053] Step 1: Acquisition and preprocessing of observation data.

[0054] 1) Acquisition of session-level normal equations

[0055] Preferably, the solution result files for each VLBI session within a predetermined time range are obtained from the International VLBI Analysis Center. Let there be S sessions within this time range, then the s-th session... The session-level normal equation is written as:

[0056] ,

[0057] In the formula, Let be the normal equation coefficient matrix for the s-th session; Let be the column vector of parameters to be estimated for the s-th session; Let be the column vector of constant terms in the normal equation for the s-th session. Where is the parameter vector. It should include at least: UTC-related parameters (such as UTC margin increments); clock error parameters for each VLBI station; tropospheric zenith delay parameters for each VLBI station; and, if necessary, station coordinate corrections.

[0058] 2) Reference frame and station coordinates are unified

[0059] To ensure consistency across different sessions, this step unifies the station coordinates of all sessions to the same Earth reference frame (e.g., ITRF2020).

[0060] 3) Universal Time Tidal Correction

[0061] To facilitate subsequent continuous processing, this invention removes short-period tidal effects from Universal Time (UTC) and defines a Universal Time Residual (UT1R). Let: For world time at time t; Let t be the world time tidal correction at time t; Let be the world time margin at time t. Then we have:

[0062] ,

[0063] In the formula, It can be calculated from a pre-given tidal model, and is generally written as the sum of multiple harmonic components:

[0064] ,

[0065] in, The number of tidal components used; for each tidal component : The amplitude of the cosine term; The amplitude of the sinusoidal term; This is the angular frequency of the component.

[0066] 4) Linearization of world time increment

[0067] To solve within a unified incremental framework, this invention selects a long-term universal time reference sequence as a priori value. Let: Let be the world time margin prior value at time t; let be the representative epoch of the s-th session. ; The world time margin observation value obtained from the solution of the s-th session; Let be the world time margin increment for the s-th session. Then the world time margin increment is defined as:

[0068] ,

[0069] Then, As one of the unknown quantities of world time, it is substituted into the conversational equation shown in formula (1) so that each conversation is represented on a unified prior basis.

[0070] Step 2: Unified modeling of session-level normal equations.

[0071] 1) Unified constraint strategy and normal equation correction

[0072] To ensure that each session solves under a unified reference frame and constraint rules, this invention applies uniform no-offset / no-rotation constraints to each session. Therefore:

[0073] ,

[0074] In the formula, Let be the constraint coefficient matrix for the s-th session; Let be the constraint matrix for the s-th session; Let be the column vector of constraint constant terms for the s-th session; The coefficient matrix of the normal equations after applying constraints; This is the vector of constant terms in the normal equation after constraints are applied.

[0075] 2) Parameter partitioning and local parameter elimination

[0076] The column vector of parameters to be estimated for the s-th session The blocks are as follows:

[0077] ,

[0078] In the formula, This is a column vector consisting of other parameters of the session, such as station clock error, tropospheric zenith delay, and station coordinate correction.

[0079] Substitute formula (6) into formula (5), and... The blocks are written as follows:

[0080] ,

[0081] In the formula, for The relevant scalar equation coefficients; This is a row vector of coupling coefficients between world time and other parameters; for The transpose (column vector); column vector The coefficient matrix of the normal equations.

[0082] Accordingly, The blocks are as follows:

[0083] ,

[0084] in, To and Related scalar constant terms; For column vectors The relevant constant term column vector.

[0085] By using Schul complement, the column vector Elimination yields only information about... The reduced equation:

[0086] ,

[0087] In the formula, The equivalent weight coefficients for world time direction after elimination are defined as follows:

[0088] ,

[0089] The equivalent constant term for world time direction after elimination is defined as:

[0090] ,

[0091] at this time, As a scalar, As scalars, they collectively represent the "equivalent observation" and "equivalent weight" of the s-th session in the world time direction.

[0092] Step 3: Continuous parameterization of world time and construction of time grid.

[0093] 1) Setting time nodes

[0094] Let the target time range of this invention be: ,in, The starting time; The end time is specified. Within this time range, a fixed time interval is selected. Set time nodes:

[0095] ,

[0096] In the formula, Let K be the time of the kth world time node; K is the total number of time nodes.

[0097] Let the increments of the world time remainder at all time points form a column vector:

[0098] ,

[0099] In the formula, For a moment The increment of the remaining amount in the world is unknown.

[0100] 2) Time interpolation and session-node relationship

[0101] At any time Within this invention, a piecewise linear function is used to approximate the world time residual increment:

[0102] ,

[0103] In the formula, The interpolation weights are defined as follows:

[0104] ,

[0105] Let the representative epoch of the s-th session be . And satisfy Then the world time margin increment for this session can be approximated as:

[0106] ,

[0107] In the formula, For session s, the interpolation coefficients relative to node k; Let be the interpolation coefficients of session s relative to node k+1.

[0108] Define the row vector corresponding to the time node of the s-th session. for:

[0109] ,

[0110] In the formula, for A 3D column vector, where all elements except the k-th and (k+1)-th elements are 0; the k-th element is... The (k+1)th element is Therefore, we can conclude that:

[0111] ,

[0112] Substituting equation (18) into the session reduction equation (9), we obtain the session time node parameters. Contributions:

[0113] ,

[0114] make The normal equation coefficient matrix represents the contribution of the s-th session to the time node parameters; The column vector of the constant terms of the normal equation contributing to the time node parameters for the s-th session can then be written as:

[0115] ,

[0116] Step 4: Fusion and joint solution of normal equations within the sliding time window.

[0117] 1) Construction of sliding time window

[0118] To improve numerical stability and facilitate time series smoothing, this invention constructs a sliding time window on the global time axis. Let: The length of the time window is ; the center time of the j-th time window is . The time interval covered by the j-th time window is ; For all representative eras The set of sessions that fall within this time interval.

[0119] 2) Superposition and weighting of normal equations within the time window

[0120] Within the j-th time window, this invention will include all elements belonging to the set We perform a weighted superposition of the conversational equations. Let: The weight factor for the s-th session (determined based on residual statistics); Let be the normal equation coefficient matrix synthesized within the j-th time window; Let be the column vector of constant terms of the normal equation synthesized within the j-th time window, then we have:

[0121] ,

[0122] To constrain the smoothness of time series and suppress ill-conditioned solutions, this invention modifies the time node parameters. A second-order difference smoothing constraint is applied. Let: The second-order difference operator matrix for time nodes; To smooth out the constraint weighting factor; This is the column vector of the incremental cosmic time balance of the time node to be estimated within the j-th time window.

[0123] The joint normal equation within the j-th time window is:

[0124] ,

[0125] Solving formula (22) yields the estimated time points within the time window:

[0126] ,

[0127] In the formula, This is the estimated vector of the world time surplus increment at time node j within the j-th time window. By sliding the time window along the time axis and repeating the above process, continuous time node estimates can be obtained over the entire target time range.

[0128] Step 5: Universal Time Continuous Sequence Recovery and Product Generation.

[0129] 1) Recovering the Universal Time Margin from Increments

[0130] After obtaining the time node increment estimate, this invention utilizes prior world time series to recover the world time residual. Let: Time node Estimation of the world time surplus increment (i.e. The k-th component); Time node Given the estimated world time surplus, we have:

[0131] ,

[0132] For any time A continuous world time residual estimate is obtained by using piecewise linear interpolation:

[0133] ,

[0134] 2) Recover the tidal term and generate a continuous sequence of Universal Time.

[0135] Finally, this invention recalculates the tidal correction at time t in universal time based on the adopted tidal model. Then add it back to the UT balance to obtain a continuous estimate of UT:

[0136] ,

[0137] In the formula, This is the world time at time t calculated by the present invention; This is the estimated world time margin at time t; Let t be the tidal correction at time t in the world.

[0138] Through the above steps, under the existing reference frame and unified constraint strategy, this invention jointly utilizes session-level normal equation information from enhanced and conventional observations to obtain a continuous Universal Time (UTC) solution sequence with high temporal resolution and high accuracy. It can also generate continuous UTC products at predetermined sampling intervals and perform corresponding accuracy assessments, as shown in Table 1 and... Figure 2 As shown.

[0139] Table 1. Weighted root mean square of world time residuals from weighted observations, conventional observations, and integrated solutions (unit: microseconds)

[0140] .

[0141] As shown in Table 1, enhanced observation is continuous but not very accurate, while conventional observation is accurate but not continuous. The solution of this invention can simultaneously meet the requirements of high accuracy and continuity.

[0142] The implementation of the various embodiments of the present invention is based on programmed processing through a device with processor functionality. Therefore, in practical engineering, the technical solutions and functions of the various embodiments of the present invention are encapsulated into various modules. Based on this reality, and building upon the above embodiments, the embodiments of the present invention provide a continuous world time (UTC) calculation apparatus, which is used to execute a UTC continuous calculation method from the above method embodiments.

[0143] The device includes: a first main module for acquiring and preprocessing the solution result files of each session within a predetermined time range; a second main module for applying uniform constraints to each preprocessed session and performing parameter block division and local parameter elimination; a third main module for constructing a sliding time window within a target time range, parameterizing Universal Time using a piecewise linear function within the window, and establishing a mapping relationship between sessions and time nodes; a fourth main module for constructing a sliding time window on a global time axis, collecting sessions within the window, superimposing the equations according to session weights, adding time node smoothing constraints, and jointly solving for the Universal Time residual increment; and a fifth main module for superimposing priors on the joint solution results, recovering the tidal term, and generating a continuous Universal Time sequence.

[0144] This invention provides a continuous Universal Time (UTC) solution device that addresses the problems of limited time resolution, insufficient utilization of session-level observation geometry and normal equation information, sensitivity to modeling differences between different session types, and difficulty in adapting to the current reference frame in existing UTC products based on daily parameter-level combinations. By employing the aforementioned modules, under the current reference frame and constraint strategy, it utilizes the normal equations of different session types in VLBI to achieve high time resolution, continuity, and high accuracy in UTC solution, thereby obtaining a continuous UTC sequence that balances timeliness and accuracy.

[0145] It should be noted that the device embodiments provided by the present invention are used not only to implement the methods in the above method embodiments, but also to implement the methods in other method embodiments provided by the present invention. The only difference is that corresponding functional modules are set. The principle is basically the same as that of the above device embodiments provided by the present invention. As long as those skilled in the art can improve the modules in the above device embodiments by referring to the specific technical solutions in other method embodiments and combining technical features to obtain corresponding technical means and technical solutions composed of these technical means, on the basis of the above device embodiments, and on the premise of ensuring the practicality of the technical solutions, they can obtain corresponding device-type embodiments for implementing the methods in other method-type embodiments.

[0146] Based on the same inventive concept as any of the foregoing embodiments, this embodiment of the invention also provides a continuous world time calculation device, including a memory and a processor. The memory stores program instructions that are executed by the processor, and the processor calls the program instructions to execute the continuous world time calculation method.

[0147] In embodiments of the present invention, the memory can be non-volatile memory, such as a hard disk drive (HDD) or a solid-state drive (SSD), or it can be volatile memory, such as random-access memory (RAM). Memory is any other medium capable of carrying or storing desired program code having an instruction or data structure form and accessible by a computer, but is not limited thereto. The memory in embodiments of the present invention can also be a circuit or any other device capable of implementing a storage function for storing program instructions and / or data.

[0148] In this embodiment of the invention, the processor may be a general-purpose processor, a digital signal processor, an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components, capable of implementing or executing the methods, steps, and logic block diagrams disclosed in this embodiment of the invention. The general-purpose processor may be a microprocessor or any conventional processor. The steps of the methods disclosed in this embodiment of the invention can be directly manifested as being executed by a hardware processor, or executed by a combination of hardware and software modules within the processor.

[0149] Based on the same inventive concept as any of the foregoing embodiments, this embodiment of the invention also provides a non-transitory computer-readable storage medium storing computer instructions that cause the computer to execute the following method for continuous world time calculation:

[0150] Obtain the solution result files for each session within a predetermined time range and perform preprocessing;

[0151] Apply uniform constraints to each session after preprocessing, and perform parameter block division and local parameter elimination;

[0152] Construct a sliding time window within the target time range, parameterize world time using a piecewise linear function within the window, and establish a mapping relationship between sessions and time nodes;

[0153] Construct a sliding time window on the global time axis, collect sessions within the window, superimpose the equations according to session weights, add time node smoothing constraints, and jointly solve the world time margin increment.

[0154] Priors are superimposed on the joint solution results, and the tidal term is recovered to generate a continuous sequence of universal time.

[0155] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the technical solutions of the embodiments of the present invention.

Claims

1. A method for continuous world time calculation, characterized in that, include: The solution result files of each session within a predetermined time range are obtained and preprocessed. Each session includes a first type of session and a second type of session. The first type of session is a short arc enhanced observation session, and the second type of session is a long arc conventional observation session. After preprocessing, uniform constraints are applied to each session, and parameter block division and local parameter elimination are performed to obtain the reduced equations for each session with respect to the world time margin increment only. Multiple time nodes with fixed intervals are set within the target time range, and the time nodes span the entire target time range. Within the window, a piecewise linear function is used to parameterize the world time surplus increment. By introducing representative epochs of each session, the world time surplus increment of the session is expressed as a linear combination of the unknown world time surplus increments at adjacent time nodes using piecewise linear interpolation. Substituting these values ​​into the reduction equation, the contribution of the session to the time node parameters is obtained, thereby establishing a mapping relationship between sessions and time nodes. A sliding time window is constructed on the global time axis, and sessions within the window are collected. The normal equations expressing the contribution of global time node parameters are weighted and superimposed according to the session weights. Second-order difference smoothing constraints are applied to the time node parameters to constrain the smoothness of the time series and suppress ill-conditioned solutions. The joint solution is used to obtain the estimated value of the world time residual increment at each time node. Using the prior world time series, the world time margin is recovered from the estimated world time margin increment value obtained from the solution, and the tidal term is also recovered, generating a continuous world time series that takes into account both time continuity and solution accuracy.

2. The method for continuous world time calculation according to claim 1, characterized in that, Obtain the solution result files for each session within a predetermined time range and perform preprocessing, including: Based on the solution result files of each VLBI session within a predetermined time range, the session normal equations are written as follows: ,in, For the first The normal equation coefficient matrix for each session; For the first A column vector of parameters to be estimated for each session; For the first A column vector of constant terms in the normal equations of each session; All station coordinates for all sessions are unified under the same Earth reference frame, and Universal Time tidal correction and incremental linearization are performed.

3. The method for continuous world time calculation according to claim 2, characterized in that, Perform Universal Time tidal correction and incremental linearization, including: Adjust the world time for tides to obtain the time. World Time Remaining for: ,in, For a moment World Time; For a moment World time tidal correction; Based on the prior value and observed value of the world time margin, the world time margin increment is obtained, and the world time margin increment is substituted into the conversational method equation as one of the world time unknowns.

4. The method for continuous world time calculation according to claim 3, characterized in that, The column vector of parameters to be estimated includes at least: the world time surplus increment, the clock error parameters of each VLBI station, and the tropospheric zenith delay parameters of each VLBI station.

5. The method for continuous world time calculation according to claim 1, characterized in that, Within the window, a piecewise linear function is used to parameterize world time, establishing a mapping relationship between sessions and time nodes, including: At any point between adjacent time nodes within the window, a piecewise linear function is used to approximate the increment of the world time balance; Introducing a representative epoch for the s-th session, where the representative epoch is located between the aforementioned adjacent time nodes, the world time residual increment is transformed using a piecewise linear function to obtain an approximate expression of the world time residual increment for the s-th session. By introducing the row vector corresponding to the time node of the s-th session, the approximate expression of the world time margin increment of the s-th session is transformed to obtain the final approximate expression of the world time margin increment. By substituting the final approximate expression of the world time margin increment into the reduction equation after local parameter elimination, a mapping relationship between sessions and time nodes is established.

6. The method for continuous world time calculation according to claim 1, characterized in that, The session weighted superposition method equations, combined with time node smoothing constraints, are used to jointly solve for the world time surplus increment, including: In the Within a time window, all items belonging to the set We perform a weighted superposition of the conversational equations to obtain the first... The normal equation coefficient matrix synthesized within each time window and the The column vector of constant terms of the normal equation synthesized within each time window ; Time node parameters A second-order difference smoothing constraint is applied, where, Let T be the second-order difference operator matrix at time node k, and let T be the transpose. To smooth out the constraint weighting factors, (j) For the first The column vector of the world time residual increment of the time node to be estimated within the n time windows, then the i-th The joint normal equations within each time window are: ， Solving the joint method equations, we obtain the first... Estimation of the increase in cosmic time residual at time points within a time window: ， By sliding the time window along the time axis and repeating the solution process, continuous time node estimates are obtained over the entire target time range.

7. The method for continuous world time calculation according to claim 6, characterized in that, After obtaining the estimate of the world time margin increment at the world node, the following is also included: By using prior world time series to recover the world time margin and combining it with piecewise linear interpolation, a continuous world time margin estimate is obtained. Based on the obtained Universal Time Residual Estimate, the tidal term is recovered, and a continuous Universal Time sequence is generated.

8. A continuous world time calculation apparatus for implementing the method according to any one of claims 1 to 7, characterized in that, include: The first main module is used to obtain the solution result files of each session within a predetermined time range and perform preprocessing. Each session includes a first type session and a second type session. The first type session is a short arc enhanced observation session, and the second type session is a long arc conventional observation session. The second main module is used to apply uniform constraints to each session after preprocessing, and to perform parameter block division and local parameter elimination. The third main module is used to construct a sliding time window within the target time range, parameterize world time using a piecewise linear function within the window, and establish a mapping relationship between sessions and time nodes. The fourth main module is used to construct a sliding time window on the global time axis, collect sessions within the window, superimpose the normal equations according to session weights, add time node smoothing constraints, and jointly solve the world time margin increment. The fifth main module is used to superimpose priors on the joint solution results and recover the tidal term to generate a continuous sequence of universal time.

9. A world time continuous calculation device, characterized in that, The device includes a memory and a processor, the memory storing program instructions that are executed by the processor, the processor invoking the program instructions to execute a continuous world time solution method as described in any one of claims 1 to 7.

10. A non-transitory computer-readable storage medium, characterized in that, The non-transitory computer-readable storage medium stores computer instructions that cause the computer to execute a continuous world time solution method as described in any one of claims 1 to 7.

Images