CN101710179A - Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method - Google Patents

Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method Download PDF

Info

Publication number
CN101710179A
CN101710179A CN200910273365A CN200910273365A CN101710179A CN 101710179 A CN101710179 A CN 101710179A CN 200910273365 A CN200910273365 A CN 200910273365A CN 200910273365 A CN200910273365 A CN 200910273365A CN 101710179 A CN101710179 A CN 101710179A
Authority
CN
China
Prior art keywords
frequency
frequently
gnss
motion
observation
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
CN200910273365A
Other languages
Chinese (zh)
Other versions
CN101710179B (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.)
Weihai Wuzhou Satellite Navigation Technology Co., Ltd.
Original Assignee
Wuhan University WHU
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 Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN 200910273365 priority Critical patent/CN101710179B/en
Publication of CN101710179A publication Critical patent/CN101710179A/en
Application granted granted Critical
Publication of CN101710179B publication Critical patent/CN101710179B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention relates to a global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method. In the original epoch, the triple-frequency precision single-point positioning technology is adopted to obtain the coordinates of two movable carriers, one of the two movable carries is selected as the reference station for the triple-frequency double-difference position, the triple-frequency double-difference positioning technology is adopted to calculate the baseline component of the two movable carriers and the integer ambiguity resolution of the double-difference carrier phase, and the baseline component of the two movable carriers and the integer ambiguity resolution of the double-difference carrier phase are used as the constraint conditions for the triple-frequency precision single-point positioning in the subsequent epochs to improve the single-point position precision and the convergence speed thereof. The geometry-Base TCAR ambiguity resolution fixation method is adopted to calculate the triple-frequency integer ambiguity resolution. The three irrelevant combination observation values of the triple-frequency non-ionizing layer and the long wave-length and low noise carrier are adopted to detect and repair the cycle slip of the original observation data.

Description

