Next Article in Journal
The Effect of Oxygen Tension on the Differentiation of Outgrowth Cells from Embryoid Bodies Produced by Mouse Induced Pluripotent Stem Cells
Next Article in Special Issue
Multitarget Search Algorithm Using Swarm Robots in an Unknown 3D Mountain Environment
Previous Article in Journal
Sunlight Bleaching of Subporphyrazine Dye Films
Previous Article in Special Issue
Heat Load Forecasting of Marine Diesel Engine Based on Long Short-Term Memory Network
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissipativity Analysis of Large-Scale Networked Systems

1
School of Automation, Qingdao University, Qingdao 266071, China
2
School of Mechanical and Electrical Engineering, Qingdao University, Qingdao 266071, China
3
Weihai Innovation Insitute, Qingdao University, Weihai 264200, China
4
Shandong Key Laboratory of Industrial Control Technology, Qingdao 266071, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(2), 1214; https://doi.org/10.3390/app13021214
Submission received: 18 December 2022 / Revised: 13 January 2023 / Accepted: 13 January 2023 / Published: 16 January 2023
(This article belongs to the Special Issue Fault Detection and State Estimation in Automatic Control)

Abstract

:
This paper investigates the dissipativity analysis of large-scale networked systems with linear time-invariant dynamics. The networked system is composed of a large number of subsystems whose connections are arbitrary, and each subsystem can have different dynamics. A sufficient and necessary condition for the strict dissipativity analysis of the networked system is derived, which takes advantage of the block-diagonal structure of the system parameter matrix and the sparseness characteristics of the subsystem interconnections. Then, a necessary condition and a sufficient condition that depend only on a single subsystem parameter are given separately. Numerical simulations illustrate that compared with the existing results, the conditions suggested in this paper have higher computational efficiency in the dissipative analysis of large-scale networked systems.

1. Introduction

In recent years, the research of large-scale networked systems has attracted great attentions [1,2,3,4]. The system can be considered as composed of a large number of subsystems with different spatial locations connected in a certain way [5]. Generally, subsystems exchange information with their neighbors straightforwardly and predictably, but the system often exhibits complicated dynamic behavior when seen as a whole. Such systems have extensive application background, including airplane formation flight [6], power network distributed system [7], automated highways [8], multi-agent formation systems [9], and so on. For such a complex system, the classic method of bringing all the subsystems together and analyzing it as a single large-scale system has very strict requirements on the calculation speed and accuracy of the computer, which will inevitably bring computational difficulties. Therefore, using the system structure to find more efficient computational conditions is of great engineering significance for the dissipativity analysis of large-scale networked systems.
There are already many results on the performance analysis of networked systems, but the research on dissipativity is not mature enough and needs to be further developed and improved. In 1972, the famous scholar Willems first put forward the concept of dissipativity [10,11]. Dissipativity describes the equilibrium relationship between the system’s internal energy, which is a vital concept in theoretical research and practical application. Its essential meaning is that there is a non-negative energy function (called storage function) so that the energy supply rate of the system is always greater than the loss of energy inside the dynamic system. Based on Willems’ work, many scholars have done in-depth studies on dissipativity and obtained rich results, which have played a significant role in the field of circuit, system, and control theory. Refs. [12,13] respectively studied the dissipative control problems of linear continuous-time and discrete systems based on linear matrix inequality (LMI) methods. A simplified mathematical model of the interconnected two-machine power system was established in [14], and its non-linear dynamic behavior such as dissipativity properties was analyzed. In [15], the analysis and improvement of the dissipativity performance of interconnected passive systems are studied. For networked control systems, Ref. [16] obtained some new sufficient conditions by utilizing Lyapunov stability theory and LMI technology to ensure that the closed-loop system is finite-time limited and dissipative. In [17], a distributed controller was created to ensure the dissipativity of a networked system made up of dynamically coupled subsystems. Its control synthesis is done locally at the subsystem level and doesn’t involve the relationship among subsystems, hence it has certain drawbacks. The linear dynamic system with the interconnection structure specified by the directed graph is studied in [18]. Based on the dissipativity inequality, an LMI for calculating system performance is established and the concept of local dissipativity is defined. Using the knowledge of graph theory to analyze large-scale networked systems has certain constraints on the dynamic characteristics and connection modes of subsystems, which has certain limitations in practical application [19,20].
Considering that the system structure of large-scale networked systems usually has sparse characteristics or specific structural forms [21,22,23], in the large-scale connected systems discussed in [24], the concept of internal input and output is introduced to represent the connections and functions among subsystems, and the connection relationship among sub-units of the entire system is described by subsystem connection matrix. This description method takes into account the situation where the dynamic characteristics of the subsystems are different and the connection relationship of the subsystems is arbitrary. The previously mentioned UAV formation flight refers to the arrangement of multiple UAVs in a certain formation so that they maintain in formation or change their relative positions within a limited scope during the flight. To maintain a certain formation shape, information interaction is required among the UAVs. In a centralized strategy, each UAV has to know information about the whole formation, demanding substantial information interaction. It is computationally intensive and requires the high performance of the airborne computer. In fact, each UAV can interact with its position, speed, attitude, and motion target with only the UAVs connected to it in the formation. In this way, the amount of computation is greatly reduced, and the system is relatively simple to implement. It is this sparse property or specific structural form among subsystem connections that we exploit to give more computationally efficient dissipativity criteria for large-scale networked systems. Dissipativity explains some of the energy losses and control problems of control systems and is a more general performance indicator of system performance.
In this paper, our objective is to reduce the computational burden of dissipativity analysis for large-scale network systems with a large number of subsystems. In general, the large-scale networked systems studied in this paper have the following characteristics. The first is that the scale of the system is large, including many subsystems. The second and most important point is that the interaction among subsystems of large-scale networked systems is usually sparse or has a specific structural form. We introduce intermediate variables in networked systems to describe the relationship among subsystems, which is more general and explicitly characterizes the structural characteristics of large-scale systems. In this regard, this paper proposes several new LMI conditions, which effectively use the block diagonal structure of the system parameter matrix and the sparsity of the subsystem connection matrix, avoiding the inverse computation of high-dimensional matrices. The simulation results indicate that the conditions proposed in this paper are more efficient than the existing results.
The following is the structure of the paper. The model and the definition of dissipativity for the networked system and some preliminary results are given in Section 2. In Section 3, some conditions for dissipative analysis of networked systems are given, and the relationship between these conditions and existing conditions is discussed. Some numerical simulation results are presented in Section 4. The research results of this paper are summarized in Section 5, and the direction of further research is also proposed here.
Notation 1.
The symbol R is used to denote the set of real numbers, and the vector space produced by real numbers of appropriate dimensions is denoted as R # . c o l Z i | i = 1 L denotes the vector/matrix stacked by Z i ( i = 1 , 2 , , L ) , and d i a g { Z i | i = 1 L } signifies a block diagonal matrix with Z i as the i-th diagonal block. Z i j i = 1 , j = 1 i = M , j = N stands a matrix with M × N blocks, and its i-th row j-th column block matrix being Z i j . 0 n and 0 n × m represent the n dimensional zero column vector and the n × m dimensional zero matrix respectively, the dimension subscript is omitted if there is no ambiguity, and the identity matrix I is the same. The superscript T represents the transpose of a matrix or vector, and ( * ) T W Z or Z W ( * ) T is shorthand for Z T W Z or Z W Z T .

