Method, device, equipment and medium for time difference calibration of low-orbit navigation enhancement transmitter

Abstract

The embodiment of the present application discloses a method, device, equipment and medium for calibrating the time difference of a low-orbit navigation enhancement transmitter, which relates to the field of satellite navigation technology. The method includes: calculating the low-orbit satellite onboard receiver clock error, ground station coordinates, station receiver clock error and other parameters through the observation data of the GNSS satellite; the ground station calculates the inter-system deviation between the low-orbit satellite system corresponding to each low-orbit satellite and the GNSS system based on the downlink navigation data of each low-orbit satellite; selecting a low-orbit satellite as a reference satellite, obtaining the inter-reference system deviation, calculating the navigation enhancement transmitter time deviation of other inter-system deviations and the inter-reference system deviation, and compensating the clock error of the corresponding low-orbit satellite based on the parameter. The present application can be used to calibrate the relative time difference of each low-orbit satellite navigation signal transmitter, realize the time synchronization of the low-orbit navigation signal, and is conducive to improving the positioning, navigation and timing capabilities of the low-orbit satellite system, and providing more accurate services to users.

Description

Technical Field

[0001] The present application relates to the field of satellite navigation technology, and in particular to a method, device, equipment and medium for calibrating time difference of a low-orbit navigation enhancement transmitter. Background Art

[0002] Because the Global Navigation Satellite System (GNSS) uses satellites in medium- to high-Earth orbits, they suffer from shortcomings such as slow carrier phase positioning convergence, weak signal strength, and susceptibility to interference. In contrast, low-Earth orbit (LEO) satellites have lower orbital altitudes, faster platform motion, less signal loss, and significantly higher signal landing power than GNSS satellites. Therefore, incorporating LEO navigation augmentation systems (LEO navigation augmentation systems) from LEO satellites into a comprehensive PNT architecture can improve positioning effectiveness.

[0003] During the positioning process, inter-satellite time synchronization is a key prerequisite for satellite navigation systems. GNSS satellites can achieve highly accurate clock error predictions through high-precision and high-stability atomic clocks. However, for low-orbit satellites, due to the limitations of satellite platform cost, size, weight and power consumption, low-orbit satellites use low-cost ultra-stable oscillators (USOs) instead of atomic clocks.

[0004] Currently, to synchronize the onboard time of low-orbit satellites with the GNSS timescale and avoid clock errors from low-cost, ultra-stable oscillators, the onboard GNSS receivers and the onboard navigation direct transmission payloads of low-orbit satellites are generally traced to the same onboard time frequency and driven by onboard GNSS technology. In this scheme, low-orbit satellites receive signals from GNSS multi-system navigation satellites in medium and high orbits, and calculate their own positions and receiver clock errors in real time to achieve inter-satellite time synchronization. At the same time, the onboard GNSS receiver's time frequency is also used to synchronously drive the navigation direct transmission payload, allowing the navigation direct transmission signal to be processed directly based on the onboard GNSS receiver clock error. By fitting the real-time onboard GNSS receiver clock error, low-orbit clock error prediction is performed, thereby providing PNT services and enhanced services to ground users.

[0005] However, this method still has some defects. Since there is a hardware link delay between the onboard navigation signal transmitter and the onboard receiver of each low-orbit satellite, and the hardware link delay of each low-orbit satellite may be different, when using the low-orbit downlink signal for navigation and positioning solution, the systematic ranging error between the observation values ​​of each low-orbit satellite cannot be ruled out, thereby affecting the user's positioning and timing accuracy.

[0006] Therefore, there is currently a lack of a method that can not only improve the positioning effect of the PNT system by utilizing the low-orbit navigation augmentation system but also eliminate the above-mentioned systematic ranging errors. Summary of the Invention

[0007] The present application provides a method, apparatus, device, and medium for calibrating the time difference of a low-orbit navigation enhancement transmitter, to address the deficiencies in the above-mentioned related technologies. The technical solution is as follows:

[0008] In a first aspect, an embodiment of the present application provides a low-orbit navigation enhancement transmitter time difference calibration method, comprising:

[0009] Acquiring first observation data of the GNSS satellite through an onboard receiver of the low-orbit satellite, and calculating a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data;

[0010] Acquire secondary observation data of the GNSS satellite through a ground station receiver, and obtain the station receiver coordinates, station receiver clock error, and tropospheric parameters through ground precise point positioning based on the secondary observation data;

[0011] Obtain downlink navigation data sent by each low-orbit satellite through the ground station, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite by combining the onboard receiver clock error, the low-orbit precise orbit, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data;

[0012] Selecting any low-orbit satellite as a reference satellite, obtaining the reference inter-system bias of the reference satellite, calculating the difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and outputting the difference as the navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite;

[0013] Based on the navigation enhancement transmitter time offset of each low-orbit satellite, the corresponding low-orbit satellite's onboard receiver clock error is compensated, and the low-orbit satellite navigation signal transmitter clock error with a unified time reference is output.

[0014] In an optional solution of the first aspect, the satellite-borne receiver clock error is calculated based on the following steps, including:

[0015] Establishing an observation equation based on the first observation data to obtain a pseudorange observation equation and a carrier phase observation equation;

[0016] Based on the pseudorange observation equation and the carrier phase observation equation, the frequency-dependent characteristics of ionospheric delay are utilized to construct a combined observation value without ionospheric delay, and the first-order term of ionospheric delay is eliminated to obtain a pseudorange observation equation without ionospheric delay and a carrier phase observation equation without ionospheric delay;

[0017] Solve the pseudorange observation equation without ionospheric delay and the carrier phase observation equation without ionospheric delay, and output the satellite receiver clock error.

[0018] In an optional solution of the first aspect, the low-orbit precise orbit is calculated based on the following steps, including:

[0019] Establishing a motion equation of the low-orbit satellite based on the first observation data;

[0020] The state transfer matrix and the orbital state information of the low-orbit satellite at each moment are calculated based on the numerical integration method;

