CN108680824A - Distributed wave recording synchronous method, device, equipment and medium - Google Patents

Distributed wave recording synchronous method, device, equipment and medium Download PDF

Info

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
Application number
CN201810473612.0A
Other languages
Chinese (zh)
Inventor
卢远宏
郭琦
李书勇
郭海平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CSG Electric Power Research Institute
China Southern Power Grid Co Ltd
Power Grid Technology Research Center of China Southern Power Grid Co Ltd
Research Institute of Southern Power Grid Co Ltd
Original Assignee
China Southern Power Grid Co Ltd
Power Grid Technology Research Center of China Southern Power Grid Co Ltd
Research Institute of Southern Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Southern Power Grid Co Ltd, Power Grid Technology Research Center of China Southern Power Grid Co Ltd, Research Institute of Southern Power Grid Co Ltd filed Critical China Southern Power Grid Co Ltd
Priority to CN201810473612.0A priority Critical patent/CN108680824A/en
Publication of CN108680824A publication Critical patent/CN108680824A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/081Locating faults in cables, transmission lines, or networks according to type of conductors
    • G01R31/086Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/08Locating faults in cables, transmission lines, or networks
    • G01R31/088Aspects of digital computing
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/50Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06KRECOGNITION OF DATA; PRESENTATION OF DATA; RECORD CARRIERS; HANDLING RECORD CARRIERS
    • G06K9/00Methods or arrangements for reading or recognising printed or written characters or for recognising patterns, e.g. fingerprints
    • G06K9/00496Recognising patterns in signals and combinations thereof
    • G06K9/00536Classification; 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

