Doppler differential positioning method, device and equipment of low-orbit satellite and medium

Abstract

The invention relates to a Doppler differential positioning method, device and equipment of a low earth orbit satellite and a medium. The method comprises the following steps: calculating an ionosphere delay rate correction value and a troposphere delay rate correction value; establishing an observation equation containing ephemeris error terms of the base station and the moving station; according to the difference between the observation value of the base station and the theoretical observation value, constructing a difference correction value; constructing the differential correction of the residual atmospheric error; calculating to obtain a corrected observation value of the moving station by using the observation value of the moving station, the differential correction value and the differential correction value of the residual atmospheric error; constructing an observation equation according to the corrected moving station observation values of the plurality of satellites, and constructing a state equation by adopting a constant acceleration model; and defining the state vector as a joint estimator of the three-dimensional motion state of the moving station and the clock drift of the receiver, and resolving the accurate position of the moving station in real time by using an unscented Kalman filtering algorithm on the basis of the observation equation and the state equation. By adopting the method, the differential positioning performance can be improved.

Description

Technical Field

[0001] This application relates to the field of satellite navigation technology, and in particular to a Doppler differential positioning method, device, equipment and medium for low-Earth orbit satellites. Background Art

[0002] In recent years, the construction of low-Earth orbit satellites is in the stage of rapid deployment. Since the Doppler frequency shift of low-satellite signals is usually used as a speed measurement observable. And due to the high-speed motion characteristics of low-Earth orbit satellites, the Doppler frequency shift of their ground signals can also be used as an effective positioning observable. Therefore, the positioning technology based on the Doppler frequency shift information of low-Earth orbit satellites has become one of the research hotspots of backup navigation technology. In 2022, Neinavaie et al. used the observations of 6 Starlink satellites' complete arcs to obtain the Doppler differential positioning accuracy of 10m horizontally for the first time. Scholars such as Wang and Schüler respectively modeled the tropospheric and ionospheric errors, improving the simulation accuracy.

[0003] Differential Doppler positioning refers to the synchronous observation of a reference station (with known precise coordinates) and a rover station (user terminal), and uses the Doppler frequency shift difference generated when both receive the same satellite signal to correct the positioning error of the rover station in real time. However, different from the differential correction of GNSS, due to the low-orbit characteristics of satellites, the differential correction amount cannot completely cancel the atmospheric Doppler error of the rover station, and there is currently a lack of research on the baseline length and the degree of atmospheric error correction. In the real-time Doppler positioning experiment of the downlink signal of low-Earth orbit satellites in the actual measurement environment, the single-point positioning accuracy is still low, and there are limitations in dealing with random errors such as the ionosphere and troposphere in the existing methods. Due to the low-orbit characteristics, the residual atmospheric error will significantly affect the differential positioning performance. Summary of the Invention

[0004] Based on this, in view of the above technical problems, it is necessary to provide a Doppler differential positioning method, device, equipment and medium for low-Earth orbit satellites that can improve the differential positioning performance.

[0005] A Doppler differential positioning method for low-Earth orbit satellites, the method includes: Construct an ionospheric error model considering the Doppler characteristics of low-Earth orbits according to the Nequick-G model, and calculate the STEC change rate at the reference epoch according to the error model; use the STEC change rate to calculate the ionospheric delay rate correction value in real time; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the correction values of the atmospheric error model; Establish the observation equations of the reference station and the rover station that include the ephemeris error term; construct the differential correction by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction of the residual atmospheric error by taking the difference between the atmospheric error model correction values of the reference station and the rover station; calculate the corrected observed value of the rover station by using the observed value of the rover station, the differential correction, and the differential correction of the residual atmospheric error; Construct the observation equation based on the corrected observed values of the rover station for multiple satellites and construct the state equation using the constant acceleration model. Define the state vector as the joint estimator of the three-dimensional motion state of the rover station and the receiver clock drift, and use the unscented Kalman filter algorithm to real-time solve the precise position of the rover station based on the observation equation and the state equation.