[0021] The orbital parameters of the low-orbit satellite are calculated based on a parameter estimation method in combination with the state transfer matrix and the orbital state information of the low-orbit satellite at each moment, and the low-orbit precise orbit is output.

[0022] In an optional solution of the first aspect, obtaining the station receiver coordinates, the station receiver clock error, and the tropospheric parameters by ground precise point positioning based on the second observation data includes:

[0023] Establishing an observation equation based on the second observation data to obtain a pseudorange observation equation and a carrier phase observation equation of the ground station;

[0024] Based on the pseudorange observation equation of the ground station and the carrier phase observation equation of the ground station, utilizing the frequency-dependent characteristics of ionospheric delay, constructing a combined observation value without ionospheric delay, eliminating the first-order term of ionospheric delay, and constructing a pseudorange observation equation without ionospheric delay and a carrier phase observation equation without ionospheric delay;

[0025] The pseudorange observation equation without ionospheric delay and the carrier phase observation equation without ionospheric delay corresponding to the ground station are solved by ground precise single point positioning, and the station receiver coordinates, the station receiver clock error and the tropospheric parameters are output.

[0026] In an optional solution of the first aspect, the calculating, by combining the onboard receiver clock error, the low-orbit precise orbit of the low-orbit satellite, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data, an inter-system bias between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite includes:

[0027] Constructing a ground station low-orbit observation equation that takes into account the hardware delay of the low-orbit satellite, including:

[0028]

[0029] The low-orbit satellite includes an onboard receiver for receiving GNSS signals and a navigation signal transmitter for sending downlink navigation data. The hardware delay is the clock difference between the onboard receiver clock error and the navigation signal transmitter clock error:

[0030]

[0031] Each time, the downlink navigation data of a low-orbit satellite is selected and substituted into the low-orbit observation equation of the ground station. Other bias terms and hardware delay are estimated as a whole to obtain the inter-system bias. The formula is applied:

[0032]

[0033] Among them, i is the low-orbit satellite number, G is the GNSS system, corresponds to the low-orbit satellite L i The satellite receiver clock error is corresponds to the low-orbit satellite L i The navigation signal transmitter clock error, ISB LEO,r is the other deviation term, corresponds to the low-orbit satellite L i Hardware latency, represents the geometric distance from the ground station r to the low-orbit satellite i, is the station receiver clock error, T r is the tropospheric parameter, is the fuzziness parameter, corresponds to the low-orbit satellite L i The inter-system deviation, ε P,IF is the measurement noise term of the pseudorange observation, ε l,IF is the measurement noise term of the carrier phase observation.

[0034] In an optional solution of the first aspect, the satellite receiver clock error is the error term due to relativistic effects and other clock error terms δt L The sum of:

[0035]

[0036] The navigation signal transmitting end clock error is the error term due to relativistic effects and other clock errors The sum of:

[0037]

[0038] The low-orbit satellite includes an onboard receiver for receiving GNSS signals and a navigation signal transmitter for sending downlink navigation data. The hardware delay is the clock error of the onboard receiver. Clock difference with the navigation signal transmitter Clock difference Application formula:

[0039]

[0040]

[0041] Eliminate the error term of low-orbit satellites affected by relativistic effects Application formula:

[0042] Δt L,ST =δt L -δt L .

[0043] In an optional solution of the first aspect, compensating the clock error of the corresponding low-orbit satellite based on the navigation augmentation transmitter timing offset of each low-orbit satellite includes:

[0044] Obtain the navigation augmentation transmitter timing offset of each low-orbit satellite respectively, and compensate the navigation augmentation transmitter timing offset to the onboard receiver clock error of the corresponding low-orbit satellite using the formula:

[0045]

[0046] Output the compensated low-orbit navigation signal clock error sequence of each low-orbit satellite

[0047] Among them, i is the low-orbit satellite number, corresponds to the low-orbit satellite L i The clock error of the navigation signal transmitter is corresponds to the low-orbit satellite L i The satellite receiver clock error before compensation.

[0048] In a second aspect, an embodiment of the present application further provides a low-orbit navigation enhancement transmitter time difference calibration device, comprising:

[0049] a low-orbit satellite-borne data acquisition unit, configured to acquire first observation data of a GNSS satellite through an onboard receiver of the low-orbit satellite, and calculate a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data;

[0050] a ground station data acquisition unit, configured to acquire second observation data of the GNSS satellite through a station receiver at the ground station, and obtain station receiver coordinates, station receiver clock errors, and tropospheric parameters through ground precise point positioning based on the second observation data;

[0051] a data processing unit, configured to obtain the downlink navigation data sent by each low-orbit satellite through the ground station data acquisition unit, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite by combining the onboard receiver clock error, the low-orbit precise orbit of the low-orbit satellite, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data;

[0052] The data processing unit is further configured to select any low-orbit satellite as a reference satellite, obtain a reference inter-system bias of the reference satellite, calculate a difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and output the difference as a navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite;

[0053] The data processing unit is further configured to compensate for the clock error of the corresponding low-orbit satellite based on the navigation enhancement transmitting end time offset of each low-orbit satellite.

[0054] In a third aspect, an embodiment of the present application further provides an electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the program, the method provided in the first aspect of the embodiment of the present application or any one of the implementations of the first aspect is implemented.

[0055] In a fourth aspect, the present application also provides a non-transitory computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, it implements the method provided by the first aspect of the embodiment of the present application or any one of the implementation methods of the first aspect.

[0056] The beneficial effects of the technical solutions provided by some embodiments of the present application include at least:

[0057] An embodiment of the present application provides a method for calibrating the time difference of a low-orbit navigation enhancement transmitter. Downlink navigation data of a low-orbit satellite is obtained through a ground station, thereby calculating the inter-system deviation between the low-orbit satellite system and the GNSS system at the station receiving end of the ground station. The obtained inter-system deviation can implicitly eliminate the influence of the relativistic effect received by the low-orbit satellite.

