CN103384145A - Method for utilizing time window to eliminate pulse per second jitter - Google Patents
Method for utilizing time window to eliminate pulse per second jitter Download PDFInfo
- Publication number
- CN103384145A CN103384145A CN2013102694489A CN201310269448A CN103384145A CN 103384145 A CN103384145 A CN 103384145A CN 2013102694489 A CN2013102694489 A CN 2013102694489A CN 201310269448 A CN201310269448 A CN 201310269448A CN 103384145 A CN103384145 A CN 103384145A
- Authority
- CN
- China
- Prior art keywords
- pulse per
- pps
- crystal oscillator
- local crystal
- time window
- 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
Images
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Measurement Of Unknown Time Intervals (AREA)
- Manipulation Of Pulses (AREA)
Abstract
The invention discloses a method for utilizing a time window to eliminate pulse per second jitter, and belongs to the field of high-precision measurement. The method comprises the following steps: step 1, the length of the time window is set as n, and a local crystal oscillator starts counting when an initial pulse per second arrives; step 2, i is made equal to 1; step 3, when an i-th pulse per second after the initial pulse per second arrives, the local crystal oscillator stops counting, and a counting value Ni is stored; step 4, the local crystal oscillator is made to count again, and the i is made equal to i plus 1; step 5, the step 3 is returned to if the i is smaller than n, and a step 6 is conducted if the i is not smaller than the n; step 6, when the i-th pulse per second after the initial pulse per second arrives, the local crystal oscillator stops counting, another counting value Ni is stored, and a calculated mean value and a mean value are averaged to obtain an Ai which is stored, wherein the mean difference of the mean value is in a set range; step 7, the local crystal oscillator is made to count again, and the i is made equal to i plus 1; step 8, the step 6 is returned to if pulse per second transmission is not finished, and a step 9 is conducted if the pulse per second transmission is finished; step 9, the local crystal oscillator is used for simulating a pulse per second according to all the stored values. The method for utilizing the time window to eliminate the pulse per second jitter is suitable for eliminating the pulse per second jitter.
Description
Technical field
The invention belongs to the high-acruracy survey field, be specifically related to a kind of method of eliminating the pulse per second (PPS) shake.
Background technology
At present, all need high-precision time measurement in the fields such as satellite navigation, target acquisition and location, traditional time reference take second as unit is practical requirement not, especially aspect the aerospace military project and in large-scale Physical Experiment, time interval measurement is as a kind of important discriminating and detection means, very strict to the required precision of measuring, sometimes need to reach the ps magnitude, thereby the precision that how to improve time measurement becomes the important topic in this field.
Split-second precision is measured and mainly focused on two aspects at present: long-time certainty of measurement and short time measurement precision wherein generally are defined as the short time to the time interval less than 1s.
Wherein long-time the measurement mainly used the GPS pulse per second (PPS), and along with the development of GPS technology, GPS pulse per second (PPS) precision can reach μ s magnitude at present.
The crystal oscillator clock that short time measurement adopts and GPS pulse per second (PPS) combination, because the GPS pulse per second (PPS) is when carrying out short time measurement, its jitter error will affect certainty of measurement greatly, therefore need constantly to eliminate the jitter error of pulse per second (PPS) in measuring process, crystal oscillator clock is revised repeatedly the raising that finally realizes accuracy, therefore the effect of pulse per second (PPS) jitter elimination has directly determined the precision of time measurement, and the disappear research of the technology of trembling of pulse per second (PPS) is the core content in this field.
In prior art, the pulse per second (PPS) jitter elimination mainly adopts following methods: because the cumulative errors of GPS pulse per second (PPS) are less, and the random error of crystal oscillator clock (ns level) is less, so traditional technology of trembling that disappears makes both reference each other according to the characteristics of pulse per second (PPS) and the complementation of crystal oscillator clock precision to both comparing to analyze.At first usage counter and comparator carry out frequency division to crystal oscillator clock, produce the crystal oscillator pps pulse per second signal, phase place to crystal oscillator pulse per second (PPS) and GPS pps pulse per second signal compares, and produces biased sequence, and this deviation comprises the cumulative errors of jitter error and the crystal oscillator of GPS pulse per second (PPS).Adopt mathematical regression model to two kinds of errors are estimated, thereby isolate error separately, and the crystal oscillator cumulative errors are carried out the online Active Compensation of least square method, realize the elimination of GPS pulse jitter error and the raising of clock accuracy.
Although can bring up to 100ns through the time precision that the above-mentioned fluttering method that disappears disappears after trembling, but related to the processing of frequency division and sequence comparison, mathematical regression model due to the method, calculation of complex, and need huge sample data amount to carry out, so the computing capability to processor CPU also has higher requirements, real-time is not high, and practical application is extremely restricted.
Summary of the invention
In view of this, the present invention proposes a kind of method of utilizing time window to eliminate the pulse per second (PPS) shake, the method realizes simple, can be good at eliminating the pulse per second (PPS) jitter error, and precision is several nanoseconds; And real-time is high.
For achieving the above object, technical scheme of the present invention comprises the steps:
Step 2, make i=1;
Step 3, after i pulse per second (PPS) after initial pulse per second (PPS) reaches, local crystal oscillator stop the counting, recording local crystal oscillator time counting is N
i, and with N
iBe stored as i the storing value that pulse per second (PPS) is corresponding;
Step 4, make local crystal oscillator restart counting, and make i from increasing 1;
Step 5, judge whether i<n sets up, return to step 3 if set up, if be false, enter step 6;
Step 6, after i pulse per second (PPS) after initial pulse per second (PPS) arrives, local crystal oscillator stops counting, recording local crystal oscillator time counting is N
i, the length take i pulse per second (PPS) as terminal point is in the time window of n, and the average of the corresponding count value of all pulse per second (PPS)s is
And calculate the deviate between each count value and average in this time window
Wherein m gets all over all integers in i-n+1 to i; Set deviation threshold, all are averaged less than count value corresponding to the deviate of deviation threshold, be designated as A
i, and with A
iBe stored as i the storing value that pulse per second (PPS) is corresponding;
Step 7, make local crystal oscillator restart counting, and make i from increasing 1;
Whether a series of pulse per second (PPS)s that step 8, judgement are launched launch end, if do not finish, return to step 6, if emission finishes, execution in step 9;
Step 9, a series of pulse per second (PPS)s of using the storing value simulation corresponding according to each pulse per second (PPS) of local crystal oscillator to launch.
Preferably, in above-mentioned steps 1 according to the precision set time window length of whole second time of required calculating.
Beneficial effect
1, the method is counted the interval of two adjacent pulse per second (PPS)s with local crystal oscillator, by count value is averaged calculating, reject with overall average and differ larger count value and again be averaged, the method realizes simple, can be good at eliminating the pulse per second (PPS) jitter error, and owing to selecting the quantity of count value according to the shake intensity of pulse per second (PPS), having realized adaptive effect, the precision that this method can reach is several nanoseconds.
2, this method is when select time window length, can carry out according to the precision of whole second time of required calculating that window is long to be adjusted, and wherein the long larger precision of window is higher, and the method arrives with pulse and be as the criterion, and real-time is high.
Description of drawings
Fig. 1 is this method flow chart.
Embodiment
Below in conjunction with the accompanying drawing specific embodiment that develops simultaneously, the present invention is described further.
There is the jitter error about 1 μ s in the GPS pulse per second (PPS) of adopting in measuring system, and local crystal oscillator clock frequency is 125MHz, and short-term stability is higher.It is the n time window of second that length is set, when arriving, each pulse per second (PPS) all use crystal oscillator to count, obtain count value, add up until after satisfying length of window Deng pulse per second (PPS), get the mean value of all pulse per second (PPS) count values in time window, only choose the count value that floats within the setting range of mean value up and down in time window, calculate the average of selected count value.Then time window moves successively backward with the arrival of new pulse per second (PPS), chooses count value and computation of mean values according to said method equally, and this moment, time window entered stable mobile running status until pulse per second (PPS) finishes.
In the present embodiment, adopt flow chart as shown in Figure 1 that one embodiment of the present invention are described:
Step 2, make i=1;
Step 3, after i pulse per second (PPS) after initial pulse per second (PPS) arrives, local crystal oscillator stops counting, recording local crystal oscillator time counting is N
i, and with N
iBe stored as i the storing value that pulse per second (PPS) is corresponding;
Step 4, local crystal oscillator restart counting, and make i from increasing 1, i.e. i=i+1;
Step 5, judge whether i<n sets up, returned to for the 3rd step if set up, repeat the 3rd step and the 4th to go on foot; If be false, entered for the 6th step;
Step 6, an i pulse per second (PPS) arrive, and local crystal oscillator stops counting, and recording local crystal oscillator time counting is N
i, this moment i>n, the length take i pulse per second (PPS) as terminal point is in the time window of n, the grand mean of all pulse per second (PPS)s correspondence count values is
Deviate corresponding to each count value is respectively
Deviate refers to each count value with respect to the deviation of average herein, and wherein m gets all over all values in i-n+1 to i.
Set deviation threshold, all are averaged less than the corresponding count value of the deviate of deviation threshold, be designated as A
i, store A
iThe deviation of calculating each count value and grand mean can show the degree of jitter of each pulse per second (PPS), and the more violent count value of shake is rejected in the setting of deviation threshold, has realized the purpose of debounce.And due at computation of mean values A
iThe time, therefore the count value of selection has reached adaptive purpose according to the difference of the shake intensity of pulse per second (PPS) in time window and quantity is different.
Step 7, local crystal oscillator restart counting, and i is from increasing 1, i.e. i=i+1;
Whether a series of pulse per second (PPS)s that step 8, judgement are launched launch end, if do not finish, return to for the 6th step, repeat for the 6th step and the 7th step, if the pulse per second (PPS) emission finishes, execution in step 9.
Step 9, a series of pulse per second (PPS)s of using the storing value simulation corresponding according to each pulse per second (PPS) of local crystal oscillator to launch.
The objective of the invention is in order to use local crystal oscillator that a series of pulse per second (PPS)s of launching are simulated, when simulation, local crystal oscillator is counted according to each pulse per second (PPS) corresponding stored value, namely uses for the pulse per second (PPS) less than the time window selected value N that stores
iSimulate, for choose greater than time window pulse per second (PPS) use the A store
iSimulate, shake has been eliminated in the pulse per second (PPS) that this moment, local crystal oscillator was simulated.
The present invention can be applicable to eliminate time error, calibration clock cycle and definite precise time benchmark in time measurement system, and the extensible high accuracy range measurement field of being applied to.
In sum, these are only preferred embodiment of the present invention, is not for limiting protection scope of the present invention.Within the spirit and principles in the present invention all, any modification of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.
Claims (2)
1. utilize time window to eliminate the method for pulse per second (PPS) shake, it is characterized in that, concrete steps are:
Step 1, setting-up time window length are n; Begin to receive pulse per second (PPS), after initial pulse per second (PPS) arrived, local crystal oscillator began counting;
Step 2, make i=1;
Step 3, after i pulse per second (PPS) after initial pulse per second (PPS) reaches, local crystal oscillator stop the counting, recording local crystal oscillator time counting is N
i, and with N
iBe stored as i the storing value that pulse per second (PPS) is corresponding;
Step 4, make local crystal oscillator restart counting, and make i from increasing 1;
Step 5, judge whether i<n sets up, return to step 3 if set up, if be false, enter step 6;
Step 6, after i pulse per second (PPS) after initial pulse per second (PPS) arrives, local crystal oscillator stops counting, recording local crystal oscillator time counting is N
i, the length take i pulse per second (PPS) as terminal point is in the time window of n, and the average of the corresponding count value of all pulse per second (PPS)s is
And calculate the deviate between each count value and average in this time window
Wherein m gets all over all integers in i-n+1 to i; Set deviation threshold, all are averaged less than count value corresponding to the deviate of deviation threshold, be designated as A
i, and with A
iBe stored as i the storing value that pulse per second (PPS) is corresponding;
Step 7, make local crystal oscillator restart counting, and make i from increasing 1;
Whether a series of pulse per second (PPS)s that step 8, judgement are launched launch end, if do not finish, return to step 6, if emission finishes, execution in step 9;
Step 9, a series of pulse per second (PPS)s of using the storing value simulation corresponding according to each pulse per second (PPS) of local crystal oscillator to launch.
2. the method for utilizing time window to eliminate the pulse per second (PPS) shake as claimed in claim 1 is characterized in that: according to the precision set time window length of whole second time of required calculating.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102694489A CN103384145A (en) | 2013-06-28 | 2013-06-28 | Method for utilizing time window to eliminate pulse per second jitter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2013102694489A CN103384145A (en) | 2013-06-28 | 2013-06-28 | Method for utilizing time window to eliminate pulse per second jitter |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103384145A true CN103384145A (en) | 2013-11-06 |
Family
ID=49491860
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2013102694489A Pending CN103384145A (en) | 2013-06-28 | 2013-06-28 | Method for utilizing time window to eliminate pulse per second jitter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103384145A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800409A (en) * | 2017-11-16 | 2018-03-13 | 湖南工业大学 | Accumulation formula disturbing pulse filter method |
CN108459239A (en) * | 2018-03-22 | 2018-08-28 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Distributed traveling wave monitoring terminal time deviation calculating and compensation method and system |
CN110782709A (en) * | 2019-11-04 | 2020-02-11 | 四川九洲空管科技有限责任公司 | High-precision clock redundancy backup method for civil aviation ADS-B ground station system |
CN114071694A (en) * | 2021-11-17 | 2022-02-18 | 成都天奥电子股份有限公司 | Time service precision improving method and device based on Beidou satellite signals and storage medium |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070127836A1 (en) * | 2005-12-07 | 2007-06-07 | Micron Technology, Inc. | Method and apparatus providing noise reduction while preserving edges for imagers |
CN101179283A (en) * | 2007-09-28 | 2008-05-14 | 杭州国芯科技有限公司 | Method of detecting and suppressing pulse noise |
CN101477801A (en) * | 2009-01-22 | 2009-07-08 | 东华大学 | Method for detecting and eliminating pulse noise in digital audio signal |
CN102538691A (en) * | 2011-12-31 | 2012-07-04 | 太原理工大学 | Method for measuring diameter of high-temperature steel pipe by using image recognition technology |
-
2013
- 2013-06-28 CN CN2013102694489A patent/CN103384145A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070127836A1 (en) * | 2005-12-07 | 2007-06-07 | Micron Technology, Inc. | Method and apparatus providing noise reduction while preserving edges for imagers |
US20080043124A1 (en) * | 2005-12-07 | 2008-02-21 | Micron Technology, Inc. | Methods and apparatuses providing noise reduction while preserving edges for imagers |
CN101179283A (en) * | 2007-09-28 | 2008-05-14 | 杭州国芯科技有限公司 | Method of detecting and suppressing pulse noise |
CN101477801A (en) * | 2009-01-22 | 2009-07-08 | 东华大学 | Method for detecting and eliminating pulse noise in digital audio signal |
CN102538691A (en) * | 2011-12-31 | 2012-07-04 | 太原理工大学 | Method for measuring diameter of high-temperature steel pipe by using image recognition technology |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800409A (en) * | 2017-11-16 | 2018-03-13 | 湖南工业大学 | Accumulation formula disturbing pulse filter method |
CN107800409B (en) * | 2017-11-16 | 2020-11-06 | 湖南工业大学 | Accumulation type interference pulse filtering method |
CN108459239A (en) * | 2018-03-22 | 2018-08-28 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | Distributed traveling wave monitoring terminal time deviation calculating and compensation method and system |
CN110782709A (en) * | 2019-11-04 | 2020-02-11 | 四川九洲空管科技有限责任公司 | High-precision clock redundancy backup method for civil aviation ADS-B ground station system |
CN114071694A (en) * | 2021-11-17 | 2022-02-18 | 成都天奥电子股份有限公司 | Time service precision improving method and device based on Beidou satellite signals and storage medium |
CN114071694B (en) * | 2021-11-17 | 2023-12-15 | 成都天奥电子股份有限公司 | Time service precision improving method and device based on Beidou satellite signals and storage medium |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Improved Constraints on D 0− D¯ 0 Mixing in D 0→ K+ π− Decays from the Belle Detector | |
CN106126961B (en) | Grey population satellite clock correction adaptive forecasting method and system | |
CN103384145A (en) | Method for utilizing time window to eliminate pulse per second jitter | |
CN107037721A (en) | One kind reviews one's lessons by oneself eurymeric time-to-digit converter | |
Rodriguez-Fuentes et al. | MediaEval 2013 spoken web search task: system performance measures | |
CN103364770B (en) | Radar target detecting system based on matrix filling and detecting method thereof | |
CN106681126B (en) | A kind of time-to-digit converter and its apparatus and method that calibrate for error | |
Cardone et al. | Testing the distance duality relation with present and future data | |
CN106597097A (en) | High-precision frequency measurement method | |
CN101788305A (en) | Method for rapid field calibration of micro inertial measurement unit | |
CN104199054A (en) | Preprocessing method for common view data of Beidou satellite navigation system | |
RU2395102C1 (en) | Method of measuring missile velocity and device for realising said method | |
CN108469609A (en) | A kind of detection information filtering method for radar target tracking | |
CN106647223A (en) | Quick stable real-time adjustment method for atomic clock timing | |
CN116321423A (en) | UWB NLOS propagation error suppression method based on deep learning | |
Guan et al. | Research on indoor positioning algorithm based on trilateral positioning and Taylor series expansion | |
Li et al. | Information theory in nonlinear error growth dynamics and its application to predictability: Taking the Lorenz system as an example | |
CN107300688A (en) | A kind of clock frequency Calibration Method in multipoint location system | |
CN103529687B (en) | Pulse interval measuring device | |
CN103575224A (en) | Principle for measuring log diameter class in double laser symmetric positioning microwave mode | |
CN116794615A (en) | Airborne radar precision evaluation analysis method | |
Shi et al. | Research on validation method for complex simulation systems | |
Shi et al. | Fmea-based control mechanism for embedded control software | |
Jansson | Improvements for IceCube Event Reconstruction through Geometry Calibration and Photon Timing Distributions | |
Chen et al. | Real-time intersection strategy of angular measuring data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20131106 |