CN108919321B - A kind of GNSS positioning Detection of Gross Errors method based on trial and error method - Google Patents

A kind of GNSS positioning Detection of Gross Errors method based on trial and error method Download PDF

Info

Publication number
CN108919321B
CN108919321B CN201810481017.1A CN201810481017A CN108919321B CN 108919321 B CN108919321 B CN 108919321B CN 201810481017 A CN201810481017 A CN 201810481017A CN 108919321 B CN108919321 B CN 108919321B
Authority
CN
China
Prior art keywords
observation
error
gnss
rough
errors
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.)
Active
Application number
CN201810481017.1A
Other languages
Chinese (zh)
Other versions
CN108919321A (en
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.)
Changan University
Original Assignee
Changan University
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 Changan University filed Critical Changan University
Priority to CN201810481017.1A priority Critical patent/CN108919321B/en
Publication of CN108919321A publication Critical patent/CN108919321A/en
Application granted granted Critical
Publication of CN108919321B publication Critical patent/CN108919321B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

Abstract

The invention discloses a kind of, and the GNSS based on trial and error method positions Detection of Gross Errors method, realizes process are as follows: construction GNSS observational equation carries out parameter Estimation using least square (LS), to obtain observation residual vector V;Error estimator in unit of account powerBeing adaptively calculated error in weight unit using excess observation number f transfinites threshold value, judges GNSS observation with the presence or absence of rough error;Rough error if it exists then utilizes the thought based on trial and error method to realize accurate Gross postionning;Rejecting is detected containing rough error observation, re-starts LS parameter Estimation.This method is not necessarily to carry out observation one by one residual error hypothesis testing, therefore can cause to reject part GPS/BDS observation excessive or very fewly because threshold value setting is unreasonable in avoiding the problem that conventional residual inspection to a certain degree.This method is simple, it is easy to accomplish, application effect is relatively good, and GNSS especially more can be significantly improved in the case where observing environment is poor and is averaged positioning accuracy.

Description

