CN110208841A - A kind of improved GNSS tight integration method towards non-overlapping frequency - Google Patents

A kind of improved GNSS tight integration method towards non-overlapping frequency Download PDF

Info

Publication number
CN110208841A
CN110208841A CN201910559647.0A CN201910559647A CN110208841A CN 110208841 A CN110208841 A CN 110208841A CN 201910559647 A CN201910559647 A CN 201910559647A CN 110208841 A CN110208841 A CN 110208841A
Authority
CN
China
Prior art keywords
disb
fuzziness
poor
ambiguity
list
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.)
Granted
Application number
CN201910559647.0A
Other languages
Chinese (zh)
Other versions
CN110208841B (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering 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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201910559647.0A priority Critical patent/CN110208841B/en
Publication of CN110208841A publication Critical patent/CN110208841A/en
Application granted granted Critical
Publication of CN110208841B publication Critical patent/CN110208841B/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
    • 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

Abstract

The invention discloses a kind of improved GNSS tight integration methods towards non-overlapping frequency, belong to satellite precise navigation and field of locating technology;The present invention passes through building GNSS pseudorange and carrier phase list difference observational equation;Establish DISB estimation model;Construct CMC combination observation amount and DISB correction model;Then traversing operation is executed to single poor Ambiguity Search Space, IFDB Orders Corrected will be generated and successively substituted into DISB correction model, execute the fixed algorithm of LAMBDA fuzziness, and count fuzziness suboptimum and optimal value rate value;Select the poor values of ambiguity of list corresponding to maximum value as optimal solution from ratio sequence again;It is finally substituted into DISB correction model using the poor fuzziness of the optimal list resolved, obtains optimal fuzziness and positioning calculation value.The present invention solves the accurate poor ambiguity issue of list for solving reference satellite, improves ambiguity resolution success rate and positioning performance, and finally effectively eliminate influence of the IFDB to RTK ambiguity resolution.

Description

A kind of improved GNSS tight integration method towards non-overlapping frequency
Technical field
The present invention relates to satellite precise navigation and positioning field, and in particular to a kind of improved towards non-overlapping frequency GNSS tight integration method.
Background technique
Carrier Phase Ambiguity Resolution is real-time dynamic (Real-time kinematic, RTK) carrier phase differential positioning method Obtain the key technology of high-precision with high reliability.In the case where star constrained environment is received in city, valley etc., because observed quantity number can be used Deficiency, Carrier Phase Ambiguity Resolution and the positioning performance of traditional RTK method will be adversely affected.With global navigational satellite system Multisystem combination is sufficiently excavated in the modernization construction of system (Global navigation satellite system, GNSS) Technical advantage can promote Carrier Phase Ambiguity Resolution success rate in the case where receiving star constrained environment for RTK method and positioning accuracy provides Opportunity.
In the research of GNSS multisystem combination, based on deviation (differential between receiver end differential system Inter inter-system bias, DISB) research with characteristic steady in a long-term, GNSS tight integration method is suggested.The party Method is directed to all GNSS systems and only chooses difference model between 1 reference satellite building star, to the DISB introduced in model using pre- First calibration mode is removed, so that the more traditional RTK method of model redundancy is increased, to help to promote fuzziness Computation success and positioning accuracy.Currently, the GNSS tight integration method towards overlaid frequency is in the case where receiving star constrained environment It is widely used in promoting fuzziness computation success and positioning performance.However, for non-overlapping frequency, except needing to remove Except DISB, it is also necessary to consider and difference on the frequency dependent deviation (inter-frequency differential bias, IFDB) Problem.IFDB depends on the product of the list poor fuzziness and difference on the frequency of reference satellite.For specified GNSS system tight integration, If BDS B1 is combined with Galileo E1, the difference on the frequency of system can not change.Therefore, the list for accurately solving reference satellite is poor Fuzziness becomes one of the effective ways for eliminating that IFDB influences.Conventional method is by subtracting pseudorange (Carrier- to carrier phase Minus-code combination, CMC) combination observation amount be rounded nearby solve reference satellite the poor fuzziness of list, resolve Fuzziness precision is limited to pseudorange observation accuracy of measurement.When pseudorange multipath error is larger, the biggish poor fuzziness deviation of list will IFDB elimination can be adversely affected.Currently, the problem of being directed to the list poor fuzziness for how accurately solving reference satellite is also It is not solved effectively, so that the GNSS tight integration method towards non-overlapping frequency can not play its advantage, is mentioned to the greatest extent Rise RTK ambiguity resolution and positioning performance.
Summary of the invention
The purpose of the present invention is to provide a kind of improved GNSS tight integration method towards non-overlapping frequency, emphasis is solved The accurate list poor ambiguity issue for solving reference satellite, effectively eliminates so that IFDB is available, improves ambiguity resolution Success rate and positioning performance.Present invention improves over traditional, and the GNSS tight integration method towards non-overlapping frequency solves reference satellite The poor fuzziness strategy of list, be first rounded method nearby by CMC and determine single poor Ambiguity Search Space, then the single poor fuzziness of building Search criteria resolves the poor fuzziness of optimal list, finally effectively eliminates influence of the IFDB to RTK ambiguity resolution.Tool of the invention Steps are as follows for body execution:
A kind of improved GNSS tight integration method towards non-overlapping frequency, executing step includes:
Step 1: building GNSS pseudorange and carrier phase list difference observational equation, and using website general location to equation into Row linearization process;
Step 2: in system initialisation phase, DISB estimation model is established, is determined using the data of specified initialization time The calibration value of DISB;
Step 3: building CMC combination observation amount determines that the poor ambiguity search of list of reference satellite is empty by being rounded method nearby Between blur center angle value determine the boundary of search space and according to the pseudorange observation precision of priori;
Step 4: the DISB calibration value determined using step 2 constructs DISB correction model;
Step 5: traversing operation is executed to single poor Ambiguity Search Space, IFDB Orders Corrected is generated and successively substitutes into DISB In correction model, the fixed algorithm of LAMBDA fuzziness is executed, and count fuzziness suboptimum and optimal value rate value;
Step 6: in the ratio sequence generated in step 5, select the poor values of ambiguity of list corresponding to maximum value as Optimal solution;
Step 7: the poor fuzziness of optimal list resolved using step 6 is substituted into the DISB correction model in step 4, is obtained Obtain optimal fuzziness and positioning calculation value.
The poor observed quantity model of list in the step 1 is expressed as:
Wherein, B refers to BDS, and E refers to Galileo;P and φ indicates that pseudorange and carrier phase observed quantity subtract the list of calculation amount Poor residual volume;X indicates that positional increment vector, corresponding geometric matrix are u;T indicates receiver clock-offsets;D and δ respectively indicate pseudorange With the hardware delay error of carrier phase;Z is values of ambiguity;λ indicates wavelength;E and ε respectively corresponds the sight of pseudorange and carrier phase Survey noise.
DISB in the step 2 estimates model are as follows:
WhereindBEIndicate the DISB of corresponding pseudorange;δBEIndicate the DISB of corresponding carrier phase;IFDB is referred to as in the present invention.
The blur center angle value of search space in the step 3 are as follows:
Wherein " []round" indicate nearest floor operation,For center values of ambiguity;It is empty for single poor ambiguity search Between, it is expressed asWhereinIndicate that set of integers, N are space search boundary.
DISB correction model in the step 4 are as follows:
Fuzziness suboptimum and optimal value rate value in the step 5 indicate are as follows:
Wherein, " | | | | " indicate quadratic form operation;It is expressed as fuzzinessCovariance.
The beneficial effects of the present invention are:
The present invention solves the accurate poor ambiguity issue of list for solving reference satellite, effectively disappears so that IFDB is available It removes, improves ambiguity resolution success rate and positioning performance;Improve traditional tight integration side GNSS towards non-overlapping frequency Method solves the poor fuzziness strategy of list of reference satellite, is first rounded method nearby by CMC and determines single poor Ambiguity Search Space, then The single poor fuzziness search criteria of building resolves the poor fuzziness of optimal list, finally effectively eliminates IFDB to RTK ambiguity resolution It influences.
Detailed description of the invention
Fig. 1 is the flow chart of the improved GNSS tight integration method towards non-overlapping frequency of the present invention.
Specific embodiment
It is an object of that present invention to provide a kind of improved GNSS tight integration method towards non-overlapping frequency, this method is main Solve traditional GNSS tight integration method towards non-overlapping frequency can not accurately solve reference satellite the poor values of ambiguity of list and It is difficult to eliminate the problem of IFDB influences RTK ambiguity resolution.This method subtracts pseudorange (CMC) method by traditional carrier phase Single poor Ambiguity Search Space has been determined, search criteria is turned to fuzziness suboptimal solution and optimal solution ratio maximum, is determined optimal The poor values of ambiguity of list.This method corrects IFDB using the poor values of ambiguity of the optimal list, so as to so that towards non-overlap frequency The positioning performance of the GNSS tight integration of rate reaches best.It, should compared to more traditional GNSS tight integration method towards non-overlapping frequency Method can obtain better satellite navigation and positioning performance in the case where receiving the limited urban canyon environment of star.It is right with reference to the accompanying drawing The present invention is described further.
Embodiment 1:
Step 1: building GNSS pseudorange and carrier phase list difference observational equation, and equation is carried out using website general location Linearization process.
Base station and movement station track BDS and Galileo satellite s simultaneously.Single poor operation can eliminate satellite end correlation between standing Error.Only consider under short baseline case simultaneously, atmosphere delay error is had ignored in model.Single poor observed quantity model can indicate For,
Wherein, B refers to BDS, and E refers to Galileo;P and φ indicates that pseudorange and carrier phase observed quantity subtract the list of calculation amount Poor residual volume;X indicates that positional increment vector, corresponding geometric matrix are u;T indicates receiver clock-offsets;D and δ respectively indicate pseudorange With the hardware delay error of carrier phase;Z is values of ambiguity;λ indicates wavelength;E and ε respectively corresponds the sight of pseudorange and carrier phase Survey noise.
Step 2: in system initialisation phase, establishing DISB estimation model, determined using the data of specified initialization time The calibration value of DISB;
In the system initial stage, need to determine the calibration value of DISB using specified time segmental arc.So according to step 1 The poor observed quantity of the list provided, construction DISB estimation model are
WhereindBEIndicate the DISB of corresponding pseudorange;δBEIndicate the DISB of corresponding carrier phase.IFDB is referred to as in the present invention.Formula (2) estimates pseudorange with carrier phase DISB as parameter. Since DISB has characteristic steady in a long-term, calibration value of the mean value of estimated sequence as DISB in segmental arc so fetching is fixed time. Specified time segmental arc is set in the present invention as 30s.
Step 3: building CMC combination observation amount determines that the poor ambiguity search of list of reference satellite is empty by being rounded method nearby Between blur center angle value determine the boundary of search space and according to the pseudorange observation precision of priori.
Pseudorange and the difference observed quantity of carrier phase list using reference satellite construct CMC combination observation amount, and by taking nearby Whole method determines that the blur center angle value of its search space is
Wherein " []round" indicate nearest floor operation.It is false in the present invention in order to cover single poor fuzziness true value If pseudorange observation error is 10m, corresponding search space boundary will be set as N=50 weeks.In this way, single differential mode of reference satellite is pasted Spend search spaceIt can be expressed asWhereinIndicate set of integers,Centered on Values of ambiguity.
Step 4: the DISB calibration value determined using step 2 constructs DISB correction model.
The DISB calibration value and candidate IFDB correction amount that step 2 is determined substitute on the left of formula as known quantity, obtain Accordingly DISB correction model is
Step 5: traversing operation being executed to single poor Ambiguity Search Space, generation IFDB Orders Corrected successively substitutes into DISB and repairs In positive model, the fixed algorithm of LAMBDA fuzziness is executed, and count fuzziness suboptimum and optimal value rate value.
By traverse fuzziness search space determine IFDB correction amount, can obtain 2N fuzziness suboptimum with it is optimal Rate value.
Wherein, " | | | | " indicate quadratic form operation;It is expressed as fuzzinessCovariance.
Step 6: in the ratio sequence generated in steps of 5, selecting the poor values of ambiguity of list corresponding to maximum value as most Excellent solution.
Step 7: the poor fuzziness of optimal list resolved using step 6 is substituted into the DISB correction model in step 4, is obtained most Excellent fuzziness and positioning calculation value.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.

Claims (6)

1. a kind of improved GNSS tight integration method towards non-overlapping frequency, which is characterized in that executing step includes:
Step 1: building GNSS pseudorange and carrier phase list difference observational equation, and line is carried out to equation using website general location Propertyization processing;
Step 2: in system initialisation phase, DISB estimation model is established, determines DISB using the data of specified initialization time Calibration value;
Step 3: building CMC combination observation amount determines the list poor Ambiguity Search Space of reference satellite by being rounded method nearby Blur center angle value, and according to the pseudorange observation precision of priori, determine the boundary of search space;
Step 4: the DISB calibration value determined using step 2 constructs DISB correction model;
Step 5: executing traversing operation to single poor Ambiguity Search Space, generates IFDB Orders Corrected and successively substitutes into DISB amendment In model, the fixed algorithm of LAMBDA fuzziness is executed, and count fuzziness suboptimum and optimal value rate value;
Step 6: in the ratio sequence generated in step 5, select the poor values of ambiguity of list corresponding to maximum value as optimal Solution;
Step 7: the poor fuzziness of optimal list resolved using step 6 is substituted into the DISB correction model in step 4, is obtained most Excellent fuzziness and positioning calculation value.
2. the improved GNSS tight integration method towards non-overlapping frequency of one kind according to claim 1, which is characterized in that The poor observed quantity model of list in the step 1 is expressed as:
Wherein, B refers to BDS, and E refers to Galileo;P and φ expression pseudorange and the list that carrier phase observed quantity subtracts calculation amount are poor residual Surplus;X indicates that positional increment vector, corresponding geometric matrix are u;T indicates receiver clock-offsets;D and δ respectively indicates pseudorange and carries The hardware delay error of wave phase;Z is values of ambiguity;λ indicates wavelength;E and ε respectively corresponds pseudorange and the observation of carrier phase is made an uproar Sound.
3. the improved GNSS tight integration method towards non-overlapping frequency of one kind according to claim 1, which is characterized in that DISB in the step 2 estimates model are as follows:
WhereindBEIndicate the DISB of corresponding pseudorange;δBEIndicate the DISB of corresponding carrier phase;IFDB is referred to as in the present invention.
4. the improved GNSS tight integration method towards non-overlapping frequency of one kind according to claim 1, which is characterized in that The blur center angle value of search space in the step 3 are as follows:
Wherein " []round" indicate nearest floor operation;For center values of ambiguity;For single poor Ambiguity Search Space, table It is shown asWhereinIndicate that set of integers, N are space search boundary.
5. the improved GNSS tight integration method towards non-overlapping frequency of one kind according to claim 1, which is characterized in that DISB correction model in the step 4 are as follows:
6. the improved GNSS tight integration method towards non-overlapping frequency of one kind according to claim 1, which is characterized in that Fuzziness suboptimum and optimal value rate value in the step 5 indicate are as follows:
Wherein, " | | | | " indicate quadratic form operation;It is expressed as fuzzinessCovariance.
CN201910559647.0A 2019-06-26 2019-06-26 Improved GNSS tight combination method facing non-overlapping frequencies Active CN110208841B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910559647.0A CN110208841B (en) 2019-06-26 2019-06-26 Improved GNSS tight combination method facing non-overlapping frequencies

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910559647.0A CN110208841B (en) 2019-06-26 2019-06-26 Improved GNSS tight combination method facing non-overlapping frequencies

Publications (2)

Publication Number Publication Date
CN110208841A true CN110208841A (en) 2019-09-06
CN110208841B CN110208841B (en) 2022-09-02

Family

ID=67794535

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910559647.0A Active CN110208841B (en) 2019-06-26 2019-06-26 Improved GNSS tight combination method facing non-overlapping frequencies

Country Status (1)

Country Link
CN (1) CN110208841B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114355421A (en) * 2021-12-22 2022-04-15 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101902195A (en) * 2009-11-27 2010-12-01 华中科技大学 Method for modeling generator excitation system
EP2477564A1 (en) * 2009-09-02 2012-07-25 Skeletal Dynamics, LLC Internal joint stabilizer device, system and method of use
CN102736094A (en) * 2012-06-19 2012-10-17 哈尔滨工程大学 Adaptive genetic algorithm-based single-frequency GNSS (Global Navigation Satellite System) integer ambiguity acquisition method
EP2541198A1 (en) * 2011-06-30 2013-01-02 Furuno Electric Company Limited Road map feedback corrections in tightly coupled gps and dead reckoning vehicle navigation
JP2014238349A (en) * 2013-06-10 2014-12-18 古野電気株式会社 Positioning device, information apparatus terminal, positioning system, positioning method, and positioning program
CN105891864A (en) * 2016-04-29 2016-08-24 辽宁工程技术大学 Method for fixing hybrid double-difference ambiguity between BDS and GPS
CN106125113A (en) * 2016-06-20 2016-11-16 武汉大学 A kind of high accuracy Baselines method utilizing multisystem GNSS observation
CN106896386A (en) * 2017-04-25 2017-06-27 武汉大学 GLONASS inter-frequency deviation precise Estimation Methods
US20170269227A1 (en) * 2016-03-18 2017-09-21 Deere & Company Estimation of inter-frequency bias for ambiguity resolution in global navigation satellite system receivers
CN108196287A (en) * 2018-02-02 2018-06-22 东南大学 A kind of tight integration RTK localization methods that need not consider reference satellite transformation
CN108519614A (en) * 2018-03-16 2018-09-11 东南大学 A kind of GPS/BDS tight integrations carrier difference localization method
CN109507696A (en) * 2018-05-25 2019-03-22 武汉大学 A kind of GNSS pseudorange multipath removing method based on adaptive-filtering
CN109597105A (en) * 2018-12-13 2019-04-09 东南大学 A kind of GPS/GLONASS tight integration localization method for taking deviation between carrier system into account

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2477564A1 (en) * 2009-09-02 2012-07-25 Skeletal Dynamics, LLC Internal joint stabilizer device, system and method of use
CN101902195A (en) * 2009-11-27 2010-12-01 华中科技大学 Method for modeling generator excitation system
EP2541198A1 (en) * 2011-06-30 2013-01-02 Furuno Electric Company Limited Road map feedback corrections in tightly coupled gps and dead reckoning vehicle navigation
CN102736094A (en) * 2012-06-19 2012-10-17 哈尔滨工程大学 Adaptive genetic algorithm-based single-frequency GNSS (Global Navigation Satellite System) integer ambiguity acquisition method
JP2014238349A (en) * 2013-06-10 2014-12-18 古野電気株式会社 Positioning device, information apparatus terminal, positioning system, positioning method, and positioning program
US20170269227A1 (en) * 2016-03-18 2017-09-21 Deere & Company Estimation of inter-frequency bias for ambiguity resolution in global navigation satellite system receivers
CN105891864A (en) * 2016-04-29 2016-08-24 辽宁工程技术大学 Method for fixing hybrid double-difference ambiguity between BDS and GPS
CN106125113A (en) * 2016-06-20 2016-11-16 武汉大学 A kind of high accuracy Baselines method utilizing multisystem GNSS observation
CN106896386A (en) * 2017-04-25 2017-06-27 武汉大学 GLONASS inter-frequency deviation precise Estimation Methods
CN108196287A (en) * 2018-02-02 2018-06-22 东南大学 A kind of tight integration RTK localization methods that need not consider reference satellite transformation
CN108519614A (en) * 2018-03-16 2018-09-11 东南大学 A kind of GPS/BDS tight integrations carrier difference localization method
CN109507696A (en) * 2018-05-25 2019-03-22 武汉大学 A kind of GNSS pseudorange multipath removing method based on adaptive-filtering
CN109597105A (en) * 2018-12-13 2019-04-09 东南大学 A kind of GPS/GLONASS tight integration localization method for taking deviation between carrier system into account

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
LIN ZHAO 等: ""Cycle slip detection and repair with triple frequency combination method"", 《2014 IEEE/ION POSITION, LOCATION AND NAVIGATION SYMPOSIUM - PLANS 2014》 *
LIU YINHUA 等: ""The analysis of Differential Code Bias of BeiDou satellite navigation system"", 《2014 IEEE INTERNATIONAL FREQUENCY CONTROL SYMPOSIUM (FCS)》 *
XIAOYAN ZHU 等: ""Using dynamic pseudo-IDs to protect privacy in location-based services"", 《2014 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS (ICC)》 *
贾春 等: ""改正接收机频间偏差的短基线北斗三频紧组合RTK方法"", 《中国科学:地球科学》 *
隋心 等: ""GPS/BDS接收机端系统偏差稳定性对整周模糊度固定的影响"", 《武汉大学学报(信息科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114355421A (en) * 2021-12-22 2022-04-15 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value
CN114355421B (en) * 2021-12-22 2023-08-18 杭州电子科技大学 Flood detection method based on Beidou satellite L4 and CMC combined observation value

Also Published As

Publication number Publication date
CN110208841B (en) 2022-09-02

Similar Documents

Publication Publication Date Title
CN110007320B (en) Network RTK resolving method
CN108415049B (en) Method for improving network RTK double-difference wide lane ambiguity fixing accuracy
CN108169774B (en) Multimode GNSS single-frequency cycle slip detection and repair method supporting RTPPP and RTK
CN111308528B (en) Positioning method for Beidou/GPS tightly-combined virtual reference station
CN109477900A (en) Estimation in global navigational satellite system receiver for the frequency offset of ambiguity resolution
CN111505685B (en) Positioning method of multisystem combination RTK model based on correcting intersystem deviation
CN111239787A (en) GNSS dynamic Kalman filtering method in cluster autonomous coordination
JPH10221109A (en) Method and apparatus for obtaining observable value for input to kalman filter process for judging rolling and pitching of platform and bearing of nose
CN111998849A (en) Differential dynamic positioning method based on inertial navigation system
CN112462397A (en) Real-time dynamic positioning method and system for full-constellation multi-frequency Beidou data
CN109683182A (en) The GNSS multimode of deviation combines Differential positioning method between a kind of introducing system
CN113325453A (en) GNSS non-differential ambiguity determination method based on parameter constraint and rapid positioning method
CN114935770B (en) Method and device for accelerating precision single-point positioning convergence speed by multiple calendars
CN110568464B (en) BDS/GNSS multimode chip-based precise positioning method and BDS/GNSS multimode chip-based precise positioning device
CN112146557A (en) GNSS-based real-time bridge deformation monitoring system and method
CN110208841A (en) A kind of improved GNSS tight integration method towards non-overlapping frequency
WO2020044129A1 (en) Attitude determination based on global navigation satellite system information
CN114114334A (en) GLONASS inter-frequency deviation calibration method and RTK method
JP4723801B2 (en) Relative positioning device
RU2624268C1 (en) Method of determining mutual position of objects by signals of global navigation satellite systems
CN110095796A (en) A kind of real-time dynamic precision navigation locating method of cross-system MW tight integration
CN113341686B (en) Single-satellite multi-calendar timing method and device, computer equipment and storage medium
CN110376629B (en) Method for determining deviation between satellite differential systems based on principle of maximum RATIO value
CN115407379A (en) Ambiguity resolution method based on pseudo satellite and UWB fusion positioning system
CN114114357A (en) Multi-constellation dynamic data high-performance quality control method

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