WO2024109466A1

WO2024109466A1 - Satellite positioning accuracy marking method and related apparatus

Info

Publication number: WO2024109466A1
Application number: PCT/CN2023/127852
Authority: WO
Inventors: 赖允斌; 王理砚; 司徒春辉; 陈亮; 韩雷晋; 施垒
Original assignee: 广州导远电子科技有限公司
Priority date: 2022-11-25
Filing date: 2023-10-30
Publication date: 2024-05-30
Also published as: CN115755113A

Abstract

A satellite positioning accuracy marking method and a related apparatus, relating to the field of satellite navigation. The method comprises: first receiving observation data of and differential data of a satellite (S20); according to the observation data and the differential data, performing ambiguity resolution processing to obtain a positioning result (S21); and subsequently, determining an ambiguity resolution processing state and, according to the ambiguity resolution processing state, determining the accuracy of the positioning result (S22). After performing ambiguity resolution processing according to the observation data and the differential data, the method can determine accuracy of positioning results according to ambiguity resolution states. Thus, the method both can determine the accuracy of the positioning results according to the ambiguity resolution states regardless of the obtained positioning results being fixed solutions or float solutions.

Description

Satellite positioning accuracy marking method and related device

CROSS-REFERENCE TO RELATED APPLICATIONS

This disclosure claims priority to a Chinese patent application with application number 2022114933850 filed with the Chinese Patent Office on November 25, 2022, and entitled “Precision Marking Method and Related Device for Satellite Positioning,” the entire contents of which are incorporated by reference in this disclosure.

Technical Field

The present disclosure relates to the field of satellite navigation, and in particular to a satellite positioning accuracy marking method and related devices.

Background technique

At present, when using satellite navigation positioning technology to calculate positioning results, it is often necessary to solve the unknown phase floating-point ambiguity and try to fix it to an integer to obtain a floating-point solution or a fixed solution. Among them, the accuracy of the fixed solution is at the centimeter level, while the accuracy of the floating-point solution often varies from decimeter level to meter level. The solutions in the prior art often cannot determine the accuracy of the floating-point solution when obtaining it.

Summary of the invention

In order to overcome the above-mentioned deficiencies in the prior art, the present invention aims to provide a method for marking the accuracy of satellite positioning, the method comprising:

Receive satellite observation data and differential data;

Performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result;

An ambiguity fixation processing state is determined, and an accuracy of the positioning result is determined according to the ambiguity fixation processing state.

Further, in the above method, the determining the ambiguity fixing processing state, and determining the accuracy of the positioning result according to the ambiguity fixing processing state, includes:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity is fixed successfully, determining that the ambiguity fixing processing state is the first state;

The accuracy of the positioning result is determined to be a first accuracy according to the first state.

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity fixation fails, determining that the ambiguity fixation processing state is the second state;

The accuracy of the positioning result is determined to be a second accuracy according to the second state.

Furthermore, in the above method, the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state includes:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than the preset number, and the narrow lane ambiguity is fixed successfully, determining that the ambiguity fixing processing state is the third state;

The accuracy of the positioning result is determined to be a third accuracy according to the third state.

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than the preset number, and the narrow lane ambiguity fixation fails, determining that the ambiguity fixation processing state is the fourth state;

The accuracy of the positioning result is determined to be a fourth accuracy according to the fourth state.

When the observation equation is solved successfully and the coordinate variance of the receiving device is greater than or equal to the preset variance, determining that the ambiguity fixing processing state is the fifth state;

The accuracy of the positioning result is determined to be a fifth accuracy according to the fifth state.

In case the observation equation solution fails, determining the ambiguity fixing processing state to be the sixth state;

The accuracy of the positioning result is determined according to the sixth state as a sixth accuracy.

In the case of pseudorange single point positioning failure, determining the ambiguity fix processing state to be the seventh state;

The positioning result is determined to be an invalid result according to the seventh state.

Furthermore, in the above method, performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result includes:

When the pseudo-range single point positioning is successful, constructing an observation equation according to the observation data and the differential data, and calculating the phase floating point ambiguity of each satellite according to the observation equation;

Performing wide lane ambiguity fixation on the phase floating point ambiguity of each of the satellites to obtain wide lane fixed ambiguity of each of the satellites;

The narrow lane ambiguity is fixed for the wide lane fixed ambiguity of each satellite to obtain a positioning result.

Furthermore, in the above method, when the pseudo-range single point positioning is successful, constructing the observation equation according to the observation data and the differential data includes:

When the pseudo-range single point positioning is successful, the observation equation is constructed according to the following formula:

Where P _s,j represents the pseudorange of satellite s at frequency j; L _s,j represents the phase observation value of satellite s at frequency j. represents the distance between the antenna phase center of the receiving device r and the phase center of the satellite s, C represents the speed of light in a vacuum, dt _r,j represents the clock error of the receiving device r at frequency j, dt ^s represents the clock error of satellite s, T represents the wet tropospheric delay, γ represents the ratio between the squares of multiple frequencies, represents the ionospheric delay of satellite s at frequency j, b _r,j represents the pseudorange hardware delay of receiving device r at frequency j, represents the pseudorange hardware delay of satellite s at frequency j, B _r,j represents the phase hardware delay of receiving device r at frequency j, Characterizes the phase hardware delay of satellite s at frequency j, Characterizes the wavelength of satellite s at frequency j, represents the phase floating point ambiguity of satellite s at frequency j, represents the ionospheric delay of the receiving device, ε _I,j represents the observation noise of the pseudorange, represents the tropospheric delay of the vehicle terminal, and ε _T,j represents the observation noise of the phase observation value.

Further, in the above method, performing wide lane ambiguity fixing on the phase floating point ambiguity of each satellite to obtain the wide lane fixed ambiguity of each satellite includes:

Selecting a reference satellite from all the satellites according to a preset rule, and using the other satellites except the reference satellite as mobile satellites;

Calculating the corresponding inter-satellite single-difference wide-lane floating ambiguity according to the phase floating ambiguity corresponding to each frequency of each mobile satellite;

If the inter-satellite single-difference wide-lane floating-point ambiguity of each of the mobile satellites meets a preset condition, then the integer wide-lane ambiguity corresponding to the mobile satellite is calculated according to the inter-satellite single-difference wide-lane floating-point ambiguity;

The constrained ambiguity parameters are updated as the inter-satellite single-difference wide-lane floating-point ambiguity of the mobile satellite according to the integer wide-lane ambiguity to obtain the wide-lane fixed ambiguity of the mobile satellite.

Furthermore, in the above method, the step of calculating the corresponding inter-satellite single-difference wide-lane floating point ambiguity according to the phase floating point ambiguity corresponding to each frequency of each mobile satellite comprises:

The corresponding inter-satellite single-difference wide-lane floating ambiguity is calculated according to the phase floating ambiguity corresponding to each frequency of each mobile satellite and the phase floating ambiguity corresponding to each frequency of the reference satellite.

Further, in the above method, the calculating corresponding inter-satellite single-difference wide-lane floating ambiguity according to the phase floating ambiguity corresponding to each frequency of each of the mobile satellites and the phase floating ambiguity corresponding to each frequency of the reference satellite comprises:

For each of the mobile satellites, the inter-satellite single-difference wide-lane floating point ambiguity is calculated by the following formula:
WL＝N _y,1 -N _y,2 -(N _c,1 -N _C,2 )

Among them, WL represents the inter-satellite single-difference wide-lane floating-point ambiguity, N _y,1 represents the phase floating-point ambiguity of the mobile satellite at the first frequency, N _y,2 represents the phase floating-point ambiguity of the mobile satellite at the second frequency, N _c,1 represents the phase floating-point ambiguity of the reference satellite at the first frequency, and N _C,2 represents the phase floating-point ambiguity of the reference satellite at the second frequency.

Further, in the above method, if the inter-satellite single-difference wide-lane floating-point ambiguity of each of the mobile satellites meets a preset condition, then calculating the integer wide-lane ambiguity corresponding to the mobile satellite according to the inter-satellite single-difference wide-lane floating-point ambiguity includes:

For each of the mobile satellites, if the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity is less than a preset threshold, the inter-satellite single-difference wide-lane floating-point ambiguity is rounded to an integer to obtain the integer wide-lane ambiguity.

Further, in the above method, updating the constrained ambiguity parameters for the inter-satellite single-difference wide-lane floating-point ambiguity of the mobile satellite according to the integer wide-lane ambiguity to obtain the wide-lane fixed ambiguity of the mobile satellite includes:

The constrained ambiguity parameters are updated for the inter-satellite single-difference wide-lane floating ambiguity of the mobile satellite by the following formula:
WL ⁰ = WL + ε

Among them, WL ⁰ represents the integer wide lane ambiguity, WL represents the inter-satellite single-difference wide lane floating point ambiguity, and ε represents the error value.

Further, in the above method, performing narrow lane ambiguity fixing on the wide lane fixed ambiguity of each satellite to obtain a positioning result includes:

Fixing narrow lane ambiguities of wide lane fixed ambiguities of the plurality of satellites according to a preset search algorithm to obtain inter-satellite single difference integer narrow lane ambiguities and ratio values;

When the ratio value does not reach a preset ratio, a target wide lane fixed ambiguity is deleted from the wide lane fixed ambiguities of the plurality of satellites, and narrow lane ambiguity is fixed again for the other wide lane fixed ambiguities according to a preset search algorithm; wherein the target wide lane fixed ambiguity is a wide lane fixed ambiguity corresponding to a satellite with the lowest satellite elevation angle among the plurality of satellites;

If the ratio value reaches a preset ratio, and the number of remaining wide lane fixed ambiguities is greater than or equal to a preset number, it is determined that the narrow lane ambiguity is fixed successfully, and the inter-satellite single-difference integer narrow lane ambiguity is used as an update constraint for the inter-satellite single-difference narrow lane floating point ambiguity to obtain the positioning result; wherein the inter-satellite single-difference narrow lane floating point ambiguity is obtained by solving the observation equation; if the ratio value does not reach the preset ratio, or the number of wide lane fixed ambiguities for narrow lane ambiguity fixation is less than a preset number, it is determined that the narrow lane ambiguity fixation fails, and the positioning result is obtained.

Another object of the present disclosure is to provide a satellite positioning accuracy marking device, the device comprising:

A receiving module, used for receiving satellite observation data and differential data;

A processing module, used for performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result;

The determination module is used to determine an ambiguity fixing processing state, and determine the accuracy of the positioning result according to the ambiguity fixing processing state.

Another object of the present disclosure is to provide a receiving device, comprising a processor and a memory, wherein the memory stores a computer program that can be executed by the processor, and the processor can execute the computer program to implement any of the methods described in the aforementioned embodiments.

Another object of the present disclosure is to provide a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method as described in any one of the aforementioned embodiments.

Compared with the prior art, the present invention has the following beneficial effects:

The accuracy marking method and related device of satellite positioning provided by the embodiment of the present disclosure first receive the observation data and differential data of the satellite, perform ambiguity fixing processing according to the observation data and differential data, thereby obtaining the positioning result, and then determine the ambiguity fixing state, and determine the accuracy of the positioning result according to the ambiguity fixing state. After the ambiguity fixing processing is performed according to the observation data and differential data, the method can determine the accuracy of the positioning result according to the ambiguity fixing state, so no matter whether the obtained positioning result is a fixed solution or a floating point solution, the accuracy of the positioning result can be determined according to the ambiguity fixing state.

BRIEF DESCRIPTION OF THE DRAWINGS

In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the drawings required for use in the embodiments will be briefly introduced below. It should be understood that the following drawings only show certain embodiments of the present disclosure and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

FIG1 shows a block diagram of a receiving device provided by an embodiment of the present disclosure;

FIG2 shows a schematic flow chart of a method for marking the accuracy of satellite positioning provided by an embodiment of the present disclosure;

FIG3 shows another schematic flow chart of a method for marking the accuracy of satellite positioning provided by an embodiment of the present disclosure;

FIG4 shows another schematic flow chart of a method for marking the accuracy of satellite positioning provided by an embodiment of the present disclosure;

FIG5 shows another schematic flow chart of the satellite positioning accuracy marking method provided by an embodiment of the present disclosure;

FIG6 is a schematic diagram showing the accuracy level of positioning results;

FIG. 7 shows a functional module diagram of a satellite positioning accuracy marking device provided in an embodiment of the present disclosure.

Icon: 10 - receiving device; 100 - memory; 110 - processor; 120 - communication module; 200 - receiving module; 210 - processing module; 220 - determination module.

Detailed ways

The following will be combined with the drawings in the embodiments of the present disclosure to clearly and completely describe the technical solutions in the embodiments of the present disclosure. Obviously, the described embodiments are only part of the embodiments of the present disclosure, not all of the embodiments. The components of the embodiments of the present disclosure described and shown in the drawings here can be arranged and designed in various different configurations.

Therefore, the following detailed description of the embodiments of the present disclosure provided in the accompanying drawings is not intended to limit the scope of the present disclosure claimed for protection, but merely represents selected embodiments of the present disclosure. Based on the embodiments of the present disclosure, all other embodiments obtained by those skilled in the art without making creative efforts are within the scope of protection of the present disclosure.