A kind of GNSS positioning Detection of Gross Errors method based on trial and error method
Technical field
The invention belongs to GNSS high accuracy positioning fields, are related to a kind of GNSS positioning Detection of Gross Errors knowledge based on trial and error method Other technology, it is actually to answer with fields such as Global Satellite Navigation System (GNSS), Surveying of Precise Control, mobile vehicle location navigations With background, it can be used for the mobile vehicle high accuracy positioning application side under the conditions of GNSS observation rough error frequently occurs in complex environment To.
Background technique
Global Satellite Navigation System (GNSS) has many advantages, such as round-the-clock, high-precision, have application to geodesic survey, navigation, The various fields such as Geological Hazards Monitoring.But GNSS signal is more fragile simultaneously, and especially under complex environment, GNSS observation is difficult Exempt to be polluted by rough error, if rough error is not dealt carefully with, will lead to final positioning result and deviate reality, give practical application band It is larger unfavorable to come, therefore Detection of Gross Errors is the important content of GNSS location data processing Quality Control Links.
At present in GNSS data processing, rough error processing can be summarized as both of which, i.e., rough error is included into stochastic model Robust filter mode and the average drifting mode that rough error is included into function model.It is external in terms of robust mode treatment rough error Mainly there are robusts scheme, the domestic scholars such as Huber weight-function method, Tukey weight-function method and Denmark's weight-function method to construct robust The robusts schemes such as IGG scheme, double factor equivalence weight and equivalent variance-covariance function.These schemes have obtained more application, Generally to power battle array deal with, by containing rough error observation weigh battle array carry out drop power, the influence of the observation containing rough error can be weakened, But it is likely to encounter the situation of matrix rank defect.The method that rough error is eliminated based on mean-shifted models during processing will not be to power Battle array makes change, will not correspondingly encounter the situation of matrix rank defect.It is usually after testing observation based on mean-shifted models Residual sequence carries out hypothesis testing, such as common normal distribution-test, χ2Distribution inspection, t distribution inspection, by underproof observation Value is rejected.Traditional residual sequence hypothesis testing has determined the threshold value whether residual error transfinites given confidence level α, In view of GNSS carrier observations residual error sensibility and may be subjected to the influence of model error, threshold size visits rough error It is extremely important to survey influential effect, it is bigger than normal may cause cannot accurately by elimination of rough difference, it is less than normal may result in some quality compared with Good observation is also removed.Real data processing in especially complex environment observation the quality of data it is poor in the case where, very Difficult matter first gives a very reasonable threshold value, therefore this Detection of Gross Errors method effect is frequently not highly desirable.
Summary of the invention
It is an object of the invention to overcome the deficiencies in the prior art, propose a kind of without carrying out residual error to observation one by one The Detection of Gross Errors method of hypothesis testing.This method passes through error estimator in weight unit first and judges whether there is rough error, then sharp Gross postionning, referred to as two-step method Detection of Gross Errors are carried out with trial and error method thought.This method is avoiding conventional residual to a certain degree Lead to the problem of rejecting part GNSS observation excessive or very fewly because threshold value setting is unreasonable in inspection, is under complex environment GNSS Detection of Gross Errors provides a kind of effective processing method.
Technical scheme is as follows:
A kind of GNSS positioning Detection of Gross Errors method based on trial and error method, includes the following steps
(1) GNSS observational equation is constructed, parameter Estimation is carried out using least square, to obtain observation residual vector V;
(2) error estimator in calculating observation value weight unit
(3) it is adaptively calculated using excess observation number fTransfinite threshold value, and according toWhether transfinite, judges GNSS observation whether there is rough error, if it does not exist rough error, then directly exports positioning result;
(4) if there are rough errors for observation, accurate Gross postionning, specific steps are carried out based on trial and error method are as follows:
(4.1) rough error judges: under the conditions of only rejecting an observation every time, obtaining the n unit rejected after observation Standard deviation is weighed, then there are rough errors for the corresponding observation of minimum therein;
(4.2) threshold decision: carrying out the threshold decision that transfinites to the minimum, if being more than the threshold value of default, rejecting should The corresponding observation of minimum, obtains n-1 observation;
(4.3) traversal is rejected: repetitive cycling carries out step (4.1)~(4.2), until the minimum of acquisition is set less than system Fixed threshold value, and reject observation corresponding to the minimum finally obtained.
(5) least-squares parameter estimation is re-started to remaining observation, obtains the GNSS high-precision influenced without rough error Positioning result.
Further, GNSS observation L described in step (1) is mainly made of pseudorange and carrier observations, corresponding residual Difference vector V is solved according to the following formula:
V=Bx-L
Wherein x is to include carrier three-dimensional coordinate, the parameter estimation of ambiguity of carrier phase;B is the corresponding design square of x Battle array, L number of observation are n;
Further, error estimator in weight unit described in step (2)Are as follows:
Wherein f is excess observation number, and P is the corresponding weight matrix of observation.
Further, the decision criteria that rough error whether there is in step (3) are as follows: if meetingThen determine Rough error is not present in GNSS observation, and otherwise there are rough errors for some or multiple observations;Wherein σ0Indicate GNSS observation noise, Fact is editor's factor,
F is excess observation number.
Further, in step (4.1) rough error the judgment method is as follows: if observation there are rough errors, successively only reject the i-th (i =1,2 ..., n) the corresponding observation L of satellitei, and unit of account weighs standard deviation againAccording to formula
It determines and corresponds to observation by kth satellite there are rough errors.
Further, the method for threshold decision is in step (4.2): determining the corresponding observation L of kth satellitekIt deposits After rough error, if corresponding weight unit standard deviation is still unsatisfactory forIllustrate that there are still thick in remaining observation Difference;
Further, the threshold value set in the rejecting of step (4.3) traversal is Fact × σ0
Compared with the prior art, the present invention has the following advantages:
The present invention first tests to GNSS Observation value error, determines whether that there are rough errors, then utilizes trial if it exists Method realizes Gross postionning, without carrying out residual error hypothesis testing to observation one by one, therefore can avoid to a certain degree conventional residual Lead to the problem of rejecting part GNSS observation excessive or very fewly because threshold value setting is unreasonable in inspection.This method is simple, easily In realization, application effect is relatively good, especially in the case where observing environment is poor, is guaranteeing the holding of data effective rate of utilization Under the premise of higher, this method more can significantly improve GNSS and be averaged positioning accuracy.
Detailed description of the invention
Fig. 1 is flow chart of the present invention;
Fig. 2 is static test positioning accuracy (X-direction) statistical chart of the present invention;
Fig. 3 is vehicle-mounted dynamic test positioning accuracy (X-direction) statistical chart of the present invention;
Fig. 4 is the vehicle-mounted dynamic test GPS/BDS plane positioning trajectory diagram of the present invention;
Specific embodiment
The present invention will be further described with reference to the accompanying drawing.
Referring to Fig.1, specific implementation step of the invention is as follows:
Step 1. constructs GNSS observational equation, parameter Estimation is carried out using least square (LS), so that it is residual to obtain observation Difference vector V.
GNSS observation mainly includes pseudorange and carrier wave, location observation equation are as follows:
Above formula,Indicate receiver-satellite s-r geometric distance,WithIt respectively indicates tropospheric delay and ionosphere is prolonged Chi Liang,WithPseudorange and the corresponding receiver-satellite clock correction of carrier wave are respectively indicated, c indicates the light velocity,WithIt respectively indicates Pseudorange and carrier observations noise,Indicate that carrier phase observation data integer ambiguity, λ indicate corresponding wavelength.
Above formula includes more error term, generallys use two receivers and two satellites constitute classical double difference observation side Journey can slacken troposphere and ionospheric error in various degree and eliminate satellite-receiver clock-offsets.Double difference observation equation is such as Under:
In above formula,Other are similar.Under normal circumstances when baseline is shorter, double difference in above formula Ionosphere and tropospheric error can be ignored.
Observation residual vector V is solved according to the following formula:
V=Bx-L
Wherein x is parameter estimation, to include carrier three-dimensional coordinate, ambiguity of carrier phase etc.;B is the corresponding design of x Matrix is determined according to different station-keeping modes, and L number of observation is denoted as n;
Error estimator in step 2. calculating observation value weight unit
Error in weight unitAre as follows:
Wherein f is excess observation number, and P is the corresponding weight matrix of observation.
Step 3. is adaptively calculated using excess observation number fTransfinite threshold value, according toWhether transfinite judgement GNSS observation whether there is rough error.
Rough error is with the presence or absence of decision criteria, if meeting following formula, determining GNSS observation, there is no rough errors, otherwise a certain There are rough errors for a or multiple observations.
Above formula, σ0It indicates GNSS observation noise, generally rule of thumb gives, for carrier observations, (whole world is fixed by GPS Position system) it can be set to 0.003, BDS (Beidou satellite system) and can be set 0.004, Pseudo-range Observations σ0It can be set to 1;Fact is editor's factor, is set as 3 under normal circumstances relatively rationally, but when observation satellite is less, even if there be no thick Difference,It may also be able to be more than 3 σ0, cause the satellite observation that reject to be removed, seriously affect positioning accuracy, or even make At can not position.Therefore the present invention gives a kind of relatively reasonable calculating Fact value empirical equation by many experiments, as follows:
When satellites in view number is less, editor's factor values can suitably be increased by above formula, i.e., accordingly relaxed in weight unit Error transfinites threshold value, and when observation is more, Fact takes 3.
Step 4. rough error if it exists then utilizes the thought based on trial and error method to realize accurate Gross postionning.
Trial and error method basic thought of the present invention is to carry out traversal one by one to all observations and delete, if a certain observation exists Rough error, after which deletes, error must have larger reduction in the weight unit that recalculates;
Thought based on trial and error method realizes that accurate Gross postionning method is as follows:
If there are rough errors for GPS/BDS observation, the corresponding observation of i-th (i=1,2..., n) satellite is successively only rejected Li, and unit of account power standard deviation (is denoted as again), after all satellites have traversed, it can determine that kth is defended according to the following formula Star corresponds to observation, and there are rough errors.
Above formula meaning is, so that the most corresponding observation of that satellite of error reduction can consider presence in weight unit Rough error.
Determining the corresponding observation L of kth satellitekThere are after rough error, if corresponding weight unit standard deviation is still unsatisfactory forIllustrate there are still rough error in remaining observation, after the observation for rejecting kth satellite, to remaining n-1 A observation continues to realize accurate Gross postionning using trial and error method thought, and circulating repetition carries out, until corresponding weight unit standard Difference is less than Fact × σ0
Step 5. rejecting is detected containing rough error observation, re-starts LS parameter Estimation.
The corresponding observation of rough error is not involved in LS parameter Estimation, and the observation without rough error then re-starts LS parameter and estimates Meter can obtain the GNSS high accuracy positioning result that no rough error influences.
Effect of the invention can be illustrated by following experiments:
1. experimental situation
Verification test of the present invention is all made of GPS/BDS measured data, respectively static observation experiment 1 and dynamic track test 2.Test 1 entirely survey section observing environment it is preferable, GPS be averaged satellites in view number be 7, GPS+BDS satellites in view number be 15, count According to sampling interval 0.5s, satellite elevation mask is set as 10 degree, and when observation is about 15 minutes;2 observation data of test are collected in There is situations such as blocking compared with multi-obstacle avoidance, signal frequently interrupts in certain urban district, base station and rover station sampling interval are 1s, are defended Star elevation mask is 10 degree, and when observation is about 15 minutes, and rover station and base station distance are within 4km.
2. experimental result
Here convenient for comparison, the method for the present invention is abbreviated as two-step method Detection of Gross Errors.
Experiment 1: since observing environment is preferable, single epoch positioning mode is carried out to this static data and is handled.Rough error Detection aspect have chosen three kinds of different schemes to handle, be respectively: no Detection of Gross Errors, it is traditional based on t distribution residual error assume inspection Test (confidence level α is set as 0.01), two-step method Detection of Gross Errors.
In order to compare influence of three kinds of different Detection of Gross Errors schemes to positioning accuracy, count that three kinds of schemes are corresponding singly to be gone through First positioning accuracy only provides X-direction position error sequence chart, as shown in Figure 2 as space is limited.Although caning be found that static observation The quality of data is integrally preferable, but there are still several rough errors for GPS and BDS observation, are distributed residual test and two-step method rough error using t Detection can be good at eliminating influence of these rough errors to positioning.In addition it can find out GPS+BDS dual system relative to GPS Single system is also promoted due to increasing observation, positioning accuracy.
Table 1 has counted the corresponding positioning result of three kinds of rough error processing schemes, including average positioning accuracy, ambiguity resolution feelings Condition and total time-consuming.Due to rough error negligible amounts, t is distributed residual test and two-step method Detection of Gross Errors effect as can be seen from Table 1 Just the same, relative to no Detection of Gross Errors, positioning accuracy and fuzziness fixed rate are promoted, but promote amplitude very little, and And overall calculation efficiency is also without the reduction big because of processing rough error.
Table 1 tests 1 static observation positioning result statistical form
Experiment 2: since GNSS observing environment is complex, recursive least square is taken to this dynamic Vehicular test data Filtering is handled, and ambiguity resolution takes float-solution mode.No Detection of Gross Errors is still taken in terms of Detection of Gross Errors, is distributed based on t Residual error hypothesis testing (confidence level α is set as 0.01), two-step method Detection of Gross Errors three kinds of schemes are handled.
Fig. 3 gives the corresponding GPS and GPS+BDS system X-direction position error sequence of three kinds of difference rough error processing schemes. As can be seen that position error is integrally larger, especially in (100,200) epoch and (600,700) epoch without rough error processing There is very big deviation in the corresponding period.After being distributed residual error hypothesis testing and two-step method Detection of Gross Errors using t, due to rejecting After some more apparent rough errors, positioning accuracy increases.Compare t distribution residual test and two-step method Detection of Gross Errors pair The position error answered, it is found that many epoch two-step method Detection of Gross Errors positioning accuracies are better than t distribution residual test, with GPS For the 45th epoch of single system, observation satellite is 5, eliminates 2 satellite carrier observations using t distribution residual test, Pseudorange One-Point Location can only be finally carried out, and a satellite carrier observation is only eliminated using two-step method Detection of Gross Errors, finally Positioning result is float-solution;The 387th epoch of GPS+BDS dual system, is not visited using t distribution residual test by several 12, satellite Rough error is measured, and two-step method Detection of Gross Errors then accurately eliminates a satellite observation, positioning accuracy is opposite to increase.
Fig. 4 gives no Detection of Gross Errors and the corresponding GPS+BDS rover station plane fortune of two-step method Detection of Gross Errors two schemes Dynamic trajectory diagram.As can be seen that the corresponding rover station motion profile of two-step method Detection of Gross Errors is relative to no Detection of Gross Errors motion profile It is smooth, the corresponding period in dotted ellipse region especially in figure, this at two corresponding location exist it is apparent high The signal blocks of layer building and trees.
This vehicle-mounted vehicle-mounted trial position of dynamic of 2 quantitative statistics of table is as a result, including three direction average localization errors, having Imitate epoch number and total time-consuming situation.In terms of positioning accuracy, in three kinds of different rough error processing schemes, two-step method rough error is visited Corresponding average positioning accuracy highest is surveyed, followed by t is distributed residual error hypothesis testing.GPS+BDS dual system is relative to GPS monosystem System, the corresponding positioning accuracy of three kinds of schemes have a distinct increment;From the point of view of effective epoch number, the corresponding number of two-step method Detection of Gross Errors According to effective rate of utilization relative to no Detection of Gross Errors without being decreased obviously, and t is distributed residual error hypothesis testing due to excessively eliminating Observation causes part epoch that can not position, and reduces data effective rate of utilization;In terms of total time-consuming, two-step method rough error is visited Survey time-consuming at most, but relatively entire experimental stage (about 900 epoch), there is no by very for two-step method Detection of Gross Errors computational efficiency Big influence is fully able to meet the real-time dynamic positioning that sample rate is 1s.
2 dynamic test positioning result statistical form of table

