Research on a Control Strategy for a Distributed Economic Dispatch System in an Active Distribution Network, Considering Communication Packet Loss

: The beneﬁts of a distributed economic dispatch system include signiﬁcant resilience, uniform processing and communication pressure, and good scalability, but, currently, there is limited in-depth research on the impact of packet loss on the communication links of a distributed economic dispatch system. To analyze how communication packet loss affects a distributed economic dispatch system and to develop appropriate scheduling approaches, in this paper, we design a fully distributed economic dispatch strategy for an active distribution network (ADN) in leader-follower mode, we use the Bernoulli stochastic process to describe the communication packet loss phenomenon, and we give sufﬁcient conditions for the system in the form of linear matrix inequality (LMI) to be able to achieve stochastic stability. The communication packet loss probability affects the convergence performance of the system, and a high packet loss probability can cause the system to lose stability. By simulating various situations in the IEEE-14 and IEEE-39 systems, in this study, we conﬁrm the validity and efﬁcacy of the distributed economic dispatch control strategy and the system stability analysis method.


Introduction
One of the most fundamental problems of power systems is economic dispatch (ED), which involves determining the output of multiple generating units to meet the total load demand and other constrains, while minimizing the total generation cost. This issue has attracted a lot of attention and impressive results have been achieved in terms of dispatch models and optimization algorithms [1,2]. Due to the increasing complexity and scale of power systems and the rapid development of distributed generations (DGs) with high-permeability distributed power, the economic dispatch of a traditional, centralized distribution network cannot meet the requirements of high system reliability and enough topology flexibility, and there is significant pressure on the dispatch center for calculations and communication [3]. An ADN is a distribution network with energy sources that are internally dispersed or decentralized, as well as active control and operational capabilities. ED in an ADN is a subject that has been studied by many scholars. At present, distribution systems have gained increasing attention in ADN economic dispatch systems due to good scalability, uniform computational and communication pressure, and greater robustness [4][5][6].
A study by [7] demonstrated that an optimal solution to the economic dispatch problem could be obtained when the cost of each generation unit increased at a consistent rate. A consensus algorithm was designed, which could solve the economic dispatch problem in a distributed way on the premise of connected communication topology. The method was independent of the initial power distribution and could automatically respond to real-time load demand changes. A study by [8] proposed a leader-follower consensus distributed economic dispatch method, with the incremental cost of each generator as the consensus function and effectively suppressed the delay and noise of information transmission between agents. However, the model did not solve the problem of obtaining the power imbalance value and the dependence on the leader performance, and no requirements were provided for the convergence of the associated management strategy.
As a result, considerable research has been carried out on distributed economic dispatch algorithms based on the consensus protocol, but, so far, there is no in-depth quantitative analysis of the impact of packet loss in the communication links of economic dispatch systems. To examine how a distributed economic dispatch system is impacted by communication packet loss and to provide suitable scheduling strategies, in this study, a distributed economic dispatch control strategy for an active distribution network in leader-follower mode is designed based on the stochastic theory using the incremental cost as the consensus variable and the Bernoulli stochastic process to describe the packet loss phenomenon, considering the coupling of the power difference between the total generated power and the load power with the consistency variable, while using the Lyapunov stability theorem and the Schur complementary theorem. Sufficient conditions are given for the LMI form to enable the system to achieve stochastic stability. Finally, the correctness and validity of the control strategy and sufficient conditions in the form of LMI to enable the system to achieve stochastic stability are verified by simulation.

