CN112305574A - Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method - Google Patents

Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method Download PDF

Info

Publication number
CN112305574A
CN112305574A CN202010434736.5A CN202010434736A CN112305574A CN 112305574 A CN112305574 A CN 112305574A CN 202010434736 A CN202010434736 A CN 202010434736A CN 112305574 A CN112305574 A CN 112305574A
Authority
CN
China
Prior art keywords
difference
observation
satellite
value
pseudo
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202010434736.5A
Other languages
Chinese (zh)
Inventor
赵光静
邵炜平
杨鸿珍
由奇林
娄竞
娄佳
李贤�
贺家乐
王甜甜
王文龙
任白杨
彭卉
龙强
蓝天
石帅
杨怀丽
李宇翔
陈端云
陈功伯
林琳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Siji Shenwang Position Service Beijing Co ltd
State Grid Corp of China SGCC
State Grid Information and Telecommunication Co Ltd
State Grid Zhejiang Electric Power Co Ltd
State Grid Fujian Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Zhejiang Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Jibei Electric Power Co Ltd
Original Assignee
State Grid Siji Shenwang Position Service Beijing Co ltd
State Grid Corp of China SGCC
State Grid Information and Telecommunication Co Ltd
State Grid Zhejiang Electric Power Co Ltd
State Grid Fujian Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Zhejiang Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Jibei Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Siji Shenwang Position Service Beijing Co ltd, State Grid Corp of China SGCC, State Grid Information and Telecommunication Co Ltd, State Grid Zhejiang Electric Power Co Ltd, State Grid Fujian Electric Power Co Ltd, Information and Telecommunication Branch of State Grid Zhejiang Electric Power Co Ltd, Information and Telecommunication Branch of State Grid Jibei Electric Power Co Ltd filed Critical State Grid Siji Shenwang Position Service Beijing Co ltd
Priority to CN202010434736.5A priority Critical patent/CN112305574A/en
Publication of CN112305574A publication Critical patent/CN112305574A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/256Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to timing, e.g. time of week, code phase, timing offset
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/24Acquisition or tracking or demodulation of signals transmitted by the system
    • G01S19/25Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS
    • G01S19/258Acquisition or tracking or demodulation of signals transmitted by the system involving aiding data received from a cooperating element, e.g. assisted GPS relating to the satellite constellation, e.g. almanac, ephemeris data, lists of satellites in view

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention relates to a Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method, which comprises the steps of inputting carrier phase observation values, pseudo-range observation values and broadcast ephemeris of a plurality of Beidou/GNSS satellites of an observation station 1 and an observation station 2; respectively calculating the coordinates of each satellite and the clock correction by using a satellite position calculation standard algorithm; and (3) comprehensively utilizing a pseudo-range positioning residual analysis method, an inter-station single-difference GF and MW combined observation value method and a high-order difference method to detect pseudo-range gross error and carrier phase cycle slip to calculate to obtain each data, and carrying out cycle slip judgment and pseudo-range positioning pre-and post-test residual judgment to finally obtain observation data with high effectiveness and reliability. The invention comprehensively utilizes a pseudo-range positioning residual error analysis method, an inter-station single-difference GF and MW combined observation value method and a high-order difference method to detect pseudo-range gross error and carrier phase cycle slip, thereby improving the effectiveness and reliability of detection.

