CN104199054A - Preprocessing method for common view data of Beidou satellite navigation system - Google Patents
Preprocessing method for common view data of Beidou satellite navigation system Download PDFInfo
- Publication number
- CN104199054A CN104199054A CN201410419574.2A CN201410419574A CN104199054A CN 104199054 A CN104199054 A CN 104199054A CN 201410419574 A CN201410419574 A CN 201410419574A CN 104199054 A CN104199054 A CN 104199054A
- Authority
- CN
- China
- Prior art keywords
- data
- delay
- altogether
- view data
- common view
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/23—Testing, monitoring, correcting or calibrating of receiver elements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a preprocessing method for common view data of a Beidou satellite navigation system. The method includes the steps of selecting proper common view data according to an observation state of a time frequency laboratory and a satellite state; correcting the delay during signal transmission of the selected common view data; performing gross error elimination on the corrected common view data; smoothing the common view data through a k-order Vondrak filter, selecting an optimal smoothing factor through iterative and least squares estimation, and eliminating random errors. The preprocessing method for the common view data of the Beidou satellite navigation system is low in calculation amount, high in computation speed and capable of quickly obtaining the proper smoothing factor and effectively eliminating the random errors in the common view data.
Description
Technical field
The present invention relates to data processing method, particularly relate to a kind of data preprocessing method for Beidou satellite navigation system common-view time Frequency Transfer.
Background technology
In Beidou satellite navigation system common-view time Frequency Transfer process, two time-frequency laboratories that participate in looking are altogether at same big-dipper satellite of synchronization observation, obtain respectively the time difference data of local clock and the spaceborne clock of this big-dipper satellite separately, through communication link, the data that record are separately exchanged mutually, through data processing, remove spaceborne clock and can obtain the time difference between two time-frequency laboratories, thereby realize the transmission of Big Dipper common-view time.
Owing to existing and disturbing when receiving Big Dipper signal, for obtaining more accurate local clock and the spaceborne clock time difference, before carrying out exchanges data, need to carry out in this locality data pre-service, comprise: altogether depending on the choosing of data, the correction of signal transmission delay, the eliminating of stochastic error.Signal transmission delay correction refers to be revised geometric delays, spatial medium time-delay, Relativistic Time Delay, receiver hardware delay etc. in time difference data, can realize by setting up the mode of model and calibration, and the weakening of stochastic error can only realize by the mode of data smoothing.
Data smoothing is very important work of Data processing, and user will obtain effective information from measurement data, must do pre-service to measurement data, reduces as far as possible the impact of stochastic error, the quality of simultaneously all right assessment of metrology data.Conventional smoothing method has curve, running mean method and Vondrak smoothing method at present.The prerequisite of curve-fitting method is to provide functional form, and great majority are polynomial expression, and its exponent number is often rule of thumb determined, and human factor impact is larger; Running mean method require data must be equidistantly and also sampled point than comparatively dense, there is significant limitation.Vondrak smoothing method essence is by selecting smoothing factor, to the absolute matching of measurement sequence with definitely smoothly, seeking a balance, the in the situation that of unknown fitting function form, Vondrak method also can be carried out reasonably smoothly survey data, and its key is to select rational smoothing factor.Conventional selection smoothing factor method has at present: observational error method, frequency response method, cross-certification method.Observational error method requires to know in advance observation data priori standard deviation, and we can not know standard deviation information accurately conventionally; Frequency response method can obtain good filter effect for the observation data of knowing frequency or cycle; Cross-certification method is to select smoothing factor by test repeatedly, and operand is large, and speed is slow, and the smoothing factor randomness of selection is also larger.
Therefore, need to provide a kind of preprocess method of looking altogether data, with fast selecting smoothing factor, reduce operand, improve travelling speed.
Summary of the invention
The technical problem to be solved in the present invention is to provide a kind of data preprocessing method for Big Dipper common-view time Frequency Transfer, solve and to look altogether at present Data processing, in the undesirable and traditional Vondrak smoothing method of conventional smoothing method smooth effect, smoothing factor is difficult to the problem of selecting.
For solving the problems of the technologies described above, the present invention adopts following technical proposals.
For Beidou satellite navigation system, look altogether a preprocess method for data, the method step comprises
S1, according to the observer state in time-frequency laboratory and satellitosis, choose the suitable data of looking altogether;
S2, to revising depending on the delay in data signal propagation process altogether after choosing;
S3, to revised, depending on data, carry out gross error altogether and reject to process;
S4, employing k rank Vondrak wave filter, to carrying out smoothing processing depending on data altogether, are chosen the postfitted orbit factor by least squares estimate iteration, eliminate stochastic error.
Preferably, in described step S1, the condition of choosing depending on data is altogether comprised
Guarantee that two time-frequency laboratories are at same satellite of synchronization observation;
The elevation angle that is observed satellite is not less than 15 °;
According to Beidou satellite navigation system feature, extend total length tracking time.
Preferably, need the delay item of revising to comprise geometric distance delay, ionosphere delay, tropospheric delay in described step S2, sagnac effect and receiver postpone.
Preferably, Deferred Correction method comprises
Utilize dual-frequency data to revise ionosphere delay;
Delay parameter and the satellite ephemeris modeling of according to Beidou satellite navigation system receiver, in row text, extracting, eliminate geometric distance delay, tropospheric delay and sagnac effect;
Mode by calibrating receiver overcomes receiver and postpones.
Preferably, in described step S3, adopt 3 σ criterions to reject and look altogether the gross error in data.
Preferably, adopt least square process of iteration to choose smoothing factor in described step S4, its step comprises
S41, the time difference sequence of making even after sliding are estimated as waiting, by Vondrak filtering ultimate principle, using time difference sequence and time difference sequence k jump divide as two class observed readings, set up respectively error equation, thereby Vondrak filtering be converted into least square adjustment problem;
S42, given initial smoothing factor, utilize least square method to solve the error equation of observed reading, obtains once level and smooth result;
S43, the method for utilizing Helmert component of variance to estimate, carry out estimation variance component according to the correction vector V last time smoothly obtaining, and revises initial smoothing factor;
S44, by iteration repeatedly, obtain the postfitted orbit factor, Simultaneous Iteration finishes to complete elimination stochastic error.
Beneficial effect of the present invention is as follows:
Technical scheme calculated amount of the present invention is little, and fast operation can obtain suitable smoothing factor fast, can effectively eliminate the stochastic error of looking altogether in data.
Accompanying drawing explanation
Below in conjunction with accompanying drawing, the specific embodiment of the present invention is described in further detail;
Fig. 1 illustrates the schematic diagram of a kind of data preprocessing method for Big Dipper common-view time Frequency Transfer of the present invention.
Embodiment
Below in conjunction with one group of embodiment and accompanying drawing, the present invention is described further.
The invention discloses a kind of data preprocessing method transmitting for Big Dipper common-view time, the method concrete steps comprise:
The first step, choosing depending on data altogether
In Big Dipper common-view time Transfer Technology, choose while looking data altogether, except what receive, be same satellite-signal, also must meet three conditions below simultaneously:
1) tracking is constantly identical, and same satellite must be observed at synchronization in two time-frequency laboratories, so just can eliminate the error of satellite clock completely;
2) elevation angle is not less than 15 °, and when elevation of satellite is too low, ionosphere delay error is difficult to correct, and Multi-Path Effects aggravation;
3) follow the tracks of length, for GPS, one time total length is tracked as 13 minutes, add preset preparation in 2 minutes and the processing of 1 minute data, for the Big Dipper, look altogether, because Beidou satellite navigation system is containing GEO and IGSO satellite, can total length tracking time of proper extension.
Second step, signal transmission delay correction
Choose and look altogether after data, the every delay of reply satellite-signal in communication process, comprises geometric distance delay, ionosphere delay, tropospheric delay, and sagnac effect and receiver delay etc. is revised, and is deducted in altogether depending on data.Wherein, ionosphere delay can utilize dual-frequency data correction, geometric distance delay, tropospheric delay, sagnac effect etc. can be eliminated according to modelings such as the delay parameter extracting from Beidou receiver navigation message, satellite ephemeris, and receiver postpones by the mode of calibrating, to obtain.
The 3rd step gross error is rejected
Because Vondrak is smoothly subject to the impact of gross error, before carrying out data smoothing, first adopt 3 σ criterions to reject and look altogether the gross error in data.Owing to being once per second depending on data altogether, twice measurement in front and back is more or less the same, therefore data mean value substitutes the rough error weeding out before and after can directly utilizing rough error.
The 4th step Vondrak is level and smooth
Adopt method that Vondrak is level and smooth to carry out data smoothing depending on data altogether after signal propagation delays correction and elimination of rough difference, eliminate stochastic error.
The ultimate principle of the level and smooth measurement data of Vondrak method is:
Q=F+λ
2S=min (1)
Wherein, F=Σ p
i(y
i'-y
i)
2, S=Σ (Δ
ky
i)
2, y
i(i=1,2 ..., N) for time argument is x
isurvey data.Y
i' be smooth value to be asked, p
iweight for measurement data.
This programme adopts k rank Vondrak wave filter as follows to looking altogether smoothed data modeling:
with
on interval between 2, with one, pass through i, i+1 ..., the i+k altogether k rank lagrange polynomial of k+1 point defines level and smooth result y '=y ' (t), and corresponding, the k jump of smooth value divides and can be expressed as
Wherein,
Make y
ifor observed reading, y
i' be to be valuated, by formula (1), set up two class error equations:
Make first kind observed reading power for P
1, the second kind observation value power is P
2.
According to least square, have:
That is:
Make P
1=ε p, P
2=I:
Order
above formula is consistent with Vondrak filtering ultimate principle, and Vondrak is smoothly converted into the described least square adjustment problem in formula (3)~(5).
Get quantity of state X=[y
1' y
2' ... y
n']
t, formula (3) can be expressed as:
Wherein,
f
1=[y
1y
2y
n]
t, f
2=[0 0 ... 0]
t
According to the principle of least square, solve, obtain the estimated value of quantity of state
Wherein,
The 5th step component of variance is estimated, right value update
From the 4th step, Vondark filtering is converted into the least square adjustment problem with the different observed readings of two classes, its key is the appropriate power of choosing all kinds of observed readings, for this reason, the method of utilizing Helmert component of variance to estimate, according to step 3 gained adjustment result, calculate correction vector V and carry out estimation variance component, revise adjustment weights used.Step is as follows:
For given initial smoothing factor ε=1/ λ
2, test accordingly front weights and be respectively P
1=ε p, P
2=I, carries out pre-adjustment according to formula (8), obtains
p
1v
1,
p
2v
2, can be calculated component of variance estimated value
and
Rear weights are tested in calculating
Wherein, C is constant.
The 6th step iteration is selected the postfitted orbit factor
Get
smoothing factor
can find out, each iteration is all to P
1, smoothing factor ε is revised.Using revised weights as testing the new adjustment of front weight, repeatedly carry out, adjustment-> component of variance estimation-> right value update is adjustment again, until twice of front and back iteration smoothing factor is consistent.
The 7th step iteration finishes, and obtains the pretreated data of looking altogether.
While utilizing the present invention to adopt k rank Vondrak wave filter to carry out data smoothing, increase along with filter order, its frequency response curve steepening, the effect of cross frequence component is better, but calculated amount is larger, and require sampling density larger, at the Big Dipper, look altogether Data processing, between local clock and the spaceborne clock of the Big Dipper, time difference data is once per second, it is 20min that one secondary tracking duration is set, each tracking period need to be carried out data smoothing to 1200 data, and calculated amount is larger, therefore select three rank Vondrak wave filters to carry out data smoothing.In addition, consider at one and follow the tracks of in the period, because receiver acquisition satellite-signal needs certain hour, before and after each comparison period, the data error of 2min is larger, for reducing rough error for the impact of level and smooth result, each can be compared to the data of period front and back 2min and remove, then do smoothing processing.
In sum, technical scheme of the present invention.
Obviously; the above embodiment of the present invention is only for example of the present invention is clearly described; and be not the restriction to embodiments of the present invention; for those of ordinary skill in the field; can also make other changes in different forms on the basis of the above description; here cannot give all embodiments exhaustive, every still row in protection scope of the present invention of apparent variation that technical scheme of the present invention extends out or change that belong to.
Claims (6)
1. for Beidou satellite navigation system, look altogether a preprocess method for data, it is characterized in that, the method step comprises
S1, according to the observer state in time-frequency laboratory and satellitosis, choose the suitable data of looking altogether;
S2, to revising depending on the delay in data signal propagation process altogether after choosing;
S3, to revised, depending on data, carry out gross error altogether and reject to process;
S4, employing k rank Vondrak wave filter, to carrying out smoothing processing depending on data altogether, are chosen the postfitted orbit factor by least squares estimate iteration, eliminate stochastic error.
2. data preprocessing method according to claim 1, is characterized in that, in described step S1, the condition of choosing depending on data is altogether comprised
Guarantee that two time-frequency laboratories are at same satellite of synchronization observation;
The elevation angle that is observed satellite is not less than 15 °;
According to Beidou satellite navigation system feature, extend total length tracking time.
3. data preprocessing method according to claim 1, is characterized in that, needs the delay item of revising to comprise geometric distance delay, ionosphere delay, tropospheric delay in described step S2, and sagnac effect and receiver postpone.
4. data preprocessing method according to claim 3, is characterized in that, Deferred Correction method comprises
Utilize dual-frequency data to revise ionosphere delay;
Delay parameter and the satellite ephemeris modeling of according to Beidou satellite navigation system receiver, in row text, extracting, eliminate geometric distance delay, tropospheric delay and sagnac effect;
Mode by calibrating receiver overcomes receiver and postpones.
5. data preprocessing method according to claim 1, is characterized in that, adopts 3 σ criterions to reject and look altogether the gross error in data in described step S3.
6. data preprocessing method according to claim 1, is characterized in that, adopts least square process of iteration to choose smoothing factor in described step S4, and its step comprises
S41, the time difference sequence of making even after sliding are estimated as waiting, by Vondrak filtering ultimate principle, using time difference sequence and time difference sequence k jump divide as two class observed readings, set up respectively error equation, thereby Vondrak filtering be converted into least square adjustment problem;
S42, given initial smoothing factor, utilize least square method to solve the error equation of observed reading, obtains once level and smooth result;
S43, the method for utilizing Helmert component of variance to estimate, carry out estimation variance component according to the correction vector V last time smoothly obtaining, and revises initial smoothing factor;
S44, by iteration repeatedly, obtain the postfitted orbit factor, Simultaneous Iteration finishes to complete elimination stochastic error.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410419574.2A CN104199054A (en) | 2014-08-22 | 2014-08-22 | Preprocessing method for common view data of Beidou satellite navigation system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410419574.2A CN104199054A (en) | 2014-08-22 | 2014-08-22 | Preprocessing method for common view data of Beidou satellite navigation system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104199054A true CN104199054A (en) | 2014-12-10 |
Family
ID=52084366
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410419574.2A Pending CN104199054A (en) | 2014-08-22 | 2014-08-22 | Preprocessing method for common view data of Beidou satellite navigation system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104199054A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104460311A (en) * | 2014-12-30 | 2015-03-25 | 四川九洲电器集团有限责任公司 | Time calibration method and device |
CN107607971A (en) * | 2017-09-08 | 2018-01-19 | 哈尔滨工程大学 | Temporal frequency transmission method and receiver based on GNSS common-view time alignment algorithms |
CN111030774A (en) * | 2019-12-19 | 2020-04-17 | 北京无线电计量测试研究所 | Real-time common-view data processing method based on Beidou satellite navigation system |
CN111368588A (en) * | 2018-12-25 | 2020-07-03 | 天津大学 | Tidal observation data preprocessing method based on Vondrak filtering |
CN113589679A (en) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | Satellite precision time transfer method based on interferometry |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2023456A1 (en) * | 2007-08-06 | 2009-02-11 | Abb Research Ltd. | Estimating a time offset between stationary clocks |
CN102624513A (en) * | 2012-03-06 | 2012-08-01 | 北京无线电计量测试研究所 | Device for verifying synchronization precision of two-way satellite time transfer modem |
-
2014
- 2014-08-22 CN CN201410419574.2A patent/CN104199054A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2023456A1 (en) * | 2007-08-06 | 2009-02-11 | Abb Research Ltd. | Estimating a time offset between stationary clocks |
CN102624513A (en) * | 2012-03-06 | 2012-08-01 | 北京无线电计量测试研究所 | Device for verifying synchronization precision of two-way satellite time transfer modem |
Non-Patent Citations (4)
Title |
---|
佘娣: "卫星高度角变化对大气折射的影响", 《全球定位系统》 * |
周晓卫: "HVF方法在GPS多路径效应研究中的应用", 《大地测量与地球动力学》 * |
杨帆: "基于北斗GEO和IGSO卫星的高精度共视时间传递", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
许国宏: "北斗共视授时技术研究与设计", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104460311A (en) * | 2014-12-30 | 2015-03-25 | 四川九洲电器集团有限责任公司 | Time calibration method and device |
CN107607971A (en) * | 2017-09-08 | 2018-01-19 | 哈尔滨工程大学 | Temporal frequency transmission method and receiver based on GNSS common-view time alignment algorithms |
CN107607971B (en) * | 2017-09-08 | 2021-01-12 | 哈尔滨工程大学 | Time frequency transmission method based on GNSS common-view time comparison algorithm and receiver |
CN111368588A (en) * | 2018-12-25 | 2020-07-03 | 天津大学 | Tidal observation data preprocessing method based on Vondrak filtering |
CN111030774A (en) * | 2019-12-19 | 2020-04-17 | 北京无线电计量测试研究所 | Real-time common-view data processing method based on Beidou satellite navigation system |
CN113589679A (en) * | 2021-06-17 | 2021-11-02 | 中国科学院国家授时中心 | Satellite precision time transfer method based on interferometry |
CN113589679B (en) * | 2021-06-17 | 2022-06-10 | 中国科学院国家授时中心 | Satellite precision time transfer method based on interferometry |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291645B (en) | The volume kalman filter method coupled deeply suitable for higher-dimension GNSS/INS | |
CN105629263B (en) | A kind of troposphere atmosphere delay estimation error correcting method and correction system | |
Psychas et al. | Assessment of ionospheric corrections for PPP-RTK using regional ionosphere modelling | |
CN108828640B (en) | Method and device for weighting satellite navigation positioning observation values | |
CN109459778B (en) | Code pseudo range/Doppler joint velocity measurement method based on robust variance component estimation and application thereof | |
CN102288978B (en) | Continuous operational reference system (CORS) base station cycle slip detection and recovering method | |
CN107678050B (en) | GLONASS phase inter-frequency deviation real-time tracking and precise estimation method based on particle filtering | |
CN104199054A (en) | Preprocessing method for common view data of Beidou satellite navigation system | |
CN104020480B (en) | A kind of satellite navigation method of the interactive multi-model UKF with adaptive factor | |
CN105842719B (en) | A kind of CORS base station net baseline Ambiguity Solution Methods for taking troposphere influence into account | |
CN108107455A (en) | A kind of satellite clock correction Real-time Forecasting Method based on phase hit | |
CN104614741B (en) | Real-time precise satellite clock error estimation method not impacted by deviation of code frequency of GLONASS | |
CN105301617B (en) | A kind of integer ambiguity validity check method in satellite navigation system | |
CN110208836B (en) | GNSS high-adaptability cycle slip detection and restoration method based on Kalman filtering | |
CN108196267B (en) | GNSS CP technology-based uninterrupted time transfer method | |
CN105807300A (en) | Method for high-precision dynamic point positioning through big dipper double frequency receiver | |
Delporte et al. | GPS Carrier-Phase Time Transfer Using Single-Difference Integer Ambiguity Resolution. | |
CN103698785A (en) | Online-sequential extreme learning machine-based satellite signal cycle slip detection and restoration method | |
CN103197326A (en) | Multi-constellation single base station receiver clock difference estimation method | |
CN104199061B (en) | A kind of method for establishing GPS system and BDS system carrier phase frequency standards | |
Shen et al. | Determination of the geopotential and orthometric height based on frequency shift equation | |
CN109521444B (en) | Self-adaptive least square fitting estimation algorithm for GPS horizontal velocity field of crustal movement | |
CN104316961A (en) | Method for obtaining geological parameters of weathered layer | |
Yin et al. | A novel cycle slips detection model for the high precision positioning | |
CN103760582A (en) | Method for optimizing satellite double-difference observation structure in occlusion environment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20141210 |
|
RJ01 | Rejection of invention patent application after publication |