A Novel Secondary Control for Microgrid Based on Synergetic Control of Multi-Agent System

Zhiwen Yu 1,*, Qian Ai 1, Jinxia Gong 2 and Longjian Piao 1 1 School of Electronic Information and Electrical Engineering, Dongchuan Road, Shanghai Jiaotong University, Shanghai 200240, China; aiqian@sjtu.edu.cn (Q.A.); longjianpiao@gmail.com (L.P.) 2 Department of Electrical Power Engineering, Changyang Road, Shanghai University of Electrical Power, Shanghai 200240 China; gongruth@foxmail.com * Correspondence: yuzhiwen@sjtu.edu.cn; Tel.: +86-21-3420-4584


Introduction
Distributed generations (DGs), especially renewable energy sources (RES), are drawing increasing attention from the world due to the lower energy consumption, high efficiency, economy and environmental benefits [1,2].Microgrids (MG) consisting of a cluster of DGs, loads and storage is a very useful way to consume the small capacity DGs, and plays an increasingly important role in power systems because of the high reliability of power supply [3][4][5][6][7][8].However, the operation of MG is very complex, including the grid-connected and islanded modes.Proper controls are necessary to keep operations stable and economically efficient.
In power systems, the frequency and voltage are needed to keep in an expectant range.For this purpose, three levels of control are adopted [9][10][11].The primary level is the local automation control to adjust the generated power for dynamic loads.In islanded MG, the droop control is usually used in this level.The secondary layer is centralized or decentralized control which restores the frequency and voltage, and activates in 30s to 15min.Tertiary control considers the economic concerns in the optimal operation of the grid, and dispatches the generators according to the load prediction.
In normal operation, the MG connects to the main grid, and most system dynamics are dictated by the main grid due to the relatively small size of DGs.In the event of disturbances, the MG disconnects from the main grid and switches to the islanded operation, and the system dynamics are dictated by the control of MG.In islanded mode, the primary control is applied to maintain the voltage and frequency stability [5][6][7][8][9][10].All DGs are responsible for maintaining the system frequency and voltage inductance and primary source.The model details of each part are described in [7,23], where the DC voltage, which is the input of inverter, is considered fixed.Since the dynamics of the voltage and current controller are much faster than that of the power controller, the combination of inverter, filter, primary sources and energy storages can be considered as a controllable voltage source.In addition, the bandwidth of the inverter voltage controller can be increased using several techniques, as described in [23].In this paper, we take an inverter as an ideal voltage source with controllable amplitude and frequency [24].The power flow in MG can be simplified to a two-node system, as shown in Figure 2, which can be written as sin sin( ( )) where  and  are the active and reactive power from node 1 to node 2, respectively; and  1 ,  1 and  2 , δ 2 are the voltage and phase of the two nodes, respectively.Z and θ are the impedance and angle of the transmission line, respectively; R and X is the resistance and inductance, respectively.From Equation (1), by multiplying Q by R and subtracting the product from multiplication of P and X, and multiplying Q by X and adding the product from multiplication of P and R, we can obtain: The power flow in MG can be simplified to a two-node system, as shown in Figure 2, which can be written as where P and Q are the active and reactive power from node 1 to node 2, respectively; and V 1 , δ 1 and V 2 , δ 2 are the voltage and phase of the two nodes, respectively.Z and θ are the impedance and angle of the transmission line, respectively; R and X is the resistance and inductance, respectively.
Energies 2016, 9, 243 3 of 14 inductance and primary source.The model details of each part are described in [7,23], where the DC voltage, which is the input of inverter, is considered fixed.Since the dynamics of the voltage and current controller are much faster than that of the power controller, the combination of inverter, filter, primary sources and energy storages can be considered as a controllable voltage source.In addition, the bandwidth of the inverter voltage controller can be increased using several techniques, as described in [23].In this paper, we take an inverter as an ideal voltage source with controllable amplitude and frequency [24].The power flow in MG can be simplified to a two-node system, as shown in Figure 2, which can be written as where  and  are the active and reactive power from node 1 to node 2, respectively; and  1 ,  1 and  2 , δ 2 are the voltage and phase of the two nodes, respectively.Z and θ are the impedance and angle of the transmission line, respectively; R and X is the resistance and inductance, respectively.From Equation (1), by multiplying Q by R and subtracting the product from multiplication of P and X, and multiplying Q by X and adding the product from multiplication of P and R, we can obtain: From Equation (1), by multiplying Q by R and subtracting the product from multiplication of P and X, and multiplying Q by X and adding the product from multiplication of P and R, we can obtain: (2) Energies 2016, 9, 243 4 of 14 When neglecting the line resistance (i.e., R = 0), Equation ( 2) can be rewritten as where δ is the power angle, i.e., δ " δ 1 ´δ2 .Therefore, the active and reactive power of DG can be controlled by drooping the frequency and voltage of the power system, respectively: where f is the output frequency; V d and V q are the output voltage in d-axis and q-axis, respectively; f 0 , V 0 are the initial frequency and voltage when the output power is 0, respectively; m and n are the droop parameters, respectively; P and Q are the active and reactive power of DG, respectively.Equation ( 4) presents the primary control of islanded MG.High gain angle droop control ensures proper load sharing, especially under weak system conditions.However, it has a negative impact on overall stability such as small signal stability [23].Furthermore, the voltage and frequency error always exist when operating.The secondary control based on the primary control is needed.In islanded MG, because the droop parameters are usually chosen based on the permitted frequency/voltage errors and DGs' power rating, the purpose of secondary control in islanded MG is to regulate f 0 and U 0 .
By differentiating the primary control characteristic in Equation ( 4), we can obtain .
where x f and x v are the differential value of frequency and voltage, respectively.Equations ( 5) and (6) show the first-order dynamic model of secondary control of islanded MG.The first-order system is a time-delay control system, and only the differential of frequency and voltage can be controlled so that the frequency and voltage are close to the target gradually.Then, the auxiliary variables, by expanding Equations ( 5) and (6), are used to precisely control the differential, that forms the secondary-order model: . .
where u c f and u cv are the auxiliary control variables of frequency and voltage, respectively.All DGs' frequency and voltage synchronize to the rated by designing u c f and u cv .Thus, the secondary-order model of MAS-based MG can be written as Energies 2016, 9, 243 x vi " u cvi , ¨¨# .
V dn " x vn .
x vn " u cvn (10) where the subscript of i and n are the sequence and total number of DGs in islanded MG.By properly designing the auxiliary variables of u c f i and u cvi , the MG can operate at the rated frequency and voltage.However, it should be noted that once the secondary frequency control is applied, the output active power of DG is distributed according to the same method used for primary control [23].After applying the primary control, the active power satisfies the following equation so that the active power of each DG is proportion to its rated value: In Equation ( 11), the droop coefficients m i are set based on the rated active power of DGs P Ni [23], which is equivalent to , where P N is the rated active power.
The secondary-order model of DG agent can be derived by differentiating Equation (11), and using the same expanding method of secondary frequency and voltage control: where u cP is the auxiliary control variable of active power control.Therefore, based on Equation ( 12), the model of active power allocation of islanded MG can be obtained as follows: x Pn " u cPn (13) Therefore, the active power allocation among DGs can be realized by designing the auxiliary control variable of u cPi .
The combination of Equations ( 9), ( 10) and ( 13) forms the secondary-order MAS model of secondary control for islanded MG.

