Distributed Operation of Microgrids Considering Secondary Frequency Restoration Based on the Di ﬀ usion Algorithm

: This paper proposes a distributed control of the microgrid (MG) system based on the di ﬀ usion algorithm. Unlike the existing decentralized strategy that focuses on the economic operation of the MG system, the proposed strategy performs secondary frequency regulation in addition to the optimization of the MG system. The hierarchical control technique is employed in this study, where the primary layer is responsible for power control and the secondary layer is responsible for the frequency control and economic operation of the MG system. A tested MG system with four distributed generations (DGs) is considered. Three types of communication topologies are evaluated in this study, which are line, ring, and full topologies. The proposed controller is compared to the conventional consensus controller to show the e ﬀ ectiveness of the proposed di ﬀ usion controller. Simulation results show that the proposed di ﬀ usion strategy improves the convergence speed of the distributed control, resulting in the improvement of power responses and frequency quality of the MG system. The tested system is implemented in the MATLAB / Simulink environment to show the feasibility of the proposed di ﬀ usion controller.

A distributed control strategy of the MG system is proposed, which maintains both economic operation and frequency restoration of the MG system in a single layer. The proposed strategy avoids the load prediction step by employing the optimization of frequency deviation, which overcomes the problem of load uncertainty in the optimization problem. A diffusion-based distributed control is proposed to improve the convergence speed of the distributed control, leading to an improvement in power response and frequency regulation.
The proposed method is based on the hierarchical strategy with two layers, in which the secondary layer is responsible for the economic operation of the MG system and the secondary frequency regulation. The typical droop controllers are used in the local controller to achieve initial power sharing among DGs. Then, the secondary controller based on the diffusion algorithm calculates the additional power of each DGs to optimize the economic operation of the MG system. The MG system consists of four DGs used to verify the proposed method. The effectiveness of the proposed diffusion-based MG controller is evaluated by testing the proposed method with three communication topologies-line, ring, and full connections. In addition, a comparison study on the proposed diffusion-based MG controller and the conventional consensus-based controller is presented in this study.
The paper organization is as follows. The hierarchical control techniques of the MG system is described in Section 2. The proposed distributed operation of the MG system based on the diffusion Energies 2020, 13, 3207 3 of 14 algorithm is presented in Section 3. Simulation results are provided in Section 4. Finally, the main conclusion is summarized in Section 5.