Description

Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method
Technical Field
The invention relates to a Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method, and belongs to the technical field of satellite positioning algorithms.
Background
With the development of technologies such as unmanned aerial vehicles and automatic driving, the application of Beidou/GNSS high-precision real-time dynamic positioning and orientation is very wide, and higher requirements are provided for guaranteeing the safety of unmanned aerial vehicles and automatic driving and the precision and reliability of real-time dynamic positioning and orientation, wherein the effect of real-time data preprocessing is to obtain the basis and the premise of high-precision and high-reliability real-time positioning and orientation. At present, in real-time high-precision dynamic positioning and orientation application, a pseudo-range coarse difference detection method mainly adopts a pseudo-range robust least square estimation method to detect a coarse difference observed value in a pseudo-range observed value, and cycle slip detection generally detects cycle slip existing in a phase observed value by using the assumed characteristic that a Geometry-free (GF) and a Melbourne-Wubbena (MW) combined observed value slowly and smoothly changes along with time. In the past, Beidou/GNSS high-precision positioning and orientation application generally requires less clearance shielding of an observation environment, and under the condition, the method can obtain a good data preprocessing effect. However, with the appearance of the application of unmanned aerial vehicles and automatic driving, the real-time dynamic positioning and orientation cannot meet the requirements of the observation environment easily, and under the complex environment with serious shielding, the method has obvious phenomena of misdetection and missed detection, which can seriously affect the precision and reliability of the real-time dynamic positioning and orientation of the Beidou/GNSS and reduce the application performance of the unmanned aerial vehicles and the automatic driving technology.
Disclosure of Invention
The invention mainly overcomes the defects in the prior art, provides a Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method, and the system and method comprehensively utilize a pseudo-range positioning residual error analysis method, an inter-station single-difference GF and MW combined observation value method and a high-order difference method to detect pseudo-range gross error and carrier phase cycle slip, thereby improving the detection effectiveness and reliability.
The technical scheme provided by the invention for solving the technical problems is as follows: the Beidou GNSS satellite real-time positioning and orientation data preprocessing method comprises the following steps:
step S1, inputting carrier phase observation values, pseudo-range observation values and broadcast ephemeris of a plurality of Beidou/GNSS satellites of the observation station 1 and the observation station 2;
step S2, calculating the coordinates and the correction of clock error of all the satellites input in the step S1 through a satellite position calculation standard algorithm according to the broadcast ephemeris, deleting the satellite data without the broadcast ephemeris product, and deleting the data of unhealthy satellites according to the satellite health identification provided by the broadcast ephemeris;
step S3, judging pseudo-range observation values of all satellites input in step S1 by using a pseudo-range positioning pre-test residual error judgment method and a pseudo-range positioning pre-test residual error judgment method to delete the pseudo-range observation values of the satellites with gross errors;
step S4, marking whether cycle slip occurs on the carrier phase observation values of all the satellites input in the step S1 by using a single difference GF and MW combined observation value method and a high-order difference value method;
and step S5, finally obtaining a pseudo-range observation value without gross error, a carrier phase observation value marked as cycle slip occurrence, and a carrier phase observation value marked as cycle slip non-occurrence.
The further technical scheme is that the specific process of the step S3 is as follows:
step S31, calculating double-difference pseudo-range observed values of all satellites by pseudo-range positioning residual error analysis method by using coordinates and clock error correction numbers of all satellites, and setting initial parameter value X0The pseudo-range double-difference observation equation is linearized to obtain the pre-test residual v of each double-difference pseudo-range observation value of all satellites0(ii) a If | v0|>30, if the gross error exists, rejecting a pseudo range observation value of the satellite; | v0If the absolute value is less than or equal to 30, considering that no gross error exists, and carrying out the next step;
step S32, obtaining parameter estimation value by least square estimation
Figure BDA0002501805820000021
The baseline distance constraints are as follows:
const=sqrt(N2+E2+U2)
wherein N, E, U respectively represent parameter vectors
Figure BDA0002501805820000031
North, middle, east, and high, const representing the baseline length;
step S33, respectively calculating the post-test residual v and the RMS value delta of the post-test residual of each satellite; if the absolute value v is larger than 4 multiplied by delta, considering that gross errors exist, rejecting the observation value of the satellite, and if the absolute value v is less than or equal to 4 multiplied by delta, considering that no gross errors exist, and keeping the pseudo-range observation value of the satellite;
in a further technical solution, the specific process of step S4 is:
step S41, calculating single difference GF and MW combined observation value delta phi through an inter-station single difference GF and MW combined observation value method according to carrier phase observation values of all satellitesGF、ΔLMW
Step S42, respectively carrying out filtering estimation on the single-difference GF and MW combined observed values obtained in step 41 by adopting static successive filtering, weakening the influence of noise, and obtaining the estimated value epsilon (phi) of the single-difference GF and MW combined observed valuesGF)、ε(LMW) And estimating the variance δ (Φ)GF)、δ(LMW);
Combining single difference GF and MW of current epoch into observed value delta phiGF、ΔLMWEstimate e (phi) of observation combined with single differences GF and MWGF)、ε(LMW) Comparing, if: l Δ ΦGF-ε(ΦGF)|>4δ(ΦGF) Or | Δ LMW-ε(LMW)|>4δ(LMW) If so, determining that cycle slip occurs, and marking the carrier phase observation value of the satellite; otherwise, considering that no cycle slip occurs, and not marking;
step S43, calculating a double-difference phase observation value and a high-order difference Δ Φ (t) of the double-difference phase observation value through a high-order difference method according to the carrier phase observation values of all satellites1,t2,t3);
Step S44, filtering and estimating the high-order difference of the double-difference phase observation value obtained in the step S43 by adopting Kalman filtering, weakening the influence of noise, and obtaining a filtering estimation value epsilon ([ delta ] phi) and a variance delta ([ delta ] phi);
comparing the high-order difference of the double-difference phase observations of the current epoch to the filtered estimate ε (Δ Φ) and the variance δ (Δ Φ), if: |. Δ Φ (t)1,t2,t3)-ε(▽ΔΦ)|>4 δ ([ Δ ] Φ), then a cycle slip is considered to have occurred, and the carrier phase observation for that satellite is marked; otherwise, the cycle slip is not considered to occur, and the marking is not carried out.
Further technical solution is that, in step S31, a calculation formula of the double-difference pseudorange observed value is as follows:
Figure BDA0002501805820000041
in the formula: Δ represents a double difference operator;
Figure BDA0002501805820000042
representing double-difference pseudorange observations;
Figure BDA0002501805820000043
representing double-difference satellite geometric distance;
Figure BDA0002501805820000044
representing pseudo-range double-difference observed value noise;
residual error v before test0The calculation formula of (a) is as follows:
Figure BDA0002501805820000045
in the formula: v. of0Representing the pre-trial residual; b denotes a coefficient matrix.
The further technical proposal is that the posterior residual error in the step S33
Figure BDA0002501805820000046
The calculation formula is as follows:
Figure BDA0002501805820000047
in the formula:
Figure BDA0002501805820000048
representing a reference estimate;
the RMS value δ of the post-test residuals is calculated as follows:
Figure BDA0002501805820000049
in the formula: delta denotes the post-test residualN represents the satellite number of the current epoch, table viShows the posterior residual, w, of the ith satelliteiAnd the pseudo range observation value weight of the ith satellite.
Further technical solution is that, in step S41, a calculation formula of the single-difference GF combination observation value is as follows:
Figure BDA00025018058200000410
in the formula: delta N1,m,k,ΔN2,m,kRespectively represent L1、L2The carrier phase observations the single-difference ambiguity between stations,
Figure BDA00025018058200000411
noise of an observation value is combined by single difference GF between stations; lambda [ alpha ]1、λ2Respectively represent L1、L2Wavelength of carrier phase observation.
The further technical solution is that the calculation formula of the single difference MW combined observed value in step S41 is as follows:
Figure BDA00025018058200000412
in the formula: Δ LMW,m,kRepresents the combined observed value of single difference MW between stations, Delta Nwid,m,kExpressing single-difference wide lane ambiguity between stations, Delta sigmam,kWhich represents the delay of the hardware, is,
Figure BDA0002501805820000051
and representing the noise of the single difference MW combined observed value between stations.
Further technical solution is that, in step S43, a calculation formula of the double-difference phase observed value is as follows:
Figure BDA0002501805820000059
in the formula: i. j represents the satellite number, k, m represent the survey station number, and λ represents the carrier phaseThe wavelength of the bits is such that,
Figure BDA0002501805820000052
representing the geometric distance of the double-differenced satellites,
Figure BDA0002501805820000053
representing a double-difference integer ambiguity,
Figure BDA0002501805820000054
representing double difference carrier phase observation noise.
Further technical solution is that, in step S43, a calculation formula of a high-order difference of the double-difference phase observed values is as follows:
Figure BDA0002501805820000055
in the formula: t is t1、t2、t3Representing the three epoch time instants that are adjacent,
Figure BDA0002501805820000056
a second order difference representing a double difference carrier phase observation,
Figure BDA0002501805820000057
represents the second difference of the geometric distance of the double-difference satellites,
Figure BDA0002501805820000058
the observation is noisy.
Big dipper GNSS satellite real time kinematic fixes a position directional data preprocessing system, includes:
the Beidou receiver I and the Beidou receiver II are used for receiving data of a satellite;
the data forwarder is electrically connected with the Beidou receiver I and the Beidou receiver II respectively and used for forwarding the received satellite data to the data preprocessing server;
the data preprocessing server is electrically connected with the data transponder and specifically comprises a positioning residual error processing module which is used for judging a pseudo-range observed value of a satellite so as to delete the pseudo-range observed value of the satellite with gross error;
the single difference GF processing module is used for marking whether cycle slip occurs to the carrier phase observed value of the satellite;
the single difference MW processing module is used for marking whether cycle slip occurs to the carrier phase observation value of the satellite;
the high-order difference processing module is used for marking whether cycle slip occurs on the carrier phase observation value of the satellite;
and the result storage server is electrically connected with the data preprocessing service and is used for storing the data after the satellite data processing.
The invention has the following beneficial effects: according to the method, the spatial correlation of the ultra-short baseline Beidou/GNSS errors is fully utilized, on one hand, a GF and MW combined observation value which is irrelevant to the geometric distance is constructed by utilizing the traditional combined observation value, the influence of ionospheric delay and multipath effect on cycle slip detection is well eliminated through single difference between stations, on the other hand, a high-order difference observation value which slowly and smoothly changes along with time is constructed, and the high-order difference observation value is complementary with a GF and MW combined observation value method, so that the reliability and the effectiveness of cycle slip detection in a complex environment can be improved. Meanwhile, aiming at the characteristics of short and fixed baseline of an automatic driving application scene, the invention provides a method for eliminating obviously larger gross errors by using the pre-test residual error, and improves the stability of least square parameter estimation and the reliability of post-test residual error distribution, thereby improving the effectiveness and reliability of gross error detection in a complex environment.
Drawings
FIG. 1 is a block flow diagram of the present invention;
fig. 2 is a block diagram of the architecture of the present invention.
Detailed Description
The present invention will be further described with reference to the following examples and the accompanying drawings.
As shown in fig. 2, the Beidou GNSS satellite real-time dynamic positioning and orientation data preprocessing system includes:
the Beidou receiver I and the Beidou receiver II are used for receiving data of a satellite;
the data forwarder is electrically connected with the Beidou receiver I and the Beidou receiver II respectively and used for forwarding the received satellite data to the data preprocessing server;
the data preprocessing server is electrically connected with the data transponder and specifically comprises a positioning residual error processing module which is used for judging a pseudo-range observed value of a satellite so as to delete the pseudo-range observed value of the satellite with gross error;
the single difference GF processing module is used for marking whether cycle slip occurs to the carrier phase observed value of the satellite;
the single difference MW processing module is used for marking whether cycle slip occurs to the carrier phase observation value of the satellite;
the high-order difference processing module is used for marking whether cycle slip occurs on the carrier phase observation value of the satellite;
and the result storage server is electrically connected with the data preprocessing service and is used for storing the data after the satellite data processing.
As shown in fig. 1, the method for preprocessing the Beidou GNSS satellite real-time positioning orientation data of the present invention comprises the following steps:
step 1, inputting time K, original non-differential observation data of all satellites of Beidou/GNSS of the observation station 1 and the observation station 2 and broadcast ephemeris.
And 2, respectively calculating all satellite coordinates and clock error correction numbers by using a satellite position calculation standard algorithm, specifically referring to section 2.3 of GPS measurement and data processing of Wuhan university Press, deleting satellite data without a broadcast ephemeris product, and deleting data of unhealthy satellites according to satellite health marks provided by the broadcast ephemeris.
Step 3, calculating non-difference GF and MW combined observed values of all satellites of the observation station 1 and the observation station 2 according to the formulas (1) and (2), calculating single-difference GF and MW combined observed values between the stations according to the formulas (3) and (4) respectively; meanwhile, pseudo-range and phase double-difference observed values are calculated according to equations (5) and (6), respectively.
Figure BDA0002501805820000071
Wherein phi1,Φ2Respectively represent L1、L2Carrier phase observation, λ1,λ2Respectively represent L1、L2Wavelength of carrier phase observation, f1,f2Respectively represent L1、L2Frequency of carrier phase observations, N1,N2Respectively represent L1、L2Ambiguity of carrier phase observation, I denotes L1Ionospheric delays on the carrier-phase observations,
Figure BDA0002501805820000081
noise representing GF combined observations.
Figure BDA0002501805820000082
Wherein λwidRepresenting a wide-lane observation wavelength, NwidRepresenting widelane ambiguity, sigma representing hardware delay,
Figure BDA0002501805820000083
representing non-difference MW combined observation noise.
Figure BDA0002501805820000084
Wherein Δ N1,m,k,ΔN2,m,kRespectively represent L1、L2The carrier phase observations the single-difference ambiguity between stations,
Figure BDA0002501805820000085
the single difference GF between stations combines the noise of the observed values.
Figure BDA0002501805820000086
Wherein, Δ LMW,m,kRepresents the combined observed value of single difference MW between stations, Delta Nwid,m,kExpressing single-difference wide lane ambiguity between stations, Delta sigmam,kWhich represents the delay of the hardware, is,
Figure BDA0002501805820000087
and representing the noise of the single difference MW combined observed value between stations.
Figure BDA0002501805820000088
Where Δ represents a double difference operator,
Figure BDA0002501805820000089
representing a double-differenced pseudorange observation,
Figure BDA00025018058200000810
representing the geometric distance of the double-differenced satellites,
Figure BDA00025018058200000811
representing pseudorange double-difference observation noise.
Figure BDA00025018058200000812
Where i, j denotes the satellite number, k, m denotes the station number, λ denotes the carrier phase wavelength,
Figure BDA00025018058200000813
representing the geometric distance of the double-differenced satellites,
Figure BDA00025018058200000814
representing a double-difference integer ambiguity,
Figure BDA00025018058200000815
representing double difference carrier phase observation noise.
Step 4, adopting static successive filtering to respectively carry out filtering estimation on the single-difference GF and MW combined observed values obtained in the step 3, and weakening noiseInfluence, resulting in an estimate ε (Φ) of the single-difference GF and MW combined observationsGF)、ε(LMW) And estimating the variance δ (Φ)GF)、δ(LMW)。
Comparing the single difference GF and MW combined observed value of the current epoch with the static successive filtering result, if: l Δ ΦGF-ε(ΦGF)|>4δ(ΦGF) Or | Δ LMW-ε(LMW)|>4δ(LMW) If not, the cycle slip is not considered to occur. And if the cycle slip of the i satellite in the time K is judged, restarting the static successive filtering process of the i satellite in the step at the time K + 1.
Step 5, simultaneously with the step 4, calculating the high-order difference of the double-difference phase observed value according to the formula (7), and performing filtering estimation on the high-order difference of the double-difference phase observed value by adopting Kalman filtering to weaken the influence of noise to obtain a filtering estimation value epsilon ([ delta ] phi) and a variance delta ([ delta ] phi),
Figure BDA0002501805820000091
wherein t is1,t2,t3Representing the three epoch time instants that are adjacent,
Figure BDA0002501805820000092
a second order difference representing a double difference carrier phase observation,
Figure BDA0002501805820000093
represents the second difference of the geometric distance of the double-difference satellites,
Figure BDA0002501805820000094
the observation is noisy.
Comparing the high-order difference of the double-difference phase observed value of the current epoch with the filtering result, if:
|▽ΔΦ(t1,t2,t3)-ε(▽ΔΦ)|>4 δ ([ Δ ] Φ), then cycle slip is considered to have occurred, otherwise cycle slip is considered to have not occurred. If the cycle slip of the i satellite in the time K is judged, at the time K +1,the kalman filtering process restarts in this step for satellite i.
Step 6, simultaneously with the step 4, selecting an initial parameter value X0The pseudo-range double-difference observation equation is linearized (1, 1, 1), and the pre-test residual V of the double-difference pseudo-range observation value is calculated according to equation (8)0If | V0|>And 30, rejecting the observation value of the satellite if the gross error exists.
V0=BX0-L (8)
Wherein V0Representing the pre-trial residuals, L representing the double-differenced pseudorange observations, and B representing the coefficient matrix.
The error equation (8) without gross error is listed in matrix form as follows:
Figure BDA0002501805820000095
obtaining parameter estimates using least squares estimation
Figure BDA0002501805820000096
Respectively calculating the post-test residual v and the RMS of the post-test residual of each satellite according to the formulas (9) and (10), if | v>4 RMS, the gross error is considered to exist, and the observation of the satellite is rejected.
Figure BDA0002501805820000101
Figure BDA0002501805820000102
Where δ represents the RMS value of the post-test residual, n represents the satellite number of the current epoch, viRepresenting the post-test residual of the ith satellite, wiPseudo range observation value weight of ith satellite
And 7, inputting the observation data at the moment K +1, and re-executing the steps 1-6 until the observation data at all the moments are preprocessed.
And 8, finally obtaining a pseudo-range observation value without gross error, a carrier phase observation value marked as cycle slip occurrence, and a carrier phase observation value marked as cycle slip non-occurrence.
Although the present invention has been described with reference to the above embodiments, it should be understood that the present invention is not limited to the above embodiments, and those skilled in the art can make various changes and modifications without departing from the scope of the present invention.

