CN109100931A - A kind of calculated using differential data user ASF carries out accurate modified method - Google Patents
A kind of calculated using differential data user ASF carries out accurate modified method Download PDFInfo
- Publication number
- CN109100931A CN109100931A CN201810825317.7A CN201810825317A CN109100931A CN 109100931 A CN109100931 A CN 109100931A CN 201810825317 A CN201810825317 A CN 201810825317A CN 109100931 A CN109100931 A CN 109100931A
- Authority
- CN
- China
- Prior art keywords
- asf
- user
- time
- value
- difference
- 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/08—Setting the time according to the time information carried or implied by the radio signal the radio signal being broadcast from a long-wave call sign, e.g. DCF77, JJY40, JJY60, MSF60 or WWVB
- G04R20/10—Tuning or receiving; Circuits therefor
-
- 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/08—Setting the time according to the time information carried or implied by the radio signal the radio signal being broadcast from a long-wave call sign, e.g. DCF77, JJY40, JJY60, MSF60 or WWVB
- G04R20/12—Decoding time data; Circuits therefor
-
- G—PHYSICS
- G04—HOROLOGY
- G04R—RADIO-CONTROLLED TIME-PIECES
- G04R40/00—Correcting the clock frequency
- G04R40/06—Correcting the clock frequency by computing the time value implied by the radio signal
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The present invention provides a kind of calculated using differential data user ASF to carry out accurate modified method, establishes differential reference station, builds long wave differential reference measuring system, acquire difference station ASF measurement data in real time, then calculates difference station ASF a reference value;It plans user's grid, in selected each user's grid points measurement ASF data and calculates ASF a reference value;User to be measured searches its four adjacent known mesh point ASF a reference value, calculates separately difference station ASF undulating value corresponding to the measurement finish time that four points are recorded, and subtract corresponding deviation with four user's ASF a reference values respectively;The accurate ASF value of user position is calculated in the way of interpolation.The present invention improves user's timing accuracy, improves the availability of long wave time service system, avoids the measurement of grid points ASF value and needs to synchronize brought Project Realization problem.
Description
Technical field
The present invention relates to a kind of long wave user time service precision modification methods, in particular to a kind of to utilize difference station ASF data
Modified method is synchronized to user's grid points ASF reference data.
Background technique
Long-wave time service is a kind of wide area radio continental rise time service mode.Before GNSS system appearance, long wave time service system
It is the Service of Timing means of China's full accuracy.Long wave time service system uses the time reference of broadcast control system on ground,
High-power emission system is used simultaneously, and compared to GNSS system, long wave system is good with system stability, signal is anti-interference
The strong feature of ability, this just compensates for the deficiency of GNSS system.Long wave system time service precision after carrying out modernization restructuring can
It is substantially suitable with GNSS time service precision to be promoted from 1 microsecond to 100 nanoseconds, it, can as the standby redundancy means of GNSS time service
With the availability and reliability of significant increase system.
Long wave user's time service precision depends on the measurement accuracy of the ASF value of its position.The user of industry research at present
The main method of ASF measurement is grid calculating method.Inserted by the four known grid points ASF data adjacent with user's point
Value is calculated.This method main problem is that the time of measuring of each grid points user ASF a reference value not can guarantee unanimously,
Each user's grid points need a longer period from addressing to rigging up and debugging, again to DATA REASONING, a grid regions
The time of measuring interval of each point is by a few hours to a couple of days etc. in domain.Certain fluctuation is had since ASF value changes over time,
When therefore will lead to calculating user ASF value, each grid points ASF a reference value has certain with the current time accurate ASF value
Deviation generates the error of dozens to hundreds of nanosecond so as to cause user's ASF value calculated result.The error will lead to user's timing
Precision is unable to reach the index request of 100 nanoseconds.Therefore, to meet long wave system high-precision time service demand, it is necessary to existing
User's ASF calculation method improves, and reduces and calculates error.
Summary of the invention
For overcome the deficiencies in the prior art, the present invention is provided a kind of calculated using differential data user ASF and carries out essence
True modified method, the consistent spy of user's ASF fluctuating change in range is functioned using long wave difference station ASF
Point repairs known grid points user ASF a reference value by calculating difference station ASF data fluctuations value near user's grid points
Just, and using ASF value of the revised grid points ASF value to user's point to be measured interpolation calculation is carried out, to overcome existing method
In each lattice user site ASF reference data it is asynchronous caused by the relatively low problem of user's ASF computational accuracy.
The technical solution adopted by the present invention to solve the technical problems the following steps are included:
Step 1, it is built using long-wave receiver, high-precision GNSS Timing Receiver, time-interval counter and industrial personal computer
Long wave reference measuring system, long-wave receiver and GNSS receiver antenna are placed in observation point known to coordinate;
Step 2, the time difference of time-interval counter real-time measurement long-wave receiver and the 1PPS of GNSS receiver output,
That is the real-time ASF measured value in difference station;Industrial personal computer continuously records 1 hour real-time ASF measured value in difference station, conduct of averaging
Difference station ASF a reference value;
Step 3, it is built using long-wave receiver, high-precision GNSS Timing Receiver, time-interval counter and industrial personal computer
Long wave reference measuring system, long-wave receiver and GNSS receiver antenna are placed in user's grid points;
Step 4, the time difference of time-interval counter real-time measurement long-wave receiver and the 1PPS of GNSS receiver output,
Industrial personal computer continuously records 1 hour time difference, averages as user's grid points ASF a reference value, and record time difference measurement to terminate
Time value;
Step 5, step 3,4 are repeated, the ASF a reference value of four or more user's grid points is obtained, establishes user's grid ASF
Benchmark database;
Step 6, it is adjacent according to the position where oneself to search itself from user's grid ASF pattern library by user U to be measured
The ASF a reference value ASF of four user's grid points C1, C2, C3, C41、ASF2、ASF3、ASF4;
Step 7, the real-time ASF measured value in difference station that query steps 2 obtain finds tetra- points of C1, C2, C3, C4 and is recorded
The real-time ASF measured value in measurement finish time corresponding difference station;
Step 8, the real-time ASF in difference station corresponding to the measurement finish time that C1, C2, C3, C4 point are recorded is calculated separately
The deviation delta 1 of measured value and difference station ASF a reference value, Δ 2, Δ 3, Δ 4;
Step 9, Δ 1, Δ 2, Δ 3, Δ 4 are subtracted with C1, C2, C3, C4 user's grid points ASF a reference value respectively, obtained
The revised a reference value ASF of C1, C2, C3, C4c1、ASFc2、ASFc3、ASFc4;
Step 10, the accurate ASF value ASF of the position user U to be measured is calculated in the way of interpolationU,
ASFU=(1- α) [(1- β) * ASFC1+β*ASFC3]+α[(1-β)*ASFC2+β*ASFC4] wherein user U to be measured
Coordinate is (x, y), and point (x (i), y (i)) indicates the mesh coordinate in the subscriber station lower left corner,
α=(x-x (i))/(x (i+1)-x (i)), p=(y-y (j))/(y (j+1)-y (j)).
The beneficial effects of the present invention are: improving user's timing accuracy using the measurement data at long wave difference station, improve
The availability of long wave time service system.Meanwhile the invention avoids the measurements of grid points ASF value to need to synchronize brought engineering reality
Existing problem.
Detailed description of the invention
Fig. 1 is ASF data measurin system equipment connecting relation figure of the invention;
Fig. 2 is user's grid coordinate relational graph.
Specific embodiment
Present invention will be further explained below with reference to the attached drawings and examples, and the present invention includes but are not limited to following implementations
Example.
The present invention includes the following steps that (for a user, step 1 to step 5 does not need to execute, it is only necessary to hold
Row step 6 is to step 10):
Step 1, differential reference station is established.Get out long-wave receiver, high-precision GNSS Timing Receiver, time interval
Counter, industrial personal computer, build long wave reference measuring system as shown in Figure 1, and long-wave receiver and GNSS receiver antenna are placed in
In observation point known to coordinate.
Step 2, counter real-time measurement long-wave receiver and GNSS receiver output 1PPS the time difference, this when difference
According to being exactly the real-time ASF measured value in difference station, and send result to industrial personal computer.Industrial personal computer continuously records ASF measured value, and counts
The average value of the 1st hour time difference data is calculated as difference station ASF a reference value.
Step 3, user's grid points ASF data are established.With step 1, long wave reference measurement system is established in selected user's point
System.
Step 4, the time difference of user's station counting unit real-time measurement long-wave receiver and the 1PPS of GNSS receiver output, and
Industrial personal computer is sent result to, is measured 1 hour altogether.Industrial personal computer calculates this 1 hour time difference data average value as user's grid
Point ASF a reference value, and record time difference data measurement end time value.
Step 5, step 3,4 are repeated, measurement obtains four or more user's grid points ASF data, and by grid points ASF number
According to deposit user's grid ASF benchmark database.
Step 6, user U to be measured searched from user's grid ASF pattern library according to the position where oneself its adjacent four
A known mesh point (being denoted as C1, C2, C3, C4) ASF a reference value, is denoted as ASF1、ASF2、ASF3、ASF4。
Step 7, in query steps 2 difference station real-time measurement ASF data, find tetra- points of C1, C2, C3, C4 and recorded
Measurement finish time corresponding difference ASF value.
Step 8, on the basis of four grid points ASF values corresponding record moment, C1, C2, C3, C4 point are calculated separately
The deviation of the difference station time difference ASF value and difference ASF a reference value corresponding to the measurement finish time recorded, is denoted as Δ 1, Δ
2、Δ3、Δ4。
Step 9, Δ 1, Δ 2, Δ 3, Δ 4 are subtracted with C1, C2, C3, C4 grid points user's ASF a reference value respectively, obtained
The revised a reference value ASF of C1, C2, C3, C4c1、ASFc2、ASFc3、ASFc4。
ASFCi=ASFi-Δi
Step 10, the accurate ASF value that the position user U is calculated in the way of interpolation, is expressed as ASFU.Such as Fig. 2 institute
User's mesh point C1, C2, C3, C4 coordinate and corresponding A SF a reference value shown, the coordinate of user U can be expressed as (x, y), (x
(i), y (i)) put the mesh coordinate for indicating the subscriber station lower left corner.Therefore, according to the position user U in grid points in Fig. 2
Coordinate relationship, the accurate ASF value ASFU of user location can be calculated.
α=(x-x (i))/(x (i+1)-x (i))
β=(y-y (j))/(y (j+1)-y (j))
According to the relationship of subscriber station and grid points, the ASF value of user's point to be measured, formula can be calculated are as follows:
ASFU=(1- α) [(1- β) * ASFC1+β*ASFC3]+α[(1-β)*ASFC2+β*ASFC4]
From above-mentioned implementation steps it is found that the present invention is described using difference station ASF data to user's grid points ASF benchmark
The method that data are modified, implementation process mainly include the real-time measurement of difference station ASF value, grid points user's ASF benchmark
Value is established, and the amendment of grid points ASF a reference value, user's ASF a reference value calculate four parts.Difference station ASF measurement and grid points
User's ASF a reference value measurement belongs to system-level operations content.User class only needs to carry out the amendment of grid points ASF a reference value and uses
Family ASF a reference value calculates work.
As can be seen from the above embodiments, the ASF for being mainly characterized by being utilized synchronization long wave difference station of the invention
The fluctuating change feature consistent with user's ASF fluctuating change, by one reduction of ASF primary system of each user's grid points to tested point
ASF calculate the moment, thus eliminate each grid points ASF variation it is inconsistent caused by user ASF the problem of calculating error.This
It invents proposed method and needs to configure long wave difference station, lasting, stable difference monitoring data are provided.User's ASF lattice are needed simultaneously
Grid database service.Differential system and ASF grid data bank are built as China's big science engineering high-precision ground time dissemination system
A part, have adequately support and ensure.For user point ASF calculating, calculating process is simple, as long as inquiry difference number
It can be obtained by accurate calculated result according to ASF grid data bank.Therefore, the present invention is for users, simple and practical, no
Additional resource overhead is needed, popularization and use are facilitated, there is preferable development prospect and application prospect.
Claims (1)
1. a kind of calculated using differential data user ASF carries out accurate modified method, it is characterised in that include the following steps:
Step 1, long wave is built using long-wave receiver, high-precision GNSS Timing Receiver, time-interval counter and industrial personal computer
Reference measuring system, long-wave receiver and GNSS receiver antenna are placed in observation point known to coordinate;
Step 2, the time difference of time-interval counter real-time measurement long-wave receiver and the 1PPS of GNSS receiver output, i.e. difference
It stands real-time ASF measured value;Industrial personal computer continuously records 1 hour real-time ASF measured value in difference station, averages as difference station
ASF a reference value;
Step 3, long wave is built using long-wave receiver, high-precision GNSS Timing Receiver, time-interval counter and industrial personal computer
Reference measuring system, long-wave receiver and GNSS receiver antenna are placed in user's grid points;
Step 4, the time difference of time-interval counter real-time measurement long-wave receiver and the 1PPS of GNSS receiver output, industrial personal computer
It continuous 1 hour time difference of record, averages as user's grid points ASF a reference value, and record time difference measurement end time value;
Step 5, step 3,4 are repeated, the ASF a reference value of four or more user's grid points is obtained, establishes user's grid ASF base value
According to library;
Step 6, user U to be measured searched from user's grid ASF pattern library according to the position where oneself itself adjacent four
The ASF a reference value ASF of user's grid points C1, C2, C3, C41、ASF2、ASF3、ASF4;
Step 7, the real-time ASF measured value in difference station that query steps 2 obtain finds the survey that tetra- points of C1, C2, C3, C4 are recorded
Measure the finish time corresponding real-time ASF measured value in difference station;
Step 8, the real-time ASF measurement in difference station corresponding to the measurement finish time that C1, C2, C3, C4 point are recorded is calculated separately
The deviation delta 1 of value and difference station ASF a reference value, Δ 2, Δ 3, Δ 4;
Step 9, Δ 1, Δ 2, Δ 3, Δ 4 are subtracted with C1, C2, C3, C4 user's grid points ASF a reference value respectively, obtain C1, C2,
The revised a reference value ASF of C3, C4c1、ASFc2、ASFc3、ASFc4;
Step 10, the accurate ASF value ASF of the position user U to be measured is calculated in the way of interpolationU,
ASFU=(1- α) [(1- β) * ASFC1+β*ASFC3]+α[(1-β)*ASFC2+β*ASFC4]
Wherein the coordinate of user U to be measured is (x, y), and point (x (i), y (i)) indicates the mesh coordinate in the subscriber station lower left corner, α=(x-
X (i))/(x (i+1)-x (i)), β=(y-y (j))/(y (j+1)-y (j)).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810825317.7A CN109100931B (en) | 2018-07-25 | 2018-07-25 | Method for accurately correcting user ASF calculation by using differential data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810825317.7A CN109100931B (en) | 2018-07-25 | 2018-07-25 | Method for accurately correcting user ASF calculation by using differential data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109100931A true CN109100931A (en) | 2018-12-28 |
CN109100931B CN109100931B (en) | 2020-07-14 |
Family
ID=64847427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810825317.7A Active CN109100931B (en) | 2018-07-25 | 2018-07-25 | Method for accurately correcting user ASF calculation by using differential data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100931B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112636893A (en) * | 2020-11-29 | 2021-04-09 | 中国科学院国家授时中心 | Method for improving eLoran system time service precision by using ASF grid and differential station |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101059543A (en) * | 2007-05-16 | 2007-10-24 | 西安理工大学 | Measuring method for time-delay correction value for long wave transmission |
CN102539939A (en) * | 2012-02-14 | 2012-07-04 | 西安理工大学 | High-precision marine ASF (Additional Secondary Factor) correcting method based on ground equivalent conductivity inversion |
KR20160000725A (en) * | 2014-06-25 | 2016-01-05 | 연세대학교 산학협력단 | Eloran asf survey system and method thereof |
CN107196716A (en) * | 2017-04-21 | 2017-09-22 | 中国科学院国家授时中心 | Calculate the difference method of long wave ground wave signals propagated time delay |
-
2018
- 2018-07-25 CN CN201810825317.7A patent/CN109100931B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101059543A (en) * | 2007-05-16 | 2007-10-24 | 西安理工大学 | Measuring method for time-delay correction value for long wave transmission |
CN102539939A (en) * | 2012-02-14 | 2012-07-04 | 西安理工大学 | High-precision marine ASF (Additional Secondary Factor) correcting method based on ground equivalent conductivity inversion |
KR20160000725A (en) * | 2014-06-25 | 2016-01-05 | 연세대학교 산학협력단 | Eloran asf survey system and method thereof |
CN107196716A (en) * | 2017-04-21 | 2017-09-22 | 中国科学院国家授时中心 | Calculate the difference method of long wave ground wave signals propagated time delay |
Non-Patent Citations (1)
Title |
---|
李瑞敏等: "长波授时附加二次相位因子修正方法", 《宇航计测技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112636893A (en) * | 2020-11-29 | 2021-04-09 | 中国科学院国家授时中心 | Method for improving eLoran system time service precision by using ASF grid and differential station |
CN112636893B (en) * | 2020-11-29 | 2022-08-02 | 中国科学院国家授时中心 | Method for improving eLoran system time service precision by using ASF grid and differential station |
Also Published As
Publication number | Publication date |
---|---|
CN109100931B (en) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102707300B (en) | Method, device and system for optimizing GPS track | |
CN103454905B (en) | A kind of method for synchronizing time | |
CN107976702A (en) | A kind of position correcting method based on CORS, positioning terminal and alignment system | |
CN110856106A (en) | Indoor high-precision three-dimensional positioning method based on UWB and barometer | |
CN114019584A (en) | VRS resolving method for high-precision CORS network in large-altitude-difference area | |
CN105182069A (en) | High resolution group quantization phase processing method under pilot frequency architecture | |
CN109100931A (en) | A kind of calculated using differential data user ASF carries out accurate modified method | |
CN103177392B (en) | A kind of correlating method of coil of strip length process | |
CN102928660A (en) | Reactive power measurement method based on finite impulse response (FIR) digital filter | |
CN114125699A (en) | Network RTK service method for reconstruction by using virtual reference station | |
CN106125044B (en) | Offline localization method based on gradient decline | |
CN107632137B (en) | Method for correcting soil radon gas instantaneous measurement data in uranium mine exploration | |
CN104614981B (en) | A kind of method of the sky frequency stability obtaining atomic clock remote calibration system | |
CN201837717U (en) | Vertical ground feature positioning error correcting system based on InSAR (interferometric synthetic aperture radar) images | |
CN115189783A (en) | Performance evaluation method of satellite navigation enhancement system reference station | |
CN109472295B (en) | Pearson correlation coefficient-based WAMS/SCADA system time synchronization sorting method | |
CN110381580B (en) | WiFi positioning method based on proportion optimization | |
CN111856530A (en) | Combined optimization method for continuously operating reference station system | |
CN108154279B (en) | Photovoltaic power station system performance ratio weighting on-line evaluation method and system | |
CN110006332A (en) | A kind of triangle deformation monitoring method of seismic precursor deformation GNSS monitoring net | |
CN116045900B (en) | Method for constructing precision retrieval table by replacing leveling measurement with GNSS (Global navigation satellite System) and geodetic level model | |
Liu et al. | Research on Reliable Long-Baseline NRTK Positioning Method Considering Ionospheric Residual Interpolation Uncertainty. Remote Sens. 2023, 15, 5353 | |
CN110794195B (en) | Three-phase synchronous sampling method and device and storage medium | |
CN106250414B (en) | Data storage system and method based on measuring point object | |
Zhiyu et al. | Indoor real-time three-side positioning technology based on low power consumption Bluetooth |
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 |