Control Strategies of the MG System
The hierarchical control of the MG system shown in Figure 1a consists of three layers, which are the primary layer, secondary layer, and tertiary layer. The primary layer is the local controller of each DG and it requires only local DG information. The primary layer is responsible for local voltage control and proper power sharing between multiple DGs. Since the rotating inertia of the MG system is low due to the lack of synchronous generators, the primary controller must be fast (i.e., operating in milliseconds) to ensure stability during operations of the MG system. Furthermore, in order to enhance the reliability of the MG, the communication system could be neglected in the primary control. The secondary and tertiary control layers support the MG operation with multiple objectives. A low-bandwidth communication link is often used in the secondary layer to achieve the global controllability. It could be based on the centralized or decentralized strategies with the objective of frequency of voltage restoration or other objectives related to the voltage control and power qualities.
DG and it requires only local DG information. The primary layer is responsible for local voltage control and proper power sharing between multiple DGs. Since the rotating inertia of the MG system is low due to the lack of synchronous generators, the primary controller must be fast (i.e., operating in milliseconds) to ensure stability during operations of the MG system. Furthermore, in order to enhance the reliability of the MG, the communication system could be neglected in the primary control. The secondary and tertiary control layers support the MG operation with multiple objectives. A low-bandwidth communication link is often used in the secondary layer to achieve the global controllability. It could be based on the centralized or decentralized strategies with the objective of frequency of voltage restoration or other objectives related to the voltage control and power qualities. In the MG system, the tertiary layer could be utilized for economic optimization considering the electricity markets and cost of power. The tertiary control layer in the conventional hierarchical control is commonly based on the centralized methods, which offer the benefit of easy implementation. However, the central controller suffers the computation burden when the system scale significantly increases. In addition, the reliability of the control system is low since all DGs in the MG system rely on a single central controller.
The problem is overcome by the proposed distributed controller shown in Figure 1b. The secondary control layer in the proposed strategy is responsible for both economic operation and frequency restoration of the MG system. The sparse communication network is used to interchange information among DGs. All DGs are coordinated to optimize the economic problem of the MG system, which results in a reduction in computation burden and an improvement in the reliability of the control system. It is well-know that the power balance is achieved when the system frequency is regulated at the nominal value. The distributed optimization in the proposed strategy is subjected to the constraint of frequency deviation, which results in the frequency restoration of the MG system. The flow diagram of the proposed distributed controller is shown in Figure 2. Optimal power indicates that the minimizing of the operation cost and frequency deviation is calculated in the secondary layer by the diffusion algorithm. Each DG receives output power information from its neighbors ( ) and locally estimates the system frequency by the phase-locked loop (PLL). Based on this information, the diffusion algorithm is used to calculate the incremental cost of each DG, and the optimal power reference ( * ) is found accordingly. The optimal power reference is sent to the local controller where the conventional droop power controller is employed to manage the output power In the MG system, the tertiary layer could be utilized for economic optimization considering the electricity markets and cost of power. The tertiary control layer in the conventional hierarchical control is commonly based on the centralized methods, which offer the benefit of easy implementation. However, the central controller suffers the computation burden when the system scale significantly increases. In addition, the reliability of the control system is low since all DGs in the MG system rely on a single central controller.
The problem is overcome by the proposed distributed controller shown in Figure 1b. The secondary control layer in the proposed strategy is responsible for both economic operation and frequency restoration of the MG system. The sparse communication network is used to interchange information among DGs. All DGs are coordinated to optimize the economic problem of the MG system, which results in a reduction in computation burden and an improvement in the reliability of the control system. It is well-know that the power balance is achieved when the system frequency is regulated at the nominal value. The distributed optimization in the proposed strategy is subjected to the constraint of frequency deviation, which results in the frequency restoration of the MG system.
The flow diagram of the proposed distributed controller is shown in Figure 2. Optimal power indicates that the minimizing of the operation cost and frequency deviation is calculated in the secondary layer by the diffusion algorithm. Each DG receives output power information from its Energies 2020, 13, 3207 4 of 14 neighbors (P j ) and locally estimates the system frequency by the phase-locked loop (PLL). Based on this information, the diffusion algorithm is used to calculate the incremental cost of each DG, and the optimal power reference (P * i ) is found accordingly. The optimal power reference is sent to the local controller where the conventional droop power controller is employed to manage the output power of each DG. Finally, the output power of each DG is measured and broadcasted to its neighbors. Iterating this process leads to the convergence of all incremental costs to a common value, resulting in the optimal condition of the MG system. The local control and secondary control layers in the proposed strategy are explained in next section.
Energies 2020, 13, x FOR PEER REVIEW 4 of 14 of each DG. Finally, the output power of each DG is measured and broadcasted to its neighbors. Iterating this process leads to the convergence of all incremental costs to a common value, resulting in the optimal condition of the MG system. The local control and secondary control layers in the proposed strategy are explained in next section.

Proposed Distributed MG Control
The MG system could include various types of distributed generations (DGs), such as photovoltaic, fuel cells, and wind generations. Such DGs are interfaced with the MG system by the power converters to improve aspects of their controllability, such as power or voltage regulation. These configurations of such DGs are known as inverter-based DGs. A typical configuration of an inverter-based DG is shown in Figure 3, which is based on the two-level voltage source converter with an output LC filter connected with a constant DC source. The proposed controller for each DG consists of two layers: a primary and secondary, where the primary layer is responsible for power control and the secondary layer is in charge of frequency restoration and economic optimization.

Proposed Distributed MG Control
The MG system could include various types of distributed generations (DGs), such as photovoltaic, fuel cells, and wind generations. Such DGs are interfaced with the MG system by the power converters to improve aspects of their controllability, such as power or voltage regulation. These configurations of such DGs are known as inverter-based DGs. A typical configuration of an inverter-based DG is shown in Figure 3, which is based on the two-level voltage source converter with an output LC filter connected with a constant DC source. The proposed controller for each DG consists of two layers: a primary and secondary, where the primary layer is responsible for power control and the secondary layer is in charge of frequency restoration and economic optimization.