Distributed Economic Dispatch Architecture
An active distribution network economic dispatch system requires the sum of the total active power output of all the DGs in the distribution network and the active power input from the external network to balance the active power demand of the total distribution network load. It also minimizes the total generation cost of the DGs in the distribution network while satisfying other constraints. In this paper, we adopt an ADN dispatch system, as shown in Figure 1, to achieve the above goals. The ADN is connected to the external network, and the DGs can communicate with each other as well as different loads. The PCC, the common connection point between the external grid and the active distribution grid, monitors the power flow exchange between the two in real time. The power difference between the total generated power and the load power is an important piece of information to enable a distributed economic dispatch control strategy, and, in this paper, we provide this information indirectly by measuring the exchanged power Ppcc at the PCC, which can be measured by the distribution network regulator.
The regulator of each DG makes independent decisions according to the control strategy based on local measurement information, such as the respective incremental cost λi The ADN is connected to the external network, and the DGs can communicate with each other as well as different loads. The PCC, the common connection point between the external grid and the active distribution grid, monitors the power flow exchange between the two in real time. The power difference between the total generated power and the load power is an important piece of information to enable a distributed economic dispatch control strategy, and, in this paper, we provide this information indirectly by measuring the exchanged power P pcc at the PCC, which can be measured by the distribution network regulator.
The regulator of each DG makes independent decisions according to the control strategy based on local measurement information, such as the respective incremental cost λ i and output power P DG,i of the respective DG, combining the information sent by other agents via the communication link and the system power difference information, etc. Then, the regulator sends the power command P DG,i to the DG that it regulates and sends its own information to the other DG regulators, so that the active power balance of the distribution network is achieved through the cooperation of the DGs and the incremental cost of each DG is equal in the end, thus, realizing the distributed economic dispatch control strategy.
In order to facilitate the incorporation of the power exchange information between the distribution network and the external network into the designed distributed economic dispatch control strategy in a uniform form, in this paper, we provide the distribution network regulator with a "virtual distributed power source" and set the number of this "virtual distributed power source" to one. The virtual DG has no actual power output and does not need to send commands to the DG.

Distributed Economic Dispatch Mode
The objective function for minimizing the total cost of generation is: where n is the total number of DGs in the distribution network; P DG,i and F(P DG,i ) are the output power and power-cost function of the ith DG, respectively. Distributed energy sources such as wind turbines and photovoltaics can be formed into dispatchable DGs by adding energy storage devices, whose power-cost functions are mostly in the form of quadratic convex functions, i.e.,: where a i , b i , and c i are the quadratic, primary, and constant term coefficients of the i-th DG generation cost function, respectively, and a i > 0. For a virtual DG, its power-cost function is constant at zero and has no real power output. The equation constraint condition of the dispatch model is power balance, namely: where P D is the total active distribution network load, P LOSS is the active power lost of the distribution network, and P PCC is the active power input from the external grid. The linear inequality constraint for the economic dispatch model is: where P min DG,i and P max DG,i are the lower limit and upper limit of output of the i-th DG, respectively.
This constrained economic dispatch problem can be solved using the Lagrange multiplier method, and the Lagrange operator L is specified as follows: where λ is the Lagrange multiplier, which also indicates the generation incremental cost for each DG; the equation constrains ∆P = P D + P LOSS − ∑ n i=1 P DG,i and indicates the power imbalance difference for this power system.
Take the partial derivatives of the economic scheduling function L(P DG,i , λ) with respect to P DG,i and λ and set them equal to zero, then, the optimal solution is solved as follows [10]: where P is the total network power demand in the ADN.