Synergetic Secondary Control of Islanded MG Based on MAS
The secondary control of islanded MGs is a tracking synchronization problem, where all DGs try to synchronize their terminal voltage amplitude and frequency to the rated values.In MAS-based MG, all agents seek to synchronize to a virtual leader with the rated frequency and voltage (see Figure 3).In a distributed synergetic control, each agent needs to exchange the information with its neighboring DGs.In MAS, the spare communication network in MG can be modelled by digraph theorem.
all agents seek to synchronize to a virtual leader with the rated frequency and voltage (see Figure 3).In a distributed synergetic control, each agent needs to exchange the information with its neighboring DGs.In MAS, the spare communication network in MG can be modelled by digraph theorem.

Graph Theorem
The communication network of MG can be modelled by a digraph [25].In MG, a DG is considered as a node of the communication digraph, and the edges of the communication network symbolize the communication links.A digraph is usually expressed as  = (, , ) with a nonempty finite set of N nodes V = { 1 ,  2 , ⋯ ,   }, a set of edges or arcs  ∈ (V × V) that is the Cartesian product of V, and the associated adjacency matrix A = [  ] ∈  × .In this paper, the digraph is assumed to be time-invariant, and then, A is constant.An edge from node j to node i is denoted by (  ,   ), which means that node i receives the information from node j.In a digraph, (  ,   ) denotes the edge from node i itself, which is usually neglected.  is the weight of edge (  ,   ), and   > 0 if (  ,   ) ∈ , otherwise   = 0. Node j is called a neighbor of node i if   > 0. The set of neighbors of node i is denoted as   = {  |(  ,   ) ∈ ε}.In MG, since all DGs exchange information through communication links, the non-directed graph is considered, that node i can receive/send the information from/to node j at the same time.The Laplacian matrix of the digraph is defined as  = [  ] × , where   = −   if  ≠ .