A kind of GNSS three is moving frequently to motion positioning method
Technical field
The invention belongs to satellite navigation positioning technical field, especially GNSS three is moving frequently to motion positioning method.
Background technology
The precision that the Global Navigation Satellite System (GNSS) carrier phase difference divides dynamic location technology can reach centimetre-sized is widely used in the high precision navigator fix field.The difference Kinematic Positioning of traditional but " quiet-moving " pattern needs the static reference station that coordinate is known, obtains the mobile observation station location thereby mobile observation station and static reference station form difference modes.And in a lot of occasions, what need care is the relative position of two motion carriers, decides appearance etc. as air refuelling, formation flight, aircraft butt joint, precision.For example shown in Figure 1, can utilize GNSS moving to moving location, the carrier signal that provides is via satellite determined two airplane accurate relative position relation.
This class does not have the dynamic phasing orientation problem of static reference station, is referred to as " moving to moving location ", moving moving location is had a wide range of applications, but domesticly for moving moving navigator fix technical research is remained in blank.The GNSS carrier phase difference divides dynamically, and the location reaches the centimetre-sized precision, need to obtain reliable ambiguity of carrier phase, to provide the satellite-signal of three frequencies for its user after the GPS of America modernization, three frequencies can increase the number of observed reading, obtain more three and make up no ionosphere frequently, the dummy observation of long wavelength and low noise carrier wave, for Carrier Phase Ambiguity Resolution provides new method, reduce the integer ambiguity set time, improve its success ratio, mainly contain for three frequency blur level study on fixation methods both at home and abroad at present: Wu Yue, the theoretical and application [D] of second generation satellite navigation system multifrequency data processing. the .2005 of Wuhan University; Yanming Feng, Chris Rizos.Geometry-Based TCAR Models and Performance Analysis[J] .IUGG 2007 .July2-13, Perugia, Italy; Wentao Zhang.Triple Frequency Cascading AmbiguityResolution for Modernized GPS and GALILEO[D] .Calgary.2005; A NewThree-Frequency, Geometry-Free Technique for Ambiguity Resolution[C] .Proc ofION GNSS 2006, P309-316,26-29 Sept, Fort Worth, TX..For above three frequency blur level fixing meanss, adopt Geometry-free TCAR or Geometry-Based TCAR technology to resolve whole integer ambiguities, but these two kinds of technology existing problems be when observation condition relatively poor, or the observation satellite number more for a long time, resolves all integer ambiguity values and can fail.
Summary of the invention
The objective of the invention is to overcome the prior art defective, provide a kind of new GNSS three moving frequently motion positioning method.
GNSS three provided by the present invention is moving frequently to motion positioning method, and technical scheme is as follows:
In initial epoch, carry out following steps,
Step 1 as the mobile observation station, gets two motion carriers three frequency GNSS data to two original observation stations in mobile observation station and carries out pre-service;
Step 2 is according to the pretreated three frequency GNSS data of step 1 gained, by the coordinate at two mobile observation stations of three Static Precise Point Positioning acquisitions frequently;
Step 3, establishing arbitrary in two mobile observation stations is the reference station of three two frequently differences location, according to the coordinate at two mobile observation stations of step 2 gained, finds the solution the baseline component at two mobile observation stations and the integer ambiguity of two poor carrier phases by three two frequently differences location;
In follow-up epoch, two original observation stations in mobile observation station are got three frequency GNSS data carry out pre-service, according to pretreated three frequency GNSS data, by the coordinate at two mobile observation stations of three Static Precise Point Positioning acquisitions frequently, and with the constraint condition of the step 3 gained solving result in initial epoch as three frequency Static Precise Point Positioning.
And, described to two original observation stations in mobile observation station get three frequently the GNSS data carry out pre-service, comprise the uncorrelated combination observation value that adopts three three not have ionosphere, long wavelength and low noise carrier wave frequently, two original observation stations in mobile observation station are got three frequency GNSS data carry out the cycle slip detection and repair.
And, when carrying out three frequency Static Precise Point Positioning, adopt three the three uncorrelated estimation of combination observation values troposphere, receiver location and the motion carrier clock correction parameters of not having ionosphere, long wavelength and low noise carrier wave frequently.And, during find the solution the two frequently differences of step 3 described three location, finding the solution of the integer ambiguity of described two difference carrier phases adopts Geometry-Based TCAR method to realize, then according to success ratio thresholding fixed part blur level, adopts assumed statistical inspection to judge the reliability of fixed part blur level.
And, describedly be according to success ratio thresholding fixed part blur level implementation, obtain the real solution and the variance and covariance battle array thereof of blur level by two poor location models, calculate the fixing success ratio of blur level according to the variance and covariance battle array behind the integer transform, get one group of blur level of its success ratio maximum in default success ratio domain value range.
Baseline component and values of ambiguity constraint Static Precise Point Positioning thereof that the present invention adopts two poor location models to obtain improve the precision of Static Precise Point Positioning and reduce its convergence time, thereby utilize the moving to moving location of three high-accuracy high-efficiency of GNSS data realization frequently rates.The present invention also adopts Geometry-Based TCAR method to determine blur level, further improves counting yield, and by the fixed part blur level, improves the baseline computational accuracy.
Description of drawings
Fig. 1 is moving to moving navigator fix synoptic diagram for GNSS;
Fig. 2 is that the GNSS three of the embodiment of the invention is moving frequently to moving positioning flow figure;
Fig. 3 is the GNSS three frequency Static Precise Point Positioning process flow diagrams of the embodiment of the invention;
Fig. 4 is three frequency observed reading linear combination observed reading wavelength graphs;
Fig. 5 is the GNSS three cycle slip detection frequently process flow diagram of the embodiment of the invention;
Fig. 6 is the GNSS three frequency ambiguity resolution process flow diagrams of the embodiment of the invention;
Fig. 7 is the blur level partial fixing process flow diagram of the embodiment of the invention.
Embodiment
Below in conjunction with drawings and Examples technical solution of the present invention is described:
Three frequencies that embodiment relates to i.e. three frequency ranges, are designated L1/L2/L5 respectively.The GNSS that the embodiment of the invention provides is moving to moving positioning flow as shown in Figure 2:
In initial epoch, carry out following three steps,
Step 1, with two motion carriers respectively as mobile observation station 1 and mobile observation station 2, under precise ephemeris and clock correction, provide original observation station to get three frequency GNSS data, then three original frequency GNSS data are carried out pre-service, comprise that three audio data cycle slips are surveyed and reparation.
Step 2, according to pretreated mobile observation station 1 gained, three frequency GNSS data, convection current in-motion viewing survey station 1 carries out three frequency Static Precise Point Positioning; According to pretreated mobile observation station 2 gained, three frequency GNSS data, convection current in-motion viewing survey station 2 carries out three frequency Static Precise Point Positioning.By existing three frequency Static Precise Point Positioning technology, can ask for obtaining two more accurate coordinates in mobile observation station.
Step 3, establishing arbitrary in two mobile observation stations is the reference station of three two frequently differences location, according to the coordinate at two mobile observation stations of step 2 gained, finds the solution the baseline component at two mobile observation stations and the integer ambiguity of two poor carrier phases by three two frequently differences location.Suggestion is concrete when implementing, with the higher mobile observation station of coordinate precision in two motion carriers as with reference to the station, promptly with the coordinate of this motion carrier as the reference station coordinates.In three two frequently differences location, at first ask for two difference observed readings according to pretreated mobile observation station 1,2 gained, three frequency GNSS data, obtain blur level real solution and variance and covariance battle array thereof according to reference station coordinates by two difference location models then, adopt blur level searching method search blur level, adopt the assumed statistical inspection of integer ambiguity optimum solution and suboptimal solution to judge the Carrier Phase Ambiguity Resolution correctness.If assumed statistical inspection is passed through, then utilize the blur level static solution to obtain the baseline vector, and with blur level static solution and integer ambiguity value thereof as the constraint condition of carrying out three frequency Static Precise Point Positioning once more; If check failure, then with blur level real solution and baseline component thereof as the constraint condition of carrying out three frequency Static Precise Point Positioning once more.
In follow-up epoch, under step 3 gained constraint condition, carry out Static Precise Point Positioning.Because the Static Precise Point Positioning in all follow-up epoch can be utilized the result of three two frequently differences location in initial epoch, therefore can be considered two mobile observation stations unites and carries out the single-point location, to carry out the precision of single-point location gained coordinate higher for the ratio of precision step 2 of single-point location gained coordinate in the case, thereby improved bearing accuracy, reduced convergence time.In follow-up epoch, the coordinate at single-point gained mobile observation station 1, location and mobile observation station 2 can directly be exported moving positioning result as moving.
Specific implementation of the present invention can for the sake of ease of implementation, provide the specific implementation of embodiment as follows referring to " GPS measures and data processing " (Li Zhenghang etc. write):
(1) original observation station gets three pre-service of GNSS data frequently
Embodiment three frequently the Static Precise Point Positioning flow processs as shown in Figure 3, process is: before carrying out the single-point location, need convection current in-motion viewing survey station 1 original observation station get three frequently GNSS data and mobile observation station 2 original observation stations get three frequently the GNSS data carry out pre-service respectively; To data after the pre-service, adopt the combination observation value estimation of no ionosphere troposphere, receiver location and motion carrier clock correction parameter then, thereby obtain three-dimensional coordinate, clock correction, the troposphere parameter at mobile observation station 1,2, realize three frequency Static Precise Point Positioning.Single-point location during initial epoch, with the single-point location in follow-up epoch, difference only is that the latter has in initial epoch two differences location find the solution the gained baseline and ambiguity of carrier phase participates in, and adopts frame of broken lines to represent among Fig. 3.As seen at first to solve the problem of pretreatment that original observation station gets three frequency GNSS data.
The detection of cycle slip is the key link in the GPS dynamic navigation location, can be correct detect cycle slip, will directly influence the precision and the reliability of navigator fix, thus the locator data pre-service all can comprise three frequently cycle slips survey and repair.And used three audio data of the present invention have bigger advantage in the detection to cycle slip with reparation.Under the situation for three frequency observed readings, can adopt wavelength longer, the more weak observating characteristic of ionosphere influence three uncorrelated combination observation values preferably carries out the detection of cycle slip, when original observation station gets carrier phase observation data and has cycle slip, to in three uncorrelated combination observation values, amplify, utilize pseudorange information to determine the size of three uncorrelated combination observations jumping on weekly duty, then can calculate the cycle slip size that obtains three original frequency carrier phases.
The GNSS three frequency cycle slip detection methods of the embodiment of the invention as shown in Figure 5, process forms three three frequencies and linear incoherent carrier phase combination observation value for get three frequency GNSS data (comprising carrier phase observation data and pseudorange raw observation) by mobile observation station 1,2 original observation stations.In order to reach good observating characteristic, the embodiment of the invention requires these uncorrelated combination observation values to have the character of no ionosphere, long wavelength's (wavelength is generally greater than 86 centimetres) and low noise carrier wave (noise is generally less than the observation noise of L1 carrier wave), and obtain corresponding combination observation value for (0,1 ,-1) according to these character conditions, (3,1,3), (1,8 ,-7).Begin then carrying out difference detecting cycle slip epoch current epoch, at first whether the judgement of cycle slip takes place, if certain linear combination observed reading is in epoch, difference surpasses certain default thresholding (for example being made as 1 phase wave length) between difference carrier phase place and the level and smooth pseudorange observed reading, then think the generation cycle slip, otherwise directly to carrying out same difference detecting cycle slip epoch next epoch.If generation cycle slip, then utilize level and smooth pseudorange observed reading to determine the cycle slip of above-mentioned three uncorrelated combination observation values, calculate original observation station and get the carrier phase observation data cycle slip, carry out the cycle slip reparation afterwards to carrying out same difference detecting cycle slip epoch next epoch according to the valuation of calculating the gained cycle slip.Described level and smooth pseudorange observed reading is handled the pseudorange raw observation with carrier phase observation data and is obtained, and described difference carrier phase place asks difference to obtain with the carrier phase observation data of adjacent epoch.
According to the GNSS standard, the reference frequency of GNSS carrier phase linear L1, L2, L5 all is 10.23MHz, and therefore frequency separately can be expressed as: f i=m if 0, i gets 1,2,5 herein.Reference frequency f wherein 0=10.23MHz (megahertz), wavelength X 0=29.31m (rice).The frequency of combination observation value can be rewritten as: f=f 0(i 1m 1+ i 2m 2+ i 5m 3)
If k=is (i 1m 1+ i 2m 2+ i 5m 3), k has integer characteristic, then f=kf 0, wavelength is: λ = λ 0 k . Coefficient i in k 1, i 2, i 5Get different values, the wavelength of combination observation value is also inequality, according to the span of k, gives the different zone of combination observation value defined, Kuan Xiang territory: 0<k<115, inter-road territory: 115≤k≤154, territory, narrow lane: k>154 respectively.Referring to accompanying drawing 4, get the wavelength of the round values combination observation value in [1-17] scope, when k=1, the maximal value of combination observation value wavelength is 29.31m (rice).
(2) three frequency Static Precise Point Positioning at first illustrate embodiment frequency three in initial epoch Static Precise Point Positioning:
According to prior art, the observation equation of GNSS single-point location is:
P i=ρ-c·δt s+c·δt t+δtrop+δion+δP mul+δrel+εp (4)
Figure G2009102733650D00072
In the formula, c is the light velocity, δ t sBe satellite clock correction, δ t tBe receiver clock correction, i represents the frequency of observed reading, P iBe the pseudorange raw observation,
Figure G2009102733650D00073
It is carrier phase observation data.ρ represents the distance of satellite to motion carrier, and δ trop represents the tropospheric delay error, and δ ion represents ionosphere delay error, δ P MulMultipath error on the expression pseudorange observed reading, δ rel represents that relativistic effect corrects, and ε p represents observation noise, and λ represents the carrier phase wavelength, and N represents phase ambiguity,
Figure G2009102733650D00074
Multipath error on the expression carrier phase observation data.
Precise ephemeris and clock correction data product (will have IGS GNSS precise ephemeris and clock correction file importing single-point positioning calculation process now gets final product) thereof that the GNSS three frequency Static Precise Point Positioning of embodiment utilize IGS to announce, adopt no ionosphere, the uncorrelated combination observation value of long wavelength and low noise carrier wave (abbreviating no ionosphere combination observation value among the present invention as) is eliminated the item influence of single order ionosphere as observed quantity, adopt Kalman filtering technique estimation troposphere parameter, location parameter and receiver clock correction, and correction solid tide, oceanic tide, antenna phase center, phase place is twined equal error, in the present embodiment in the uncorrelated combination observation value of GNSS three employing of Static Precise Point Positioning frequently and the data pre-service cycle slip survey used identical with modification, the combination coefficient of these combination observation values is (0,24 ,-23).
Then, illustrate that embodiment carries out three frequency Static Precise Point Positioning in follow-up epoch:
After additional baseline and the blur level constraint condition, the function model of two mobile observation station three frequency Static Precise Point Positioning the following is:
Figure G2009102733650D00081
In the formula, p IfA i,
Figure G2009102733650D00082
p IfB i,
Figure G2009102733650D00083
The pseudorange and the carrier phase that are respectively two mobile observation stations do not have ionosphere combination observation value, i=1, and 2...n, i ≠ j, n are number of satellites, p IfA j,
Figure G2009102733650D00084
p IfB j,
Figure G2009102733650D00085
Being respectively two selected pseudorange and carrier phases with reference to satellite in mobile observation station does not have ionosphere combination observation value, and (Δ X, Δ Y, Δ Z),
Figure G2009102733650D00086
Be respectively baseline component and two poor values of ambiguity thereof at two mobile observation stations, M A i,
Figure G2009102733650D00088
M B iBe respectively the surplus profound and troposphere projection function of satellite direction at two mobile observation stations,
Figure G2009102733650D000810
M A j,
Figure G2009102733650D000811
M B jBe respectively satellite direction cosine and the troposphere projection function of two mobile observation stations, (X with reference to satellite A, Y A, Z A), (X B, Y B, Z B) be respectively the coordinate at two mobile observation stations, δ T A, δ T BBe respectively the receiver clock correction at two mobile observation stations, δ trop A, δ trop BBe respectively the troposphere parameter at two mobile observation stations, N A j, N B jBe respectively the satellite values of ambiguity at two mobile observation stations, N A j, N B jBe respectively the values of ambiguity of the reference star at two mobile observation stations, ξ is an observation noise.
After carrying out three two frequently differences location initial epoch, just can be in follow-up epoch all with baseline and carrier phase ambiguity as constraint condition, according to the parameter estimation of above-mentioned function model, obtain three-dimensional coordinate, clock correction, the troposphere parameter at two mobile observation stations to no ionosphere combination observation value.Other with initial epoch in three identical treatment steps in the Static Precise Point Positioning frequently, the present invention will not give unnecessary details.
(3) three in initial epoch, two frequently differences were located
According to prior art, the two difference of GNSS observation equation is:
▿ Δ P i = ▿ Δρ + ▿ ΔT + ▿ ΔI f i 2 + ϵ ▿ Δ P i - - - ( 1 )
λ i ▿ Δ φ i = ▿ Δρ + ▿ ΔT - ▿ ΔI f i 2 + λ i ▿ Δ N i + ϵ ▿ Δ φ i - - - ( 2 )
I represents the frequency of observed reading, λ iBe the wavelength of correspondence, P iBe the pseudorange raw observation, φ iBe carrier phase observation data, ρ is the distance of satellite to motion carrier, Expression is two poor,
Figure G2009102733650D00094
Be two poor tropospheres residual error,
Figure G2009102733650D00095
Be two poor Ionosphere Residual Error, Be two poor blur leveles,
Figure G2009102733650D00097
Expression pseudorange observation noise,
Figure G2009102733650D00098
Expression carrier phase observation noise, f iThe expression signal frequency.When base length during less than 20km (kilometer), two poor tropospheres residual error
Figure G2009102733650D00099
With two poor Ionosphere Residual Error Can ignore.
In the GPS relative positioning, the influence that starting point coordinate changes can be divided into two parts, and one is to make the baseline vector produce translation, and in rectangular coordinate system in space, this influence relation is simple, and its value only depends on the variation of starting point coordinate; Its two, be the deviation of starting point coordinate, by GPS relative positioning model, to the influence of asking baseline vector.This influence is a base length, orientation, initial point position and all multifactor complicated functions such as the satellite geometry distribution of surveying.
Starting point coordinate changes, to the influence of surveys baseline, mainly and base length closely related, under worst situation, the size of its influence, can estimate by following approximation relation formula:
δs=0.60·10 -4·D·δX (3)
D is a base length, with km (km) is unit, and δ X is for starting at the point coordinate error, and δ s starts at the influence of point coordinate error to baseline, when baseline accuracy requires to be 3mm ± 0.01ppmD when (mm represents that millimeter, ppm represent micrometre), starting at the point coordinate error should be less than 0.2m (rice).
From the above, the moving relative position of three two frequently differences location in the moving location accurately being determined carrier of GNSS, need accurate reference station coordinates, because the position of two motion carriers is all in motion, therefore adopt three frequency Static Precise Point Positioning to determine the initial position of motion carrier, can obtain inferior decimeter grade bearing accuracy, and with this coordinate as the reference station, carry out the difference location and find the solution dynamic baseline, and the integer ambiguity constraint Static Precise Point Positioning of utilizing difference to resolve the baseline component of acquisition and solve, the precision and the convergence time thereof of raising Static Precise Point Positioning.
In the GNSS navigator fix, often utilize low noise, the light current absciss layer, long wavelength's carrier phase combination observation value improves the success ratio that blur level is calculated, and three provide new method for the GNSS ambiguity resolution frequently.In order to improve the success ratio of counting yield and ambiguity resolution, can adopt existing Geometry-based TCAR method to consider the geometry intensity information of satellite constellation during concrete enforcement, adopt existing ambiguity resolution technology Lambda method to find the solution super wide lane ambiguity earlier, find the solution wide lane ambiguity again, find the solution the L1 values of ambiguity at last, that is: the first step, utilize super wide lane combination observation value and pseudorange observed reading to form observation equation, obtain the real solution of super wide lane combination observation value, according to the success ratio thresholding, adopt the super wide lane ambiguity of Lambda method partial fixing, second step is similar with the first step, super wide lane combination observation value and fixing integer ambiguity value and wide lane combination observation value substitution observation equation, find the solution the real solution and the variance and covariance battle array thereof of wide lane ambiguity, adopt the wide lane of Lambda method partial fixing integer ambiguity.The 3rd step was adopted the integer ambiguity of the direct fixed L 1 of similar methods.In the linear combination observed reading that adopts.Super Kuan Xiang and wide lane all belong to the Kuan Xiang territory.The GNSS three of embodiment Ambiguity Solution Methods frequently may further comprise the steps as shown in Figure 6:
Steps A is calculated super wide lane combination observation value WL2 (combination coefficient is (0,1 ,-1)) blur level real solution and variance and covariance thereof, LAMBDA blur level partial fixing, WL2 ambiguity resolution result (comprising fixed value and real number value); Judge whether then all to fix, be then to enter step B, otherwise enter step D;
Step B calculates wide lane combination observation value WL1 (combination coefficient is (1,0 ,-1)) blur level real solution and variance and covariance thereof, LAMBDA blur level partial fixing, WL1 ambiguity resolution result (comprising fixed value and real number value); Judge whether then all to fix, be then to enter step C, otherwise enter step D;
Step C calculates L1 blur level real solution and variance and covariance thereof, LAMBDA blur level partial fixing, and L1 ambiguity resolution result (comprising fixed value and real number value) enters step D;
Step D utilizes the blur level that obtains to calculate the baseline component, output result and corresponding integer ambiguity thereof.
For the GNSS ambiguity resolution, when the blur level number more for a long time, can not correctly fix all blur leveles single epoch, the integer ambiguity according to certain success ratio fixed part can be fixed can improve the precision that baseline resolves.Its specific implementation method is as follows: at first the variance and covariance battle array of real number blur level (being the real solution of blur level) is carried out decorrelation, obtain the integer transition matrix, utilize the integer transition matrix that the real number blur level is carried out integer transform, real number blur level and variance and covariance battle array thereof behind the acquisition integer transform, calculate the fixing success ratio of blur level according to the variance and covariance battle array behind the integer transform, get one group of blur level of its success ratio maximum in domain value range, utilize the search of Lambda method, obtain the optimum solution and the suboptimal solution (when differing greatly) of blur level, utilize existing assumed statistical inspection theory, the fixing correctness of test blur level.If upcheck, think that then the integer ambiguity that optimum solution obtains is exactly the blur level of asking, blur level is fixing correct, adopts the integer ambiguity value to find the solution baseline, if the real number values of ambiguity is then adopted in the check failure.The part blur level fixing means of embodiment is summed up its process and is as shown in Figure 7: for real number blur level and variance and covariance battle array thereof, with real number blur level behind the Lambda integer transform and variance and covariance battle array thereof, by the ordering of variance size; Select to treat fixedly blur level group and variance and covariance battle array thereof according to given success ratio thresholding (suggestion value 96%), utilize the search blur level to obtain to treat the fixedly optimal value and time figure of merit of blur level group, upcheck then as static solution, the check failure is then as real solution.

