CN102565813A - Method and device for performing pseudorange observation estimation by carrier smoothing - Google Patents
Method and device for performing pseudorange observation estimation by carrier smoothing Download PDFInfo
- Publication number
- CN102565813A CN102565813A CN2010106197744A CN201010619774A CN102565813A CN 102565813 A CN102565813 A CN 102565813A CN 2010106197744 A CN2010106197744 A CN 2010106197744A CN 201010619774 A CN201010619774 A CN 201010619774A CN 102565813 A CN102565813 A CN 102565813A
- Authority
- CN
- China
- Prior art keywords
- tec
- pseudorange
- frequency
- eff
- smoothing
- 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
Links
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a method and a device for performing pseudorange observation estimation by carrier smoothing. The method comprises the following steps: respectively getting a pseudorange and carrier phase observation quantity of each of n frequency points, wherein n is an integer which is not less than 2; simultaneously performing carrier smoothing pseudorange operation on all the frequency points to get the pseudoranges after smoothing; enabling the pseudoranges after smoothing of the multiple frequency points to constitute a measurement equation; and performing value estimation according to the measurement equation and taking the obtained estimated value of r as the estimated value of pseudorange observation quantity. According to the method and the device, on the basis of greatly reducing multi-path and measuring noise, the influence of an ionized layer is eliminated, and the pseudorange estimation precision is further significantly improved.
Description
Technical field
The present invention relates to navigation field, be specifically related to a kind ofly carry out pseudo range observed quantity estimation approach and device through carrier smoothing.
Background technology
The speed that ripple is propagated in a kind of medium depends on the refractive index of this medium.Refractive index n be defined as ripple free space velocity of propagation c with at the ratio of medium medium velocity v.For the electromagnetic wave of GNSS (GlobalNavigation Satellite System, GLONASS) emission, with light velocity propagation, the velocity of propagation in ionosphere and troposphere is not equal to the light velocity in free space.
Refractive index: n=c/v (1)
If the refractive index of ripple in a kind of medium is unfixing, but the function of frequency f, then this medium is chromatic dispersion for this ripple.For the electromagnetic wave of L-band, ionosphere is dispersive medium, and the troposphere is a nondispersive medium.In dispersive medium, velocity of propagation (the phase velocity v of signal carrier phase
p) velocity of propagation (the group velocity v of the signal message that carried with ripple
g) unequal.
n
p=c/v
p v
g=c/v
g (2)
Group index n
gWith refractive index n mutually
pRelation (can be referring to Hofmann-Wellenhof, B., H.Lichtenegger, and J.Collins, GPS Theory and Practice, New York:Springer-Verlag, 1993) as follows:
The expression formula of ionosphere phase refractive index is:
Can know that by relational expression 3 group index is:
Obtain after ignoring higher order term:
n
eBe the ionospheric refraction rate.
Be respectively apart from l and measured value (pseudorange) S thereof between satellite SV and receiver user User:
Anaclastic measuring error is:
Therefore, reflecting the error of on carrier phase and pseudo range observed quantity, introducing is respectively:
Electron density along the electromagnetic wave propagation path is called total electron content, is defined as:
Formula 9 can be rewritten as:
Can draw to draw a conclusion based on above quantitative test: the error that ionospheric refraction is introduced carrier phase measurement is leading, then is to lag behind to pseudorange, but its absolute value equates, with square being inversely proportional to of frequency, is directly proportional with TEC.Theoretically, the scattering effect of utilizing frequency to cause possibly estimate TEC, utilizes TEC and frequency can obtain the range error at each frequency again.
Can know that by the analysis of preceding text ionosphere time-delay is the function of frequency, theoretically, the GPS-L1 of dual-frequency receiver output and L2 frequency pseudo range observed quantity carried out calculus of differences can eliminate its influence fully.The expression formula of pseudorange P is following:
P=ρ+dρ+c(dt-dT)+ΔS
ion,g+ΔS
trop,g+ε
P (12)
For the pseudo range observed quantity of L1 and two frequencies of L2, very apart from ρ, satellite orbital error d ρ, satellite clock correction dt, receiver clock correction dT and troposphere time-delay Δ S
Trop, gEquate relation below the ionosphere time-delay is satisfied:
f
L1 2ΔS
ion,g,L1=f
L2 2ΔS
ion,g,L2 (13)
Measure noise and multipath ε
PBetween two frequencies, do not have definite relation, therefore be modeled as noise jointly.According to relational expression 13, can obtain ionosphere time-delay estimated value
Can know according to aforementioned analysis
Finally
Clearly, the accuracy of ionosphere time-delay estimated value depends on and the size of noise under the very big situation of multipath error, can't guarantee the validity of estimated value.In addition, if the too small f that causes of the frequency difference between frequency
L2 2-f
L1 2Too small, noise is exaggerated, and also can reduce the accuracy of estimated value.
Strictly speaking, the pseudorange of different frequent points has only receiver clock correction to equate, because ionospheric effect of dispersion, if the electromagnetic wave of two frequencies arrives receiving antenna at synchronization, then together, its travel path also has minute differences to its x time scarcely.For no other reason than that the time-delay maximum that effect of dispersion causes hundreds of rice only; Then between two frequencies x time at interval below the microsecond level; In the so short mistiming; The motion of satellite and receiver user can be ignored, and the change of satellite clock correction is also insignificant, and the correlativity in time and space makes that the difference of two frequency troposphere time-delays and ionized layer TEC is equally very small.
For the higher GNSS user of positioning accuracy request, ionospheric propagation time-delay and multipath effect are pseudo range measurements and based on two main sources of pseudorange positioning error.At the beginning of design, the Service Signal of GPS just provides two frequencies, hopes to utilize the electromagnetic effect of dispersion of L-band, and estimation is also proofreaied and correct the pseudorange delay time error that is caused by ionosphere.Excessive, possibly cause the time-delay of double frequency pseudorange ionosphere to estimate failure but if multipath effect acts on the error that observed quantity causes.Multipath error is closely related with the surrounding environment of receiving antenna, although industry has been developed multiple sign indicating number ring (DLL) multipath inhibition technology, still is difficult to thoroughly solve the pseudorange multipath error problem in the harsh and unforgiving environments such as city, mountain region.Commercial receiver thereby generally adopted carrier smoothing pseudorange technology utilizes carrier wave observed quantity multipath and measures noise much smaller than pseudo range observed quantity, reduces multipath and noise, satisfies the requirement of pseudorange locating accuracy.
Because ionosphere is time-delay to the influence of pseudorange, be leading to the influence of carrier phase, the carrier smoothing pseudorange can be introduced a system deviation (ionosphericdivergence) again on time-delay basis, the original ionosphere of pseudorange.And this system deviation can accumulate with the continuous service of smoothing algorithm, finally to cause the size of level and smooth back pseudorange error to surpass original/and the error of smoothing pseudo range not.Therefore, smoothing algorithm must reset after lasting computing a period of time, restarted level and smooth.The length of duration depends on the estimation to ionosphere systematic error parameter, if ionosphere is active, the ionosphere system deviation can cause level and smooth back pseudorange error bigger, and frequent reset operation makes multipath and measurement Noise Suppression effect descend.
Summary of the invention
The technical matters that the present invention will solve provides a kind ofly carries out pseudo range observed quantity estimation approach and device through carrier smoothing, the multipath that greatly reduces with measure on the effective noise floor, rejected ionospheric influence, thereby significantly improved the estimated accuracy of pseudorange.
In order to address the above problem, the invention provides and a kind ofly carry out the pseudo range observed quantity estimation approach through carrier smoothing, comprising:
Obtain the pseudorange and the carrier phase observed quantity of n frequency respectively; N is the integer more than or equal to 2;
Carry out carrier smoothing pseudorange operation simultaneously at each frequency respectively, obtain the pseudorange
after level and smooth;
Pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is since 1, and the integer till the n comprises 1 and n; d
1And d
2Be constant;
According to above-mentioned measurement equation to TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
Preferably,
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
Preferably, d
1Span be 40.3 ± 10%, d
2Span be 1 ± 10%.
Preferably, d
1Be 40.3, d
2Be 1.
Preferably, said according to measurement equation to TEC
EffCarry out in the step of valuation with r, adopt least-squares estimation, weighted least-squares estimation or Kalman filtering to carry out valuation.
The present invention also provides a kind of and has carried out the device that pseudo range observed quantity is estimated through carrier smoothing, comprising:
Measurement module is used for obtaining respectively the pseudorange and the carrier phase observed quantity of n frequency; N is the integer more than or equal to 2;
Level and smooth module; Be used for respectively carrying out the operation of carrier smoothing pseudorange simultaneously, obtain the pseudorange
after level and smooth at each frequency
Constructing module is used for the pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is since 1, and the integer till the n comprises 1 and n; d
1And d
2Be constant;
Estimator module, according to above-mentioned measurement equation to TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
Preferably,
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
Preferably, d
1Span be 40.3 ± 10%, d
2Span be 1 ± 10%.
Preferably, d
1Be 40.3, d
2Be 1.
Preferably, said estimator module adopts least-squares estimation, weighted least-squares estimation or Kalman filtering to carry out valuation.
Technical scheme of the present invention can significantly reduce the multipath effect influence, significantly reduces and measures noise intensity, effectively rejects the influence of ionosphere to pseudo-range measurements; Owing to can carry out whole valuation correction to the ionosphere system deviation that ionospheric error itself is introduced when level and smooth, so the selection of time constant need not be considered ionospheric situation again, duration that also need not limit algorithm.
Description of drawings
Fig. 1 is the schematic block diagram of the pseudo range observed quantity estimation unit of embodiment two.
Embodiment
To combine accompanying drawing and embodiment that technical scheme of the present invention is explained in more detail below.
Need to prove that if do not conflict, each characteristic among the embodiment of the invention and the embodiment can mutually combine, all within protection scope of the present invention.In addition; Can in computer system, carry out in the step shown in the process flow diagram of accompanying drawing such as a set of computer-executable instructions, and, though logical order has been shown in process flow diagram; But in some cases, can carry out step shown or that describe with the order that is different from here.
At first, the inventor is following to the analysis of carrier smoothing pseudorange algorithm:
Multipath in the pseudorange disturbs and measures noise and is significantly higher than the carrier phase observed quantity, but the latter has comprised a fixing integral circumference ambiguity.Bearing accuracy based on the raw pseudo range observed quantity is lower, though carrier phase difference location can provide the precision of centimetre-sized, finding the solution integral circumference ambiguity with the observed quantity that provides by means of the difference base station is prerequisite.Utilize accurate relatively carrier phase observed quantity smoothing pseudo range can under the condition that need not find the solution integral circumference ambiguity, improve the receiver-autonomous bearing accuracy of GNSS.The expression formula of carrier phase:
Φ=ρ+dρ+c(dt-dT)+ΔS
ion,p+ΔS
trop,p+N
Φ+M
Φ+n
Φ (17)
Comparison expression (12) and (17), carrier phase are very apart from ρ, satellite orbital error d ρ, satellite clock correction dt, receiver clock correction dT and troposphere time-delay Δ S
Trop, pIdentical with in the expression formula of pseudorange, visible from formula (11), ionosphere time-delay equal and opposite in direction, opposite in sign.N
ΦThe expression integral circumference ambiguity, the size of carrier phase multipath error and measurement noise is much smaller than pseudo range observed quantity.Carrier smoothing pseudorange algorithm comes down to a regressive filter, progressively increases the weight of carrier phase and reduces the weight of pseudorange and keep up to resetting.
K is the discrete time constant.Under the original state:
When the m increase and above K, weight no longer changes:
Therefore, thereafter algorithm mathematics model can be described by constant coefficient recursion equation:
Error analysis to carrier smoothing pseudorange algorithm is following:
Because carrier phase is very apart from ρ, satellite orbital error d ρ, satellite clock correction dt, receiver clock correction dT and troposphere time-delay Δ S
Trop, pIdentical with in the expression formula of pseudorange, and also equal basically between different frequencies, can rewrite the expression formula of pseudorange and carrier phase:
P
i=r
i+I
i+M
P_i+n
P_i (23)
Φ
i=r
i-I
i+N
Φ+M
Φ_i+n
Φ_i
Wherein
r
i=ρ
i+dρ
i+c(dt
i-dT
i)+ΔS
trop_i (24)
I
i=ΔS
ion,g_i ΔS
trop_i=ΔS
trop,g_i=ΔS
trop,p_i
Definition is dispersed the error that (ionospheric divergence), multipath and measurement noise cause by ionosphere sign indicating number-carrier wave:
The Δ t of definition smoothingtime constant τ=(K-1), Δ t is the minimum SI of smoothing pseudo range.Ignore carrier wave multipath and carrier wave and measure noise, under Δ t → 0 condition, can represent the Algorithm Error of following formula definition by following integral equation:
Multipath error depends on the geometry site between reflecting body and antenna fully; Its Changing Pattern is similar to sine function; Cycle is greatly about a minute level (can be referring to M.Elizabeth Cannon, Lecture Notes ofENGO 561 Satellite Positioning, Winter 2005):
M
P=A
Pcos(ω
Pt+α
P) (27)
Multipath signal strength is A
P, the cycle is ω
PThe inverse of/2 π.Pseudorange multipath effect steady-state error behind carrier smoothing is:
The level and smooth back multipath error cycle that is all sine function is constant, and intensity weakens.Clearly, the smoothingtime constant is big more, and effect is good more.
Suppose that the pseudo range measurement noise is a white noise, time correlation function is:
Wherein δ (t) is impulse function/Dirac function, σ
P 2It is the pseudorange noise power.After level and smooth, noise power is:
σ
SP 2=σ
P 2/2τ (30)
Thus it is clear that,, also hope to adopt bigger time constant from suppressing the angle of noise.
Suppose that ionosphere time-delay Changing Pattern is as follows:
I=a+bt+ct
2 (31)
Since pseudorange and carrier phase ionosphere time-delay opposite in sign, the steady-state system deviation that smoothing algorithm is introduced:
ε
I=4cτt (32)
The once existence of a 4c τ t about the algorithm duration means the constantly accumulation of ionosphere system deviation, and the speed of accumulating rate is directly proportional with time constant.In order to reduce the ionosphere system deviation; The algorithm time constant is dwindled in requirement; And the algorithm continuous working period can not be long, otherwise the benefit that complete cancellation algorithms inhibition multipath of system deviation meeting and measurement noise bring, even the pseudorange total error after causing smoothly is greater than raw pseudo range.The selection of time constant and duration is all depended on (can be referring to Euiho Kim to the estimation of ionosphere and multipath effect parameter; Todd Walter; And J.D.Powell, " Adaptive CarrierSmoothing Using Code and Carrier Divergence ", ION NTM 2007; San Diego, CA).
Embodiment one, a kind ofly carries out the pseudo range observed quantity estimation approach through carrier smoothing, comprising:
Obtain the pseudorange and the carrier phase observed quantity of n frequency respectively; N is the integer more than or equal to 2;
Carry out carrier smoothing pseudorange operation simultaneously at each frequency respectively, obtain the pseudorange
after level and smooth;
Pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;(33)
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is 1,2 ..., n, just be since 1, the integer till the n comprises 1 and n; d
1And d
2Be constant;
According to above-mentioned measurement equation to TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
In the present embodiment, the smoothing pseudo range equation structure measurement equation through a plurality of frequencies obtains the pseudo range observed quantity estimated value with the mode of valuation, can reject the influence that bring in ionosphere, improves result's precision.
Wherein, d
1And d
2Can confirm according to derivation or simulation result.
In the present embodiment, d
2Span can but be not limited to 1 ± 10%, optimum is 1.
Obviously, the Z in the formula (33) is the measurement of state X, and V is an error in measurement.Treat estimator TEC for two
EffAnd r, can adopt multiple optimum estimation theory, criterion and method to estimate.
Further analyze from the principle angle below.
Select the reasonable time constant can effectively suppress to measure noise and multipath; Under the negative exponential function effect, initial deviation item ε
0+Be tending towards 0 gradually with the algorithm accumulated time, thereby can in steady-state analysis, ignore.Be illustrated in the smoothing process by the pseudo range measurement noise and the multipath that fully suppress with
; And the noise like item of the measurement noise of carrier phase and multipath interference composition; Then according to formula (25), the pseudorange after level and smooth can be expressed as:
Bringing formula (11) into following formula can get:
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
In the present embodiment, can use TEC
EffAs the ionospheric equivalent TEC of carrier smoothing pseudorange:
Can find out TEC
Eff' be the equivalent TEC that has taken all factors into consideration behind the population effect of the time-delay of raw pseudo range ionosphere and ionosphere system deviation two error terms relevant of carrier smoothing pseudorange algorithm introducing itself with ionosphere; Present embodiment through type (36) can be merged into a unknown quantity to the time-delay of raw pseudo range ionosphere and two unknown quantitys of additional ionospheric error that carrier smoothing pseudorange algorithm is introduced; Add r two unknown quantitys altogether, only need the smoothing pseudo range equation of two frequencies just can carry out valuation; Usually, the GPS receiver can provide three frequencies, and therefore the redundant equation of a frequency will be arranged, and can make that the valuation result is more accurate.
In a kind of embodiment of present embodiment, can be directly with the TEC in the formula (37)
Eff' as the TEC in the formula (34)
Eff, at this moment, d
1Span can but be not limited to 40.3 ± 10%, optimum is 40.3.
Theoretically, TEC
EffWith r and frequency-independent, if multifrequency receiver carries out the operation of carrier smoothing pseudorange synchronously for each frequency (in time), then aforementioned two variablees of synchronization equate in the value of each frequency.Two unknown quantitys need only plural equation/frequency can be found the solution and obtained about TEC
EffWith the valuation of r, wherein r is carried out PVT resolve and just can accomplish receiver-autonomous location and test the speed.Because present embodiment can carry out whole valuation correction to the ionosphere system deviation that ionospheric error itself and carrier smoothing pseudorange algorithm are introduced, so the selection of time constant need not be considered ionospheric situation again, duration that also need not limit algorithm.
Behind formula (37) substitution (36) substitution following formula, the pseudorange after level and smooth can be expressed as:
Can find out and since in the valuation of r multipath with measure noise and received effective inhibition of carrier smoothing, ionospheric influence is rejected through above-mentioned Least Squares Estimating again.Therefore, the valuation of r is accurate more than original pseudorange and single-frequency carrier smoothing pseudorange, and the PVT that carries out on this basis resolves and can access higher receiver-autonomous bearing accuracy.
Suppose that receiver generates the pseudorange and the carrier phase observed quantity of n frequency, and carry out the operation of carrier smoothing pseudorange synchronously, then can obtain following system of equations at each frequency:
...
The matrix form of following formula is formula (33): Z=HX+V; Z, H, X and V see formula (34).
In another embodiment of present embodiment, also can use TEC
Eff"=40.3 * TEC
Eff' as the TEC in the formula (34)
Eff, at this moment,
d
1Span can but be not limited to 1 ± 10%, optimum is 1.
During practical application, can use d
3* TEC
Eff' be used as as the TEC in the formula (34)
Eff, d
3Be constant; At this moment, corresponding adjustment d
1Just can reach the same valuation effect, optimum is to make d
1And d
3Amassing is 40.3; Further, because the TEC in the formula (34)
EffBe to treat estimated value, therefore also can and TEC
EffAs long as ' disproportional relation is adjustment d
1Get final product.
In the present embodiment, said according to above-mentioned measurement equation to TEC
EffCarry out in the step of valuation with r, can adopt estimation methods such as least-squares estimation, weighted least-squares estimation or Kalman filtering to carry out valuation.
The criterion of least-squares estimation is following:
The measurement noise power of considering each frequency is variant, and multipath intensity maybe be also not quite identical, can adopt accurate more weighted least-squares valuation.Its estimation criterion:
W=R in the formula
-1, R is the variance battle array of error in measurement V.
If hope to utilize temporal correlation to obtain higher estimation accuracy, also can use the bigger Kalman filter of operand.Its expulsion mechanism is described by following state equation:
X
k=φ
k,k-1X
k-1+Γ
k-1W
k-1 (41)
In the formula: φ
K, k-1Be t
K-1A step transfer matrix constantly; Γ
K-1For system noise drives battle array; W
kBe the system incentive noise sequence.
Embodiment two, a kind of device that carries out the pseudo range observed quantity estimation through carrier smoothing, as shown in Figure 1, comprising:
Measurement module is used for obtaining respectively the pseudorange and the carrier phase observed quantity of n frequency; N is the integer more than or equal to 2;
Level and smooth module; Be used for respectively carrying out the operation of carrier smoothing pseudorange simultaneously, obtain the pseudorange
after level and smooth at each frequency;
Constructing module is used for the pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;(33)
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is since 1, and the integer till the n comprises 1 and n; d
1And d
2Be constant;
Estimator module is used for according to above-mentioned measurement equation TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
In the present embodiment,
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
In the present embodiment, d
1Span can but be not limited to 40.3 ± 10%, d
2Span can but be not limited to 1 ± 10%.
In the present embodiment, d
1Optimum is 40.3, d
2Optimum is 1.
In the present embodiment, said estimator module can adopt estimation methods such as least-squares estimation, weighted least-squares estimation or Kalman filtering to carry out valuation.
Other variation and realization details can be with embodiment one.
One of ordinary skill in the art will appreciate that all or part of step in the said method can instruct related hardware to accomplish through program, said program can be stored in the computer-readable recording medium, like ROM (read-only memory), disk or CD etc.Alternatively, all or part of step of the foregoing description also can use one or more integrated circuit to realize.Correspondingly, each the module/unit in the foregoing description can adopt the form of hardware to realize, also can adopt the form of software function module to realize.The present invention is not restricted to the combination of the hardware and software of any particular form.
Certainly; The present invention also can have other various embodiments; Under the situation that does not deviate from spirit of the present invention and essence thereof; Those of ordinary skill in the art work as can make various corresponding changes and distortion according to the present invention, but these corresponding changes and distortion all should belong to the protection domain of claim of the present invention.
Claims (10)
1. one kind is carried out the pseudo range observed quantity estimation approach through carrier smoothing, comprising:
Obtain the pseudorange and the carrier phase observed quantity of n frequency respectively; N is the integer more than or equal to 2;
Carry out carrier smoothing pseudorange operation simultaneously at each frequency respectively, obtain the pseudorange
after level and smooth;
Pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is since 1, and the integer till the n comprises 1 and n; d
1And d
2Be constant;
According to above-mentioned measurement equation to TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
2. the method for claim 1 is characterized in that:
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
3. method as claimed in claim 2 is characterized in that:
d
1Span be 40.3 ± 10%, d
2Span be 1 ± 10%.
4. like claim 2 or 3 described methods, it is characterized in that:
d
1Be 40.3, d
2Be 1.
5. the method for claim 1 is characterized in that:
Said according to measurement equation to TEC
EffCarry out in the step of valuation with r, adopt least-squares estimation, weighted least-squares estimation or Kalman filtering to carry out valuation.
6. a device that carries out the pseudo range observed quantity estimation through carrier smoothing is characterized in that, comprising:
Measurement module is used for obtaining respectively the pseudorange and the carrier phase observed quantity of n frequency; N is the integer more than or equal to 2;
Level and smooth module; Be used for respectively carrying out the operation of carrier smoothing pseudorange simultaneously, obtain the pseudorange
after level and smooth at each frequency
Constructing module is used for the pseudorange behind a plurality of frequencies level and smooth is formed measurement equation:
Z=HX+V;
Wherein:
TEC
EffWith r for treating estimator, f
iBe the frequency of frequency i,
For frequency i goes up by the pseudo range measurement noise and the multipath that fully suppress; I is since 1, and the integer till the n comprises 1 and n; d
1And d
2Be constant;
Estimator module, according to above-mentioned measurement equation to TEC
EffCarry out valuation with r, with the valuation of resulting r as the pseudo range observed quantity estimated value.
7. device as claimed in claim 6 is characterized in that:
Wherein, K is the discrete time constant, and t is a time point, and e is the truth of a matter of natural logarithm; τ is the smoothingtime constant, the Δ t of τ=(K-1), and Δ t is the minimum SI of smoothing pseudo range, TEC is a total electron content.
8. device as claimed in claim 7 is characterized in that:
d
1Span be 40.3 ± 10%, d
2Span be 1 ± 10%.
9. like claim 7 or 8 described devices, it is characterized in that:
d
1Be 40.3, d
2Be 1.
10. method as claimed in claim 6 is characterized in that:
Said estimator module adopts least-squares estimation, weighted least-squares to estimate or Kalman filtering is carried out valuation.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010619774 CN102565813B (en) | 2010-12-31 | 2010-12-31 | Method and device for performing pseudorange observation estimation by carrier smoothing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010619774 CN102565813B (en) | 2010-12-31 | 2010-12-31 | Method and device for performing pseudorange observation estimation by carrier smoothing |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102565813A true CN102565813A (en) | 2012-07-11 |
CN102565813B CN102565813B (en) | 2013-10-16 |
Family
ID=46411688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010619774 Active CN102565813B (en) | 2010-12-31 | 2010-12-31 | Method and device for performing pseudorange observation estimation by carrier smoothing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102565813B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792558A (en) * | 2014-01-10 | 2014-05-14 | 中国人民解放军63921部队 | GNSS carrier phase smoothness pseudo-range processing method |
CN103926609A (en) * | 2014-04-30 | 2014-07-16 | 深圳市星耀航威通讯技术有限公司 | Observed quantity smoothing system and method |
CN104516006A (en) * | 2014-11-07 | 2015-04-15 | 中国电子科技集团公司第二十研究所 | Carrier phase smoothing pseudorange algorithm based on improved Kalman filtering |
CN105068097A (en) * | 2015-09-01 | 2015-11-18 | 中国电子科技集团公司第二十研究所 | Self-adaptive filtering method for carrier smoothed code pseudorange |
CN107193026A (en) * | 2017-05-06 | 2017-09-22 | 千寻位置网络有限公司 | Pseudorange positioning smooth method and system, positioning terminal |
CN108267135A (en) * | 2017-12-25 | 2018-07-10 | 中铁第四勘察设计院集团有限公司 | For the accurate positioning method and system of track automatic measurement vehicle |
CN108845338A (en) * | 2018-07-26 | 2018-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Star ground united carrier smoothing pseudo range distance measuring method |
CN109709579A (en) * | 2019-02-21 | 2019-05-03 | 哈尔滨工程大学 | A kind of GNSS satellite ephemeris fault detection method based on user ranging errors real-time estimation |
CN111856512A (en) * | 2020-07-06 | 2020-10-30 | 上海交通大学 | Dual-frequency ionosphere error estimation method and system based on portable terminal |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050024263A1 (en) * | 2003-07-30 | 2005-02-03 | Sharpe Richard T. | Method for generating clock corrections for a wide-area or global differential GPS system |
WO2010050433A1 (en) * | 2008-10-28 | 2010-05-06 | 古野電気株式会社 | Satellite navigation device |
CN101825717A (en) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | Carrier smoothing code pseudorange technology-based dynamic attitude positioning method |
-
2010
- 2010-12-31 CN CN 201010619774 patent/CN102565813B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050024263A1 (en) * | 2003-07-30 | 2005-02-03 | Sharpe Richard T. | Method for generating clock corrections for a wide-area or global differential GPS system |
WO2010050433A1 (en) * | 2008-10-28 | 2010-05-06 | 古野電気株式会社 | Satellite navigation device |
CN101825717A (en) * | 2010-04-16 | 2010-09-08 | 北京航空航天大学 | Carrier smoothing code pseudorange technology-based dynamic attitude positioning method |
Non-Patent Citations (2)
Title |
---|
赵胜 等: "GNSS接收机的伪距平滑技术研究", 《无线电工程》 * |
陈怡 等: "载波相位平滑在GPS伪距提取中的应用", 《航天控制》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792558B (en) * | 2014-01-10 | 2019-05-31 | 中国人民解放军63921部队 | GNSS carrier phase smoothing pseudo-range processing method |
CN103792558A (en) * | 2014-01-10 | 2014-05-14 | 中国人民解放军63921部队 | GNSS carrier phase smoothness pseudo-range processing method |
CN103926609A (en) * | 2014-04-30 | 2014-07-16 | 深圳市星耀航威通讯技术有限公司 | Observed quantity smoothing system and method |
CN104516006A (en) * | 2014-11-07 | 2015-04-15 | 中国电子科技集团公司第二十研究所 | Carrier phase smoothing pseudorange algorithm based on improved Kalman filtering |
CN104516006B (en) * | 2014-11-07 | 2017-04-12 | 中国电子科技集团公司第二十研究所 | Carrier phase smoothing pseudorange algorithm based on improved Kalman filtering |
CN105068097A (en) * | 2015-09-01 | 2015-11-18 | 中国电子科技集团公司第二十研究所 | Self-adaptive filtering method for carrier smoothed code pseudorange |
CN107193026A (en) * | 2017-05-06 | 2017-09-22 | 千寻位置网络有限公司 | Pseudorange positioning smooth method and system, positioning terminal |
CN108267135A (en) * | 2017-12-25 | 2018-07-10 | 中铁第四勘察设计院集团有限公司 | For the accurate positioning method and system of track automatic measurement vehicle |
CN108845338A (en) * | 2018-07-26 | 2018-11-20 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Star ground united carrier smoothing pseudo range distance measuring method |
CN108845338B (en) * | 2018-07-26 | 2022-04-01 | 西南电子技术研究所(中国电子科技集团公司第十研究所) | Satellite-ground combined carrier smoothing pseudorange ranging method |
CN109709579A (en) * | 2019-02-21 | 2019-05-03 | 哈尔滨工程大学 | A kind of GNSS satellite ephemeris fault detection method based on user ranging errors real-time estimation |
CN109709579B (en) * | 2019-02-21 | 2023-02-14 | 哈尔滨工程大学 | GNSS satellite ephemeris fault detection method based on user ranging error real-time estimation |
CN111856512A (en) * | 2020-07-06 | 2020-10-30 | 上海交通大学 | Dual-frequency ionosphere error estimation method and system based on portable terminal |
CN111856512B (en) * | 2020-07-06 | 2022-03-11 | 上海交通大学 | Dual-frequency ionosphere error estimation method and system based on portable terminal |
Also Published As
Publication number | Publication date |
---|---|
CN102565813B (en) | 2013-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102565813B (en) | Method and device for performing pseudorange observation estimation by carrier smoothing | |
RU2354991C2 (en) | Method of using three gps frequencies for resolving integral uncertainties of carrier phases | |
CN103344978B (en) | Area enhanced precision positioning service method suitable for large-scale users | |
US7423583B2 (en) | Method for refining state estimates of a global positioning system (GPS) receiver based upon sequential GPS carrier phase measurements | |
JP2013019893A (en) | Error detection for satellite navigation system based on biased measurement | |
JP2010528320A (en) | Reduction of distance-dependent error in real-time kinematic (RTK) positioning | |
TW200821613A (en) | System and/or method for acquisition of GNSS signals | |
Krypiak-Gregorczyk et al. | Carrier phase bias estimation of geometry-free linear combination of GNSS signals for ionospheric TEC modeling | |
Rabbou et al. | Precise point positioning using multi-constellation GNSS observations for kinematic applications | |
Rahemi et al. | Accurate solution of navigation equations in GPS receivers for very high velocities using pseudorange measurements | |
JP7153842B2 (en) | Positioning method and positioning terminal | |
Gowdayyanadoddi et al. | A Ray‐Tracing Technique to Characterize GPS Multipath in the Frequency Domain | |
Zhang et al. | Performance assessment of GNSS diffraction models in urban areas | |
US20050203702A1 (en) | Method for backup dual-frequency navigation during brief periods when measurement data is unavailable on one of two frequencies | |
Liu | A new approach for cycle slip detection and fix using single GPS receiver’s single satellite dual frequency data containing arbitrarily large pseudorange errors | |
JP5670649B2 (en) | How to determine the location of a mobile device | |
García Fernández | Contributions to the 3D ionospheric sounding with GPS data | |
Blanch et al. | Application of spatial statistics to ionosphere estimation for WAAS | |
O’Keefe et al. | Availability, accuracy, reliability, and carrier-phase ambiguity resolution with Galileo and GPS | |
US9086479B2 (en) | Convergence zone | |
He et al. | A comprehensive method for multipath performance analysis of GNSS navigation signals | |
JP7088846B2 (en) | Methods and systems for satellite signal processing | |
Avram et al. | Experimental results of multipath behavior for GPS L1-L2 and Galileo E1-E5b in static and kinematic scenarios | |
Hoque et al. | Transionospheric Microwave Propagation: Higher-Order Effects up to 100 GHz | |
Yang et al. | An Enhanced Method for Detecting and Repairing the Cycle Slips of Dual‐Frequency Onboard GPS Receivers of LEO Satellites |
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 |