Primary Control of DGs
The primary control of DG is shown in Figure 4, which includes the inner current/voltage control and outer power control loops. The proportional integral (PI) regulators are used to regulate the converter current and voltage. The current controller is based on the decoupling of a dq-reference frame using Equations (1) and (2) [25].
where and are PI control coefficients of the inner current control loop of d axes, and are PI control coefficients of the inner current control loop of q axes, * and * are the stator currents of the and reference of the d and q axes, and * and * are the output voltage references.
The outer control loop is based on the droop control mechanism to manage the power sharing between multiple DGs in the MG system. The droop controller is used to share load power demand through proportional constants and ensure angular frequency/voltage stability. The droop control loop generates the references of angular frequency and voltage to the inner voltage control loop, as given by Equations (3) and (4) where * and ∆ * are the angular frequency and amplitude deviation of voltage, is the nominal frequency, and and are v, respectively. * and ∆ * are the reference values for the DG's reactive power and deviation of active power, respectively, and and and represent the measurements. Thus, the angular frequency is dependent on active power and the difference in voltage amplitude is dependent on reactive power.

Primary Control of DGs
The primary control of DG is shown in Figure 4, which includes the inner current/voltage control and outer power control loops. The proportional integral (PI) regulators are used to regulate the converter current and voltage. The current controller is based on the decoupling of a dq-reference frame using Equations (1) and (2) [25].
where k P1 and k I1 are PI control coefficients of the inner current control loop of d axes, k P2 and k I2 are PI control coefficients of the inner current control loop of q axes, i * gd and i * gq are the stator currents of the i gd and i gq reference of the d and q axes, and v * gd and v * gq are the output voltage references.

Optimal Operation of the MG System
It is assumed that there are n distributed generators in the MG system. The output power of i-th DGi is equal to , and is the load demand of the MG system. The cost function for each DG can be expressed by the quadratic function ( ), as shown in Equation (5) [26].
where , and are cost coefficients. The optimal output powers of DGs are found by minimizing the objective function in Equation (6). The outer control loop is based on the droop control mechanism to manage the power sharing between multiple DGs in the MG system. The droop controller is used to share load power demand through proportional constants and ensure angular frequency/voltage stability. The droop control loop generates the references of angular frequency and voltage to the inner voltage control loop, as given by Equations (3) and (4).
Energies 2020, 13, 3207 6 of 14 where θ * and ∆V * gd are the angular frequency and amplitude deviation of voltage, f is the nominal frequency, and k p and k q are v, respectively. Q * and ∆P * are the reference values for the DG's reactive power and deviation of active power, respectively, and P and Q and represent the measurements. Thus, the angular frequency is dependent on active power and the difference in voltage amplitude is dependent on reactive power.

Optimal Operation of the MG System
It is assumed that there are n distributed generators in the MG system. The output power of i-th DG i is equal to P i , and P D is the load demand of the MG system. The cost function for each DG can be expressed by the quadratic function C i (P i ), as shown in Equation (5) [26].
where α i , β i and γ i are cost coefficients.
The optimal output powers of DGs are found by minimizing the objective function C in Equation (6).
Subject to the constraint's power balance and power capacity.
The Lagrange multiplier method is applied to solve the objective function C. The Lagrange function L is given by the objective function C and the constraint function multiplied by the undetermined multiplier λ, as shown in Equation (9).
Considering the derivative of L with respect to each variable, we have: Setting the derivatives of Equations (10) and (11) to 0 gives the optimal solution. Thus, the cost operation of the MG system reaches to minimum when the IC of each DG is equal to λ, as given by Equation (12) [27].

