CN104181550B - A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system - Google Patents

A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system Download PDF

Info

Publication number
CN104181550B
CN104181550B CN201410409512.3A CN201410409512A CN104181550B CN 104181550 B CN104181550 B CN 104181550B CN 201410409512 A CN201410409512 A CN 201410409512A CN 104181550 B CN104181550 B CN 104181550B
Authority
CN
China
Prior art keywords
time
satellite
difference
data
station
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.)
Active
Application number
CN201410409512.3A
Other languages
Chinese (zh)
Other versions
CN104181550A (en
Inventor
王宏博
张升康
易航
王海峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Radio Metrology and Measurement
Original Assignee
Beijing Institute of Radio Metrology and Measurement
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Radio Metrology and Measurement filed Critical Beijing Institute of Radio Metrology and Measurement
Priority to CN201410409512.3A priority Critical patent/CN104181550B/en
Publication of CN104181550A publication Critical patent/CN104181550A/en
Application granted granted Critical
Publication of CN104181550B publication Critical patent/CN104181550B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • GPHYSICS
    • G04HOROLOGY
    • G04RRADIO-CONTROLLED TIME-PIECES
    • G04R20/00Setting the time according to the time information carried or implied by the radio signal
    • G04R20/02Setting 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
    • G04R20/04Tuning or receiving; Circuits therefor

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electric Clocks (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

The invention discloses a kind of Beidou satellite navigation system common-view time Frequency Transfer method, the method extracts the time difference of local zone time and each tracking satellite using the measurement of timing type Beidou receiver, clock correction parameter in navigation message, ionosphere time delay and troposphere delay parameter and height angular data, calculate the clock correction of local zone time and tracked big-dipper satellite star clock;Through the pretreatment of initial data, and the elevation angle weighting of each tracking satellite, calculate the difference of local zone time and dipper system time.With Big Dipper system time as intermediary, two comparison results for regarding station altogether are exchanged by the Big Dipper short message communication function of No., the time difference at two stations can be obtained in real time, and then calculate the relative frequency difference at two stations.Technical scheme of the present invention can overcome the dependence using GPS passing time, realize accurately Time Transmission by Beidou satellite navigation system;This method do not need high cost hardware supported, it is easy to operate, can long-time automatic running, greatly reduce labour and use cost.

Description

A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system
Technical field
The present invention relates to satellite navigation application, more particularly to a kind of common apparent time based on Beidou satellite navigation system Between Frequency Transfer method.
Background technology
Time, as one of seven fundamental physical quantities, is remotely propagated have the characteristics that passing through means of electromagnetic waves.High-precision Degree Time Transmission is the important research direction of time-frequency metering field.Existing remote time Frequency Transfer mode mainly has GPS to regard altogether Method, gps carrier phase method and satellite two dimensional method.
GPS regards method altogether with gps satellite as medium, in two websites for carrying out Time Transmission using GPS collection ratio To data, afterwards by two places data exchange processing, Time Transmission result is obtained.This method is that countries in the world time-frequency laboratory should With a kind of remote time transmission method that popularizes the most, low cost, data processing technique are simply ripe, but real-time is poor.
Gps carrier phase time TRANSFER METHOD is also required to two websites for carrying out Time Transmission by measurement type GPS The gps carrier phase data of collection certain time length, using issued after IGS (International GNSS Service) 15 days Precise ephemeris and star clock data, and accurately antenna phase center coordinate, carry out the Data Post of complexity, when obtaining each station Between time difference with IGST, so as to complete the Time Transmission between standing.This method data processing algorithm is complicated, and needs to wait IGS issues precise ephemeris, therefore cannot also complete real time data transmission.
Satellite two-way time transfer method is stood in carry out Time Transmission two typically using geostationary orbit telecommunication satellite Point carries out the long-range transmission of time by ground satellite station and Two Way Satellite Time Transfer to modem.As electromagnetic wave is double To the symmetry of propagation path, main Time transfer receiver error can be eliminated, obtain the Time Transmission result of two places in real time.This Method real-time is good, and index is high, but relatively costly, needs two websites for carrying out Time Transmission to set up ground satellite station, merges Special two-way Time transfer receiver modem is bought, while leased satellite channel.
In sum, GPS cannot all complete real-time Time Transmission depending on method and gps carrier phase method altogether, and satellite is two-way Although Time Transmission method can complete real-time time transmission, equipment complex and expensive, satellite channel lease expenses are very high.
Accordingly, it is desirable to provide a kind of both can guarantee that real-time time transmission quality, and can cost-effective, using simple when Between transmission method.
Content of the invention
The technical problem to be solved in the present invention is to provide a kind of common-view time frequency based on Beidou satellite navigation system to pass Method is passed, to solve the problems, such as real-time time transmission cannot be carried out using GPS in prior art.
For solving above-mentioned technical problem, the present invention adopts following technical proposals.
A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system, the method step include
S1, using measurement type Beidou receiver receive satellite-signal, gather receiver output satellite ephemeris, pseudorange, star Clock parameter and atmospheric propagation delay parameter etc.;
S2, the data obtained according to step S1, try to achieve the data of each Satellite clock and the difference of local zone time, and as former Beginning data;
S3, the initial data utilized in step S2, and linear fit is carried out to data using least square method, obtain each Satellite time is with local zone time in the time difference for comparing period midpoint;
The time of S4, calculating local zone time and Beidou satellite navigation system is in the time difference for comparing period midpoint, and conduct is originally Ground regards result altogether;
Local data exchange is carried out by Beidou satellite communication system between S5, user station, and calculate between two users station Clock correction and relative frequency difference.
Preferably, in step S1, clock correction parameter includes telegraph text data, delayed data and satellite altitude angular data.
Preferably, elimination of rough difference and smothing filtering are carried out to initial data in step S3.
Preferably, in step S4 using the data of each satellite are weighted by satellite species and elevation of satellite Average method obtains the difference of the time of the local zone time of midpoint and Beidou satellite navigation system in comparison time.
Preferably, using formulaCalculate relative weighting, wherein, WPRNiFor No. i-th satellite time difference The relative weighting of data, EPRNiFor No. i-th elevation of satellite.
6th, common-view time Frequency Transfer method according to claim 1, it is characterised in that using formulaAverage relative frequency departure of the two station frequency standards within the τ time period is calculated, wherein ΔTAB(ti) at a time tiMeasure the two station time difference of A, B, Δ TAB(ti+ τ) it is ti+ τ the moment measures the time difference at two station of A, B.
Preferably, the method further includes to accumulate geostationary orbit satellite comparison data sample, and Time transfer receiver is entered Row optimizes.
Beneficial effects of the present invention are as follows:
Technical scheme of the present invention can overcome the dependence using GPS passing time, by Beidou satellite navigation system Realize accurately Time Transmission;This method do not need high cost hardware supported, it is easy to operate, can long-time automatic running, pole Big reduces labour and use cost.
Description of the drawings
Specific embodiment to the present invention is described in further detail below in conjunction with the accompanying drawings;
Fig. 1 illustrates a kind of showing for common-view time Frequency Transfer method based on Beidou satellite navigation system of the present invention It is intended to;
Fig. 2 illustrates the operation principle schematic diagram of the present embodiment.
Specific embodiment
The present invention is described further with reference to one group of embodiment and accompanying drawing.
As shown in figure 1, the invention discloses a kind of based on the concrete of Beidou satellite navigation system common-view time transmission method Step includes:
The first step, the pseudorange of collection timing type Beidou receiver output, ephemeris, star clock parameter and atmospheric propagation delay parameter Etc. data
During Time Transmission is carried out using Beidou navigation satellite, need to select timing type Big Dipper receiving device to receive Big Dipper satellite signal, the PRN numbering of the tracking satellite of collection Beidou receiver output, time, pseudorange, ephemeris orbital parameters, star Clock parameter a0, a1, a2 and t0, ionosphere delay parameter and tropospheric delay parameter.
The initial data of the difference of second step, each tracking satellite star clock of calculating and local zone time
According to all data that extracts in the first step, positioning equation is resolved, local zone time and Big Dipper time can be calculated Difference initial data.
3rd step, each tracking satellite of calculating and local zone time are in the time difference for comparing period midpoint
Initial data to obtaining in second step carries out elimination of rough difference and smothing filtering, carries out Linear Quasi with least square method Close, obtain corresponding to the time difference at comparison time section midpoint.
The time difference of the 4th step, calculating local zone time and dipper system Time transfer receiver period midpoint, result is regarded altogether as local
The data of each tracking satellite are weighted averagely, obtaining in comparison time section by satellite species and elevation of satellite The difference of local zone time and Big Dipper time at point, used as local Big Dipper common-view mode result.The following institute of relative weighting computing formula Show
Wherein, WPRNiRelative weighting for No. i-th satellite time difference data;
EPRNiFor No. i-th elevation of satellite.
Its absolute weight can be obtained by following formula
Pass through Big Dipper short message exchange data between the 5th step, comparison station and calculate clock correction and relative frequency difference
Geostationary geosynchronous satellite in dipper system has short message communication capacity, between common-view mode earth station Can be by its local Big Dipper common-view mode result of exchange.Our station comparison result is sent to distant station by local station through satellite, with When receive the comparison result that sends through satellite of distant station.Clock between our station and distant station is tried to achieve by formula (3) and (4) Difference.
As shown in Fig. 2 A station and B station are respectively two earth stations, two places temporal frequency is carried out depending on method by dipper system altogether Compare and transmission.A station is utilized respectively Research on passive BeiDou navigation position method with B station, calculates the difference at two stations and dipper system time.
Wherein, TAStand the time for A;
TBStand the time for B;
TBDFor the dipper system time;
ΔTAFor A station time and the time deviation of dipper system;
ΔTBFor B station time and the time deviation of dipper system.
Two results obtained in identical epoch or comparison time section are made the difference, you can obtained for two station time differences
ΔTAB=TA-TB
=(TA-TBD)-(TB-TBD) (4)
=Δ TA-ΔTB
If at a time tiThe two station time difference of A, B is measured for Δ TAB(ti), after τ after a while, i.e. ti+ τ the moment measures The time difference at two station of A, B is Δ TAB(ti+ τ), then average relative frequency of the two station frequency standards within the τ time period can be obtained with following formula Deviation:
6th step, accumulation GEO satellite comparison data sample, long-time are compared and obtain optimum results
Due to GEO satellite be continuously visible, therefore respectively stand in tracking GEO satellite during will not occur interrupt feelings Condition.Therefore, it can constantly accumulate the comparison data of GEO satellite receiving channel.Often increase the comparison data of one group of GEO satellite, enter One secondary data of row is fitted.After its data sample is more than 1h, replace the linear fit in the 3rd step by quadratic fit, permissible The common-view time comparison result for optimizing is obtained by the observation of longer time.
Below by one group of specific embodiment, the present invention will be further described:
During Big Dipper common-view mode is carried out, each station that compares is needed to be compared according to the prior comparison time table that works out Right.The period is compared for the ease of each station is unified, each the integral point hour for planning universal time (UTC time) is divided into three sections, Every section of comparison time is set to 20min.So, the comparison period of one day is 72.
Within each Big Dipper common-view mode period of upper step, the time difference data for gathering local zone time and Big Dipper time per second. For the satellite of continuous tracking, 1200 in each time data can make full use of, i.e., 1200 comparison data are one group. For needing to switch the passage of tracking satellite, as receiver acquisition satellite-signal needs the regular hour, accordingly, it would be desirable to will be per Before and after the individual comparison period, the data of 2min are removed, and there are 960 data so each period.
For the data of each period, it is required for calculating the comparison result at period midpoint, corresponding time point is respectively UTC 10min, 30min and 50min hourly.
In sum, the hardware device low cost that this method of the present invention is related to, it is easy to operate, can long-time automatic running, Be easy to promote the use of, and by the Beidou satellite navigation system of China carry out Time Transmission can break away to GPS of America according to Rely.
Obviously, the above embodiment of the present invention is only intended to clearly illustrate example of the present invention, and is not right The restriction of embodiments of the present invention, for those of ordinary skill in the field, may be used also on the basis of the above description To make other changes in different forms, all of embodiment cannot be exhaustive here, every belong to this Obvious change that bright technical scheme is extended out changes row still in protection scope of the present invention.

Claims (4)

1. a kind of common-view time Frequency Transfer method based on Beidou satellite navigation system, it is characterised in that the method step bag Include:
S1, using measurement type Beidou receiver receive satellite-signal, gather receiver output satellite ephemeris, pseudorange, star clock ginseng Number and atmospheric propagation delay parameter;
S2, the data obtained according to step S1, try to achieve the data of each Satellite clock and the difference of local zone time, and as original number According to;
S3, the initial data utilized in step S2, and linear fit is carried out to data using least square method, obtain each satellite Time is with local zone time in the time difference for comparing period midpoint;
The time of S4, calculating local zone time and Beidou satellite navigation system is in the time difference for comparing period midpoint, and conduct is local altogether Depending on result;Average side is weighted using by the data of each satellite by satellite species and elevation of satellite in step S4 Method obtains the difference of the time of the local zone time of midpoint and Beidou satellite navigation system in comparison time, wherein, using formulaCalculate relative weighting, wherein, WPRNiFor the relative weighting of No. i-th satellite time difference data, EPRNiFor i-th Number elevation of satellite;
Local data exchange is carried out by Beidou satellite communication system between S5, user station, and calculate the clock between two users station Difference and relative frequency difference.
2. common-view time Frequency Transfer method according to claim 1, it is characterised in that first to original in step S3 Data carry out elimination of rough difference and the disposal of gentle filter, recycle least square method to carry out linear fit to data.
3. common-view time Frequency Transfer method according to claim 1, it is characterised in that using formulaAverage relative frequency departure of the two station frequency standards within the τ time period is calculated, wherein ΔTAB(ti) at a time tiMeasure the two station time difference of A, B, Δ TAB(ti+ τ) it is ti+ τ the moment measures the time difference at two station of A, B.
4. common-view time Frequency Transfer method according to claim 1, it is characterised in that the method further include S6, Accumulation geostationary orbit satellite comparison data sample, is optimized to Time transfer receiver.
CN201410409512.3A 2014-08-19 2014-08-19 A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system Active CN104181550B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410409512.3A CN104181550B (en) 2014-08-19 2014-08-19 A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410409512.3A CN104181550B (en) 2014-08-19 2014-08-19 A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system

Publications (2)

Publication Number Publication Date
CN104181550A CN104181550A (en) 2014-12-03
CN104181550B true CN104181550B (en) 2017-03-08

Family

ID=51962747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410409512.3A Active CN104181550B (en) 2014-08-19 2014-08-19 A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system

Country Status (1)

Country Link
CN (1) CN104181550B (en)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104702401B (en) * 2015-03-26 2017-03-01 中国人民解放军国防科学技术大学 A kind of synchronous method of ground base navigation system time
CN104950301A (en) * 2015-04-09 2015-09-30 北京寰亚翔宇科技有限责任公司 Exterior trajectory measurement method
CN105158771B (en) * 2015-08-25 2018-01-05 中国科学院国家授时中心 A kind of model parameter production method of tracing to the source in satellite navigation system
CN105044747B (en) * 2015-08-27 2017-09-15 交通信息通信技术研究发展中心 Time synchronization device and method based on multi-satellite common view and filtering
CN105578587B (en) * 2015-12-17 2019-06-04 北京无线电计量测试研究所 A kind of temporal frequency comparison method and equipment
CN106533529B (en) * 2016-09-27 2020-03-03 北京无线电计量测试研究所 Satellite two-way time comparison modulation and demodulation system and method
CN106850177A (en) * 2016-12-15 2017-06-13 北京无线电计量测试研究所 A kind of real-time time transmission method based on Static Precise Point Positioning
CN107643527B (en) * 2017-09-08 2021-01-05 桂林电子科技大学 GPS (global positioning system) common-view signal simulation generation method and GPS common-view signal simulator
CN108196267B (en) * 2017-12-20 2021-06-04 中国科学院国家授时中心 GNSS CP technology-based uninterrupted time transfer method
CN108680930A (en) * 2018-05-17 2018-10-19 中国电子科技集团公司第五十四研究所 A kind of method that embedded remote real time GPS common-view time compares
CN109743776B (en) * 2018-11-27 2020-12-15 中国科学院光电研究院 Base station networking time synchronization method based on GNSS
CN109765582A (en) * 2019-02-13 2019-05-17 重庆市计量质量检测研究院 A kind of temporal frequency calibration system based on GNSS
CN109902882A (en) * 2019-03-21 2019-06-18 北京工业大学 Atomic clock clock deviation prediction model training method and device
CN109991837B (en) * 2019-03-22 2020-10-27 中国科学技术大学 System and method for comparing two-place clocks by using laser common view
CN111030774A (en) * 2019-12-19 2020-04-17 北京无线电计量测试研究所 Real-time common-view data processing method based on Beidou satellite navigation system
CN111552171A (en) * 2020-04-26 2020-08-18 北京无线电计量测试研究所 Atomic frequency standard remote time frequency calibration method, equipment and system
CN112286038A (en) * 2020-07-03 2021-01-29 国网思极神往位置服务(北京)有限公司 Common-view time synchronization method and device for Beidou satellite navigation
CN111884748B (en) * 2020-07-24 2023-02-03 北京卓越航导科技有限责任公司 Coding method applied to real-time data interaction in satellite common-view time synchronization
CN114637186A (en) * 2020-12-15 2022-06-17 ***通信有限公司研究院 Time difference processing method, equipment and storage medium
CN113325444B (en) * 2021-05-28 2023-02-24 北京卓越航导科技有限责任公司 Beidou satellite system time monitoring method and device based on satellite common view
CN113341686B (en) * 2021-07-02 2022-05-13 长沙学院 Single-satellite multi-calendar timing method and device, computer equipment and storage medium
CN114995098B (en) * 2022-04-22 2024-05-03 中国科学院国家授时中心 Picosecond-level precision satellite interference time transfer method based on phase delay measurement
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

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5577122A (en) * 1994-12-29 1996-11-19 Trimble Navigation Limited Secure communication of information
JP3545640B2 (en) * 1999-04-30 2004-07-21 株式会社東芝 Local area integrated positioning system
ES2350270T3 (en) * 2007-12-18 2011-01-20 Thales Alenia Space Italia S.P.A. SYNCHRONIZATION PROCEDURE FOR NODES OF A CORRESPONDING NETWORK, SYSTEM AND DEVICE.
CN201294513Y (en) * 2008-09-26 2009-08-19 华东电网有限公司 System for synchronizing whole network time based on satellite common vision
JP2011164902A (en) * 2010-02-09 2011-08-25 National Institute Of Advanced Industrial Science & Technology Traceable time standardization method
CN102023290B (en) * 2010-11-04 2013-04-10 中国民用航空总局第二研究所 High-precision distributed pulse signal time difference of arrival detection system

Also Published As

Publication number Publication date
CN104181550A (en) 2014-12-03

Similar Documents

Publication Publication Date Title
CN104181550B (en) A kind of common-view time Frequency Transfer method based on Beidou satellite navigation system
CN105629263B (en) A kind of troposphere atmosphere delay estimation error correcting method and correction system
Strasser et al. Processing of GNSS constellations and ground station networks using the raw observation approach
CN103197340B (en) Gridding real-time monitoring method for total electron content of ionized layer
EP3410144B1 (en) High-precision, real-time satellite positioning device and method thereof
CN105158780B (en) One kind is based on the interchangeable navigation locating method of a variety of aeronautical satellites
CN109154668A (en) Accurate low latency GNSS satellite clock estimation
CN108196272A (en) A kind of satellite navigation positioning device and method based on real-time accurate One-Point Location
CN106850177A (en) A kind of real-time time transmission method based on Static Precise Point Positioning
CN103728643B (en) With the Big Dipper three network RTK blur level single epoch fixing means frequently that wide lane retrains
CN103163533B (en) Seamless fusion expression and correction method of global navigation satellite system (GNSS) global and regional ionospheric delay
CN103760572A (en) Single-frequency PPP ionosphere weighting method based on region CORS
CN108919305A (en) Beidou ground enhances band-like method of servicing and system in communications and transportation
CN103516457A (en) High-precision remote time synchronization method
CN101887128A (en) Method for determining inter-frequency deviation of navigation satellite of global satellite navigation system
CN103576166B (en) Based on the method and system that single satellite positions interference source
CN105717784B (en) Electric power time service method
CN109709584A (en) A kind of GNSS Timing Receiver and time service method
CN107332637A (en) Two-way satellite time and frequency transfer method based on dual carrier phase real-Time Compensation
CN106292265A (en) A kind of many ground based on aeronautical satellite method for synchronizing time
Banville et al. Defining the basis of an integer-levelling procedure for estimating slant total electron content
CN113031036B (en) Ionosphere phase flicker factor construction method based on GNSS 30s sampling frequency data
CN107196716A (en) Calculate the difference method of long wave ground wave signals propagated time delay
CN104536023B (en) Non-delayed sub-meter differential positioning method of high and low frequency error shunting prediction
CN114088080A (en) Positioning device and method based on multi-sensor data fusion

Legal Events

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