CN109743776A - A kind of base station networking method for synchronizing time based on GNSS - Google Patents
A kind of base station networking method for synchronizing time based on GNSS Download PDFInfo
- Publication number
- CN109743776A CN109743776A CN201811422469.9A CN201811422469A CN109743776A CN 109743776 A CN109743776 A CN 109743776A CN 201811422469 A CN201811422469 A CN 201811422469A CN 109743776 A CN109743776 A CN 109743776A
- Authority
- CN
- China
- Prior art keywords
- satellite
- time
- station
- base station
- value
- 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
Landscapes
- Synchronisation In Digital Transmission Systems (AREA)
- Mobile Radio Communication Systems (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The present invention provides a kind of base station networking method for synchronizing time based on GNSS, can be realized the time synchronization of large coverage, high synchronization accuracy.Base station networking method for synchronizing time based on GNSS of the invention, compared with tradition is based on RTK method for synchronizing time, construction base station is not needed, it will be able to realize that the chronometer time of each website nanosecond in any region 20km is synchronous, has the characteristics that non-stop layer website, survival ability is strong.Base station networking method for synchronizing time based on GNSS of the invention is by way of the interactive information of each website networking and data interaction, it can not only make have the high-precision time synchronization of region website holding for possessing total view satellite condition originally, and high-precision time synchronization is also able to achieve for the region of no identical observation satellite or the only identical observation satellite in part, not only the high-precise synchronization between adjoining website had been realized, but also has realized the continuous covering of whole network.
Description
Technical field
The invention belongs to Time synchronization technique fields, and in particular to a kind of base station networking time synchronization side based on GNSS
Method.
Background technique
High-precision time synchronization (two base station times within any 10 kilometers are required in current many technical fields
Synchronization accuracy is less than or equal to 3ns), such as electric power networks, mobile communications network system, very big region is covered, is often required that
Time synchronization is realized in entire area.Meanwhile in such systems, the synchronization accuracy between mutually adjacent website closes very much again
Key.For example, the base station of mobile communications network will realize the functions such as more base station collaborations, Internet of Things positioning, it is desirable to mobile network's letter
Timing tracking accuracy number mutually between overlapping base station was better than for 10 nanoseconds, even higher.Global Navigation Satellite System (GNSS)
With round-the-clock, Global coverage, high-precision characteristic, thus become the important means of time synchronization between realization multistation.
Mainly there are satellite common vision, RTK, PPP etc. currently based on the high-precision Time synchronization technique of GNSS.Satellite common vision skill
Art can be realized wide area covering, but its to base station underlying hardware facility requirements height, generally will website arrange atomic clock, cost compared with
Greatly；The principle for the Time Synchronization Network that RTK is realized is mainly based upon central station to navigation signal transmission path error and and satellite
Amendment in relation to error, coverage area are localities, and the cost for establishing central station is often also higher；Based on PPP technology
Time synchronization coverage area is big between the station of realization, and precision is high, but needs to rely on third-party enhancing data, and operating cost is high,
Reliability is uncontrolled.Therefore, it is badly in need of a kind of low cost, the time synchronization of high reliablity of a kind of covering of wide area, high-precise synchronization
Method.
Summary of the invention
In view of this, the present invention provides a kind of base station networking method for synchronizing time based on GNSS, can be realized and cover greatly
The time synchronization of lid range, high synchronization accuracy.
To achieve the above object, a kind of base station networking method for synchronizing time based on GNSS of the invention, including walk as follows
It is rapid:
Step 1, for each base station, single star clock bias estimation value average to the time of all satellites in the base station carries out satellite weighting
It is average, obtain the base station initial time adjustment amount；The list star clock bias estimation value is current time satellite pseudorange and geometric distance
Difference, the time average single star clock bias estimation value is the average value of each moment list star clock bias estimation value；
Step 2, i-th satellite observed for our station, (i=1,2,3......N), N is observed by our station website
The population of satellite, calculate the time deviation of the satellite of other base stations communicated with our station；Wherein, this is defended under j-th of base station
The time deviation at star current time is the single star clock bias estimation value and j-th of base station initial time tune of the satellite under j-th of base station
The difference of whole amount, wherein j=1,2,3,4 ... n, n are the sum at every other station received by our station；
Average statistical under each base station that acquisition and our station communicate in present period to this satellite time deviation utilizes
Each average statistical obtains the time deviation correction value for i-th satellite that present period our station is observed；
Using the time deviation correction value for i-th satellite that present period and upper period our station are observed, our station is constructed
The final time deviation correction value for i-th satellite observed；
The optimal time in our station website region is calculated according to the final time deviation correction value of all satellites of our station
Adjustment amount correction value；
Step 3, our station initial time adjustment amount is added to obtain our station with the optimal time adjustment amount correction value final
Time error Δ T；
Step 4, time error Δ T is sent into time error and corrects circuit, correct local zone time, output sync pulse signal
And reference frequency signal.
Wherein, the GNSS satellite signal observed quantity is original observed quantity or statistical information.
Wherein, each out by the pseudo range observed quantity of GNSS satellite signal and the positional information calculation of base station in the step 1
Single star clock bias estimation value of a satellite；In the satellite weighted average, satellite weighted value has with corresponding satellite elevation angle and signal-to-noise ratio
It closes.
Wherein, the step 2 includes following sub-step:
Step 2.1, i-th satellite time deviation of j-th of base station of current time is calculated:
Wherein, Δ TjFor the initial time adjustment amount of j-th of base station,For single star of i-th of satellite of j-th of base station
Clock bias estimation value；
J-th of base station is obtained to the average statistical and variance of all time deviations under i-th satellite present period；
Step 2.2, it calculates present period our station and observes i-th satellite time drift correction value
Wherein, ΩkFor the set at every other station received by our station；Mj iIt is current to i-th satellite for j-th of base station
The average statistical of time deviation under period；Uj iFor Mj iWeighting coefficient, the weighting coefficient by with a distance from this website, observation
Number of satellite and the variance of number and time deviation is observed to determine；
Step 2.3, i-th satellite time drift correction value was observed according to the our station that a upper period calculatesWith work as
Preceding period our station observes i-th satellite time drift correction valueThe final time deviation of construction i-th satellite of our station is repaired
Positive value Fi:
Wherein, v is Dynamic gene, value range 0~1；
Step 2.4, the final time deviation correction value F of i-th satellite of our station resolved according to step 2.3iCalculate our station
Optimal time adjustment amount correction value Δ T in point regioncorrection:
Wherein, WiFor the weighted value of i-th satellite.
The utility model has the advantages that
1. the base station networking method for synchronizing time of the invention based on GNSS is based on RTK method for synchronizing time phase with tradition
Than not needing construction base station, it will be able to realize that the chronometer time of each website nanosecond in any region 20km is synchronous, there is nothing
The characteristics of central site, survival ability are strong.
2. the friendship that the base station networking method for synchronizing time of the invention based on GNSS passes through each website networking and data interaction
The mode of mutual information, can not only make to have the region website for possessing total view satellite condition originally and kept for high-precision time same
Step, and for the region of no identical observation satellite or the only identical observation satellite in part, (as shown in Figure 3 two face domain station
" north and south wall " phenomenon figure between point) it is also able to achieve high-precision time synchronization, both realize the height between adjoining website
Accurate synchronization, and realize the continuous covering of whole network.
3. compared with satellite common vision technology, the base station networking method for synchronizing time of the invention based on GNSS does not depend on height
Steady atomic clock, only need to configure general crystal oscillator for time synchronism equipment can realize high-precision time synchronization, and between website
Data interaction link can not depend on wired optical fiber networking, by means of current mobile communications network or data radio station broadcast just
It may be implemented, therefore construction cost is low.
4. the base station networking method for synchronizing time of the invention based on GNSS only needs to utilize compared with traditional PPP technology
The moonscope information of other base stations of region does not need to enhance data by means of third party, split-second precision between station can be realized
It is synchronous, have and relies on infrastructure low, and independently controllable advantage.
Detailed description of the invention
Fig. 1 is the precise synchronization networking topological diagram of the invention based on GNSS.
Fig. 2 is that GNSS signal reception device single machine information handles frame diagram in the present invention.
Fig. 3 is two " north and south wall " phenomenon figures faced between the website of domain.
Specific embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
The present invention provides a kind of neighborhood high-precision time synchronization method based on GNSS, the method achieve go to center
Change, realizes that the precise synchronization of region multistation, each website are not only available by each website networking and data interaction
Our station GNSS satellite signal observed quantity information, the GNSS satellite signal that other base stations can also be received by relevant interface agreement are seen
Metrical information, and the GNSS satellite signal observed quantity information of our station is gone out by relevant interface agreement broadcast, realize big covering
The time synchronization of range.Precise synchronization networking topological diagram based on GNSS is as shown in Figure 1.
Wherein, time error statistical disposition strategy thinking is as follows:
Our station initial time adjustment amount is calculated according to our station GNSS satellite signal observed quantity first, then according to other bases
GNSS satellite of standing signal observed quantity calculates the optimal time adjustment amount correction value in this website region, our station initial time
Adjustment amount is final time error Δ T plus time adjustment amount correction value.
Specifically, GNSS signal reception device single machine information processing frame in, GNSS signal reception device single machine when
Between in error statistics processing module, set a base station as our station, the pseudo range observed quantity using our station GNSS satellite signal is (original
Observed quantity or statistical information) and the positional information calculation of our station go out the pseudorange of each satellite of our station and the difference of geometric distance, i.e., singly
Star clock bias estimation value；Our station initial time tune is obtained according to single star clock bias estimation value of all satellites of our station and corresponding weighted information
Whole amount；Obtain the initial time adjustment amount of each base station；
Step 1, for each base station, single star clock bias estimation value average to the time of all satellites in the base station carries out satellite weighting
It is average, obtain the base station initial time adjustment amount；The list star clock bias estimation value is current time satellite pseudorange and geometric distance
Difference, the time average single star clock bias estimation value is the average value of each moment list star clock bias estimation value；
Step 2, i-th satellite observed for our station, (i=1,2,3......N), N is observed by our station website
The population of satellite, calculate the time deviation of the satellite of other base stations communicated with our station；Wherein, this is defended under j-th of base station
The time deviation at star current time is the single star clock bias estimation value and j-th of base station initial time tune of the satellite under j-th of base station
The difference of whole amount, wherein j=1,2,3,4 ... n, n are the sum at every other station received by our station；
Average statistical under each base station that acquisition and our station communicate in present period to this satellite time deviation utilizes
Each average statistical obtains the time deviation correction value for i-th satellite that present period our station is observed；
Using the time deviation correction value for i-th satellite that present period and upper period our station are observed, our station is constructed
The final time deviation correction value for i-th satellite observed；
The optimal time in our station website region is calculated according to the final time deviation correction value of all satellites of our station
Adjustment amount correction value；
Step 3, our station initial time adjustment amount is added to obtain our station with the optimal time adjustment amount correction value final
Time error Δ T；
Step 4, time error Δ T is sent into time error and corrects circuit, correct local zone time, output sync pulse signal
And reference frequency signal, such as PPS time synchronization pulse signal, 10MHz reference frequency signal, temporal information.
In time error statistical disposition module, different time error statistical disposition strategies can be formulated to obtain our station
Final time error Δ T, the specific implementation steps are as follows:
Step 1, single star clock bias estimation value of i-th satellite is calculated:
Wherein c is the light velocity, (xi,yi,zi) it is i-th coordinate of the satellite position resolved by navigation message, (x0,y0,z0)
For the position coordinates of the base station, ρiFor the pseudo range observed quantity of i-th satellite, (i=1,2,3......N), N is our station website institute
The population of satellite observed；
Our station initial time adjustment amount is calculated according to single star clock bias estimation value of N number of satellite and corresponding weighted information are as follows:
Wherein, WiFor the weighted value of i-th satellite, the satellite that the usual elevation angle is high, signal-to-noise ratio is high gives bigger weighted value.
Step 2, the optimal time in this website region is calculated according to the GNSS satellite signal observed quantity of other websites
Adjustment amount correction value, including following sub-step:
Step 2.1, i-th satellite time deviation of j-th of base station of current time is calculated:
Wherein, Δ TjFor the initial time adjustment amount of j-th of base station,For single star of i-th of satellite of j-th of base station
Clock bias estimation value, wherein j=1,2,3,4 ... n, n are the sum at every other station received by our station；
J-th of base station is obtained to the average statistical and variance of all time deviations under i-th satellite present period；
Step 2.2, it calculates present period our station and observes i-th satellite time drift correction value
Wherein, ΩkFor the set at every other station received by our station；Mj iIt is current to i-th satellite for j-th of base station
The average statistical of time deviation under period；Uj iFor Mj iWeighting coefficient, the weighting coefficient by with a distance from this website, observation
Number of satellite and the variance of number and time deviation is observed to determine；
Step 2.3, i-th satellite time drift correction value was observed according to the our station that a upper period calculatesWith work as
Preceding period our station observes i-th satellite time drift correction valueThe final time deviation of construction i-th satellite of our station is repaired
Positive value Fi:
Wherein, v is Dynamic gene, value range 0~1.
Step 2.4, the final time deviation correction value F of i-th satellite of our station resolved according to step 2.3iCalculate our station
Optimal time adjustment amount correction value Δ T in point regioncorrection, calculation method are as follows:
Step 3, the final time error Δ T of our station is calculated:
Δ T=Δ Torigin+ΔTcorrection (6)
In conclusion the above is merely preferred embodiments of the present invention, being not intended to limit the scope of the present invention.
All within the spirits and principles of the present invention, any modification, equivalent replacement, improvement and so on should be included in of the invention
Within protection scope.
Claims (4)
1. a kind of base station networking method for synchronizing time based on GNSS, which comprises the steps of:
Step 1, for each base station, it is flat that single star clock bias estimation value average to the time of all satellites in the base station carries out satellite weighting
, the base station initial time adjustment amount is obtained；The list star clock bias estimation value be current time satellite pseudorange and geometric distance it
Difference, the time average single star clock bias estimation value are the average value of each moment list star clock bias estimation value；
Step 2, i-th satellite observed for our station, (i=1,2,3......N), N is defended by what our station website observed
Star sum, calculates the time deviation of the satellite of other base stations communicated with our station；Wherein, the satellite is worked as under j-th of base station
The time deviation at preceding moment is the single star clock bias estimation value and j-th of base station initial time adjustment amount of the satellite under j-th of base station
Difference, wherein j=1,2,3,4 ... n, n are the sum at every other station received by our station；
Average statistical under each base station that acquisition and our station communicate in present period to this satellite time deviation, utilizes each system
Meter mean value obtains the time deviation correction value for i-th satellite that present period our station is observed；
Using the time deviation correction value for i-th satellite that present period and upper period our station are observed, constructs our station and seen
The final time deviation correction value for i-th satellite surveyed；
The optimal time adjustment in our station website region is calculated according to the final time deviation correction value of all satellites of our station
Measure correction value；
Step 3, our station initial time adjustment amount is added to obtain the our station final time with the optimal time adjustment amount correction value
Error delta T；
Step 4, time error Δ T is sent into time error and corrects circuit, correct local zone time, output sync pulse signal and ginseng
Examine frequency signal.
2. a kind of base station networking method for synchronizing time based on GNSS as described in claim 1, which is characterized in that the GNSS
Satellite-signal observed quantity is original observed quantity or statistical information.
3. a kind of base station networking method for synchronizing time based on GNSS as described in claim 1, which is characterized in that the step
In 1, single star clock bias estimation of each satellite is gone out by the pseudo range observed quantity of GNSS satellite signal and the positional information calculation of base station
Value；In the satellite weighted average, satellite weighted value is related with corresponding satellite elevation angle and signal-to-noise ratio.
4. a kind of base station networking method for synchronizing time based on GNSS as described in claim 1, which is characterized in that the step
2 include following sub-step:
Step 2.1, i-th satellite time deviation of j-th of base station of current time is calculated:
Wherein, Δ TjFor the initial time adjustment amount of j-th of base station,For single star clock deviation of i-th of satellite of j-th of base station
Estimated value；
J-th of base station is obtained to the average statistical and variance of all time deviations under i-th satellite present period；
Step 2.2, it calculates present period our station and observes i-th satellite time drift correction value
Wherein, ΩkFor the set at every other station received by our station；Mj iIt is j-th of base station to i-th satellite present period
Under time deviation average statistical；Uj iFor Mj iWeighting coefficient, the weighting coefficient is by with a distance from this website, observation satellite
Number and the variance of number and time deviation is observed to determine；
Step 2.3, i-th satellite time drift correction value was observed according to the our station that a upper period calculatesAnd present period
Our station observes i-th satellite time drift correction valueConstruct the final time deviation correction value F of i-th satellite of our stationi:
Wherein, v is Dynamic gene, value range 0~1；
Step 2.4, the final time deviation correction value F of i-th satellite of our station resolved according to step 2.3iCalculate this website place
Optimal time adjustment amount correction value Δ T in regioncorrection:
Wherein, WiFor the weighted value of i-th satellite.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811422469.9A CN109743776B (en) | 2018-11-27 | 2018-11-27 | Base station networking time synchronization method based on GNSS |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811422469.9A CN109743776B (en) | 2018-11-27 | 2018-11-27 | Base station networking time synchronization method based on GNSS |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109743776A true CN109743776A (en) | 2019-05-10 |
| CN109743776B CN109743776B (en) | 2020-12-15 |
Family
ID=66358257
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811422469.9A Active CN109743776B (en) | 2018-11-27 | 2018-11-27 | Base station networking time synchronization method based on GNSS |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109743776B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110505685A (en) * | 2019-08-29 | 2019-11-26 | 成都精位科技有限公司 | Base station calibration method and device based on global position system |
| WO2021253916A1 (en) * | 2020-06-16 | 2021-12-23 | 华为技术有限公司 | Data processing method and related device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007026098A1 (en) * | 2005-08-30 | 2007-03-08 | Alcatel Lucent | Device for checking the consistency of broadcast information, with no constraints on the clocks of the monitoring stations, for a satellite navigation system |
| CN101881817A (en) * | 2009-05-06 | 2010-11-10 | 中国科学院微电子研究所 | Device and method for detecting feeble signals of global navigation satellite system |
| CN102279404A (en) * | 2010-06-13 | 2011-12-14 | 上海伽利略导航有限公司 | Seamless positioning method and device |
| CN103048671A (en) * | 2011-11-16 | 2013-04-17 | 郭杭 | New RTD (Relative Time Difference), MEMS-IMU (Micro-electromechanical Systems-Inertial Measurement Unit) and ZigBee combined high-precision seamless navigation method for global navigation satellite system |
| CN104181550A (en) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | Common view time and frequency transmitting method based on BeiDou navigation satellite system |
| CN106292265A (en) * | 2016-09-09 | 2017-01-04 | 北京航空航天大学 | A kind of many ground based on aeronautical satellite method for synchronizing time |
-
2018
- 2018-11-27 CN CN201811422469.9A patent/CN109743776B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007026098A1 (en) * | 2005-08-30 | 2007-03-08 | Alcatel Lucent | Device for checking the consistency of broadcast information, with no constraints on the clocks of the monitoring stations, for a satellite navigation system |
| CN101881817A (en) * | 2009-05-06 | 2010-11-10 | 中国科学院微电子研究所 | Device and method for detecting feeble signals of global navigation satellite system |
| CN102279404A (en) * | 2010-06-13 | 2011-12-14 | 上海伽利略导航有限公司 | Seamless positioning method and device |
| CN103048671A (en) * | 2011-11-16 | 2013-04-17 | 郭杭 | New RTD (Relative Time Difference), MEMS-IMU (Micro-electromechanical Systems-Inertial Measurement Unit) and ZigBee combined high-precision seamless navigation method for global navigation satellite system |
| CN104181550A (en) * | 2014-08-19 | 2014-12-03 | 北京无线电计量测试研究所 | Common view time and frequency transmitting method based on BeiDou navigation satellite system |
| CN106292265A (en) * | 2016-09-09 | 2017-01-04 | 北京航空航天大学 | A kind of many ground based on aeronautical satellite method for synchronizing time |
Non-Patent Citations (1)
| Title |
|---|
| 李东新: "伪卫星共视法时间同步模型及其精度分析", 《导航定位与授时》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110505685A (en) * | 2019-08-29 | 2019-11-26 | 成都精位科技有限公司 | Base station calibration method and device based on global position system |
| CN110505685B (en) * | 2019-08-29 | 2021-01-15 | 成都精位科技有限公司 | Base station calibration method and device based on satellite positioning system |
| WO2021253916A1 (en) * | 2020-06-16 | 2021-12-23 | 华为技术有限公司 | Data processing method and related device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109743776B (en) | 2020-12-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1709460B1 (en) | Transfer of calibrated time information in a mobile terminal | |
| JP5791968B2 (en) | Method and apparatus for determining position using a small number of GPS satellites and synchronous and asynchronous base stations | |
| US8837530B2 (en) | Method and system for adaptive synchronization of timing information generated by independently clocked communication nodes | |
| CN105867108A (en) | National standard time remote recurrence method | |
| CN102023290A (en) | High-precision distributed pulse signal time difference of arrival detection system | |
| JP2008547340A (en) | Method and apparatus for determining whether a mobile terminal has moved outside a predetermined location | |
| TW200807008A (en) | Associating a universal time with a received signal | |
| US20020167934A1 (en) | Method and system for timebase synchronization | |
| CN108882356A (en) | The synchronous method of clock, time reference source device and clock reappear equipment | |
| CN106292264A (en) | A kind of wide area accurate time transmission system strengthening system based on GNSS high accuracy | |
| CN105717784B (en) | Electric power time service method | |
| CN109743776A (en) | A kind of base station networking method for synchronizing time based on GNSS | |
| Liu et al. | Correction of fractional cycle bias of pseudolite system for user integer ambiguity resolution | |
| Liu et al. | Time synchronization in communication networks based on the Beidou foundation enhancement system | |
| CN201945685U (en) | High-accuracy time difference of arrival (TDOA) measuring system for distribution type pulse signals | |
| CN111366960A (en) | Navigation positioning method based on communication time service | |
| Lei et al. | Hybsync: nanosecond wireless position and clock synchronization based on UWB communication with multisensors | |
| EP1323317A4 (en) | Distributed location system | |
| KR101099175B1 (en) | Transfer of calibrated time information in a mobile terminal | |
| CN109728868A (en) | A kind of GNSS base station networking method for synchronizing time examined based on multiple integrity | |
| CN110286580A (en) | A kind of time service method of optical fiber and the punctual synchronous seamless interfacing of wireless time service | |
| CN109412692A (en) | A kind of distribution type fiber-optic time dissemination system and its sagnac effect compensation method | |
| CN112034697B (en) | High-precision RDSS time service method | |
| CN113075702B (en) | Unmanned aerial vehicle defense high-precision time service method and system based on pseudolite technology | |
| CN106886035A (en) | A kind of signal synchronizing method and device based on adaptive loop circuit correction |
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 |