It should be noted that relational terms such as "first" and "second" are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any actual relationship or order between these entities or operations. Moreover, the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device that includes a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements for such process, method, article or device. Elements inherent to a product or device. In the absence of further restrictions, an element defined by the sentence "comprising a ..." does not exclude the presence of other identical elements in the process, method, article or device comprising the element. At present, satellite navigation positioning technology can be used for positioning, especially in the field of autonomous driving, which often uses the fusion positioning technology of vehicle-mounted combined navigation, which generally uses GNSS (Global Navigation Satellite System) satellite positioning technology to provide the absolute position coordinates of the vehicle.

When using satellite navigation positioning technology to calculate positioning results, it is often necessary to solve the unknown phase floating-point ambiguity and try to fix it to an integer. Since the phase floating-point ambiguity has integer characteristics in theory, when the phase floating-point ambiguity is successfully fixed to an integer, it is called a fixed solution in GNSS satellite positioning technology. When the phase floating-point ambiguity is not successfully fixed to an integer, it is called a floating-point solution in GNSS satellite positioning technology. Among them, the accuracy of the fixed solution is at the centimeter level, while the accuracy of the floating-point solution often varies from decimeters to meters.

In the prior art, since the accuracy of the floating-point solution is unreliable, the accuracy of the floating-point solution is often assumed to be at the meter level, that is, even if the floating-point solution is actually at the decimeter level, the electronic device will process it as a meter-level accuracy. At present, there is a general tendency to obtain a fixed solution and use the positioning result of the fixed solution for positioning. However, in actual application, the positioning accuracy requirement of some application scenarios is only at the decimeter level, and the solutions in the prior art cannot distinguish which floating-point solutions are at the decimeter level and which are at the meter level. In application scenarios where the positioning accuracy requirement is only at the decimeter level, a centimeter-level fixed solution must also be used for positioning. As a result, the value of the decimeter-level floating-point solution cannot be fully utilized, which greatly reduces the availability of GNSS satellite positioning results.

Therefore, the present disclosure provides a satellite positioning accuracy marking method to solve the above problems. Please refer to FIG1 , which is a block diagram of a receiving device 10 .

Optionally, the receiving device 10 can be set on a device that needs satellite positioning. The receiving device 10 can receive data sent by the satellite and calculate the coordinate information of the receiving device. It can be understood that the coordinate information of the receiving device can be used as the coordinate information of the device that needs satellite positioning.

Optionally, the device that requires satellite positioning may be a vehicle equipped with autonomous driving technology, that is, the receiving device 10 may be set on the vehicle equipped with autonomous driving technology to provide the absolute position coordinates of the vehicle through GNSS satellite positioning technology.

Optionally, the receiving device 10 includes a memory 100, a processor 110, and a communication module 120. The memory 100, the processor 110, and the communication module 120 are electrically connected to each other directly or indirectly to achieve data transmission or interaction. For example, these components can be electrically connected to each other through one or more communication buses or signal lines.

The memory 100 is used to store programs or data. The memory 100 may be, but is not limited to, a random access memory (RAM), a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable read-only memory (EEPROM), etc.

The processor 110 is used to read/write data or programs stored in the memory and execute corresponding functions.

The communication module 120 is used to establish a communication connection between the server and other communication terminals through the network, and to send and receive data through the network.

It should be understood that the structure shown in FIG1 is only a schematic diagram of the structure of the receiving device 10, and the receiving device 10 may also include more or fewer components than those shown in FIG1, or have a configuration different from that shown in FIG1. Each component shown in FIG1 may be implemented by hardware, software, or a combination thereof.

The embodiments of the present disclosure also provide a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the accuracy marking method for satellite positioning provided by the embodiments of the present disclosure can be implemented.

Next, the receiving device 10 in FIG. 1 is used as the execution subject, and the accuracy marking method of satellite positioning provided by the embodiment of the present disclosure is exemplarily described in combination with the flowchart. Specifically, FIG. 2 is a flowchart of the accuracy marking method of satellite positioning provided by the embodiment of the present disclosure. Referring to FIG. 2, the method includes:

Step S20, receiving satellite observation data and differential data;

Optionally, the observation data may be observation data of a GNSS navigation satellite, and the differential data may be differential data of a geostationary communication satellite.

Optionally, the GNSS navigation satellite may include a GPS satellite navigation system, a Galileo satellite navigation system, a BDS Beidou satellite navigation system, etc., and the receiving device may obtain observation data sent by one or more of these satellite navigation systems.

Optionally, the differential data may be provided by a satellite-based augmentation service provider.

Optionally, the observation data may include pseudorange and carrier phase observations and Doppler observations, and the differential data may include precise orbit clock error data, satellite-side pseudorange and phase hardware delay data, and ionospheric and tropospheric delay data on the satellite signal propagation path.

Step S21, performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result;

Step S22, determining the ambiguity fixing processing state, and determining the accuracy of the positioning result according to the ambiguity fixing processing state.

Optionally, the positioning result is coordinate information of the receiving device.

Optionally, the ambiguity fixation processing state may characterize the processing status of each processing process when ambiguity fixation processing is performed based on observation data and differential data.

Optionally, the accuracy of the positioning result may include the accuracy of a fixed solution and the accuracy of a floating-point solution. It can be understood that if the positioning result obtained is a fixed solution, the accuracy of the positioning result is at the centimeter level; if the positioning result obtained is a floating-point solution, the method can be used to determine whether the accuracy of the positioning result is specifically at the decimeter level or the meter level.

Obviously, this method can determine whether the accuracy of the floating-point solution is at the decimeter level or the meter level. Therefore, in actual application, for some application scenarios where the positioning accuracy requirement is only at the decimeter level, the floating-point solution with decimeter accuracy can be directly used for positioning without having to obtain a fixed solution. This method can make full use of the value of the decimeter-level floating-point solution and improve the availability of GNSS satellite positioning results.

The accuracy marking method of satellite positioning provided by the embodiment of the present disclosure first receives the observation data and differential data of the satellite, performs ambiguity fixing processing according to the observation data and differential data, thereby obtaining a positioning result, and then determines the ambiguity fixing state, and determines the accuracy of the positioning result according to the ambiguity fixing state. After performing ambiguity fixing processing according to the observation data and differential data, the method can determine the accuracy of the positioning result according to the ambiguity fixing state, so no matter whether the obtained positioning result is a fixed solution or a floating point solution, the accuracy of the positioning result can be determined according to the ambiguity fixing state.

Optionally, the ambiguity fixation processing process may include an observation equation solving process, a wide lane ambiguity fixation process, and a narrow lane ambiguity fixation process, wherein the phase floating-point ambiguity of each satellite and the initial coordinates of the receiving device can be obtained by solving the observation equation, and the wide lane ambiguity fixation process and the narrow lane ambiguity fixation process can improve the accuracy of the initial coordinates of the receiving device, thereby obtaining centimeter-level, decimeter-level, or meter-level coordinates of the receiving device.

Specifically, based on FIG. 2 , FIG. 3 is another flow chart of the satellite positioning accuracy marking method provided by the embodiment of the present disclosure. Referring to FIG. 3 , the above step S21 can be implemented by the following steps:

Step S21-1, when the pseudo-range single point positioning is successful, construct an observation equation based on the observation data and the differential data, and calculate the phase floating point ambiguity of each satellite according to the observation equation;

Optionally, before constructing the observation equation, pseudo-range single-point positioning needs to be performed in advance. If the pseudo-range single-point positioning is successful, the observation equation can be constructed and solved. If the pseudo-range single-point positioning fails, the positioning result is an invalid solution and there is no need to perform subsequent steps.

Alternatively, the observation equation can be constructed by:

in, (x ^s ,y ^s ,z ^s ) are the coordinates of satellite s, and (x _r ,y _r ,z _r ) are the coordinates of receiving device r.

It can be understood that the coordinates of the satellite s are known parameters, and the coordinates of the receiving device r are unknown parameters. By solving the observation equation, the unknown coordinates of the receiving device r can be obtained.

Optionally, the observation equation also includes other unknown parameters, and the receiving device can obtain the unknown data therein by solving the equation, including the coordinates of the receiving device r, the clock error of the receiving device r at the j frequency, the wet tropospheric delay, the ionospheric delay of the receiving device, and the phase floating point ambiguity of the satellite at each frequency.

It can be understood that the receiving device can construct the observation equation for each satellite, thereby obtaining the phase floating point ambiguity of each satellite at each frequency.

Optionally, in this equation, the hardware delay of the receiving device is absorbed by the clock error of the receiving device, and the orbit error and clock error of the satellite and the hardware delay of the satellite end can be corrected by differential data. In addition, in order to improve the accuracy of the solution result, errors such as tropospheric dry delay, phase winding, antenna phase center offset, and relativistic effect can be corrected in advance through the set error correction model.

Optionally, in this formula, ε _I,j and ε _T,j are small, so It can be regarded as approximately equal to It can be considered as approximately equal to T.

Step S21-2, performing wide lane ambiguity fixation on the phase floating point ambiguity of each satellite to obtain the wide lane fixed ambiguity of each satellite;

Step S21-3, fix the narrow lane ambiguity of the wide lane fixed ambiguity of each satellite to obtain a positioning result.

Optionally, after solving the unknown data, the wide lane ambiguity of the phase floating point ambiguity of each satellite may be fixed to obtain the wide lane fixed ambiguity. Specifically, based on FIG. 3 , FIG. 4 is another flow chart of the accuracy marking method for satellite positioning provided by an embodiment of the present disclosure. Please refer to FIG. 4 . The above step S21-2 may be implemented by the following steps:

Step S21-2-1, selecting a reference satellite from all satellites according to a preset rule, and using other satellites other than the reference satellite as mobile satellites;

Step S21-2-2, calculating the corresponding inter-satellite single-difference wide-lane floating ambiguity according to the phase floating ambiguity corresponding to each frequency of each mobile satellite;

Optionally, the preset rule may be to randomly select a reference satellite, or to select a reference satellite according to a satellite altitude angle.

In this embodiment, the receiving device can select one from all satellites as a reference satellite and use the other satellites as mobile satellites, so that for each mobile satellite, the inter-satellite single-difference wide-lane floating-point ambiguity of the satellite is calculated according to the phase floating-point ambiguity corresponding to each frequency of the satellite.

Optionally, the corresponding inter-satellite single-difference wide-lane floating ambiguity may be calculated according to the phase floating ambiguity corresponding to each frequency of each mobile satellite and the phase floating ambiguity corresponding to each frequency of the reference satellite.

Specifically, for each mobile satellite, the inter-satellite single-difference wide-lane floating point ambiguity is calculated by the following formula:
WL＝N _y,1 -N _y,2 -(N _c,1 -N _C,2 )

Among them, WL represents the inter-satellite single-difference wide-lane floating ambiguity, N _y,1 represents the phase floating ambiguity of the mobile satellite at the first frequency, N _y,2 represents the phase floating ambiguity of the mobile satellite at the second frequency, N _c,1 represents the phase floating ambiguity of the reference satellite at the first frequency, and N _C,2 represents the phase floating ambiguity of the reference satellite at the second frequency.

Step S21-2-3, if the inter-satellite single-difference wide-lane floating-point ambiguity of each mobile satellite meets the preset conditions, then the integer wide-lane ambiguity corresponding to the mobile satellite is calculated according to the inter-satellite single-difference wide-lane floating-point ambiguity;

Optionally, the preset condition may be that the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity is less than a preset threshold. On this basis, for each mobile satellite, if the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity is less than the preset threshold, the inter-satellite single-difference wide-lane floating-point ambiguity is rounded to obtain an integer wide-lane ambiguity.

Optionally, the preset threshold can be 0.3. It can be understood that if the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity of a mobile satellite is greater than or equal to 0.3, it is determined that the mobile satellite does not meet the preset conditions, and there is no need to round off the satellite that does not meet the preset conditions.

Optionally, the inter-satellite single-difference wide-lane float ambiguity may be rounded by the following formula: WL ⁰ =round(N _y,1 -N _y,2 -(N _c,1 -N _C,2 )).

Among them, WL ⁰ represents the integer wide lane ambiguity, and round() represents the rounding operation on the data in the brackets.

Step S21-2-4, updating the constrained ambiguity parameters according to the integer wide lane ambiguity as the inter-satellite single-difference wide lane floating point ambiguity of the mobile satellite, and obtaining the wide lane fixed ambiguity of the mobile satellite.

Optionally, the constrained ambiguity parameters can be updated for the inter-satellite single-difference wide-lane float ambiguities of the mobile satellite by the following formula:
WL ⁰ = WL + ε

Optionally, after obtaining the wide lane fixed ambiguity of each satellite, it is also necessary to perform narrow lane ambiguity fixation on the wide lane fixed ambiguity of each satellite to obtain a positioning result. Specifically, based on FIG. 3 , FIG. 5 is another flow chart of the accuracy marking method for satellite positioning provided by an embodiment of the present disclosure. Please refer to FIG. 5 . The above step S21-3 can be implemented by the following steps:

Step S21-3-1, performing narrow lane ambiguity fixation on wide lane fixed ambiguities of multiple satellites according to a preset search algorithm to obtain inter-satellite single difference integer narrow lane ambiguities and ratio values;

Optionally, the preset search algorithm may be a lambda search algorithm.

In one possible implementation manner, the wide lane fixed ambiguities of the multiple satellites may include only the wide lane fixed ambiguities of satellites that meet a preset condition in the wide lane ambiguity fixation process; in another possible implementation manner, the wide lane fixed ambiguities of the multiple satellites may include the wide lane fixed ambiguities of satellites that meet the preset condition in the wide lane ambiguity fixation process, and the wide lane fixed ambiguities of satellites that do not meet the preset condition in the wide lane ambiguity fixation process.

Step S21-3-2, when the ratio value does not reach the preset ratio, the target wide lane fixed ambiguity is deleted from the wide lane fixed ambiguities of the multiple satellites, and the narrow lane ambiguity is fixed for the other wide lane fixed ambiguities again according to the preset search algorithm;

The target wide lane fixed ambiguity is the wide lane fixed ambiguity corresponding to the satellite with the lowest satellite elevation angle among multiple satellites. Fuzziness;

Optionally, the preset ratio may be 2.5.

In this embodiment, after the receiving device performs narrow lane ambiguity fixation on the wide lane fixed ambiguities of multiple satellites through the lambda search algorithm, the receiving device can obtain the inter-satellite single difference integer narrow lane ambiguity and the ratio value, and the receiving device can determine whether the ratio value reaches a preset ratio. If it does not reach the preset ratio, the receiving device can delete the wide lane fixed ambiguity corresponding to the satellite with the lowest satellite elevation angle from all wide lane fixed ambiguities, and perform a lambda search again based on the other wide lane fixed ambiguities from which the wide lane fixed ambiguity corresponding to the satellite with the lowest satellite elevation angle has been deleted.

Step S21-3-3, if the ratio value reaches a preset ratio, and the number of remaining wide lane fixed ambiguities is greater than or equal to a preset number, it is determined that the narrow lane ambiguity is fixed successfully, and the inter-satellite single difference integer narrow lane ambiguity is constrained to update the inter-satellite single difference narrow lane floating point ambiguity to obtain a positioning result;

Among them, the inter-satellite single-difference narrow-lane floating-point ambiguity is obtained by solving the observation equation.

Optionally, the receiving device can obtain the inter-satellite single-difference narrow-lane floating-point ambiguity when solving the observation equation. Therefore, when performing narrow-lane ambiguity fixation on the wide-lane fixed ambiguity, if it is determined that the narrow-lane ambiguity fixation is successful based on the ratio value and the number of remaining wide-lane fixed ambiguities, the receiving device can update the constraints for the inter-satellite single-difference narrow-lane floating-point ambiguity based on the obtained inter-satellite single-difference integer narrow-lane ambiguity to obtain the positioning result.

Optionally, the preset number can be set according to actual conditions. In one possible implementation, the preset number can be 4.

It can be understood that the updating of constraints for inter-satellite single-difference narrow-lane floating-point ambiguities according to the integer narrow-lane ambiguities is the same as the above-mentioned step of updating the constraint ambiguity parameters for the inter-satellite single-difference wide-lane floating-point ambiguities of the mobile satellite according to the integer wide-lane ambiguities, and will not be repeated here.

Step S21-3-4, if the ratio value does not reach the preset ratio, or the number of wide lane fixed ambiguities for narrow lane ambiguity fixation is less than the preset number, it is determined that the narrow lane ambiguity fixation has failed, and the positioning result is obtained.

In this embodiment, if the ratio value does not reach the preset ratio, or the number of wide lane fixed ambiguities for narrow lane ambiguity fixation is less than the preset number, it means that the narrow lane ambiguity fixation fails at this time, so there is no need to update the constraints for the inter-satellite single-difference narrow lane floating point ambiguities according to the inter-satellite single-difference integer narrow lane ambiguities.

Optionally, after obtaining the positioning result, the ambiguity fixing processing state of each step in the above ambiguity fixing process can be determined, so as to determine the accuracy of the positioning result finally obtained. Optionally, the receiving device can first determine whether the observation equation is solved successfully. If the observation equation fails to be solved, the ambiguity fixing processing state can be determined to be the sixth state, and the accuracy of the positioning result can be determined to be the sixth accuracy according to the sixth state.

In one possible implementation, it can be determined whether the observation equation is successfully solved by determining whether all unknown parameters can be completely solved. It can be understood that if all unknown parameters cannot be completely solved, it means that the observation equation has failed to be solved.

Optionally, the accuracy of the positioning result calculated in the sixth state is the worst. In a possible implementation, the sixth accuracy can represent an accuracy level of 10.0 meters.

Optionally, if the receiving device determines that the observation equation is solved successfully, it can further determine whether the coordinate variance of the receiving device solved by the observation equation is greater than or equal to a preset variance. That is, when the observation equation is solved successfully and the coordinate variance of the receiving device is greater than or equal to the preset variance, the ambiguity fixed processing state is determined to be the fifth state, and the accuracy of the positioning result is determined to be the fifth accuracy based on the fifth state.

Optionally, the coordinates of the receiving device may be obtained by solving the observation equation, and the receiving device may calculate the coordinate variance of the receiving device according to the coordinates, and determine whether the coordinate variance is greater than or equal to a preset variance.

In a possible implementation, the preset variance may be 0.5.

Optionally, the accuracy of the positioning result obtained in the fifth state is better than the accuracy of the positioning result obtained in the sixth state. In a possible implementation, the fifth accuracy can represent an accuracy level of 5.0 rice.

It can be understood that the fifth precision is the default precision of the floating-point solution in the prior art.

Optionally, if the coordinate variance of the receiving device is less than a preset variance, it can be further determined whether the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than or equal to a preset number, and whether the narrow lane ambiguity is fixed successfully.

Specifically, when the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide-lane ambiguity is fixed is less than the preset number, and the narrow-lane ambiguity fixation fails, the ambiguity fixation processing state is determined to be the fourth state; according to the fourth state, the accuracy of the positioning result is determined to be the fourth accuracy.

Optionally, the accuracy of the positioning result obtained by solving in the fourth state is better than the accuracy of the positioning result obtained by solving in the fifth state. In a possible implementation, the fourth accuracy can represent an accuracy level of 3.0 meters.

When the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than the preset number, and the narrow lane ambiguity is fixed successfully, the ambiguity fixing processing state is determined to be the third state; according to the third state, the accuracy of the positioning result is determined to be the third accuracy.

Optionally, the accuracy of the positioning result obtained by solving in the third state is better than the accuracy of the positioning result obtained by solving in the fourth state. In a possible implementation, the third accuracy can represent an accuracy level of 0.5 meters.

When the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity fixation fails, the ambiguity fixation processing state is determined to be the second state; and the accuracy of the positioning result is determined to be the second accuracy according to the second state.

Optionally, the accuracy of the positioning result obtained by solving in the second state is better than the accuracy of the positioning result obtained by solving in the third state. In a possible implementation, the second accuracy can represent an accuracy level of 0.3 meters.

When the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity is fixed successfully, the ambiguity fixing processing state is determined to be the first state; and the accuracy of the positioning result is determined to be the first accuracy according to the first state.

Optionally, the accuracy of the positioning result obtained by solving in the first state is better than the accuracy of the positioning result obtained by solving in the second state. In a possible implementation, the first accuracy can represent an accuracy level of 0.1 meters.

It can be understood that the first accuracy is the accuracy of the fixed solution.

In this embodiment, if it is determined in advance that the pseudo-range single point positioning fails, the positioning result is an invalid solution.

Optionally, referring to Figure 6, the present disclosure can divide the positioning results into the following 7 levels, namely level 0, level 1, level 2, level 3, level 4, level 5, level 6 and level 7, wherein the accuracy level corresponding to level 0 is an invalid solution, the accuracy level corresponding to level 1 is 0.1 meter, the accuracy level corresponding to level 2 is 0.3 meter, the accuracy level corresponding to level 3 is 0.5 meter, the accuracy level corresponding to level 4 is 3.0 meters, the accuracy level corresponding to level 5 is 5.0 meters, and the accuracy level corresponding to level 6 is 10.0 meters.

Optionally, in order to facilitate the user to determine whether the accuracy of the obtained positioning results meets the positioning accuracy requirements of the current application scenario, the receiving device can present the accuracy grade and accuracy level of the obtained positioning results to the user, so that the user can choose whether to select the positioning results at this accuracy level.

For example, if it is determined that the obtained positioning result is of the second accuracy, the receiving device may display that the accuracy level of the current positioning result is level 2, and the accuracy level is 0.3 meters.

In order to execute the corresponding steps in the above embodiments and various possible methods, a method for implementing a satellite positioning accuracy marking device is provided below. Functional module diagram of the positioning accuracy marking device. It should be noted that the basic principle and technical effect of the satellite positioning accuracy marking device provided in this embodiment are the same as those of the above embodiments. For the sake of brief description, for the parts not mentioned in this embodiment, reference can be made to the corresponding contents in the above embodiments. The satellite positioning accuracy marking device includes: a receiving module 200, a processing module 210 and a determining module 220.

The receiving module 200 is used to receive satellite observation data and differential data;

It can be understood that the receiving module 200 can also be used to perform the above step S20;

The processing module 210 is used to perform ambiguity fixation processing according to the observation data and the differential data to obtain a positioning result;

It can be understood that the processing module 210 can also be used to perform the above step S21;

The determination module 220 is used to determine the ambiguity fixation processing state, and determine the accuracy of the positioning result according to the ambiguity fixation processing state.

It can be understood that the determination module 220 can also be used to execute the above step S22.

Optionally, the determination module 220 is also used to determine that the ambiguity fixing processing state is the first state when the observation equation is solved successfully, the coordinate variance of the receiving device is less than a preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to a preset number, and the narrow lane ambiguity is fixed successfully; and determine the accuracy of the positioning result as the first accuracy according to the first state.

Optionally, the determination module 220 is also used to determine that the ambiguity fixing processing state is the second state when the observation equation is solved successfully, the coordinate variance of the receiving device is less than a preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to a preset number, and the narrow lane ambiguity fixation fails; and determine the accuracy of the positioning result as the second accuracy according to the second state.

Optionally, the determination module 220 is also used to determine that the ambiguity fixing processing state is a third state when the observation equation is solved successfully, the coordinate variance of the receiving device is less than a preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than a preset number, and the narrow lane ambiguity is fixed successfully; and determine the accuracy of the positioning result as the third accuracy according to the third state.

Optionally, the determination module 220 is also used to determine that the ambiguity fixing processing state is a fourth state when the observation equation is solved successfully, the coordinate variance of the receiving device is less than a preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than a preset number, and the narrow lane ambiguity fixation fails; and determine the accuracy of the positioning result as the fourth accuracy according to the fourth state.

Optionally, the determination module 220 is also used to determine that the ambiguity fixing processing state is the fifth state when the observation equation is solved successfully and the coordinate variance of the receiving device is greater than or equal to the preset variance; and determine the accuracy of the positioning result as the fifth accuracy according to the fifth state.

Optionally, the determination module 220 is further configured to, when the observation equation fails to be solved, determine that the ambiguity fixing processing state is a sixth state; and determine that the accuracy of the positioning result is a sixth accuracy according to the sixth state.

Optionally, the determination module 220 is further configured to, when the pseudorange single-point positioning fails, determine that the ambiguity fixing processing state is a seventh state; and determine that the positioning result is an invalid solution according to the seventh state.

Optionally, the processing module 210 is also used to, when the pseudorange single-point positioning is successful, construct an observation equation based on the observation data and the differential data, and calculate the phase floating-point ambiguity of each satellite based on the observation equation; perform wide-lane ambiguity fixation on the phase floating-point ambiguity of each satellite to obtain the wide-lane fixed ambiguity of each satellite; and perform narrow-lane ambiguity fixation on the wide-lane fixed ambiguity of each satellite to obtain a positioning result.

It can be understood that the processing module 210 can also be used to execute the above steps S21 - 1 to S21 - 3.

Optionally, the processing module 210 is further configured to construct an observation equation according to the following formula when the pseudorange single point positioning is successful:

Where P _s,j represents the pseudorange of satellite s at frequency j; L _s,j represents the phase observation value of satellite s at frequency j. represents the distance between the antenna phase center of the receiving device r and the phase center of the satellite s, C represents the speed of light in a vacuum, dt _r,j represents the clock error of the receiving device r at frequency j, dt ^s represents the clock error of the satellite s, T represents the wet tropospheric delay, γ represents the ratio between the squares of multiple frequencies, represents the ionospheric delay of satellite s at frequency j, b _r,j represents the pseudorange hardware delay of receiving device r at frequency j, represents the pseudorange hardware delay of satellite s at frequency j, B _r,j represents the phase hardware delay of receiving device r at frequency j, Characterizes the phase hardware delay of satellite s at frequency j, Characterizes the wavelength of satellite s at frequency j, represents the phase floating point ambiguity of satellite s at frequency j, represents the ionospheric delay of the receiving device, ε _I,j represents the observation noise of the pseudorange, represents the tropospheric delay of the vehicle terminal, and ε _T,j represents the observation noise of the phase observation value.

Optionally, the processing module 210 is further used to select a reference satellite from all satellites according to a preset rule, and use other satellites other than the reference satellite as mobile satellites; calculate the corresponding inter-satellite single-difference wide-lane floating-point ambiguity according to the phase floating-point ambiguity corresponding to each frequency of each mobile satellite; if the inter-satellite single-difference wide-lane floating-point ambiguity of each mobile satellite meets a preset condition, calculate the integer wide-lane ambiguity corresponding to the mobile satellite according to the inter-satellite single-difference wide-lane floating-point ambiguity; update the constrained ambiguity parameters for the inter-satellite single-difference wide-lane floating-point ambiguity of the mobile satellite according to the integer wide-lane ambiguity, and obtain the wide-lane fixed ambiguity of the mobile satellite.

It can be understood that the processing module 210 can also be used to execute the above steps S21-2-1 to S21-2-4.