[0006] A Doppler differential positioning device for low Earth orbit satellites, the device includes: An atmospheric error model correction value calculation module, which is used to construct an ionospheric error model considering the low Earth orbit Doppler characteristics according to the Nequick-G model, calculate the STEC change rate at the reference epoch according to the error model; use the STEC change rate to calculate the ionospheric delay rate correction value in real time; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the atmospheric error model correction values; A differential correction module, which is used to establish the observation equations of the reference station and the rover station that include the ephemeris error term; construct the differential correction by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction of the residual atmospheric error by taking the difference between the atmospheric error model correction values of the reference station and the rover station; calculate the corrected observed value of the rover station by using the observed value of the rover station, the differential correction, and the differential correction of the residual atmospheric error; A positioning module, which is used to construct the observation equation based on the corrected observed values of the rover station for multiple satellites and construct the state equation using the constant acceleration model. Define the state vector as the joint estimator of the three-dimensional motion state of the rover station and the receiver clock drift, and use the unscented Kalman filter algorithm to real-time solve the precise position of the rover station based on the observation equation and the state equation.

[0007] A computer device, including a memory and a processor, the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented: Construct an ionospheric error model considering the low Earth orbit Doppler characteristics according to the Nequick-G model, calculate the STEC change rate at the reference epoch according to the error model; use the STEC change rate to calculate the ionospheric delay rate correction value in real time; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the atmospheric error model correction values; Establish the observation equations of the reference station and the rover station that include the ephemeris error term; construct the differential correction by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction of the residual atmospheric error by taking the difference between the corrected values of the atmospheric error models of the reference station and the rover station; calculate the corrected observed value of the rover station by using the observed value of the rover station, the differential correction, and the differential correction of the residual atmospheric error; Construct the observation equation based on the corrected observed values of multiple satellites for the rover station and construct the state equation by using the constant acceleration model. Define the state vector as the combined estimator of the three-dimensional motion state of the rover station and the receiver clock drift, and use the unscented Kalman filter algorithm to real-time solve the precise position of the rover station based on the observation equation and the state equation.

[0008] A computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, the following steps are implemented: Construct the ionospheric error model considering the low-orbit Doppler characteristics according to the Nequick-G model, and calculate the STEC change rate at the reference epoch according to the error model; calculate the ionospheric delay rate correction value in real time by using the STEC change rate; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the corrected values of the atmospheric error model; Establish the observation equations of the reference station and the rover station that include the ephemeris error term; construct the differential correction by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction of the residual atmospheric error by taking the difference between the corrected values of the atmospheric error models of the reference station and the rover station; calculate the corrected observed value of the rover station by using the observed value of the rover station, the differential correction, and the differential correction of the residual atmospheric error; Construct the observation equation based on the corrected observed values of multiple satellites for the rover station and construct the state equation by using the constant acceleration model. Define the state vector as the combined estimator of the three-dimensional motion state of the rover station and the receiver clock drift, and use the unscented Kalman filter algorithm to real-time solve the precise position of the rover station based on the observation equation and the state equation.

[0009] The above-mentioned Doppler differential positioning method, device, equipment and medium for low-orbit satellites. In this application, an ionospheric error model is constructed based on the Nequick-G model, and the low-orbit Doppler characteristics are fully considered to accurately calculate the STEC change rate at the reference epoch moment, and then the ionospheric delay rate correction value is obtained in real time to accurately compensate for the ionospheric delay error. At the same time, combining the zenith tropospheric delay and the satellite elevation angle, the tropospheric delay rate correction value is calculated to effectively reduce the tropospheric delay error. Secondly, the observation equations including the ephemeris error terms of the reference station and the roving station are constructed to comprehensively cover the error factors. By taking the difference between the reference station observation value and the theoretical observation value, the differential correction amount is constructed to reduce the common errors between the reference station and the roving station. Then, by taking the difference between the corrected values of the atmospheric error models of the reference station and the roving station, the differential correction amount of the residual atmospheric error is obtained to further correct the deviation of the roving station caused by the residual atmospheric error. Finally, according to the corrected roving station observation values of multiple satellites, the observation equation is constructed, and the state equation is constructed using the constant acceleration model. The state vector is set as the joint estimator of the three-dimensional motion state of the roving station and the receiver clock drift. According to the unscented Kalman filter algorithm, integrating the observation information of multiple satellites and considering the motion and clock drift of the roving station, high-precision real-time positioning is achieved, effectively overcoming the adverse effects of the residual atmospheric error under the low-orbit characteristics on the differential positioning performance, and significantly improving the accuracy and reliability of positioning. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] Figure 1 It is a schematic flowchart of a Doppler differential positioning method for a low-orbit satellite in an embodiment; Figure 2 It is a schematic diagram of an inter-station differential model in an embodiment; Figure 3 It is a structural block diagram of a Doppler differential positioning device for a low-orbit satellite in an embodiment; Figure 4 It is an internal structure diagram of a computer device in an embodiment. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0011] In order to make the objectives, technical solutions and advantages of this application clearer, the following further details this application in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not used to limit this application.