2. System Description and Some Preliminaries

The networked system Γ is consisted of N linear time-invariant subsystems, and the i-th subsystem Γ i is defined by the state-space model below,
x ˙ ( t , i ) z ( t , i ) y ( t , i ) = A TT ( i ) A TS ( i ) B T ( i ) A ST ( i ) A SS ( i ) B S ( i ) C T ( i ) C S ( i ) D T ( i ) x ( t , i ) v ( t , i ) u ( t , i ) ,
in which t and i denote respectively for the temporal variable and the index number of a subsystem, i = 1 , 2 , , N . Moreover, x ( t , i ) is the state vector of the i-th subsystem Γ i at time t. y ( t , i ) and u ( t , i ) represent the external output vector and external input vector of the Γ i , respectively. z ( t , i ) and v ( t , i ) are the output vector to other subsystems and input vector from others, which is also called internal output vector and input vector. The connection relationship among subsystems can be expressed as
v ( t ) = Φ z ( t ) ,
here, v ( t ) = c o l v ( t , i ) i = 1 N and z ( t ) = c o l z ( t , i ) i = 1 N . Φ is called the subsystem connectivity matrix. We assume that each row of the matrix Φ has only one non-zero element equal to one and there are no columns in which all of the items are equal to zero. That means the internal output channels of a subsystem can affect the internal inputs channels of other subsystems, and some subsystem internal input channels depend on the internal output of multiple subsystems. This assumption, as explained in [23], does not jeopardize the generality of the adopted system model. Approximate power-law degree distribution widely exists in engineering systems, such as protein interaction networks, gene regulatory networks, power systems, the Internet, etc. [23]. In these systems, the dimension of the subsystem connection matrix Φ is usually much smaller than the state dimension of the system, and the interactions among subsystems are sparse.
In this paper, we assume that the dimensions of vectors x ( t , i ) , v ( t , i ) , z ( t , i ) , u ( t , i ) and y ( t , i ) are m x i , m v i , m z i , m u i and m y i , respectively. Based on the above assumptions and Equation (2), the dimension of the matrix Φ is i = 1 N m v i × i = 1 N m z i . Then we can get Φ T Φ = Σ 2 , where Σ 2 = diag Σ j 2 j = 1 N , Σ j 2 = diag m ( i ) i = M z , j 1 + 1 M z , j , M z , i = k = 1 i m z k , m ( i ) indicates the number of subsystems directly affected by the i-th element of the vector z ( t ) , i = 1 , , k = 1 N m v k , j = 1 , , N .
To simplify the mathematical derivation, we define the following matrix, A * # = d i a g A * # ( i ) i = 1 N , B * = d i a g B * ( i ) i = 1 N , C * = d i a g C * ( i ) i = 1 N and D * = d i a g D * ( i ) i = 1 N , in which * , # = T , S . By exploiting the connection relationship among subsystems, the dynamic system Γ may be expressed equivalently in the following state-space form,
x ˙ ( t ) y ( t ) = A B C D x ( t ) u ( t ) ,
where
A = A TT + A TS Φ ( I A SS Φ ) 1 A ST ,
B = B T + A TS Φ ( I A SS Φ ) 1 B S ,
C = C T + C S Φ I A SS Φ 1 A ST ,
D = D T + C S Φ ( I A SS Φ ) 1 B S .
Note: Well -posedness is very important in system design, and ill-posed systems are usually difficult to control or impossible to estimate [25,26,27]. Therefore, this paper assumes that System Γ is well-posed, which means that ( I A SS Φ ) 1 exists.
This paper intends to establish computationally effective conditions for the dissipativity analysis of large-scale networked systems Γ . The concept of dissipativity is very important in the system, whether from the perspective of theoretical research or the perspective of practical application. Roughly speaking, dissipative systems can be described as such properties. At any time, the energy that the system may provide cannot exceed the energy already supplied. We first describe the definition of dissipativity for System Γ .
The definition is related to the supply function. For the i-th subsystem Γ i , its supply function is defined as
s i ( u ( t , i ) , y ( t , i ) ) = y ( t , i ) u ( t , i ) T Q ( i ) y ( t , i ) u ( t , i ) ,
where Q ( i ) is a symmetric matrix of suitable dimensions.
Definition 1.
The large-scale networked system (1) and (2) with x ( 0 , i ) = 0 is said to be dissipative with supply function s i ( u ( t , i ) , y ( t , i ) ) if and only if there is a matrix P ( i ) 0 , such that,
t 0 t 1 i = 1 N s i ( u ( t , i ) , y ( t , i ) ) d t i = 1 N x T ( t 1 , i ) P ( i ) x ( t 1 , i ) i = 1 N x T ( t 0 , i ) P ( i ) x ( t 0 , i )
holds for all t 0 t 1 .
According to the definition, the supply function can be interpreted as the energy transferred to the system, which means that within a period of time [ t 0 , t 1 ] , as long as t 0 t 1 i = 1 N s i ( u ( t , i ) , y ( t , i ) ) d t is positive, the system will work normally, otherwise, the system will not work. i = 1 N x T ( t 1 , i ) P ( i ) x ( t 1 , i ) i = 1 N x T ( t 0 , i ) P ( i ) x ( t 0 , i ) represents the actual energy consumption of the system after the time interval t 1 t 0 . Therefore, Equation (5) shows that in any time period [ t 0 , t 1 ] , the energy change inside the system will not exceed the energy supplied by the outside.
It can be seen from the following derivation that the definition of dissipativity for the networked system (1) and (2) are consistent with the one in [28] based on (3).
The supply function based on the large-scale networked system (1) and (2) is as follows,
s ( u ( t ) , y ( t ) ) = i = 1 N s i ( u ( t , i ) , y ( t , i ) ) = y ( t , 1 ) u ( t , 1 ) y ( t , 2 ) u ( t , 2 ) y ( t , N ) u ( t , N ) T Q y ( t , 1 ) u ( t , 1 ) y ( t , 2 ) u ( t , 2 ) y ( t , N ) u ( t , N ) .
The supply function of System (3) is
s ( u ( t ) , y ( t ) ) = y ( t ) u ( t ) T Q 1 y ( t ) u ( t ) = y ( t , 1 ) y ( t , 2 ) y ( t , N ) u ( t , 1 ) u ( t , 2 ) u ( t , N ) T Q 1 y ( t , 1 ) y ( t , 2 ) y ( t , N ) u ( t , 1 ) u ( t , 2 ) u ( t , N ) = y ( t , 1 ) u ( t , 1 ) y ( t , 2 ) u ( t , 2 ) y ( t , N ) u ( t , N ) T Q y ( t , 1 ) u ( t , 1 ) y ( t , 2 ) u ( t , 2 ) y ( t , N ) u ( t , N ) ,
in which y ( t ) = c o l { y t , i | i = 1 N } , u ( t ) = c o l { u t , i | i = 1 N } , Q = d i a g { Q ( i ) | i = 1 N } and
Q = I 0 0 0 0 I 0 0 0 0 0 0 I 0 0 I 0 0 0 0 0 0 0 I 0 0 0 0 0 I Q 1 ( * ) T .
For the convenience of the following discussion, we introduce the following preliminary results that need to be used.
Lemma 1
([29]). For matrices L and U with compatible dimensions, there is a scalar α > 0 such that,
L U + U T L T α L L T + α 1 U T U .
Lemma 2
([30]). Given symmetric matrices F and G with appropriate dimensions, if v T F v > 0 can be obtained for every non-zero vector v satisfying v T G v = 0 , then there must be a real number r such that F + r G is positive definite, and vice versa.
Lemma 3
([29]). For an LMI in the form of an M × M ( M 1 ) block matrix: G P < 0 , except for the symmetric independent variable matrix P, other known coefficient matrices or constant matrices are all block diagonal matrices of appropriate dimensions, and all have N ( N > 1 ) diagonals. If it is divided into blocks, there is a full block feasible solution P for this LMI, and there must be a feasible solution for the diagonal division of the appropriate dimension.