Packet Loss Modeling for Multi-Intelligence Communications
In this paper, we use a distributed consensus protocol to complete the exchange of information between the DG agents.
The topology of a multi-intelligence network containing n agents can be represented by a graph G = (V, E), where V = {1,2, . . . , n} is the set of vertices; E ⊆ V × V is the set of edges; and the ordered node pair (i, j)∈E denotes a directed edge starting at node i and ending at node j, i.e., agent j can receive messages from agent i. Define N + i = {j∈V:(i, j)∈E} and N − i = {j∈V:(j, i)∈E} as the set of information output and input neighbor nodes, respectively. For an undirected graph N + i = N − i , the set of neighbors of vertex i is uniformly denoted by N i = {j∈V:(i,j)∈E} and |N i | denotes the number of elements in the set, also known as the degree of node i, denoted as D i .
For a first-order linear system, the dynamics of the multi-agent system consensus protocol (control strategy of node i) can be expressed as: where ω is the feedback factor; N i (k) denotes the input type domain of agent i at moment k, agent i receives information from the agent in that domain at moment k; x i (k) is the state of the i-th agent at moment k; a ij > 0 is the weight set by node i for the state information received from neighboring node j. In this paper, we select the following scheme for setting information weights a ij : The normal operation of a multi-agent system depends on the state of the communication network. In practice, the communication links between the agents suffer from packet loss, delays, noise, and restricted bandwidth.
In this paper, we conduct research on the communication packet loss situation. There are currently three control design approaches for communication packet loss, namely the zero-setting method, the zero-order retainer method, and the predictive control method [20]. When a communication packet loss occurs: The zero-setting method sets the information received on that communication line to zero; using the zero-order retainer method, the information received from that communication line remains at the last received value until new data are received; and with the predictive control method, the agent uses a type of predictive algorithm to generate data to replace the information received from that communication line.
Since the zero-setting method has the advantages of simple implementation and system analysis, no loss of communication efficiency, and small information storage space requirements, in this paper, we chose to use the zero-setting method.
The communication between intelligences is propagated through communication links, and the information exchange between different DGs is often independent of each other. In distributed power systems, external information attacks, aging failure of communication links, or interface instability can lead to communication packet loss, which can be regarded as a random occurrence; therefore, we consider a multi-intelligent system under asynchronous communication packet loss, where the packet loss of all communication lines is a random event with independent identical distribution (i.i.d.), and we adopt a Bernoulli distribution communication packet loss model based on the zero-setting method.
Taking agent j sending a message to agent i as an example, the stability of the communication line v j , v i is represented by γ ij . When the line is stable, γ ij = 1 and, when the line Electronics 2022, 11, 3288 6 of 18 experiences packet loss, γ ij = 0. γ ij uses a Bernoulli distribution with an expectation of p ij , i.e.,: In this paper, assuming that the model shown in Equation (11) is used for packet loss on each communication channel of the communication network, Equation (9) can be rewritten under the zero-setting method as: Assumption 1. For any pair of agents i and j, the communication channel between them exists or disappears at the same time.
Assumption 1 guarantees that the communication topology is always biaxially symmetric, in a dynamic process such that the Laplace matrix L(k) of the system is a symmetric matrix.
The elements of the system Laplacian matrix L(k) = l ij (k) ∈ R n×n corresponding to Equation (12) are:

Distributed Economic Dispatch Strategy
With the incremental cost λ of the distributed power as the consensus variable, the distributed economic dispatch strategy considering the virtual DG is: (14) where λ i (k) and u i (k) are the incremental cost and control input of DG i at moment k, respectively; d i is the leader identifier, d i = 1 means that this virtual DG is the leader and for other distributed power supplies there is d i = 0; ω takes values in the range 0 < ω <1/∆ and ∆ = max{|N 1 |, |N 2 |, · · · , |N n |} is the maximum degree of the graph; ε > 0 is the power deviation factor. In general, too large a value of ε will lead to non-convergence and too small will lead to slow convergence.
Equation (14) is written in matrix form as follows: The control strategy shown in Equation (15) consists of two components: (1) Regarding the share component, λ(k + 1) = (I − ωL)λ(k) can reduce the difference between each DG incremental cost λ and achieve an equal incremental rate for each generation unit through multiple iterations. (2) Regarding the total adjustment component, the diag(d i )ε∆ f p (k) signal received by the leader node is used to adjust the total generation of each DG to satisfy the limit of zero power deviation.
In this paper, the economic dispatch performance evaluation function f p is defined as [12]: Take ∆ f p (k) to be the negative value of the derivative of the performance evaluation function f p with respect to the leader's output force, i.e.,: Based on the updated value of the incremental cost obtained in Equation (15), each generating unit can obtain an instruction for its own output: Considering the power constraint for each DG output, Equation (18) can be changed to: Consider ∆ f p (k) is coupled to λ, and substitute Equation (19) into Equation (15), Equation (15) can be rewritten in the following matrix form: is the leader selection matrix after coupling and K = [k 1 , k 2 , . . . , k n ] is the system power difference coupling matrix, The distributed economic dispatch control strategy starts by initializing ∆ f p (k) and λ 0 , then updates the incremental cost of each DG after the iteration according to Equation (14), and then obtains the output of each DG at that time by calculating the incremental cost of each DG obtained, and finally substitutes the updated λ and the output of each DG into Equation (14) to continue the iteration until the system converges.

System Stability Analysis
The stability analysis of a distributed economic dispatch system under the probabilistic communication packet loss model is a stochastic stability analysis problem, and the necessary reasoning is first given below.
Lemma 1 [21]. For symmetric positive definite matrices P ∈ R n×n and matrices E ∈ R r×n , F ∈ R r×n , there are the following relations: The following three conditions are equivalent: Lemma 3 [22]. For a symmetric graph G with a given Laplacian matrix L(k) and symmetric matrix P, the expectation E{L(k)PL(k)} can be derived from the following equation: where L (0) is the Laplacian matrix of the graph G corresponding to no consideration of communication packet loss, and Ξ(P) is a function of the matrix P defined as follows: where the matrix E (m,q) ∈ R n×n is 1 only at (q, m) and 0 at all other positions. Each E T (m1,q1) PE (m2,q2) will produce a new matrix. This matrix has a unique non-zero element in row m1 and column m2 only, and the value of this non-zero element is equal to P (q1, q2) .
The Lyapunov function is designed as V(k) = Z T (k)PZ(k) and Z(k + 1) = G(k)Z(k). P is a real symmetric positive definite matrix, and the mathematical expectation of the difference between the functions V(k+1) and V(k) is: where N = I − εM.
The mathematical expectation of the matrix L(k) can be calculated from E{L(k)} = pL (0) , from which the expectation of the difference of Lyapunov functions can be obtained as: The mathematical expectation of the difference between the functions V(k + 1) and V(k), E{∆V(k)} < 0, is equivalent to the matrix inequality: By treating −L (0) as E T and N as F in Lemma 1, we can obtain: Substituting Equation (27) into Equation (26), obtain: Then, a sufficient condition for Equation (26) to hold is obtained as: Repeated use of Lemma 2 for Equation (29) yields: where Θ = p(1 − p)Ξ(P), ϑ = sym L (0) PN , = 1 2 ωp + 1N T P, and Ψ = 1 2 ωp + ω 2 p 2 L (0) P. Combining the above equations, if for given positive real numbers ε, ω, and p, there exists a real symmetric positive definite matrix P such that Equation (30) holds, it follows from the Lyapunov stability theory that the system Z(k+1) = G(k)Z(k) is stochastically asymptotically stable, and hence, the distributed economic dispatch system Equation (20) is stochastically asymptotically stable.
According to the LMI in Equation (30), when p → 0, Equation (30) is approximated by the following linear matrix inequality: It is verified in MATLAB that Equation (31) is not strict (marginal infeasibility), and when the communication packet loss probability is high, i.e., p → 0 , the system cannot remain stable. In addition, as the communication packet loss probability increases, the efficiency of the information interaction between the DGs decreases and the equalization part of the control strategy (15) reduces the difference between the incremental costs λ of the DGs slower, which can lead to the convergence performance of the system (20) being affected.

Simulation Verification
The IEEE-14 node distribution network system was used for simulation calculations. The system structure and local communication network are shown in Figure 2 [8]. Five DGs are connected to the distribution network, located at nodes 2, 3, 7, 10, and 12. The total active load of the entire network is 28.7 MW and the specific information for each DG is given in Table 1. The adjudication condition for convergence of the system in the simulation is that the difference between any two DG incremental costs |∆λ| < 0.0001 $/(MWh).

Scenario 1: Ideal Communication Conditions
Taking the feedback coefficient ω = 0.2 and the power deviation coefficient ε = 0.1, the eigenvalues of the control system shown in Equation (20) are 0.7611 ± j 0.3707, 0.0593, 0.3208, 0.6978, and 0.6000, all within the unit circle, indicating that the system is stable.
The total active load of the distribution network increases from an initial 28.7 MW to 33.7 MW at time t = 0.7 s. The response of the distributed economic dispatch system is

Scenario 1: Ideal Communication Conditions
Taking the feedback coefficient ω = 0.2 and the power deviation coefficient ε = 0.1, the eigenvalues of the control system shown in Equation (20) are 0.7611 ± j 0.3707, 0.0593, 0.3208, 0.6978, and 0.6000, all within the unit circle, indicating that the system is stable.
The total active load of the distribution network increases from an initial 28.7 MW to 33.7 MW at time t = 0.7 s. The response of the distributed economic dispatch system is shown in Figure 3.

Scenario 1: Ideal Communication Conditions
Taking the feedback coefficient ω = 0.2 and the power deviation coefficient ε = 0.1, the eigenvalues of the control system shown in Equation (20) are 0.7611 ± j 0.3707, 0.0593, 0.3208, 0.6978, and 0.6000, all within the unit circle, indicating that the system is stable.
The total active load of the distribution network increases from an initial 28.7 MW to 33.7 MW at time t = 0.7 s. The response of the distributed economic dispatch system is shown in Figure 3.  As can be seen from Figure 3, when the load is 28.7 MW, the DG incremental cost converges to the same optimal value λ* = 1.497$/(MWh) from each initial value, the sum of each DG output is balanced with the load in the distribution network; when the load As can be seen from Figure 3, when the load is 28.7 MW, the DG incremental cost converges to the same optimal value λ* = 1.497$/(MWh) from each initial value, the sum of each DG output is balanced with the load in the distribution network; when the load increases to 33.7 MW, the DG output fluctuates briefly, and the DG output meets the new load demand. As the total load power increases, the incremental cost of each DG stabilization increases to the new optimal consensus value λ* = 1.667 $/(MWh). It is shown that the method in this paper can achieve the objective of economic dispatch of the active distribution network in a distributed manner during load changes, when the system converges to the optimal value in 0.63 s and 0.55 s

Scenario 2: Communication Packet Loss
(1) Setting p = 0.9, i.e., packet loss probability is 0.1, Equation (29) holds and the system is able to maintain stability. The total active load of the distribution network increases from the initial 28.7 MW to 33.7 MW at t = 1 s. Simulations are performed using the control strategies of [4] and this paper, respectively, and the responses of the distributed economic scheduling system are shown in Figures 4 and 5.

Scenario 2: Communication Packet Loss
(1) Setting p = 0.9, i.e., packet loss probability is 0.1, Equation (29) holds and the system is able to maintain stability. The total active load of the distribution network increases from the initial 28.7 MW to 33.7 MW at t = 1 s. Simulations are performed using the control strategies of [4] and this paper, respectively, and the responses of the distributed economic scheduling system are shown in Figures 4 and 5.