Claims (10)

1. The Beidou GNSS satellite real-time positioning and orientation data preprocessing method is characterized by comprising the following steps:
step S1, inputting carrier phase observation values, pseudo-range observation values and broadcast ephemeris of a plurality of Beidou/GNSS satellites of the observation station 1 and the observation station 2;
step S2, calculating the coordinates and the correction of clock error of all the satellites input in the step S1 through a satellite position calculation standard algorithm according to the broadcast ephemeris, deleting the satellite data without the broadcast ephemeris product, and deleting the data of unhealthy satellites according to the satellite health identification provided by the broadcast ephemeris;
step S3, judging pseudo-range observation values of all satellites input in step S1 by using a pseudo-range positioning pre-test residual error judgment method and a pseudo-range positioning pre-test residual error judgment method to delete the pseudo-range observation values of the satellites with gross errors;
step S4, marking whether cycle slip occurs on the carrier phase observation values of all the satellites input in the step S1 by using a single difference GF and MW combined observation value method and a high-order difference value method;
and step S5, finally obtaining a pseudo-range observation value without gross error, a carrier phase observation value marked as cycle slip occurrence, and a carrier phase observation value marked as cycle slip non-occurrence.
2. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data according to claim 1, wherein the specific process of the step S3 is as follows:
step S31, calculating double-difference pseudo-range observed values of all satellites by a pseudo-range positioning residual error analysis method according to coordinates and clock error correction numbers of all satellites, and setting an initial parameter value X0(1, 1, 1), linearizing pseudo-range double-difference observation equation, and calculating to obtain pre-test residual v of double-difference pseudo-range observation value0(ii) a If | v0|>30, if the gross error exists, rejecting a pseudo range observation value of the satellite; | v0If the absolute value is less than or equal to 30, considering that no gross error exists, and carrying out the next step;
step S32, obtaining parameter estimation value by least square estimation
Figure FDA0002501805810000011
The baseline distance constraints are as follows:
const=sqrt(N2+E2+U2)
wherein N, E, U respectively represent parameter vectors
Figure FDA0002501805810000021
North, middle, east, and high, const representing the baseline length;
step S5, and estimating the parameters according to step 32
Figure FDA0002501805810000022
Calculating the post-test residual v and RMS (root mean square) values delta of the post-test residual of all satellites; if | v $>And if the absolute value is less than or equal to 4 multiplied by delta, the pseudorange observed value of the satellite is considered to be not present, and the pseudorange observed value of the satellite is reserved.
3. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data according to claim 1, wherein the specific process of the step S4 is as follows:
step S41, calculating single difference GF and MW combined observation value delta phi through an inter-station single difference GF and MW combined observation value method according to carrier phase observation values of all satellitesGF、ΔLMW
Step S42, respectively carrying out filtering estimation on the single-difference GF and MW combined observed values obtained in step 41 by adopting static successive filtering, weakening the influence of noise, and obtaining the estimated value epsilon (phi) of the single-difference GF and MW combined observed valuesGF)、ε(LMW) And estimating the variance δ (Φ)GF)、δ(LMW);
Combining single difference GF and MW of current epoch into observed value delta phiGF、ΔLMWEstimate e (phi) of observation combined with single differences GF and MWGF)、ε(LMW) Comparing, if: l Δ ΦGF-ε(ΦGF)|>4δ(ΦGF) Or | Δ LMW-ε(LMW)|>4δ(LMW) If so, determining that cycle slip occurs, and marking the carrier phase observation value of the satellite; otherwise, considering that no cycle slip occurs, and not marking;
step S43, calculating the double-difference phase observed value and the high-order difference of the double-difference phase observed value by a high-order difference method according to the carrier phase observed values of all satellites
Figure FDA0002501805810000023
Step S44, filtering and estimating the high-order difference of the double-difference phase observed value obtained in the step S5 by Kalman filtering, weakening the influence of noise and obtaining a filtering estimated value
Figure FDA0002501805810000024
Sum variance
Figure FDA0002501805810000025
The high-order difference of the double-difference phase observed value of the current epoch and the filtering estimation value
Figure FDA0002501805810000026
Sum variance
Figure FDA0002501805810000027
Make a comparison if
Figure FDA0002501805810000028
Considering that cycle slip occurs, and marking the carrier phase observed value of the satellite; otherwise, the cycle slip is not considered to occur, and the marking is not carried out.
4. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data of claim 2, wherein the double-difference pseudorange observation in step S31 is calculated as follows:
Figure FDA0002501805810000031
in the formula:
Figure FDA0002501805810000032
representing a double difference operator;
Figure FDA0002501805810000033
representing double-difference pseudorange observations;
Figure FDA0002501805810000034
representing double-difference satellite geometric distance;
Figure FDA0002501805810000035
representing pseudo-range double-difference observed value noise;
residual error v before test0The calculation formula of (a) is as follows:
Figure FDA0002501805810000036
in the formula: v. of0Representing the pre-trial residual; b denotes a coefficient matrix.
5. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data of claim 2, wherein the posterior residual error in the step S33 is
Figure FDA0002501805810000037
The calculation formula is as follows:
Figure FDA0002501805810000038
in the formula:
Figure FDA0002501805810000039
representing a reference estimate;
the RMS value δ of the post-test residuals is calculated as follows:
Figure FDA00025018058100000310
in the formula: delta represents the RMS value of the post-test residual, n represents the satellite number of the current epoch, table viShows the posterior residual, w, of the ith satelliteiAnd the pseudo range observation value weight of the ith satellite.
6. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data according to claim 3, wherein the single-difference GF combined observation value in the step S41 is calculated according to the following formula:
Figure FDA00025018058100000311
in the formula: delta N1,m,k,ΔN2,m,kRespectively represent L1、L2The carrier phase observations the single-difference ambiguity between stations,
Figure FDA00025018058100000312
GF combined observation of single difference between stationsNoise of the value; lambda [ alpha ]1、λ2Respectively represent L1、L2Wavelength of carrier phase observation.
7. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data of claim 3, wherein the calculation formula of the single difference MW combined observation value in the step S41 is as follows:
Figure FDA0002501805810000041
in the formula: Δ LMW,m,kRepresents the combined observed value of single difference MW between stations, Delta Nwid,m,kExpressing single-difference wide lane ambiguity between stations, Delta sigmam,kWhich represents the delay of the hardware, is,
Figure FDA0002501805810000042
and representing the noise of the single difference MW combined observed value between stations.
8. The method for preprocessing the Beidou GNSS satellite real-time positioning and orientation data of claim 1, wherein the calculation formula of the double-difference phase observation in the step S43 is as follows:
Figure FDA0002501805810000043
in the formula: i. j denotes the satellite number, k, m denote the station number, λ denotes the carrier phase wavelength,
Figure FDA0002501805810000044
representing the geometric distance of the double-differenced satellites,
Figure FDA0002501805810000045
representing a double-difference integer ambiguity,
Figure FDA0002501805810000046
representing double difference carrier phase observation noise.
9. The method as claimed in claim 6, wherein the calculation formula of the high order difference of the double difference phase observations in the step S43 is as follows:
Figure FDA0002501805810000047
in the formula: t is t1、t2、t3Representing the three epoch time instants that are adjacent,
Figure FDA0002501805810000048
a second order difference representing a double difference carrier phase observation,
Figure FDA0002501805810000049
represents the second difference of the geometric distance of the double-difference satellites,
Figure FDA00025018058100000410
the observation is noisy.
10. Big dipper GNSS satellite real time kinematic fixes a position directional data preprocessing system, its characterized in that includes:
the Beidou receiver I and the Beidou receiver II are used for receiving data of a satellite;
the data forwarder is electrically connected with the Beidou receiver I and the Beidou receiver II respectively and used for forwarding the received satellite data to the data preprocessing server;
the data preprocessing server is electrically connected with the data transponder and specifically comprises a positioning residual error processing module which is used for judging a pseudo-range observed value of a satellite so as to delete the pseudo-range observed value of the satellite with gross error;
the single difference GF processing module is used for marking whether cycle slip occurs to the carrier phase observed value of the satellite;
the single difference MW processing module is used for marking whether cycle slip occurs to the carrier phase observation value of the satellite;
the high-order difference processing module is used for marking whether cycle slip occurs on the carrier phase observation value of the satellite;
and the result storage server is electrically connected with the data preprocessing service and is used for storing the data after the satellite data processing.
CN202010434736.5A 2020-05-21 2020-05-21 Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method Pending CN112305574A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010434736.5A CN112305574A (en) 2020-05-21 2020-05-21 Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010434736.5A CN112305574A (en) 2020-05-21 2020-05-21 Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method