[0058] Furthermore, by selecting a reference star, the reference system deviation of the reference star is subtracted from the system deviation of other low-orbit satellites, and the navigation enhancement transmitter time deviation obtained by the difference is compensated to the clock error of the corresponding low-orbit satellite, thereby obtaining the clock error of the navigation signal transmitter of the accurately synchronized low-orbit satellite. The time synchronization calibration from the ground station receiver to the low-orbit satellite navigation signal transmitter can be completed using only the data from the ground station, avoiding the defect of different receiving and transmitting delays caused by hardware delays of different low-orbit satellites, and greatly simplifying the calculation model, which can improve the solution efficiency. Low-orbit satellites that compensate for hardware delays can provide users with higher-precision positioning, navigation and timing services. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] In order to more clearly illustrate the technical solutions in this application or related technologies, the following briefly introduces the drawings required for use in the embodiments or related technical descriptions. Obviously, the drawings described below are some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0060] Figure 1 This is a schematic diagram of an application scenario of a time difference calibration method for a low-orbit navigation enhancement transmitter provided in an embodiment of the present application;

[0061] Figure 2 This is a flow chart of a method for calibrating the time difference of a low-orbit navigation enhancement transmitter provided in an embodiment of the present application;

[0062] Figure 3 This is a structural diagram of a time difference calibration device for a low-orbit navigation enhancement transmitter provided in an embodiment of the present application;

[0063] Figure 4 It is a structural diagram of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0064] To make the objectives, technical solutions, and advantages of this application more clear, the technical solutions of this application will be clearly and completely described below in conjunction with the accompanying drawings. Obviously, the embodiments described are only part of the embodiments of this application, not all of them. Based on the embodiments of this application, all other embodiments obtained by ordinary technicians in this field without making any creative efforts are within the scope of protection of this application.

[0065] The terms "including" and "having," and any variations thereof, in the specification and claims of this application and the accompanying drawings are intended to cover non-exclusive inclusions. For example, a process, method, system, product, or apparatus comprising a series of steps or modules is not limited to the listed steps or modules, but may optionally include steps or modules not listed, or may optionally include other steps or modules inherent to the process, method, product, or apparatus.

[0066] It should be noted that the terms "first" and "second" used in this application are merely used to distinguish similar objects and do not represent a specific ordering of the objects. It is understood that the terms "first" and "second" may interchangeably represent a specific order or precedence, where permitted. It should be understood that the objects distinguished by "first" and "second" may interchangeably represent a specific order or precedence, where appropriate, such that the embodiments of the present application described herein can be implemented in an order other than that described or illustrated herein.

[0067] Next, combine Figure 1 , introduces the application scenario of a low-orbit navigation enhanced transmitter time difference calibration method provided in an embodiment of the present application.

[0068] This application scenario includes a GNSS satellite 100, a low-orbit satellite 110, and a ground station 120. GNSS satellite 100 is used to transmit GNSS navigation signals, which can be acquired by both low-orbit satellite 110 and ground station 120. This includes observation data from low-orbit satellite 110 on GNSS satellites and observation data from ground station 120 on GNSS satellites. Low-orbit satellite 110 can enhance navigation, receiving GNSS navigation signals transmitted by GNSS satellite 100 through its onboard receiver and transmitting navigation enhancement signals to ground station 120 through its onboard navigation signal transmitter. Ground station 120 can receive both GNSS navigation signals transmitted by GNSS satellite 100 and navigation enhancement signals transmitted by the navigation signal transmitter of low-orbit satellite 110.

[0069] It should be noted that Figure 1 The GNSS satellites 100 shown may include a plurality of different GNSS satellites, not limited to Figure 1 GNSS satellites 101, 102, 103, 104, and 105 are shown. In actual applications, GNSS satellite 100 can be selected from one of the GNSS systems such as GPS, Beidou Navigation, and Galileo, or can be selected from a combination of multiple GNSS systems. This embodiment of the present application is not limited to this. Figure 1 The low-orbit satellite 110 shown can actually include multiple different low-orbit satellites, not limited to Figure 1 The low-orbit satellite 111 and the low-orbit satellite 112 are shown; similarly, Figure 1 Only one ground station 120 is shown. In actual application, multiple ground stations may be included. The ground stations may include base stations, GNSS receivers, and terminal devices with the function of receiving GNSS signals. The embodiments of the present application are not limited to this. In some application scenarios, the mobile terminal held by the user can also be regarded as a ground station.

[0070] The present application is described in detail below with reference to specific embodiments.

[0071] Next, combine Figure 2 , introduces a time difference calibration method for the low-orbit navigation enhancement transmitter provided by the embodiment of this application. For details, please refer to Figure 2 , Figure 2 The figure shows a flow chart of a method for calibrating the time difference of a low-orbit navigation enhancement transmitter provided by an embodiment of the present application. Figure 2 As shown, the method includes the following steps:

[0072] S201: Acquire first observation data of a GNSS satellite through an onboard receiver of a low-orbit satellite, and calculate a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data.

[0073] S202: Acquire, by a station receiver at a ground station, second observation data of the GNSS satellite, and obtain station receiver coordinates, station receiver clock errors, and tropospheric parameters by ground precise point positioning based on the second observation data.

[0074] S203: Obtain downlink navigation data sent by each LEO satellite through the ground station, and calculate the inter-system deviation between the LEO satellite system and the GNSS system corresponding to each LEO satellite by combining the onboard receiver clock error, the LEO precise orbit of the LEO satellite, the station receiver coordinates, the station receiver clock error, tropospheric parameters, and the downlink navigation data;

[0075] S204: Select any low-orbit satellite as a reference satellite, obtain the reference inter-system bias of the reference satellite, calculate the difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and output the difference as the navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite;

[0076] S205 , based on the navigation enhancement transmitter time offset of each low-orbit satellite, the onboard receiver clock error of the corresponding low-orbit satellite is compensated, and the low-orbit satellite navigation signal transmitter clock error with a unified time reference is output.

[0077] Specifically, S201 includes the following sub-steps:

[0078] In step S2011, a GNSS observation data, namely, first observation data, is received by an onboard receiver of a low-orbit satellite, and an observation equation is established based on the first observation data:

[0079] For example, pseudorange observation values ​​have the following pseudorange observation equation:

[0080]

[0081] Wherein, the superscript G indicates that the GNSS satellite is the transmitter, the subscript L indicates that the low Earth orbit satellite (LEO) is the receiver, that is, the low Earth orbit satellite receives the GNSS signal transmitted by the GNSS satellite through the onboard receiver, the subscript j indicates the observation value frequency; P indicates the pseudorange observation value, the unit is meter; Represents the geometric distance between the low-orbit satellite and the GNSS satellite. If the coordinates of the low-orbit satellite are marked as (x L ,y L , z L ), the coordinates of the GNSS satellite are marked as (x G ,y G , z G ), then:

[0082]

[0083] Among them, δt L is the receiver clock error of the onboard receiver of the low-orbit satellite; δt G is the clock error of the GNSS satellite. Both are defined on the gravity potential W0 at the mean sea level. The relativistic effect of the low-orbit satellite is The relativistic effect of GNSS satellites is represents the ionospheric error at frequency j; d L,j is the pseudo-range hardware delay of the onboard receiver of the low-orbit satellite representing the channel correlation, is the pseudorange hardware delay of the channel-related GNSS satellite; ε P,j represents the unmodeled errors such as pseudorange observation noise and multipath error.

[0084] For carrier observations, there are the following observation equations:

[0085]

[0086] Where, l represents the carrier phase observation value, in meters; b L,j is the channel-related phase delay of the onboard receiver of the low-orbit satellite, is the phase delay of the GNSS satellite associated with the channel; is the carrier phase integer ambiguity; λ j is the wavelength of frequency j; εl,j represents the carrier observation noise.

[0087] It should be noted that the orbital parameters and clock errors of GNSS satellites are known values, which are then used to estimate the orbital parameters and clock errors of low-orbit satellites.

[0088] S2012, using the frequency-dependent characteristics of ionospheric delay, constructs a combined observation value without ionospheric delay, eliminates the influence of the first-order ionospheric delay, reduces unknown parameters, and only estimates the low-orbit orbit, satellite clock error, and ambiguity parameters. The specific combined model constructed is:

[0089]

[0090] in,

[0091]

[0092] λ IF =A·λ j1 +B·λ j2 ;

[0093]

[0094] Among them, M P is the error in the pseudorange observation that can be modeled, M L is the error in the carrier phase observation that can be modeled; IF means ionosphere-free.

[0095] Among them, the tropospheric error can be corrected by the Saastamoinen model and the GMF projection function model; the GNSS satellite antenna phase center deviation PCO and the antenna phase center variation PCV can be corrected by the antenna products released by IGS; the phase winding error δφ on the GNSS satellite side can be corrected by the following model:

[0096]

[0097] where k is the unit vector from the GNSS satellite to the LEO satellite-borne receiver; D′ and D are the effective dipole vectors of the GNSS satellite and onboard receiver antennas, respectively, which are calculated from the unit vectors (x′, y′, z′) in the satellite-fixed system and the unit vector (x, y, z) in the Earth-fixed system.

[0098] It should be noted that for GNSS systems, the correction method for relativistic effects can apply a unified protocol specification, and only correct its periodic variation term, using the formula:

[0099]

[0100] Among them, RG is the position vector of the GNSS satellite, V G is the velocity vector of the GNSS satellite, and c represents the speed of light.

[0101] It should be noted that in actual processing, as long as the relativistic effect correction model used by the server for processing GNSS data and the user end using GNSS data is consistent, the relativistic effect of GNSS satellites can be basically eliminated. The influence of the low-orbit satellite, and the low-orbit satellite is different from the GNSS satellite. The low-orbit satellite is in a lower orbit and the gravity potential changes faster. The relativistic effect is more significantly affected by the high-order gravity field and the tide than the GNSS satellite. The impact on the low-orbit satellite clock is more significant, and the correction method is more complicated. Therefore, the above method cannot eliminate the influence of the relativistic effect of the low-orbit satellite. The relativistic effect on the low-orbit satellite Requires special consideration.

[0102] Specifically, the onboard receiver clock error of the low-orbit satellite has absorbed the influence of the low-orbit relativistic effect. It can be understood that the actual onboard receiver clock error is The error term of the relativistic effect is included in Should be the error term of relativistic effect and other clock error terms δt L The combination can be expressed as:

[0103]

[0104] Based on the steps of S2011-S2012, the total clock error including the relativistic effect can be obtained. That is, the satellite receiver clock error finally obtained in S201

[0105] Further, the satellite orbit of the low-orbit satellite is determined by the following steps:

[0106] S2013, Establishing the equations of motion for low-orbit satellites.

[0107] It should be noted that the satellite orbit of a low-orbit satellite is a series of state quantities characterized by initial orbital parameters and a dynamic model, which can be specifically expressed as:

[0108]

[0109] The differential equation is as follows:

[0110]

[0111] Among them, function F is an n-dimensional nonlinear functional, and x0 is the initial state.

[0112] Let X = xx * Then we have:

[0113]

[0114] in, Assume that the solution of the above formula is:

[0115] X=ψ(t,t0)X0;

[0116] Then we have:

[0117]

[0118] Where I is the identity matrix and ψ(t, t0) is the state transition matrix, which is specifically expressed as:

[0119]

[0120] Finally, the state transfer matrix and the orbital state information of the low-orbit satellite at each moment are calculated through numerical integration methods such as RKF and ADAMS, which can be used to linearize the subsequent observation equations and solve the satellite orbit of the low-orbit satellite.

[0121] S2014, parameter estimation.

[0122] Optionally, the parameter estimation method includes least squares estimation, Kalman filtering, etc.

[0123] Based on the state transition matrix calculated in S2013 and the orbital state information of the low-orbit satellite at each moment, the following formula is applied:

[0124]

[0125] Where x * (t) is the reference orbital parameter obtained by integration under the initial conditions. * (t) is x * (t) is the observed value calculated. * (t) and y * Expanding at (t) and taking the first term yields:

[0126]

[0127] in:

[0128]

[0129] Solving the differential equation yields:

[0130] X=ψ(t,t0)X0;

[0131] Among them, the transfer matrix ψ(t,t0) satisfies:

[0132]

[0133] The solution of the transfer matrix is ​​realized by orbital integration, which can be obtained:

[0134] Y i =H i ψ(t i ,t0)X0+V i ;

[0135] is the partial derivative of the observed quantity with respect to each parameter to be estimated, based on the above formula:

[0136] L=BX0+V;

[0137] Subsequently, the satellite orbit parameters of the low-orbit satellite can be obtained by least squares estimation or filtering estimation methods.

[0138] Based on this, through the above steps S2011-S2012, the onboard receiver clock error of each low-orbit satellite can be determined. Through the above steps S2013-S2014, the orbital parameters of each low-orbit satellite can be determined, that is, the satellite orbit of the low-orbit satellite can be determined.

[0139] In some embodiments, in S202, each ground station receives a GNSS signal transmitted by a GNSS satellite through a station receiver, that is, obtains second observation data of the GNSS satellite through the station receiver, and obtains ground precise point positioning parameters through ground precise point positioning solution, which specifically includes the following steps:

[0140] S2021: Establish observation equations based on GNSS observation data obtained by the station receiver.

[0141] The pseudorange observation equation is as follows:

[0142]

[0143] Where r is the number of the ground station, G is the number of the GNSS satellite, and P represents the pseudorange observation value in meters; Represents the geometric distance from the ground station r to the GNSS satellite G. If the position coordinates of the ground station are (x r ,y r , z r ), the position coordinates of the GNSS satellite are (x G ,y G , z G ), then is the station receiver clock error of the ground station, specifically the clock error of the ground station relative to the GNSS reference system, δt G represents the clock error of GNSS satellite; T r is the tropospheric delay error; represents the ionospheric delay error at frequency j; is the hardware delay of the ground station receiver at frequency j, is the hardware delay of the GNSS satellite at frequency j; ε P,j represents the measurement noise term of the pseudorange observation.

[0144] For carrier observations, there are the following observation equations:

[0145]

[0146] Where, l represents the carrier phase observation value in meters. is the wavelength of the signal with frequency number j transmitted by the GNSS system satellite. For a given system and signal frequency, this value is a constant. is the carrier phase hardware delay corresponding to the signal of frequency j at the ground station receiver, is the carrier phase hardware delay corresponding to the signal with frequency j at the GNSS satellite end, and the unit is week; is the integer ambiguity of the carrier phase observation value, in weeks; ε l,j represents the measurement noise term of the carrier phase observation.

[0147] The hardware delay of pseudorange and carrier phase is often absorbed into the clock error parameter, so the observation equation can be further written as follows:

[0148]

[0149] in, is the ground station receiver clock error, δt G It is the clock difference including hardware delay of GNSS satellite.

[0150] S2022, construct ionospheric-free delay combined observations.

[0151] By utilizing the frequency-dependent characteristics of ionospheric delay, a combined observation value without ionospheric delay is constructed to eliminate the influence of the first-order term of ionospheric delay and reduce unknown parameters. The specific combined model is:

[0152]

[0153] S2023, perform parameter estimation.

[0154] Optionally, parameter estimation can be performed by least square estimation, Kalman filtering, etc. During the calculation process, filtering can be performed by a suitable stochastic model.

[0155] Specifically, the ionospheric-free delay combination model of S2022 is a linear combination of the original observations, and the initial variance can be calculated according to the error propagation law.

[0156] Specifically, weighting can be performed using the satellite elevation angle, wherein the ground station receiver can be considered to be stationary. Strong constraints can be placed on the position coordinates of the station receiver during filtering, and the station receiver can be considered to be static.

[0157] Specifically, after parameter estimation, the ground precise point positioning parameters are output. The ground precise point positioning parameters calculated based on S2021-S2023 include the coordinates of the station receiver, the station receiver clock error, and the ground station receiver clock error. and the tropospheric parameter T r .

[0158] It should be noted that the embodiment of the present application does not limit the execution order of S201 and the corresponding sub-steps S2011-S2014 and S202 and the corresponding sub-steps S2021-S2023. S201 and the corresponding sub-steps S2011-S2014 can be executed first, and then S202 and the corresponding sub-steps S2021-S2023 can be executed; or S202 and the corresponding sub-steps S2021-S2023 can be executed first, and then S201 and the corresponding sub-steps S2011-S2014 can be executed; the two can also be executed simultaneously.

[0159] Furthermore, S201 calculates the satellite orbit of the low-orbit satellite and the onboard receiver clock error of the low-orbit satellite based on the observation data of the GNSS satellite by the receiver of the low-orbit satellite, and S202 calculates the ground precise single-point positioning parameters obtained after performing ground precise single-point positioning based on the observation data of the GNSS satellite by the ground station, including: station receiver coordinates, station receiver clock error and tropospheric parameters, and then executes step S203.

[0160] Specifically, in S203, the ground station can obtain the downlink navigation data sent by each low-orbit satellite, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite based on the downlink navigation data of each low-orbit satellite.

[0161] Satellite receiver clock error calculated based on S201 S202 calculates the station receiver clock error and the tropospheric parameter T r, the ground station low-orbit observation equation considering the time difference of the low-orbit navigation enhancement transmitter is constructed, specifically:

[0162]

[0163] Among them, i is the low-orbit satellite number, G is the GNSS system, corresponds to the low-orbit satellite L i The onboard receiver clock error, ISB LEO,r represents other deviation terms, corresponds to the low-orbit satellite L i Low Earth Orbit satellite hardware latency; Represents the geometric distance from the ground station r to the low-orbit satellite i. If the position coordinates of the ground station are (x r ,y r , z r ), low-orbit satellite L i The position coordinates are Then there are:

[0164]

[0165] Specifically, the above equation can be solved by least squares or Kalman filtering methods.

[0166] Among the known parameters, the station receiver clock error and the tropospheric parameter T r It has been determined in step S202; low-orbit satellite L i The onboard receiver clock error In step S201, it is determined that

[0167] Thus, the unknown parameters include other bias terms ISB LEO,r , low-orbit satellite hardware delay

[0168] It should be noted that since only one low-orbit satellite data is used each time, it is impossible to distinguish other bias terms and low-orbit satellite hardware delays. LEO,r , low-orbit satellite hardware delay That is, the inter-system deviation calculated by the above-mentioned ground station low-orbit observation equation based on the downlink navigation data of each low-orbit satellite in S203 is actually the sum of other deviation terms and low-orbit satellite hardware delay, thereby calculating each low-orbit satellite L i The corresponding inter-system deviation between the low-orbit satellite system and the GNSS system is recorded as As shown in the following formula:

[0169]

[0170] It should be noted that low-orbit satellites not only have onboard receivers for receiving GNSS data, but also have navigation signal transmitters for sending downlink navigation data. Ground stations can obtain the downlink navigation data sent by low-orbit satellites, that is, the ground stations' observation data of low-orbit satellites. They can also obtain GNSS navigation data sent by GNSS satellites, that is, the ground stations' observation data of GNSS satellites. It is understandable that there is a hardware delay between each low-orbit satellite's onboard receiver for receiving GNSS signals and the navigation signal transmitter for sending downlink navigation data. Different low-orbit satellites have different hardware delays. The above formula That is, each low-orbit satellite L i The corresponding hardware delay, that is, the clock difference between the onboard receiver and the navigation signal transmitter, is as follows:

[0171]

[0172] in, is the clock error of the navigation signal transmitter, is the satellite receiver clock error.

[0173] It should be noted that the station coordinates, station receiver clock errors, and troposphere parameters calculated based on the GNSS system can be regarded as fixed values, and only the observation data of the low-orbit satellite tracked by the ground station at a certain moment are selected for calculation, and the deviation value due to the ground station directly receiving the GNSS satellite signal and receiving the downlink navigation data sent by the low-orbit satellite (that is, the other deviation term ISB in the above formula) is taken into account. LEO,r Considered as a constant.

[0174] It is understandable that the other bias terms ISB for each LEO satellite LEO,r is a satellite-independent term, so different low-orbit satellites can be considered to have the same other bias terms ISB LEO,r .

[0175] Furthermore, in S204, any low-orbit satellite can be selected as a reference star to obtain the inter-reference system deviation of this reference star, which is recorded as Calculate other low-orbit satellites L i Inter-system deviation Deviation from the reference system The difference is the low-orbit satellite L i The navigation enhancement transmitter time deviation.

[0176] Specifically, the formula is applied:

[0177]

[0178] It can be seen from the above formula that other deviation terms ISB LEO,rIt is a satellite-independent term and is eliminated in the process of calculating the difference, leaving only the low-orbit satellite L i Hardware latency Hardware delay with low-orbit satellite L0 The relative difference of the low-orbit satellite L i Navigation enhancement transmitter time offset

[0179] Furthermore, in S205, the low-orbit satellite L calculated in S204 can be i Navigation enhancement transmitter time offset Compensation to low-orbit satellite L i The formula is applied to the satellite receiver clock error:

[0180]

[0181] This ensures that the clock differences of the satellite navigation signal transmitters have the same time base after compensation.

[0182] It should be noted that when navigation positioning enhancement is performed based on the compensated clock error of the navigation signal transmitter, the influence of hardware delay between the onboard receiver and the navigation signal transmitter of the low-orbit satellite can be avoided, which is conducive to improving the positioning accuracy of low-orbit navigation enhancement through low-orbit satellites.

[0183] In some embodiments, each low-orbit satellite L i Navigation enhancement transmitter time offset Compensation is performed to obtain the low-orbit navigation signal clock error sequence with time delay correction between the onboard receiver and the navigation signal transmitter.

[0184] In this way, the impact of hardware delays of multiple low-orbit satellites can be eliminated, allowing users to combine multiple low-orbit satellites for navigation and positioning enhancement, solving the problem of relative time deviation between navigation enhancement signals broadcast by different low-orbit satellites.

[0185] In some embodiments, taking into account the influence of relativistic effects, the actual onboard receiver clock error is Should be the error term of relativistic effect and other clock error terms δt L In the process of obtaining the navigation enhancement transmitter time deviation by difference in S203, the error term of the relativistic effect will also be eliminated, as shown below:

[0186]

[0187] Eliminate the error term of low-orbit satellites affected by relativistic effects Application formula:

[0188] ΔtL,ST =δt L -δt L .

[0189] Therefore, when calculating the inter-system deviation between the LEO satellite system and the GNSS system corresponding to each LEO satellite in step S203, the relativistic effect be eliminated.

[0190] In this way, the method provided in the embodiment of the present application can implicitly correct the low-orbit relativistic effects, is highly efficient and simple to calculate, and can provide a fast, simple and efficient relativistic effect correction method, which is conducive to improving the efficiency of low-orbit satellites in providing navigation enhancement services.

[0191] The following are device embodiments of the present application, which can be used to implement the method embodiments of the present application. For details not disclosed in the device embodiments of the present application, please refer to the method embodiments of the present application.

[0192] See next Figure 3 , is a structural diagram of a low-orbit navigation enhancement transmitter time difference calibration device provided by an exemplary embodiment of the present application. The device can be implemented as all or part of a terminal through software, hardware, or a combination of both, and can also be integrated on a server as an independent module. A low-orbit navigation enhancement transmitter time difference calibration device in an embodiment of the present application can be applied to a terminal or the cloud. The device 30 includes a low-orbit satellite-borne data acquisition unit 301, a ground station data acquisition unit 302, and a data processing unit 303, wherein:

[0193] The low-orbit satellite-borne data acquisition unit 301 is configured to acquire first observation data of a GNSS satellite through an onboard receiver of a low-orbit satellite, and calculate a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data;

[0194] The ground station data acquisition unit 302 acquires second observation data of the GNSS satellite through a station receiver of the ground station, and obtains station receiver coordinates, station receiver clock errors, and tropospheric parameters through ground precise point positioning based on the second observation data;

[0195] The data processing unit 303 is configured to obtain the downlink navigation data sent by each low-orbit satellite through the ground station data acquisition unit 302, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite by combining the onboard receiver clock error, the low-orbit precise orbit of the low-orbit satellite, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data;

[0196] The data processing unit 303 is further configured to select any low-orbit satellite as a reference satellite, obtain a reference inter-system bias of the reference satellite, calculate the difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and output the difference as a navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite;

[0197] The data processing unit is further configured to compensate for the clock error of the onboard receiver of the corresponding low-orbit satellite based on the navigation enhancement transmitter time offset of each low-orbit satellite, and output the low-orbit satellite navigation signal transmitter clock error with a unified time reference.

[0198] It should be noted that the device 30 provided in the above embodiment, when performing a method for calibrating the time difference between a transmitter and a low-orbit navigation enhancement device, only uses the division of the above-mentioned functional modules as an example. In actual applications, the above-mentioned functions can be assigned to different functional modules as needed, that is, the internal structure of the device can be divided into different functional modules to complete all or part of the functions described above. In addition, the device provided in the above embodiment and the embodiment of the method for calibrating the time difference between a transmitter and a low-orbit navigation enhancement device are based on the same concept. The implementation process is detailed in the method embodiment and will not be repeated here.

[0199] An embodiment of the present application also provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the program, the steps of the method of any of the above embodiments are implemented.

[0200] See Figure 4 , is a structural block diagram of an electronic device provided in an embodiment of the present application.

[0201] like Figure 4 As shown, the electronic device 400 includes a processor 401 and a memory 402 .

[0202] In the embodiment of the present application, the processor 401 is the control center of the computer system and can be the processor of a physical machine or the processor of a virtual machine. The processor 401 may include one or more processing cores, such as a 4-core processor, an 8-core processor, etc. The processor 401 can be implemented in the form of at least one hardware selected from the group consisting of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array).

[0203] The processor 401 may also include a main processor and a coprocessor. The main processor is a processor for processing data in the awake state, also known as a CPU (Central Processing Unit); the coprocessor is a low-power processor for processing data in the standby state.

[0204] The memory 402 may include one or more computer-readable storage media, which may be non-transitory. The memory 402 may also include high-speed random access memory and non-volatile memory, such as one or more disk storage devices and flash memory storage devices. In some embodiments of the present application, the non-transitory computer-readable storage medium in the memory 402 is used to store at least one instruction, which is used to be executed by the processor 401 to implement the method in the embodiment of the present application.

[0205] In some embodiments, the electronic device 400 further includes: a peripheral device interface 403 and at least one peripheral device 404. The processor 401, the memory 402, and the peripheral device interface 403 can be connected via a bus or signal lines. Each peripheral device 404 can be connected to the peripheral device interface 403 via a bus, signal lines, or a circuit board. Specifically, the peripheral devices 404 include: a display screen, a camera, and an audio circuit. The peripheral device interface 403 can be used to connect at least one I / O (Input / Output)-related peripheral device to the processor 401 and the memory 402.

[0206] In some embodiments of the present application, the processor 401, the memory 402, and the peripheral device interface 403 are integrated on the same chip or circuit board; in some other embodiments of the present application, any one or two of the processor 401, the memory 402, and the peripheral device interface 403 may be implemented on separate chips or circuit boards. This embodiment of the present application is not specifically limited to this.

[0207] The electronic device structure block diagram shown in the embodiment of the present application does not constitute a limitation on the electronic device 400. The electronic device 400 may include more or fewer components than shown in the figure, or combine certain components, or adopt a different component arrangement.

[0208] The present application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the method of any of the aforementioned embodiments. The computer-readable storage medium may include, but is not limited to, any type of disk, including a floppy disk, an optical disk, a DVD, a CD-ROM, a microdrive, a magneto-optical disk, a ROM, a RAM, an EPROM, an EEPROM, a DRAM, a VRAM, a flash memory device, a magnetic card or an optical card, a nanosystem (including a molecular memory IC), or any type of medium or device suitable for storing instructions and / or data.

[0209] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, or of course, by hardware. Based on this understanding, the above technical solution, in essence, or the part that contributes to the relevant technology, can be embodied in the form of a software product. The computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, a magnetic disk, an optical disk, etc., and includes a number of instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or certain parts of the embodiments.

[0210] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present application, rather than to limit them. Although the present application has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for calibrating the time difference of a low-orbit navigation enhancement transmitter, characterized in that: include: Acquiring first observation data of the GNSS satellite through an onboard receiver of the low-orbit satellite, and calculating a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data; Acquire secondary observation data of the GNSS satellite through a ground station receiver, and obtain the station receiver coordinates, station receiver clock error, and tropospheric parameters through ground precise point positioning based on the secondary observation data; Obtain downlink navigation data sent by each low-orbit satellite through the ground station, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite by combining the onboard receiver clock error, the low-orbit precise orbit, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data; Selecting any low-orbit satellite as a reference satellite, obtaining the reference inter-system bias of the reference satellite, calculating the difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and outputting the difference as the navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite; Based on the navigation enhancement transmitter time offset of each low-orbit satellite, the corresponding low-orbit satellite's onboard receiver clock error is compensated, and the low-orbit satellite navigation signal transmitter clock error with a unified time reference is output.