Distributed wave recording synchronous method, device, equipment and medium
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 wave-recording for being distributed in each device, Each device recording has different a time span and sample rate, operations staff need manual sorting analysis failure wave-recording 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 above-mentioned 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, tb0It is the left end point of the recording curve b to be synchronized, tbmIt is the record to be synchronized The right endpoint of wave profile b, taNIt is the right endpoint of the benchmark recording curve a, ta0It 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 time-domain With effective domain;Wherein, the Optimum Matching time-domain is:
τ∈(ta0-tbM,taN-tb0)=(τminmax)=Tτ
Effective domain is:t∈min(Ta,Tb)=(tmin,tmax)=Tm;Wherein, TaFor benchmark recording song The domain of line a, TbFor 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 anti-time 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 anti-time 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, tmaxFor the maximum value of the benchmark recording curve and the effective domain of recording curve to be synchronized, tminThe benchmark recording is bent The minimum value of line and the effective domain of recording curve to be synchronized;τ∈Tτ, t ∈ Tm
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 low-density 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 low-density 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 (ti)=a (ti)-b(ti+ τ), τ ∈ Tτ,t∈Tm
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 (ti) it is on trapezoidal Bottom;f(ti+1) it is trapezoidal bottom;(ti+1-ti) 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 time-obtaining 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 above-described 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 computer-readable to deposit Storage media includes the computer program of storage, wherein the computer-readable storage is controlled when the computer program is run Equipment where medium executes the distributed wave recording synchronous method described in above-mentioned any one.
One technical solution of above-mentioned 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 sub-fraction 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 synchronized1, b2, b3..., bnIt 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, tb0It is the left end point of the recording curve b to be synchronized, tbmIt is the record to be synchronized The right endpoint of wave profile b, taNIt is the right endpoint of the benchmark recording curve a, ta0It 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 time-domain With effective domain;Wherein, the Optimum Matching time-domain is:
τ∈(ta0-tbM,taN-tb0)=(τminmax)=Tτ
Effective domain is:t∈min(Ta,Tb)=(tmin,tmax)=Tm;Wherein, TaFor benchmark recording song The domain of line a, TbFor 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 TaAnd Tb:A=a (t), t ∈ Ta, b=b (t), t ∈ Tb;Construct b'(t), indicate b (t) records that offset THS is formed on a timeline Wave profile:B'(t)=b (t+ τ), t ∈ Tb;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:
τ∈(ta0-tbM,taN-tb0)=(τminmax)=Tτ
So the corresponding effective domains of curve a and b are:
t∈min(Ta,Tb)=(tmin,tmax)=Tm
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 τ ∈ TfWhen 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 co-insides, 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 anti-time 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 anti-time 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, tmaxFor the benchmark recording curve and the effective domain of recording curve to be synchronized Maximum value, tminThe minimum value of the benchmark recording curve and the effective domain of recording curve to be synchronized;τ∈Tτ, t ∈ Tm
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 sub-fraction intersect, from And assert difference minimum, this is unreasonable, so to construct an amount related with time of coincidence length, i.e., divided by tmax- tmin, 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 presentm=(tmin,tmax) 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, (tmax-tmin)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 (ti)=a (ti)-b(ti+ τ), τ ∈ Tτ,t∈Tm
Preferably, the discrete integration formula is trapezoidal integration formula;
The trapezoidal integration formula isWherein, f (ti) it is on trapezoidal Bottom;f(ti+1) it is trapezoidal bottom;(ti+1-ti) 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 low-density 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 low-density 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=(i-1) (τmaxmin)/n, i=1,2,3 ... n.τ=τ is enabled successivelyi, 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 τ0With 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 τ0With inceptive direction s=1, temperature T ∈ [Tmin,Tmax] it is initial temperature T0.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;TmaxFor maximum temperature;TiFor 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 low-density traversal curve of gamut is obtained, herein In the range of low-density 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 non-fixed 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 TyIt is interior, a'(t), b'(t) and it can be expressed as:With Wherein, fInter(a, t) indicates to carry out a the interpolation of moment t, fInter(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 ∈ [tai,ta(i+1)], then: Domain TyInterior, 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 time-obtaining 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, tb0It is the record to be synchronized The left end point of wave profile b, tbmIt is the right endpoint of the recording curve b to be synchronized, taNIt is the right end of the benchmark recording curve a Point, ta0It 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 time-domain With effective domain;Wherein, the Optimum Matching time-domain is:
τ∈(ta0-tbM,taN-tb0)=(τminmax)=Tτ
Effective domain is:t∈min(Ta,Tb)=(tmin,tmax)=Tm;Wherein, TaFor benchmark recording song The domain of line a, TbFor 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 anti-time 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 anti-time 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, tmaxFor the benchmark recording curve and the recording curve to be synchronized The maximum value of effective domain, tminThe minimum value of the benchmark recording curve and the effective domain of recording curve to be synchronized; τ∈Tτ, t ∈ Tm
Preferably, the Optimum Matching time-obtaining 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 low-density 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 (ti)=a (ti)-b (ti+ τ), τ ∈ Tτ,t∈Tm
Preferably, the discrete adaptation function acquiring unit includes:
The discrete integration formula is trapezoidal integration formula;
The trapezoidal integration formula isWherein, f (ti) it is on trapezoidal Bottom;f(ti+1) it is trapezoidal bottom;(ti+1-ti) 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 high-speed RAM memory, it is also possible to further include non-labile deposit Reservoir (non-volatilememory), 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 above-described 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), application-specific integrated circuit (Application Specific Integrated Circuit, ASIC), ready-made 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 sound-playing 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 high-speed random access memory, can also include nonvolatile memory, example Such as hard disk, memory, plug-in 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 above-described 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 above-mentioned 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 computer-readable 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, read-only memory (ROM, Read-Only 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 computer-readable 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, computer-readable 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 above-described 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, tb0It is the left end point of the recording curve b to be synchronized, tbmIt is that the recording to be synchronized is bent The right endpoint of line b, taNIt is the right endpoint of the benchmark recording curve a, ta0It 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 time-domain and has Imitate domain;Wherein, the Optimum Matching time-domain is:
τ∈(ta0-tbM,taN-tb0)=(τminmax)=Tτ
Effective domain is:t∈min(Ta,Tb)=(tmin,tmax)=Tm;Wherein, TaFor the benchmark recording curve a's Domain, TbFor 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 anti-time 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 anti-time 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, tmaxFor institute State the maximum value of benchmark recording curve and the effective domain of recording curve to be synchronized, tminThe benchmark recording curve and institute State the minimum value of the effective domain of recording curve to be synchronized;τ∈Tτ, t ∈ Tm
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 low-density 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 low-density 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 (ti)=a (ti)-b(ti+ τ), τ ∈ Tτ,t∈Tm
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 (ti) it is trapezoidal upper bottom;f (ti+1) it is trapezoidal bottom;(ti+1-ti) 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 time-obtaining 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.
CN201810473612.0A 2018-05-16 2018-05-16 Distributed wave recording synchronous method, device, equipment and medium Pending CN108680824A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810473612.0A CN108680824A (en) 2018-05-16 2018-05-16 Distributed wave recording synchronous method, device, equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810473612.0A CN108680824A (en) 2018-05-16 2018-05-16 Distributed wave recording synchronous method, device, equipment and medium

Publications (1)

Publication Number Publication Date
CN108680824A true CN108680824A (en) 2018-10-19

Family

ID=63805636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810473612.0A Pending CN108680824A (en) 2018-05-16 2018-05-16 Distributed wave recording synchronous method, device, equipment and medium

Country Status (1)

Country Link
CN (1) CN108680824A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004132746A (en) * 2002-10-08 2004-04-30 Nippon Kouatsu Electric Co Fault locating method and fault locating system
CN101900777A (en) * 2009-05-27 2010-12-01 河南省电力勘测设计院 Monitoring method and device of power system
CN102243272A (en) * 2011-04-01 2011-11-16 重庆大学 High-precision harmonic analysis method of sampling data synchronization
CN102288170A (en) * 2011-07-14 2011-12-21 浙江大学 Correction method of electronic compass in underwater vehicle
CN105606894A (en) * 2016-01-28 2016-05-25 南京信息工程大学 Instantaneous frequency estimation method based on simulated annealing (SA)
CN107959969A (en) * 2017-12-11 2018-04-24 武汉中原电子信息有限公司 A kind of method for synchronizing time applied to transient state recording type fault detector

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004132746A (en) * 2002-10-08 2004-04-30 Nippon Kouatsu Electric Co Fault locating method and fault locating system
CN101900777A (en) * 2009-05-27 2010-12-01 河南省电力勘测设计院 Monitoring method and device of power system
CN102243272A (en) * 2011-04-01 2011-11-16 重庆大学 High-precision harmonic analysis method of sampling data synchronization
CN102288170A (en) * 2011-07-14 2011-12-21 浙江大学 Correction method of electronic compass in underwater vehicle
CN105606894A (en) * 2016-01-28 2016-05-25 南京信息工程大学 Instantaneous frequency estimation method based on simulated annealing (SA)
CN107959969A (en) * 2017-12-11 2018-04-24 武汉中原电子信息有限公司 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) Cohort-based 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 heap-hoarding 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) Real-time 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 re-positioning 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