Claims (5)

1. a GNSS three is moving frequently to motion positioning method, it is characterized in that:
In initial epoch, carry out following steps,
Step 1 as the mobile observation station, gets two motion carriers three frequency GNSS data to two original observation stations in mobile observation station and carries out pre-service;
Step 2 is according to the pretreated three frequency GNSS data of step 1 gained, by the coordinate at two mobile observation stations of three Static Precise Point Positioning acquisitions frequently;
Step 3, establishing arbitrary in two mobile observation stations is the reference station of three two frequently differences location, according to the coordinate at two mobile observation stations of step 2 gained, finds the solution the baseline component at two mobile observation stations and the integer ambiguity of two poor carrier phases by three two frequently differences location;
In follow-up epoch, two original observation stations in mobile observation station are got three frequency GNSS data carry out pre-service, according to pretreated three frequency GNSS data, by the coordinate at two mobile observation stations of three Static Precise Point Positioning acquisitions frequently, and with the constraint condition of the step 3 gained solving result in initial epoch as three frequency Static Precise Point Positioning.
2. GNSS three as claimed in claim 1 is moving frequently to motion positioning method, it is characterized in that: described to two original observation stations in mobile observation station get three frequently the GNSS data carry out pre-service, comprise the uncorrelated combination observation value that adopts three three not have ionosphere, long wavelength and low noise carrier wave frequently, two original observation stations in mobile observation station are got three frequency GNSS data carry out the cycle slip detection and repair.
3. GNSS as claimed in claim 1 is moving to motion positioning method, it is characterized in that: when carrying out three frequency Static Precise Point Positioning, adopt three the three uncorrelated estimation of combination observation values troposphere, receiver location and the motion carrier clock correction parameters of not having ionosphere, long wavelength and low noise carrier wave frequently.
4. moving frequently to motion positioning method as claim 1 or 2 or 3 described GNSS three, it is characterized in that: during find the solution the two frequently differences of step 3 described three location, finding the solution of the integer ambiguity of described two difference carrier phases adopts Geometry-Based TCAR method to realize, according to success ratio thresholding fixed part blur level, adopt assumed statistical inspection to judge the reliability of fixed part blur level then.
5. GNSS three as claimed in claim 4 is moving frequently to motion positioning method, it is characterized in that: describedly be according to success ratio thresholding fixed part blur level implementation, obtain the real solution and the variance and covariance battle array thereof of blur level by two poor location models, calculate the fixing success ratio of blur level according to the variance and covariance battle array behind the integer transform, get one group of blur level of its success ratio maximum in default success ratio domain value range.
CN 200910273365 2009-12-23 2009-12-23 Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method Active CN101710179B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200910273365 CN101710179B (en) 2009-12-23 2009-12-23 Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200910273365 CN101710179B (en) 2009-12-23 2009-12-23 Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method

Publications (2)

Publication Number Publication Date
CN101710179A true CN101710179A (en) 2010-05-19
CN101710179B CN101710179B (en) 2012-12-05

Family

ID=42402976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200910273365 Active CN101710179B (en) 2009-12-23 2009-12-23 Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method

Country Status (1)

Country Link
CN (1) CN101710179B (en)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102288978A (en) * 2011-07-20 2011-12-21 东南大学 Continuous operational reference system (CORS) base station cycle slip detection and recovering method
CN102353969A (en) * 2011-09-02 2012-02-15 东南大学 Method for estimating phase deviation in precise single-point positioning technology
CN102636800A (en) * 2012-04-28 2012-08-15 北京航空航天大学 Decorrelation method of global navigation satellite system (GNSS) integer ambiguity based on finishing pretreatment
CN102650692A (en) * 2011-02-25 2012-08-29 中国人民解放军61081部队 Method for detecting and repairing cycle slip by utilizing BeiDou three-frequency observed quantity
CN103376454A (en) * 2012-04-12 2013-10-30 天宝导航有限公司 Advanced global navigation satellite systems (gnss) positioning using precise satellite information
CN103487821A (en) * 2012-06-12 2014-01-01 广州中海达卫星导航技术股份有限公司 Baseline vector solution method
CN103675874A (en) * 2013-12-20 2014-03-26 北京遥测技术研究所 Method for determining whole-cycle ambiguity of three-frequency carrier phase of BeiDou navigation system
CN104155668A (en) * 2014-08-28 2014-11-19 哈尔滨工程大学 Underway locating method for long-base-line carrier phase difference
CN104483691A (en) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 GNSS combined precise single-point positioning method
CN104502943A (en) * 2014-06-03 2015-04-08 西北工业大学 Indoor pseudo-satellite differential relative positioning method
CN104536027A (en) * 2015-01-28 2015-04-22 中国人民解放军国防科学技术大学 Real-time bayou precision relative positioning method
CN104536026A (en) * 2015-01-08 2015-04-22 中国航空无线电电子研究所 Dynamic-to-dynamic real-time measurement system
CN105093251A (en) * 2015-07-13 2015-11-25 中国电子科技集团公司第十研究所 High-precision relative positioning method of GNSS receiver in static state
CN105242293A (en) * 2014-07-08 2016-01-13 成都国星通信有限公司 High-precision centimeter-level positioning method of global navigation satellite system
CN106772478A (en) * 2016-11-11 2017-05-31 哈尔滨工程大学 The localization method of difference constraint between a kind of star based on epoch
CN108267135A (en) * 2017-12-25 2018-07-10 中铁第四勘察设计院集团有限公司 For the accurate positioning method and system of track automatic measurement vehicle
CN108508470A (en) * 2018-05-16 2018-09-07 武汉大学 Towards the instantaneous decimeter grade navigation locating method that the whole world is seamless
CN108548479A (en) * 2018-04-16 2018-09-18 武汉大学 Bridge based on GNSS monitors fast initializing method in real time
CN109507703A (en) * 2018-12-17 2019-03-22 成都国星通信有限公司 A kind of GNSS carrier phase calculation method
CN110568464A (en) * 2019-06-19 2019-12-13 航天信息股份有限公司 BDS/GNSS (broadband navigation satellite system/global navigation satellite system) multi-mode chip-based precision positioning method and device
CN110687558A (en) * 2019-12-11 2020-01-14 深圳微品致远信息科技有限公司 Independent double-baseline calculation verification method based on Beidou dynamic positioning
CN112051598A (en) * 2020-06-24 2020-12-08 中铁第四勘察设计院集团有限公司 Vehicle-mounted GNSS/INS integrated navigation method based on double correction
CN112731502A (en) * 2020-11-23 2021-04-30 南京航空航天大学 Unmanned aerial vehicle aerial refueling inertia-assisted Beidou three-frequency precise relative navigation method
CN114675314A (en) * 2022-05-30 2022-06-28 深圳市智联时空科技有限公司 Re-convergence precise point positioning method
CN115902981A (en) * 2022-11-14 2023-04-04 中南大学 Train positioning optimization method and system and rail transit vehicle

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2849209B1 (en) * 2002-12-19 2007-04-06 Agence Spatiale Europeenne METHOD AND SYSTEM FOR REAL-TIME NAVIGATION USING THREE-CARRIER SATELLITE-TYPE RADIO ELECTRIC SIGNALS AND IONOSPHERIC CORRECTIONS
US6934632B2 (en) * 2003-10-08 2005-08-23 Navcom Technology, Inc. Method for using three GPS frequencies to resolve carrier-phase integer ambiguities

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102650692A (en) * 2011-02-25 2012-08-29 中国人民解放军61081部队 Method for detecting and repairing cycle slip by utilizing BeiDou three-frequency observed quantity
CN102650692B (en) * 2011-02-25 2014-07-23 中国人民解放军61081部队 Method for detecting and repairing cycle slip by utilizing BeiDou three-frequency observed quantity
CN102288978B (en) * 2011-07-20 2013-09-18 东南大学 Continuous operational reference system (CORS) base station cycle slip detection and recovering method
CN102288978A (en) * 2011-07-20 2011-12-21 东南大学 Continuous operational reference system (CORS) base station cycle slip detection and recovering method
CN102353969A (en) * 2011-09-02 2012-02-15 东南大学 Method for estimating phase deviation in precise single-point positioning technology
CN103376454A (en) * 2012-04-12 2013-10-30 天宝导航有限公司 Advanced global navigation satellite systems (gnss) positioning using precise satellite information
CN103376454B (en) * 2012-04-12 2017-03-01 天宝导航有限公司 Advanced GLONASS using precise satellite Information locating(GNSS)
CN102636800A (en) * 2012-04-28 2012-08-15 北京航空航天大学 Decorrelation method of global navigation satellite system (GNSS) integer ambiguity based on finishing pretreatment
CN103487821A (en) * 2012-06-12 2014-01-01 广州中海达卫星导航技术股份有限公司 Baseline vector solution method
CN103487821B (en) * 2012-06-12 2015-10-28 广州中海达卫星导航技术股份有限公司 A kind of baseline vector solution method
CN103675874B (en) * 2013-12-20 2015-11-25 北京遥测技术研究所 A kind of triones navigation system three ambiguity of carrier phase defining method frequently
CN103675874A (en) * 2013-12-20 2014-03-26 北京遥测技术研究所 Method for determining whole-cycle ambiguity of three-frequency carrier phase of BeiDou navigation system
CN104502943B (en) * 2014-06-03 2017-07-04 西北工业大学 A kind of indoor pseudo satellite, pseudolite difference relative positioning method
CN104502943A (en) * 2014-06-03 2015-04-08 西北工业大学 Indoor pseudo-satellite differential relative positioning method
CN105242293B (en) * 2014-07-08 2017-10-17 成都国星通信有限公司 A kind of high-precision centimeter-level positioning method of GLONASS
CN105242293A (en) * 2014-07-08 2016-01-13 成都国星通信有限公司 High-precision centimeter-level positioning method of global navigation satellite system
CN104155668A (en) * 2014-08-28 2014-11-19 哈尔滨工程大学 Underway locating method for long-base-line carrier phase difference
CN104483691B (en) * 2015-01-05 2017-05-31 中国科学院重庆绿色智能技术研究院 A kind of GNSS combines accurate one-point positioning method
CN104483691A (en) * 2015-01-05 2015-04-01 中国科学院重庆绿色智能技术研究院 GNSS combined precise single-point positioning method
CN104536026A (en) * 2015-01-08 2015-04-22 中国航空无线电电子研究所 Dynamic-to-dynamic real-time measurement system
CN104536027B (en) * 2015-01-28 2015-09-23 中国人民解放军国防科学技术大学 A kind of method of real-time Big Dipper Precise Relative Positioning
CN104536027A (en) * 2015-01-28 2015-04-22 中国人民解放军国防科学技术大学 Real-time bayou precision relative positioning method
CN105093251A (en) * 2015-07-13 2015-11-25 中国电子科技集团公司第十研究所 High-precision relative positioning method of GNSS receiver in static state
CN106772478A (en) * 2016-11-11 2017-05-31 哈尔滨工程大学 The localization method of difference constraint between a kind of star based on epoch
CN108267135A (en) * 2017-12-25 2018-07-10 中铁第四勘察设计院集团有限公司 For the accurate positioning method and system of track automatic measurement vehicle
CN108548479B (en) * 2018-04-16 2019-09-10 武汉大学 Bridge real-time monitoring fast initializing method based on GNSS
CN108548479A (en) * 2018-04-16 2018-09-18 武汉大学 Bridge based on GNSS monitors fast initializing method in real time
CN108508470A (en) * 2018-05-16 2018-09-07 武汉大学 Towards the instantaneous decimeter grade navigation locating method that the whole world is seamless
CN109507703A (en) * 2018-12-17 2019-03-22 成都国星通信有限公司 A kind of GNSS carrier phase calculation method
CN110568464B (en) * 2019-06-19 2023-10-10 航天信息股份有限公司 BDS/GNSS multimode chip-based precise positioning method and BDS/GNSS multimode chip-based precise positioning device
CN110568464A (en) * 2019-06-19 2019-12-13 航天信息股份有限公司 BDS/GNSS (broadband navigation satellite system/global navigation satellite system) multi-mode chip-based precision positioning method and device
CN110687558A (en) * 2019-12-11 2020-01-14 深圳微品致远信息科技有限公司 Independent double-baseline calculation verification method based on Beidou dynamic positioning
CN110687558B (en) * 2019-12-11 2020-05-12 深圳微品致远信息科技有限公司 Independent double-baseline calculation verification method based on Beidou dynamic positioning
CN112051598A (en) * 2020-06-24 2020-12-08 中铁第四勘察设计院集团有限公司 Vehicle-mounted GNSS/INS integrated navigation method based on double correction
CN112051598B (en) * 2020-06-24 2023-09-29 中铁第四勘察设计院集团有限公司 Dual-correction-based vehicle-mounted GNSS/INS integrated navigation method
CN112731502A (en) * 2020-11-23 2021-04-30 南京航空航天大学 Unmanned aerial vehicle aerial refueling inertia-assisted Beidou three-frequency precise relative navigation method
CN114675314A (en) * 2022-05-30 2022-06-28 深圳市智联时空科技有限公司 Re-convergence precise point positioning method
CN114675314B (en) * 2022-05-30 2022-08-26 深圳市智联时空科技有限公司 Re-convergence precise point positioning method
CN115902981A (en) * 2022-11-14 2023-04-04 中南大学 Train positioning optimization method and system and rail transit vehicle
CN115902981B (en) * 2022-11-14 2024-01-30 中南大学 Train positioning optimization method and system and rail transit vehicle

Also Published As

Publication number Publication date
CN101710179B (en) 2012-12-05

Similar Documents

Publication Publication Date Title
CN101710179B (en) Global navigation satellite system (GNSS) triple-frequency motion-to-motion positioning method
CN109709591B (en) GNSS high-precision positioning method for intelligent terminal
CN110531392A (en) A kind of high-precision locating method and system based on PPP algorithm
US7292185B2 (en) Attitude determination exploiting geometry constraints
US8035552B2 (en) Distance dependant error mitigation in real-time kinematic (RTK) positioning
AU2008260578B2 (en) Distance dependant error mitigation in real-time kinematic (RTK) positioning
CN101770033B (en) Fixing method of ambiguity network between CORS and system station
US8120527B2 (en) Satellite differential positioning receiver using multiple base-rover antennas
CN108802782B (en) Inertial navigation assisted Beidou three-frequency carrier phase integer ambiguity solving method
Li et al. Modeling and quality control for reliable precise point positioning integer ambiguity resolution with GNSS modernization
CN103837879B (en) The method of hi-Fix is realized based on the civilian combination carrier phase observation of dipper system
CN105223598B (en) A kind of GNSS ambiguity of carrier phase simple epoch solution method
CN104483690B (en) A kind of frequency Static Precise Point Positioning fuzziness fixing means of GNSS tri-
CN102016624B (en) Inter-moving body interferometric positioning system, device and method thereof
US10739471B2 (en) GNSS receiver with a capability to resolve ambiguities using an uncombined formulation
Su et al. Triple-frequency carrier phase precise time and frequency transfer models for BDS-3
CN110058282A (en) A kind of PPP high-precision locating method based on double frequency GNSS smart phone
Morales et al. GNSS vertical dilution of precision reduction using terrestrial signals of opportunity
Mosavi et al. Position estimation in single-frequency GPS receivers using Kalman filter with pseudo-range and carrier phase measurements
Zhang et al. A novel GNSS based V2V cooperative localization to exclude multipath effect using consistency checks
CN105510945A (en) PPP positioning method applied to satellite navigation landing outfield detection
CN111983650B (en) High-precision time transfer method based on GNSS
Seepersad Reduction of initial convergence period in GPS PPP data processing
Tolman et al. Absolute precise kinematic positioning with GPS and GLONASS
CN103543454A (en) Satellite orbit determination system inserted in mobile communication network

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180726

Address after: 264200 Shandong Weihai Nanhai New District Binhai Road North, Longhai Road East

Patentee after: Weihai Wuzhou Satellite Navigation Technology Co., Ltd.

Address before: 430072 Hubei city of Wuhan province Wuchang Luojiashan

Patentee before: Wuhan University

TR01 Transfer of patent right