[0012] In one embodiment, as Figure 1 shown, a Doppler differential positioning method for a low-orbit satellite is provided. Under the inter-station differential model shown in Figure 2 , the method includes the following steps: Step 102: Construct an ionospheric error model considering the low-orbit Doppler characteristics according to the Nequick-G model, and calculate the STEC change rate at the reference epoch according to the error model; use the STEC change rate to calculate the ionospheric delay rate correction value in real time; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the atmospheric error model correction values.

[0013] First, a real-time correction method based on the Nequick-G (Galileo single-frequency ionospheric correction algorithm) model was designed. The Nequick-G model calculates the effective ionization factor through the ionospheric coefficients in the corresponding broadcast ephemeris (solar activity parameters ), and the calculation formula is as follows: ; ; where is the transmitted ionospheric coefficient, is the modified magnetic declination (MODIP) at the receiver position.

[0014] The conversion formula between F10.7 (the average value of the solar radiation flux at a wavelength of 10.7 cm) and R12 (the 12-month moving average sunspot number) is as follows: ; ; The calculation process of ; ; ; ; ; ; ; In the formula, the [] symbol represents rounding; is the geographical latitude, is the longitude; is the MODIP grid attribute value ; is the third-order interpolation function; is the interpolation function at x.

[0015] Specifically, using the known satellite velocity at the reference epoch, recursively calculate The satellite position after that is regarded as a virtual epoch. Then, the satellite position at this virtual epoch is substituted into the Nequick-G model to obtain the correction value at the virtual epoch moment. Next, the difference is taken with the correction value at the reference epoch moment, so as to obtain the STEC change rate at the reference epoch moment. The specific mathematical expression is as follows: ; where is the satellite coordinate, is the Nequick-G model, and the expression adds and two adaptation parameters on the basis of the Nequick-G model, so it is called the Adaptive Nequick-G (A-Nequick) model. can effectively reduce the jump of the correction value caused by the difference between epochs; while can adjust the systematic deviation caused by the movement of low-orbit satellites in the ionosphere.

[0016] By calculating the STEC (slant total electron content) change rate at the reference epoch moment, the dynamic change of the ionospheric electron content can be grasped in real time, and then the ionospheric delay rate correction value can be calculated by using it to effectively compensate for the ionospheric delay error. The process of obtaining the ionospheric delay rate correction value in real time according to the STEC change rate at the reference epoch moment is the prior art and will not be elaborated too much in this application.

[0017] At the same time, a simplified tropospheric Doppler model considering the wet component, receiver height, and seasonal temperature is proposed: ; ; where is the tropospheric delay rate correction value in m / s. ZTD(m) is the zenith tropospheric delay empirical value, which is a function of the month value and is obtained from Table 2. and are the random angular quantities of the dry and wet components. is the zenith tropospheric delay at altitude 0; is the altitude in km. is the satellite elevation angle in radians. E is the base of the natural logarithm (approximately equal to 2.71828).

[0018] Table 1

[0019] As can be seen from Table 1, the wide range from high latitudes (Ulan Bator, 47.865°) to low latitudes (Bangkok, 13.736°) can reflect the tropospheric delay characteristics under different climate conditions. In high-latitude regions, the temperature is usually low and the water vapor is less, resulting in smaller tropospheric delays; in mid-latitude regions (such as Beijing and Wuhan), the four seasons change, and the water vapor and temperature fluctuate greatly; while in low-latitude regions (such as Hong Kong and Bangkok), it is high temperature and high humidity, which may lead to larger tropospheric delays.

[0020] Table 2 shows the mean values of ZTD measurement data at five stations in different months (January - February, March - April, May - June, July - August, September - October, and November - December), reflecting the climate characteristics of each region in different seasons. The measured values in Ulan Bator are relatively stable, while Beijing reaches the highest value in spring and summer, showing obvious seasonal variations. The values in Wuhan are relatively balanced, and Hong Kong and Bangkok show higher measured values in summer, reflecting the characteristics of tropical climate. Considering the influence of seasonal temperature, the model takes the mean values in Table 2 as the values of ZTD(m).

[0021] Table 2