Claims (3)

1. a kind of GNSS based on trial and error method positions Detection of Gross Errors method, which comprises the steps of:
(1) GNSS observational equation is constructed, parameter Estimation is carried out using least square, to obtain observation residual vector V;
(2) error estimator in calculating observation value weight unit
(3) it is adaptively calculated using excess observation number fTransfinite threshold value, and according toWhether transfinite, judges that GNSS is seen Measured value whether there is rough error, if it does not exist rough error, then directly exports positioning result;
The wherein decision criteria that rough error whether there is are as follows: if meetingThen determining GNSS observation, there is no thick Difference, otherwise there are rough errors for some or multiple observations;Wherein σ0Indicate that GNSS observation noise, Fact are editor's factor,
F is excess observation number;
(4) if there are rough errors for observation, accurate Gross postionning, specific steps are carried out based on trial and error method are as follows:
(4.1) Gross postionning: under the conditions of only rejecting an observation every time, the n unit token rejected after observation is obtained Quasi- poor, then there are rough errors for the corresponding observation of minimum therein;
Wherein Gross postionning comprises the concrete steps that: if there are rough errors for observation, successively only rejecting i-th, (i=1 2 ..., n) is defended The corresponding observation L of stari, and unit of account weighs standard deviation againAccording to formula
Kth satellite can be positioned to correspond to observation there are rough errors;
(4.2) threshold decision: to (4.1) minimumThe threshold decision that transfinites is carried out, it, can be true if being more than the threshold value of default There are still rough errors in remaining fixed observation;
The method of threshold decision is: determining the corresponding observation L of kth satellitekThere are after rough error, if corresponding unit token Quasi- difference is still unsatisfactory forIllustrate that there are still rough errors in remaining observation;
(4.3) traversal is rejected: repetitive cycling carries out step (4.1)~(4.2), until error minimum is small in the weight unit of acquisition In threshold value Fact × σ of default0, and reject observation corresponding to the minimum finally obtained;
(5) least-squares parameter estimation is re-started to remaining observation, obtains the GNSS high accuracy positioning influenced without rough error As a result.
2. the GNSS according to claim 1 based on trial and error method positions Detection of Gross Errors method, wherein described in step (1) GNSS observation L is mainly made of pseudorange and carrier observations, and corresponding residual vector V is solved according to the following formula:
V=Bx-L
Wherein x is to include carrier three-dimensional coordinate, the parameter estimation of ambiguity of carrier phase;B is that x corresponds to design matrix, is seen L number of measured value is n.
3. the GNSS according to claim 1 based on trial and error method positions Detection of Gross Errors method, wherein list described in step (2) Error estimator in the power of positionAre as follows:
Wherein f is excess observation number, and P is the corresponding weight matrix of observation.
CN201810481017.1A 2018-05-18 2018-05-18 A kind of GNSS positioning Detection of Gross Errors method based on trial and error method Active CN108919321B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810481017.1A CN108919321B (en) 2018-05-18 2018-05-18 A kind of GNSS positioning Detection of Gross Errors method based on trial and error method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810481017.1A CN108919321B (en) 2018-05-18 2018-05-18 A kind of GNSS positioning Detection of Gross Errors method based on trial and error method