Proposed Diffusion-Based Distributed Control of the MG System
It can be observed in previous section that the optimal operation of the MG system is maintained if the incremental cost of each DG in the MG system converges to a common value. In the secondary control layer, the communication network is utilized to inform the neighbor-to-neighbor information trade and to find the optimal value of the incremental cost. In the mechanism, each DG will send their respective output powers P i to neighborhoods, and receive information from those neighborhoods with DG j (P j ) with j ∈ N i , where N i is the neighborhood of DG i . The undirected graph G N = (N N , E N ), with the set of nodes N N = {n 1 , n 2 , . . . , n N } representing the N controllers of DGs, and the communication links between controllers represented by the set of edges E N ⊆ N N × N N express the communication network established between the DGs. The edge n i , n j in the aforementioned edge set E N shows that node n receives information from node n j . The undirected graph G N defined as as A N = a ij ∈ R N×N is also associated with the adjacency matrix, where a ij is a positive value supposing that if n i , n j ∈ E N and a ij = 0 otherwise. The Laplacian matrix of the undirected graph G N is noted as L N = l ij ∈ R N×N , where l ii = N j=1, j i a ij and l ij = −a ij is for j i. The global dynamic of the combine-then-adapt diffusion strategy is as shown in Equation (14) [28], where X k = x 1,k ;  (14) is used in each DG to find the optimal incremental cost. The secondary control objective of frequency is added in the diffusion algorithm to improve the frequency of the MG system. Unlike the conventional distributed operation strategy, the constraint of power balance is maintained by regulating the system frequency at the nominal value. The proposed diffusion algorithm is given by the following: where are the incremental cost of a DG at time k; f is the measured system frequency; and f * is the nominal frequency of the MG system. The optimal power of each DG is found according to the corresponding incremental cost, as given by In the secondary diffusion control, the control block diagram is supposed as Figure 5, where the inputs are the power information of the adjacent DGs and the output is the optimal power reference for the local controller.

Simulation Results
In this study, the tested MG system is illustrated in Figure 6, consisting of four DGs and three loads. The distribution lines are modelled by the resistance/inductance components. Parameters of the DGs and the tested MG system are shown in Tables 1 and 2. Table 3 shows the converter parameters used in this study. Since this study focuses on the economic operation and secondary frequency regulation of the MG system, the converter modelling is based on the average model. The designs of the PI gains of current and voltage control are based on the pole-placement method, as discussed in [29]. The three communication topologies (line, ring and full topology) of this study are given in Figure 7. In the MATLAB/Simulink environment, the communication network is modelled using From/Goto blocks, in which the Goto block is used to transmit information and the From block is used to receive information. The scenario of step load change is used to verify the performance of the proposed control strategy. It is assumed that the initial load is 8 kW and an additional load of 4 kW is connected to the MG system at 5 s.  Figure 5. Secondary diffusion controller considering frequency restoration.

Simulation Results
In this study, the tested MG system is illustrated in Figure 6, consisting of four DGs and three loads. The distribution lines are modelled by the resistance/inductance components. Parameters of the DGs and the tested MG system are shown in Tables 1 and 2. Table 3 shows the converter parameters used in this study. Since this study focuses on the economic operation and secondary frequency regulation of the MG system, the converter modelling is based on the average model. The designs of the PI gains of current and voltage control are based on the pole-placement method, as discussed in [29]. The three communication topologies (line, ring and full topology) of this study are given in Figure 7. In the MATLAB/Simulink environment, the communication network is modelled using From/Goto blocks, in which the Goto block is used to transmit information and the From block is used to receive information. The scenario of step load change is used to verify the performance of the proposed control strategy. It is assumed that the initial load is 8 kW and an additional load of 4 kW is connected to the MG system at 5 s.  When solving the incremental cost (IC) of each DG using the Lagrange multiplier method, the optimal IC (λ * ) with a power demand of 12 kW is 17.5 USD/kW. The corresponding powers of the four DGs are 4.35 kW, 3.52 kW, 1.4 kW, and 2.73 kW, respectively. Figures 8-10 show the optimal IC found by the distributed controllers, and it is consistent with the Lagrange multiplier method. The optimal IC found by the distributed control is 17.5 USD/kW after several iterations. The convergence of ICs depends on the communication topologies and the control methods. The optimal IC is found after 1 s with the full connection, whereas it takes more than 5 s in the case of line connection. The advantage of the proposed distributed control based on the diffusion algorithm can be seen in Figure  8, which shows a faster convergence speed compared to the conventional consensus method. The     When solving the incremental cost (IC) of each DG using the Lagrange multiplier method, the optimal IC (λ * ) with a power demand of 12 kW is 17.5 USD/kW. The corresponding powers of the four DGs are 4.35 kW, 3.52 kW, 1.4 kW, and 2.73 kW, respectively. Figures 8-10 show the optimal IC found by the distributed controllers, and it is consistent with the Lagrange multiplier method. The optimal IC found by the distributed control is 17.5 USD/kW after several iterations. The convergence of ICs depends on the communication topologies and the control methods. The optimal IC is found after 1 s with the full connection, whereas it takes more than 5 s in the case of line connection. The advantage of the proposed distributed control based on the diffusion algorithm can be seen in Figure  8, which shows a faster convergence speed compared to the conventional consensus method. The optimal IC is found at 8.9 s, whereas it is found at more than 5 s in the case of the conventional controller.  When solving the incremental cost (IC) of each DG using the Lagrange multiplier method, the optimal IC (λ * ) with a power demand of 12 kW is 17.5 USD/kW. The corresponding powers of the four DGs are 4.35 kW, 3.52 kW, 1.4 kW, and 2.73 kW, respectively. Figures 8-10 show the optimal IC found by the distributed controllers, and it is consistent with the Lagrange multiplier method. The optimal IC found by the distributed control is 17.5 USD/kW after several iterations. The convergence of ICs depends on the communication topologies and the control methods. The optimal IC is found after 1 s with the full connection, whereas it takes more than 5 s in the case of line connection. The advantage of the proposed distributed control based on the diffusion algorithm can be seen in Figure 8, which shows a faster convergence speed compared to the conventional consensus method. The optimal IC is found at 8.9 s, whereas it is found at more than 5 s in the case of the conventional controller.  When solving the incremental cost (IC) of each DG using the Lagrange multiplier method, the optimal IC (λ * ) with a power demand of 12 kW is 17.5 USD/kW. The corresponding powers of the four DGs are 4.35 kW, 3.52 kW, 1.4 kW, and 2.73 kW, respectively. Figures 8-10 show the optimal IC found by the distributed controllers, and it is consistent with the Lagrange multiplier method. The optimal IC found by the distributed control is 17.5 USD/kW after several iterations. The convergence of ICs depends on the communication topologies and the control methods. The optimal IC is found after 1 s with the full connection, whereas it takes more than 5 s in the case of line connection. The advantage of the proposed distributed control based on the diffusion algorithm can be seen in Figure  8, which shows a faster convergence speed compared to the conventional consensus method. The optimal IC is found at 8.9 s, whereas it is found at more than 5 s in the case of the conventional controller.   The convergence speed of the ICs has an impact on the response of the output powers, which is represented in Figures 11-13. Load changing at 5 s causes the transient output power. The steadystate output powers of DGs shown in these figures are consistent with the Lagrange multiplier method. It can be observed that fast convergence speed results in a fast response of the output powers. The communication topology of full connection provides the fastest power responses compared to the line and ring connections. With the proposed diffusion-based controller, the responses of output powers are much quicker than the conventional consensus method, particularly in the case of the line connection.  The convergence speed of the ICs has an impact on the response of the output powers, which is represented in Figures 11-13. Load changing at 5 s causes the transient output power. The steadystate output powers of DGs shown in these figures are consistent with the Lagrange multiplier method. It can be observed that fast convergence speed results in a fast response of the output powers. The communication topology of full connection provides the fastest power responses compared to the line and ring connections. With the proposed diffusion-based controller, the responses of output powers are much quicker than the conventional consensus method, particularly in the case of the line connection. The convergence speed of the ICs has an impact on the response of the output powers, which is represented in Figures 11-13. Load changing at 5 s causes the transient output power. The steady-state output powers of DGs shown in these figures are consistent with the Lagrange multiplier method. It can be observed that fast convergence speed results in a fast response of the output powers. The communication topology of full connection provides the fastest power responses compared to the line and ring connections. With the proposed diffusion-based controller, the responses of output powers are much quicker than the conventional consensus method, particularly in the case of the line connection.
The system frequency under the two distributed controllers is shown in Figure 14. With the conventional consensus controller, there is an existence of a frequency drop below the nominal frequency. When the load is increased, the frequency drops to 59.9 Hz. However, with the proposed diffusion controller, the frequency is maintained stably at 60 Hz. When additional load connects to the MG system, the system frequency temporarily shows a slight reduction, and then it gradually recovers to the recorded value of 5.58 s.   The system frequency under the two distributed controllers is shown in Figure 14. With the conventional consensus controller, there is an existence of a frequency drop below the nominal frequency. When the load is increased, the frequency drops to 59.9 Hz. However, with the proposed diffusion controller, the frequency is maintained stably at 60 Hz. When additional load connects to the MG system, the system frequency temporarily shows a slight reduction, and then it gradually recovers to the recorded value of 5.58 s.   The system frequency under the two distributed controllers is shown in Figure 14. With the conventional consensus controller, there is an existence of a frequency drop below the nominal frequency. When the load is increased, the frequency drops to 59.9 Hz. However, with the proposed diffusion controller, the frequency is maintained stably at 60 Hz. When additional load connects to the MG system, the system frequency temporarily shows a slight reduction, and then it gradually recovers to the recorded value of 5.58 s.   The system frequency under the two distributed controllers is shown in Figure 14. With the conventional consensus controller, there is an existence of a frequency drop below the nominal frequency. When the load is increased, the frequency drops to 59.9 Hz. However, with the proposed diffusion controller, the frequency is maintained stably at 60 Hz. When additional load connects to the MG system, the system frequency temporarily shows a slight reduction, and then it gradually recovers to the recorded value of 5.58 s.  Figure 15 shows the cost comparison study with the secondary frequency control schemes in [30,31]. With the proposed control strategies in [30,31], all DGs equally supply power to restore the system frequency, which leads to a higher operation cost of the MG system. By comparison, the proposed strategy manages different power generations among DGs to minimize the operation cost while restoring the system frequency to the nominal value. It can be seen that the operation cost of the MG system with the proposed diffusion controller is smaller than that of the existing controller in [30,31].

Conclusions
This paper has proposed a distributed control of the MG system based on the diffusion strategy. The proposed controller not only optimized the economic operation of the MG system but also maintained the system frequency at the nominal value. By considering the secondary frequency regulation, the load prediction could be avoided in the proposed distributed control since the power balance constraint is kept by regulating the system frequency. The proposed method is in the form of a diffusion algorithm added to the existing primary control system. The proposed strategy was tested in three communication topologies and was compared to the conventional method. The proposed diffusion-based controller showed the advantage of a fast convergence speed compared to the conventional consensus method, particularly in the case of line connection, which has a positive impact on the power response of DGs. The proposed method can be considered as a promising solution for the distributed operation of the MG system.   Figure 15 shows the cost comparison study with the secondary frequency control schemes in [30,31]. With the proposed control strategies in [30,31], all DGs equally supply power to restore the system frequency, which leads to a higher operation cost of the MG system. By comparison, the proposed strategy manages different power generations among DGs to minimize the operation cost while restoring the system frequency to the nominal value. It can be seen that the operation cost of the MG system with the proposed diffusion controller is smaller than that of the existing controller in [30,31].  Figure 15 shows the cost comparison study with the secondary frequency control schemes in [30,31]. With the proposed control strategies in [30,31], all DGs equally supply power to restore the system frequency, which leads to a higher operation cost of the MG system. By comparison, the proposed strategy manages different power generations among DGs to minimize the operation cost while restoring the system frequency to the nominal value. It can be seen that the operation cost of the MG system with the proposed diffusion controller is smaller than that of the existing controller in [30,31].

Conclusions
This paper has proposed a distributed control of the MG system based on the diffusion strategy. The proposed controller not only optimized the economic operation of the MG system but also maintained the system frequency at the nominal value. By considering the secondary frequency regulation, the load prediction could be avoided in the proposed distributed control since the power balance constraint is kept by regulating the system frequency. The proposed method is in the form of a diffusion algorithm added to the existing primary control system. The proposed strategy was tested in three communication topologies and was compared to the conventional method. The proposed diffusion-based controller showed the advantage of a fast convergence speed compared to the conventional consensus method, particularly in the case of line connection, which has a positive impact on the power response of DGs. The proposed method can be considered as a promising solution for the distributed operation of the MG system.

Conclusions
This paper has proposed a distributed control of the MG system based on the diffusion strategy. The proposed controller not only optimized the economic operation of the MG system but also maintained the system frequency at the nominal value. By considering the secondary frequency regulation, the load prediction could be avoided in the proposed distributed control since the power balance constraint is kept by regulating the system frequency. The proposed method is in the form of a diffusion algorithm added to the existing primary control system. The proposed strategy was tested in three communication topologies and was compared to the conventional method. The proposed diffusion-based controller showed the advantage of a fast convergence speed compared to the conventional consensus method, particularly in the case of line connection, which has a positive impact on the power response of DGs. The proposed method can be considered as a promising solution for the distributed operation of the MG system.