Optionally, the processing module 210 is further used to calculate corresponding inter-satellite single-difference wide-lane floating-point ambiguities according to the phase floating-point ambiguities corresponding to each frequency of each mobile satellite and the phase floating-point ambiguities corresponding to each frequency of the reference satellite.

Optionally, the processing module 210 is further configured to calculate, for each mobile satellite, an inter-satellite single-difference wide-lane floating point ambiguity by the following formula:
WL＝N _y,1 -N _y,2 -(N _c,1 -N _C,2 )

Optionally, the processing module 210 is further configured to, for each mobile satellite, round off the inter-satellite single-difference wide-lane floating-point ambiguity to obtain an integer wide-lane ambiguity if the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity is less than a preset threshold.

Optionally, the processing module 210 is further configured to update the constrained ambiguity parameter for the inter-satellite single-difference wide-lane floating point ambiguity of the mobile satellite by using the following formula:
WL ⁰ = WL + ε

Optionally, the processing module 210 is further configured to perform narrow lane ambiguity fixation on the wide lane fixed ambiguities of the plurality of satellites according to a preset search algorithm to obtain inter-satellite single difference integer narrow lane ambiguity and a ratio value; if the ratio value does not reach a preset ratio, delete the target wide lane fixed ambiguity from the wide lane fixed ambiguities of the plurality of satellites, and perform narrow lane ambiguity fixation on the other wide lane fixed ambiguities again according to a preset search algorithm; wherein the target wide lane fixed ambiguity is the wide lane fixed ambiguity corresponding to the satellite with the lowest satellite elevation angle among the plurality of satellites. The method comprises the steps of: determining an ambiguity by fixing the narrow lane ambiguity; if the ratio value reaches a preset ratio, and the number of remaining wide lane fixed ambiguities is greater than or equal to a preset number, determining that the narrow lane ambiguity is fixed successfully, and updating the inter-satellite single-difference narrow lane floating point ambiguity as a constraint for inter-satellite single-difference narrow lane floating point ambiguity according to the inter-satellite single-difference integer narrow lane ambiguity, to obtain the positioning result; wherein the inter-satellite single-difference narrow lane floating point ambiguity is obtained by solving the observation equation; if the ratio value does not reach the preset ratio, or the number of wide lane fixed ambiguities for narrow lane ambiguity fixation is less than a preset number, determining that the narrow lane ambiguity fixation fails, and obtaining the positioning result.

It can be understood that the processing module 210 can also be used to execute the above steps S21-3-1 to S21-3-4.

The precision marking device for satellite positioning provided by the embodiment of the present disclosure receives the observation data and differential data of the satellite through the receiving module; performs ambiguity fixing processing according to the observation data and differential data through the processing module to obtain the positioning result; determines the ambiguity fixing processing state through the determining module, and determines the precision of the positioning result according to the ambiguity fixing processing state. After performing ambiguity fixing processing according to the observation data and differential data, the device can determine the precision of the positioning result according to the ambiguity fixing state, so no matter whether the obtained positioning result is a fixed solution or a floating point solution, the precision of the positioning result can be determined according to the ambiguity fixing state.

In several embodiments provided in the present disclosure, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely schematic. For example, the flowcharts and block diagrams in the accompanying drawings show the possible architecture, functions and operations of the devices, methods and computer program products according to multiple embodiments of the present disclosure. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of the code, and a module, a program segment or a part of the code contains one or more executable instructions for implementing the specified logical functions. It should also be noted that in some alternative implementations, the functions marked in the box can also occur in a different order from the order marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram and/or flowchart, and the combination of boxes in the block diagram and/or flowchart can be implemented with a dedicated hardware-based system that performs a specified function or action, or can be implemented with a combination of dedicated hardware and computer instructions.

In addition, the functional modules in the various embodiments of the present disclosure may be integrated together to form an independent part, or each module may exist separately, or two or more modules may be integrated to form an independent part.

If the function is implemented in the form of a software function module and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present disclosure, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present disclosure. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, and other media that can store program codes.

The above description is only a preferred embodiment of the present disclosure and is not intended to limit the present disclosure. For those skilled in the art, the present disclosure may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present disclosure shall be included in the protection scope of the present disclosure.

Industrial Applicability:

The present disclosure provides a satellite positioning accuracy marking method and related devices, which can determine the accuracy of the positioning result according to the ambiguity fixing state after performing ambiguity fixing processing according to observation data and differential data. Therefore, no matter whether the obtained positioning result is a fixed solution or a floating point solution, the accuracy of the positioning result can be determined according to the ambiguity fixing state.

Claims

A satellite positioning accuracy marking method, characterized in that the method comprises:

Receive satellite observation data and differential data;

Performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result;

An ambiguity fixation processing state is determined, and an accuracy of the positioning result is determined according to the ambiguity fixation processing state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity is fixed successfully, determining that the ambiguity fixing processing state is the first state;

The accuracy of the positioning result is determined to be a first accuracy according to the first state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is greater than or equal to the preset number, and the narrow lane ambiguity fixation fails, determining that the ambiguity fixation processing state is the second state;

The accuracy of the positioning result is determined to be a second accuracy according to the second state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than the preset number, and the narrow lane ambiguity is fixed successfully, determining that the ambiguity fixing processing state is the third state;

The accuracy of the positioning result is determined to be a third accuracy according to the third state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In the case where the observation equation is solved successfully, the coordinate variance of the receiving device is less than the preset variance, the number of satellites that meet the preset conditions after the wide lane ambiguity is fixed is less than the preset number, and the narrow lane ambiguity fixation fails, determining that the ambiguity fixation processing state is the fourth state;

The accuracy of the positioning result is determined to be a fourth accuracy according to the fourth state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

When the observation equation is solved successfully and the coordinate variance of the receiving device is greater than or equal to the preset variance, determining that the ambiguity fixing processing state is the fifth state;

The accuracy of the positioning result is determined to be a fifth accuracy according to the fifth state.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In case of failure in solving the observation equation, determining the ambiguity fixing processing state to be the sixth state;

The accuracy of the positioning result is determined according to the sixth state as a sixth accuracy.
The method according to claim 1, characterized in that the determining of the ambiguity fixing processing state and determining the accuracy of the positioning result according to the ambiguity fixing processing state comprises:

In the case of pseudorange single point positioning failure, determining the ambiguity fix processing state to be the seventh state;

The positioning result is determined to be an invalid solution according to the seventh state.
The method according to claim 1, characterized in that the performing ambiguity fixing processing according to the observation data and the differential data to obtain the positioning result comprises:

When the pseudo-range single point positioning is successful, constructing an observation equation according to the observation data and the differential data, and calculating the phase floating point ambiguity of each satellite according to the observation equation;

Performing wide lane ambiguity fixation on the phase floating point ambiguity of each of the satellites to obtain wide lane fixed ambiguity of each of the satellites;

The narrow lane ambiguity is fixed for the wide lane fixed ambiguity of each satellite to obtain a positioning result.
The method according to claim 9 is characterized in that, when the pseudorange single point positioning is successful, constructing the observation equation according to the observation data and the differential data comprises:

When the pseudo-range single point positioning is successful, the observation equation is constructed according to the following formula:

Where P s,j represents the pseudorange of satellite s at frequency j; L s,j represents the phase observation value of satellite s at frequency j. represents the distance between the antenna phase center of the receiving device r and the phase center of the satellite s, C represents the speed of light in a vacuum, dt r,j represents the clock error of the receiving device r at frequency j, dt s represents the clock error of satellite s, T represents the wet tropospheric delay, γ represents the ratio between the squares of multiple frequencies, represents the ionospheric delay of satellite s at frequency j, b r,j represents the pseudorange hardware delay of receiving device r at frequency j, represents the pseudorange hardware delay of satellite s at frequency j, B r,j represents the phase hardware delay of receiving device r at frequency j, Characterizes the phase hardware delay of satellite s at frequency j, Characterizes the wavelength of satellite s at frequency j, represents the phase floating point ambiguity of satellite s at frequency j, represents the ionospheric delay of the receiving device, ε I,j represents the observation noise of the pseudorange, represents the tropospheric delay of the vehicle terminal, and ε T,j represents the observation noise of the phase observation value.
The method according to claim 9, characterized in that the performing wide lane ambiguity fixing on the phase floating point ambiguity of each satellite to obtain the wide lane fixed ambiguity of each satellite comprises:

Selecting a reference satellite from all the satellites according to a preset rule, and using the other satellites except the reference satellite as mobile satellites;

Calculating the corresponding inter-satellite single-difference wide-lane floating ambiguity according to the phase floating ambiguity corresponding to each frequency of each mobile satellite;

If the inter-satellite single-difference wide-lane floating-point ambiguity of each of the mobile satellites meets a preset condition, then the integer wide-lane ambiguity corresponding to the mobile satellite is calculated according to the inter-satellite single-difference wide-lane floating-point ambiguity;

The constrained ambiguity parameters are updated as the inter-satellite single-difference wide-lane floating-point ambiguity of the mobile satellite according to the integer wide-lane ambiguity to obtain the wide-lane fixed ambiguity of the mobile satellite.
The method according to claim 11 is characterized in that the step of calculating the corresponding inter-satellite single-difference wide-lane floating-point ambiguity according to the phase floating-point ambiguity corresponding to each frequency of each of the mobile satellites comprises:

The corresponding inter-satellite single-difference wide-lane floating ambiguity is calculated according to the phase floating ambiguity corresponding to each frequency of each mobile satellite and the phase floating ambiguity corresponding to each frequency of the reference satellite.
The method according to claim 12, characterized in that the calculating the corresponding inter-satellite single-difference wide-lane floating ambiguity according to the phase floating ambiguity corresponding to each frequency of each of the mobile satellites and the phase floating ambiguity corresponding to each frequency of the reference satellite comprises:

For each of the mobile satellites, the inter-satellite single-difference wide-lane floating point ambiguity is calculated by the following formula:
WL＝N y,1 -N y,2 -(N c,1 -N c,2 )

Among them, WL represents the inter-satellite single-difference wide-lane floating-point ambiguity, N y,1 represents the phase floating-point ambiguity of the mobile satellite at the first frequency, N y,2 represents the phase floating-point ambiguity of the mobile satellite at the second frequency, N c,1 represents the phase floating-point ambiguity of the reference satellite at the first frequency, and N c,2 represents the phase floating-point ambiguity of the reference satellite at the second frequency.
The method according to claim 11 is characterized in that if the inter-satellite single-difference wide-lane floating-point ambiguity of each of the mobile satellites meets a preset condition, then calculating the integer wide-lane ambiguity corresponding to the mobile satellite according to the inter-satellite single-difference wide-lane floating-point ambiguity comprises:

For each of the mobile satellites, if the decimal part of the inter-satellite single-difference wide-lane floating-point ambiguity is less than a preset threshold, the inter-satellite single-difference wide-lane floating-point ambiguity is rounded to an integer to obtain the integer wide-lane ambiguity.
The method according to claim 11 is characterized in that the updating of constrained ambiguity parameters for the inter-satellite single-difference widelane floating-point ambiguities of the mobile satellite according to the integer widelane ambiguities to obtain the widelane fixed ambiguities of the mobile satellite comprises:

The constrained ambiguity parameters are updated for the inter-satellite single-difference wide-lane floating point ambiguity of the mobile satellite by the following formula:
WL 0 = WL + ε

Among them, WL 0 represents the integer wide lane ambiguity, WL represents the inter-satellite single-difference wide lane floating point ambiguity, and ε represents the error value.
The method according to claim 9, characterized in that the step of performing narrow lane ambiguity fixing on the wide lane fixed ambiguity of each satellite to obtain a positioning result comprises:

Fixing narrow lane ambiguities of wide lane fixed ambiguities of the plurality of satellites according to a preset search algorithm to obtain inter-satellite single difference integer narrow lane ambiguities and ratio values;

When the ratio value does not reach a preset ratio, a target wide lane fixed ambiguity is deleted from the wide lane fixed ambiguities of the plurality of satellites, and narrow lane ambiguity is fixed again for the other wide lane fixed ambiguities according to a preset search algorithm; wherein the target wide lane fixed ambiguity is a wide lane fixed ambiguity corresponding to a satellite with the lowest satellite elevation angle among the plurality of satellites;

If the ratio value reaches a preset ratio, and the number of remaining wide lane fixed ambiguities is greater than or equal to a preset number, it is determined that the narrow lane ambiguity is fixed successfully, and the inter-satellite single-difference integer narrow lane ambiguity is used as an inter-satellite single-difference narrow lane floating point ambiguity update constraint to obtain the positioning result; wherein the inter-satellite single-difference narrow lane floating point ambiguity is obtained by solving the observation equation;

If the ratio value does not reach a preset ratio, or the number of wide lane fixed ambiguities for narrow lane ambiguity fixation is less than a preset number, it is determined that the narrow lane ambiguity fixation fails, and the positioning result is obtained.
A satellite positioning accuracy marking device, characterized in that the device comprises:

A receiving module, used for receiving satellite observation data and differential data;

A processing module, used for performing ambiguity fixing processing according to the observation data and the differential data to obtain a positioning result;

The determination module is used to determine an ambiguity fixing processing state, and determine the accuracy of the positioning result according to the ambiguity fixing processing state.
A receiving device, characterized in that it includes a processor and a memory, wherein the memory stores a computer program that can be executed by the processor, and the processor can execute the computer program to implement any method described in claims 1-16.
A computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the method according to any one of claims 1 to 16.