CN107064646A - A kind of electric network impedance real-time identification method for multiple-input and multiple-output asymmetric system - Google Patents
A kind of electric network impedance real-time identification method for multiple-input and multiple-output asymmetric system Download PDFInfo
- Publication number
- CN107064646A CN107064646A CN201710367062.XA CN201710367062A CN107064646A CN 107064646 A CN107064646 A CN 107064646A CN 201710367062 A CN201710367062 A CN 201710367062A CN 107064646 A CN107064646 A CN 107064646A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- mtable
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/08—Measuring resistance by measuring both voltage and current
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
- G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
- G01R27/16—Measuring impedance of element or network through which a current is passing from another source, e.g. cable, power line
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
The invention discloses a kind of electric network impedance real-time identification method of multiple-input and multiple-output asymmetric system, including：PCC point three-phase networking electric currents are obtained by current sensor；In dsp respectively by the three-phase current i collecteda,ib,icCoordinate transform is carried out, the electric current under two cordic phase rotators is obtained；Produce the orthogonal-binary sequence with two signals；In the reference current that orthogonal-binary sequence is injected to power network current d axles and q by digital obtaining unit；Measurement line voltage and power frequency respond and are broken down into positive sequence and negative sequence component；Calculate and obtain electric network impedance.The problem of off-line identification stability is not enough in the prior art is solved, and shortens time of measuring.
Description
Technical field
The invention belongs to electric and electronic technical field, it is related to a kind of real-time impedance identification side of multi-input multi-output system
Method.
Background technology
Real-time identification impedance can be roughly divided into parametrization and imparametrization method, and nonparametric frequency domain method is in electric network impedance
Extensive use in identification and stability analysis.Different with parametric method, nonparametric technique needs to assume a system model, and
And only select that an excitation.By excitation injection power network and recognize electric network impedance using it.In stable state, measurement electric current and voltage
Respond and carry out Fourier transformation analysis.
With the change and the limitation of fossil fuel of weather, distributed power generation is becoming increasingly popular, inverter output impedance
The problem of harmonic response that mismatch between electric network impedance causes therebetween is one very main.Research in recent years is all
It has been shown that, the impedance by measuring power network and inverter can avoid unstable situation, and inverse to change based on these parameters
Become the design parameter of device.According to identification technique, off-line identification and on-line identification can be divided into, offline impedance is recognized it is impossible to ensure that being
The stability of system, because the impedance of power network and inverter is real-time change.Research in recent years is also increasingly tended in real time
Impedance is recognized, and real-time impedance identification is compared to offline impedance identification speed faster.Accordingly, it would be desirable to be improved to prior art.
The content of the invention
The present invention proposes a kind of method of the real-time impedance identification based on multi-input multi-output system, solves existing skill
The problem of off-line identification stability is not enough in art, and shorten time of measuring.
To realize above-mentioned technical purpose and the technique effect, the present invention is achieved through the following technical solutions.
Step 1, PCC point three-phase networking electric currents i is obtained respectively by current sensora,ib,ic；
Step 2, the three-phase current i collected in dspa, ib, icCoordinate transform is carried out, is obtained under two cordic phase rotators
Electric current, i.e. id,iq：
Step 3, the orthogonal-binary sequence with two signals is produced：
The generation (d axles Injection Signal) of primary signal：DIBS signal amplitudes switch between being worth at 1 and -1 two, make Ab (t)
For periodicity real number binary sequence.A amplitudes, then the Fourier transform expression formula of its k subharmonic be
Designed DIBS sequence signals are defined as：
By Matlab simulation softwares, optimization problem is solved, and by the sequence under specific frequency with the shape of code table
Formula is stored in dsp, is called during to inject.
Q axles inject the generation of reference signal sequence：The addition sequence 01010 ... of mould 2 forms the input of q axles into primary signal
The disturbance of reference signal, storage in dsp, is called during to inject in the form of code table.
Step 4, while orthogonal-binary sequence is injected into power network current d axles and q reference electricity by digital obtaining unit
In stream；
Step 5, measurement line voltage and power frequency respond and are broken down into positive sequence and negative sequence component, and decomposition formula is such as
Under：
The decomposition of current-responsive：
Wherein p, n represent positive sequence and negative phase-sequence, a=e respectivelyj2π/3
Because ic=-ia-ib, introduce matrix：
Then have：
The decomposition of voltage responsive：
By vca=-vab-vbc, introduce matrix：
Then, the positive-negative sequence of voltage component can be reduced to：
Step 6, the impedance components coupled for positive-negative sequence, have：
Produce electric network impedance：
Wherein, x, y represent the response of two orthogonal-binary sequence injections respectively.
Beneficial effects of the present invention are：(1) existing off-line identification is solved it cannot be guaranteed that the problem of the stability of a system；(2)
Reasonable in design, usability is strong, simple to operate；(3) relative to traditional single input/mono- output e measurement technology, system dynamic is overcome
The shortcoming of characteristic variations；Shorten time of measuring.
Brief description of the drawings
Fig. 1 is the block diagram of real-time impedance device for identifying.
Fig. 2 is impedance discrimination method flow chart of steps.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with accompanying drawing and example, to this
Invention is further elaborated.
As shown in figure 1, direct voltage source passes through three-phase grid-connected inverter, direct current is converted into three-phase alternating current through filtering
After wave apparatus and in power network, current sensor collection networking current information, voltage sensor collection networking information of voltage.
The discrimination method of the present invention is transformed to the magnitude of current under dq coordinates by the way that the networking collected electric current is carried out into coordinate,
Orthogonal-binary sequence is injected at the reference current of d, q axle simultaneously by DSP, analysis power network current and voltage responsive, respectively
Positive-negative sequence composition is broken down into, the positive and negative order components responded according to Current Voltage obtain electric network impedance.
It is the identification step of the present invention, details are as follows for identification process with reference to shown in Fig. 2：
Step 1, PCC point three-phase networking electric currents i is obtained by current sensora,ib,ic；
Step 2, in dsp by the three-phase current i collecteda,ib,icCoordinate transform is carried out, is obtained under two cordic phase rotators
Electric current, i.e. id,iq；
Step 3, the orthogonal-binary sequence with two signals is produced：
The generation (d axles Injection Signal) of primary signal：DIBS signal amplitudes switch between being worth at 1 and -1 two, make Ab (t)
For periodicity real number binary sequence.A is amplitude, then the Fourier transform expression formula of its k subharmonic is：
Designed DIBS sequence signals are defined as：
By Matlab simulation softwares, optimization problem is solved, and by the sequence under specific frequency with the shape of code table
Formula is stored in dsp, is called during to inject.
Q axles inject the generation of reference signal sequence：The addition sequence 01010 ... of mould 2 forms the input of q axles into primary signal
The disturbance of reference signal, storage in dsp, is called during to inject in the form of code table.
Step 4, while orthogonal-binary sequence to be injected to the reference current of power network d axles and q axles by digital obtaining unit；
Step 5, measurement line voltage and power frequency respond and are broken down into positive sequence and negative sequence component is (assuming that three contraries
It is infinite to become device zero-sequence component, i.e., power network connected mode is three-wire system), decomposition formula is as follows：
The decomposition of current-responsive：
P, n represent positive sequence and negative phase-sequence respectively；
Because ic=-ia-ib, introduce matrix：
Then have：
The decomposition of voltage responsive：
By vca=-vab-vbc, introduce matrix：
Then, the positive-negative sequence of voltage component can be reduced to：
Step 6, the impedance components coupled for positive-negative sequence, have：
Produce electric network impedance：
Wherein, x, y represent the response of two orthogonal-binary sequence injections respectively.
Merely illustrating the principles of the invention described in above-described embodiment and specification and most preferred embodiment, are not departing from this
On the premise of spirit and scope, the present invention can also various changes and modifications, these changes and improvements both fall within claimed
The scope of the invention in.
Claims (1)
1. a kind of electric network impedance real-time identification method for multiple-input and multiple-output asymmetric system, direct voltage source passes through inversion
Device, then PCC points are incorporated to after device after filtering, the electric network impedance real-time identification method includes：
Step 1, PCC point three-phase networking electric currents i is obtained by current sensora,ib,ic；
Step 2, in dsp respectively by the three-phase current i collecteda,ib,icCoordinate transform is carried out, is obtained under two cordic phase rotators
Electric current be id,iq；
Step 3, the orthogonal-binary sequence with two signals is produced：The generation (d axles Injection Signal) of primary signal：DIBS believes
Number amplitude switches between being worth at 1 and -1 two, and it is periodicity real number binary sequence to make Ab (t).A amplitudes, then its k subharmonic
Fourier transform expression formula is：
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>b</mi>
<mi>k</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>A</mi>
<mi>T</mi>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</msubsup>
<mi>b</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mi>&pi;</mi>
<mi>j</mi>
<mi>k</mi>
<mi>t</mi>
<mo>/</mo>
<mi>T</mi>
</mrow>
</msup>
<mi>d</mi>
<mi>t</mi>
</mrow>
Designed DIBS sequence signals are defined as：
<mrow>
<mover>
<mi>A</mi>
<mo>^</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mrow>
<mo>|</mo>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
</mrow>
</mrow>
By Matlab simulation softwares, optimization problem is solved, and the sequence under specific frequency is deposited in the form of code table
Storage in dsp, is called during to inject.
Q axles inject the generation of reference signal sequence：The addition sequence 01010 ... of mould 2 forms the input reference of q axles into primary signal
The disturbance of signal, storage in dsp, is called during to inject in the form of code table.
Step 4, while orthogonal-binary sequence to be injected to the reference current of power network current d axles and q axles by digital obtaining unit
In；
Step 5, measurement line voltage and power frequency respond and are broken down into positive sequence and negative sequence component decomposition formula is as follows：
The decomposition of current-responsive：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>a</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>b</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>c</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Wherein p, n represent positive sequence and negative phase-sequence, a=e respectivelyj2π/3
Because ic=-ia-ib, introduce matrix：
<mrow>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>3</mn>
</msqrt>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>6</mn>
</mrow>
</msup>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>6</mn>
</mrow>
</msup>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Then have：
The decomposition of voltage responsive：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>3</mn>
<msqrt>
<mn>3</mn>
</msqrt>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>6</mn>
</mrow>
</msup>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>6</mn>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
</mtd>
<mtd>
<mi>a</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>c</mi>
<mi>a</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
By vca=-vab-vbc, introduce matrix
<mrow>
<msub>
<mi>S</mi>
<mi>v</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msqrt>
<mn>3</mn>
</msqrt>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mo>/</mo>
<mn>3</mn>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Then, the positive and negative order components of voltage component can be reduced to：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msub>
<mi>S</mi>
<mi>v</mi>
</msub>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>a</mi>
<mi>b</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mrow>
<mi>b</mi>
<mi>c</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Step 6, the impedance components coupled for positive-negative sequence, have
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>V</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>p</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>I</mi>
<mi>n</mi>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Produce electric network impedance：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>p</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>n</mi>
<mi>p</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mi>n</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>p</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>p</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>V</mi>
<mrow>
<mi>n</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>n</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>n</mi>
<mi>y</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
Wherein, x, y represent the response of two orthogonal-binary sequence injections respectively.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710367062.XA CN107064646B (en) | 2017-05-23 | 2017-05-23 | Power grid impedance real-time identification method for multi-input multi-output asymmetric system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710367062.XA CN107064646B (en) | 2017-05-23 | 2017-05-23 | Power grid impedance real-time identification method for multi-input multi-output asymmetric system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107064646A true CN107064646A (en) | 2017-08-18 |
| CN107064646B CN107064646B (en) | 2020-02-21 |
Family
ID=59610363
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710367062.XA Active CN107064646B (en) | 2017-05-23 | 2017-05-23 | Power grid impedance real-time identification method for multi-input multi-output asymmetric system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107064646B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109412178A (en) * | 2018-10-19 | 2019-03-01 | 湖南大学 | A kind of switching control method of high-voltage large-capacity impedance measurement equipment |
| CN112567584A (en) * | 2018-08-20 | 2021-03-26 | 西门子股份公司 | Grid impedance estimation for power converter controller |
| CN112861315A (en) * | 2021-01-11 | 2021-05-28 | 广西大学 | One-dimensional descent search method for non-convex single-target optimal power flow global solution of power system |
| CN113743806A (en) * | 2021-09-09 | 2021-12-03 | 广西康耐肯特科技有限公司 | Method for searching non-convex dual-target optimal power flow global solution of power system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007205934A (en) * | 2006-02-02 | 2007-08-16 | Hioki Ee Corp | Inspecting device |
| CN103376363A (en) * | 2013-07-04 | 2013-10-30 | 西安交通大学 | Impedance measurement method of three-phase alternating current system |
| CN103630748A (en) * | 2013-11-11 | 2014-03-12 | 西安交通大学 | Device and method for harmonic impedance measurement of micro-grid |
| CN105158592A (en) * | 2015-07-09 | 2015-12-16 | 中国矿业大学 | Real-time three-phase electrical network impedance detection method on the asymmetric electrical network condition |
| CN105259414A (en) * | 2015-11-23 | 2016-01-20 | 上海电力学院 | Power grid impedance on-line detection method on the basis of inverter |
-
2017
- 2017-05-23 CN CN201710367062.XA patent/CN107064646B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2007205934A (en) * | 2006-02-02 | 2007-08-16 | Hioki Ee Corp | Inspecting device |
| CN103376363A (en) * | 2013-07-04 | 2013-10-30 | 西安交通大学 | Impedance measurement method of three-phase alternating current system |
| CN103630748A (en) * | 2013-11-11 | 2014-03-12 | 西安交通大学 | Device and method for harmonic impedance measurement of micro-grid |
| CN105158592A (en) * | 2015-07-09 | 2015-12-16 | 中国矿业大学 | Real-time three-phase electrical network impedance detection method on the asymmetric electrical network condition |
| CN105259414A (en) * | 2015-11-23 | 2016-01-20 | 上海电力学院 | Power grid impedance on-line detection method on the basis of inverter |
Non-Patent Citations (1)
| Title |
|---|
| TOMI ROINILA 等: "Online Grid Impedance Measurement Using Discrete-Interval Binary Sequence Injection", 《CONTROL AND MODELING FOR POWER ELECTRONICS》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112567584A (en) * | 2018-08-20 | 2021-03-26 | 西门子股份公司 | Grid impedance estimation for power converter controller |
| CN112567584B (en) * | 2018-08-20 | 2024-03-05 | 西门子股份公司 | Grid impedance estimation for a power converter controller |
| CN109412178A (en) * | 2018-10-19 | 2019-03-01 | 湖南大学 | A kind of switching control method of high-voltage large-capacity impedance measurement equipment |
| CN109412178B (en) * | 2018-10-19 | 2021-04-30 | 湖南大学 | Switching control method of high-voltage high-capacity impedance measurement equipment |
| CN112861315A (en) * | 2021-01-11 | 2021-05-28 | 广西大学 | One-dimensional descent search method for non-convex single-target optimal power flow global solution of power system |
| CN113743806A (en) * | 2021-09-09 | 2021-12-03 | 广西康耐肯特科技有限公司 | Method for searching non-convex dual-target optimal power flow global solution of power system |
| CN113743806B (en) * | 2021-09-09 | 2024-03-19 | 广西康耐肯特科技有限公司 | Searching method for non-convex double-target optimal power flow global solution of power system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107064646B (en) | 2020-02-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107064646A (en) | A kind of electric network impedance real-time identification method for multiple-input and multiple-output asymmetric system | |
| CN107121609A (en) | A kind of electric network impedance on-line identification method and device that injection is disturbed based on PRBS | |
| CN104716668B (en) | Improve feed forward control method of the LCL type combining inverter to grid adaptability | |
| CN105158592B (en) | Three phase network impedance real-time detection method under power network asymmetrical | |
| CN107608933B (en) | Inverter generalized impedance measurement method based on secondary side disturbance | |
| CN103376363B (en) | Impedance measurement method of three-phase alternating current system | |
| CN102723722B (en) | Parallel inverter droop control based reactive power control method | |
| CN204425298U (en) | A kind of faint AC signal amplitude lock-in amplifier | |
| CN111239491B (en) | Generalized impedance real-time experimental measurement method adopting physical controller disturbance injection | |
| CN111506869A (en) | Grid-connected inverter frequency coupling impedance characteristic extraction method based on multi-sine-wave signal injection | |
| CN108155643B (en) | A kind of robust estimation method of the single-phase mains voltage parameter based on sliding mode observer | |
| CN107632200B (en) | Harmonic detection method based on discrete Fourier transform | |
| CN110323745B (en) | Analysis method for AC-DC side harmonic transmission characteristics of modular multilevel converter | |
| CN104158222A (en) | Grid-connected inverter direct power control method with voltage compensation | |
| CN107576851A (en) | System harmonic impedance measuring method based on rotatory current phasor | |
| Wei et al. | Extended Park's vector method in early inter‐turn short circuit fault detection for the stator windings of offshore wind doubly‐fed induction generators | |
| CN105914789B (en) | Inverter type distributed generation resource simplifies modeling method | |
| CN110138241A (en) | A kind of method for analyzing stability of the single-phase electricity potential source based on Jacobi's theory | |
| CN112964935B (en) | Controllable broadband power grid impedance measuring device and measuring method thereof | |
| CN207472983U (en) | A kind of electric network impedance on-line identification device based on PRBS disturbance injections | |
| CN206117616U (en) | High -power adjustable high frequency fractional order electric capacity that exponent number is greater than 1 | |
| CN106249586B (en) | A kind of phase-locked loop method for supporting the tracking of single-phase and three-phase voltage | |
| CN207148289U (en) | A kind of wireless mobile oil machine collector | |
| CN108984847B (en) | Real-time digital hybrid simulation interface method based on frequency division impedance compensation | |
| CN103592984A (en) | Method for decomposing and reconstructing current instantaneous sequence component of triangular connection current transformer |
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 | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Yanming Inventor after: Wu Huifang Inventor after: Li Hong Inventor after: Ding Xin Inventor after: Ren Xiaotong Inventor after: Han Chao Inventor after: Yang Zaixue Inventor after: Yang Liulin Inventor after: Huang Hongquan Inventor after: Mo Like Inventor before: Li Hong Inventor before: Wu Huifang Inventor before: Ding Xin Inventor before: Ren Xiaotong Inventor before: Han Chao Inventor before: Yang Zaixue Inventor before: Yang Liulin Inventor before: Huang Hongquan Inventor before: Chen Yanming Inventor before: Mo Like |
|
| CB03 | Change of inventor or designer information |