Publications (1)

Publication Number Publication Date
CN112305574A true CN112305574A (en) 2021-02-02

Family

ID=74336589

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010434736.5A Pending CN112305574A (en) 2020-05-21 2020-05-21 Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method

Country Status (1)

Country Link
CN (1) CN112305574A (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113009517A (en) * 2021-03-02 2021-06-22 中国铁路设计集团有限公司 Beidou multi-antenna array-based high-speed railway infrastructure deformation monitoring method
CN113281796A (en) * 2021-07-23 2021-08-20 腾讯科技(深圳)有限公司 Position determining method, speed determining method, device, equipment and storage medium
CN113865592A (en) * 2021-09-09 2021-12-31 河海大学 Multi-path parameterization method and storage medium suitable for multi-frequency GNSS precision navigation positioning
CN113970773A (en) * 2021-10-29 2022-01-25 北京百度网讯科技有限公司 Positioning method and device and electronic equipment
CN114167456A (en) * 2021-10-21 2022-03-11 国网新源控股有限公司 Method, system and terminal for detecting and repairing clock jumps of receiver
CN115201864A (en) * 2022-07-13 2022-10-18 涟漪位置(广州)科技有限公司 Method, device, storage medium and equipment for detecting clock difference jump of satellite
CN115877414A (en) * 2023-02-20 2023-03-31 广州导远电子科技有限公司 Fixed solution coordinate checking method and device, storage medium and electronic equipment
CN115993623A (en) * 2023-03-24 2023-04-21 武汉大学 Adaptive star selection method, device, equipment and readable storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581452A (en) * 2018-12-18 2019-04-05 辽宁工程技术大学 A kind of GNSS reference station ambiguity of carrier phase calculation method
CN109655854A (en) * 2019-02-21 2019-04-19 哈尔滨工程大学 It is a kind of based on zero base line constraint multi-receiver PPP quickly restrain technology again
CN110907960A (en) * 2018-09-17 2020-03-24 千寻位置网络有限公司 Cycle slip detection method and device based on K-Means dynamic clustering analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110907960A (en) * 2018-09-17 2020-03-24 千寻位置网络有限公司 Cycle slip detection method and device based on K-Means dynamic clustering analysis
CN109581452A (en) * 2018-12-18 2019-04-05 辽宁工程技术大学 A kind of GNSS reference station ambiguity of carrier phase calculation method
CN109655854A (en) * 2019-02-21 2019-04-19 哈尔滨工程大学 It is a kind of based on zero base line constraint multi-receiver PPP quickly restrain technology again

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
YANG YUANXI 等: "Preliminary assessment of the navigation and positioning performance of BeiDou regional navigation satellite system", 《SCI CHINA EARTH SCI》 *
曾琪: "BDS/GPS组合动态相对定位及其质量控制方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
蔡华等: "GNSS实时数据质量控制", 《武汉大学学报(信息科学版)》 *
谢瑞煜 等: "一种基于双差观测的BDS周跳探测与修复方", 《火力与指挥控制》 *
郭向等: "基于单频星载GPS数据的低轨卫星精密定轨", 《中国空间科学技术》 *
陈向东 等: "双源车载差分卫星导航定位稳定性分析", 《导航定位学报》 *
陈玉龙: "GPS/BDS基线解算方法研究及程序实现", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113009517A (en) * 2021-03-02 2021-06-22 中国铁路设计集团有限公司 Beidou multi-antenna array-based high-speed railway infrastructure deformation monitoring method
CN113281796A (en) * 2021-07-23 2021-08-20 腾讯科技(深圳)有限公司 Position determining method, speed determining method, device, equipment and storage medium
CN113281796B (en) * 2021-07-23 2021-10-15 腾讯科技(深圳)有限公司 Position determining method, speed determining method, device, equipment and storage medium
CN113865592A (en) * 2021-09-09 2021-12-31 河海大学 Multi-path parameterization method and storage medium suitable for multi-frequency GNSS precision navigation positioning
CN113865592B (en) * 2021-09-09 2024-05-10 河海大学 Multipath parameterization method and storage medium suitable for multi-frequency GNSS precise navigation positioning
CN114167456A (en) * 2021-10-21 2022-03-11 国网新源控股有限公司 Method, system and terminal for detecting and repairing clock jumps of receiver
CN113970773A (en) * 2021-10-29 2022-01-25 北京百度网讯科技有限公司 Positioning method and device and electronic equipment
CN113970773B (en) * 2021-10-29 2024-04-16 北京百度网讯科技有限公司 Positioning method and device and electronic equipment
CN115201864A (en) * 2022-07-13 2022-10-18 涟漪位置(广州)科技有限公司 Method, device, storage medium and equipment for detecting clock difference jump of satellite
CN115877414A (en) * 2023-02-20 2023-03-31 广州导远电子科技有限公司 Fixed solution coordinate checking method and device, storage medium and electronic equipment
CN115877414B (en) * 2023-02-20 2023-04-28 广州导远电子科技有限公司 Fixed solution coordinate checking method and device, storage medium and electronic equipment
CN115993623A (en) * 2023-03-24 2023-04-21 武汉大学 Adaptive star selection method, device, equipment and readable storage medium

Similar Documents

Publication Publication Date Title
CN112305574A (en) Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method
CN108415049B (en) Method for improving network RTK double-difference wide lane ambiguity fixing accuracy
CN108802782B (en) Inertial navigation assisted Beidou three-frequency carrier phase integer ambiguity solving method
CN111751853B (en) GNSS dual-frequency carrier phase integer ambiguity resolution method
CN111965673B (en) Time-frequency transfer method of single-frequency precise single-point positioning algorithm based on multiple GNSS
CN104181562B (en) The satellite of a kind of GLONASS is preferably and localization method
CN110045407A (en) A kind of distribution pseudo satellite, pseudolite/GNSS optimum position method
CN103529462A (en) Probing and repairing method for dynamic cycle slip of global navigation satellite system
CN104459722B (en) A kind of integer ambiguity certificate authenticity method based on redundant obser ration part
CN112285745B (en) Three-frequency ambiguity fixing method and system based on Beidou third satellite navigation system
CN109613582B (en) Vehicle-mounted real-time single-frequency meter-level pseudo-range positioning method
CN113358017B (en) Multi-station cooperative processing GNSS high-precision deformation monitoring method
CN111983641A (en) Method for generating Beidou satellite-based augmentation system integrity parameters in real time
CN107966722A (en) A kind of GNSS satellite clock solutions method
CN105510945A (en) PPP positioning method applied to satellite navigation landing outfield detection
CN114895330A (en) Single-station displacement monitoring method, equipment and storage medium based on broadcast ephemeris
CN114935770B (en) Method and device for accelerating precision single-point positioning convergence speed by multiple calendars
CN116009042A (en) Method and system for detecting relative deformation in real time by difference between single-station carrier epochs
CN115220078A (en) GNSS high-precision positioning method and navigation method based on carrier phase difference
CN109143289B (en) GNSS single-station displacement monitoring method
CN113819863B (en) Deformation monitoring method and system
Zheng et al. Multipath mitigation for improving GPS narrow-lane uncalibrated phase delay estimation and speeding up PPP ambiguity resolution
Chen et al. A new cycle slip detection and repair method for single-frequency GNSS data
Bisnath Relative Positioning and Real‐Time Kinematic (RTK)
CN110058274B (en) Method and system for monitoring time difference between satellite navigation systems

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210202