CN108680824A  Distributed wave recording synchronous method, device, equipment and medium  Google Patents
Distributed wave recording synchronous method, device, equipment and medium Download PDFInfo
 Publication number
 CN108680824A CN108680824A CN201810473612.0A CN201810473612A CN108680824A CN 108680824 A CN108680824 A CN 108680824A CN 201810473612 A CN201810473612 A CN 201810473612A CN 108680824 A CN108680824 A CN 108680824A
 Authority
 CN
 China
 Prior art keywords
 recording
 curve
 synchronized
 domain
 benchmark
 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
 230000001360 synchronised Effects 0.000 title claims abstract description 114
 230000001808 coupling Effects 0.000 claims abstract description 62
 238000010168 coupling process Methods 0.000 claims abstract description 62
 238000005859 coupling reaction Methods 0.000 claims abstract description 62
 230000004301 light adaptation Effects 0.000 claims abstract description 58
 238000002922 simulated annealing Methods 0.000 claims abstract description 25
 238000004590 computer program Methods 0.000 claims description 16
 238000005070 sampling Methods 0.000 claims description 16
 238000003860 storage Methods 0.000 claims description 14
 230000035852 Tmax Effects 0.000 claims description 12
 230000000875 corresponding Effects 0.000 claims description 12
 230000002441 reversible Effects 0.000 claims description 9
 230000001276 controlling effect Effects 0.000 claims 1
 238000004364 calculation method Methods 0.000 abstract description 11
 238000010586 diagram Methods 0.000 description 17
 238000004422 calculation algorithm Methods 0.000 description 7
 238000004891 communication Methods 0.000 description 5
 230000002493 climbing Effects 0.000 description 4
 238000009826 distribution Methods 0.000 description 3
 238000000034 method Methods 0.000 description 3
 238000001514 detection method Methods 0.000 description 2
 230000005611 electricity Effects 0.000 description 2
 238000005516 engineering process Methods 0.000 description 2
 239000004744 fabric Substances 0.000 description 2
 230000004048 modification Effects 0.000 description 2
 238000006011 modification reaction Methods 0.000 description 2
 238000004088 simulation Methods 0.000 description 2
 238000004458 analytical method Methods 0.000 description 1
 230000006399 behavior Effects 0.000 description 1
 230000005540 biological transmission Effects 0.000 description 1
 239000000969 carrier Substances 0.000 description 1
 238000006243 chemical reaction Methods 0.000 description 1
 235000019800 disodium phosphate Nutrition 0.000 description 1
 230000000694 effects Effects 0.000 description 1
 238000002156 mixing Methods 0.000 description 1
 230000035772 mutation Effects 0.000 description 1
 238000005457 optimization Methods 0.000 description 1
 238000002360 preparation method Methods 0.000 description 1
 239000002965 rope Substances 0.000 description 1
 239000007787 solid Substances 0.000 description 1
 230000001131 transforming Effects 0.000 description 1
 230000001052 transient Effects 0.000 description 1
 230000001960 triggered Effects 0.000 description 1
Classifications

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/08—Locating faults in cables, transmission lines, or networks
 G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
 G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/08—Locating faults in cables, transmission lines, or networks
 G01R31/088—Aspects of digital computing

 G—PHYSICS
 G01—MEASURING; TESTING
 G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
 G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
 G01R31/50—Testing of electric apparatus, lines, cables or components for shortcircuits, continuity, leakage current or incorrect line connections

 G—PHYSICS
 G06—COMPUTING; CALCULATING; COUNTING
 G06K—RECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
 G06K9/00—Methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
 G06K9/00496—Recognising patterns in signals and combinations thereof
 G06K9/00536—Classification; Matching