Publications (2)

Publication Number Publication Date
CN108919321A CN108919321A (en) 2018-11-30
CN108919321B true CN108919321B (en) 2019-05-10

Family

ID=64404380

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810481017.1A Active CN108919321B (en) 2018-05-18 2018-05-18 A kind of GNSS positioning Detection of Gross Errors method based on trial and error method

Country Status (1)

Country Link
CN (1) CN108919321B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109597841B (en) * 2018-12-12 2019-10-08 中国人民解放军32021部队 A kind of target location accuracy optimization method based on many types of cartographic satellite repeated measures
CN110398758A (en) * 2019-07-24 2019-11-01 广州中海达卫星导航技术股份有限公司 Detection of Gross Errors method, apparatus, equipment and storage medium in real-time clock bias estimation
CN111077550A (en) * 2019-12-26 2020-04-28 广东星舆科技有限公司 Gross error detection method and system applied to RTD positioning of intelligent terminal
CN111505687B (en) * 2020-04-17 2021-12-21 中国科学院国家授时中心 Original observation value gross error rejection method based on GPS satellite navigation system
CN112130174B (en) * 2020-09-30 2022-10-18 长安大学 Improved GNSS-IR snow depth extraction method
CN112327340B (en) * 2021-01-06 2021-04-13 腾讯科技(深圳)有限公司 Terminal positioning accuracy evaluation method, device, equipment and medium
CN113281796B (en) * 2021-07-23 2021-10-15 腾讯科技(深圳)有限公司 Position determining method, speed determining method, device, equipment and storage medium
CN115840241B (en) * 2022-12-15 2024-02-23 北京六分科技有限公司 Partial ambiguity fixing method, apparatus, storage medium, and program product
CN116719073B (en) * 2023-08-09 2023-10-20 深圳华大北斗科技股份有限公司 GNSS (Global navigation satellite System) solution domain-oriented coarse difference detection and rejection method

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103576175B (en) * 2013-11-06 2016-01-20 西南交通大学 A kind of double frequency many constellations GNSS OTF Ambiguity Resolution method
CN106168672B (en) * 2016-01-01 2019-05-21 广州中海达卫星导航技术股份有限公司 A kind of GNSS multimode single-frequency RTK Cycle Slips Detection and device
EP3226034A1 (en) * 2016-04-01 2017-10-04 Centre National d'Etudes Spatiales Improved gnss receiver using velocity integration
CN107102346B (en) * 2017-06-08 2020-02-07 中国电子科技集团公司第五十四研究所 Multi-antenna attitude measurement method based on Beidou system
CN107728180B (en) * 2017-09-05 2021-01-29 西南交通大学 GNSS precision positioning method based on multi-dimensional particle filter deviation estimation
CN107678050B (en) * 2017-09-05 2020-09-18 西南交通大学 GLONASS phase inter-frequency deviation real-time tracking and precise estimation method based on particle filtering