[0022] Step 104: Establish the observation equations including ephemeris error terms for the reference station and the rover station; construct the differential correction amount by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction amount of the residual atmospheric error by taking the difference between the corrected values of the atmospheric error models of the reference station and the rover station; calculate the corrected observed value of the rover station using the observed value of the rover station, the differential correction amount, and the differential correction amount of the residual atmospheric error.

[0023] In view of the influence of the ephemeris error of low-earth orbit communication satellites, establish the observation equations including ephemeris error terms. The complete observation equation can be expressed as: ; In the formula, the subscript is a wildcard (replaced by when it is the reference station and replaced by when it is the rover station), represents the geometric distance from the satellite to the receiver, is the observation noise. The formula contains four physical terms: the theoretical Doppler observation value term characterizes the Doppler observation of the receiver and the true satellite state, and for the rover station, it contains the position and velocity parameters to be estimated; the ephemeris error term reflects the influence of satellite position and velocity errors; the clock drift term contains the satellite clock drift and the receiver clock drift ; the propagation delay term mainly considers the time-varying effects of the ionosphere and the troposphere . The relativistic effect is not listed separately because of its small magnitude and easy coupling to the clock error term.

[0024] Since the coordinates of the reference station are known, the theoretical observation values of the reference station can be defined as: ; Thus, the differential correction amount is constructed . In addition, according to the error model, the differential correction amount of the residual atmospheric error is constructed: ; wherein is the correction amount of the residual atmospheric error. and are the correction values of the atmospheric error models of the reference station and the roving station respectively. The receiver coordinates required to calculate the roving station correction value can be approximately replaced by the initial positioning result. Compensating the correction amount to the roving station observation value, we can get: ; is the roving station observation value corrected by the differential correction information. Based on the "parallel" assumption in the previous text, includes errors such as ephemeris, atmosphere, and satellite clock drift that are approximately equal to the various errors of the roving station. And includes a part of the residual atmospheric error caused by spatio-temporal correlation. Therefore, this inter-station differential operation can effectively suppress the common error terms including ephemeris error of the roving station and significantly improve the positioning accuracy. The corrected roving station observation value only contains the difference between the receiver clock drift and the uncorrected observation noise, and the specific expression is ; wherein are the observation noise values of the roving station and the reference station respectively. While eliminating the systematic error, this model retains the Doppler observables directly related to the motion state of the roving station, providing an observation basis for subsequent UKF filtering.

[0025] Step 106, construct the observation equation based on the corrected roving station observation values of multiple satellites and construct the state equation using the constant acceleration model. Define the state vector as the joint estimator of the three-dimensional motion state of the roving station and the receiver clock drift, and use the unscented Kalman filter algorithm to solve the accurate position of the roving station in real time based on the observation equation and the state equation.

[0026] Define the state vector at the kth epoch as the joint estimator of the three-dimensional motion state of the roving station and the receiver clock drift : ; The observation equation is composed of the differential correction Doppler frequency shifts of N satellites output by the satellite selection algorithm: ; The state transition model adopts the constant acceleration model, and its discrete state equation is as follows: ; where is the epoch interval, and the process noise covariance matrix needs to be adjusted according to the dynamic characteristics of the carrier.

[0027] The UKF accurately propagates the state distribution through the unscented transform. Its core steps are as follows: First, generate the Sigma point set to capture the statistical characteristics of the state mean and covariance; propagate the Sigma points through the nonlinear observation model, and calculate the predicted observation and its covariance; update the state estimate by combining the actual observation value. To improve the calculation efficiency, a symmetric sampling strategy is adopted to reduce the number of Sigma points, and a dynamic adjustment mechanism is introduced to optimize the scaling parameter , to balance the high-order moment approximation accuracy and numerical stability. Regarding the design problem of the observation noise covariance matrix, since the atmospheric error is strongly correlated with the elevation angle, the elevation angle weighted model commonly used in GNSS is adopted. The specific formula for the weight is as follows: ; where is the elevation angle of the th satellite.

[0028] The process of using the unscented Kalman filter algorithm to real-time calculate the accurate position of the rover is shown in Table 3.

[0029] Table 3

