CN106788092B - A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition - Google Patents
A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition Download PDFInfo
- Publication number
- CN106788092B CN106788092B CN201710112598.7A CN201710112598A CN106788092B CN 106788092 B CN106788092 B CN 106788092B CN 201710112598 A CN201710112598 A CN 201710112598A CN 106788092 B CN106788092 B CN 106788092B
- Authority
- CN
- China
- Prior art keywords
- atom
- signal
- current
- parameter
- decomposition
- 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.)
- Active
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 45
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000013016 damping Methods 0.000 claims description 13
- 230000001052 transient effect Effects 0.000 claims description 12
- 238000004804 winding Methods 0.000 claims description 8
- 238000005457 optimization Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000005284 excitation Effects 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000012360 testing method Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 6
- 230000007812 deficiency Effects 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P23/00—Arrangements or methods for the control of AC motors characterised by a control method other than vector control
- H02P23/14—Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P9/00—Arrangements for controlling electric generators for the purpose of obtaining a desired output
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2103/00—Controlling arrangements characterised by the type of generator
- H02P2103/20—Controlling arrangements characterised by the type of generator of the synchronous type
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P2207/00—Indexing scheme relating to controlling arrangements characterised by the type of motor
- H02P2207/05—Synchronous machines, e.g. with permanent magnets or DC excitation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
The PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition that the invention discloses a kind of, characterized in that the following steps are included: 1) constructing relevant atomic library and by atom discretization continuous parameters;2) evolution matching pursuit algorithm is used, fundamental current, DC current, frequency multiplication electric current are quickly extracted from synchronous motor sudden short-circuit current;3) parameter of synchronous machine is extrapolated from each atom signals characteristic parameter of decomposition.Advantageous effects of the invention: being directed to the deficiency of traditional algorithm, in particular by empirical mode decomposition (EMD) and Prony algorithm it is not good enough to the discomposing effect of Noise signal shorts electric current, noiseproof feature is bad, this algorithm can accurately extract parameter of synchronous machine, and have preferable noiseproof feature.
Description
Technical field
The PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition that the present invention relates to a kind of belongs to signal processing technology neck
Domain.
Background technique
Accurate PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE has important meaning to Operation of Electric Systems and design control system is researched and analysed
Justice.Wherein, reflect the transient parameter of synchronous motor transient process and power equipment selection, stability of power system calculate and event
Barrier calculating etc. is closely related.
The parameter of synchronous machine measurement method that GB/T1029-2005 is recommended to use is three-phase suddenly-applied short circuit method, passes through short circuit
Envelope extracting cycle component and aperiodic component above and below electric current.This data processing method obtains the resultant error of transient parameter
It is larger.Short circuit current can be represented by exponential functions, and the parameter identification of Prony algorithm tailored index function, the parameter of extraction
Effect is preferable, but Prony algorithm there are problems that determining noise-sensitive and order.
The influence of limitation and noise to PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE precision for conventional method proposes many
Improved method, but also there are still some drawbacks simultaneously for these improved methods:
The method that wavelet transformation and Prony are combined reduces information collection noise using wavelet transformation, obtained parameter
Error is smaller.By the array signal processing method total least square-constant technology of rotating vector (total least square-
Estimation of signal parameters via rotational invariance technology, TLS-
ESPRIT it) is applied in PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE, signal is subjected to Subspace partition and total least square (total
Least square, TLS) dual denoising, anti-noise ability is improved.
PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE side based on Hilbert-Huang transform (Hilbert-Huang transform, HHT)
Method constitutes a kind of new time domain filter based on empirical mode decomposition (empirical mode decomposition, EMD)
Wave method obtains preferable identification result under strong noise background.But there are insoluble " endpoint effects for the EMD signal of HHT
Answer " problem.
The discrimination method of (local mean decomposition, LMD) is decomposed based on local mean value, identification result is made an uproar
Sound shadow sound is smaller, but the sliding average span selection of LMD also needs further to study.
Summary of the invention
To solve the deficiencies in the prior art, the purpose of the present invention is to provide a kind of synchronous motors based on atom decomposition
Parameter identification method completely shows each local feature of a non-stationary signal, can cancelling noise well in parsing
Information improves the precision of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE.
In order to achieve the above objectives, the present invention adopts the following technical scheme that:
A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, characterized in that include the following steps:
1) relevant atomic library is constructed, by atom discretization continuous parameters;
2) evolution matching pursuit algorithm is used, fundamental current, direct current are quickly extracted from synchronous motor sudden short-circuit current
Electric current and frequency multiplication electric current;
3) parameter of synchronous machine is extrapolated from each atom signals characteristic parameter of decomposition.
Further, the particular content of the step 1) are as follows:
11) Fundamental atom library is constructedIn formula:f1For fundamental frequency,For fundamental phase, ts1At the beginning of indicating fundamental wave, te1Indicate the knot of fundamental wave
The beam moment;T indicates time variable;
U (t) is unit jump function;kr1To make | | gr1(t) | |=1 coefficient;
12) using the fundamental wave sinusoidal signal exponentially to decay in short circuit current and decaying frequency multiplication sinusoidal signal as damping
Sinusoidal atom, the damped sine atom of construction are as follows:
In formula:f2For signal frequency,For phase, ρ2For attenuation parameter, ts2At the beginning of indicating disturbance
It carves, te2Indicate the finish time of disturbance, u (t) is unit jump function, kr2To make | | gr2(t) | |=1 coefficient;
13) attenuating dc component in current signal, structural damping direct current atom are extractedIn formula: r3=[ρ3,ts3,te3], ρ3For attenuation coefficient, ts3Expression is opened
Begin the time, te3Indicate the end time, u (t) is unit jump function, kr3To make | | gr3(t) | |=1 coefficient;
Further, the coefficient value range in Fundamental atom library function are as follows: 49HZ≤f1≤ 51HZ, Wherein, N indicates the sampled data length of signal i (t) to be analyzed, fsIndicate signal i to be analyzed
(t) sample frequency;
Coefficient value range in damped sine atom library function are as follows:By finger
The value range of the decaying fundamental frequency of the fundamental wave sinusoidal signal of number rule decaying is 49HZ≤f2≤ 51HZ, decaying frequency multiplication are sinusoidal
The value range of the component frequencies of signal is 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Range it is true according to signal oscillating type
It is fixed;
Coefficient value range in decaying dc atom library function isAttenuation coefficient ρ3Range according to
Concrete signal feature determines.
Further, the step 2) includes the following steps:
21) it to time-frequency Atomic Decomposition, and is reconstructed;
22) by evolution matching pursuit algorithm, the parameter numerical value of each atom is obtained;
23) Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current is calculated separately;
After short circuit occurs after about tens cycles, stable state fundamental current and noise electricity are contained only in short circuit current
Stream, stable state fundamental current value is by synchronous motor d-axis reactance xdIt determines;
Atom decomposition extraction stable state fundametal compoment atomic parameter is utilized after intercepting short circuit current signal.
Further, the step 21) particular content is as follows:
Signal i (t) ∈ H to be analyzed, wherein H indicates the space Hilbert, and D is over-complete dictionary of atoms space, and D ∈ H;
If the corresponding discrete system signal of signal is i (n), grFor the atom of D, r is atomic combination, and Γ is coefficient sets
Set, r ∈ Γ;Atom is normalized, i.e., | | gr| |=1;It is selected from D and signal i (n) original the most matched
Sub- gr(0), i.e. gr(0) it is the atom in D with i (n) inner product maximum absolute value:< in formula,
> indicates the inner product of the two;
Signal i (n) is decomposed into gr(0) component and residual components two parts on: i=< i, gr(0) > gr(0)+
ri 1, r in formulai 1To carry out residual components after the 1st subatom decomposes to signal i (n);
Decomposition, iterative formula r are iterated to the residual components after each decomposei m=< ri m,gr(m) > gr(m)+ri m +1, gr (m) satisfaction
After carrying out n times iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikIt is denoted as reconstruction signal, the similarity with original signal i is
Further, the step 22) particular content is as follows:
221) initialization population: initial population { x is definedij| i=1,2 ..., NP;J=1,2 ..., E1, wherein NP is kind
Group's scale, E1To optimize space dimensionality;
It pressesEach individual is randomly generated, rand is in formulaIt is interior with
Machine number function,For xijLower limit,For xijUpper limit value;
222) it calculates fitness value: each population value being substituted into optimization fitness function, calculates each population fitness value;
223) mutation operation: three different individuals are chosen from current pth for population, variance component is obtained by difference strategy
ForIn formula: i ≠ r1 ≠ r2 ≠ r3, F are scale factor;
224) crossover operation: to current pth for populationAnd its variation vector Vi p+1Crossover operation is carried out, test is obtained
Body In formula: CR is the crossover probability factor, jrandFor [1,
NP] interior random parameter;
225) selection operation:WithBetween next-generation individual is selected by competing In formula: O () is the fitness function for minimizing optimization problem;
226) step 222) is repeated to 225), when reaching maximum number of iterations or similarity reaches certain value, stops changing
In generation, exports the parameter numerical value of each atom.
Further, the formula of calculating Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current is in the step 23)
Wherein, x " d is d axis subtranient reactance, and x'd is the transient state reactance of d axis, x "qFor q axis subtranient reactance, T "dFor d axis
Transient state damping time constant, T'dFor d axis transient state damping time constant, TaFor the time constant of stator winding, f0For electric voltage frequency,For initial phase angle;For fundametal compoment,For harmonic, i01, i02, i03For direct current
Component, e1(t), e2(t), e3It (t) is noise current component, E indicates excitation electric gesture.
Further, F=0.5 in the step 223).
Further, CR=0.9 in the step 224).
Further, the particular content of the step 3) are as follows: enableWherein, it is matched with evolving
Tracing algorithm can carry out Atomic Decomposition to short circuit current, obtain the amplitude A of each atom00, A01, A02, A1, A2, damping time constant
T″d,T′d,Ta, frequency f0And initial phase angle
Simultaneously x " can be obtained using the method being averaged in simultaneous solutiond、x′d、x″q、xd, each time constant passes through atom point
Atomic parameter after solution directly acquires.
Advantageous effects of the invention: be directed to traditional algorithm deficiency, in particular by empirical mode decomposition
(EMD) and Prony algorithm is not good enough to the discomposing effect of Noise signal shorts electric current, and noiseproof feature is bad, this algorithm can be accurate
Ground extracts parameter of synchronous machine, and has preferable noiseproof feature.
Detailed description of the invention
Fig. 1 is the algorithm flow chart of the embodiment of the present invention;
Fig. 2 (a) (b) (c) (d) (e) (f) (g) is that the three short circuit current waveform of the not Noise of the embodiment of the present invention is former
Sub- exploded view respectively corresponds original signal figure, DC component figure, steady-state short-circuit fundamental wave figure, decaying 1 figure of fundamental wave, decaying fundamental wave 2
Figure, harmonic figure and reconstruction signal figure;
Fig. 3 (a) (b) (c) (d) (e) (f) (g) is the three short circuit current waveform atom of the Noise of the embodiment of the present invention
Exploded view respectively corresponds original signal figure, DC component figure, steady-state short-circuit fundamental wave schematic diagram, decaying 1 figure of fundamental wave, decaying fundamental wave 2
Figure, harmonic figure and reconstruction signal figure.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following embodiment is only used for clearly illustrating the present invention
Technical solution, and not intended to limit the protection scope of the present invention.
Specific step is as follows for this algorithm:
Step 1) constructs relevant atomic library, by atom discretization continuous parameters, particular content are as follows:
Step 11) fundamental wave is that the main component of sudden short-circuit current needs to accurately analyze its parameter of synchronous machine
Fundametal compoment is extracted, Fundamental atom library is constructedFormula
In:
f1For fundamental frequency,For fundamental phase, ts1At the beginning of indicating fundamental wave, te1It indicates
The finish time of fundamental wave;U (t) is unit jump function;kr1To make | | gr1(t) | |=1 coefficient.
Coefficient value range in Fundamental atom library function are as follows: 49HZ≤f1≤ 51HZ, Wherein, N indicates the sampled data length of signal i (t) to be analyzed, fsIndicate signal i to be analyzed
(t) sample frequency.
Step 12) is using the fundamental wave sinusoidal signal exponentially to decay in short circuit current and decaying frequency multiplication sinusoidal signal as resistance
Buddhist nun's sine atom constructs damped sine atom are as follows:
In formula:f2For signal frequency,For phase, ρ2For attenuation parameter, ts2At the beginning of indicating disturbance
It carves, te2Indicate the finish time of disturbance, u (t) is unit jump function, kr2To make | | gr2(t) | |=1 coefficient.
Coefficient value range in damped sine atom library function are as follows:By finger
The value range of the decaying fundamental frequency of the fundamental wave sinusoidal signal of number rule decaying is 49HZ≤f2≤ 51HZ, decaying frequency multiplication are sinusoidal
The value range of the component frequencies of signal is 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Range it is true according to signal oscillating type
It is fixed.
Step 13) extracts the attenuating dc component in current signal, structural damping direct current atomIn formula: r3=[ρ3,ts3,te3], ρ3For attenuation coefficient, ts3Expression is opened
Begin the time, te3Indicate the end time, u (t) is unit jump function, kr3To make | | gr3(t) | |=1 coefficient.
Coefficient value range in decaying dc atom library function isAttenuation coefficient ρ3Range according to
Concrete signal feature determines.
Step 2) use evolution matching pursuit algorithm, quickly from synchronous motor sudden short-circuit current extract fundamental current,
DC current and frequency multiplication electric current, specifically comprise the following steps:
21) to time-frequency Atomic Decomposition, and be reconstructed: signal i (t) ∈ H to be analyzed, wherein H indicates that Hilbert is empty
Between, D is over-complete dictionary of atoms space, and D ∈ H;
If the corresponding discrete system signal of signal is i (n), grFor the atom of D, r is atomic combination, and Γ is coefficient sets
Set, r ∈ Γ;Atom is normalized, i.e., | | gr| |=1;It is selected from D and signal i (n) original the most matched
Sub- gr(0), i.e. gr(0) it is the atom in D with i (n) inner product maximum absolute value:< in formula,
> indicates the inner product of the two;
Signal i (n) is decomposed into gr(0) component and residual components two parts on: i=< i, gr(0) > gr(0)+
ri 1, r in formulai 1To carry out residual components after the 1st subatom decomposes to signal i (n);
Decomposition, iterative formula r are iterated to the residual components after each decomposei m=< ri m,gr(m) > gr(m)+ri m +1, gr(m) meet
After carrying out n times iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikIt is denoted as reconstruction signal, the similarity with original signal i is
22) evolution matching pursuit algorithm obtains the parameter numerical value of each atom, comprises steps that:
221) initialization population: initial population { x is definedij| i=1,2 ..., NP;J=1,2 ..., E1, wherein NP is kind
Group's scale, E1To optimize space dimensionality;
It pressesEach individual is randomly generated, rand is in formulaIt is interior with
Machine number,For xijLower and upper limit value;
222) it calculates fitness value: each population value being substituted into optimization fitness function, calculates each population fitness value;
223) mutation operation: three different individuals are chosen from current pth for population, variance component is obtained by difference strategy
ForIn formula: i ≠ r1 ≠ r2 ≠ r3, F are scale factor, F=0.5 in the present embodiment.
224) crossover operation: to current pth for populationAnd its variation vector Vi p+1Crossover operation is carried out, is tested
Individual In formula: CR is the crossover probability factor, jrandFor
[1, NP] interior random parameter, CR=0.9 in the present embodiment.
225) selection operation:WithBetween next-generation individual is selected by competing In formula: O () is the fitness function for minimizing optimization problem;
226) step 222) is repeated to 225), when reaching maximum number of iterations or similarity reaches certain value, stops changing
In generation, exports the parameter numerical value of each atom.
23) Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current is calculated separately, formula is
Wherein, x "dFor d axis subtranient reactance, x'dFor d axis transient state reactance, x "qFor q axis subtranient reactance, T "dFor d axis
Transient state damping time constant, T'dFor d axis transient state damping time constant, TaFor the time constant of stator winding, f0For electric voltage frequency,For initial phase angle.For fundametal compoment,For harmonic, i01, i02, i03For direct current
Component, e1(t), e2(t), e3It (t) is noise current component, E indicates excitation electric gesture.
For above-mentioned formula, as t →+∞,
After short circuit occurs after about tens cycles, stable state fundamental current and noise electricity are contained only in short circuit current
Stream, stable state fundamental current value is by synchronous motor d-axis reactance xdIt determines;
Atom decomposition extraction stable state fundametal compoment atomic parameter is utilized after intercepting short circuit current signal.Short circuit current signal
Energy value shared by middle DC component is big, can preferentially decompose during Atomic Decomposition and extract its characteristic value, and harmonic is wrapped
Energy content is smaller, then can decompose extraction.
Step 3) extrapolates parameter of synchronous machine, particular content from each atom signals characteristic parameter of decomposition are as follows: enablesWherein, Atomic Decomposition can be carried out to short circuit current with evolution matching pursuit algorithm, obtain each atom
Amplitude A00, A01, A02, A1, A2, damping time constant Td ", Td ', Ta, frequency f0And initial phase angleSimultaneous solution and using taking
Xd ", xd ', xq ", xd can be obtained in the method for average value, and each time constant is directly acquired by the atomic parameter after Atomic Decomposition.
According to fig. 2 with the embodiment of Fig. 3, compared respectively with EMD+Prony method in the prior art:
Not parameter identification result comparison sheet when Noise
Parameter identification result comparison sheet when Noise
The present invention constructs relevant atomic library and first by atom discretization continuous parameters, calculates with evolution match tracing
Method quickly extracts fundamental current, DC current, frequency multiplication electric current from synchronous motor sudden short-circuit current;Then from each of decomposition
Parameter of synchronous machine is extrapolated in atom signals characteristic parameter.By intercepting the sampled signal of steady-state shortcircuit current, pick out same
Walk the synchronous reactance and short-circuit first phase angle value of motor.Atomic Decomposition method can accurately extract parameter of synchronous machine, and have preferably
Noiseproof feature.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations
Also it should be regarded as protection scope of the present invention.
Claims (9)
1. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, characterized in that include the following steps:
1) relevant atomic library is constructed, by atom discretization continuous parameters;
2) evolution matching pursuit algorithm is used, fundamental current, DC current are quickly extracted from synchronous motor sudden short-circuit current
With frequency multiplication electric current;
3) parameter of synchronous machine, particular content are extrapolated from each atom signals characteristic parameter of decomposition are as follows:
It enablesWherein, Atomic Decomposition can be carried out to short circuit current with evolution matching pursuit algorithm, obtain
The amplitude A of each atom00, A01, A02, A1, A2, damping time constant T "d,T′d,Ta, frequency f0And initial phase angle
Simultaneously x " can be obtained using the method being averaged in simultaneous solutiond、x′d、x″q、xd, after each time constant passes through Atomic Decomposition
Atomic parameter directly acquires.
2. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 1, characterized in that institute
State the particular content of step 1) are as follows:
11) Fundamental atom library is constructedIn formula: f1For fundamental wave frequency
Rate,For fundamental phase, ts1At the beginning of indicating fundamental wave, te1Indicate the finish time of fundamental wave;T indicates time variable;
U (t) is unit jump function;kr1To make | | gr1(t) | |=1 coefficient;
12) using the fundamental wave sinusoidal signal exponentially to decay in short circuit current and decaying frequency multiplication sinusoidal signal as damped sine
Atom, the damped sine atom of construction are as follows:
In formula: f2For signal frequency,For phase, ρ2For attenuation parameter, ts2At the beginning of indicating disturbance, te2Indicate the end of disturbance
Moment, u (t) are unit jump function, kr2To make | | gr2(t) | |=1 coefficient;
13) attenuating dc component in current signal, structural damping direct current atom are extracted
In formula: ρ3For attenuation coefficient, ts3Indicate time started, te3Indicate the end time, u (t) is unit jump function, kr3To make | |
gr3(t) | |=1 coefficient.
3. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 2, characterized in that base
Coefficient value range in wave atom library function are as follows: 49HZ≤f1≤ 51HZ, Wherein, N
Indicate the sampled data length of signal i (t) to be analyzed, fsIndicate the sample frequency of signal i (t) to be analyzed;
Coefficient value range in damped sine atom library function are as follows: Exponentially
The value range of the decaying fundamental frequency of the fundamental wave sinusoidal signal of decaying is 49HZ≤f2≤ 51HZ, decaying frequency multiplication sinusoidal signal
The value range of component frequencies is 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Range according to signal oscillating type determine;
Coefficient value range in decaying dc atom library function isAttenuation coefficient ρ3Range is according to specific
Signal characteristic determines.
4. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 1, characterized in that institute
Step 2) is stated to include the following steps:
21) it to time-frequency Atomic Decomposition, and is reconstructed;
22) by evolution matching pursuit algorithm, the parameter numerical value of each atom is obtained;
23) Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current is calculated separately;
After short circuit occurs after about tens cycles, stable state fundamental current and noise current are contained only in short circuit current,
Stable state fundamental current value is by synchronous motor d-axis reactance xdIt determines;
Atom decomposition extraction stable state fundametal compoment atomic parameter is utilized after intercepting short circuit current signal.
5. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 4, characterized in that institute
It is as follows to state step 21) particular content:
Signal i (t) ∈ H to be analyzed, wherein H indicates the space Hilbert, and D is over-complete dictionary of atoms space, and D ∈ H;
If the corresponding discrete system signal of signal is i (n), grFor the atom of D, r is atomic combination, and Γ is the collection of coefficient sets
It closes, r ∈ Γ;Atom is normalized, i.e., | | gr| |=1;It is selected from D and signal i (n) atom g the most matchedr
(0), i.e. gr(0) it is the atom in D with i (n) inner product maximum absolute value:< in formula, > table
Show the inner product of the two;
Signal i (n) is decomposed into gr(0) component and residual components two parts on: i=< i, gr(0) > gr(0)+ri 1, formula
Middle ri 1To carry out residual components after the 1st subatom decomposes to signal i (n);
Decomposition, iterative formula r are iterated to the residual components after each decomposei m=< ri m,gr(m) > gr(m)+ri m+1, gr
(m) meet
After carrying out n times iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikIt is denoted as reconstruction signal, the similarity with original signal i is
6. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 5, characterized in that institute
It is as follows to state step 22) particular content:
221) initialization population: initial population { x is definedij| i=1,2 ..., NP;J=1,2 ..., E1, wherein NP is population rule
Mould, E1To optimize space dimensionality;
It pressesEach individual is randomly generated, rand is in formulaInterior random number letter
Number,For xijLower limit,For xijUpper limit value;
222) it calculates fitness value: each population value being substituted into optimization fitness function, calculates each population fitness value;
223) mutation operation: three different individuals are chosen from current pth for population, obtaining variance component by difference strategy isIn formula: i ≠ r1 ≠ r2 ≠ r3, F are scale factor;
224) crossover operation: to current pth for populationAnd its variation vector Vi p+1Crossover operation is carried out, test individual is obtained In formula: CR is the crossover probability factor, jrandFor [1, NP]
Interior random parameter;
225) selection operation:WithBetween next-generation individual is selected by competing In formula: O () is the fitness function for minimizing optimization problem;
226) step 222) is repeated to 225), when reaching maximum number of iterations or similarity reaches certain value, stops iteration, it is defeated
The parameter numerical value of each atom out.
7. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, characterized in that institute
Stating the middle formula for calculating Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current of step 23) is
Wherein, x "dFor d axis subtranient reactance, x'dFor d axis transient state reactance, x "qFor q axis subtranient reactance, T "dFor d axis time transient state
Damping time constant, T 'dFor d axis transient state damping time constant, TaFor the time constant of stator winding, f0For electric voltage frequency,For
Initial phase angle;It is corresponding fundametal compoment, It is corresponding harmonic, i01,
i02, i03It is corresponding DC component, e1(t), e2(t), e3It (t) is corresponding noise current component, E indicates excitation electric
Gesture.
8. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, characterized in that institute
State F=0.5 in step 223).
9. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, characterized in that institute
State CR=0.9 in step 224).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710112598.7A CN106788092B (en) | 2017-02-28 | 2017-02-28 | A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710112598.7A CN106788092B (en) | 2017-02-28 | 2017-02-28 | A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106788092A CN106788092A (en) | 2017-05-31 |
CN106788092B true CN106788092B (en) | 2019-03-05 |
Family
ID=58959978
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710112598.7A Active CN106788092B (en) | 2017-02-28 | 2017-02-28 | A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106788092B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110069744B (en) * | 2018-01-22 | 2020-11-03 | 北京航空航天大学 | Method for estimating stable value of step response signal of pressure sensor |
CN109379011B (en) * | 2018-09-21 | 2021-11-02 | 武汉科技大学 | Ripple thrust compensation method of permanent magnet synchronous linear servo system based on MP algorithm |
CN109981011B (en) * | 2019-04-03 | 2021-03-09 | 中国水利水电科学研究院 | Generator parameter identification method |
CN110071667B (en) * | 2019-05-10 | 2020-09-25 | 福建省宏闽电力工程监理有限公司 | Doubly-fed wind generator parameter identification method |
CN110376517A (en) * | 2019-07-31 | 2019-10-25 | 国网安徽省电力有限公司电力科学研究院 | Determine the method, system and storage medium of Current Waveform Characteristics parameter |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5271699B2 (en) * | 2005-04-19 | 2013-08-21 | ディ.イー.ショー リサーチ,エルエルシー | A zone method for the calculation of particle interactions. |
CN102854437B (en) * | 2012-08-22 | 2015-02-11 | 广东电网公司电力科学研究院 | Fault line selection method of low current grounding system using time-frequency atom decomposition theory |
CN103646146B (en) * | 2013-12-20 | 2017-03-29 | 武汉大学 | Based on the SVC controller method for designing for improving Atomic Decomposition parameter identification |
-
2017
- 2017-02-28 CN CN201710112598.7A patent/CN106788092B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN106788092A (en) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106788092B (en) | A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition | |
CN107451557B (en) | Power transmission line short-circuit fault diagnosis method based on empirical wavelet transform and local energy | |
CN108469560B (en) | Electromagnetic interference objective complexity evaluation method based on rapid S-transform time-frequency space model | |
Jain et al. | Harmonics estimation in emerging power system: Key issues and challenges | |
Han et al. | A recursive sparse representation strategy for bearing fault diagnosis | |
CN105548739B (en) | A kind of arrester operating state signal processing method | |
CN109100814B (en) | A kind of audio magnetotelluric method signal antinoise method based on noise classification | |
Panigrahi et al. | Hybrid signal processing and machine intelligence techniques for detection, quantification and classification of power quality disturbances | |
Li et al. | Damage localization of stacker’s track based on EEMD-EMD and DBSCAN cluster algorithms | |
CN101919695A (en) | Electrocardiosignal QRS complex detection method based on wavelet transform | |
Du et al. | Weighted low-rank sparse model via nuclear norm minimization for bearing fault detection | |
CN106501602B (en) | A kind of fundamental wave measurement method of parameters based on the separation of sliding window frequency spectrum | |
CN110726875A (en) | New energy flexible direct-current grid-connected transient harmonic detection method and system | |
CN115656817A (en) | Permanent magnet synchronous motor turn-to-turn short circuit fault detection method based on neural network technology | |
Xu et al. | Match-extracting chirplet transform with application to bearing fault diagnosis | |
CN117290788B (en) | Power distribution network fault identification method and system based on improved wavelet transformation algorithm | |
CN103245830B (en) | A kind of harmonic detection method in conjunction with AR Power estimation and nonlinear optimization | |
CN106137184A (en) | Electrocardiosignal QRS complex detection method based on wavelet transformation | |
CN116613751B (en) | Small interference stability analysis method and system for new energy grid-connected system | |
CN108008262B (en) | Rapid partial discharge detection method | |
Srivastava et al. | Harmonic and interharmonic estimation using poincaré filtering assisted ESPRIT method | |
Li et al. | Detection of interharmonics using sparse signal decomposition based on ICA-MP | |
Liu et al. | Classification and identification of electric shock current for safety operation in power distribution network | |
CN109617051A (en) | A kind of New-energy power system low-frequency oscillation parameter identification method | |
Li et al. | Based on short-time fourier transform impulse frequency response analysis in the application of the transformer winding deformation |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20170531 Assignee: Nanjing Tangyi Information Technology Co.,Ltd. Assignor: NANJING INSTITUTE OF TECHNOLOGY Contract record no.: X2024980002000 Denomination of invention: A parameter identification method for synchronous motors based on atomic decomposition method Granted publication date: 20190305 License type: Common License Record date: 20240206 |
|
EE01 | Entry into force of recordation of patent licensing contract |