Abstract
The invention discloses a kind of distributed wave recording synchronous method, including：Obtain benchmark recording curve and recording curve to be synchronized；According to the benchmark recording curve and the recording curve acquisition domain to be synchronized；Continuous coupling function is obtained according to the domain and least square method；The Optimum Matching time of the continuous coupling function is obtained according to the simulated annealing traversed in advance based on the overall situation, under the premise of ensureing precision, reduce calculation amount, accelerate the speed of search optimal value, the acquisition of Optimum Matching time can be fast implemented, discretization is carried out to the continuous coupling function according to discrete integration formula, obtains discrete adaptation function；Realize that the benchmark recording curve is synchronous with the recording curve to be synchronized with the discrete adaptation function according to the Optimum Matching time.With minimum operand and space hold, the calculating of discrete adaptation function is realized, finally realizes the synchronization of different durations, sample frequency curve.
Description
Technical field
The present invention relates to field of power more particularly to distributed wave recording synchronous method, device, equipment and media.
Background technology
As intelligent power grid technology advances in big strides, power department wants the Distributed Application of fault wave recording device
Ask higher and higher, the requirement synchronized to recording in a distributed system is also increasingly tighter.
Fault oscillograph is failure forward and backward process can be automatically recorded when electric system is broken down and vibrated various electrical
Measure a kind of device of variation.It can be recorded because the large disturbances such as short trouble, system oscillation, collapse of frequency, collapse of voltage cause
System power, voltage and its derived quantity (such as active and reactive and system frequency) overall process variation.Mainly for detection of
The action behavior of relay protection and automatic safety device understands the changing rule of each electric parameter in system during system transient modelling,
And check the correctness etc. of power system computation program and model parameter.Currently, the recording of fault wave recording device is the result is that divide
Analyse the important evidence of electric power system fault.
When DC transmission system breaks down, concentration is passed to operations staff by the failure waverecording for being distributed in each device,
Each device recording has different a time span and sample rate, operations staff need manual sorting analysis failure waverecording etc. numerous from
Information is dissipated, can not realize quick positioning and solution to failure in a short time, to influence to restore electricity speed and power grid peace
Row for the national games, so needing to synchronize processing to distributed recording.
About the method that recording curve synchronizes, a kind of thinking is the Characteristics of Mutation amount of detection recording curve, this method profit
With the Partial Feature of curve rather than whole feature, in fact it could happen that the problem of timing error.
Invention content
In view of the abovementioned problems, the purpose of the present invention is to provide a kind of distributed wave recording synchronous method, record can be accurately realized
Curve line locking.
In a first aspect, the present invention provides a kind of distributed wave recording synchronous method, including：
Obtain benchmark recording curve and recording curve to be synchronized；
According to the benchmark recording curve and the recording curve acquisition domain to be synchronized；
Continuous coupling function is obtained according to the domain and least square method；
The Optimum Matching time of the continuous coupling function is obtained according to the simulated annealing traversed in advance based on the overall situation；
Discretization is carried out to the continuous coupling function according to discrete integration formula, obtains discrete adaptation function；
According to the Optimum Matching time and the discrete adaptation function realize the benchmark recording curve and it is described wait for it is same
Recording curve is walked to synchronize.
It is described according to the benchmark recording curve and described to be synchronized in the first possible realization method of first aspect
Recording curve acquisition domain includes：
Obtaining the benchmark recording curve and the recording curve to be synchronized has the condition of intersection；Wherein, described to have intersection
Condition beWherein, t_{b0}It is the left end point of the recording curve b to be synchronized, t_{bm}It is the record to be synchronized
The right endpoint of wave profile b, t_{aN}It is the right endpoint of the benchmark recording curve a, t_{a0}It is the left end point of the benchmark recording curve a,
τ is the Optimum Matching time；
Domain is obtained according to the condition for having intersection；Wherein, the domain includes Optimum Matching timedomain
With effective domain；Wherein, the Optimum Matching timedomain is：
τ∈(t_{a0}t_{bM},t_{aN}t_{b0})=(τ_{min},τ_{max})=T_{τ}
Effective domain is：t∈min(T_{a},T_{b})=(t_{min},t_{max})=T_{m}；Wherein, T_{a}For benchmark recording song
The domain of line a, T_{b}For the domain of the benchmark recording curve b.
The possible realization method of with reference to first aspect the first, in second of possible realization method of first aspect, institute
It states and includes according to the domain and least square method acquisition continuous coupling function：
Time reversal power is formed according to the time of coincidence length of the benchmark recording curve and the recording curve to be synchronized
Value；
The least square method with antitime weight is obtained according to the time reversal weights and least square method；
Continuous coupling function is obtained according to domain and the least square method of the antitime weight of the band
Wherein, τ is the Optimum Matching time, and a (t) is the benchmark recording curve, and b (t) is the recording curve to be synchronized,
t_{max}For the maximum value of the benchmark recording curve and the effective domain of recording curve to be synchronized, t_{min}The benchmark recording is bent
The minimum value of line and the effective domain of recording curve to be synchronized；τ∈T_{τ}, t ∈ T_{m}。
In the third possible realization method of first aspect, the basis is based on the global simulated annealing traversed in advance
The Optimum Matching time for obtaining the continuous coupling function includes:
The gamut lowdensity that the continuous coupling function is obtained according to the traversal of small decile n traverses curve；
When the continuous coupling functional value is smaller, the range of the corresponding lowdensity traversal curve is obtained；
In the range, the continuous coupling function is obtained most according to the temperature T faster simulated annealings that decay
Excellent match time.
In the 4th kind of possible realization method of first aspect, it is described according to discrete integration formula to the continuous coupling letter
Number carries out discretization, and obtaining discrete adaptation function includes：
Discretization is carried out to the continuous coupling function according to discrete integration formula, obtains the discrete adaptation functionWherein, g (t_{i})=a (t_{i})b(t_{i}+ τ), τ ∈ T_{τ},t∈T_{m}。
The 4th kind of possible realization method with reference to first aspect, in the 5th kind of possible realization method of first aspect, institute
It is trapezoidal integration formula to state discrete integration formula；
The trapezoidal integration formula isWherein, f (t_{i}) it is on trapezoidal
Bottom；f(t_{i+1}) it is trapezoidal bottom；(t_{i+1}t_{i}) it is trapezoidal height.
It is described according to the Optimum Matching time and discrete described in the 6th kind of possible realization method of first aspect
With function realize the benchmark recording curve and the recording curve to be synchronized synchronization include：
Obtain the former sampled point of the benchmark recording curve and the recording curve to be synchronized；
The discrete adaptation function is calculated respectively according to the former sampled point, obtains the discrete matching letter after operation
Number；
According to the discrete adaptation function after the Optimum Matching time and the operation realize the benchmark recording curve and
The recording curve to be synchronized synchronizes.
Second aspect, the present invention also provides a kind of distributed wave recording synchronizing devices, including：
Curve acquisition module, for obtaining benchmark recording curve and recording curve to be synchronized；
Domain acquisition module, for being defined according to the benchmark recording curve and the recording curve acquisition to be synchronized
Domain；
Continuous coupling function acquisition module, for obtaining continuous coupling function according to the domain and least square method；
Optimum Matching timeobtaining module, it is described continuous for being obtained according to the simulated annealing traversed in advance based on the overall situation
The Optimum Matching time of adaptation function；
Discrete adaptation function acquisition module, it is discrete for being carried out to the continuous coupling function according to discrete integration formula
Change, obtains discrete adaptation function；
Synchronization module, for realizing that the benchmark recording is bent according to the Optimum Matching time and the discrete adaptation function
Line is synchronous with the recording curve to be synchronized.
The third aspect, the embodiment of the present invention additionally provide a kind of distributed wave recording synchronizer, including processor, memory
And it is stored in the memory and is configured as the computer program executed by the processor, the processor executes institute
The distributed wave recording synchronous method as described in any of the abovedescribed one is realized when stating computer program.
Fourth aspect, the embodiment of the present invention additionally provide a kind of computer readable storage medium, described computerreadable to deposit
Storage media includes the computer program of storage, wherein the computerreadable storage is controlled when the computer program is run
Equipment where medium executes the distributed wave recording synchronous method described in abovementioned any one.
One technical solution of abovementioned technical proposal has the following advantages that：It obtains benchmark recording curve and recording to be synchronized is bent
Line；According to the benchmark recording curve and the recording curve acquisition domain to be synchronized；According to the domain and minimum two
Multiplication obtains continuous coupling function；The continuous coupling function is obtained according to the simulated annealing traversed in advance based on the overall situation most
Excellent match time reduces calculation amount under the premise of ensureing precision, accelerates the speed of search optimal value, can fast implement
The Optimum Matching time obtains, and carries out discretization to the continuous coupling function according to discrete integration formula, obtains discrete matching letter
Number；The benchmark recording curve and the recording to be synchronized are realized according to the Optimum Matching time and the discrete adaptation function
Curve synchronizes.With minimum operand and space hold, the calculating of discrete adaptation function is realized, when finally realizing different
Long, sample frequency curve synchronization can accurately realize that recording curve synchronizes.
Description of the drawings
Fig. 1 is the distributed wave recording synchronous method flow diagram that first embodiment of the invention provides；
Fig. 2 is the typical case schematic diagram that first embodiment of the invention provides effective domain；
Fig. 3 is the flow diagram for the acquisition continuous coupling functional based method that first embodiment of the invention provides；
Fig. 4 is the upper curve head that first embodiment of the invention provides and lower curve tail subfraction intersection schematic diagram；
Fig. 5 is the adaptation function curve synoptic diagram that first embodiment of the invention provides；
Fig. 6 is the adaptation function curve magnification schematic diagram that first embodiment of the invention provides；
Fig. 7 be first embodiment of the invention provide synchronize before recording curve comparison schematic diagram；
Fig. 8 be first embodiment of the invention provide synchronize rear recording curve comparison schematic diagram；
Fig. 9 is another distributed wave recording synchronous method flow diagram that second embodiment of the invention provides；
Figure 10 is another distributed wave recording synchronous method flow diagram that third embodiment of the invention provides；
Figure 11 is the sampling algorithm schematic diagram that third embodiment of the invention provides；
Figure 12 is a kind of distributed wave recording synchronizing device structural schematic diagram that fourth embodiment of the invention provides；
Figure 13 is the schematic diagram for the distributed wave recording synchronizer that fifth embodiment of the invention provides.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other without creative efforts
Embodiment shall fall within the protection scope of the present invention.
An embodiment of the present invention provides a kind of distributed wave recording synchronous method, for different durations, sample frequency curve
It is synchronous, it is described in detail separately below.
Embodiment one
The distributed wave recording synchronous method flow diagram provided referring to Fig. 1, first embodiment of the invention；
S11, benchmark recording curve and recording curve to be synchronized are obtained；
It should be noted that distributed wave recording includes the multiple recorded wave files for protecting equipment from multiple controls, each control, which is protected, to be set
Standby correspondence has corresponding recorded wave file, and each recording curve beginning and ending time having the same inside each recorded wave file adopts
Sample rate, that is to say, that each each recording of file internal is synchronous, but each recording for being different its inside of recorded wave file is bent
Line not necessarily has identical beginning and ending time, sample rate, to carry out the synchronization of multiple and different recorded wave files, as long as with wherein
Recorded wave file on the basis of one recorded wave file, certain recording curve a in the benchmark recorded wave file as benchmark recording curve,
A recording curve b is respectively chosen from other recorded wave files to be synchronized_{1}, b_{2}, b_{3}..., b_{n}It is synchronous with recording curve a respectively, you can
Realize the synchronization of all recorded wave files.
S12, according to the benchmark recording curve and the recording curve acquisition domain to be synchronized；
Preferably, described to include according to the benchmark recording curve and the recording curve acquisition domain to be synchronized：
Obtaining the benchmark recording curve and the recording curve to be synchronized has the condition of intersection；Wherein, described to have intersection
Condition beWherein, t_{b0}It is the left end point of the recording curve b to be synchronized, t_{bm}It is the record to be synchronized
The right endpoint of wave profile b, t_{aN}It is the right endpoint of the benchmark recording curve a, t_{a0}It is the left end point of the benchmark recording curve a,
τ is the Optimum Matching time；
Domain is obtained according to the condition for having intersection；Wherein, the domain includes Optimum Matching timedomain
With effective domain；Wherein, the Optimum Matching timedomain is：
τ∈(t_{a0}t_{bM},t_{aN}t_{b0})=(τ_{min},τ_{max})=T_{τ}
Effective domain is：t∈min(T_{a},T_{b})=(t_{min},t_{max})=T_{m}；Wherein, T_{a}For benchmark recording song
The domain of line a, T_{b}For the domain of the benchmark recording curve b.
Specifically, in the present embodiment, for two recording curves that needs synchronize, use a respectively, b is indicated, a
For the benchmark recording curve, b is the recording curve to be synchronized, its general domain is not consistent, if its domain is distinguished
For T_{a}And T_{b}：A=a (t), t ∈ T_{a}, b=b (t), t ∈ T_{b}；Construct b'(t), indicate b (t) records that offset THS is formed on a timeline
Wave profile：B'(t)=b (t+ τ), t ∈ T_{b}；Only in the case where curve a and curve b domains have intersection, adaptation function just has
Meaning, so it is m (τ)=m (a, b'), τ ∈ T to define adaptation function_{τ}；Curve a and curve b domains have 4 kinds of typical cases of intersection
Situation is as follows, and as shown in Fig. 2, the right endpoint of (1) curve b is overlapped with the left end point of curve a；(2) right endpoint of curve b and song
The right endpoint of line a overlaps；(3) left end point of curve b is overlapped with the left end point of curve a；(4) left end point of curve b and curve a
Right endpoint overlap；It follows that curve a and curve b domains have the condition of intersection to be exactly
It is arranged according to the condition for having intersection：
τ∈(t_{a0}t_{bM},t_{aN}t_{b0})=(τ_{min},τ_{max})=T_{τ}
So the corresponding effective domains of curve a and b are：
t∈min(T_{a},T_{b})=(t_{min},t_{max})=T_{m}。
S13, continuous coupling function is obtained according to the domain and least square method；
It should be noted that judging that the matching degree of two recording curves often uses least square method in the prior art：Meet τ ∈ T_{f}When the corresponding τ of m acquirements minimum value is Optimum Matching in range
Between, least square method refers to square integrating to the difference of two curves, finally integrates bigger, indicates two curves more not
Unanimously.Obviously, when two curve coinsides, it is 0 to integrate.And difference is bigger and time of coincidence length is bigger, integral is bigger.
In most cases least square method can meet the requirement that the corresponding τ of m acquirements minimum value is the Optimum Matching time.And in reality
In, the inventors discovered that it sometimes appear that the almost the same situation of a recording curve head and the tail waveform, such as in circuit therefore
Barrier restarts busbar voltage and line current in situation, and waveform, only meeting during failure are all up before and after failure
Significant change, if at this time only with least square method can misjudgment, obtain be before and after failure the Optimum Matching time mistake
Conclusion exists and judges mistaken ideas.
Preferably, described to include according to the domain and least square method acquisition continuous coupling function referring to Fig. 3：
S21, according to the time of coincidence length of the benchmark recording curve and the recording curve to be synchronized formed the time it is anti
To weights；
S22, the least square method with antitime weight is obtained according to the time reversal weights and least square method；
S23, continuous coupling function is obtained according to domain and the least square method of the antitime weight of the bandWherein, τ is the Optimum Matching time, and a (t) is the benchmark recording curve, b
(t) it is the recording curve to be synchronized, t_{max}For the benchmark recording curve and the effective domain of recording curve to be synchronized
Maximum value, t_{min}The minimum value of the benchmark recording curve and the effective domain of recording curve to be synchronized；τ∈T_{τ}, t ∈ T_{m}。
In the present embodiment, least square method difference is bigger and time of coincidence length is bigger, and integral is bigger.When coincidence
Between length be variation, it is possible to cause least square method integral mainly by time of coincidence effect length without by differentia influence.
For example, with reference to Fig. 4, possible least square method judges that integral is minimum when upper curve head and lower curve tail subfraction intersect, from
And assert difference minimum, this is unreasonable, so to construct an amount related with time of coincidence length, i.e., divided by t_{max}
t_{min}, formed time reversal weights, it is ensured that when comparing weight mainly the time of coincidence length zone of reasonableness in.Dead zone is constructed,
T ∈ T are thought roughly at present_{m}=(t_{min},t_{max}) be all reasonable, this section can be actually reduced, to exclude the apparent small model of head and the tail
Enclose the case where connecting.Actually recording, which is substantially, is triggered in its 10%~30% range, that is to say, that finally right
Neat point should be located in every 10~30% range of recording, if recording all meets this requirement, can specify that each
Curve is aligned in 5%~35% range and is reasonable, unreasonable more than the value of its range.
It should be noted that changing time reversal weights number, (t_{max}t_{min})^{h}, h is the real number more than 0, can basis
Different occasions optimize setting, and the present invention is not especially limited this.
When S14, basis obtain the Optimum Matching of the continuous coupling function based on the simulated annealing that the overall situation traverses in advance
Between；
In the present embodiment, first with traversal, the traversal curve of gamut is obtained, in the traversal curve smaller value
In range, then into the faster simulated annealing of trip temperature T decaying, to improve simulated annealing, in guarantee precision
Under the premise of, calculation amount is reduced, the speed of search optimal value is accelerated.
S15, discretization is carried out to the continuous coupling function according to discrete integration formula, obtains discrete adaptation function；
It should be noted that practical recording curve is discrete function, so needing the song of obtained continuous discrete function
Line has much into discretization, discrete integration formula, such as Rectangular Method, trapezoidal method, Simpson's integration method etc., and principle is similar, essence
Degree is different, and the present invention is not especially limited this.
Preferably, described that discretization is carried out to the continuous coupling function according to discrete integration formula, obtain discrete matching
Function includes：
Discretization is carried out to the continuous coupling function according to discrete integration formula, obtains the discrete adaptation functionWherein, g (t_{i})=a (t_{i})b(t_{i}+ τ), τ ∈ T_{τ},t∈T_{m}。
Preferably, the discrete integration formula is trapezoidal integration formula；
The trapezoidal integration formula isWherein, f (t_{i}) it is on trapezoidal
Bottom；f(t_{i+1}) it is trapezoidal bottom；(t_{i+1}t_{i}) it is trapezoidal height.
S16, the benchmark recording curve and described is realized according to the Optimum Matching time and the discrete adaptation function
Recording curve to be synchronized synchronizes.
Specifically, the pole bus current IDLH of the S1P1PPR obtained using in NW YUNNAN engineering experiment is bent as benchmark recording
The pole bus current IDLH of line a, S1P1CCP1 are as recording curve b to be synchronized, wherein the recording curve b to be synchronized compares institute
State benchmark recording curve a lag 1s, referring to Fig. 5 and Fig. 6, the high form in intermediate low both sides is presented in adaptation function m (τ), τ=
M is minimum when 1s, and τ is Optimum Matching point at this time, and referring to Fig. 7 and Fig. 8, curve difference 1s, synchronizes rear curve almost before synchronizing
It overlaps.
Implement the present embodiment to have the advantages that：Obtain benchmark recording curve and recording curve to be synchronized；According to institute
State benchmark recording curve and the recording curve acquisition domain to be synchronized；It is obtained and is connected according to the domain and least square method
Continuous adaptation function；When obtaining the Optimum Matching of the continuous coupling function according to the simulated annealing traversed in advance based on the overall situation
Between, under the premise of ensureing precision, calculation amount is reduced, the speed of search optimal value is accelerated, Optimum Matching can be fast implemented
Time obtains, and carries out discretization to the continuous coupling function according to discrete integration formula, obtains discrete adaptation function；According to institute
It states the Optimum Matching time and realizes that the benchmark recording curve is synchronous with the recording curve to be synchronized with the discrete adaptation function.
With minimum operand and space hold, the calculating of discrete adaptation function is realized, finally realizes different durations, sample frequency
The synchronization of curve.
Embodiment two is another distribution that second embodiment of the invention provides referring to Fig. 9 on the basis of embodiment one
Formula recording synchronous method flow diagram；
Preferably, the simulated annealing that the basis is traversed in advance based on the overall situation obtains the optimal of the continuous coupling function
Match time includes:
S31, the gamut lowdensity traversal curve that the continuous coupling function is obtained according to the traversal of small decile n；
S32, when the continuous coupling functional value is smaller, obtain the range of the corresponding lowdensity traversal curve；
S33, in the range, the continuous coupling function is obtained according to the temperature T faster simulated annealing that decays
The Optimum Matching time.
It should be noted that in adaptation function domain range T_{τ}It is optimal that m, which obtains the corresponding τ of minimum value, in range
With the time, common search optimization algorithm is as follows：
(1) traversal
By the domain T of adaptation function domain, that is, Optimum Matching time_{τ}It is divided into n minizone, all interval endpoints
For：τ_{i}=(i1) (τ_{max}τ_{min})/n, i=1,2,3 ... n.τ=τ is enabled successively_{i}, i=1,2,3 ..., n+1, and ask corresponding
With function m (τ), to obtain Optimum Matching time τ.The advantages of traversal is can to traverse entire effective range, fitting precision
Depending on the size of decile n, n is bigger, and decile is finer, and matching is more accurate, but the disadvantage is that as precision increases, calculation amount
It can significantly increase.
(2) hill climbing method
Set the initial value τ of τ_{0}With inceptive direction s=1, increase increment s Δs τ to τ every time, wherein s indicates direction, only
1 indicates τ augment directions with 1 two value, 1, and 1 indicates that τ reduces direction.Δ τ indicates augmental interval, is a positive number.sΔτ
Indicate the increment of τ.Calculate corresponding m (τ_{0}+ s Δ τ), if m (τ_{0}+ s Δ τ) increment Delta m<0, then it keeps s constant and continues growing
sΔτ；If m (τ_{0}+ s Δ τ) increment Delta m>0, then make s involutions ( 1) and continue growing s Δ τ, repeatedly, that is, first, s is not
Become, second is that s involutions ( 1), refer to that the deterministic process of this two kinds of situations carries out repeatedly repeatedly.Finally make m (τ_{0}+sΔ
τ) reach minimum value and mean that all amounts being related to finally think to find minimum value by several wheels cycle, algorithm terminates.
The advantages of hill climbing method, makes that for arbitrary recording curve its extreme value can be found, the disadvantage is that it is easily trapped into local optimum trap,
And global optimum can not be searched.
In embodiments of the present invention, it in order to improve hill climbing method, avoids being absorbed in local optimum trap, simulated annealing calculation can be used
Method.Simulated annealing has the ability for jumping out local optimum trap, i.e., random quantity is added on the often step increment of hill climbing method, steps on
Mountain has probability to carry out " direction of mistake " according to random quantity, even if falling into local optimum trap, after a period of time, algorithm
Trap can also be jumped out again, will finally be restrained toward the direction of global optimum.Simulated annealing is described in detail as follows：
Set the initial value τ of τ_{0}With inceptive direction s=1, temperature T ∈ [T_{min},T_{max}] it is initial temperature T_{0}.Increase every time to τ
Increment s Δ τ, calculate corresponding m (τ_{0}+ s Δ τ), if m (τ_{0}+ s Δ τ) increment Delta m<0, then it keeps s constant and continues growing s Δs
τ；If m (τ_{0}+ s Δ τ) increment Delta m>0, then it keeps s constant according to probability P and continues growing s Δ τ, otherwise make s involutions ( 1)
And s Δ τ are continued growing, repeatedly, finally make m (τ_{0}+ s Δ τ) reach minimum value.The expression formula of probability P is as follows：Wherein, Δ m is the increment of adaptation function；T_{max}For maximum temperature；T_{i}For Current Temperatures, subtract as i increases
It is small；I is current calculating number of steps；Simulated annealing receives the solution of the current solution difference of a ratio with certain probability, therefore
It is possible that the optimal solution of this part can be jumped out, reach global optimal solution.Simulated annealing is searching locally optimal solution
Afterwards, movements of the m toward higher value can be received with certain probability, perhaps after as not several times being the movement of local optimum
Globe optimum can be reached, local optimum has then just been jumped out.Ti becomes smaller with i in formula, specific formula be T (i+1)=
R*Ti, r are a fixed real number and 0<r<1, r is bigger, and search is slower, but is more easily found optimal value.
It should be noted that the calculation amount of simulated annealing is related with precision, when the T decaying of simulated annealing is very fast
When, still precision is relatively low for speed, and in turn when T decaying is slower, precision is higher but speed is slower.
In the present embodiment, first with traversal smaller decile n, the lowdensity traversal curve of gamut is obtained, herein
In the range of lowdensity traverses curve m smaller values, then the faster simulated annealing of T decaying is carried out, is moved back to improve simulation
Fiery algorithm reduces calculation amount under the premise of ensureing precision, accelerates the speed of search optimal value.
Implement the present embodiment to have the advantages that：Under the premise of ensureing precision, calculation amount is reduced, accelerates and searches
The speed of rope optimal value is absorbed in local optimum trap convenient for more rapidly finding the Optimum Matching time, and being avoided that, can not search
To global optimum.
Embodiment three is the another kind point that third embodiment of the invention provides referring to Figure 10 on the basis of embodiment one
Cloth recording synchronous method flow diagram；
It is described according to the Optimum Matching time and the discrete adaptation function realizes the benchmark recording curve and described
Recording curve synchronization to be synchronized includes：
S41, the former sampled point for obtaining the benchmark recording curve and the recording curve to be synchronized；
S42, the discrete adaptation function is calculated respectively according to the former sampled point, obtains discrete after operation
With function；
S43, realize that the benchmark recording is bent according to the discrete adaptation function after the Optimum Matching time and the operation
Line is synchronous with the recording curve to be synchronized.
It should be noted that for two curves to be synchronized, recording duration, sample frequency are generally inconsistent, do not deposit
The case where sampled point corresponds, can not directly it calculate.Previous common method is constant frequency Lagrange's interpolation algorithm,
Resampling is carried out with two curves of constant sample rate pair and obtain new curve, then calculated.But there are one for this algorithm
A problem, in order to ensure that sampling can cover all the points of two curves, the sampling period can be less than or equal to two recording curves from
The minimum interval of scatterplot, sample frequency is generally higher, has larger impact for memory space and arithmetic speed；And if adopting
With relatively low sample frequency, then whole original recording points can not be covered, sample waveform is caused to be distorted.
In the present embodiment, in order to which while retaining whole raw data points, reduction sample frequency is to reduce in terms of subsequently
The complexity of calculation, set forth herein nonfixed frequency sampling curve blending algorithms, and two recording curve original sampled points are only covered in sampling
The position of appearance is described in detail based on the following Figure 11 for convenience of description：
The the 1st, 2 article of curve indicates that the distribution of recording curve a, b on a timeline, curve a expression formulas are as follows respectively in Figure 11
It is shown：
Curve b expression formulas are as follows：
Wherein, the duration of recording curve a, b, sample rate are all inconsistent, so can not directly be calculated, so herein will
The the 3rd, 4 article of curve indicates that recording is bent respectively in the consistent recording curve a` and b` in recording curve a, b transformation definition domain, Figure 11
Line a` and b`.The 5th article of curve indicates the distributions of recording y on a timeline after a` and b` operations in Figure 11：
In order to calculate the domain of recording curve a` and b` and y, first with recording curve a, the period that b is overlapped is as y
Effective period of time, i.e.,In effective period of time, the domain of recording curve a` and b` and y is a, b
The union of domain：
In domain T_{y}It is interior, a'(t), b'(t) and it can be expressed as：With
Wherein, f_{Inter}(a, t) indicates to carry out a the interpolation of moment t, f_{Inter}(b, t) indicates to carry out b the interpolation of moment t, and interpolation can be used
A variety of interpolation algorithms such as Lagrangian linear interpolation method, cubic spline interpolation, in order to describe for the sake of simplicity, using Lagrangian line
Property interpolation method.By taking curve a as an example, if the section where t is t ∈ [t_{ai},t_{a(i+1)}], then：
Domain T_{y}Interior, y curves are represented by
Implement the present embodiment to have the advantages that：With minimum operand and space hold, realize different durations,
The synchronization of sample frequency curve finally realizes the calculating of discrete adaptation function.
It is a kind of distributed wave recording synchronizing device that fourth embodiment of the invention provides referring to Figure 12, Figure 12, including：
Curve acquisition module 11, for obtaining benchmark recording curve and recording curve to be synchronized；
Domain acquisition module 12, for being defined according to the benchmark recording curve and the recording curve acquisition to be synchronized
Domain；
Continuous coupling function acquisition module 13, for obtaining continuous coupling letter according to the domain and least square method
Number；
Optimum Matching timeobtaining module 14, for obtaining the company according to the simulated annealing traversed in advance based on the overall situation
The Optimum Matching time of continuous adaptation function；
Discrete adaptation function acquisition module 15, it is discrete for being carried out to the continuous coupling function according to discrete integration formula
Change, obtains discrete adaptation function；
Synchronization module 16, for realizing the benchmark recording according to the Optimum Matching time and the discrete adaptation function
Curve is synchronous with the recording curve to be synchronized.
Preferably, the domain acquisition module 12 includes：Obtain the benchmark recording curve and the recording to be synchronized
Curve has the condition of intersection；Wherein, described to there is the condition of intersection to beWherein, t_{b0}It is the record to be synchronized
The left end point of wave profile b, t_{bm}It is the right endpoint of the recording curve b to be synchronized, t_{aN}It is the right end of the benchmark recording curve a
Point, t_{a0}It is the left end point of the benchmark recording curve a, τ is the Optimum Matching time；
Domain is obtained according to the condition for having intersection；Wherein, the domain includes Optimum Matching timedomain
With effective domain；Wherein, the Optimum Matching timedomain is：
τ∈(t_{a0}t_{bM},t_{aN}t_{b0})=(τ_{min},τ_{max})=T_{τ}
Effective domain is：t∈min(T_{a},T_{b})=(t_{min},t_{max})=T_{m}；Wherein, T_{a}For benchmark recording song
The domain of line a, T_{b}For the domain of the benchmark recording curve b.
Preferably, the continuous coupling function acquisition module 13 includes：
Reversed weights acquiring unit, when being used for the coincidence according to the benchmark recording curve and the recording curve to be synchronized
Between length formed time reversal weights；
Square law acquiring unit, for being obtained with antitime weight according to the time reversal weights and least square method
Least square method；
Continuous coupling function acquiring unit, for being obtained according to domain and the least square method of the antitime weight of the band
Continuous coupling functionWherein, τ is the Optimum Matching time, and a (t) is the base
Quasi recording curve, b (t) are the recording curve to be synchronized, t_{max}For the benchmark recording curve and the recording curve to be synchronized
The maximum value of effective domain, t_{min}The minimum value of the benchmark recording curve and the effective domain of recording curve to be synchronized；
τ∈T_{τ}, t ∈ T_{m}。
Preferably, the Optimum Matching timeobtaining module 14 includes：
Curve acquisition unit, the gamut for obtaining the continuous coupling function according to the traversal of small decile n are low close
Degree traversal curve；
Range acquiring unit, for when the continuous coupling functional value is smaller, obtaining the corresponding lowdensity traversal
The range of curve；
Time acquisition unit, in the range, being obtained according to the faster simulated annealing of temperature T decaying described
The Optimum Matching time of continuous coupling function.
Preferably, the discrete adaptation function acquisition module 15 includes：
Discrete adaptation function acquiring unit, it is discrete for being carried out to the continuous coupling function according to discrete integration formula
Change, obtains the discrete adaptation functionWherein, g (t_{i})=a (t_{i})b
(t_{i}+ τ), τ ∈ T_{τ},t∈T_{m}。
Preferably, the discrete adaptation function acquiring unit includes：
The discrete integration formula is trapezoidal integration formula；
The trapezoidal integration formula isWherein, f (t_{i}) it is on trapezoidal
Bottom；f(t_{i+1}) it is trapezoidal bottom；(t_{i+1}t_{i}) it is trapezoidal height.
Preferably, the synchronization module 16 includes：
Former sampled point acquiring unit, the former sampling for obtaining the benchmark recording curve and the recording curve to be synchronized
Point；
Operation function acquiring unit, for being calculated respectively the discrete adaptation function according to the former sampled point,
Obtain the discrete adaptation function after operation；
Synchronization unit, for realizing the base according to the discrete adaptation function after the Optimum Matching time and the operation
Quasi recording curve is synchronous with the recording curve to be synchronized.
Implement the present embodiment to have the advantages that：
Obtain benchmark recording curve and recording curve to be synchronized；According to the benchmark recording curve and the recording to be synchronized
Curve acquisition domain；Continuous coupling function is obtained according to the domain and least square method；According to based on global pre traversal
Simulated annealing obtain Optimum Matching time of the continuous coupling function and reduce meter under the premise of ensureing precision
Calculation amount accelerates the speed of search optimal value, the acquisition of Optimum Matching time can be fast implemented, according to discrete integration formula to described
Continuous coupling function carries out discretization, obtains discrete adaptation function；According to the Optimum Matching time and the discrete matching letter
Number realizes that the benchmark recording curve is synchronous with the recording curve to be synchronized.With minimum operand and space hold, realize
The calculating of discrete adaptation function, finally realizes the synchronization of different durations, sample frequency curve.
Figure 13 is referred to, Figure 13 is the schematic diagram for the distributed wave recording synchronizer that fifth embodiment of the invention provides, and is used
In executing distributed wave recording synchronous method provided in an embodiment of the present invention, as shown in figure 13, the terminal which synchronizes
Equipment includes：At least one processor 11, such as CPU, at least one network interface 14 or other users interface 13, memory
15, at least one communication bus 12, communication bus 12 is for realizing the connection communication between these components.Wherein, user interface
13 may include optionally USB interface and other standards interface, wireline interface.Network interface 14 may include optionally Wi
Fi interfaces and other wireless interfaces.Memory 15 may include highspeed RAM memory, it is also possible to further include nonlabile deposit
Reservoir (nonvolatilememory), for example, at least a magnetic disk storage.Memory 15 can include optionally at least one
A storage device for being located remotely from aforementioned processor 11.
In some embodiments, memory 15 stores following element, executable modules or data structures, or
Their subset or their superset:
Operating system 151, including various system programs, for realizing various basic businesses and hardware based of processing
Business；
Program 152.
Specifically, processor 11 executes point described in abovedescribed embodiment for calling the program 152 stored in memory 15
Cloth recording synchronous method.
Alleged processor can be central processing unit (Central Processing Unit, CPU), can also be it
His general processor, digital signal processor (Digital Signal Processor, DSP), applicationspecific integrated circuit
(Application Specific Integrated Circuit, ASIC), readymade programmable gate array (Field
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor can also be any conventional processor
Deng the processor is the control centre of the distributed wave recording synchronous method, utilizes the entire institute of various interfaces and connection
State the various pieces of distributed wave recording synchronous method.
The memory can be used for storing the computer program and/or module, and the processor is by running or executing
Computer program in the memory and/or module are stored, and calls the data being stored in memory, realizes distribution
The various functions for the electronic device that formula recording synchronizes.The memory can include mainly storing program area and storage data field,
In, storing program area can storage program area, application program (such as soundplaying function, word needed at least one function
Conversion function etc.) etc.；Storage data field can store that (for example audio data, word disappear according to the created data that use of mobile phone
Cease data etc.) etc..In addition, memory may include highspeed random access memory, can also include nonvolatile memory, example
Such as hard disk, memory, plugin type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital (Secure
Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or other volatibility are solid
State memory device.
Wherein, if the distributed wave recording synchronizes integrated module and is realized in the form of SFU software functional unit and as only
Vertical product is sold or in use, can be stored in a computer read/write memory medium.Based on this understanding, this hair
All or part of flow in bright realization abovedescribed embodiment method, can also be instructed by computer program relevant hardware come
It completes, the computer program can be stored in a computer readable storage medium, and the computer program is by processor
When execution, it can be achieved that the step of abovementioned each embodiment of the method.Wherein, the computer program includes computer program code,
The computer program code can be source code form, object identification code form, executable file or certain intermediate forms etc..Institute
Stating computerreadable medium may include：Can carry the computer program code any entity or device, recording medium,
USB flash disk, mobile hard disk, magnetic disc, CD, computer storage, readonly memory (ROM, ReadOnly Memory), arbitrary access
Memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and software distribution medium etc..It needs
It is bright, the content that the computerreadable medium includes can according in jurisdiction legislation and patent practice requirement into
Row increase and decrease appropriate, such as in certain jurisdictions, according to legislation and patent practice, computerreadable medium does not include that electricity carries
Wave signal and telecommunication signal.
It should be noted that the apparatus embodiments described above are merely exemplary, wherein described be used as separating component
The unit of explanation may or may not be physically separated, and the component shown as unit can be or can also
It is not physical unit, you can be located at a place, or may be distributed over multiple network units.It can be according to actual
It needs that some or all of module therein is selected to achieve the purpose of the solution of this embodiment.In addition, device provided by the invention
In embodiment attached drawing, the connection relation between module indicates there is communication connection between them, specifically can be implemented as one or
A plurality of communication bus or signal wire.Those of ordinary skill in the art are without creative efforts, you can to understand
And implement.
The above is the preferred embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
It should be noted that in the abovedescribed embodiments, all emphasizing particularly on different fields to the description of each embodiment, in some embodiment
In the part that is not described in, may refer to the associated description of other embodiment.Secondly, those skilled in the art should also know
It knows, embodiment described in this description belongs to preferred embodiment, and involved action and simulation must be that the present invention must
Must.
Claims (10)
1. a kind of distributed wave recording synchronous method, which is characterized in that including：
Obtain benchmark recording curve and recording curve to be synchronized；
According to the benchmark recording curve and the recording curve acquisition domain to be synchronized；
Continuous coupling function is obtained according to the domain and least square method；
The Optimum Matching time of the continuous coupling function is obtained according to the simulated annealing traversed in advance based on the overall situation；
Discretization is carried out to the continuous coupling function according to discrete integration formula, obtains discrete adaptation function；
The benchmark recording curve and the record to be synchronized are realized according to the Optimum Matching time and the discrete adaptation function
Curve line locking.
2. distributed wave recording synchronous method according to claim 1, which is characterized in that described according to benchmark recording song
Line and the recording curve acquisition domain to be synchronized include：
Obtaining the benchmark recording curve and the recording curve to be synchronized has the condition of intersection；Wherein, the item for having intersection
Part isWherein, t_{b0}It is the left end point of the recording curve b to be synchronized, t_{bm}It is that the recording to be synchronized is bent
The right endpoint of line b, t_{aN}It is the right endpoint of the benchmark recording curve a, t_{a0}It is the left end point of the benchmark recording curve a, τ is
The Optimum Matching time；
Domain is obtained according to the condition for having intersection；Wherein, the domain includes Optimum Matching timedomain and has
Imitate domain；Wherein, the Optimum Matching timedomain is：
τ∈(t_{a0}t_{bM},t_{aN}t_{b0})=(τ_{min},τ_{max})=T_{τ}
Effective domain is：t∈min(T_{a},T_{b})=(t_{min},t_{max})=T_{m}；Wherein, T_{a}For the benchmark recording curve a's
Domain, T_{b}For the domain of the benchmark recording curve b.
3. distributed wave recording synchronous method according to claim 2, which is characterized in that it is described according to the domain and most
Small square law obtains continuous coupling function：
Time reversal weights are formed according to the time of coincidence length of the benchmark recording curve and the recording curve to be synchronized；
The least square method with antitime weight is obtained according to the time reversal weights and least square method；
Continuous coupling function is obtained according to domain and the least square method of the antitime weight of the band
Wherein, τ is the Optimum Matching time, and a (t) is the benchmark recording curve, and b (t) is the recording curve to be synchronized, t_{max}For institute
State the maximum value of benchmark recording curve and the effective domain of recording curve to be synchronized, t_{min}The benchmark recording curve and institute
State the minimum value of the effective domain of recording curve to be synchronized；τ∈T_{τ}, t ∈ T_{m}。
4. distributed wave recording synchronous method according to claim 1, which is characterized in that the basis is based on global pre traversal
Simulated annealing obtain Optimum Matching time of the continuous coupling function and include:
The gamut lowdensity that the continuous coupling function is obtained according to the traversal of small decile n traverses curve；
When the continuous coupling functional value is smaller, the range of the corresponding lowdensity traversal curve is obtained；
In the range, optimal of the continuous coupling function is obtained according to the temperature T faster simulated annealings that decay
With the time.
5. distributed wave recording synchronous method according to claim 1, which is characterized in that described according to discrete integration formula pair
The continuous coupling function carries out discretization, and obtaining discrete adaptation function includes：
Discretization is carried out to the continuous coupling function according to discrete integration formula, obtains the discrete adaptation functionWherein, g (t_{i})=a (t_{i})b(t_{i}+ τ), τ ∈ T_{τ},t∈T_{m}。
6. distributed wave recording synchronous method according to claim 5, which is characterized in that the discrete integration formula is trapezoidal
Integral formula；
The trapezoidal integration formula isWherein, f (t_{i}) it is trapezoidal upper bottom；f
(t_{i+1}) it is trapezoidal bottom；(t_{i+1}t_{i}) it is trapezoidal height.
7. distributed wave recording synchronous method according to claim 1, which is characterized in that it is described according to the Optimum Matching when
Between with the discrete adaptation function realize the benchmark recording curve and the recording curve to be synchronized synchronization include：
Obtain the former sampled point of the benchmark recording curve and the recording curve to be synchronized；
The discrete adaptation function is calculated respectively according to the former sampled point, obtains the discrete adaptation function after operation；
The benchmark recording curve and described is realized according to the discrete adaptation function after the Optimum Matching time and the operation
Recording curve to be synchronized synchronizes.
8. a kind of distributed wave recording synchronizing device, which is characterized in that including：
Curve acquisition module, for obtaining benchmark recording curve and recording curve to be synchronized；
Domain acquisition module, for according to the benchmark recording curve and the recording curve acquisition domain to be synchronized；
Continuous coupling function acquisition module, for obtaining continuous coupling function according to the domain and least square method；
Optimum Matching timeobtaining module, for obtaining the continuous coupling according to the simulated annealing traversed in advance based on the overall situation
The Optimum Matching time of function；
Discrete adaptation function acquisition module is obtained for carrying out discretization to the continuous coupling function according to discrete integration formula
Take discrete adaptation function；
Synchronization module, for according to the Optimum Matching time and the discrete adaptation function realize the benchmark recording curve and
The recording curve to be synchronized synchronizes.
9. a kind of distributed wave recording synchronizer, including processor, memory and it is stored in the memory and is configured
For the computer program executed by the processor, the processor realizes such as claim 1 when executing the computer program
To the distributed wave recording synchronous method described in any one of 7.
10. a kind of computer readable storage medium, which is characterized in that the computer readable storage medium includes the calculating of storage
Machine program, wherein equipment where controlling the computer readable storage medium when the computer program is run is executed as weighed
Profit requires the distributed wave recording synchronous method described in any one of 1 to 7.
Priority Applications (1)
Application Number  Priority Date  Filing Date  Title 

CN201810473612.0A CN108680824A (en)  20180516  20180516  Distributed wave recording synchronous method, device, equipment and medium 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

CN201810473612.0A CN108680824A (en)  20180516  20180516  Distributed wave recording synchronous method, device, equipment and medium 
Publications (1)
Publication Number  Publication Date 

CN108680824A true CN108680824A (en)  20181019 
Family
ID=63805636
Family Applications (1)
Application Number  Title  Priority Date  Filing Date 

CN201810473612.0A Pending CN108680824A (en)  20180516  20180516  Distributed wave recording synchronous method, device, equipment and medium 
Country Status (1)
Country  Link 

CN (1)  CN108680824A (en) 
Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

JP2004132746A (en) *  20021008  20040430  Nippon Kouatsu Electric Co  Fault locating method and fault locating system 
CN101900777A (en) *  20090527  20101201  河南省电力勘测设计院  Monitoring method and device of power system 
CN102243272A (en) *  20110401  20111116  重庆大学  Highprecision harmonic analysis method of sampling data synchronization 
CN102288170A (en) *  20110714  20111221  浙江大学  Correction method of electronic compass in underwater vehicle 
CN105606894A (en) *  20160128  20160525  南京信息工程大学  Instantaneous frequency estimation method based on simulated annealing (SA) 
CN107959969A (en) *  20171211  20180424  武汉中原电子信息有限公司  A kind of method for synchronizing time applied to transient state recording type fault detector 

2018
 20180516 CN CN201810473612.0A patent/CN108680824A/en active Pending
Patent Citations (6)
Publication number  Priority date  Publication date  Assignee  Title 

JP2004132746A (en) *  20021008  20040430  Nippon Kouatsu Electric Co  Fault locating method and fault locating system 
CN101900777A (en) *  20090527  20101201  河南省电力勘测设计院  Monitoring method and device of power system 
CN102243272A (en) *  20110401  20111116  重庆大学  Highprecision harmonic analysis method of sampling data synchronization 
CN102288170A (en) *  20110714  20111221  浙江大学  Correction method of electronic compass in underwater vehicle 
CN105606894A (en) *  20160128  20160525  南京信息工程大学  Instantaneous frequency estimation method based on simulated annealing (SA) 
CN107959969A (en) *  20171211  20180424  武汉中原电子信息有限公司  A kind of method for synchronizing time applied to transient state recording type fault detector 
Similar Documents
Publication  Publication Date  Title 

US20170329462A1 (en)  Graphical user interface for static and adaptive thresholds  
CN106067080B (en)  Configurable workflow capabilities are provided  
US10333799B2 (en)  Monitoring IT services at an individual overall level from machine data  
US9864493B2 (en)  Cohortbased presentation of user interaction data  
US9256908B2 (en)  Utility consumption disaggregation using low sample rate smart meters  
Melamed  TES: A class of methods for generating autocorrelated uniform variates  
US7421460B2 (en)  Method for determining execution of backup on a database  
CN104252452B (en)  The method and device of data management  
CN107438836A (en)  Complex event processing device for the data of history/in real time/replay  
US9519692B2 (en)  Visualizing related events within a timeline  
US8041652B2 (en)  Measuring web site satisfaction of information needs using page traffic profile  
US10534643B2 (en)  Correlation of thread intensity and heap usage to identify heaphoarding stack traces  
US9285242B2 (en)  Determining a connectivity model in smart grids  
US9247502B2 (en)  Method and apparatus for power management on mobile devices and associated mobile devices  
CN103499806B (en)  For portable detector and the method for electric energy measuring equipment  
CN106096905A (en)  A kind of project management method and system  
Zhang et al.  Stochastic model predictive control using a combination of randomized and robust optimization  
CN105488130A (en)  Associated data visualization display method  
US10692048B2 (en)  Apparatus and method for generating a chronological timesheet  
US20160187864A1 (en)  Apparatus and System for Providing Energy Information  
US20130191052A1 (en)  Realtime simulation of power grid disruption  
CN1225713C (en)  Design method for specific chip of intelligent electric appliance  
Brunst et al.  Comprehensive performance tracking with vampir 7  
US8539367B2 (en)  Automatic repositioning of graphical program nodes during node placement or node movement  
CN105069690B (en)  A kind of grid operation data processing method 
Legal Events
Date  Code  Title  Description 

PB01  Publication  
PB01  Publication  
SE01  Entry into force of request for substantive examination  
SE01  Entry into force of request for substantive examination 