[0030] For the above Doppler differential positioning method of low-earth orbit satellites, the present application constructs an ionospheric error model based on the Nequick-G model, fully considers the Doppler characteristics of low-earth orbits, accurately calculates the STEC change rate at the reference epoch, and then obtains the ionospheric delay rate correction value in real time to accurately compensate for the ionospheric delay error. At the same time, in combination with the zenith tropospheric delay and the satellite elevation angle, the tropospheric delay rate correction value is calculated to effectively reduce the tropospheric delay error. Secondly, observation equations including ephemeris error terms for the reference station and the rover are constructed to comprehensively cover error factors. By taking the difference between the observed values at the reference station and the theoretical observed values, a differential correction quantity is constructed to reduce the common errors between the reference station and the rover. Then, by taking the difference between the correction values of the atmospheric error models at the reference station and the rover, the differential correction quantity of the residual atmospheric error is obtained to further correct the deviation of the rover caused by the residual atmospheric error. Finally, based on the observed values of the rover corrected by multiple satellites, an observation equation is constructed, and a state equation is constructed using the constant acceleration model. The state vector is set as the joint estimator of the three-dimensional motion state of the rover and the receiver clock drift. According to the unscented Kalman filter algorithm, integrating the observed information of multiple satellites and considering the motion and clock drift of the rover, high-precision real-time positioning is achieved, effectively overcoming the adverse effects of the residual atmospheric error under low-orbit characteristics on the differential positioning performance, and significantly improving the accuracy and reliability of positioning.

[0031] In one embodiment, calculating the STEC change rate at the reference epoch according to the ionospheric error model includes: Using the known satellite velocity at the reference epoch, recursively calculate the satellite position after ; where is the satellite coordinate, is the Nequick-G model, represents a value that can reduce the jump of the correction value caused by the difference between epochs; represents the time interval.

[0032] In one embodiment, calculating the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle includes: Calculating the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle is ; ; where is the tropospheric delay rate correction value in m / s, and ZTD(m) is the zenith tropospheric delay empirical value. and are the random angular quantities of the dry and wet components. is the zenith tropospheric delay at altitude 0. is the altitude in km. is the satellite elevation angle in radians.

[0033] In one embodiment, the observation equations including the ephemeris error terms for the reference station and the rover are established, including: ; where the subscript is a wildcard, which is replaced by for the reference station and for the rover. represents the geometric distance from the satellite to the receiver. is the observation noise. The theoretical Doppler observation value term characterizes the Doppler observation of the receiver and the true satellite state. represents the velocity vector of the satellite. represents the velocity vector of the reference station or the rover. represents the unit direction vector from the satellite to the reference station or the rover. represents the satellite position. represents the speed of light, and for the rover, it includes the position and velocity parameters to be estimated. The ephemeris error term reflects the influence of the satellite position and velocity errors. The clock drift term includes the satellite clock drift and the receiver clock drift ; represents the propagation delay term of the ionosphere. represents the propagation delay term of the troposphere.

[0034] In one embodiment, the corrected rover observation values are calculated using the rover's observation values, differential correction amounts, and differential correction amounts of the residual atmospheric errors, including: The corrected rover observation values calculated using the rover's observation values, differential correction amounts, and differential correction amounts of the residual atmospheric errors are ; where represents the rover's observation value. represents the differential correction amount. represents the theoretical observation value. represents the velocity vector of the reference station. represents the unit direction vector of the reference station. represents the carrier frequency. represents the velocity vector of the rover. represents the unit direction vector of the rover station, represents the receiver clock drift term, represents the observation value of the reference station, represents the differential correction of the residual atmospheric error, and are the correction values of the atmospheric error models of the reference station and the rover station respectively, are the observation noise values of the rover station and the reference station respectively.

[0035] In one embodiment, the state vector is defined as the joint estimator of the three-dimensional motion state of the rover station and the receiver clock drift, including: Define the k epoch state vector as the joint estimator of the three-dimensional motion state of the rover station and the receiver clock drift is ; where, T represents the transpose operation.

[0036] In one embodiment, an observation equation is constructed based on the observed values of the rover station corrected by multiple satellites, and a state equation is constructed using a constant acceleration model, including: The observation equation constructed based on the observed values of the rover station corrected by multiple satellites is: ; The state equation constructed using a constant acceleration model is: ; where, N represents the satellite number, represents the observed value of the rover station after correction, is the epoch interval, represents the process noise covariance matrix.

[0037] It should be understood that although Figure 1 the steps in the flowchart of Figure 1 are shown in sequence according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless there is a clear indication in this article, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover,