2. The method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to claim 1, characterized in that: The satellite receiver clock error is calculated based on the following steps, including: Establishing an observation equation based on the first observation data to obtain a pseudorange observation equation and a carrier phase observation equation; Based on the pseudorange observation equation and the carrier phase observation equation, the frequency-dependent characteristics of ionospheric delay are utilized to construct a combined observation value without ionospheric delay, and the first-order term of ionospheric delay is eliminated to obtain a pseudorange observation equation without ionospheric delay and a carrier phase observation equation without ionospheric delay; Solve the pseudorange observation equation without ionospheric delay and the carrier phase observation equation without ionospheric delay, and output the satellite receiver clock error.

3. The method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to claim 1, characterized in that: The LEO precise orbit is calculated based on the following steps, including: Establishing a motion equation of the low-orbit satellite based on the first observation data; The state transfer matrix and the orbital state information of the low-orbit satellite at each moment are calculated based on the numerical integration method; The orbital parameters of the low-orbit satellite are calculated based on a parameter estimation method in combination with the state transfer matrix and the orbital state information of the low-orbit satellite at each moment, and the low-orbit precise orbit is output.

4. The method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to claim 1, wherein: The obtaining of the station receiver coordinates, the station receiver clock error, and the tropospheric parameters by ground precise point positioning based on the second observation data includes: Establishing an observation equation based on the second observation data to obtain a pseudorange observation equation and a carrier phase observation equation of the ground station; Based on the pseudorange observation equation of the ground station and the carrier phase observation equation of the ground station, utilizing the frequency-dependent characteristics of ionospheric delay, constructing a combined observation value without ionospheric delay, eliminating the first-order term of ionospheric delay, and constructing a pseudorange observation equation without ionospheric delay and a carrier phase observation equation without ionospheric delay; The pseudorange observation equation without ionospheric delay and the carrier phase observation equation without ionospheric delay corresponding to the ground station are solved by ground precise single point positioning, and the station receiver coordinates, the station receiver clock error and the tropospheric parameters are output.

5. The method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to claim 1, characterized in that: The calculating, by combining the onboard receiver clock error, the low-orbit precise orbit of the low-orbit satellite, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data, an inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite includes: Constructing a ground station low-orbit observation equation that takes into account the hardware delay of the low-orbit satellite, including: The low-orbit satellite includes an onboard receiver for receiving GNSS signals and a navigation signal transmitter for sending downlink navigation data. The hardware delay is the clock difference between the onboard receiver clock error and the navigation signal transmitter clock error: Each time, the downlink navigation data of a low-orbit satellite is selected and substituted into the low-orbit observation equation of the ground station. Other bias terms and hardware delay are estimated as a whole to obtain the inter-system bias. The formula is applied: Among them, i is the low-orbit satellite number, G is the GNSS system, corresponds to the low-orbit satellite L i The satellite receiver clock error is corresponds to the low-orbit satellite L i The navigation signal transmitter clock error, ISB LEO,r is the other deviation term, corresponds to the low-orbit satellite L i Hardware latency, represents the geometric distance from the ground station r to the low-orbit satellite i, is the station receiver clock error, T r is the tropospheric parameter, is the wavelength corresponding to the ionosphere-free combination, is the fuzziness parameter, corresponds to the low-orbit satellite L i The inter-system deviation, ε P,IF is the measurement noise term of the pseudorange observation, ε l,IF is the measurement noise term of the carrier phase observation.

6. The method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to claim 5, characterized in that: The satellite receiver clock error is the error term due to relativistic effects The low-orbit satellite includes an onboard receiver for receiving GNSS signals and a navigation signal transmitter for sending downlink navigation data. The hardware delay is the clock error of the onboard receiver. Clock difference with the navigation signal transmitter Clock difference Application formula: Eliminate the error term of low-orbit satellites affected by relativistic effects Application formula: Δt L,ST =δt L -δt L 。 7. A method for calibrating the time difference of a low-orbit navigation enhancement transmitter according to any one of claims 1 to 6, characterized in that: The compensating the clock error of the corresponding low-orbit satellite based on the navigation enhancement transmitting end time offset of each low-orbit satellite includes: Obtain the navigation augmentation transmitter timing offset of each low-orbit satellite respectively, and compensate the navigation augmentation transmitter timing offset to the onboard receiver clock error of the corresponding low-orbit satellite using the formula: Output the compensated low-orbit navigation signal clock error sequence of each low-orbit satellite Among them, i is the low-orbit satellite number, corresponds to the low-orbit satellite L i The clock error of the navigation signal transmitter is corresponds to the low-orbit satellite L i The onboard receiver clock error before compensation.

8. A low-orbit navigation enhanced transmitter time difference calibration device, characterized in that: include: a low-orbit satellite-borne data acquisition unit, configured to acquire first observation data of a GNSS satellite through an onboard receiver of the low-orbit satellite, and calculate a clock error of the onboard receiver and a low-orbit precise orbit of the low-orbit satellite based on the first observation data; a ground station data acquisition unit, configured to acquire second observation data of the GNSS satellite through a station receiver at the ground station, and obtain station receiver coordinates, station receiver clock errors, and tropospheric parameters through ground precise point positioning based on the second observation data; a data processing unit, configured to obtain the downlink navigation data sent by each low-orbit satellite through the ground station data acquisition unit, and calculate the inter-system deviation between the low-orbit satellite system and the GNSS system corresponding to each low-orbit satellite by combining the onboard receiver clock error, the low-orbit precise orbit of the low-orbit satellite, the station receiver coordinates, the station receiver clock error, the tropospheric parameters, and the downlink navigation data; The data processing unit is further configured to select any low-orbit satellite as a reference satellite, obtain a reference inter-system bias of the reference satellite, calculate a difference between the inter-system bias corresponding to each low-orbit satellite other than the reference satellite and the reference inter-system bias, and output the difference as a navigation augmentation transmitting end timing bias of the corresponding low-orbit satellite; The data processing unit is further configured to compensate for the clock error of the onboard receiver of the corresponding low-orbit satellite based on the navigation enhancement transmitter time offset of each low-orbit satellite, and output the low-orbit satellite navigation signal transmitter clock error with a unified time reference.

9. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images