CN106788092A - 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
- CN106788092A CN106788092A CN201710112598.7A CN201710112598A CN106788092A CN 106788092 A CN106788092 A CN 106788092A CN 201710112598 A CN201710112598 A CN 201710112598A CN 106788092 A CN106788092 A CN 106788092A
- Authority
- CN
- China
- Prior art keywords
- atom
- signal
- current
- parameter
- synchronous machine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 46
- 238000000354 decomposition reaction Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000013016 damping Methods 0.000 claims description 13
- 230000001052 transient effect Effects 0.000 claims description 12
- 238000005457 optimization Methods 0.000 claims description 9
- 238000004804 winding Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 230000035772 mutation Effects 0.000 claims description 3
- 238000002474 experimental method Methods 0.000 claims description 2
- 230000000452 restraining effect Effects 0.000 claims 1
- 239000000284 extract Substances 0.000 abstract description 7
- 230000000694 effects Effects 0.000 abstract description 4
- 230000007812 deficiency Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 2
- 238000012545 processing Methods 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
- 230000006872 improvement Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011160 research 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 invention discloses a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, it is characterized in that, comprise the following steps：1) relevant atomic storehouse is built and by atom discretization continuous parameters；2) evolution matching pursuit algorithm is used, fundamental current, DC current, frequency multiplication electric current is quickly extracted from synchronous motor sudden short-circuit current；3) parameter of synchronous machine is extrapolated from each atom signals characteristic parameter for decomposing.The beneficial effect that the present invention is reached：For the deficiency of traditional algorithm, in particular by empirical mode decomposition (EMD) and Prony algorithms it is not good enough to the discomposing effect of Noise signal shorts electric current, noiseproof feature is bad, and this algorithm can exactly extract parameter of synchronous machine, and have preferable noiseproof feature.
Description
Technical field
The present invention relates to a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, belong to signal processing technology neck
Domain.
Background technology
Accurate PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE has important meaning to researching and analysing Operation of Electric Systems and design control system
Justice.Wherein, the transient parameter of reflection synchronous motor transient process is calculated and event with power equipment selection, stability of power system
Barrier calculating etc. is closely related.
The parameter of synchronous machine measuring method that GB/T1029-2005 is recommended is three-phase suddenly-applied short circuit method, by short circuit
The upper and lower envelope extracting cycle component of electric current and aperiodic component.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 functions, the parameter of extraction
Effect is preferable, but Prony algorithms there are problems that determining noise-sensitive and exponent number.
For the limitation of conventional method, and influence of the noise to PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE precision, it is proposed that many
Improved method, but these improved methods are simultaneously also there are still some drawbacks：
The method that wavelet transformation and Prony are combined, information gathering noise is reduced using wavelet transformation, the parameter for obtaining
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 carried out into Subspace partition and total least square (total
Least square, TLS) dual denoising Processing, anti-noise ability is improved.
PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE side based on Hilbert-Huang transform (Hilbert-Huang transform, HHT)
Method, based on empirical mode decomposition (empirical mode decomposition, EMD), constitutes a kind of new time domain filter
Wave method, obtains preferable identification result under strong noise background.But there is insoluble " end points effect in the EMD signals of HHT
Should " problem.
The discrimination method of (local mean decomposition, LMD) is decomposed based on local mean value, identification result is made an uproar
Sound shadow rings smaller, but the moving average span selection of LMD also needs further research.
The content of the invention
To solve the deficiencies in the prior art, it is an object of the invention to provide a kind of synchronous motor based on atom decomposition
Parameter identification method, intactly shows each local feature of non-stationary signal, can cancelling noise well in parsing
Information, improves the precision of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE.
In order to realize above-mentioned target, the present invention is adopted the following technical scheme that：
A kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, it is characterized in that, comprise the following steps：
1) relevant atomic storehouse is built, by atom discretization continuous parameters；
2) evolution matching pursuit algorithm is used, fundamental current, direct current is 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 for decomposing.
Further, the step 1) particular content be：
11) construction Fundamental atom storehouseIn formula：
f1It is fundamental frequency,It is fundamental phase, ts1Carved at the beginning of representing fundamental wave, te1Represent the finish time of fundamental wave；T represents the time
Variable；
U (t) is unit jump function；kr1To make | | gr1(t) | |=1 coefficient；
12) using the fundamental wave sinusoidal signal exponentially decayed in short circuit current and the frequency multiplication sinusoidal signal as damping that decays
Sinusoidal atom, the damped sine atom of construction is：
In formula：f2It is signal frequency,It is phase, ρ2It is attenuation parameter, ts2At the beginning of representing disturbance
Carve, te2The finish time of disturbance is represented, 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], ρ3It is attenuation coefficient, ts3Expression is opened
Time beginning, te3The end time is represented, u (t) is unit jump function, kr3To make | | gr3(t) | |=1 coefficient；
Further, the coefficient value scope in Fundamental atom built-in function is：49HZ≤f1≤ 51HZ, Wherein, N represents the sampled data length of signal i (t) to be analyzed, fsRepresent signal i to be analyzed
The sample frequency of (t)；
Coefficient value scope in damped sine atom built-in function is：By finger
The span of the decay fundamental frequency of the fundamental wave sinusoidal signal of number rule decay is 49HZ≤f2≤ 51HZ, decay frequency multiplication is sinusoidal
The span of the component frequencies of signal is 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Scope it is true according to signal oscillating type
It is fixed；
Coefficient value scope in decaying dc atom built-in function isAttenuation coefficient ρ3Scope according to
Concrete signal feature determines.
Further, the step 2) comprise the following steps：
21) time-frequency atom is decomposed, 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 respectively；
By after about tens cycles, comprising only stable state fundamental current and noise electricity in short circuit current after short circuit occurs
Stream, stable state fundamental current value is by synchronous motor d-axis reactance xdDetermine；
After interception short circuit current signal stable state fundametal compoment atomic parameter is extracted using atom decomposition.
Further, the step 21) particular content is as follows：
Signal i (t) ∈ H to be analyzed, wherein H represent Hilbert spaces, and D is over-complete dictionary of atoms space, and D ∈ H；
If the corresponding discrete system signal of signal is i (n), grIt is the atom of D, r is combined for atomic, and Γ is coefficient sets
Set, r ∈ Γ；Atom is normalized, i.e., | | gr| |=1；The original matched the most with signal i (n) is selected from D
Sub- gr(0), i.e. gr(0)It is the atom in D with i (n) inner product maximum absolute values：＜ in formula, ＞
Represent both inner products；
Signal i (n) is decomposed into gr(0)On component and residual components two parts：I=＜ i, gr(0)＞ gr(0)+ri 1, formula
Middle ri 1It is that residual components after the 1st subatom is decomposed are carried out to signal i (n)；
Decomposition is iterated to the residual components after decomposing every time, iterative formula is ri m=＜ ri m,gr(m)＞ gr(m)+ri m +1, gr(m)Meet
After carrying out n iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikReconstruction signal is designated as, the similarity with primary signal i is
Further, the step 22) particular content is as follows：
221) population is initialized：Define initial population { xij| i=1,2 ..., NP；J=1,2 ..., E1, wherein NP is to plant
Group's scale, E1It is optimization space dimensionality；
PressRandomly generate each individual, rand is in formulaIt is interior with
Machine number function,It is xijLower limit,It is xijHigher limit；
222) fitness value is calculated：Each population value is substituted into optimization fitness function, each population fitness value is calculated；
223) mutation operation：Three different individualities are chosen for population from current pth, 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, experiment is obtained
Body In formula：CR is the crossover probability factor, jrandFor [1,
NP] interior random parameter；
225) selection operation：WithBetween individuality of future generation is selected by competing In formula：O () is the fitness function for minimizing optimization problem；
226) repeat step 222) to 225), when maximum iteration is reached or similarity reaches certain value, stop changing
In generation, export the parameter numerical value of each atom.
Further, the step 23) in calculate Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current formula be
Wherein, x "dIt is d axle subtranient reactances, x'dIt is d axle transient state reactance, x "qIt is q axle subtranient reactances, Td" it is d axles
Transient state damping time constant, Td' it is d axle transient state damping time constants, TaIt is the time constant of stator winding, f0It is electric voltage frequency,It is initial phase angle；It is fundametal compoment,It is harmonic, i01, i02, i03It is direct current
Component, e1(t), e2(t), e3T () is noise current component, E represents excitation electric gesture.
Further, the step 223) in F=0.5.
Further, the step 224) in CR=0.9.
Further, the step 3) particular content be：OrderWherein, 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
Td″,Td′,Ta, frequency f0And initial phase angle
Simultaneous solution simultaneously can obtain x using the method averaged "d、x′d、xq″、xd, each time constant is by atom point
Atomic parameter direct access after solution.
The beneficial effect that the present invention is reached：For the deficiency of traditional algorithm, in particular by empirical mode decomposition
(EMD) and Prony algorithms are not good enough to the discomposing effect of Noise signal shorts electric current, noiseproof feature is bad, and this algorithm can be accurate
Ground extracts parameter of synchronous machine, and has preferable noiseproof feature.
Brief description of the drawings
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, corresponds to primary signal figure, DC component figure, steady-state short-circuit fundamental wave figure, decay fundamental wave 1 figure, decay fundamental wave 2 respectively
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, corresponds to primary signal figure, DC component figure, steady-state short-circuit fundamental wave schematic diagram, decay fundamental wave 1 figure, decay fundamental wave 2 respectively
Figure, harmonic figure and reconstruction signal figure.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following examples are only used for clearly illustrating the present invention
Technical scheme, and can not be limited the scope of the invention with this.
This algorithm is comprised the following steps that：
Step 1) relevant atomic storehouse is built, by atom discretization continuous parameters, particular content is：
Step 11) fundamental wave be sudden short-circuit current main component, in order to analyze its parameter of synchronous machine exactly, it is necessary to
Fundametal compoment is extracted, construction Fundamental atom storehouseFormula
In：
f1It is fundamental frequency,It is fundamental phase, ts1Carved at the beginning of representing fundamental wave, te1Represent
The finish time of fundamental wave；U (t) is unit jump function；kr1To make | | gr1(t) | |=1 coefficient.
Coefficient value scope in Fundamental atom built-in function is：49HZ≤f1≤ 51HZ, Wherein, N represents the sampled data length of signal i (t) to be analyzed, fsRepresent signal i to be analyzed
The sample frequency of (t).
Step 12) using the fundamental wave sinusoidal signal exponentially decayed in short circuit current and decay frequency multiplication sinusoidal signal as
Damped sine atom, constructing damped sine atom is：
In formula：f2It is signal frequency,It is phase, ρ2It is attenuation parameter, ts2At the beginning of representing disturbance
Carve, te2The finish time of disturbance is represented, u (t) is unit jump function, kr2To make | | gr2(t) | |=1 coefficient.
Coefficient value scope in damped sine atom built-in function is：By finger
The span of the decay fundamental frequency of the fundamental wave sinusoidal signal of number rule decay is 49HZ≤f2≤ 51HZ, decay frequency multiplication is sinusoidal
The span of the component frequencies of signal is 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Scope it is true according to signal oscillating type
It is fixed.
Step 13) extract current signal in attenuating dc component, structural damping direct current atomIn formula：R3=[ρ3,ts3,te3], ρ3It is attenuation coefficient, ts3Expression is opened
Time beginning, te3The end time is represented, u (t) is unit jump function, kr3To make | | gr3(t) | |=1 coefficient.
Coefficient value scope in decaying dc atom built-in function isAttenuation coefficient ρ3Scope 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 include following steps：
21) time-frequency atom is decomposed, and is reconstructed：Signal i (t) ∈ H, wherein H to be analyzed represents 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), grIt is the atom of D, r is combined for atomic, and Γ is coefficient sets
Set, r ∈ Γ；Atom is normalized, i.e., | | gr| |=1；The original matched the most with signal i (n) is selected from D
Sub- gr(0), i.e. gr(0)It is the atom in D with i (n) inner product maximum absolute values：＜ in formula, ＞
Represent both inner products；
Signal i (n) is decomposed into gr(0)On component and residual components two parts：I=＜ i, gr(0)＞ gr(0)+ri 1, formula
Middle ri 1It is that residual components after the 1st subatom is decomposed are carried out to signal i (n)；
Decomposition is iterated to the residual components after decomposing every time, iterative formula is ri m=＜ ri m,gr(m)＞ gr(m)+ri m +1, gr(m)Meet
After carrying out n iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikReconstruction signal is designated as, the similarity with primary signal i is
22) evolution matching pursuit algorithm, obtains the parameter numerical value of each atom, comprises steps that：
221) population is initialized：Define initial population { xij| i=1,2 ..., NP；J=1,2 ..., E1, wherein NP is to plant
Group's scale, E1It is optimization space dimensionality；
PressRandomly generate each individual, rand is in formulaIt is interior with
Machine number,It is xijLower and upper limit value；
222) fitness value is calculated：Each population value is substituted into optimization fitness function, each population fitness value is calculated；
223) mutation operation：Three different individualities are chosen for population from current pth, 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
It is 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 individuality of future generation is selected by competing In formula：O () is the fitness function for minimizing optimization problem；
226) repeat step 222) to 225), when maximum iteration is reached or similarity reaches certain value, stop changing
In generation, export the parameter numerical value of each atom.
23) Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current is calculated respectively, and formula is
Wherein, x "dIt is d axle subtranient reactances, x'dIt is d axle transient state reactance, x "qIt is q axle subtranient reactances, Td" it is d axles
Transient state damping time constant, Td' it is d axle transient state damping time constants, TaIt is the time constant of stator winding, f0It is electric voltage frequency,It is initial phase angle.It is fundametal compoment, It is harmonic, i01, i02, i03It is direct current
Component, e1(t), e2(t), e3T () is noise current component, E represents excitation electric gesture.
For above-mentioned formula, as t →+∞,
By after about tens cycles, comprising only stable state fundamental current and noise electricity in short circuit current after short circuit occurs
Stream, stable state fundamental current value is by synchronous motor d-axis reactance xdDetermine；
After interception short circuit current signal stable state fundametal compoment atomic parameter is extracted using atom decomposition.Short circuit current signal
Energy value shared by middle DC component is big, can preferentially be decomposed during Atomic Decomposition and extract its characteristic value, and harmonic is wrapped
Energy content is smaller, can then decompose extraction.
Step 3) parameter of synchronous machine is extrapolated from each atom signals characteristic parameter for decomposing, particular content is：OrderWherein, Atomic Decomposition can be carried out to short circuit current, obtain each atom with evolution matching pursuit algorithm
Amplitude A00, A01, A02, A1, A2, damping time constant Td ", Td ', Ta, frequency f0And initial phase angleSimultaneous solution and using taking
The method of average value can obtain xd ", xd ', xq ", xd, and each time constant is by the atomic parameter direct access after Atomic Decomposition.
According to the embodiment of Fig. 2 and Fig. 3, contrasted with EMD+Prony methods of the prior art respectively：
Not parameter identification result comparison sheet during Noise
Parameter identification result comparison sheet during Noise
The present invention builds relevant atomic storehouse and by atom discretization continuous parameters first, is calculated with evolution match tracing
Method, quickly extracts fundamental current, DC current, frequency multiplication electric current from synchronous motor sudden short-circuit current；Then it is each from what is decomposed
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 exactly extract parameter of synchronous machine, and have preferably
Noiseproof feature.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, on the premise of the technology of the present invention principle is not departed from, some improvement and deformation can also be made, these improve and deform
Also should be regarded as protection scope of the present invention.
Claims (10)
1. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition, it is characterized in that, comprise the following steps：
1) relevant atomic storehouse is built, by atom discretization continuous parameters；
2) evolution matching pursuit algorithm is used, fundamental current, DC current is quickly extracted from synchronous motor sudden short-circuit current
With frequency multiplication electric current；
3) parameter of synchronous machine is extrapolated from each atom signals characteristic parameter for decomposing.
2. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 1, it is characterized in that, institute
State step 1) particular content be：
11) construction Fundamental atom storehouseIn formula：
f1It is fundamental frequency,It is fundamental phase, ts1Carved at the beginning of representing fundamental wave, te1Represent the finish time of fundamental wave；T represents the time
Variable；
U (t) is unit jump function；kr1To make | | gr1(t) | |=1 coefficient；
12) using the fundamental wave sinusoidal signal exponentially decayed in short circuit current and the frequency multiplication sinusoidal signal as damped sine that decays
Atom, the damped sine atom of construction is：
In formula：f2It is signal frequency,It is phase, ρ2It is attenuation parameter, ts2At the beginning of representing disturbance
Carve, te2The finish time of disturbance is represented, 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 extracted
In formula：R3=[ρ3,ts3,te3], ρ3It is attenuation coefficient, ts3Represent time started, te3The end time is represented, u (t) is unit rank
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, it is characterized in that, base
Coefficient value scope in ripple atom built-in function is：49HZ≤f1≤ 51HZ, Wherein, N
Represent the sampled data length of signal i (t) to be analyzed, fsRepresent the sample frequency of signal i (t) to be analyzed；
Coefficient value scope in damped sine atom built-in function is：Advised by index
The span for restraining the decay fundamental frequency of the fundamental wave sinusoidal signal of decay is 49HZ≤f2≤ 51HZ, decay frequency multiplication sinusoidal signal
Component frequencies span be 98HZ≤f2≤ 101HZ, attenuation parameter ρ2Scope according to signal oscillating type determine；
Coefficient value scope in decaying dc atom built-in function isAttenuation coefficient ρ3Scope 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, it is characterized in that, institute
State step 2) comprise the following steps：
21) time-frequency atom is decomposed, 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 respectively；
By after about tens cycles, stable state fundamental current and noise current being comprised only in short circuit current after short circuit occurs,
Stable state fundamental current value is by synchronous motor d-axis reactance xdDetermine；
After interception short circuit current signal stable state fundametal compoment atomic parameter is extracted using atom decomposition.
5. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 4, it is characterized in that, institute
State step 21) particular content is as follows：
Signal i (t) ∈ H to be analyzed, wherein H represent Hilbert spaces, and D is over-complete dictionary of atoms space, and D ∈ H；
If the corresponding discrete system signal of signal is i (n), grIt is the atom of D, r is combined for atomic, Γ is the collection of coefficient sets
Close, r ∈ Γ；Atom is normalized, i.e., | | gr| |=1；The atom matched the most with signal i (n) is selected from D
gr(0), i.e. gr(0)It is the atom in D with i (n) inner product maximum absolute values：＜ in formula, ＞ are represented
Both inner products；
Signal i (n) is decomposed into gr(0)On component and residual components two parts：I=＜ i, gr(0)＞ gr(0)+ri 1, r in formulai 1
It is that residual components after the 1st subatom is decomposed are carried out to signal i (n)；
Decomposition is iterated to the residual components after decomposing every time, iterative formula is ri m=＜ ri m,gr(m)＞ gr(m)+ri m+1, gr(m)
Meet
After carrying out n iteration, current residual components | | ri m| | it is sufficiently small or decay to 0, then this signal i can approximate representation beikReconstruction signal is designated as, the similarity with primary signal i is
6. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 5, it is characterized in that, institute
State step 22) particular content is as follows：
221) population is initialized：Define initial population { xij| i=1,2 ..., NP；J=1,2 ..., E1, wherein NP is advised for population
Mould, E1It is optimization space dimensionality；
PressRandomly generate each individual, rand is in formulaInterior random number letter
Number,It is xijLower limit,It is xijHigher limit；
222) fitness value is calculated：Each population value is substituted into optimization fitness function, each population fitness value is calculated；
223) mutation operation：Three different individualities are chosen for population from current pth, 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, experiment is obtained individual In formula：CR is the crossover probability factor, jrandIt is [1, NP]
Interior random parameter；
225) selection operation：WithBetween individuality of future generation is selected by competing In formula：O () is the fitness function for minimizing optimization problem；
226) repeat step 222) to 225), when maximum iteration is reached or similarity reaches certain value, stop iteration, it is defeated
Go out the parameter numerical value of each atom.
7. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, it is characterized in that, institute
State step 23) in calculate Damper Winding synchronous motor Sudden Three-phase Short Circuit electric current formula be
Wherein, x "dIt is d axle subtranient reactances, x'dIt is d axle transient state reactance, x "qIt is q axle subtranient reactances, T "dIt is d axles time transient state
Damping time constant, T 'dIt is d axle transient state damping time constants, TaIt is the time constant of stator winding, f0It is electric voltage frequency,For
Initial phase angle；Corresponding fundametal compoment is, It is corresponding harmonic, i01,
i02, i03It is corresponding DC component, e1(t), e2(t), e3T () is corresponding noise current component, E represents excitation electric
Gesture.
8. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, it is characterized in that, institute
State step 223) in F=0.5.
9. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 6, it is characterized in that, institute
State step 224) in CR=0.9.
10. a kind of PARAMETER IDENTIFICATION OF SYNCHRONOUS MACHINE method based on atom decomposition according to claim 1, it is characterized in that,
The step 3) particular content be：OrderWherein, can be to short-circuit electricity with evolution matching pursuit algorithm
Stream carries out Atomic Decomposition, obtains the amplitude A of each atom00, A01, A02, A1, A2, damping time constant Td″,Td′,Ta, frequency f0With
Initial phase angle
Simultaneous solution simultaneously can obtain x using the method averaged "d、x′d、x″q、xd, each time constant is by after Atomic Decomposition
Atomic parameter direct access.
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 true CN106788092A (en) | 2017-05-31 |
| CN106788092B 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) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109379011A (en) * | 2018-09-21 | 2019-02-22 | 武汉科技大学 | Permanent magnet synchronous linear motor servo ripple thrust-compensating method based on MP algorithm |
| CN109981011A (en) * | 2019-04-03 | 2019-07-05 | 中国水利水电科学研究院 | A kind of generator parameter identification method |
| CN110069744A (en) * | 2018-01-22 | 2019-07-30 | 北京航空航天大学 | A kind of estimation method of pressure sensor phase step response signals stationary value |
| CN110071667A (en) * | 2019-05-10 | 2019-07-30 | 福建省宏闽电力工程监理有限公司 | 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 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120116737A1 (en) * | 2005-04-19 | 2012-05-10 | D.E. Shaw Research LLC | Zonal methods for computation of particle interactions |
| CN102854437A (en) * | 2012-08-22 | 2013-01-02 | 广东电网公司电力科学研究院 | Fault line selection method of low current grounding system using time-frequency atom decomposition theory |
| CN103646146A (en) * | 2013-12-20 | 2014-03-19 | 武汉大学 | Design method for SVC controller based on improved atomic decomposition parameter identification |
-
2017
- 2017-02-28 CN CN201710112598.7A patent/CN106788092B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120116737A1 (en) * | 2005-04-19 | 2012-05-10 | D.E. Shaw Research LLC | Zonal methods for computation of particle interactions |
| CN102854437A (en) * | 2012-08-22 | 2013-01-02 | 广东电网公司电力科学研究院 | Fault line selection method of low current grounding system using time-frequency atom decomposition theory |
| CN103646146A (en) * | 2013-12-20 | 2014-03-19 | 武汉大学 | Design method for SVC controller based on improved atomic decomposition parameter identification |
Non-Patent Citations (1)
| Title |
|---|
| 杨宗军等: "《表贴式永磁同步电机的多参数在线辨识》", 《电工技术学报》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110069744A (en) * | 2018-01-22 | 2019-07-30 | 北京航空航天大学 | A kind of estimation method of pressure sensor phase step response signals stationary value |
| CN110069744B (en) * | 2018-01-22 | 2020-11-03 | 北京航空航天大学 | Method for estimating stable value of step response signal of pressure sensor |
| CN109379011A (en) * | 2018-09-21 | 2019-02-22 | 武汉科技大学 | Permanent magnet synchronous linear motor servo ripple thrust-compensating method based on MP algorithm |
| CN109379011B (en) * | 2018-09-21 | 2021-11-02 | 武汉科技大学 | Ripple thrust compensation method of permanent magnet synchronous linear servo system based on MP algorithm |
| CN109981011A (en) * | 2019-04-03 | 2019-07-05 | 中国水利水电科学研究院 | A kind of generator parameter identification method |
| CN110071667A (en) * | 2019-05-10 | 2019-07-30 | 福建省宏闽电力工程监理有限公司 | Doubly-fed wind 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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106788092B (en) | 2019-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106788092A (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 | |
| Rosenhaus et al. | Entanglement entropy: a perturbative calculation | |
| Frei et al. | Intrinsic time-scale decomposition: time–frequency–energy analysis and real-time filtering of non-stationary signals | |
| Youssef et al. | Disturbance classification utilizing dynamic time warping classifier | |
| CN111160798B (en) | Non-invasive household appliance load identification method based on bee colony algorithm | |
| CN109100814A (en) | A kind of audio magnetotelluric method signal antinoise method based on noise classification | |
| CN105760347A (en) | HHT end effect restraining method based on data/extreme value joint symmetric prolongation | |
| CN109783919A (en) | Frequency domain multi-conductor transmission lines electromagnetic pulse response quickly modeling method based on waveform relaxation iteration | |
| CN109446975A (en) | Multiple dimensioned noise adjusts the Detection of Weak Signals of accidental resonance | |
| Peng et al. | Short-term wind power prediction based on wavelet feature arrangement and convolutional neural networks deep learning | |
| Guesmi et al. | Fingerprint verification system based on curvelet transform and possibility theory | |
| CN117290788B (en) | Power distribution network fault identification method and system based on improved wavelet transformation algorithm | |
| Hafiz et al. | An approach for classification of power quality disturbances based on Hilbert Huang transform and Relevance vector machine | |
| CN114325072B (en) | Ferromagnetic resonance overvoltage identification method and device based on gram angular field coding | |
| CN109617051A (en) | A kind of New-energy power system low-frequency oscillation parameter identification method | |
| Li et al. | Detection of interharmonics using sparse signal decomposition based on ICA-MP | |
| CN115329820A (en) | Power grid disturbance identification method and system based on improved FastICA | |
| Peng et al. | Modulated broadband mode decomposition for the feature extraction of double pulse metal inert gas welding | |
| Wang et al. | Mill load identification method for ball milling process based on grinding signal | |
| Behzadi et al. | Dual identification of multi-complex and non-stationary power quality disturbances using variational mode decomposition in hybrid modern power systems | |
| Fan et al. | Voltage disturbance signals identification based on ILMD and neural network | |
| CN117390508B (en) | Hydroelectric generating set signal state identification method based on time-shifting multi-scale cosine similarity entropy | |
| Zhuo et al. | The evaluation of wind power time series multi-scale modeling methods | |
| Liao et al. | Identification of power quality disturbances based on sample entropy and weighted optimization random forest |
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 |