[0038] In one embodiment, as Figure 3 shown, a Doppler differential positioning device for a low-earth orbit satellite is provided, including: an atmospheric error model correction value calculation module 302, a differential correction amount module 304, and a positioning module 306, where: The atmospheric error model correction value calculation module 302 is configured to construct an ionospheric error model considering the Doppler characteristics of the low-earth orbit according to the Nequick-G model, calculate the STEC change rate at the reference epoch according to the error model; calculate the ionospheric delay rate correction value in real time by using the STEC change rate; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the atmospheric error model correction values; The differential correction amount module 304 is configured to establish an observation equation including the ephemeris error term for the reference station and the rover; construct a differential correction amount by taking the difference between the observed value of the reference station and the theoretical observed value; construct a differential correction amount for the residual atmospheric error by taking the difference between the atmospheric error model correction values of the reference station and the rover; calculate the corrected observed value of the rover by using the observed value of the rover, the differential correction amount, and the differential correction amount of the residual atmospheric error; The positioning module 306 is configured to construct an observation equation based on the corrected observed values of multiple satellites and construct a state equation by using a constant acceleration model, define the state vector as a joint estimator of the three-dimensional motion state of the rover and the receiver clock drift, and use the unscented Kalman filter algorithm to solve the exact position of the rover in real time based on the observation equation and the state equation.

[0039] For the specific limitations of the Doppler differential positioning device for the low-earth orbit satellite, reference can be made to the limitations of the Doppler differential positioning method for the low-earth orbit satellite in the above text, which will not be elaborated here. Each module in the above Doppler differential positioning device for the low-earth orbit satellite can be implemented in whole or in part by software, hardware, and their combination. The above modules can be embedded in the processor in the computer device in the form of hardware or be independent of it, or can be stored in the memory in the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0040] In one embodiment, a computer device is provided. The computer device can be a terminal, and its internal structure diagram can be as Figure 4As shown in the figure. The computer device includes a processor, a memory, a network interface, a display screen, and an input device connected via a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external terminal via a network connection. When the computer program is executed by the processor, it implements a Doppler differential positioning method for low-earth orbit satellites. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen. The input device of the computer device can be a touch layer covering the display screen, or a button, a trackball, or a touchpad provided on the housing of the computer device, or an external keyboard, touchpad, or mouse, etc.

[0041] Those skilled in the art can understand that Figure 4 the structure shown in the figure is only a block diagram of some structures related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0042] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, storage, database, or other medium used in the various embodiments provided in the present application can include non-volatile and / or volatile memories. Non-volatile memories can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memories can include random access memory (RAM) or external cache memory. By way of illustration and not limitation, RAM is available in various forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM), etc.

[0043] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope recorded in this specification.

[0044] The above-described embodiments merely represent several implementation manners of the present application. The description is relatively specific and detailed, but it should not be construed as a limitation on the scope of the present application. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the appended claims.

Claims

1. A Doppler differential positioning method for low-earth orbit satellites, characterized in that, The method includes: Constructing an ionospheric error model considering the characteristics of low-orbit Doppler according to the Nequick-G model, and calculating the STEC change rate at the reference epoch according to the error model; calculating the ionospheric delay rate correction value in real time by using the STEC change rate; calculating the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the correction values of the atmospheric error model; Establishing the observation equations of the reference station and the roving station including the ephemeris error term; constructing the differential correction amount by taking the difference between the observed value of the reference station and the theoretical observed value; constructing the differential correction amount of the residual atmospheric error by taking the difference between the correction values of the atmospheric error model of the reference station and the roving station; calculating the corrected observed value of the roving station by using the observed value of the roving station, the differential correction amount and the differential correction amount of the residual atmospheric error; Constructing the observation equations according to the corrected observed values of the roving station of multiple satellites and constructing the state equations by using the constant acceleration model, defining the state vector as the joint estimator of the three-dimensional motion state of the roving station and the receiver clock drift, and solving the precise position of the roving station in real time by using the unscented Kalman filter algorithm based on the observation equations and the state equations.

2. The method according to claim 1, characterized in that, Calculating the STEC change rate at the reference epoch according to the ionospheric error model, including: Using the known satellite velocity at the reference epoch, recursively calculate the satellite position after and regard it as a virtual epoch. Substitute the satellite position of this virtual epoch into the Nequick-G model. After obtaining the correction value at the virtual epoch, subtract it from the correction value at the reference epoch to obtain the STEC change rate at the reference epoch. The specific formula is as follows: ； Among them, is the satellite coordinate, is the Nequick-G model, represents a value that can reduce the jump of the correction value caused by the epoch difference; represents the time interval.

