CN103346585A - Grid-connected inverter control method based on state observer - Google Patents

Grid-connected inverter control method based on state observer Download PDF

Info

Publication number
CN103346585A
CN103346585A CN2013102847484A CN201310284748A CN103346585A CN 103346585 A CN103346585 A CN 103346585A CN 2013102847484 A CN2013102847484 A CN 2013102847484A CN 201310284748 A CN201310284748 A CN 201310284748A CN 103346585 A CN103346585 A CN 103346585A
Authority
CN
China
Prior art keywords
state
state observer
observer
current
combining inverter
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
Application number
CN2013102847484A
Other languages
Chinese (zh)
Other versions
CN103346585B (en
Inventor
张兴
乔彩霞
余畅舟
徐海珍
汪杨俊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201310284748.4A priority Critical patent/CN103346585B/en
Publication of CN103346585A publication Critical patent/CN103346585A/en
Application granted granted Critical
Publication of CN103346585B publication Critical patent/CN103346585B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Inverter Devices (AREA)

Abstract

The invention discloses a grid-connected inverter control method based on a state observer. The grid-connected inverter control method comprises the steps that (1) a bridge arm current iL1 and network voltage eg are detected, and a phase angle theta is obtained through phase locking of a phase-locked loop; (2) an observation state quantity is obtained through the state observer based on an internal model; (3) a state quantity feedback signal Xfd is obtained through a state feedback matrix K; (4) grid currents of an observational network are processed, and then grid current error signals ed and eq are obtained; (5) closed-loop processing is carried out on the grid current error signals ed and eq through a PI controller, and then through coordinate inverse transformation, wave generation voltage ui1 is obtained; (6) the wave generation voltage ui1 and the state quantity feedback signal Xfd are subtracted from each other to generate a SVPWM control signal ui of an inverter bridge switching tube, output of a three-phase full bridge inverter is controlled, and therefore a distributed power generation system is controlled to be combined to the grid. Practices show that the control method obtains good dynamic and steady-state performance for a grid-connected current.

Description

A kind of combining inverter control method based on state observer
Technical field
The invention belongs to grid converter control technology field, be specifically related to a kind of combining inverter control method based on state observer.
Background technology
Along with the continuous application of renewable resources such as photovoltaic generation, wind power generation, fuel cell, distributed generation system becomes the emphasis of Recent study, is subjected to extensive concern.Combining inverter is as the core of distributed generation system, and the LCL filter construction is often adopted in its filter design.With respect to the L filter, the LCL filter can more effectively suppress the electric current high order harmonic component, and reduces the total inductance amount.Yet the LCL filter has low resistance high-order characteristic, easily causes system resonance.
For suppressing LCL combining inverter resonance, there is the scholar to propose to utilize inverter brachium pontis electric current to control the method for its networking electric current indirectly, this method has been improved the stability of system, but can not guarantee the current waveform quality that networks.Subsequently, a lot of scholars begin one's study and adopt two closed loop control methods of direct networking Current Control, and outer shroud adopts the networking current closed-loop, and interior ring adopts damping schemes such as capacitor current feedback, capacitance voltage feedback, brachium pontis current closed-loop respectively.Above-mentioned pair of close-loop control scheme suppresses system oscillation by interior ring, improves the stability of a system; Yet there is coupling in the inner and outer ring design in each scheme, influences system control performance; And these schemes need increase extra transducer, have increased the system hardware cost.
There is the scholar to propose to transform controlled device by state feedback, and then the control scheme of simplified control system design of Regulator, guaranteed the control performance of grid-connected current.Yet the realization of state feedback needs a plurality of quantity of states of detection system, has increased extra transducer.For this reason, scholars begin one's study and come the observer state feedback variable by designing traditional state observer, the additional sensors that replaces active damping link and state feedback to increase, but the method does not consider that the variation of state observer eigenmatrix parameter and state observer input deviation to the influence of observation effect, can not guarantee system control performance.
Summary of the invention
The objective of the invention is to overcome traditional state observer does not consider that state observer eigenmatrix parameter changes and the state observer input deviation to the deficiency of the influence of observation effect, a kind of combining inverter control method and state observer thereof based on state observer proposed, this method adopts a kind of state observer based on internal mold to obtain observer state amount accurately, transform controlled device by state feedback, and then carry out combining inverter networking Current Control.
To achieve these goals, the present invention is by the following technical solutions:
A kind of combining inverter control method based on state observer, the topological structure of the combining inverter that this control method relates to comprises DC source U Dc, DC side filter capacitor C Dc, three phase full bridge inverter circuit, LCL filter, DC side filter capacitor C DcBe connected in parallel on DC source U DcTwo ends, DC source U DcTwo power output ends link to each other with two inputs of three phase full bridge inverter circuit respectively, the three-phase input end of the three-phase output end of three phase full bridge inverter circuit and LCL filter is corresponding linking to each other one by one, the three-phase output end of LCL filter respectively with three phase network e a, e b, e cLink to each other; Described LCL filter is by the inverter side inductance L 1, net side filter inductance L 2With filter capacitor C 1Form; This control method comprises the steps:
Step 1, utilize the current sensor of brachium pontis electric current to detect the brachium pontis current i L1Utilize line voltage transducer detection of grid voltage e gThe horizontal lock of going forward side by side obtains phase angle theta;
Step 2, according to bridge arm voltage u iWith line voltage e g, obtain the observer state amount through the state observer based on internal mold: observation brachium pontis electric current
Figure BDA00003478918800021
The observation capacitance voltage
Figure BDA00003478918800022
And observation current on line side
Figure BDA00003478918800023
Step 3, with the described observer state amount in the step 2, obtain quantity of state feedback signal X through state feedback matrix K Fd, specific as follows:
1. try to achieve state feedback matrix K, K=[k by expectation POLE PLACEMENT USING method 1, k 2, k 3], k 1, k 2, k 3Be the state feedback coefficient;
2. with the described observer state amount in the step 2, namely observe the brachium pontis electric current
Figure BDA00003478918800024
The observation capacitance voltage
Figure BDA00003478918800025
The observation current on line side Multiply by the state feedback coefficient k respectively 1, k 2, k 3, can obtain described quantity of state feedback signal X fd = [ k 1 × i ^ L 1 ; k 2 × u ^ C ; k 3 × i ^ L 2 ] ;
Step 4, will observe current on line side Obtain current on line side real component i through the abc-dq coordinate transform D_fdbWith current on line side idle component i Q_fdb, with current on line side is given function signal i arranged respectively D_fref, given no function signal i Q_refDiffer from, obtain the current on line side error signal e d, e q, that is:
e d=i d_fref-i d_fdb
e q=i q_ref-i q_fdb
Step 5, by the current on line side error signal e of PI controller to obtaining in the step 4 d, e qCarry out closed-loop process, then through dq-abc coordinate inverse transformation, obtain sending out wave voltage u I1
Step 6, with described wave voltage u I1With described quantity of state feedback signal X FdStack can generate the SVPWM control signal u of three phase full bridge inverter circuit i
Further, the state observer based on internal mold in the described step 2, realize by following steps:
1. select the brachium pontis current i L1, capacitance voltage u C, current on line side i L2As system state amount, with bridge arm voltage u iWith line voltage e gAs input variable, obtain combining inverter system discrete state space equation by following formula:
x ( k + 1 ) = A d x ( k ) + B d u ( k ) y ( k ) = C d 1 x ( k ) + D d u ( k ) - - - ( 1 )
In the formula (1), A dBe combining inverter system features matrix, B dBe combining inverter system input matrix, C D1Be combining inverter system output matrix, D dBe the direct transmission matrix of combining inverter system;
X (k) represents described system state amount: the brachium pontis current i L1, capacitance voltage u CAnd current on line side i L2
X (k+1) claps the system state amount of delaying time through one;
U (k) is combining inverter system input variable;
Y (k) is combining inverter system output variable;
2. according to described combining inverter system discrete state space equation, consider that respectively state observer eigenmatrix parameter changes
Figure BDA00003478918800032
And state observer input deviation Δ u (k), by the principle of duality, utilize following formula to obtain state observer state space equation based on internal mold:
x ^ ( k + 1 ) = A ^ d x ^ ( k ) + B ^ d u ′ ( k ) + G ( 1 + φ e - 1 ( z ) ) [ y ( k ) - y ^ ( k ) ] y ^ ( k ) = C ^ d 1 x ^ ( k ) + D ^ d u ( k ) - - - ( 2 )
In the formula (2),
Figure BDA00003478918800034
Be the state observer eigenmatrix,
Figure BDA00003478918800035
Be the state observer input matrix,
Figure BDA00003478918800036
Be the state observer output matrix,
Figure BDA00003478918800037
Be the direct transmission matrix of state observer;
Figure BDA00003478918800038
Represent described observer state amount: observation brachium pontis electric current The observation capacitance voltage
Figure BDA000034789188000310
And observation current on line side
Figure BDA000034789188000312
For clapping the observer state amount of delaying time through one;
U ' is the state observer input variable (k), and G is the feedback matrix of state observer;
Figure BDA000034789188000313
Be the state observer output variable, y (k) is combining inverter system output variable,
Figure BDA000034789188000314
Be the state observer output error;
Figure BDA000034789188000315
Be the internal mold item.
3. described state observer eigenmatrix When parameter changes, the observation error Δ x that obtains by described formula (1) and formula (2) 1Be expressed from the next for:
Figure BDA00003478918800041
In the formula (3),
Figure BDA00003478918800042
Be state observer eigenmatrix variable quantity,
Figure BDA00003478918800043
For considering the state observer eigenmatrix
Figure BDA00003478918800044
The internal mold item that arranges when parameter changes is according to described state observer output error
Figure BDA00003478918800045
Be set to single order integral element or ratio resonance link;
(a) part in the formula (3) reaches 0 value by expectation POLE PLACEMENT USING method, and formula (3) is by (b) part internal mold
Figure BDA00003478918800046
Setting, make (b) part and (c) part the two and reach 0 value, final described observation error Δ x 1 Reach 0 value.
4. described state observer input variable u ' is (k) when existing deviation, the observation error Δ x that through type (1) and formula (2) obtain 2Be expressed from the next for:
Figure BDA00003478918800047
In the formula (4), Δ u (k) is the variable quantity of state observer input variable u (k),
Figure BDA00003478918800048
For the internal mold item of considering that state observer input u ' (k) arranges when having deviation, according to described state observer output error
Figure BDA00003478918800049
Be set to single order integral element or ratio resonance link;
(d) part in the formula (4) reaches 0 value by expectation POLE PLACEMENT USING method, and formula (4) is by (e) part internal mold
Figure BDA000034789188000410
Setting, make (e) part and (f) part the two and reach 0 value, final described observation error Δ x 2 Reach 0 value.
The present invention considered that state observer eigenmatrix parameter changes and the state observer input deviation to the influence of observation effect, obtain observer state amount accurately by the state observer based on internal mold, transform controlled device through state feedback, adopt the observation current on line side to carry out the combining inverter current loop control, obtained good dynamic and steady-state behaviour.
Description of drawings
Fig. 1 is the combining inverter control structure figure that the present invention is based on state observer.
Fig. 2 is traditional state observer structure chart.
Fig. 3 is based on the state observer structure chart of internal mold.
Fig. 4 a, 4b are respectively traditional state observers and dynamic observe effect contrast simulation figure based on the state observer of internal mold at rest frame brachium pontis electric current.
Fig. 5 a, 5b are respectively traditional state observers and dynamic observe effect contrast simulation figure based on the state observer of internal mold at rotating coordinate system brachium pontis electric current.
Fig. 6 a, 6b be respectively traditional state observer and based on the state observer of internal mold at rest frame brachium pontis electric current stable state observation effect contrast simulation figure.
Fig. 7 a, 7b be respectively traditional state observer and based on the state observer of internal mold at rest frame brachium pontis electric current dynamic contrast experimental waveform.
Fig. 8 a, 8b be respectively traditional state observer and based on the state observer of internal mold at rest frame brachium pontis electric current stable state contrast experiment waveform.
Fig. 9 a, 9b be respectively traditional state observer and based on the state observer of internal mold at rotating coordinate system brachium pontis active current contrast experiment waveform.
Figure 10 a, 10b be respectively traditional state observer and based on the state observer of internal mold at rest frame capacitance voltage dynamic contrast experimental waveform.
Figure 11 a, 11b be respectively traditional state observer and based on the state observer of internal mold at the current on line side dynamic experiment waveform based on the observer state feedback.
Embodiment
In order more specifically to describe the present invention, below in conjunction with the drawings and the specific embodiments combining inverter control method and the state observer thereof that the present invention is based on state observer is elaborated.
In the present embodiment, as shown in Figure 1, three-phase LCL combining inverter converts the 260V direct current to the 220V alternating current, and through importing electrical network behind the LCL filter filtering.
Described LCL filter parameter: brachium pontis filter inductance L 1=1mH, brachium pontis filter inductance dead resistance r 1=0.001 Ω, filter capacitor C 1=20 μ F, net side filter inductance L 2=0.5mH, net side filter inductance dead resistance r 2=0.001 Ω, electric network impedance L gExcursion 0~1.5mH.
A kind of combining inverter control method based on state observer may further comprise the steps:
1, utilizes the current sensor H of brachium pontis electric current L1a, H L1bDetect the brachium pontis current i L1Utilize line voltage transducer H Ea, H EbDetection of grid voltage e a, e b, e cThe horizontal lock of going forward side by side obtains phase angle theta;
2, according to bridge arm voltage u iWith line voltage e gObtain following observer state amount through the state observer based on internal mold: observation brachium pontis electric current
Figure BDA00003478918800051
The observation capacitance voltage
Figure BDA00003478918800052
And observation current on line side
Figure BDA00003478918800053
Present embodiment is as follows:
Choose state variable: the brachium pontis current i L1, capacitance voltage u C, current on line side i L2, be designated as x=[i L1, u C, i L2] TChoose input variable: inverter brachium pontis output voltage u iWith electrical network electromotive force e g, be designated as u=[u i, e g] T,
Set up LCL combining inverter system state space model:
x · = Ax + Bu y = Cx + Du
In the following formula, A is combining inverter system features matrix, A = - r 1 + R p L 1 - 1 L 1 R p L 1 1 C 0 - 1 C R p L 2 ′ 1 L 2 ′ - r 2 + R P L 2 ′ ;
Wherein, R pBe filter capacitor C 1The damping resistance of last series connection;
L 2' be net side filter inductance L 2With electric network impedance equivalent inductance L gSum;
B is combining inverter system input matrix, B = . B 1 . B 2 . = 1 / L 1 0 0 0 0 - 1 / L 2 ′ ;
C is combining inverter system output matrix, when selecting the brachium pontis current i L1During as output, C 1=[1,0,0]; When selecting current on line side i L2During as output, C 2=[0,0,1];
D is the direct transmission matrix of combining inverter system, D=[0,0].
For the ease of digital control, select sample frequency f s=10kHz, utilize the discretization method of zero-order holder: C D1=C 1, D d=D with continuous state equation discretization, obtains LCL combining inverter discrete state space equation:
x ( k + 1 ) = A d x ( k ) + B d u ( k ) y ( k ) = C d 1 x ( k ) + D d u ( k )
In the following formula, C D1=[1,0,0] namely selects the brachium pontis electric current as output;
X (k) is discretization quantity of state x (k)=[i L1(k), u C(k), i L2(k)] T, wherein: i L1(k) be discrete domain brachium pontis electric current, u C(k) be the discrete domain capacitance voltage, i L2(k) be the discrete domain current on line side;
U (k) is the discrete domain input variable, u (k)=[u i(k), e g(k)] T, wherein: u i(k) be discrete domain inverter brachium pontis output voltage, e g(k) be discrete domain electrical network electromotive force.
According to LCL combining inverter discrete state spatial model, on traditional state observer basis of (as shown in Figure 2), consider that state observer eigenmatrix parameter changes and the variation of observer input, obtain the state observer based on internal mold, its state space is expressed as:
x ^ ( k + 1 ) = A ^ d x ^ ( k ) + B ^ d u ′ ( k ) + G ( 1 + φ e - 1 ( z ) ) [ y ( k ) - y ^ ( k ) ] y ^ ( k ) = C ^ d 1 x ^ ( k ) + D ^ d u ( k )
Wherein,
Figure BDA00003478918800072
Be the state observer eigenmatrix,
Figure BDA00003478918800073
Be the state observer input matrix,
Figure BDA00003478918800074
Be the state observer output matrix, Be the direct transmission matrix of state observer,
Figure BDA00003478918800076
Be state observer eigenmatrix variable quantity;
Figure BDA00003478918800077
Represent described observer state amount: observation brachium pontis electric current
Figure BDA00003478918800078
The observation capacitance voltage And observation current on line side
Figure BDA000034789188000710
U (k) is the state observer input variable, and Δ u (k) is the variable quantity of state observer input variable u (k);
G is the feedback matrix of state observer;
Be the state observer output variable, y (k) is the combining inverter output variable;
Figure BDA000034789188000712
Expression internal mold item.
With reference to the accompanying drawings 3, described state observer based on internal mold, on the basis of LCL discrete state space equation, with combining inverter output variable y (k), i.e. combining inverter brachium pontis side filter inductance current i L1With the state observer output variable
Figure BDA000034789188000713
Namely observe the brachium pontis electric current Differ from, obtain output error e (k); Because in different coordinate systems, the expression-form difference of output error for ease of analyzing, is considered the general type of output error, is designated as E;
According to described output error E, obtain its proper polynomial:
φ e(z)=det(zI-E)
In the formula, φ e(z) expression output error E proper polynomial, z is the discrete transform operator;
Described output error E is thought of as the input of internal mold, represents φ by following formula e(z) unstable limit multinomial:
φ(z)=z m+a m-1z m-1+a m-2z m-2+…+a 1z+a 0
In the formula, in the formula, φ (z) is the monic polynomial of z, a iThe expression multinomial coefficient, described a i∈ [0, m-1], m are discrete transform factor order;
Internal mold φ -1(z) realize by following state equation:
x c ( k + 1 ) = G c x c ( k ) + b c e ( k ) y c ( k ) = x c ( k )
In the formula, G cFor internal mold is realized the formulation character matrix, G c = 0 1 0 0 . . . 0 0 0 1 0 . . . 0 . . . . . . 0 . . . 0 1 - a 0 - a 1 . . . - a m - 1 ,
b cFor internal mold is realized input matrix, b c=[0,0 ... 0,1] T
x c(k) be the internal mold quantity of state, x c(k+1) be that e (k) is actual output error, is expressed as through the internal mold state of a bat time-delay e ( k ) = y ( k ) - y ^ ( k ) ;
In the three-phase rotating coordinate system, error and the quantity of state of output variable are constant, internal mold Usually can be set to the single order integral controller; In rest frame, error and the quantity of state of output variable are sinusoidal quantity, and internal mold can be set to the resonance adjuster.
As Fig. 3, with the output item x of output error e (k) with the process internal mold c(k) addition, G can obtain observer state amount accurately through the state observation matrix: observation brachium pontis electric current The observation capacitance voltage
Figure BDA00003478918800086
The observation current on line side
Figure BDA00003478918800087
Among Fig. 3, I gBe net side active current, Zoh is zero-order holder, 1/ (sL 1+ r 1) be brachium pontis inductance transfer function, 1/ (sC 1) be the filter capacitor transfer function, 1/ (sL 2'+r 2) be net side filter inductance transfer function.
3, according to the described observer state amount of step 2, obtain quantity of state feedback signal X through state feedback matrix K FdIntroduce described state feedback, the limit of system is reconfigured, make the pole distribution of system in the unit circle of z territory, thereby increase the damping of control system.
In the present embodiment, separate to obtain according to the described state observer based on internal mold of separation principle state feedback and step 2.
Discrete state space expression formula based on state feedback is
x(k+1)=(A d-B dK)x(k)+B du(k)
Comprise damping resistance R according to capacitive branch pThe system features matrix A d, A dIn get filter capacitor series resistance R p=0.5 Ω obtains the system features multinomial that comprises passive damping; According to the system features matrix A through state feedback d-B dK utilizes the respective items coefficient to equate, tries to achieve state feedback matrix K, namely
det(zI-A d)=det(zI-(A d-B dK))
Wherein, z is the discrete transform operator, and I is 3 * 3 unit matrixs, det () representation feature multinomial, and K is state feedback matrix, is designated as K=[k 1, k 2, k 3], k 1, k 2, k 3Be called the state feedback coefficient;
The state feedback coefficient of trying to achieve: k 1=2.8108609696657359180517147320431,
k 2=-0.42233855315199843186755694448278,
k 3=-2.8075596495029625451005704996724;
With the described observer state amount of step 2: observation brachium pontis electric current
Figure BDA00003478918800091
The observation capacitance voltage
Figure BDA00003478918800092
The observation current on line side
Figure BDA00003478918800093
Multiply by the state feedback coefficient k respectively 1, k 2, k 3, obtain the quantity of state feedback signal X fd = [ k 1 × i ^ L 1 ; k 2 × u ^ C ; k 3 × i ^ L 2 ] .
4, will observe current on line side
Figure BDA00003478918800095
Obtain current on line side real component i through coordinate transform D_fdbWith current on line side idle component i Q_fdb, with current on line side is given function signal i arranged respectively D_fref, given no function signal i Q_refDiffer from, obtain the current on line side error signal e d, e q
5, pass through the PI controller to the current on line side error signal e d, e qCarry out closed-loop process, wherein, the PI controller is got k p=15, k i=1800, through the coordinate inverse transformation, obtain sending out wave voltage u I1
6, with u I1With state feedback signal X FdIt is the SVPWM control signal u of three phase full bridge inverter circuit that stack generates the converter bridge switching parts pipe i, and then control combining inverter networking electric current.
Fig. 4 a~Fig. 6 b be traditional state observer and based on the state observer of internal mold when the net side impedance variation, observation effect dynamically, stable state contrast simulation waveform.Fig. 4 a, 4b rest frame brachium pontis electric current dynamic observe effect contrast simulation figure, and number in the figure ' 1 ' represents brachium pontis electric current actual measurement waveform, and label ' 2 ' represents brachium pontis electric current observation waveform; Fig. 5 a, 5b rotating coordinate system brachium pontis electric current dynamic observe effect contrast simulation figure, and number in the figure ' 1 ' represents brachium pontis active current actual measurement waveform, and label ' 2 ' represents brachium pontis active current observation waveform; Label ' 3 ' represents brachium pontis reactive current actual measurement waveform, and label ' 4 ' represents brachium pontis reactive current observation waveform.
Fig. 7 a~Figure 10 b be traditional state observer and based on the internal mold state observer when the net side impedance variation, observation effect dynamically, stable state contrast experiment's waveform.From analogous diagram and lab diagram as can be seen, based on the state observer of internal mold, still be that steady state effect all is better than traditional state observer from dynamic effect no matter, and the observer state amount is used for state feedback, and then being used for grid-connected current control, the current on line side dynamic waveform is good.
The present invention can apply in the multiple combining inverter control structure, reduces the control system cost, improves the control system reliability simultaneously.This method has versatility, is not limited to hardware configuration and control structure.

Claims (7)

1. combining inverter control method based on state observer, the topological structure of the combining inverter that this control method relates to comprises DC source U Dc, DC side filter capacitor C Dc, three phase full bridge inverter circuit, LCL filter, described DC side filter capacitor C DcBe connected in parallel on DC source U DcTwo ends, DC source U DcTwo power output ends link to each other with two inputs of three phase full bridge inverter circuit respectively, the three-phase input end of the three-phase output end of three phase full bridge inverter circuit and LCL filter is corresponding linking to each other one by one, the three-phase output end of LCL filter respectively with three phase network e a, e b, e cLink to each other; Described LCL filter is by the inverter side inductance L 1, net side filter inductance L 2With filter capacitor C 1Form; It is characterized in that comprising the steps:
Step 1, utilize the current sensor of brachium pontis electric current to detect the brachium pontis current i L1Utilize line voltage transducer detection of grid voltage e gThe horizontal lock of going forward side by side obtains phase angle theta;
Step 2, according to bridge arm voltage u iWith line voltage e g, obtain the observer state amount through the state observer based on internal mold: observation brachium pontis electric current
Figure FDA00003478918700011
The observation capacitance voltage
Figure FDA00003478918700012
And observation current on line side
Figure FDA00003478918700013
Step 3, with the described observer state amount of step 2, obtain quantity of state feedback signal X through state feedback matrix K Fd
Step 4, will observe current on line side
Figure FDA00003478918700014
Obtain current on line side real component i through the abc-dq coordinate transform D_fdbWith current on line side idle component i Q_fdb, with current on line side is given function signal i arranged respectively D_fref, given no function signal i Q_refDiffer from, obtain the current on line side error signal e d, e q, that is:
e d=i d_fref-i d_fdb
e q=i q_ref-i q_fdb
Step 5, by the current on line side error signal e of PI controller to obtaining in the step 4 d, e qCarry out closed-loop process, abc-dq coordinate inverse transformation obtains sending out wave voltage u then I1
Step 6, with described wave voltage u I1With described quantity of state feedback signal X FdStack can generate the SVPWM control signal u of three phase full bridge inverter circuit i
2. a kind of combining inverter control method based on state observer according to claim 1 is characterized in that obtaining quantity of state feedback signal X through state feedback matrix K described in the step 3 FdStep as follows:
1. try to achieve state feedback matrix K, K=[k by expectation POLE PLACEMENT USING method 1, k 2, k 3], k 1, k 2, k 3Be the state feedback coefficient;
2. with the observer state amount described in the step 2, namely observe the brachium pontis electric current
Figure FDA00003478918700015
The observation capacitance voltage The observation current on line side
Figure FDA00003478918700021
Multiply by the state feedback coefficient k respectively 1, k 2, k 3, can obtain the quantity of state feedback signal
3. a kind of combining inverter control method based on state observer according to claim 1 is characterized in that the described state observer based on internal mold of step 2 realizes by following steps:
1. select the brachium pontis current i L1, capacitance voltage u C, current on line side i L2As system state amount, with bridge arm voltage u iWith line voltage e gAs input variable, obtain combining inverter system discrete state space equation by following formula:
Figure RE-FDA00003597886500013
In the formula (1), A dBe combining inverter system features matrix, B dBe combining inverter system input matrix, C D1Be combining inverter system output matrix, D dBe the direct transmission matrix of combining inverter system;
X (k) represents described system state amount: the brachium pontis current i L1, capacitance voltage u CAnd current on line side i L2
X (k+1) claps the system state amount of delaying time through one;
U (k) is combining inverter system input variable;
Y (k) is combining inverter system output variable;
2. according to described combining inverter system discrete state space equation, consider that respectively state observer eigenmatrix parameter changes
Figure RE-FDA00003597886500014
And state observer input deviation Δ u (k), by the principle of duality, utilize following formula to obtain state observer state space equation based on internal mold:
Figure RE-FDA00003597886500015
In the formula (2),
Figure RE-FDA00003597886500016
Be the state observer eigenmatrix,
Figure RE-FDA00003597886500017
Be the state observer input matrix,
Figure RE-FDA00003597886500018
Be the state observer output matrix,
Figure RE-FDA00003597886500019
Be the direct transmission matrix of state observer;
Figure RE-FDA000035978865000110
Represent described observer state amount: observation brachium pontis electric current
Figure RE-FDA000035978865000111
The observation capacitance voltage
Figure RE-FDA000035978865000112
And observation current on line side
Figure RE-FDA000035978865000113
Figure RE-FDA000035978865000114
For clapping the observer state amount of delaying time through one;
U ' is the state observer input variable (k), and G is the feedback matrix of state observer;
Figure RE-FDA000035978865000115
Be the state observer output variable, y (k) is combining inverter system output variable,
Figure RE-FDA000035978865000116
Be the state observer output error;
Be the internal mold item.
4. a kind of combining inverter control method based on state observer according to claim 3 is characterized in that: described state observer eigenmatrix
Figure FDA00003478918700031
When parameter changes, the observation error Δ x that obtains by described formula (1) and formula (2) 1Be expressed from the next for:
Figure FDA00003478918700032
In the formula (3),
Figure FDA00003478918700033
Be state observer eigenmatrix variable quantity, For considering the state observer eigenmatrix
Figure FDA00003478918700035
The internal mold item that arranges when parameter changes,
(a) part in the formula (3) reaches 0 value by expectation POLE PLACEMENT USING method, and formula (3) is by (b) part internal mold
Figure FDA00003478918700036
Setting, make (b) part and (c) part the two and reach 0 value, final described observation error Δ x 1Reach 0 value.
5. a kind of combining inverter control method and state observer thereof based on state observer according to claim 3, it is characterized in that: described state observer input variable u ' is (k) when existing deviation, the observation error Δ x that through type (1) and formula (2) obtain 2Be expressed from the next for:
In the formula (4), Δ u (k) is the variable quantity of state observer input variable u (k),
Figure RE-FDA00003597886500028
Be the internal mold item of considering that state observer input u ' (k) arranges when having deviation,
(d) part in the formula (4) reaches 0 value by expectation POLE PLACEMENT USING method, and formula (4) is by (e) part internal mold Setting, make (e) part and (f) part the two and reach 0 value, final described observation error Δ x 2Reach 0 value.
6. a kind of combining inverter control method based on state observer according to claim 4 is characterized in that: the internal mold item in the described formula (3)
Figure FDA000034789187000310
According to described state observer output error
Figure FDA000034789187000311
Be set to single order integral element or ratio resonance link.
7. a kind of combining inverter control method based on state observer according to claim 5 is characterized in that: the internal mold item in the described formula (4)
Figure FDA000034789187000312
According to described state observer output error
Figure FDA000034789187000313
Be set to single order integral element or ratio resonance link.
CN201310284748.4A 2013-07-08 2013-07-08 Grid-connected inverter control method based on state observer Expired - Fee Related CN103346585B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310284748.4A CN103346585B (en) 2013-07-08 2013-07-08 Grid-connected inverter control method based on state observer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310284748.4A CN103346585B (en) 2013-07-08 2013-07-08 Grid-connected inverter control method based on state observer

Publications (2)

Publication Number Publication Date
CN103346585A true CN103346585A (en) 2013-10-09
CN103346585B CN103346585B (en) 2015-05-27

Family

ID=49281368

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310284748.4A Expired - Fee Related CN103346585B (en) 2013-07-08 2013-07-08 Grid-connected inverter control method based on state observer

Country Status (1)

Country Link
CN (1) CN103346585B (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104022511A (en) * 2014-05-28 2014-09-03 国家电网公司 Dimension reduction observer used for LCL grid-connected inversion damping control
CN104158222A (en) * 2014-08-29 2014-11-19 东南大学 Grid-connected inverter direct power control method with voltage compensation
CN104898071A (en) * 2015-06-12 2015-09-09 东南大学 Fault diagnosis method for modularized multi-level inverter based on state observation
CN104950201A (en) * 2015-06-12 2015-09-30 东南大学 Fault diagnosis method of modular multilevel inverter based on state observation
CN106301049A (en) * 2015-05-28 2017-01-04 中国科学院沈阳自动化研究所 Current source inverter mixing H2/ H∞optimum guaranteed cost control method
CN106374765A (en) * 2016-11-04 2017-02-01 龙岩学院 Inverter control system of sensor without back electromotive force and control method thereof
CN107134936A (en) * 2017-05-22 2017-09-05 天津威瀚电气股份有限公司 A kind of inverter active damping control method based on overall-finished housing
CN107181415A (en) * 2017-06-22 2017-09-19 同济大学 A kind of solid-state transformer model predictive control method
CN107749641A (en) * 2017-11-07 2018-03-02 国网重庆市电力公司电力科学研究院 Consider that the combining inverter grid-connected current of control delay improves repetitive control
CN107769251A (en) * 2017-11-07 2018-03-06 国网重庆市电力公司电力科学研究院 Parallel inverter current control method based on modified state observer
CN108155643A (en) * 2017-12-22 2018-06-12 上海交通大学 A kind of robust estimation method of the single-phase mains voltage parameter based on sliding mode observer
CN109193622A (en) * 2018-09-06 2019-01-11 上海电力学院 It is a kind of for exchanging the Parameters design of micro-capacitance sensor power quality controlling device
CN109193793A (en) * 2018-09-13 2019-01-11 深圳市禾望电气股份有限公司 A kind of current transformer exempts from the grid-connection control system and method for voltage detecting
CN110311416A (en) * 2019-06-14 2019-10-08 合肥工业大学 The adaptive control method of grid-connected inverter of bandwidth of phase lock loop based on state feedback
CN111435141A (en) * 2019-01-15 2020-07-21 中国电力科学研究院有限公司 Power grid impedance measuring device and method
CN111740675A (en) * 2020-07-02 2020-10-02 合肥工业大学 Two-degree-of-freedom control method for discrete domain current loop high robustness of permanent magnet synchronous motor
CN112491275A (en) * 2020-11-20 2021-03-12 西安热工研究院有限公司 DC-DC conversion circuit control method based on active disturbance rejection observer
CN112564133A (en) * 2020-12-14 2021-03-26 广西电网有限责任公司北海供电局 Intelligent power generation control method based on deep learning full-state optimal feedback and application
CN113517710A (en) * 2021-06-16 2021-10-19 武汉理工大学 Electric vehicle V2G inversion control method and system based on state observer
CN114142774A (en) * 2021-12-07 2022-03-04 杭州电子科技大学 PMSM phase current reconstruction method based on sine curve fitting observer
US11757375B2 (en) 2021-12-27 2023-09-12 Industrial Technology Research Institute DC-AC inverter system using state observer and control method thereof

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6198248B1 (en) * 1998-11-04 2001-03-06 Alstom Enterprise Sa Method of controlling a rotary electrical machine, a servo-control system for implementing the method, and a rotary machine fitted with such a system
CN101557174A (en) * 2009-04-22 2009-10-14 华中科技大学 A state-tracked digitally-controlled inverter power supply

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6198248B1 (en) * 1998-11-04 2001-03-06 Alstom Enterprise Sa Method of controlling a rotary electrical machine, a servo-control system for implementing the method, and a rotary machine fitted with such a system
CN101557174A (en) * 2009-04-22 2009-10-14 华中科技大学 A state-tracked digitally-controlled inverter power supply

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
侯朝勇等: "基于离散状态空间模型的LCL 滤波并网变换器控制策略", 《中国电机工程学报》 *
孔雪娟等: "基于内模原理的三相电压源型逆变电源的波形控制技术", 《中国电机工程学报》 *
张凯等: "基于状态反馈与重复控制的逆变器控制技术", 《中国电机工程学报》 *
杨淑英等: "LCL 滤波电压源并网逆变器多环控制策略设计", 《电力系统自动化》 *

Cited By (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104022511A (en) * 2014-05-28 2014-09-03 国家电网公司 Dimension reduction observer used for LCL grid-connected inversion damping control
CN104158222A (en) * 2014-08-29 2014-11-19 东南大学 Grid-connected inverter direct power control method with voltage compensation
CN106301049A (en) * 2015-05-28 2017-01-04 中国科学院沈阳自动化研究所 Current source inverter mixing H2/ H∞optimum guaranteed cost control method
CN106301049B (en) * 2015-05-28 2018-09-04 中国科学院沈阳自动化研究所 The optimal guaranteed cost control methods of current source inverter mixing H2/H ∞
CN104950201B (en) * 2015-06-12 2017-11-28 东南大学 A kind of modular multilevel inverter method for diagnosing faults based on state observation
CN104898071A (en) * 2015-06-12 2015-09-09 东南大学 Fault diagnosis method for modularized multi-level inverter based on state observation
CN104950201A (en) * 2015-06-12 2015-09-30 东南大学 Fault diagnosis method of modular multilevel inverter based on state observation
CN104898071B (en) * 2015-06-12 2018-02-02 东南大学 Modular multilevel inverter method for diagnosing faults based on state observation
CN106374765B (en) * 2016-11-04 2019-05-17 龙岩学院 A kind of control method of the inverter control system of no inverse electromotive force sensor
CN106374765A (en) * 2016-11-04 2017-02-01 龙岩学院 Inverter control system of sensor without back electromotive force and control method thereof
CN107134936A (en) * 2017-05-22 2017-09-05 天津威瀚电气股份有限公司 A kind of inverter active damping control method based on overall-finished housing
CN107181415A (en) * 2017-06-22 2017-09-19 同济大学 A kind of solid-state transformer model predictive control method
CN107749641A (en) * 2017-11-07 2018-03-02 国网重庆市电力公司电力科学研究院 Consider that the combining inverter grid-connected current of control delay improves repetitive control
CN107769251A (en) * 2017-11-07 2018-03-06 国网重庆市电力公司电力科学研究院 Parallel inverter current control method based on modified state observer
CN108155643B (en) * 2017-12-22 2019-12-03 上海交通大学 A kind of robust estimation method of the single-phase mains voltage parameter based on sliding mode observer
CN108155643A (en) * 2017-12-22 2018-06-12 上海交通大学 A kind of robust estimation method of the single-phase mains voltage parameter based on sliding mode observer
CN109193622A (en) * 2018-09-06 2019-01-11 上海电力学院 It is a kind of for exchanging the Parameters design of micro-capacitance sensor power quality controlling device
CN109193793B (en) * 2018-09-13 2021-10-19 深圳市禾望电气股份有限公司 Converter voltage detection-free grid-connected control system and method
CN109193793A (en) * 2018-09-13 2019-01-11 深圳市禾望电气股份有限公司 A kind of current transformer exempts from the grid-connection control system and method for voltage detecting
CN111435141A (en) * 2019-01-15 2020-07-21 中国电力科学研究院有限公司 Power grid impedance measuring device and method
CN111435141B (en) * 2019-01-15 2023-12-15 中国电力科学研究院有限公司 Power grid impedance measurement device and method
CN110311416B (en) * 2019-06-14 2020-09-01 合肥工业大学 Phase-locked loop bandwidth self-adaptive grid-connected inverter control method based on state feedback
CN110311416A (en) * 2019-06-14 2019-10-08 合肥工业大学 The adaptive control method of grid-connected inverter of bandwidth of phase lock loop based on state feedback
CN111740675A (en) * 2020-07-02 2020-10-02 合肥工业大学 Two-degree-of-freedom control method for discrete domain current loop high robustness of permanent magnet synchronous motor
CN111740675B (en) * 2020-07-02 2022-07-19 合肥工业大学 Two-degree-of-freedom control method for discrete domain current loop high robustness of permanent magnet synchronous motor
CN112491275A (en) * 2020-11-20 2021-03-12 西安热工研究院有限公司 DC-DC conversion circuit control method based on active disturbance rejection observer
CN112564133A (en) * 2020-12-14 2021-03-26 广西电网有限责任公司北海供电局 Intelligent power generation control method based on deep learning full-state optimal feedback and application
CN112564133B (en) * 2020-12-14 2023-03-21 广西电网有限责任公司北海供电局 Intelligent power generation control method based on deep learning full-state optimal feedback and application
CN113517710A (en) * 2021-06-16 2021-10-19 武汉理工大学 Electric vehicle V2G inversion control method and system based on state observer
CN113517710B (en) * 2021-06-16 2022-05-03 武汉理工大学 Electric vehicle V2G inversion control method and system based on state observer
CN114142774A (en) * 2021-12-07 2022-03-04 杭州电子科技大学 PMSM phase current reconstruction method based on sine curve fitting observer
US11757375B2 (en) 2021-12-27 2023-09-12 Industrial Technology Research Institute DC-AC inverter system using state observer and control method thereof

Also Published As

Publication number Publication date
CN103346585B (en) 2015-05-27

Similar Documents

Publication Publication Date Title
CN103346585B (en) Grid-connected inverter control method based on state observer
Zhang et al. Harmonic transfer-function-based impedance modeling of a three-phase VSC for asymmetric AC grid stability analysis
CN103475029B (en) Three-phase LCL type grid-connected inverter control system and method based on pole assignment
CN103326386B (en) Capacitor-voltage-based grid-connected inverter active damping method
Zhang et al. Direct grid‐side current model predictive control for grid‐connected inverter with LCL filter
CN104135003B (en) APF control method based on active disturbance rejection and repetitive control
CN111555339B (en) Converter grid-connected general sequence impedance model for stability analysis and modeling method
CN110020510B (en) A kind of method for analyzing stability of the EMU net side Pulse rectifier based on Floquet theory
CN105790306B (en) LCL type control method of grid-connected inverter based on inverter current single-sensor
Shah et al. Sequence domain transfer matrix model of three-phase voltage source converters
CN112803461B (en) Active disturbance rejection control method of battery energy storage converter based on extended state observer
Lu et al. Resonance propagation of parallel-operated DC-AC converters with LCL filters
CN107332244A (en) The level SAPF of three-phase four-wire system NPC types three nonlinear control method
Givaki et al. Stability analysis of large wind farms connected to weak AC networks incorporating PLL dynamics
Mahamat et al. Optimized design of an LCL filter for grid connected photovoltaic system and analysis of the impact of neighbors' consumption on the system
CN113595430B (en) Three-loop controller of LCL type grid-connected inverter and parameter design method thereof
CN111064225A (en) Grid-connected inverter current loop control device based on linear active disturbance rejection control
Gonzatti et al. Virtual impedance-based grid synchronization for converters connected through long cables
CN104377721B (en) VSC-HVDC optimal control method during a kind of unbalanced source voltage
CN105305392A (en) Symmetrical component method for short circuit calculation of voltage-controlled type IIDG included power distribution network
Yu et al. Impedance modeling and stability analysis of LCL-type grid-connected inverters with different current sampling schemes
Tang et al. Grid-current-sensorless control of grid-forming inverter with LCL filter
Chou et al. Reflection coefficient stability criterion for multi-bus multi-VSC power systems
CN105915094A (en) LCL type grid-connected inverter control method based on power grid current single sensor
D'Arco et al. Eigenvalue-based analysis of small-signal dynamics and stability in DC grids

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150527

Termination date: 20190708

CF01 Termination of patent right due to non-payment of annual fee