Active Power Control
The active power of DG is allocated according to the droop coefficients used for primary control of islanded MG.Therefore, the active power allocation for secondary control is designed firstly.In this section, the synergetic control is used to design the auxiliary control variables.
To achieve precise power allocation among DGs, it is assumed that DGs communicate with their neighbors through a prescribed communication digraph .Since the active power of DG can not be predicted, the synergetic control with leaderless is adopted in this paper.The auxiliary controls   are chosen based on the own information of each DG and the information of its neighbors in Figure 3 as follows: ()

Graph Theorem
The communication network of MG can be modelled by a digraph [25].In MG, a DG is considered as a node of the communication digraph, and the edges of the communication network symbolize the communication links.A digraph is usually expressed as G " pV, ε, Aq with a non-empty finite set of N nodes V " tv 1 , v 2 , ¨¨¨, v N u, a set of edges or arcs " P pV ˆVq that is the Cartesian product of V, and the associated adjacency matrix A " " a ij ‰ P R NˆN .In this paper, the digraph is assumed to be time-invariant, and then, A is constant.An edge from node j to node i is denoted by `vj , v i ˘, which means that node i receives the information from node j.In a digraph, pv i , v i q denotes the edge from node i itself, which is usually neglected.a ij is the weight of edge `vj , v i ˘, and In MG, since all DGs exchange information through communication links, the non-directed graph is considered, that node i can receive/send the information from/to node j at the same time.The Laplacian matrix of the digraph is defined as