3. The method according to claim 1, characterized in that, Calculating the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle, including: Calculating the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle as ； ； Among them, is the correction value of the tropospheric delay rate in m / s, and ZTD(m) is the empirical value of the zenith tropospheric delay, and are the random angular quantities of the dry and wet components, is the zenith tropospheric delay at an altitude of 0; is the altitude in km, is the satellite elevation angle in radians.

4. The method according to claim 1, wherein Establishing the observation equations of the reference station and the roving station including the ephemeris error term, including: ； Among them, the subscript is a wildcard, which is replaced by at the reference station and by at the rover station. represents the geometric distance from the satellite to the receiver. is the observation noise. The theoretical Doppler observation term characterizes the Doppler observation of the receiver and the true satellite state. represents the velocity vector of the satellite. represents the velocity vector of the reference station or the rover station. represents the unit direction vector from the satellite to the reference station or the rover station. represents the satellite position. represents the speed of light, which includes the position and velocity parameters to be estimated for the rover station; the ephemeris error term reflects the influence of satellite position and velocity errors; the clock drift term includes the satellite clock drift and the receiver clock drift . represents the propagation delay term of the ionosphere. represents the propagation delay term of the troposphere.

5. The method according to claim 1, wherein Calculating the corrected observed value of the roving station by using the observed value of the roving station, the differential correction amount and the differential correction amount of the residual atmospheric error, including: Calculating the corrected observed value of the roving station by using the observed value of the roving station, the differential correction amount and the differential correction amount of the residual atmospheric error as ； Among them, represents the observation value of the rover station, represents the differential correction amount, represents the theoretical observation value, represents the velocity vector of the reference station, represents the unit direction vector of the reference station, represents the carrier frequency, represents the velocity vector of the rover station, represents the unit direction vector of the rover station, represents the receiver clock drift term, represents the observation value of the reference station, represents the differential correction amount of the residual atmospheric error, and are the correction values of the atmospheric error models of the reference station and the rover station respectively, are the observation noise values of the rover station and the reference station respectively.

6. The method according to claim 1, wherein Defining the state vector as the joint estimator of the three-dimensional motion state of the roving station and the receiver clock drift, including: Define the k epoch state vector as the combined estimator of the three-dimensional motion state of the rover station and the receiver clock drift is ； Where, T represents the transpose operation.

7. The method according to claim 1, characterized in that, Constructing the observation equations according to the corrected observed values of the roving station of multiple satellites and constructing the state equations by using the constant acceleration model, including: Constructing the observation equations according to the corrected observed values of the roving station of multiple satellites as: ； Constructing the state equations by using the constant acceleration model as: ； Among them, N represents the serial number of the satellite, represents the corrected rover observation value, is the epoch interval, represents the process noise covariance matrix.

8. A Doppler differential positioning device for low-earth orbit satellites, characterized in that, The device includes: An atmospheric error model correction value calculation module, which is used to construct an ionospheric error model considering the characteristics of low-orbit Doppler according to the Nequick-G model, calculate the STEC change rate at the reference epoch according to the error model; calculate the ionospheric delay rate correction value in real time by using the STEC change rate; calculate the tropospheric delay rate correction value according to the zenith tropospheric delay and the satellite elevation angle; the ionospheric delay rate correction value and the tropospheric delay rate correction value are the correction values of the atmospheric error model; The differential correction module is used to establish the observation equations of the reference station and the rover station including the ephemeris error terms; construct the differential correction by taking the difference between the observed value and the theoretical observed value of the reference station; construct the differential correction of the residual atmospheric error by taking the difference between the corrected values of the atmospheric error models of the reference station and the rover station; calculate the corrected observed value of the rover station by using the observed value of the rover station, the differential correction, and the differential correction of the residual atmospheric error. The positioning module is used to construct the observation equation based on the corrected observed values of multiple satellites of the rover station and construct the state equation by using the constant acceleration model, define the state vector as the joint estimator of the three-dimensional motion state of the rover station and the receiver clock drift, and use the unscented Kalman filter algorithm to solve the exact position of the rover station in real time based on the observation equation and the state equation.

9. A computer device, comprising a memory and a processor, the memory storing a computer program, characterized in that, When the processor executes the computer program, it realizes the steps of the method according to any one of claims 1 to 7.

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

Images