Scenario 2: Communication Packet Loss
(1) Setting p = 0.9, i.e., packet loss probability is 0.1, Equation (29) holds and the system is able to maintain stability. The total active load of the distribution network increases from the initial 28.7 MW to 33.7 MW at t = 1 s. Simulations are performed using the control strategies of [4] and this paper, respectively, and the responses of the distributed economic scheduling system are shown in Figures 4 and 5.  From the comparison of Figures 4 and 5a, the control strategy of [4] is not applicable at all when there is communication packet loss, and the packet loss condition will have a wide range of fluctuations and even break the steady state that has been achieved, while the control strategy of this paper can adapt well to the communication packet loss condition and achieve the control goal. As can be seen from Figure 5, for a load of 28.7 MW and a packet loss probability of 0.1, the DG incremental costs converge from their respective initial values to the same optimal value of λ* = 1.497 ($/MWh), and the sum of each DG's output is balanced with the load in the distribution network. The output of each DG briefly fluctuates when the demand reaches 33.7 MW, but soon stabilizes at a new ideal consensus value of λ* = 1.667 ($/MWh).
It is shown that the system (20) considering communication packet loss can achieve the goal of economic dispatch of an active distribution network in a distributed manner, but due to the effect of communication packet loss, each DG will not be able to accurately receive the increment cost from neighboring DGs or pass out relevant information about itself in a timely manner, and the time to converge to the optimal value before and after the load change increases to 0.65 s and 0.62 s, respectively.
(2) Holding other parameters constant, the response of the system at p = 0.7 and p = 0.5 is shown in Figure 6.
output is balanced with the load in the distribution network. The output of each DG briefly fluctuates when the demand reaches 33.7 MW, but soon stabilizes at a new ideal consensus value of λ* = 1.667 ($ /MWh). It is shown that the system (20) considering communication packet loss can achieve the goal of economic dispatch of an active distribution network in a distributed manner, but due to the effect of communication packet loss, each DG will not be able to accurately receive the increment cost from neighboring DGs or pass out relevant information about itself in a timely manner, and the time to converge to the optimal value before and after the load change increases to 0.65 s and 0.62 s, respectively.
(2) Holding other parameters constant, the response of the system at p = 0.7 and p = 0.5 is shown in Figure 6. As can be seen from Figure 6, after the load increases to 33.7 MW, the incremental costs eventually stabilize to a new optimal consistent value of λ* = 1.667 ($ /MWh), and the convergence time of the system to the final consistent value is 0.67 s, 0.92 s, and 1.09 s when p is 0.9, 0.7, and 0.5 respectively, indicating that the convergence speed of the distributed economic scheduling system becomes slower and slower as the communication packet loss probability increases. In addition, the variation in the incremental cost of each DG as the load changes increases with the probability of packet loss.
(3) For an initial load of 28.7 MW, the time for the system to converge from the initial value to the optimum value for different packet loss probabilities is shown in Figure  7.
As can be seen from Figure 7, the time required for the system to converge to the optimum value increases rapidly with the probability of packet loss. In the communication topology shown in Figure 2, the system convergence time is <5 s when p > 0.2, i.e., the packet loss probability is less than 0.8; when the packet loss probability is higher than 0.8, As can be seen from Figure 6, after the load increases to 33.7 MW, the incremental costs eventually stabilize to a new optimal consistent value of λ* = 1.667 ($ /MWh), and the convergence time of the system to the final consistent value is 0.67 s, 0.92 s, and 1.09 s when p is 0.9, 0.7, and 0.5 respectively, indicating that the convergence speed of the distributed economic scheduling system becomes slower and slower as the communication packet loss probability increases. In addition, the variation in the incremental cost of each DG as the load changes increases with the probability of packet loss.
(3) For an initial load of 28.7 MW, the time for the system to converge from the initial value to the optimum value for different packet loss probabilities is shown in Figure 7.
As can be seen from Figure 7, the time required for the system to converge to the optimum value increases rapidly with the probability of packet loss. In the communication topology shown in Figure 2, the system convergence time is <5 s when p > 0.2, i.e., the packet loss probability is less than 0.8; when the packet loss probability is higher than 0.8, the system convergence time increases significantly; when the packet loss probability exceeds 0.9, the convergence time is already close to 50 s and the system can no longer converge. The maximum packet loss probability p max = 0.904, which holds by solving Equation (29) through the LMI toolbox in MATLAB, is consistent with the simulation results, verifying the validity of the sufficient condition (29) for the stochastic asymptotic stability of system (20). the system convergence time increases significantly; when the packet loss probability exceeds 0.9, the convergence time is already close to 50 s and the system can no longer converge. The maximum packet loss probability pmax = 0.904, which holds by solving Equation (29) through the LMI toolbox in MATLAB, is consistent with the simulation results, verifying the validity of the sufficient condition (29) for the stochastic asymptotic stability of system (20).

Scenario 3: IEEE-39 Node System
To further verify the effectiveness of the strategy proposed in this paper, a larger scale IEEE-39 node 10 DG system [23] was used for simulation analysis. Table 2 lists the basic parameters of each DG, G1 is a virtual distributed power supply, and the system communication network topology is shown in Figure 8