3. Dissipativity Analysis

In [28], the dissipativity criterion of System (3) is proposed.
Lemma 4.
Assume that the networked system Γ is controllable. Then, System (3) is strictly dissipative with the supply function s ( u ( t ) , y ( t ) ) if and only if there exists a matrix P > 0 such that,
A T P + P A P B B T P 0 C D 0 I T Q C D 0 I < 0 .
Note that the matrices A, B, C and D in the condition of Lemma 4 all contain ( I A SS Φ ) 1 terms. Although the subsystem connection matrix Φ is sparse and the system parameter A * # , B * , C * and D * with * , # = T, S are block diagonal, the matrix ( I A SS Φ ) 1 is generally dense. When there are a large number of subsystems in large-scale networked systems, the calculation of matrices A, B, C, and D involves the inversion of high-dimensional matrices. Therefore, when the scale of the networked system increases, the computational complexity of Equation (9) will become very high.
Lemma 4 is a dissipative analysis condition based on lumped networked model. Due to the establishment of the lumped model, the connection relationship among subsystems is hidden inside the parameters, and its structural information is not effectively utilized. As a result, for large-scale networked systems, the use of this condition for dissipative testing will inevitably bring computational difficulties and even cannot be calculated.
Then, to reduce the computational difficulty caused by the increase of system scale, we establish a computationally efficient sufficient, and necessary condition for the strict dissipativity analysis of large-scale networked systems. This condition effectively utilizes the sparse structure of the subsystem connection matrix Φ in the networked system, that is, each subsystem is only connected to a limited number of other subsystems.
Theorem 1.
Assume that the networked system Γ is controllable. Then, System Γ is strictly dissipative with the supply function s i ( u ( t , i ) , y ( t , i ) ) if and only if there exists a symmetric positive definite matrix P and a positive scalar h such that,
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I h × ( * ) T I Φ Φ T Σ 2 I 0 0 A SS A ST B S < 0 .
Proof of Theorem 1.
Equation (9) can be expressed equivalently as follows,
( * ) T 0 P P 0 Q I 0 A B C D 0 I < 0 .
We express Equation (11) in the following equivalent form,
( * ) T 0 P P 0 Q I 0 0 0 0 0 I 0 0 0 0 I 0 I 0 0 I 0 0 I A B C D < 0 .
Matrices A, B, C, and D can be written as follows,
A B C D = A TT B T C T D T + A TS C S Φ ( I A SS Φ ) 1 A ST B S .
Substituting the above formula into Equation (12), we can get that,
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I × Φ ( I A SS Φ ) 1 A ST Φ ( I A SS Φ ) 1 B S I 0 0 I < 0 .
Then we define matrices F, M, and K as follows,
F = ( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I ,
M = Φ ( I A SS Φ ) 1 A ST B S I 0 0 I ,
K = I Φ I 0 0 A SS A ST B S .
Obviously, M T ( F ) M < 0 . When v = M ζ , ζ R # , for any v 0 , we can get that K v = 0 , which means v T F v > 0 . According to Lemma 2, there must be a real number h such that,
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I + h × ( * ) T I Φ Φ T Σ 2 I 0 0 A SS A ST B S > 0 .
That is,
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I h × ( * ) T I Φ Φ T Σ 2 I 0 0 A SS A ST B S < 0 .
The characterization of the left term of Equation (19) shows that if the inequality has a solution, then there must be h > 0 . So far, the necessity has been proved. Then, multiply the left and right sides of Equation (10) by the matrices M and M T respectively, and direct algebraic operations can complete the sufficiency proof.
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I Φ ( I A SS Φ ) 1 A ST B S I 0 0 I h × ( * ) T I Φ Φ T Σ 2 I 0 0 A SS A ST B S Φ ( I A SS Φ ) 1 A ST B S I 0 0 I < 0 .
The proof is completed. □
It can be seen that the condition in Lemma 4 hides the connection relationship among subsystems inside the parameters, while the left side of Equation (10) in Theorem 1 linearly depends on the symmetric matrix P, and the structure of the system is specifically reflected in it, which can effectively make use of the sparse structure of the subsystem connection matrix. Furthermore, the matrices A * # , B * , C * and D * with * , # = T , S are all block diagonal, and large-scale networked systems are sparse. Combined with the research on sparse semi-definite programming problems [31,32,33], when the system is relatively large, the computational complexity of solving the above sparse LMI is frequently lower than the condition in Lemma 4. This aspect can also be explained in subsequent numerical simulations. It is worth noting that the condition of Theorem 1 does not bring conservatism compared with Lemma 4, a dissipative criterion based on the lumped description.
When there are a huge number of subsystems, the strict dissipativity analysis using the condition in Theorem 1 may still encounter computational difficulties. To overcome this difficulty, we further explore the structural characteristics of the subsystem connection matrix Φ , and put forward the conditions for strictly dissipative analysis based on the parameters of each subsystem.
A simple derivation leads to the following relationship,
I Φ Φ T Φ T Φ 2 I 0 I 0 + 0 Φ T 0 Φ .
Combined with Lemma 4 and the properties of the subsystem connection matrix, on the basis of Equation (10), the necessary condition for the strict dissipativity analysis that only depends on the parameters of a single subsystem can be obtained.
Theorem 2.
Assume that the networked system Γ is controllable. A necessary condition for the strict dissipativity of System Γ with the supply function s i ( u ( t , i ) , y ( t , i ) ) is that each subsystem has a symmetric positive definite matrix P ( i ) and a positive scalar h such that,
( * ) T 0 P ( i ) P ( i ) 0 Q ( i ) 0 A TS ( i ) C S ( i ) 0 I A TT ( i ) C T ( i ) 0 0 B T ( i ) D T ( i ) I h × ( * ) T I Σ i 2 I 0 0 A SS ( i ) A ST ( i ) B S ( i ) < 0 .
For large-scale networked systems, sometimes the parameters of multiple subsystems are the same. In this case, using Theorem 2 is more efficient. A sufficient condition for strict dissipativity analysis which only depends on the parameters of a single subsystem is given below.
Theorem 3.
Assume that the networked system Γ is controllable. Then, System Γ is strictly dissipative with the supply function s i ( u ( t , i ) , y ( t , i ) ) if there exists a symmetric positive definite matrix P ( i ) and real number h 2 h 1 0 (or h 1 h 2 0 ) for each subsystem such that,
( * ) T 0 P ( i ) P ( i ) 0 Q ( i ) 0 A TS ( i ) C S ( i ) 0 I A TT ( i ) C T ( i ) 0 0 B T ( i ) D T ( i ) I ( * ) T h 1 I h 2 Σ i 2 I 0 0 A SS ( i ) A ST ( i ) B S ( i ) < 0 .
Proof of Theorem 3.
From Lemma 1, we can get
I Φ Φ T Φ T Φ ( 1 α ) I 0 ( * ) T + 1 1 α 0 Φ T ( * ) T .
Using the above formula and the conclusion in Theorem 1, one can obtain a sufficient condition for System Γ to be strictly dissipative is the existence of a symmetric positive definite matrix P and two real numbers h > 0 , α > 0 , such that,
( * ) T 0 P P 0 Q 0 A TS C S 0 I A TT C T 0 0 B T D T I h × ( * ) T ( 1 α ) I 0 ( * ) T + 1 1 α 0 Φ T ( * ) T I 0 0 A SS A ST B S < 0 .
Let h 1 = ( 1 α ) h , h 2 = 1 α 1 h , we can get h 2 = α 1 h 1 . Therefore, when α 1 , h 2 h 1 0 ; when α 1 , h 1 h 2 0 . The proof can be completed by combining Lemma 3. □
Compared with Theorem 1, the left side of Equation (23) in Theorem 3 is linearly related to the matrix P ( i ) , and its dimension is entirely governed by the dimension of the subsystem Γ i . When the state dimension of each subsystem is fixed, the computational complexity of Equation (23) only linearly depends on the number of subsystems N. Therefore, Theorem 3 has a substantially higher computing efficiency than Theorem 1 for large-scale networked systems. However, it should be noted that Theorems 2 and 3 are conservative.

4. Numerical Simulations

Several numerical simulations are employed in this section to demonstrate the efficacy of the strict dissipativity conditions presented in this paper. The simulation experiments are performed on a laptop computer with an Intel(R) Core(TM) i5-3230M CPU @ 2.60 GHz 2.60 GHz and 6 G RAM. In these simulations, we assume that m u i = m x i = m v i = m z i = m y i = 2 . Furthermore, all the parameters of the subsystem are independent of each other, and the parameters of each subsystem are randomly generated according to a continuous uniform distribution with an interval of [−0.9, 0.9]. The subsystem connection matrix is randomly generated, but there is only one non-zero element 1 in each row and column.
The conditions in Lemma 4, Theorem 1, and Theorem 3 are used to verify the strict dissipativity of the system. Among them, the conditions in Lemma 4 and Theorem 3 are calculated by the LMI toolbox provided by MATLAB, and the condition in Theorem 1 is calculated by the sparse solvers DSDP. For System Γ introduced in this paper, we generate 10 systems for calculation, and the average value and standard deviation of system dissipativity analysis calculation time are obtained. Table 1 and Table 2 give some results when the number of subsystems is among 2 and 45.
The tables show that the calculation time of the above three methods all increases with the increase of the number of subsystems. When the number of subsystems is 10 or less, the computational efficiency based on Lemma 4 is comparable to that of Theorem 1 and Theorem 3. This is because the dimensionality of the matrix inequality in Theorem 1 is higher than that in Lemma 4, and Theorem 3 requires several inequalities to be verified. With the expansion of the number of subsystems, when the number of subsystems is 20, 30, 40, 45, the ratio of calculation time based on the conditions in Lemma 4 and Theorem 1 is 1.0957, 1.3414, 1.3911, 1.4581. The ratio of the average computation time becomes larger and larger, which means that the computational efficiency of Theorem 1 is improved to some extent. Because the condition in Lemma 4 requires operations such as inversion of high-dimensional matrices. Clearly, Theorem 3 is more computationally efficient than both Lemma 4 and Theorem 1. This is due to the fact that the conditions of Theorem 3 are tested based on individual subsystem parameters, and their computational complexity only increases linearly with the number of subsystems N. In addition, due to the limitation of computer memory, the conditions in Lemma 4 and Theorem 1 may not be calculated, but Theorem 3, which is tested independently for each subsystem, can still be calculated. Therefore, Theorem 3 has more computational advantages in the dissipativity analysis of large-scale networked systems. It should be noted that Theorem 3 is conservative compared to Lemma 4 and Theorem 1.

5. Conclusions

This paper investigates the strict dissipativity of networked systems composed of a large number of subsystems. At first, according to the model of large-scale networked systems, the definition of the dissipativity of networked systems is given in this paper. Then, we study the dissipative criteria of networked systems. For large-scale networked systems, when the number of subsystems is large, the performance analysis using the existing linear system theory will encounter computational difficulties. Some LMI-form conditions for dissipativity analysis of large-scale networked systems are derived. Among them, Theorem 1 is a necessary and sufficient condition, which effectively utilizes the block diagonal structure of the system parameter matrix and the sparsity of the subsystem connection matrix. Combined with the use of sparse semidefinite programming tools, it is more efficient than the lumped analysis method for medium-scale networked systems. In addition, the proposed sufficient condition, and necessary condition only depend on the parameters of a single subsystem, which are more suitable for the dissipative analysis of networked systems with a large number of subsystems, but they are conservative compared with other conditions.
Regarding large-scale networked systems, the design of distributed controllers to ensure the dissipativity of large-scale networked systems will be investigated in further research. For instance, in UAV formation flight, relying on a centralized controller to observe the entire formation and control all UAVs at once is both impractical and increases operational costs in engineering applications. A more reasonable option would be to decentralize the controller to each UAV platform and achieve the overall objective by interacting and sharing information between platforms. Therefore, if the structural information of the network topology is capable of being fully utilized and a distributed control strategy that relies on local information sharing is adopted, the amount of data transmission in the network will be greatly reduced and the computational efficiency will be raised. Furthermore, the presence of quantization errors, time delays, data packet loss, and other phenomena when communicating networked among subsystems or among subsystems and their local controllers will be explored. In practical engineering applications, network connections would possibly be non-idealized, and the arrival of information delivered is frequently unable to be achieved immediately.

Author Contributions

Conceptualization, Y.S. and H.L.; methodology, Y.S.; software, Y.S.; validation, H.L.; formal analysis, Y.S.; investigation, Y.S.; resources, J.W. and H.L.; data curation, J.W.; writing—original draft preparation, Y.S.; writing—review and editing, Y.S.; visuConceptualizationalization, J.W. and H.L.; supervision, J.W. and H.L.; project administration, J.W. and H.L.; funding acquisition, J.W. and H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Natural Science Foundation of China under grant number 62273189; by the Shandong Provincial Natural Science Foundation under grant number ZR2019MF063 and ZR2020MF064.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, H.; Wang, D. Stability and stabilisation of a class of networked dynamic systems. Int. J. Syst. Sci. 2018, 49, 964–973. [Google Scholar] [CrossRef]
  2. Montanari, A.N.; Duan, C.; Aguirre, L.A.; Motter, A.E. Functional observability and target state estimation in large-scale networks. Proc. Natl. Acad. Sci. USA 2022, 119, e2113750119. [Google Scholar] [CrossRef] [PubMed]
  3. Yu, T.; Yu, L.; Xiong, J. Observer-based distributed control of large-scale systems under gossip communication protocol. Asian J. Control 2022, 24, 956–972. [Google Scholar] [CrossRef]
  4. Rahman, M.H.; Al-Naeem, M.; Ali, M.N.; Sufian, A. TFBN: A cost effective high performance hierarchical interconnection network. Appl. Sci. 2020, 10, 8252. [Google Scholar] [CrossRef]
  5. Bamieh, B.; Paganini, F.; Dahleh, M.A. Distributed control of spatially invariant systems. IEEE Trans. Autom. Control 2002, 47, 1091–1107. [Google Scholar] [CrossRef] [Green Version]
  6. Zaeri Amirani, M.; Bigdeli, N.; Haeri, M. Distributed fault detection and isolation in time-varying formation tracking UAV multi-agent systems. Asian J. Control. 2022, 25, 604–622. [Google Scholar] [CrossRef]
  7. Venkat, A.N.; Hiskens, I.A.; Rawlings, J.B.; Wright, S.J. Distributed output feedback MPC for power system control. In Proceedings of the 45th IEEE Conference on Decision and Control, San Diego, CA, USA, 13–15 December 2006; pp. 4038–4045. [Google Scholar]
  8. Knies, C.; Diermeyer, F. Data-driven test scenario generation for cooperative maneuver planning on highways. Appl. Sci. 2020, 10, 8154. [Google Scholar] [CrossRef]
  9. Lai, J.; Chen, S.; Lu, X.; Zhou, H. Formation tracking for nonlinear multi-agent systems with delays and noise disturbance. Asian J. Control 2015, 17, 879–891. [Google Scholar] [CrossRef]
  10. Willems, J.C. Dissipative dynamical systems part I: General theory. Arch. Ration. Mech. Anal. 1972, 45, 321–351. [Google Scholar] [CrossRef]
  11. Willems, J.C. Dissipative dynamical systems part II: Linear systems with quadratic supply rates. Arch. Ration. Mech. Anal. 1972, 45, 352–393. [Google Scholar] [CrossRef]
  12. Xie, S.; Xie, L.; De Souza, C.E. Robust dissipative control for linear systems with dissipative uncertainty. Int. J. Control 1998, 70, 169–191. [Google Scholar] [CrossRef]
  13. Tan, Z.; Soh, Y.C.; Xie, L. Dissipative control for linear discrete-time systems. Automatica 1999, 35, 1557–1564. [Google Scholar] [CrossRef]
  14. Zhu, D.; Zhang, W.; Liu, C.; Duan, J. Fractional-order hyperbolic tangent sliding mode control for chaotic oscillation in power system. Math. Probl. Eng. 2021, 2021. [Google Scholar] [CrossRef]
  15. Urata, K.; Inoue, M. Dissipativity reinforcement in feedback systems and its application to expanding power systems. Int. J. Robust Nonlinear Control 2018, 28, 1528–1546. [Google Scholar] [CrossRef]
  16. Zhang, Q.; Yan, H.; Chen, S.; Zhan, X.; Jiang, X. Finite-time dissipative control for networked control systems with hybrid-triggered scheme. Trans. Inst. Meas. Control 2021, 43, 891–901. [Google Scholar] [CrossRef]
  17. Agarwal, E.; Sivaranjani, S.; Gupta, V.; Antsaklis, P.J. Distributed synthesis of local controllers for networked systems with arbitrary interconnection topologies. IEEE Trans. Autom. Control 2020, 66, 683–698. [Google Scholar] [CrossRef]
  18. Chen, B.D.; Lall, S. Dissipation inequalities for distributed systems on graphs. In Proceedings of the 42nd IEEE International Conference on Decision and Control (IEEE Cat. No. 03CH37475), Maui, HI, USA, 9–12 December 2003; Volume 3, pp. 3084–3090. [Google Scholar]
  19. Ji, Z.; Lin, H.; Cao, S.; Qi, Q.; Ma, H. The complexity in complete graphic characterizations of multiagent controllability. IEEE Trans. Cybern. 2020, 51, 64–76. [Google Scholar] [CrossRef]
  20. Ji, Z.; Yu, H. A new perspective to graphical characterization of multiagent controllability. IEEE Trans. Cybern. 2016, 47, 1471–1483. [Google Scholar] [CrossRef]
  21. Yu-Geng, X. Large-scale systems control and complex networks-exploration and thinking. Acta Autom. Sin. 2013, 39, 1758–1768. [Google Scholar]
  22. Liu, H.; Yu, H. Stability and stabilization of spatially interconnected systems. Control Decis. 2020, 35, 749–756. [Google Scholar]
  23. Zhou, T.; Zhang, Y. On the stability and robust stability of networked dynamic systems. IEEE Trans. Autom. Control 2015, 61, 1595–1600. [Google Scholar] [CrossRef]
  24. Liu, H.; Yu, H. Finite-time control of continuous-time networked dynamical systems. IEEE Trans. Syst. Man, Cybern. Syst. 2018, 50, 4623–4632. [Google Scholar] [CrossRef]
  25. Zhou, T. Coordinated one-step optimal distributed state prediction for a networked dynamical system. IEEE Trans. Autom. Control 2013, 58, 2756–2771. [Google Scholar] [CrossRef]
  26. van Schuppen, J.H.; Boutin, O.; Kempker, P.L.; Komenda, J.; Masopust, T.; Pambakian, N.; Ran, A.C. Control of distributed systems: Tutorial and overview. Eur. J. Control 2011, 17, 579–602. [Google Scholar] [CrossRef] [Green Version]
  27. Khalil, I.; Doyle, J.; Glover, K. Robust and Optimal Control; Prentice Hall: Hoboken, NJ, USA, 1996. [Google Scholar]
  28. Duan, G.R.; Yu, H.H. LMIs in Control Systems: Analysis, Design and Applications; CRC Press: Boca Raton, FL, USA, 2013. [Google Scholar]
  29. Liu, H.; Yu, H. Decentralized state estimation for a large-scale spatially interconnected system. ISA Trans. 2018, 74, 67–76. [Google Scholar] [CrossRef]
  30. Pólik, I.; Terlaky, T. A survey of the S-lemma. SIAM Rev. 2007, 49, 371–418. [Google Scholar] [CrossRef] [Green Version]
  31. Andersen, M.S.; Pakazad, S.K.; Hansson, A.; Rantzer, A. Robust stability analysis of sparsely interconnected uncertain systems. IEEE Trans. Autom. Control 2014, 59, 2151–2156. [Google Scholar] [CrossRef]
  32. Benson, S.J.; Ye, Y. DSDP5 User Guide-Software for Semidefinite Programming; Technical Repor; Argonne National Lab.(ANL): Argonne, IL, USA, 2006. [Google Scholar]
  33. George, A.; Gilbert, J.R.; Liu, J.W. Graph Theory and Sparse Matrix Computation; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012; Volume 56. [Google Scholar]
Table 1. Average of calculation time.
Table 1. Average of calculation time.
Subsystem NumberLemma 4 (s)Theorem 1 (s)Theorem 3 (s)
20.2568090.1136300.280457
100.5065720.3573340.332980
204.7183641.8251620.550354
3049.84341310.1020701.132749
38160.80455125.1942601.990903
40201.48465631.9422832.235822
Table 2. Standard deviation of calculation time.
Table 2. Standard deviation of calculation time.
Subsystem NumberLemma 4 (s)Theorem 1 (s)Theorem 3 (s)
20.0051450.0198870.003355
100.0174580.0387210.005182
200.0355180.0512970.005435
300.5851510.1607430.011004
381.3052060.3171350.017758
401.8061960.3896990.092795
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, Y.; Wang, J.; Liu, H. Dissipativity Analysis of Large-Scale Networked Systems. Appl. Sci. 2023, 13, 1214. https://doi.org/10.3390/app13021214

AMA Style

Sun Y, Wang J, Liu H. Dissipativity Analysis of Large-Scale Networked Systems. Applied Sciences. 2023; 13(2):1214. https://doi.org/10.3390/app13021214

Chicago/Turabian Style

Sun, Yuanfei, Jirong Wang, and Huabo Liu. 2023. "Dissipativity Analysis of Large-Scale Networked Systems" Applied Sciences 13, no. 2: 1214. https://doi.org/10.3390/app13021214

APA Style

Sun, Y., Wang, J., & Liu, H. (2023). Dissipativity Analysis of Large-Scale Networked Systems. Applied Sciences, 13(2), 1214. https://doi.org/10.3390/app13021214

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop