CN108919634A - A kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou - Google Patents
A kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou Download PDFInfo
- Publication number
- CN108919634A CN108919634A CN201810914774.3A CN201810914774A CN108919634A CN 108919634 A CN108919634 A CN 108919634A CN 201810914774 A CN201810914774 A CN 201810914774A CN 108919634 A CN108919634 A CN 108919634A
- Authority
- CN
- China
- Prior art keywords
- time
- frequency
- ionosphere
- vtec
- model
- 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
Classifications
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R20/00—Setting the time according to the time information carried or implied by the radio signal
- G04R20/02—Setting the time according to the time information carried or implied by the radio signal the radio signal being sent by a satellite, e.g. GPS
Abstract
The invention discloses a kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou, this method includes:Obtain the parameters such as carrier phase observable, Pseudo-range Observations, satellite ephemeris, earth rotation, the antenna phase center of Chinese Beidou satellite navigation system at two passing time stations;Carry out data inspection, elimination of rough difference and Detection of Cycle-slip;Carry out ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna phase deviation equal error model correction;Construct the ionosphere virtual observation model of additional ionosphere prior information constraint, space region constraint and time region constraint;The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, One-Point Location time service resolving is carried out to the data after correction, it obtains the comprehensive clock deviation of two passing time station receivers and is compared, it is poor to obtain Time Transmission, it is compared with the standard time at one of passing time station, obtains the chronometer time at another passing time station.The invention can reduce observation noise, improve Time Transmission precision and reliability.
Description
Technical field
The present invention relates to big-dipper satellite passing time technical field, especially a kind of three non-non-combined observation of difference of frequency of Beidou
Time Transmission system and method.
Background technique
Time Transmission is that the retention time is synchronous between time laboratory and establishes and maintain the important set of standard time scale
At part.Time Transmission precision is an important indicator of measure time net synchronization capability, is affected synchronous between each time reference
Precision.Time delivering method has total view method Time Transmission and full view method Time Transmission, obtains in many time laboratories extensively
Using still, they are all based on Pseudo-range Observations, it is also necessary to be based on satellite common vision condition, sphere of action is limited, and traditional is complete
No ionospheric combination is usually used depending on method, is exaggerated observation noise, and whole positioning time service precision and reliability are lower.
Summary of the invention
The object of the present invention is to provide a kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou, energy
Observation noise is reduced, Time Transmission precision and reliability are improved.
To achieve the above object, the present invention provides following schemes:
A kind of three non-non-combined observation time delivering method of difference of frequency of Beidou, including:
Obtain carrier phase observable, Pseudo-range Observations, the satellite of the Chinese Beidou satellite navigation system at two passing time stations
Ephemeris, earth rotation, antenna phase center with etc. parameters;
Data inspection, elimination of rough difference and Detection of Cycle-slip are carried out to the carrier phase observable, Pseudo-range Observations, pre-processed
Data afterwards;
Ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna phase deviation etc. are carried out to the pretreated data
Error model correction;
Construct the ionosphere virtual observation mould of additional ionosphere prior information constraint, space region constraint and time region constraint
Type;
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and One-Point Location time service solution is carried out to the data after correction
It calculates, obtains the comprehensive clock deviation of the respective receiver in two passing time stations;
The comprehensive clock deviation of the respective receiver in two passing time stations is compared, it is poor to obtain Time Transmission;
It is compared according to the standard time at the Time Transmission difference and one of passing time station, when obtaining another transmitting
Between the chronometer time stood.
Optionally, described that ephemeris, tide, the theory of relativity, earth rotation, atmosphere, day are carried out to the pretreated data
Phase of line deviation equal error model correction, specifically includes:
It is modified just based on Precise Orbit product to orbit error, precise clock correction is repaired based on precise clock correction product
Just;
By additional ionosphere prior information, time-domain and the related constraint of spatial domain and ionosphere delay error as not
Know that parameter carries out real-time estimation;
Earth tide, tide, Ghandler motion, relativistic effect, antenna phase center variation variation, phase are turned around using corresponding
Model is corrected.
Optionally, the virtual observation mould for constructing the constraint of ionosphere prior information, space region constraint and time region constraint
Type specifically includes:
Using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error,
σpriorValue is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,With
λ0For the latitude and longitude of survey station, EijThe coefficient of representative model.
Using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value model
It encloses for 0.009-0.025m2。
Optionally, the three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and the data after correction are carried out single
Point location time service resolves, and obtains the comprehensive clock deviation of the respective receiver in two passing time stations, specifically includes:
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction,
Unit is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
It is as follows to establish stochastic model:
One-Point Location time service resolving is carried out in conjunction with the non-combined Static Precise Point Positioning model of the non-difference of three frequencies of foundation, obtains two stations
The comprehensive clock deviation of reception;
Wherein, I is unit matrix;σ0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian.
A kind of three non-non-combined observation Time Transmission system of difference of frequency of Beidou, including:
Observation obtains module, the phase observations of the Chinese Beidou satellite navigation system for obtaining two passing time stations
The parameters such as value, Pseudo-range Observations, satellite ephemeris, earth rotation, antenna phase center;
Preprocessing module, for carrying out data inspection, elimination of rough difference and cycle slip to the carrier phase observable, Pseudo-range Observations
Detection, obtains pretreated data;
Correct module, for the pretreated data carry out ephemeris, tide, the theory of relativity, earth rotation, atmosphere,
Antenna phase deviation equal error model correction;
Virtual observation model construction module, for constructing additional ionosphere prior information constraint, space region constraint and time
The ionosphere virtual observation model of region constraint;
Module is resolved, for constructing the three non-non-combined Static Precise Point Positioning models of difference of frequency, and the data after correction are carried out
One-Point Location time service resolves, and obtains the comprehensive clock deviation of the respective receiver in two passing time stations;
Time Transmission difference computing module, for the comprehensive clock deviation of the respective receiver in two passing time stations to be compared,
It is poor to obtain Time Transmission;
Passing time module is compared according to the standard time at the Time Transmission difference and one of passing time station, is obtained
To the chronometer time at another passing time station.
Optionally, the correction module, specifically includes:
Orbit parameter amending unit is modified orbit error based on Precise Orbit product;
Satellite ephemeris error amending unit is modified clock deviation error based on precise clock correction product;
Ionosphere delay error amending unit, it is related for ionosphere prior information, time-domain and spatial domain will to be added
Constraint and ionosphere delay error are as unknown parameter progress real-time estimation;
Satellite navigation relevant parameter amending unit, for in earth tide, tide, Ghandler motion, relativistic effect, antenna phase
Heart change of error, phase are turned around to be corrected using corresponding model.
Optionally, the virtual observation model construction module, including:
Based on grid ionospheric model, prior information constraint is carried out, prior information constraint equation is obtained;According to ionosphere sky
Between characteristic, carry out single station modeling, solving model coefficient forms spatial domain constraint equation;It is right according to ionosphere time behavior
Adjacent epoch is constrained, and is formed time-domain constraint equation, is specifically included:
Ionosphere prior information restricted model construction unit, for using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error,
σpriorValue is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Ionosphere spatial domain restricted model construction unit, for using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,With
λ0For the latitude and longitude of survey station, EijThe coefficient of representative model.
Ionosphere time-domain restricted model construction unit, for using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value model
It encloses for 0.009-0.025m2。
Optionally, the resolving module, including:
It is linearized by model correction and parameter, constructs observational equation, and determine stochastic model, carry out parameter Estimation, obtain
The comprehensive clock deviation of receiver to two stations, specifically includes:
The three non-non-combined Static Precise Point Positioning model construction units of difference of frequency, for constructing the three non-non-combined accurate one-points of difference of frequency
Location model is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction,
Unit is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
Stochastic model establishes unit, as follows for establishing stochastic model:
One-Point Location time service resolving is carried out in conjunction with the non-combined Static Precise Point Positioning model of the non-difference of three frequencies of foundation, obtains two stations
The comprehensive clock deviation of receiver;
Wherein, I is unit matrix;σ0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
The present invention provides a kind of three non-non-combined observation Time Transmission system and methods of difference of frequency of Beidou, based on GPS essence
Close One-Point Location (GPS Precise Point Positioning, PPP) technology is combined the three non-non-combined sight of difference of frequency of Beidou
The time delivering method of measured value passes through additional ionosphere prior information constraint, space region constraint and time region constraint virtual observation
Equation constructs the non-combined PPP model of non-difference of three frequency signals, by the shared dipper system time, realizes the precision at two stations
Time Transmission service, effectively reduces observation noise, improves Time Transmission precision, and available ionosphere and hardware delay
Product.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain his attached drawing.
Fig. 1 is the flow diagram of 1 Beidou of the embodiment of the present invention, the three non-non-combined observation time delivering method of difference of frequency;
Fig. 2 is the structural schematic diagram of 2 Beidou of the embodiment of the present invention, the three non-non-combined observation Time Transmission system of difference of frequency;
Fig. 3 is the flow diagram of 3 Beidou of the embodiment of the present invention, the three non-non-combined observation time delivering method of difference of frequency.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, obtained by those of ordinary skill in the art without making creative efforts all he is real
Example is applied, shall fall within the protection scope of the present invention.
The object of the present invention is to provide a kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou, energy
Observation noise is reduced, Time Transmission precision and reliability are improved.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Applying mode, the present invention is described in further detail.
Fig. 1 is the flow diagram of the three non-non-combined observation time delivering method of difference of frequency of Beidou of the present invention.
As shown in Figure 1, a kind of three non-non-combined observation time delivering method of difference of frequency of Beidou, including:
Step 101:Obtain carrier phase observable, the pseudorange observation of the Chinese Beidou satellite navigation system at two passing time stations
The parameters such as value, satellite ephemeris, earth rotation, antenna phase center;
Step 102:Data inspection, elimination of rough difference and Detection of Cycle-slip are carried out to the carrier phase observable, Pseudo-range Observations, obtained
To pretreated data;
Step 103:Ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna phase are carried out to the pretreated data
Position deviation equal error model correction;
Step 104:The ionosphere for constructing additional ionosphere prior information constraint, space region constraint and time region constraint is virtual
Observation model;
Step 105:The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and it is fixed to carry out single-point to the data after correction
Position time service resolves, and obtains the comprehensive clock deviation of the respective receiver in two passing time stations;
Step 106:The comprehensive clock deviation of the respective receiver in two passing time stations is compared, it is poor to obtain Time Transmission;
Step 107:It is compared, is obtained another according to the standard time at the Time Transmission difference and one of passing time station
The chronometer time at a passing time station.
The step 103:Ephemeris, tide, the theory of relativity, earth rotation, atmosphere, day are carried out to the pretreated data
Phase of line deviation equal error model correction, specifically includes:
Orbit error is modified by Precise Orbit product, and clock deviation error is modified by precise clock correction product;
By additional ionosphere prior information, time-domain and the related constraint of spatial domain and ionosphere delay error as not
Know that parameter carries out real-time estimation;
Earth tide, tide, Ghandler motion, relativistic effect, antenna phase center variation variation, phase are turned around using corresponding
Model is corrected.
The step 104:Construct the ionosphere of additional ionosphere prior information constraint, space region constraint and time region constraint
Virtual observation model, including:
Based on grid ionospheric model, prior information constraint is carried out, prior information constraint equation is obtained;According to ionosphere sky
Between characteristic, carry out single station modeling, solving model coefficient forms spatial domain constraint equation;It is right according to ionosphere time behavior
Adjacent epoch is constrained, and is formed time-domain constraint equation, is specifically included:
Using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error,
σpriorValue is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,With
λ0For the latitude and longitude of survey station, EijThe coefficient of representative model.
Using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value model
It encloses for 0.009-0.025m2。
The step 105:The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and the data after correction are carried out single
Point location time service resolves, and obtains the comprehensive clock deviation of the respective receiver in two passing time stations, including:
It is linearized by model correction and parameter, constructs observational equation, and determine stochastic model, carry out parameter Estimation, obtain
To the synthesis clock deviation at two stations, specifically include:
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction, single
Position is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
It is as follows to establish stochastic model:
One-Point Location time service resolving is carried out in conjunction with the non-combined Static Precise Point Positioning model of the non-difference of three frequencies of foundation, obtains two stations
The comprehensive clock deviation of reception;
Wherein, I is unit matrix;σ0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian.
Fig. 2 is the structural schematic diagram of the three non-non-combined observation Time Transmission system of difference of frequency of Beidou of the present invention.
As shown in Fig. 2, a kind of three non-non-combined observation Time Transmission system of difference of frequency of Beidou, including:
Observation obtains module 201, the phase of the Chinese Beidou satellite navigation system for obtaining two passing time stations
The parameters such as observation, Pseudo-range Observations, satellite ephemeris, earth rotation, antenna phase center;
Preprocessing module 202, for the carrier phase observable, Pseudo-range Observations carry out data inspection, elimination of rough difference and
Detection of Cycle-slip obtains pretreated data;
Module 203 is corrected, for carrying out ephemeris, tide, the theory of relativity, earth rotation, big to the pretreated data
Gas, antenna phase deviation equal error model correction;
Virtual observation model construction module 204, for construct the constraint of additional ionosphere prior information, space region constraint and when
Between region constraint ionosphere virtual observation model;
Module 205, the three non-non-combined Static Precise Point Positioning model of difference of frequency of building are resolved, and the data after correction are carried out single
Point location time service resolves, and obtains the comprehensive clock deviation of the respective receiver in two passing time stations;
Time Transmission difference computing module 206, for comparing the comprehensive clock deviation of the respective receiver in two passing time stations
It is right, it is poor to obtain Time Transmission;
Passing time module 207 is compared according to the standard time at the Time Transmission difference and one of passing time station,
Obtain the chronometer time at another passing time station.
The correction module 203, specifically includes:
Orbit parameter amending unit carries out orbit error amendment based on Precise Orbit product;
Satellite ephemeris error amending unit carries out clock deviation error correction based on precise clock correction product;
Ionosphere delay error amending unit, it is related for ionosphere prior information, time-domain and spatial domain will to be added
Constraint and ionosphere delay error are as unknown parameter progress real-time estimation;
Satellite navigation relevant parameter amending unit, for in earth tide, tide, Ghandler motion, relativistic effect, antenna phase
Heart change of error, phase are turned around to be corrected using corresponding model.
The virtual observation model construction module 204, including:
Based on grid ionospheric model, prior information constraint is carried out, prior information constraint equation is obtained;According to ionosphere sky
Between characteristic, carry out single station modeling, solving model coefficient forms spatial domain constraint equation;It is right according to ionosphere time behavior
Adjacent epoch is constrained, and is formed time-domain constraint equation, is specifically included:
Ionosphere prior information restricted model construction unit, for using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error,
σpriorValue is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Ionosphere spatial domain restricted model construction unit, for using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,With
λ0For the latitude and longitude of survey station, EijThe coefficient of representative model.
Ionosphere time-domain restricted model construction unit, for using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value model
It encloses for 0.009-0.025m2。
The resolving module 205, including:
It is linearized by model correction and parameter, constructs observational equation, and determine stochastic model, carry out parameter Estimation, obtain
To the synthesis clock deviation at two stations, specifically include:
It is the foundation of function model first for model foundation, ionosphere delay uses the parameter with constraint condition to estimate
Meter method, tropospheric delay error are corrected using empirical model, and nubbin uses piecewise constant or random walk mould
Type is estimated.
The three non-non-combined Static Precise Point Positioning model construction units of difference of frequency, for constructing the three non-non-combined accurate one-points of difference of frequency
Location model is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction, single
Position is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
Stochastic model establishes unit, as follows for establishing stochastic model:
One-Point Location time service resolving is carried out in conjunction with the non-combined Static Precise Point Positioning model of the non-difference of three frequencies of foundation, obtains two stations
The comprehensive clock deviation of reception;
Wherein, I is unit matrix;σ0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian;
The hardware delay deviation D CB of satellite end is corrected using the product that IGS is announced, and the DCB of receiver end is joined
Number estimation;Satellite ephemeris error is corrected using the correction of real-time broadcasting, and receiver clock-offsets first carry out time-difference parameter correction, then when
Estimated as white Gaussian noise to get local zone time of the user relative to the NTSC reference time is arrived.Multipath effect does not have temporarily
There are reliable model or method correction, can be handled as observation noise.Its stochastic model can be according to pseudorange, carrier phase observable and mould
Type precision, and elevation of satellite is cooperated to carry out comprehensive determination.
Fig. 3 is the flow diagram of 3 Beidou of the embodiment of the present invention, the three non-non-combined observation time delivering method of difference of frequency.Such as
Shown in Fig. 3, the observation data acquired on survey station are collected, what Precise Orbit, clock deviation, ionosphere product and data processing needed
Assist product (earth rotation parameter (ERP), DCB correction, day line file etc.);
Based on grid ionospheric model, prior information constraint is carried out, prior information constraint equation is obtained;According to ionosphere sky
Between characteristic, carry out single station modeling, solving model coefficient forms spatial domain constraint equation;It is right according to ionosphere time behavior
Adjacent epoch is constrained, and time-domain constraint equation is formed.
It is linearized by model correction and parameter, constructs observational equation, and determine stochastic model, carry out parameter Estimation, obtain
The comprehensive clock deviation of receiver to each station.
The comparison that comprehensive clock deviation is carried out between each station, obtains that Time Transmission is poor, relies on the standard time of reference station, Ji Keshi
The Time Transmission function of existing subscriber station.
Using high-precision phase position observation, the precision of Time Transmission is improved.
The three non-non-combined observation Time Transmission system and method for difference of frequency of a kind of Beidou of the invention, high-precision phase are seen
Measured value carries out Time Transmission calculating, compensates for routine and regards method altogether, the limitation of low precision Pseudo-range Observations is used only in full view method,
The precision of Time Transmission will be greatly promoted;Using three frequency observations, whole Time Transmission precision and reliability can be improved, made
With three frequency observations, not only there is combination of frequency abundant, but also Time Transmission precision and reliability can be improved;Using non-difference
Non-combined model, not only greatly reduces observation noise, improves the precision of Time Transmission;Ionization abundant is also obtained simultaneously
Layer and hardware delay product.
Used herein a specific example illustrates the principle and implementation of the invention, and above embodiments are said
It is bright to be merely used to help understand method and core concept of the invention;At the same time, for those skilled in the art, according to this
The thought of invention, there will be changes in the specific implementation manner and application range.In conclusion the content of the present specification is not answered
It is interpreted as limitation of the present invention.
Claims (8)
1. a kind of three non-non-combined observation time delivering method of difference of frequency of Beidou, which is characterized in that including:
Obtain the carrier phase observable of the Chinese Beidou satellite navigation system at two passing time stations, Pseudo-range Observations, satellite ephemeris,
Earth rotation, antenna phase center parameter;
Data inspection, elimination of rough difference and Detection of Cycle-slip are carried out to the carrier phase observable, Pseudo-range Observations, obtained pretreated
Data;
Ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna phase deviation equal error are carried out to the pretreated data
Model correction;
Construct the ionosphere virtual observation model of additional ionosphere prior information constraint, space region constraint and time region constraint;
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and One-Point Location time service resolving is carried out to the data after correction,
Obtain the comprehensive clock deviation of the respective receiver in two passing time stations;
The comprehensive clock deviation of the respective receiver in two passing time stations is compared, it is poor to obtain Time Transmission;
It is compared according to the standard time at the Time Transmission difference and one of passing time station, obtains another passing time station
Chronometer time.
2. the three non-non-combined observation time delivering method of difference of frequency of a kind of Beidou according to claim 1, which is characterized in that
It is described that ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna phase deviation equal error are carried out to the pretreated data
Model correction specifically includes:
Orbit error is corrected by Precise Orbit, clock deviation error is corrected by precise clock correction;
By additional ionosphere prior information, time-domain and the related constraint of spatial domain and ionosphere delay error as unknown ginseng
Number carries out real-time estimation;
Earth tide, tide, Ghandler motion, relativistic effect, antenna phase center variation variation, phase are turned around using corresponding model
It is corrected.
3. the three non-non-combined observation time delivering method of difference of frequency of a kind of Beidou according to claim 1, which is characterized in that
The virtual observation model for constructing the constraint of ionosphere prior information, space region constraint and time region constraint, specifically includes:
Using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error, σpriorIt takes
Value is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,And λ0For
The latitude and longitude of survey station, EijThe coefficient of representative model.
Using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value range is
0.009-0.025m2。
4. the three non-non-combined observation time delivering method of difference of frequency of a kind of Beidou according to claim 1, which is characterized in that
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building, and One-Point Location time service resolving is carried out to the data after correction,
The comprehensive clock deviation of the respective receiver in two passing time stations is obtained, is specifically included:
The three non-non-combined Static Precise Point Positioning model of difference of frequency of building is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction,
Unit is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
It is as follows to establish stochastic model:
One-Point Location time service solution is carried out to the amendment data in conjunction with the three non-non-combined Static Precise Point Positioning models of difference of frequency of foundation
It calculates, obtains the comprehensive clock deviation of reception at two stations;
Wherein, I is unit matrix;σ0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian.
5. a kind of three non-non-combined observation Time Transmission system of difference of frequency of Beidou, which is characterized in that including:
Observation acquisition module, the carrier phase observable of the Chinese Beidou satellite navigation system for obtaining two passing time stations,
Pseudo-range Observations, satellite ephemeris, earth rotation, antenna phase center parameter;
Preprocessing module is visited for carrying out data inspection, elimination of rough difference and cycle slip to the carrier phase observable, Pseudo-range Observations number
It surveys, obtains pretreated data;
Module is corrected, for carrying out ephemeris, tide, the theory of relativity, earth rotation, atmosphere, antenna to the pretreated data
Phase deviation equal error model correction;
Virtual observation model construction module, for constructing additional ionosphere prior information constraint, space region constraint and time-domain about
The ionosphere virtual observation model of beam;
Module is resolved, for constructing the three non-non-combined Static Precise Point Positioning models of difference of frequency, and single-point is carried out to the data after correction
It positions time service to resolve, obtains the comprehensive clock deviation of the respective receiver in two passing time stations;
Time Transmission difference computing module is obtained for the comprehensive clock deviation of the respective receiver in two passing time stations to be compared
Time Transmission is poor;
Passing time module is obtained for being compared according to the standard time of the Time Transmission difference and one of passing time station
To the chronometer time at another passing time station.
6. the three non-non-combined observation Time Transmission system of difference of frequency of a kind of Beidou according to claim 5, which is characterized in that
The correction module, specifically includes:
Orbit parameter amending unit, for the amendment to orbit error;
Satellite ephemeris error amending unit, for the amendment to clock deviation error;
Ionosphere delay error amending unit, for ionosphere prior information, time-domain and the related constraint of spatial domain will to be added
And ionosphere delay error carries out real-time estimation as unknown parameter;
Satellite navigation relevant parameter amending unit, for inclined to earth tide, tide, Ghandler motion, relativistic effect, antenna phase center
Difference variation, phase are turned around to be corrected using corresponding model.
7. the three non-non-combined observation Time Transmission system of difference of frequency of a kind of Beidou according to claim 5, which is characterized in that
The virtual observation model construction module, specifically includes:
Ionosphere prior information restricted model construction unit, for using such as drag:VTEC=VTECGIM+εGIM;
Construct the constraint of ionosphere prior information;
Wherein, VTECGIMIndicate the total electron content that grid ionospheric model extracts, εGIMIndicate prior model error, σpriorIt takes
Value is 0.3-0.6 meters, and B indicates that the latitude of ionosphere point of puncture, t indicate the local time as unit of hour;
Ionosphere spatial domain restricted model construction unit, for using such as drag:
VTEC=VTECspace+εspace;
Establish ionosphere space region constraint;
Wherein, m and n is the end of surface model, and usual value is 2,With λ be point of puncture latitude and longitude,And λ0For
The latitude and longitude of survey station, EijThe coefficient of representative model.
Ionosphere time-domain restricted model construction unit, for using such as drag:
VTEC=VTEClast+ΔVTEC+εtemp;
Establish ionosphere time region constraint;
Wherein, the epoch that Δ VTEC is ionosphere VTEC changes,For the variance of Δ VTEC, unit m2, value range is
0.009-0.025m2。
8. the three non-non-combined observation Time Transmission system of difference of frequency of a kind of Beidou according to claim 5, which is characterized in that
The resolving module, specifically includes:
The three non-non-combined Static Precise Point Positioning model construction units of difference of frequency, for constructing the three non-non-combined Static Precise Point Positionings of difference of frequency
Model is as follows:
Wherein, P and φ is respectively pseudorange and carrier phase observation data, and f is frequency;Coefficient of frequencyρ is geometric distance of the satellite to survey station, c
For the light velocity, dt12For P1And P2The receiver clock-offsets of no ionospheric combination observation, dtsFor satellite clock correction, dtropFor tropospheric delay;
M is represented and earth tide, tide, earth rotation, relativistic effect, antenna phase center variation, turns around and multipath composition error;
N is fuzziness item, and the hardware delay deviation comprising satellite and receiver end, ε is observation noise, and IFB is receiver inter-frequency deviation;
DCB is pseudorange code inter-frequency deviation, and the code deviation of satellite end is provided by international GNSS service centre, and is possessed between different frequencyRelationship, dionFor the ionosphere delay in L1 frequency, dion=f (Z) *
VTEC, footmark behalf satellite,Z is satellite zenith angle, and VTEC is the total electron content of zenith direction,
Unit is TECU, and f (Z) is the mapping function for converting the VTEC of zenith direction to oblique path direction;
Stochastic model establishes unit, as follows for establishing stochastic model:
One-Point Location time service solution is carried out to the amendment data in conjunction with the three non-non-combined Static Precise Point Positioning models of difference of frequency of foundation
It calculates, obtains the comprehensive clock deviation of receiver at two stations;
Wherein, I is unit matrix;σ 0=a/sin (E) is error in unit power;A is constant, and phase value is 0.002-
0.003m, pseudorange value are 0.2-2.0m;E is elevation of satellite, and unit is radian.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810914774.3A CN108919634B (en) | 2018-08-13 | 2018-08-13 | Beidou three-frequency non-differential non-combined observation value time transmission system and method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810914774.3A CN108919634B (en) | 2018-08-13 | 2018-08-13 | Beidou three-frequency non-differential non-combined observation value time transmission system and method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108919634A true CN108919634A (en) | 2018-11-30 |
CN108919634B CN108919634B (en) | 2020-08-07 |
Family
ID=64404483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810914774.3A Active CN108919634B (en) | 2018-08-13 | 2018-08-13 | Beidou three-frequency non-differential non-combined observation value time transmission system and method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108919634B (en) |
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109856656A (en) * | 2019-02-14 | 2019-06-07 | 上海华测导航技术股份有限公司 | A kind of navigation locating method, device, electronic equipment and storage medium |
CN110007328A (en) * | 2019-05-10 | 2019-07-12 | 国网浙江省电力有限公司信息通信分公司 | Non-combined RTK localization method based on four frequency observation of No. three satellites of Beidou |
CN110018507A (en) * | 2019-05-08 | 2019-07-16 | 中国科学院国家授时中心 | It is a kind of based on make between constellation difference combination accurate one-point positioning method and system |
CN110221320A (en) * | 2019-05-07 | 2019-09-10 | 中国科学院国家授时中心 | A kind of Beidou inter-frequency deviation measuring method based on parabola antenna observation |
CN110286396A (en) * | 2019-05-20 | 2019-09-27 | 苏州时空复弦网络科技有限公司 | A kind of non-combined PPP method of non-difference based on ionosphere delay prior information and the double constraints of change in time and space information |
CN110780588A (en) * | 2019-10-16 | 2020-02-11 | 北京航空航天大学 | Wide-area accurate time service WPT system and method |
CN111323748A (en) * | 2018-12-13 | 2020-06-23 | 千寻位置网络有限公司 | Differential positioning method and system |
CN111367160A (en) * | 2020-03-30 | 2020-07-03 | 武汉大学 | GNSS time service receiver precision clock regulation and control method and system |
CN111447030A (en) * | 2020-03-31 | 2020-07-24 | 中国科学院地质与地球物理研究所 | Network time service precision control method and system of earthquake collector and terminal equipment |
CN111766616A (en) * | 2020-06-15 | 2020-10-13 | 中国人民解放军61081部队 | Beidou second-order time transfer satellite-side multipath error correction method |
CN111831965A (en) * | 2019-04-23 | 2020-10-27 | 千寻位置网络有限公司 | Integrity monitoring method and device for UPD correction number in SSR |
CN112859120A (en) * | 2021-01-11 | 2021-05-28 | 中国科学院国家授时中心 | Continuous GNSS carrier phase time and frequency transfer method |
CN113099381A (en) * | 2021-04-06 | 2021-07-09 | 苏州迭慧智能科技有限公司 | Antenna industrial parameter intelligent sensor and multifunctional intelligent sensing network |
CN113325446A (en) * | 2021-06-01 | 2021-08-31 | 中国科学院国家授时中心 | Multi-mode common-frequency GNSS carrier phase time transfer method and system |
WO2021174482A1 (en) * | 2020-03-05 | 2021-09-10 | 上海华测导航技术股份有限公司 | Method and device for converting state space reputation into observation space reputation |
CN113568020A (en) * | 2021-09-27 | 2021-10-29 | 长沙学院 | Satellite navigation positioning error correction method and device considering hardware inter-frequency difference |
CN113608248A (en) * | 2021-06-25 | 2021-11-05 | 北京建筑大学 | Beidou 5G fused high-precision routing inspection personnel positioning method and related equipment |
CN113671545A (en) * | 2021-08-27 | 2021-11-19 | 中国科学院国家授时中心 | Satellite-ground precise time synchronization and carrier phase cycle slip detection method |
CN114002723A (en) * | 2021-09-14 | 2022-02-01 | 中国科学院国家授时中心 | Beidou GEO satellite precise time service system and method based on directional antenna |
CN114280644A (en) * | 2021-12-22 | 2022-04-05 | 广州南方卫星导航仪器有限公司 | PPP-B2B service-based precise point positioning system and method |
CN114355419A (en) * | 2021-12-15 | 2022-04-15 | 中国科学院国家授时中心 | Distributed Beidou position service center RTK product positioning method and positioning device |
CN114397683A (en) * | 2021-12-30 | 2022-04-26 | 深圳市华芯云创科技有限公司 | Beidou non-differential non-combination PPP-RTK positioning method |
CN114740513A (en) * | 2022-05-20 | 2022-07-12 | 航天宏图信息技术股份有限公司 | GNSS positioning method and device |
CN114911153A (en) * | 2022-03-31 | 2022-08-16 | 中国科学院国家授时中心 | Universal time UT1 determination method based on satellite-ground ultra-long baseline interferometry |
CN114995098A (en) * | 2022-04-22 | 2022-09-02 | 中国科学院国家授时中心 | Satellite interference time transfer method based on picosecond-level precision of phase delay measurement |
CN115201870A (en) * | 2022-07-04 | 2022-10-18 | 中国科学院精密测量科学与技术创新研究院 | Multi-frequency multi-mode GNSS non-differential non-combination time transfer method with prior constraint |
CN115225245A (en) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | Non-differential non-combination PPP-RTK long baseline real-time transfer method |
CN115407371A (en) * | 2022-09-02 | 2022-11-29 | 中国人民解放军国防科技大学 | PPP-B2B-based real-time high-precision time transfer method and device |
CN115453851A (en) * | 2022-09-01 | 2022-12-09 | 广东电网有限责任公司 | BD 3-based multi-satellite comprehensive time service method, system, device and medium |
CN115856945A (en) * | 2023-02-03 | 2023-03-28 | 北京航空航天大学 | Precise single-point time service method, device and medium suitable for GNSS multi-frequency point |
CN116184441A (en) * | 2022-12-20 | 2023-05-30 | 北京航空航天大学 | Precise atmospheric correction enhanced intelligent terminal PPP instantaneous sub-meter positioning method |
CN116243591A (en) * | 2023-01-28 | 2023-06-09 | 北京航空航天大学 | Subnanosecond time service method integrating UTC (k) and Beidou broadcast ephemeris |
CN116299618A (en) * | 2023-03-24 | 2023-06-23 | 中国科学院精密测量科学与技术创新研究院 | Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters |
CN116540282A (en) * | 2023-05-08 | 2023-08-04 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
CN117031507A (en) * | 2023-10-08 | 2023-11-10 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Precise single point positioning method suitable for BDS-3B 1CB2a double-frequency signals |
CN117310770A (en) * | 2023-11-28 | 2023-12-29 | 武汉理工大学三亚科教创新园 | Sea level height inversion method based on Beidou ionosphere-free combined precise single-point positioning |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100278220A1 (en) * | 2007-12-14 | 2010-11-04 | Magellan Systems Japan, Inc. | Process for sub-microsecond time transfer using weak gps/gnss signals |
CN103728876A (en) * | 2014-01-02 | 2014-04-16 | 东南大学 | Satellite clock error real-time estimation method based on area multi-reference-station joint resolving |
CN104503223A (en) * | 2014-12-17 | 2015-04-08 | 同济大学 | GNSS (Global Navigation Satellite System) three-frequency high-precision satellite clock correction estimating and service method |
CN106405576A (en) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | Ionosphere monitoring and modeling method based on carrier phase observed value of GNSS |
CN107390510A (en) * | 2017-08-22 | 2017-11-24 | 中国科学院国家授时中心 | Beidou navigation satellite precise time service unified approach based on carrier phase observable |
-
2018
- 2018-08-13 CN CN201810914774.3A patent/CN108919634B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100278220A1 (en) * | 2007-12-14 | 2010-11-04 | Magellan Systems Japan, Inc. | Process for sub-microsecond time transfer using weak gps/gnss signals |
CN103728876A (en) * | 2014-01-02 | 2014-04-16 | 东南大学 | Satellite clock error real-time estimation method based on area multi-reference-station joint resolving |
CN104503223A (en) * | 2014-12-17 | 2015-04-08 | 同济大学 | GNSS (Global Navigation Satellite System) three-frequency high-precision satellite clock correction estimating and service method |
CN106405576A (en) * | 2016-08-25 | 2017-02-15 | 中国科学院国家授时中心 | Ionosphere monitoring and modeling method based on carrier phase observed value of GNSS |
CN107390510A (en) * | 2017-08-22 | 2017-11-24 | 中国科学院国家授时中心 | Beidou navigation satellite precise time service unified approach based on carrier phase observable |
Non-Patent Citations (2)
Title |
---|
张如飞 等: "BDS PPP高精度时间传递方法与精度比较", 《大地测量与地球动力学》 * |
辜声峰: "多频GNSS非差非组合精密数据处理理论及其应用", 《中国博士学位论文全文数据库 基础科学辑》 * |
Cited By (52)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111323748A (en) * | 2018-12-13 | 2020-06-23 | 千寻位置网络有限公司 | Differential positioning method and system |
CN111323748B (en) * | 2018-12-13 | 2023-07-21 | 千寻位置网络有限公司 | Differential positioning method and system |
CN109856656A (en) * | 2019-02-14 | 2019-06-07 | 上海华测导航技术股份有限公司 | A kind of navigation locating method, device, electronic equipment and storage medium |
CN111831965A (en) * | 2019-04-23 | 2020-10-27 | 千寻位置网络有限公司 | Integrity monitoring method and device for UPD correction number in SSR |
CN111831965B (en) * | 2019-04-23 | 2023-09-12 | 千寻位置网络有限公司 | Integrity monitoring method and device for UPD correction in SSR |
CN110221320A (en) * | 2019-05-07 | 2019-09-10 | 中国科学院国家授时中心 | A kind of Beidou inter-frequency deviation measuring method based on parabola antenna observation |
CN110221320B (en) * | 2019-05-07 | 2023-01-17 | 中国科学院国家授时中心 | Beidou inter-frequency deviation measurement method based on parabolic antenna observation |
CN110018507A (en) * | 2019-05-08 | 2019-07-16 | 中国科学院国家授时中心 | It is a kind of based on make between constellation difference combination accurate one-point positioning method and system |
CN110018507B (en) * | 2019-05-08 | 2020-11-20 | 中国科学院国家授时中心 | Combined precise point positioning method and system based on constellation intercropping difference |
CN110007328A (en) * | 2019-05-10 | 2019-07-12 | 国网浙江省电力有限公司信息通信分公司 | Non-combined RTK localization method based on four frequency observation of No. three satellites of Beidou |
CN110286396A (en) * | 2019-05-20 | 2019-09-27 | 苏州时空复弦网络科技有限公司 | A kind of non-combined PPP method of non-difference based on ionosphere delay prior information and the double constraints of change in time and space information |
CN110780588A (en) * | 2019-10-16 | 2020-02-11 | 北京航空航天大学 | Wide-area accurate time service WPT system and method |
CN110780588B (en) * | 2019-10-16 | 2020-08-11 | 北京航空航天大学 | Wide-area accurate time service WPT system and method |
WO2021174482A1 (en) * | 2020-03-05 | 2021-09-10 | 上海华测导航技术股份有限公司 | Method and device for converting state space reputation into observation space reputation |
CN111367160A (en) * | 2020-03-30 | 2020-07-03 | 武汉大学 | GNSS time service receiver precision clock regulation and control method and system |
CN111447030A (en) * | 2020-03-31 | 2020-07-24 | 中国科学院地质与地球物理研究所 | Network time service precision control method and system of earthquake collector and terminal equipment |
CN111766616A (en) * | 2020-06-15 | 2020-10-13 | 中国人民解放军61081部队 | Beidou second-order time transfer satellite-side multipath error correction method |
CN112859120A (en) * | 2021-01-11 | 2021-05-28 | 中国科学院国家授时中心 | Continuous GNSS carrier phase time and frequency transfer method |
CN113099381A (en) * | 2021-04-06 | 2021-07-09 | 苏州迭慧智能科技有限公司 | Antenna industrial parameter intelligent sensor and multifunctional intelligent sensing network |
CN113325446A (en) * | 2021-06-01 | 2021-08-31 | 中国科学院国家授时中心 | Multi-mode common-frequency GNSS carrier phase time transfer method and system |
CN113325446B (en) * | 2021-06-01 | 2023-08-08 | 中国科学院国家授时中心 | Multimode common-frequency GNSS carrier phase time transfer method and system |
CN113608248A (en) * | 2021-06-25 | 2021-11-05 | 北京建筑大学 | Beidou 5G fused high-precision routing inspection personnel positioning method and related equipment |
CN113608248B (en) * | 2021-06-25 | 2023-06-13 | 北京建筑大学 | Beidou 5G fusion high-precision patrol personnel positioning method and related equipment |
CN113671545A (en) * | 2021-08-27 | 2021-11-19 | 中国科学院国家授时中心 | Satellite-ground precise time synchronization and carrier phase cycle slip detection method |
CN114002723A (en) * | 2021-09-14 | 2022-02-01 | 中国科学院国家授时中心 | Beidou GEO satellite precise time service system and method based on directional antenna |
CN113568020A (en) * | 2021-09-27 | 2021-10-29 | 长沙学院 | Satellite navigation positioning error correction method and device considering hardware inter-frequency difference |
CN114355419A (en) * | 2021-12-15 | 2022-04-15 | 中国科学院国家授时中心 | Distributed Beidou position service center RTK product positioning method and positioning device |
CN114280644A (en) * | 2021-12-22 | 2022-04-05 | 广州南方卫星导航仪器有限公司 | PPP-B2B service-based precise point positioning system and method |
CN114397683A (en) * | 2021-12-30 | 2022-04-26 | 深圳市华芯云创科技有限公司 | Beidou non-differential non-combination PPP-RTK positioning method |
CN114397683B (en) * | 2021-12-30 | 2023-01-13 | 深圳市华芯云创科技有限公司 | Beidou non-differential non-combination PPP-RTK positioning method |
CN114911153A (en) * | 2022-03-31 | 2022-08-16 | 中国科学院国家授时中心 | Universal time UT1 determination method based on satellite-ground ultra-long baseline interferometry |
CN114911153B (en) * | 2022-03-31 | 2023-02-21 | 中国科学院国家授时中心 | Universal time UT1 determination method based on satellite-ground overlength baseline interferometry |
CN114995098A (en) * | 2022-04-22 | 2022-09-02 | 中国科学院国家授时中心 | Satellite interference time transfer method based on picosecond-level precision of phase delay measurement |
CN114740513A (en) * | 2022-05-20 | 2022-07-12 | 航天宏图信息技术股份有限公司 | GNSS positioning method and device |
CN115201870A (en) * | 2022-07-04 | 2022-10-18 | 中国科学院精密测量科学与技术创新研究院 | Multi-frequency multi-mode GNSS non-differential non-combination time transfer method with prior constraint |
CN115453851A (en) * | 2022-09-01 | 2022-12-09 | 广东电网有限责任公司 | BD 3-based multi-satellite comprehensive time service method, system, device and medium |
CN115453851B (en) * | 2022-09-01 | 2024-02-13 | 广东电网有限责任公司 | BD 3-based multi-satellite comprehensive time service method, system, equipment and medium |
CN115407371A (en) * | 2022-09-02 | 2022-11-29 | 中国人民解放军国防科技大学 | PPP-B2B-based real-time high-precision time transfer method and device |
CN115407371B (en) * | 2022-09-02 | 2023-08-15 | 中国人民解放军国防科技大学 | PPP-B2B-based real-time high-precision time transmission method and device |
CN115225245A (en) * | 2022-09-20 | 2022-10-21 | 中国科学院国家授时中心 | Non-differential non-combination PPP-RTK long baseline real-time transfer method |
CN116184441B (en) * | 2022-12-20 | 2023-10-20 | 北京航空航天大学 | Precise atmospheric correction enhanced intelligent terminal PPP instantaneous sub-meter positioning method |
CN116184441A (en) * | 2022-12-20 | 2023-05-30 | 北京航空航天大学 | Precise atmospheric correction enhanced intelligent terminal PPP instantaneous sub-meter positioning method |
CN116243591B (en) * | 2023-01-28 | 2023-09-29 | 北京航空航天大学 | Subnanosecond time service method integrating UTC (k) and Beidou broadcast ephemeris |
CN116243591A (en) * | 2023-01-28 | 2023-06-09 | 北京航空航天大学 | Subnanosecond time service method integrating UTC (k) and Beidou broadcast ephemeris |
CN115856945A (en) * | 2023-02-03 | 2023-03-28 | 北京航空航天大学 | Precise single-point time service method, device and medium suitable for GNSS multi-frequency point |
CN116299618A (en) * | 2023-03-24 | 2023-06-23 | 中国科学院精密测量科学与技术创新研究院 | Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters |
CN116299618B (en) * | 2023-03-24 | 2024-03-19 | 中国科学院精密测量科学与技术创新研究院 | Carrier phase satellite common view time transfer method based on PPP (point-to-point protocol) calculation parameters |
CN116540282B (en) * | 2023-05-08 | 2023-09-29 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
CN116540282A (en) * | 2023-05-08 | 2023-08-04 | 北京交通大学 | Inertial sensor assisted GNSS-based mobile carrier accurate time transfer method |
CN117031507A (en) * | 2023-10-08 | 2023-11-10 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Precise single point positioning method suitable for BDS-3B 1CB2a double-frequency signals |
CN117310770A (en) * | 2023-11-28 | 2023-12-29 | 武汉理工大学三亚科教创新园 | Sea level height inversion method based on Beidou ionosphere-free combined precise single-point positioning |
CN117310770B (en) * | 2023-11-28 | 2024-02-13 | 武汉理工大学三亚科教创新园 | Sea level height inversion method based on Beidou ionosphere-free combined precise single-point positioning |
Also Published As
Publication number | Publication date |
---|---|
CN108919634B (en) | 2020-08-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108919634A (en) | A kind of three non-non-combined observation Time Transmission system and method for difference of frequency of Beidou | |
CN107367716B (en) | A kind of high-precision satellite-borne SAR geometric calibration method | |
CN107153209B (en) | A kind of low rail navigation satellite real-time accurate orbit determination method of short arc segments | |
US10078140B2 (en) | Navigation satellite system positioning involving the generation of advanced correction information | |
Chen et al. | Real-time precise point positioning using single frequency data | |
RU2565386C2 (en) | Method, apparatus and system for determining position of object, having global navigation satellite system receiver, by processing non-differential data, similar to carrier phase measurements, and external data similar to ionospheric data | |
CN106970398B (en) | Satellite visibility analysis and ephemeris forecasting method considering satellite shielding condition | |
CN108196272A (en) | A kind of satellite navigation positioning device and method based on real-time accurate One-Point Location | |
CN106324629A (en) | BDS_GPS_GLONASS precise point positioning integration method | |
CN110174104A (en) | A kind of Combinated navigation method, device, electronic equipment and readable storage medium storing program for executing | |
CN113253314A (en) | Time synchronization method and system between low-earth-orbit satellites | |
Ning et al. | Single-frequency precise point positioning enhanced with multi-GNSS observations and global ionosphere maps | |
CN110146904B (en) | Accurate modeling method suitable for regional ionized layer TEC | |
WO2016141030A1 (en) | Gnss cooperative receiver system | |
Li et al. | Precise orbit determination for the Haiyang-2D satellite using new onboard BDS-3 B1C/B2a signal measurements | |
JP3598372B2 (en) | Measurement method and apparatus for positioning satellite receiver bias | |
Wielgosz et al. | Research on GNSS positioning and applications in Poland in 2015–2018 | |
Zhang et al. | Study on PPP time comparison based on BeiDou-3 new signal | |
Ji et al. | First preliminary fast static ambiguity resolution results of medium-baseline with triple-frequency BeiDou wavebands | |
Liu | Positioning performance of single-frequency GNSS receiver using Australian regional ionospheric corrections | |
Caojun et al. | BeiDou-GPS integrated dual-system with multi-satellites for positioning and navigating farm vehicles | |
Din12 et al. | Monitoring vertical land motion in Malaysia using global positioning system (GPS) | |
Nishimoto et al. | Determination of the parameters of a Japanese-regional Klobuchar ionospheric model based on GR models and SCH analysis by GEONET data | |
Liu et al. | Mapping the regional ionospheric TEC using a spherical cap harmonic model and IGS products in high latitudes and the arctic region | |
Vankadara et al. | Performance Analysis of Various Ionospheric Delay Corrections in Single-frequency GPS Positioning solution at a low latitude Indian location |
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 |