Scenario 3: IEEE-39 Node System
To further verify the effectiveness of the strategy proposed in this paper, a larger scale IEEE-39 node 10 DG system [23] was used for simulation analysis. Table 2 lists the basic parameters of each DG, G 1 is a virtual distributed power supply, and the system communication network topology is shown in Figure 8 G7 G8 Figure 8. Communication network topology of the IEEE-39 system.
(1) At p = 0.9 and t = 2 s, the system load increases by 200 MW from an initial 3420 MW, and decreases by 250 MW at t = 4 s. The response of the distributed economic dispatch system is shown in Figure 9. (1) At p = 0.9 and t = 2 s, the system load increases by 200 MW from an initial 3420 MW, and decreases by 250 MW at t = 4 s. The response of the distributed economic dispatch system is shown in Figure 9. (1) At p = 0.9 and t = 2 s, the system load increases by 200 MW from an initial 3420 MW, and decreases by 250 MW at t = 4 s. The response of the distributed economic dispatch system is shown in Figure 9.  Figure 9. Incremental cost at p = 0.9.
As can be seen from Figure 9, when the load is 3420 MW and there is a packet loss probability of 0.1, the DG incremental costs converge to the same optimal value λ* = 17.903 $/(MWh) from their respective initial values, and the sum of each DG output is balanced with the load in the distribution network. When the load fluctuates at t = 2 s and t = 4 s, there are brief fluctuations in each DG output, and the incremental costs stabilize to the new optimal consistent values λ* = 18.9075 $/(MWh) and λ* = 17.6519 $/(MWh). It is further shown that system (20) can achieve the goal of economic dispatch of the active distribution network in the event of packet loss in communication.
(2) Regarding the impact of network connectivity, the distributed control systems can not achieve control objectives when the communication network is not connected, and communication packet loss will affect network connectivity to some extent. The degrees of each DG node in Figure 8 are: 3, 4, 5, 4, 3, 5, 5, 3, 3, 4, and 3, where the size of D1 of the leader node has a greater impact on the convergence performance of the system. The convergence times of the system from the initial value to the optimal value for different values of p and the number of connected communication links between the leader and other DGs are shown in Table 3 and Figure 10. As can be seen from Figure 9, when the load is 3420 MW and there is a packet loss probability of 0.1, the DG incremental costs converge to the same optimal value λ* = 17.903 $/(MWh) from their respective initial values, and the sum of each DG output is balanced with the load in the distribution network. When the load fluctuates at t = 2 s and t = 4 s, there are brief fluctuations in each DG output, and the incremental costs stabilize to the new optimal consistent values λ* = 18.9075 $/(MWh) and λ* = 17.6519 $/(MWh). It is further shown that system (20) can achieve the goal of economic dispatch of the active distribution network in the event of packet loss in communication.
(2) Regarding the impact of network connectivity, the distributed control systems can not achieve control objectives when the communication network is not connected, and communication packet loss will affect network connectivity to some extent. The degrees of each DG node in Figure 8 are: 3, 4, 5, 4, 3, 5, 5, 3, 3, 4, and 3, where the size of D 1 of the leader node has a greater impact on the convergence performance of the system. The convergence times of the system from the initial value to the optimal value for different values of p and the number of connected communication links between the leader and other DGs are shown in Table 3 and Figure 10.  As can be seen from Table 3 and Figure 10, for the same packet loss probability, the system converges slower for both D1 = 1 and D1 = 3 than for D1 = 5, and the smaller the value of D1, the slower the system converges. The main reason is that, when the degree of the leader node is small, the DGs cannot quickly obtain information about the power difference of the system through the leader node, and thus, quickly adjust their own output. When the degree of the leader increases, the convergence speeds up. However, the increase in D1 leads to an increase in communication costs; therefore, it is important to choose the right D1 value to ensure the quality of communication and convergence performance while ensuring economic efficiency.
(3) Regarding individual DG communication failure, let DG3, DG6, and DG9 each have communication failures and packet loss with their respective adjacent DGs. The system's convergence time is illustrated in Figure 11 below, and shows that the communication packet loss probability increases between the failing DG and the neighboring DGs.  As can be seen from Table 3 and Figure 10, for the same packet loss probability, the system converges slower for both D 1 = 1 and D 1 = 3 than for D 1 = 5, and the smaller the value of D 1 , the slower the system converges. The main reason is that, when the degree of the leader node is small, the DGs cannot quickly obtain information about the power difference of the system through the leader node, and thus, quickly adjust their own output. When the degree of the leader increases, the convergence speeds up. However, the increase in D 1 leads to an increase in communication costs; therefore, it is important to choose the right D 1 value to ensure the quality of communication and convergence performance while ensuring economic efficiency.
(3) Regarding individual DG communication failure, let DG 3 , DG 6 , and DG 9 each have communication failures and packet loss with their respective adjacent DGs. The system's convergence time is illustrated in Figure 11 below, and shows that the communication packet loss probability increases between the failing DG and the neighboring DGs.  As can be seen from Table 3 and Figure 10, for the same packet loss probability, the system converges slower for both D1 = 1 and D1 = 3 than for D1 = 5, and the smaller the value of D1, the slower the system converges. The main reason is that, when the degree of the leader node is small, the DGs cannot quickly obtain information about the power difference of the system through the leader node, and thus, quickly adjust their own output. When the degree of the leader increases, the convergence speeds up. However, the increase in D1 leads to an increase in communication costs; therefore, it is important to choose the right D1 value to ensure the quality of communication and convergence performance while ensuring economic efficiency.
(3) Regarding individual DG communication failure, let DG3, DG6, and DG9 each have communication failures and packet loss with their respective adjacent DGs. The system's convergence time is illustrated in Figure 11 below, and shows that the communication packet loss probability increases between the failing DG and the neighboring DGs.  Figure 11. System convergence time for DG 3 , DG 6 , and DG 9 failures.
As can be seen from Figure 11, the system convergence performance is most affected when DG 3 experiences a communication failure.

Conclusions
With incremental cost as the consensus variable, the Bernoulli's stochastic packet loss process is used to describe the packet loss phenomenon, and based on the stochastic theory, a distributed economic dispatch strategy for the ADN in the leader-follower mode is designed, considering the coupling of the power difference between the total generation power and the load power with the consensus variable, while using Lyapunov's stability theorem and the Schur complementary theorem. Sufficient conditions are given in the form of LMI that enable the system to achieve stochastic stability, which were not considered and are given in the references. Simulation calculations were carried out for several scenarios in different-sized ADNs and the results showed that: The distributed economic dispatch control strategy based on communication packet loss can effectively meet the requirements of ADN power balance and minimum generation cost at the same time, and the ADN can maintain stable operation, which improves the practicality and reliability of the distributed economic dispatch control strategy of the ADN.
(2) The validity of the sufficient conditions for the stochastic stability of the system in the form of LMI is verified by Scenario 2, and the convergence of the system under different communication packet loss probabilities can be judged by Equation (29) and the maximum packet loss probability p max of the system can be obtained. (3) The communication packet loss probability will affect the convergence performance of the system, and too high a packet loss probability will cause the system to lose stability; the selection of degree D affects the convergence performance of the system, and a suitable degree should be selected to ensure that each DG is not disturbed by too much information, while it collects enough information from neighboring DGs.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

ADN
Active distribution network LMI Linear matrix inequality ED Economic dispatch DG distributed generation PCC Point of common coupling i, j Indices of the DGs n Total number of DGs ω Feedback factor in consensus protocol P DG,i Output power of the i-th DG v j Node j F(P DG,i ) Power-cost function of the ith DG a i , b i , c i Quadratic, primary, and constant term coefficients of the i-th DG generation cost function P D Total active distribution network load P LOSS Active power lost of the distribution network P PCC Active power input from the external grid P min DG,i , P max DG,i Lower limit and upper limit of output of the i-th DG G A gragh V Set of vertices E Set of edges N i (k) Input type domain of agent i at moment k x i (k) State of the ith agent at moment k