Also Published As

Publication number Publication date
CN108919321A (en) 2018-11-30

Similar Documents

Publication Publication Date Title
CN108919321B (en) A kind of GNSS positioning Detection of Gross Errors method based on trial and error method
CN108169774B (en) Multimode GNSS single-frequency cycle slip detection and repair method supporting RTPPP and RTK
CN106168672B (en) A kind of GNSS multimode single-frequency RTK Cycle Slips Detection and device
US11709281B2 (en) High-precision point positioning method and device based on smartphone
CN104714244B (en) A kind of multisystem dynamic PPP calculation methods based on robust adaptable Kalman filter
Matosevic et al. A comparison of accuracy using a GPS and a low-cost DGPS
CN108363084A (en) Utilize the method and apparatus of satellite positioning, satellite navigation receiver, storage medium
CN105629279B (en) A kind of wide lane ambiguity fixing means between Fiducial station of the network
CN105301617B (en) A kind of integer ambiguity validity check method in satellite navigation system
CN102116867A (en) Method for detecting and restoring cycle slip of GPS (Global Positioning System) carrier phase under dynamic environment
CN109143298B (en) Beidou and GPS observation value cycle slip detection and restoration method, equipment and storage equipment
CN109116394A (en) A kind of real-time dynamic positioning method suitable for different length baseline
CN102998681A (en) High-frequency clock error estimation method of satellite navigation system
CN105158783A (en) Real-time dynamic differential positioning method and device thereof
CN105158778B (en) Multisystem combined implementation carrier phase difference fault satellites elimination method and its system
TW201445168A (en) A receiver and method for satellite positioning and speed measuring
CN105388496B (en) Traffic application vulnerability checking system and method based on GPS
CN115792980A (en) Model and data dual-driven GNSS RTK positioning satellite selection method and system
CN114994727A (en) Equipment for realizing high-precision time calibration and satellite positioning
Obst et al. Probabilistic non-line-of-sight detection in reliable urban GNSS vehicle localization based on an empirical sensor model
CN112540393A (en) Offshore precise single-point positioning protection horizontal quality checking method and system
CN115616637B (en) Urban complex environment navigation positioning method based on three-dimensional grid multipath modeling
CN106443729A (en) Global navigation satellite system (GNSS) adaptive pseudo range-phase weight ratio determination method
Yin et al. Modelling and assessment of Galileo and Galileo/GPS velocity determination with stand-alone receiver
CN111999750B (en) Real-time single-station cycle slip detection improvement method aiming at inaccurate lever arm

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
GR01 Patent grant
GR01 Patent grant