CN104142624A - Time synchronization method and system based on waveform matching - Google Patents
Time synchronization method and system based on waveform matching Download PDFInfo
- Publication number
- CN104142624A CN104142624A CN201410403983.3A CN201410403983A CN104142624A CN 104142624 A CN104142624 A CN 104142624A CN 201410403983 A CN201410403983 A CN 201410403983A CN 104142624 A CN104142624 A CN 104142624A
- Authority
- CN
- China
- Prior art keywords
- data
- time
- matched
- bias
- reference data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a time synchronization method and system based on waveform matching. According to the method and system, based on the correlation principle, accurate time information is synchronized for a set of data not containing time scales or containing wrong time scales by utilizing a set of data containing reference time. The method comprises the steps that firstly, the data to be matched are segmented, correlation solution is performed on all segments and different sub-segments of reference data, when the correlation coefficient is maximum, the time difference between the segment of data and the reference data is calculated, and therefore time difference sequences of all the segments are obtained; secondly, the time difference sequences serve as observed values, time differences and time-scale drift factors serve as parameters to be estimated, and least square adjustment is performed to figure out the parameters to be estimated; finally, accurate time is synchronized for the data to be matched through the time differences and the time-scale drift factors. The time synchronization method and system can be widely applied to time synchronization during multi-sensor data fusion and can be further applied to correctness detection for time synchronization based on other methods.
Description
Technical field
The present invention relates to sensor application field, concrete is a kind of method for synchronizing time and system based on Waveform Matching.
Background technology
At present, sensor Data Fusion research is in the ascendant, and some proven technique has been applied in Practical Project, and has obtained good effect.The application of sensor Data Fusion not only can improve precision and the reliability of system, can also improve the measurement range of system, increases the confidence level of system, shortens the reaction time of system.But Fusion is a complicated uncertain information processing procedure, the precondition that can merge is that the information obtaining from each sensor is must be the description to the synchronization of same target.This comprises two aspects, first will guarantee that the information that each sensor obtains is the description that the same parameter of same target is carried out.The temporal information that secondly, guarantee the data that merge is synchronous.In dynamic duty environment, it is particularly outstanding that time synchronization problem shows.
In different engineering practices, there is the method for synchronizing time for particular problem.For example, utilize the method for curve to carry out the algorithm of time synchronized, while utilizing the method for sequence to solve, become the Time synchronization algorithm of observation, and adopt smothing filtering algorithm, by the methods such as time synchronized corresponding to the measurement data between each sensor.Before these algorithm application, each sensor is all stamped identical markers when starting to measure, and carries out timing afterwards according to sampling rate separately.Yet, for whole measuring process, all there is not identical markers in each sensor, or the processor that itself carries out timing the phenomenon of frequency marking drift occurs because of the impact of the factors such as temperature characterisitic, and above algorithm all can not correctly be realized the time synchronized between multisensor.
The more typical example of Fusion is exactly the integrated navigation system that GPS (Global Positioning System, GPS) and INS (Inertial Navigation System, inertial navigation system) form.The atomic clock that wherein GPS is equipped with degree of precision provides precise time for it, and INS can only carry out counter timing by the frequency of setting, and drift phenomenon often occurs INS frequency marking.GPS/INS integrated navigation system adopts Kalman filtering to carry out data fusion conventionally, and while only having GPS and INS subsystem data to put at one time, integrated navigation just has actual meaning.The most general, the effective practice is to utilize GPS pps pulse per second signal as the time synchronized benchmark of integrated navigation system at present, and integrated structure design realizes time synchronized.Yet when a IMU newly going out (Inertial Measurement Unit, Inertial Measurement Unit) need carry out integrated navigation test, though said method precision is high, complexity is high, cost is large, needs software and hardware support, inapplicable.In actual test process, the integrated navigation system that often carries one group of degree of precision is as with reference to system, and INS has wherein stamped high-precision GPS markers.
INS markers in integrated navigation system can be used as reference time system, only need the suitable algorithm of design, can realize the time synchronized between IMU to be tested and frame of reference, this algorithm is not only applicable to the time synchronized between inertial sensor, also can be applicable to the time synchronized between other sensors with reference time scale.
Summary of the invention
The present invention is directed to the problems referred to above, designed a kind of method for synchronizing time and system based on Waveform Matching; The method and system be to the whether identical not restriction of the original time system of different sensors data, and simple to operate, computing is quick, cost is low;
Technical scheme of the present invention is: a kind of method for synchronizing time based on Waveform Matching, comprises the steps:
Step 1: read reference data R
0and data T to be matched
0;
Step 2: unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1;
Step 3: respectively to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
Step 4: the thick mistiming int_bias according to trying to achieve, makes reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2; Enter step 5;
Step 5: unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3;
Step 6: try to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k.
Step 7: according to best bias and k combination, upgrade data T to be matched
0time, and storage;
Step 8: finish.
Unified reference data R in described step 2
0and data T to be matched
0the formula of data processing form as follows:
v=Δs/Δt
s
Wherein, Δ s is gyro or the measured increment size of accelerometer between former and later two samplings, Δ t
sfor the time interval between former and later two samplings, v is the rate value that Δ s is corresponding;
Reference data R after described conversion
1and data T to be matched
1that the time row after changing and an axis data of carrying out Waveform Matching form;
In described step 3, comprise the steps:
Step 3.1: respectively to R
1and T
1down-sampled, obtain the reference data R after down-sampled
2and data T to be matched
2;
Step 3.2: to down-sampled rear R
2and T
2among shorter sequence trailing zero, until both are equal in length;
Step 3.3: calculate R
2with T
2simple crosscorrelation sequence; Simple crosscorrelation sequence computing formula is as follows:
Wherein, x and y are the sequence of carrying out correlativity processing, and x and y length are N, and K and m are respectively used to represent the ordinal number of the value in sequence x and y, m value is-(N-1) to+(N-1);
Step 3.4: get simple crosscorrelation sequence
corresponding interval during maximal value m, obtains thick mistiming int_bias now, and the thick mistiming can utilize following formula to try to achieve:
int_bias=R
2(t
0)-T
2(t
0)-T
2_interval×m
Wherein, R
2(t
0) be R
2start time, T
2(t
0) be T
2start time, T
2_ interval is the time span of the interval representative that staggers, is T
2sampling interval, T
2_ interval * m is the time that two column datas stagger;
In described step 6.15, comprise following sub-step:
Step 6.1: calculate sub-sampling T_inter interval time;
Step 6.2: take T_inter as factor, calculate the data to be matched δ s of translation distance to the right;
Step 6.3: judge whether δ s is the positive integer times of sampling interval; If so, forward step S6.15 to, if not, forward step S6.4 to;
Step 6.4: by data data T to be matched
3translation δ s obtains T
4;
Step 6.5: take Δ t as length of window, by T
4be divided into limited n section, get present segment data T
4(i), 1≤i≤n;
Step 6.6: judge whether data to be matched extract n section; If so, forward step S6.13 to; If not, forward step S6.7 to;
Step 6.7: set moving window W
1, window is long is Δ t, starting point is the present segment zero hour;
Step 6.8: judge whether window sliding distance is greater than preseting length σ
s; As if so, forwarded step S6.11 to; Otherwise, with W
1for window, get reference data R
3(W
1);
Step 6.9: calculate T
4and R (i)
3(W
1) between Pearson's related coefficient; Pearson's Calculation of correlation factor as shown in the formula:
Wherein
for x serial mean,
for y serial mean, N is the length of sequence x and y, and r is the Pearson's related coefficient between sequence x and y;
Step 6.10: judge that whether reference data reads last position, if so, forwards step S6.11 to; Otherwise, by window W
1slide backward after the sampling interval of a reference data, forward step S6.8 to;
Step 6.11: calculate this section of mistiming Bias (i) and corresponding time time (i);
Step 6.12: extract next section of data to be matched.Take off one section of T
4(i), forward step 6.6 to;
Step 6.13: with the poor bias of least square adjustment seeking time and markers drift parameter k; Using bias and k as solve for parameter, and Bias (i) is observed reading, and observation equation is as follows:
Bias(i)=k×[time(i)-t
0]+bias
Wherein the initial value Bias0 of bias is Bias (i) sequence intermediate value, and the corresponding time is t
0, the initial value of k is 0;
Step 6.14: judge that whether δ s is greater than a sampling interval, if so, forwards step S6.15 to; Otherwise, by the δ s T_inter that doubles, forward step S6.3 to;
Step 6.15: select best bias and k combination; Utilize not bias and k on the same group, obtain the to be matched data of different time after synchronous;
In described step 6.15, if data are single shaft, that group bias while selecting covariance maximum and k are as best bias and k combination; If data are multiaxis, select the most believable axis data to carry out covariance and solve, select best bias and k combination.
In described step 7, utilize following formula to calculate and upgrade T
0time:
t′=(t-t
0)×k+bias+t
Wherein, t is T
0time row, t
0for time corresponding to mistiming initial value of selecting in step 6.12, k is markers drift parameter, and bias is the mistiming, and t ' is the rear time for t upgrades.
A clock synchronization system based on Waveform Matching, is characterized in that, comprises as lower module:
Read module: for reading reference data R
0and data T to be matched
0;
Data processing format module: for unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1;
Down-sampled module: for respectively to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
Thick mistiming module: for according to the thick mistiming int_bias trying to achieve, make reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2;
Sampling rate computing module: for unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3;
Computing module: for trying to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k;
Update module: for according to best bias and k combination, upgrade data T to be matched
0time, and storage.
Accompanying drawing explanation
Fig. 1 is schematic flow sheet of the present invention;
Fig. 2 is the schematic flow sheet of step 6 of the present invention;
Fig. 3-1st, the vertical deflection speed of EPSON gyro waveform schematic diagram;
Fig. 3-2nd, the vertical deflection speed of FSAS gyro waveform schematic diagram;
Fig. 4-1st, the maximum correlation coefficient schematic diagram that three axle segmentations of EPSON gyro are tried to achieve;
Fig. 4-2nd, the maximum correlation coefficient schematic diagram that three axle segmentations of EPSON accelerometer are tried to achieve;
Fig. 5-1st, the time scale difference sequence schematic diagram that three axle segmentations of EPSON gyro solve;
Fig. 5-2nd, the time scale difference sequence schematic diagram that three axle segmentations of EPSON accelerometer solve;
Fig. 6 is the vertical direction time match of gyro result schematic diagram;
Fig. 7 is the vertical direction time match of gyro result partial schematic diagram;
Fig. 8 is system architecture schematic diagram of the present invention.
Embodiment
Below in conjunction with drawings and Examples, the present invention is described in further detail.
Embodiment describes with the data of inertial sensor collection, and inertial sensor can gather multiaxis data, and the present embodiment take multiaxis data processing as example.
Data Source in embodiment: when vehicle-mounted test, carry low precision IMU EPSON and the GPS/INS integrated navigation system that contains high precision IMU FSAS simultaneously, in the data that wherein EPSON collects with counter timing, in integrated navigation system, FSAS process and GPS integrated design, stamped accurate gps time.
Data layout in embodiment: EPSON sampling rate is 125Hz, data are containing 7 row, wherein first classify time row as, 2nd~4 classify the three row angular speeds that gyroscope gathers as, be respectively forward direction, dextrad and vertical, 5th~7 classify the three row acceleration that accelerometer gathers as, are respectively forward direction, dextrad and vertical; FSAS sampling rate is 200Hz, and first classifies time row as, and 2nd~4 classify the three row angular velocity increments that gyroscope gathers as, are respectively forward direction, dextrad and vertical, and 5th~7 classify the three row acceleration increments that accelerometer gathers as, are respectively forward direction, dextrad and vertical.
The present invention can be without time mark or time and marks the synchronous precise time information of inaccurate data.As shown in Figure 1, Fig. 2 is the poor and markers drift parameter solution procedure of optimal time in Fig. 1 to implementation step.
S1: read reference data R
0and data T to be matched
0.Wherein FSAS data are as with reference to data R
0, EPSON data are as data T to be matched
0.
S2: unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1(useful data row); According to existing prior imformation, judgement R
0with T
0whether represented parameter is identical.If different, according to the mutual relationship between parameter, change, finally with frequently-used data form, be as the criterion.In this example, EPSON is rate form, and FSAS is increment form, and conventional data processing form is rate format.It is as follows that increment is converted to the formula of speed:
v=Δs/Δt
s
Wherein, Δ s is gyro or the measured increment size of accelerometer between former and later two samplings, Δ t
sfor the time interval between former and later two samplings, v is the rate value that Δ s is corresponding.
Utilize above formula to R
0change the average that the corresponding time is former and later two sampling times.If carry out Waveform Matching with an axis data wherein, extract time row after conversion and form new reference data R when axis data
1.In this example, the data of 6 axles are all carried out to Waveform Matching, final FSAS data R
1for R
0data after transforming, EPSON data T
1for T
0.
As Fig. 3-1 and Fig. 3-2 are depicted as R
1with T
1the vertical gyro data waveform of the 4th row, as can be seen from Figure, two waveforms have high consistency, and other axles also have same characteristic features, can carry out Waveform Matching according to the related coefficient of two data like this.
S3: to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
First, treat matched data and reference data and carry out down-sampled processing, down-sampled rear frequency is the common factor of two groups of data original frequencies, meeting on the basis of thick mistiming precision, the smaller the better.In embodiment by R
1and T
1down-sampled is 5Hz, to down-sampled rear FSAS data R
2and EPSON data T
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now.Concrete steps are as follows:
S3.1: to R
1and T
1down-sampled, obtain the reference data R after down-sampled
2and data T to be matched
2;
S3.2: at R
2and T
2shorter sequence trailing zero is until both are equal in length;
S3.3: calculate R
2with T
2simple crosscorrelation sequence.It is as follows that simple crosscorrelation sequence is calculated General Principle formula:
Wherein x and y are the sequence of carrying out correlativity processing, and x and y length are N, and subscript K and m be for representing the ordinal number of value of sequence x and y, m value be-(N-1) arrive+(N-1).
for simple crosscorrelation sequence,
while being less than zero for m, utilize m to be greater than the simple crosscorrelation sequence that zero formula calculates.When m get-(N-1) to+in (N-1) during different value, all can utilize above formula to calculate and try to achieve a cross correlation value
finally, obtain the cross correlation function sequential value that vector length is 2 * N-1.
Step 3.4: get simple crosscorrelation sequence
corresponding interval during maximal value m, obtains thick mistiming int_bias now, in this example, and the R after end zero-padded length equates
2with T
2the 4th row be sequence x and the y that need carry out correlativity processing, utilize above formula to calculate and try to achieve simple crosscorrelation sequential value
corresponding m when value is maximum, is two column data correlativitys when maximum, the interval of staggering.The thick mistiming can utilize following formula to try to achieve:
int_bias=R
2(t
0)-T
2(t
0)-T
2_interval×m
R
2(t
0) be R
2start time, T
2(t
0) be T
2start time, T
2_ interval is the time span of the interval representative that staggers, is T in this example
2sampling interval, T
2_ interval * m is the time that two column datas stagger.
S4: the thick mistiming int_bias according to trying to achieve, makes reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2.The object of this step is according to the thick mistiming, finds data T to be matched
2with reference data R
2between the corresponding time period, make reference data R
2comprise data T to be matched
2, and front and back respectively exceed data T to be matched
2the principle of certain length, intercepting reference data R
2and data T to be matched
2.In the present embodiment, thick match time of the precision of trying to achieve is 1 second level, therefore when getting FSAS data, make respectively to exceed FSAS data 0.5 second before and after it, enters S5.
S5: unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3; FSAS sampling rate is 200Hz, and EPSON sampling rate is 125Hz, finally to T
0carry out linear interpolation, making its sampling rate is 200Hz.
Through above-mentioned in steps after, FSAS data are R
3, EPSON data are T
3.
S6: try to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k; Solving of mistiming and markers drift parameter is the core of this algorithm, and calculated amount is maximum, is also a most complicated step simultaneously.
S6.1: calculate sub-sampling T_inter interval time, with sampling interval 1/10~1/2 between be advisable.Sub-sampling interval is defined as the arbitrary value lower than sampling interval, in concrete enforcement, take user's accuracy requirement as criterion, suitably value.T_inter is too small, and precision is improved, but not raising of accuracy, and calculated amount increases; T_inter is excessive, does not reach and proposes high-precision requirement.For example, R in this example
3with T
3the sampling interval time be 0.005 second, sub-sampling interval T _ inter is defined as 1/5 second of sampling interval.
S6.2: take T_inter as factor, calculate the data to be matched δ s of translation distance to the right.At circulation time for the first time, δ s is 0 times of T_inter, i.e. δ s=0 * T_inter=0.
S6.3: judge whether δ s is the positive integer times of sampling interval.If so, forward step S6.15 to if not, forward step S6.4 to.
S6.4: by EPSON data translation δ s.First by T
3time row carry out to right translation, translation distance is δ s, after translation, the time be listed as interpolated point, the numerical value T that other column datas linearity interpolations are made new advances
4.For example there is original time t
1, t
2, corresponding numerical value is val
1, val
2, for new time t
3(t
1< t
3< t
2), corresponding numerical value val
3by following formula linear interpolation out:
S6.5: take Δ t as length of window, by T
4be divided into limited n section, get present segment data T
4(i), 1≤i≤n.In Δ t length of window, should be able to comprise abundant multidate information, to restrain the impact of random noise.In this example, Δ t is taken as 60 seconds, T
4be divided into 61 sections, get first paragraph T
4(1).
S6.6: judge whether data to be matched extract n section.If so, forward step S6.13 to.If not, forward step S6.7 to.
S6.7: set moving window W
1, window is long is Δ t, starting point is the present segment zero hour.It is 60 seconds that window length is grown the same with the window in S6.5, W
1at every turn along time shaft backward sampling interval of translation get reference data, window sliding distance is W
1backward translation T.T., account form is: W
1the number of translation sampling interval is backward multiplied by the sampling interval time.
S6.8: judge whether window sliding distance is greater than preseting length σ
s(σ
saccording to thick match time of precision and determining, be about twice thick match time of precision), as if so, forwarded step S6.11 to.Otherwise, with W
1for window, get reference data R
3(W
1).In this example, σ
sit is 1 second.
S6.9: calculate T
4and R (i)
3(W
1) between Pearson's related coefficient.Pearson's Calculation of correlation factor principle as shown in the formula:
Wherein
for x serial mean,
for y serial mean, N is the length of sequence x and y, and r is the Pearson's related coefficient between sequence x and y.In this example, T
4and R (i)
3(W
1) respectively containing six column datas, each respective column data are carried out separately related coefficient according to above formula and solved, each respective shaft all obtains a related coefficient, totally six.
S6.10: judge that whether reference data reads last position, if so, forwards step S6.11 to.Otherwise, by window W
1slide backward after the sampling interval of a reference data, forward step S6.8 to.
S6.11: calculate this section of mistiming Bias (i) and corresponding time time (i).For this section of EPSON data, W
1sliding process in, T
4and R (i)
3(W
1) all can obtain corresponding facies relationship ordered series of numbers between each respective column data.For one group of respective column data wherein, when related coefficient is maximum, two column data waveform optimum matching, maximum correlation coefficient max_corr now, W
1the sampling interval number sliding is to the right num.This axle can calculate a Bias (i) with following formula, and other axles in like manner.
Bias(i)=R
3(t
0)-T
4(t
0)+R
3_interval×num
R
3(t
0) be R
3start time, T
4(t
0) be T
4start time, R
3_ interval is W
1dextrad moves the time of a sampling interval representative, and this example is R
3sampling interval, R
3_ interval * num is W
1total duration slides to the right.
Calculate after the Bias (i) of six axles, each axle max_corr of take is weight, and weighted mean obtains T
4(i) unique Bias (i), the max_corr of this section (i) is the arithmetic mean of each axle max_corr, time (i) is T
4(i) time in the moment in the middle of.
Mistiming between this section of data to be matched and reference data has calculated complete.
S6.12: extract next section of data to be matched.Take off one section of T
4(i), forward step 6.6 to.
As shown in Fig. 4-1 and 4-2, be all sections of maximum correlation coefficients of trying to achieve of six axles of EPSON.Related coefficient is larger, and data waveform matching degree is higher.From figure, also can find out, after the 16th section, related coefficient, close to 1, show that EPSON and FSAS data waveform consistance are higher, and method identifies optimum matching section thus.
As shown in Fig. 5-1 and 5-2, the mistiming sequence being calculated by each section of maximum correlation coefficient, chooses the period of related coefficient when larger.In figure, can find out that the method can not only calculate the mistiming of EPSON, also can find that the frequency marking of EPSON exists drift.
S6.13: the poor bias of least square adjustment seeking time and markers drift parameter k.Using bias and k as solve for parameter, and Bias (i) is observed reading, and observation equation is as follows:
Bias(i)=k×[time(i)-t
0]+bias
Wherein the initial value Bias0 of bias is Bias (i) sequence intermediate value, and the corresponding time is t
0, the initial value of k is 0.
The power of Bias (i) is max_corr (i).For removing the rough error in observed reading, set a correlation coefficient threshold σ
r, in this example, σ
rbe 0.9.Only have when max_corr (i) is greater than 0.9, observed reading is substitution observation equation, carries out least square adjustment and solves, and tries to achieve one group of bias and k.
S6.14: judge that whether δ s is greater than a sampling interval, if so, forwards step S6.15 to.Otherwise, by the δ s T_inter that doubles, forward step S6.3 to.
S6.15: select best bias and k combination.Utilize not bias and k on the same group, obtain the to be matched data of different time after synchronous.Calculate the covariance between data to be matched and reference data.If data are single shaft, that group bias while selecting covariance maximum and k are as best bias and k combination; If when data are multiaxis, be as the criterion with practical experience, select the most believable axis data to carry out covariance and solve, select best bias and k combination.
S7: according to best bias and k combination, upgrade data T to be matched
0time, and storage.In this example, the most believable is the 4th column data, and after time synchronized, the 3rd group of covariance that data obtain is maximum.Utilize following formula to calculate and upgrade T
0time:
t′=(t-t
0)×k+bias+t
Wherein, t is T
0time row, t
0for time corresponding to mistiming initial value of selecting in step 6.12, k is markers drift parameter, and bias is the mistiming, and t ' is the rear time for t upgrades.
By T
0in t ' replacement for time row, and data are stored as and T again
0the data file that data layout is consistent.
S8: finish.
As shown in Figure 6, for utilizing algorithm of the present invention to carry out after time synchronized EPSON and FSAS, vertical gyrobearing is marked lower oscillogram at one time.From then on figure can find out, this algorithm can correctly carry out time synchronized computing.
As shown in Figure 7, be the Waveform Matching partial, detailed view in certain period in Fig. 6.In figure, can find out, this algorithm can not only correctly carry out the time synchronized between data to be matched and reference data, and synchronization accuracy is higher.
A clock synchronization system based on Waveform Matching, is characterized in that, comprises as lower module:
Read module: for reading reference data R
0and data T to be matched
0;
Data processing format module: for unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1;
Down-sampled module: for respectively to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
Thick mistiming module: for according to the thick mistiming int_bias trying to achieve, make reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2;
Sampling rate computing module: for unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3;
Computing module: for trying to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k;
Update module: for according to best bias and k combination, upgrade data T to be matched
0time, and storage.
Claims (7)
1. the method for synchronizing time based on Waveform Matching, is characterized in that, comprises the steps:
Step 1: read reference data R
0and data T to be matched
0;
Step 2: unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1;
Step 3: respectively to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
Step 4: the thick mistiming int_bias according to trying to achieve, makes reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2, enter step 5;
Step 5: unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3;
Step 6: try to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k;
Step 7: according to best bias and k combination, upgrade data T to be matched
0time, and storage;
Step 8: finish.
2. a kind of method for synchronizing time based on Waveform Matching according to claim 1, is characterized in that: unified reference data R in described step 2
0and data T to be matched
0the formula of data processing form as follows:
v=Δs/Δt
s
Wherein, Δ s is gyro or the measured increment size of accelerometer between former and later two samplings, Δ t
sfor the time interval between former and later two samplings, v is the rate value that Δ s is corresponding;
Reference data R after described conversion
1and data T to be matched
1that the time row after changing and an axis data of carrying out Waveform Matching form.
3. a kind of method for synchronizing time based on Waveform Matching according to claim 1, is characterized in that: in described step 3, comprise the steps:
Step 3.1: respectively to R
1and T
1down-sampled, obtain the reference data R after down-sampled
2and data T to be matched
2;
Step 3.2: to down-sampled rear R
2and T
2among shorter sequence trailing zero, until both are equal in length;
Step 3.3: calculate R
2with T
2simple crosscorrelation sequence; Simple crosscorrelation sequence computing formula is as follows:
Wherein, x and y are the sequence of carrying out correlativity processing, and x and y length are N, and K and m are respectively used to represent the ordinal number of the value in sequence x and y, m value is-(N-1) to+(N-1);
Step 3.4: get simple crosscorrelation sequence
corresponding interval during maximal value m, obtains thick mistiming int_bias now, and the thick mistiming can utilize following formula to try to achieve:
int_bias=R
2(t
0)-T
2(t
0)-T
2_interval×m
Wherein, R
2(t
0) be R
2start time, T
2(t
0) be T
2start time, T
2_ interval is the time span of the interval representative that staggers, is T
2sampling interval, T
2_ interval * m is the time that two column datas stagger.
4. a kind of method for synchronizing time based on Waveform Matching according to claim 1, is characterized in that: in described step 6.15, comprise following sub-step:
Step 6.1: calculate sub-sampling T_inter interval time;
Step 6.2: take T_inter as factor, calculate the data to be matched δ s of translation distance to the right;
Step 6.3: judge whether δ s is the positive integer times of sampling interval; If so, forward step S6.15 to, if not, forward step S6.4 to;
Step 6.4: by data data T to be matched
3translation δ s obtains T
4;
Step 6.5: take Δ t as length of window, by T
4be divided into limited n section, get present segment data T
4(i), 1≤i≤n;
Step 6.6: judge whether data to be matched extract n section; If so, forward step S6.13 to; If not, forward step S6.7 to;
Step 6.7: set moving window W
1, window is long is Δ t, starting point is the present segment zero hour;
Step 6.8: judge whether window sliding distance is greater than preseting length σ
s; As if so, forwarded step S6.11 to;
Otherwise, with W
1for window, get reference data R
3(W
1);
Step 6.9: calculate T
4and R (i)
3(W
1) between Pearson's related coefficient; Pearson's Calculation of correlation factor as shown in the formula:
Wherein
for x serial mean,
for y serial mean, N is the length of sequence x and y, and r is the Pearson's related coefficient between sequence x and y;
Step 6.10: judge that whether reference data reads last position, if so, forwards step S6.11 to; Otherwise,
By window W
1slide backward after the sampling interval of a reference data, forward step S6.8 to;
Step 6.11: calculate this section of mistiming Bias (i) and corresponding time time (i);
Step 6.12: extract next section of data to be matched.Take off one section of T
4(i), forward step 6.6 to;
Step 6.13: with the poor bias of least square adjustment seeking time and markers drift parameter k; Using bias and k as solve for parameter, and Bias (i) is observed reading, and observation equation is as follows:
Bias(i)=k×[time(i)-t
0]+bias
Wherein the initial value Bias0 of bias is Bias (i) sequence intermediate value, and the corresponding time is t
0, the initial value of k is 0;
Step 6.14: judge that whether δ s is greater than a sampling interval, if so, forwards step S6.15 to; Otherwise, by the δ s T_inter that doubles, forward step S6.3 to;
Step 6.15: select best bias and k combination; Utilize not bias and k on the same group, obtain the to be matched data of different time after synchronous.
5. a kind of method for synchronizing time based on Waveform Matching according to claim 4, is characterized in that: in described step 6.15, if data are single shaft, that group bias while selecting covariance maximum and k are as best bias and k combination; If data are multiaxis, select the most believable axis data to carry out covariance and solve, select best bias and k combination.
6. a kind of method for synchronizing time based on Waveform Matching according to claim 1, is characterized in that: in described step 7, utilize following formula to calculate and upgrade T
0time:
t′=(t-t
0)×k+bias+t
Wherein, t is T
0time row, t
0for time corresponding to mistiming initial value of selecting in step 6.12, k is markers drift parameter, and bias is the mistiming, and t ' is the rear time for t upgrades.
7. the clock synchronization system based on Waveform Matching, is characterized in that, comprises as lower module:
Read module: for reading reference data R
0and data T to be matched
0;
Data processing format module: for unified reference data R
0and data T to be matched
0data processing form, extract the reference data R after conversion
1and data T to be matched
1;
Down-sampled module: for respectively to R
1and T
1down-sampled, to the reference data R after down-sampled
2and data T to be matched
2carry out simple crosscorrelation sequence and solve, corresponding interval while getting simple crosscorrelation sequence maximal value, obtains thick mistiming int_bias now;
Interception module: for according to the thick mistiming int_bias trying to achieve, make reference data R
2time period cover data T to be matched
2time period, and front and back respectively exceed data T to be matched
2certain length, intercepting reference data R
2and data T to be matched
2;
Sampling rate computing module: for unified R
2and T
2sampling rate, obtains the reference data R of identical sampling rate
3with data T to be matched
3;
Computing module: for trying to achieve R
3and T
3the poor bias of optimal time and markers drift parameter k;
Update module: for according to best bias and k combination, upgrade data T to be matched
0time, and storage.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403983.3A CN104142624B (en) | 2014-08-15 | 2014-08-15 | Time synchronization method and system based on waveform matching |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410403983.3A CN104142624B (en) | 2014-08-15 | 2014-08-15 | Time synchronization method and system based on waveform matching |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104142624A true CN104142624A (en) | 2014-11-12 |
CN104142624B CN104142624B (en) | 2015-04-01 |
Family
ID=51851840
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410403983.3A Active CN104142624B (en) | 2014-08-15 | 2014-08-15 | Time synchronization method and system based on waveform matching |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104142624B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226836A (en) * | 2016-12-31 | 2018-06-29 | 科大智能电气技术有限公司 | A kind of calibration method of fault detector synchronous waveform |
CN110477895A (en) * | 2019-07-24 | 2019-11-22 | 苏州国科医疗科技发展有限公司 | The continuous method for measuring heart rate of multiple light courcess detector based on volumetric blood wave |
CN110703586A (en) * | 2019-09-09 | 2020-01-17 | 广州市中海达测绘仪器有限公司 | Time synchronization method, data synchronization method, device, system, equipment and medium |
CN110895321A (en) * | 2019-12-06 | 2020-03-20 | 南京南瑞继保电气有限公司 | Secondary equipment time mark alignment method based on recording file reference channel |
CN112181906A (en) * | 2019-07-02 | 2021-01-05 | 广州汽车集团股份有限公司 | Test data synchronization method, device and system based on combined data acquisition instrument |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7752483B1 (en) * | 2006-12-13 | 2010-07-06 | Science Applications International Corporation | Process and system for three-dimensional urban modeling |
CN102183253A (en) * | 2010-12-31 | 2011-09-14 | 北京航空航天大学 | Software time synchronization method for position and orientation system |
CN203204161U (en) * | 2013-05-03 | 2013-09-18 | 武汉大学 | Knocking-type pulse-trigger INS/GPS time synchronizer |
-
2014
- 2014-08-15 CN CN201410403983.3A patent/CN104142624B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7752483B1 (en) * | 2006-12-13 | 2010-07-06 | Science Applications International Corporation | Process and system for three-dimensional urban modeling |
CN102183253A (en) * | 2010-12-31 | 2011-09-14 | 北京航空航天大学 | Software time synchronization method for position and orientation system |
CN203204161U (en) * | 2013-05-03 | 2013-09-18 | 武汉大学 | Knocking-type pulse-trigger INS/GPS time synchronizer |
Non-Patent Citations (3)
Title |
---|
刘占超 等: "一种改进的高精度POS时间同步方法", 《仪器仪表学报》, vol. 32, no. 10, 31 October 2011 (2011-10-31) * |
张涛 等: "不同动态条件下组合导航系统的时间同步", 《中国惯性技术学报》, vol. 20, no. 3, 30 June 2012 (2012-06-30), pages 320 - 325 * |
李倩 等: "GPS/INS组合导航系统时间同步系统设计", 《传感技术学报》, vol. 22, no. 12, 31 December 2009 (2009-12-31), pages 1752 - 1756 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108226836A (en) * | 2016-12-31 | 2018-06-29 | 科大智能电气技术有限公司 | A kind of calibration method of fault detector synchronous waveform |
CN108226836B (en) * | 2016-12-31 | 2021-06-01 | 科大智能电气技术有限公司 | Calibration method for waveform synchronization of fault indicator |
CN112181906A (en) * | 2019-07-02 | 2021-01-05 | 广州汽车集团股份有限公司 | Test data synchronization method, device and system based on combined data acquisition instrument |
CN110477895A (en) * | 2019-07-24 | 2019-11-22 | 苏州国科医疗科技发展有限公司 | The continuous method for measuring heart rate of multiple light courcess detector based on volumetric blood wave |
CN110703586A (en) * | 2019-09-09 | 2020-01-17 | 广州市中海达测绘仪器有限公司 | Time synchronization method, data synchronization method, device, system, equipment and medium |
CN110703586B (en) * | 2019-09-09 | 2022-02-18 | 广州市中海达测绘仪器有限公司 | Time synchronization method, data synchronization method, device, system, equipment and medium |
CN110895321A (en) * | 2019-12-06 | 2020-03-20 | 南京南瑞继保电气有限公司 | Secondary equipment time mark alignment method based on recording file reference channel |
CN110895321B (en) * | 2019-12-06 | 2021-12-10 | 南京南瑞继保电气有限公司 | Secondary equipment time mark alignment method based on recording file reference channel |
Also Published As
Publication number | Publication date |
---|---|
CN104142624B (en) | 2015-04-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104142624B (en) | Time synchronization method and system based on waveform matching | |
CN107797129B (en) | Point cloud data acquisition method and device under no GNSS signal | |
CN103630137B (en) | A kind of for the attitude of navigational system and the bearing calibration of course angle | |
CN102435763B (en) | Measuring method for attitude angular velocity of spacecraft based on star sensor | |
CN105509738B (en) | Vehicle positioning orientation method based on inertial navigation/Doppler radar combination | |
JP6019504B2 (en) | Traveling direction estimation apparatus and traveling direction estimation method for moving body | |
CN105910602B (en) | A kind of Combinated navigation method | |
US20150153178A1 (en) | Car navigation system and method in which global navigation satellite system (gnss) and dead reckoning (dr) are merged | |
CN104061934A (en) | Pedestrian indoor position tracking method based on inertial sensor | |
CN109507706B (en) | GPS signal loss prediction positioning method | |
CN103245797A (en) | State detection device, electronic apparatus, and program | |
CN107024212A (en) | A kind of astronomical doppler combined navigation method of deep space probe X-ray pulsar/time difference | |
CN106500693A (en) | A kind of AHRS algorithms based on adaptive extended kalman filtering | |
CN107747953A (en) | A kind of multi-sensor data and orbit information method for synchronizing time | |
CN103604430A (en) | Marginalized cubature Kalman filter (CKF)-based gravity aided navigation method | |
CN103900574A (en) | Attitude estimation method based on iteration volume Kalman filter | |
CN102878989A (en) | Triaxial angular vibration measuring method through adopting satellite-borne linear accelerometers | |
CN104613965A (en) | Stepping type pedestrian navigation method based on bidirectional filtering smoothing technology | |
CN104359496A (en) | High-precision attitude correction method based on vertical deviation compensation | |
CN103395435A (en) | Real-time positioning system method of high-precision and high-speed train | |
CN102967868A (en) | Positioning device and signal processing method thereof | |
JP5664059B2 (en) | Vehicle trajectory estimation device | |
CN103364842B (en) | A kind of error separation method of strapdown airborne gravitometer | |
CN105547291A (en) | Self-adaptive static detection method for indoor personnel self-localization system | |
CN109484434A (en) | Train integrated positioning device and method |
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 |