Active Power Control
The active power of DG is allocated according to the droop coefficients used for primary control of islanded MG.Therefore, the active power allocation for secondary control is designed firstly.In this section, the synergetic control is used to design the auxiliary control variables.
To achieve precise power allocation among DGs, it is assumed that DGs communicate with their neighbors through a prescribed communication digraph G. Since the active power of DG can not be predicted, the synergetic control with leaderless is adopted in this paper.The auxiliary controls u cP are chosen based on the own information of each DG and the information of its neighbors in Figure 3 as follows: where u cPαi and u cPβi are two parts of the auxiliary control u cPi , which are given in Equations (14a,b), respectively.In the synergetic control, the first part is to directly control the active power of DGs to allocate according to the droop control and the last one is to control the errors of pmPq in each DG The solutions of Equation ( 13) px Pi " 0q show that the equilibrium point of the control is u cPi = 0. Thus, the output active power of DG is allocated in inverse proportion to the droop coefficients.Since the active power allocation is realized, the system stability is also needed to be taken into consideration.
The small signal stability based on the linearization model [23,24] and bifurcation method [26][27][28] based on nonlinear analysis method are two common-used methods for MG.In this paper, the extended Lyapunov theorem is adopted to analyze the globally asymptotic stability.
In control system, the behavior of a system about its equilibrium point can be studied by Lyapunov's stability theorem [29,30].Lyapunov's stability theorem allows us to determine the stability for the MG system by analyzing an energy function, which is also called as the Lyapunov function V pxq.According to Lyapunov's direct method, the equilibrium point is globally asymptotically stable if V pxq satisfies the following four properties: V pxq ą 0, for all x ‰ 0; V pxq Ñ 8 as x Ñ 8; V pxq ă 0 for all x ‰ 0.
To analyze the stability of the active power control in Equation ( 14), the following Lyapunov function is considered.
Vpxq " 1 2 Obviously, the first three conditions in Lyapunov's theorem are always satisfied.Then, only the last condition is needed to analyze.
By differentiating the Lyapunov function in Equation ( 17), we obtain the following equation: . Vpxq " x Pi .
x Pi `1 2 Then, add the secondary-order MAS model in Equation (13) to Equation (17) .
Vpxq " The first and third parts in Equation ( 18) are always satisfied by the following equation, i.e., Then, Equation ( 18) can be simplified to .
Vpxq " Energies 2016, 9, 243 8 of 14 Consider the signum function in Equation ( 15), we can obtain the result that the control system in Equation ( 14) satisfies the following function: . V pxq ď 0. ( Thus, all the conditions of Lyapunov's theorem are satisfied, i.e., the stability of synergetic control in Equation ( 14) leads to globally asymptotic stability.

Secondary Voltage Control
In this section, a distributed synergetic control is designed to control all the voltage magnitudes of DGs V di to the reference voltage V re f .Because the reference of secondary voltage is given at first, a virtual leader with the rated value is set to design the synergetic control.The auxiliary controls u cv are chosen based on the own information of each DG and the information of its neighbors in Figure 3 as follows: where the subscript L means the virtual leader; u cvαi , u cvβi and u cvγi are three parts of the auxiliary control u cvi that are given from Equations (22a) to (22c), respectively.In Equation (22), Equation (22a) is to directly control the DG voltage to work to the same value by using its own and the neighbors' voltage; Equation (22b) is to control the voltage errors to be the same, and Equation (22c) is to control the DG to synchronize to the virtual leader.The non-negative coefficients of k v1 , k v2 and k v3 are three control gains, i.e., k v1 , k v2 , k v3 ě 0. Similarly, we use the following Lyapunov function to analyze the system stability: According to the analysis results in Equation ( 16), the above function satisfies the first three conditions in Lyapunov's theorem.By differentiating Equation ( 23), we obtain .Vpxq " Then, adding the secondary-order MAS model in Equation (10) to Equation (24), and utilizing the analysis results in Equation ( 19), we can obtain .Vpxq " Suppose that xvi " x vi ´xvL .Then, we can obtain Energies 2016, 9, 243 9 of 14 Adding Equations ( 19) and (26) to Equation (25) yields: . Vpxq " According to the analysis results in Equation ( 20), the secondary voltage control always satisfies the conditions in Lyapunov's theorem.Therefore, the synergetic secondary voltage control has globally asymptotic stability.
The block diagram of the secondary voltage control based on the distributed synergetic control is shown in Figure 4.The control input V oi is where V oi0 is the initial value of V oi in the primary control, and Q i is the reactive power in the primary control [23].
Energies 2016, 9, 243 9 of 14 According to the analysis results in Equation ( 20), the secondary voltage control always satisfies the conditions in Lyapunov's theorem.Therefore, the synergetic secondary voltage control has globally asymptotic stability.
The block diagram of the secondary voltage control based on the distributed synergetic control is shown in Figure 4.The control input   is 0 () where  0 is the initial value of   in the primary control, and   is the reactive power in the primary control [23].

Secondary Frequency Control
In this section, a distributed synergetic control is designed to control the frequency in each DG to the reference value of   .To complete the target, it is assumed that DGs can communicate with each other through a prescribed communication digraph .The auxiliary variable   is chosen based on its own information and its neighbors in Figure 3 as follows: where   ,      are three parts of the auxiliary control   that are given from Equation (29a) to (29c), respectively.In Equation ( 29), Equation (29a) is to directly control the frequency to work to the same value by using its own and the neighbors' information; Equation (29b) is to control the frequency error to be the same, and Equation ( 29c) is to control the DG to synchronize with the virtual leader.The non-negative coefficients  1 ,  2 and  3 are three control gains, i.e.  1 ,  2 ,  3 ≥ 0.

Secondary Frequency Control
In this section, a distributed synergetic control is designed to control the frequency in each DG to the reference value of f N .To complete the target, it is assumed that DGs can communicate with each other through a prescribed communication digraph G.The auxiliary variable u c f i is chosen based on its own information and its neighbors in Figure 3 as follows: where u c f αi , u c f βi and u c f γi are three parts of the auxiliary control u c f i that are given from Equation (29a) to (29c), respectively.In Equation ( 29), Equation (29a) is to directly control the frequency to work to the same value by using its own and the neighbors' information; Equation (29b) is to control the frequency error to be the same, and Equation (29c) is to control the DG to synchronize with the virtual leader.The non-negative coefficients k f 1 , k f 2 and k f 3 are three control gains, i.e. k f 1 , k f 2 , k f 3 ě 0.
According to the analysis results in the secondary voltage control, we can obtain that the secondary frequency control in Equation ( 29) also has globally asymptotic stability.Therefore, the secondary frequency control input f oi is where f oi0 is the initial value of f oi in the primary control.
Because the active power of DGs needs to be allocated according to Equation ( 11), the active power P i in secondary frequency control is obtained from Equation ( 14), but not the value calculated from the DG's output current and voltage.Then, the whole secondary frequency is Based on the distributed synergetic control, the block diagram of the secondary frequency control is illustrated in Figure 5.
Energies 2016, 9, 243 10 of 14 According to the analysis results in the secondary voltage control, we can obtain that the secondary frequency control in Equation ( 29) also has globally asymptotic stability.Therefore, the secondary frequency control input   is 0 () where  0 is the initial value of   in the primary control.
Because the active power of DGs needs to be allocated according to Equation ( 11), the active power   in secondary frequency control is obtained from Equation ( 14), but not the value calculated from the DG's output current and voltage.Then, the whole secondary frequency is 0 () Based on the distributed synergetic control, the block diagram of the secondary frequency control is illustrated in Figure 5.

Simulation
An islanded MG illustrated in Figure 6 is used to verify the effectiveness of the synergetic secondary control [22], which consists of four DGs.The DG is considered as a controllable voltage source.The rated power of DG1 and DG2 is 100 kVA, and that of DG3 and DG4 is 75 kVA.The lines are modelled as series RL branches.The parameters of the DGs, lines, and loads are given in Table 1.The dynamic simulation is done in EMTP, and the results are analyzed with Matlab.

Simulation
An islanded MG illustrated in Figure 6 is used to verify the effectiveness of the synergetic secondary control [22], which consists of four DGs.The DG is considered as a controllable voltage source.The rated power of DG1 and DG2 is 100 kVA, and that of DG3 and DG4 is 75 kVA.The lines are modelled as series RL branches.The parameters of the DGs, lines, and loads are given in Table 1.The dynamic simulation is done in EMTP, and the results are analyzed with Matlab.
Energies 2016, 9, 243 10 of 14 According to the analysis results in the secondary voltage control, we can obtain that the secondary frequency control in Equation ( 29) also has globally asymptotic stability.Therefore, the secondary frequency control input   is 0 () where  0 is the initial value of   in the primary control.
Because the active power of DGs needs to be allocated according to Equation ( 11), the active power   in secondary frequency control is obtained from Equation ( 14), but not the value calculated from the DG's output current and voltage.Then, the whole secondary frequency is 0 () Based on the distributed synergetic control, the block diagram of the secondary frequency control is illustrated in Figure 5.

Simulation
An islanded MG illustrated in Figure 6 is used to verify the effectiveness of the synergetic secondary control [22], which consists of four DGs.The DG is considered as a controllable voltage source.The rated power of DG1 and DG2 is 100 kVA, and that of DG3 and DG4 is 75 kVA.The lines are modelled as series RL branches.The parameters of the DGs, lines, and loads are given in Table 1.The dynamic simulation is done in EMTP, and the results are analyzed with Matlab.It is assumed that the DGs communicate with each other as depicted in Figure 3.The Laplace matrix of the digraph is The reference value for the MG frequency and terminal amplitude voltage of DGs are f re f " 50Hz and V dire f " 311V.The control gains for secondary control are given in Table 2.It is assumed that the islanded MG starts to work with the primary control, and the secondary control works at t = 1.2 s.The secondary frequency and voltage control results are shown in Figures 7 and 8 respectively.It is assumed that the DGs communicate with each other as depicted in Figure 3.The Laplace matrix of the digraph is The reference value for the MG frequency and terminal amplitude voltage of DGs are   = 50 and   = 311.The control gains for secondary control are given in Table 2.It is assumed that the islanded MG starts to work with the primary control, and the secondary control works at t = 1.2 s.The secondary frequency and voltage control results are shown in Figures 7 and 8, respectively.

Discussion
As seen in Figure 7, the islanded MG starts to work with primary control and goes to stability state after 0.8 s.In Figure 7a, the operating frequencies all go to a common value (49.3 Hz) less than the rated, which is the frequency of MG.Due to the primary control, the output active power of the four DGs is allocated according to their droop parameters.After 1.2 s, the secondary frequency control in Equation ( 31) is applied, and the operating frequency of the islanded MG returns to its rated value after about 0.3 s, shown in Figure 7a.As seen in Figure 7b, the output active power of DGs is allocated in inverse proportion to the droop parameters because the distributed active power control in Equation ( 14) is applied.
Figure 8 presents the voltage and reactive power of DGs.Similar to the frequency control, the terminal voltages of DGs works with a different value less than the rated   .Once the synergetic secondary control is applied, all the DGs' voltage   returns to the rated value, as shown in Figure 8a.Because of the application of primary in Equation ( 4), the quadratic term of DG voltage is 0. Therefore, the active and reactive power of DGs is re-allocated according to the primary control, which corresponds to the results shown in Figure 7b and Figure 8b, respectively.
In power systems, since the system frequency error can be used to measure the active power deviation, the centralized control based on the secondary control in traditional power grid can be used to restore the system frequency, which has been verified in [9][10][11].However, the system with centralized control may be unstable based on the analysis results in [23], which is globally asymptotic stability under the proposed control.Meanwhile, the secondary voltage control is very complex through the centralized control for the following reasons.Firstly, the reactive power deviation of the grid cannot be derived directly through the same way because the voltage at each node is not the same.In addition, even though the reactive power deviation is obtained, its allocation is still a complex problem because the power flow is different to the traditional power grid while the balance point of the grid does not exist in droop-controlled MG.

Conclusions
In this paper, the distributed synergetic control based on MAS is proposed to design the secondary control of MG.In MAS, DGs communicate with their nearest neighbors through a communication network modelled by digraph theorem.Based on primary controls, the secondaryorder differential-algebraic equations are used to describe the dynamics of agents.The synergetic control is used to design the distributed secondary frequency and voltage control.In the proposed control, the auxiliary variables are designed so that each DG only requires its own and its neighbors' information in the communication network.The Lyapunov theorem is used to verify the globally asymptotic stability.In the proposed secondary control, since the complex communication network is avoided, the system is more reliable.The distributed control structure provides a robust secondary

Discussion
As seen in Figure 7, the islanded MG starts to work with primary control and goes to stability state after 0.8 s.In Figure 7a, the operating frequencies all go to a common value (49.3 Hz) less than the rated, which is the frequency of MG.Due to the primary control, the output active power of the four DGs is allocated according to their droop parameters.After 1.2 s, the secondary frequency control in Equation ( 31) is applied, and the operating frequency of the islanded MG returns to its rated value after about 0.3 s, shown in Figure 7a.As seen in Figure 7b, the output active power of DGs is allocated in inverse proportion to the droop parameters because the distributed active power control in Equation ( 14) is applied.
Figure 8 presents the voltage and reactive power of DGs.Similar to the frequency control, the terminal voltages of DGs works with a different value less than the rated V re f .Once the synergetic secondary control is applied, all the DGs' voltage V odi returns to the rated value, as shown in Figure 8a.Because of the application of primary in Equation ( 4), the quadratic term of DG voltage is 0. Therefore, the active and reactive power of DGs is re-allocated according to the primary control, which corresponds to the results shown in Figures 7b and 8b, respectively.
In power systems, since the system frequency error can be used to measure the active power deviation, the centralized control based on the secondary control in traditional power grid can be used to restore the system frequency, which has been verified in [9][10][11].However, the system with centralized control may be unstable based on the analysis results in [23], which is globally asymptotic stability under the proposed control.Meanwhile, the secondary voltage control is very complex through the centralized control for the following reasons.Firstly, the reactive power deviation of the grid cannot be derived directly through the same way because the voltage at each node is not the same.In addition, even though the reactive power deviation is obtained, its allocation is still a complex problem because the power flow is different to the traditional power grid while the balance point of the grid does not exist in droop-controlled MG.

Conclusions
In this paper, the distributed synergetic control based on MAS is proposed to design the secondary control of MG.In MAS, DGs communicate with their nearest neighbors through a communication network modelled by digraph theorem.Based on primary controls, the secondary-order differential-algebraic equations are used to describe the dynamics of agents.The synergetic control is used to design the distributed secondary frequency and voltage control.In the proposed control, the auxiliary variables are designed so that each DG only requires its own and its neighbors' information in the communication network.The Lyapunov theorem is used to verify the globally asymptotic stability.In the proposed secondary control, since the complex communication network is avoided, the system is more reliable.The distributed control structure provides a robust secondary control framework, which is able to appropriately operate in time varying and unreliable communication networks.

Figure 2 .
Figure 2. Power flow in MGs.(a) voltage vector and (b) power flow

Figure 2 .
Figure 2. Power flow in MGs.(a) voltage vector and (b) power flow

Figure 2 .
Figure 2. Power flow in MGs.(a) voltage vector and (b) power flow.

Figure 4 .
Figure 4. Diagram of the synergetic secondary voltage control.

Figure 4 .
Figure 4. Diagram of the synergetic secondary voltage control.

Figure 5 .
Figure 5. Diagram of synergetic secondary frequency control.

Figure 5 .
Figure 5. Diagram of synergetic secondary frequency control.

Figure 5 .
Figure 5. Diagram of synergetic secondary frequency control.

Figure 7 .
Figure 7. Simulation results of the secondary frequency control.(a) the islanded MG frequency and (b) the active power of DGs.Figure 7. Simulation results of the secondary frequency control.(a) the islanded MG frequency and (b) the active power of DGs.

Figure 7 .Figure 8 .
Figure 7. Simulation results of the secondary frequency control.(a) the islanded MG frequency and (b) the active power of DGs.Figure 7. Simulation results of the secondary frequency control.(a) the islanded MG frequency and (b) the active power of DGs.

Figure 8 .
Figure 8. Simulation results of the secondary voltage control.(a) the amplitude voltage and (b) the reactive power of DGs.
to be the same.The non-negative coefficients k P1 and k P2 are two control gains, i.e., k P1 , k P2 ě 0.Moreover, sign means the signum function, i.e., is

Table 1 .
Parameters of the MG system.

Table 2 .
Parameters of the primary and secondary control.

Table 1 .
Parameters of the MG system.

Table 2 .
Parameters of the primary and secondary control.