Next Article in Journal
Ramping Up Customer-Centric Modular Design Projects: Mobile App Development for Pandemic Relief
Previous Article in Journal
Cognitive Network Science Reconstructs How Experts, News Outlets and Social Media Perceived the COVID-19 Pandemic
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Determining Asymptotic Stability and Robustness of Networked Systems

Department of Engineering, Reykjavik University, Menntavegur 1, 102 Reykjavik, Iceland
Submission received: 27 August 2020 / Revised: 5 October 2020 / Accepted: 26 October 2020 / Published: 29 October 2020
(This article belongs to the Section Complex Systems)

Abstract

:
This paper is motivated by the notion that coupling systems allows for mitigating the failure of individual ones. We present a novel approach to determining asymptotic stability and robustness of a network consisting of coupled dynamical systems, where individual system dynamics are represented through polynomial or rational functions. The analysis relies on a local analysis; thus, making it computationally implementable. We present an efficient computational method that relies on semidefinite programming. Importantly, for cases where multiple equilibrium points exist, we show how to determine regions around an asymptotically stable equilibrium point that bounds solutions. These regions increase when systems are coupled as we observe when applying the presented analysis framework to a mathematical model of a continuous stirred tank reactor. Importantly, the presented work has implications to other fields as well.

1. Introduction

We are currently experiencing rapid digitalisation and automatisation in industry and society, which are culminating in the notions of smart factories, smart cities, the industry of things, and the internet of things. This makes networked systems and their control increasingly gaining in importance [1,2,3]. Significantly, networked sensor and actuator systems are also already pervasive in many traditional industries, for instance, the petroleum industry [4]. Additionally, networks of coupled dynamical systems appear also in other areas such as modelling of the heart [5] and of flocking behaviour [2]. An important question in network science, which links the dynamics of individual nodes with the emergent properties of macroscopic networks [6], is how much failure a networked system can tolerate. From the perspective of control engineering, a networked system can be viewed as a system of coupled dynamical systems and the question can be formulated in terms of asymptotic stability, the ability to return to a steady state, and robustness, the ability to withstand failures and perturbations.
In recent years, tools have been developed to determine whether a given complex network, consisting of connected systems, is robustly controllable, that is, whether any desired state can be reached from any initial state for almost all values of network edge-weights. For linear systems, conditions based on structural controllability, which identifies the minimal number of inputs or controllable nodes for system controllability of a complex network, were developed in [7,8]. For nonlinear systems, an equivalent approach based on feedback vertex set control, which is an attractor-based control method and requires thorough system knowledge, has been presented in [9]. However, these works and the majority of related ones rather consider systems that do not represent dynamical processes, as highlighted in [10], where the authors studied networks with node behaviour, which better represents dynamical behaviour that is typical for many systems, and examined how this influences the properties of the network. Using their analysis framework, they generated more realistic, albeit still simplified, results.
In this work, we provide certificates for asymptotic stability and robustness of complex networked systems consisting of coupled realistic system models through local analysis, which makes the approach computationally implementable and does not require simplifying the individual node dynamics. For this, we use so-called barrier certificates [11] and semidefinite programming [12,13]. For dynamical systems, the authors of [11] introduced a safety verification method that relies on state functions termed barrier certificates. Such functions form contour levels and, thus, create different regions in the state space. For a given set of initial conditions, the zero-level set of a barrier certificate bounds corresponding system trajectories and, by doing so, separates safe from unsafe regions. This means that one obtains certificates or proofs of system safety without the need to compute system behaviour explicitly. The use of different types of barrier certificates is popular in fields as diverse as biological model-testing, using the noise inherent in such systems [14], and avoiding collisions in multi-robot systems [15]. For our results, barrier certificates provide means to assess the set of initial conditions that will converge towards a desired operating point.
Related work has been presented in [16], where global stability for coupled differential equations was investigated assuming a Lyapunov function [17] of the form V ( x ) = x x * ln x , where x * is a positive equilibrium point. This approach was later extended to coupled systems with time-varying coupling [18]. As a novelty, the approach presented in this paper provides not only means to search for a Lyapunov function that guarantees asymptotic stability of the coupled system, for systems with multiple equilibria, it also uses Lyapunov theory to provide regions of attraction, even when model parameters are uncertain (robustness).
In many areas of research, realistic dynamical system models are based on ordinary differential equations described by polynomial or rational functions. Assuming such models, in Section 2.1 and Section 2.2, we provide results on testing for asymptotic stability and robustness of a networked system through a local analysis of the system. Locality makes the analysis feasible and allows for its computational implementation, as we show in Section 2.3. Meanwhile, for simplicity, we first assume all-to-all coupling, in Section 2.4, we extend the results to arbitrary coupling topologies. Section 3 provides an application to process engineering, where we analyse the stability and robustness of a high biomass concentration steady state for connected continuous stirred tank reactors. Finally, we discuss our results and conclude the paper in Section 4.

2. Asymptotic Stability of Coupled Dynamical Systems

In this paper, we consider bidirectional coupling only. Furthermore, let us consider N dynamical systems that are linearly all-to-all coupled. Then, we represent the i-th system through
x ˙ i = f i ( x i ) + j = 1 N k D ( x j x i ) ,
where f i : R n R n describes the system’s dynamics for all i, i { 1 , 2 , , N } , x = [ x 1 x N ] T , x i ( t ) R n , and D R n × n . Vector x is the state vector. Matrix D is the matrix that denotes the variables that are used in the coupling; that is, D ( k , l ) 0 if and only if states x i ( k ) and x j ( l ) are coupled, for all i , j and i j . Positive constant k corresponds to the coupling strength. Next, we provide results on asymptotic stability of the coupled system based on Lyapunov theory.

2.1. Unique Equilibrium Point

First, let us assume that the networked system described by N coupled systems given by Equation (1) has a unique global equilibrium point. Without loss of generality, let it be the origin; otherwise, the equilibrium point can be shifted to the origin through a change of variables. Moreover, we assume that x ˙ 1 = f 1 ( x 1 ) is asymptotically stable. We do this, because, without proof, the networked system can only be asymptotically stable if at least one individual system is asymptotically stable, whose state we denote by x 1 , without loss of generality. If we consider all other systems to be similar to this system but to possess a certain amount of bounded parameter uncertainty then condition (4), which we derive in the following, guarantees that the coupled system is asymptotically stable. This allows us to assess the robustness of a networked system.
Now, consider the following Lyapunov function,
V ( x ) = i = 2 N 1 2 1 N 1 x 1 T M x 1 + V ¯ ( X i ) = i = 2 N V ˜ ( x 1 , X i ) ,
where X i = x i x 1 , M = M T 0 , and V ˜ ( x 1 , X i ) 0 . Moreover, V ˜ ( x 1 , X i ) > 0 if neither x 1 = 0 nor X i = 0 , for all i. Thus, V ( x ) > 0 unless x = 0 . Then,
V ˙ ( x ) = i = 2 N 1 2 1 N 1 x ˙ 1 T M x 1 + x 1 T M x ˙ 1 + V ¯ ˙ ( X i ) = i = 2 N 1 N 1 x 1 T M x ˙ 1 + V ¯ ˙ ( X i ) = i = 2 N 1 N 1 x 1 T M f 1 ( x 1 ) + k D j = 2 N X j + V ¯ ( X i ) X i ( f i ( x i ) f 1 ( x 1 ) k N D X i ) = i = 2 N 1 N 1 x 1 T M f 1 ( x 1 ) + x 1 T k M D X i + V ¯ ( X i ) X i ( f i ( x i ) f 1 ( x 1 ) k N D X i ) .
The second equality in the first line of Equation (3) follows from the symmetry of M. Thus, if
1 N 1 x 1 T M f 1 ( x 1 ) + x 1 T k M D X i + V ¯ ( X i ) X i ( f i ( x i ) f 1 ( x 1 ) k N D X i ) < 0
holds for all i unless x 1 = 0 and X i = 0 for all i, and, thus, unless x = 0 , then V ˙ < 0 unless x = 0 , which guarantees asymptotic stability of the coupled system. Finally, for all i, we consider the differences in function f i ( x i ) to depend on r real parameters that are uncertain but whose values are bounded. Let these parameters be elements of vector q Q R r . Then, we can rewrite condition (4) by letting g ( q 1 , x 1 ) f 1 ( x 1 ) and g ( q i , x i ) f i ( x i ) . In this case, if condition (4) holds, for all x 1 , x i R n and all q 1 , q i Q , then asymptotic stability of the coupled system is guaranteed.

2.2. Multiple Equilibrium Points

Second, if the networked system possesses multiple equilibrium points, then inequality (4) can only hold in a certain region of the state space. In order to estimate this region, we assume, again, that the equilibrium point, whose stability properties we wish to investigate, is the origin. Then, if inequality (4) holds for V ˜ ( x 1 , X i ) < c , c R , c > 0 , and all i then V ˜ ( x 1 , X i ) decreases with time for all i and, thus, so does V ( x ) , which means that the origin is asymptotically stable in the region bounded by V ( x ) = c . As we will show in Section 3, in addition to above robustness analysis, determining this region allows us to assess how far a system can deviate from a particular steady state and still remain attracted by the latter.

2.3. Computational Implementation

2.3.1. Coupled Linear Systems

Consider N coupled linear dynamical systems of the form x ˙ 1 = A x and x ˙ i = B x for all i > 1 , where A , B R n × n , matrix A is stable, and matrix B is a matrix that is similar to A, however, with some uncertain but bounded entries. Then, for V ¯ ( X i ) = X i T P X i , where P = P T 0 , inequality (4) becomes
x 1 T M 1 N 1 A k D + P ( A k N D ) x 1 + x 1 T k M D + P ( k N D B ) x i + x i T P ( k N D A ) x 1 + x i T P ( B k N D ) x i < 0 .
Given k, the problem of determining matrices M and P such that Equation (2) and inequality (5) hold can be cast as a semidefinite programme of the following form
given A , B , D , N , k search for M , P subject to M = M T 0 , P = P T 0 , ( 5 ) .
Problem (6) is a classic case of linear matrix inequalities [19] that can be solved efficiently using, for example, YALMIP, a free, third-party MATLAB toolbox [20], which can handle bounded uncertainties in B [21] to obtain robustness certificates.

2.3.2. Coupled Nonlinear Systems

For nonlinear systems, we investigate those ones whose dynamics are modelled using polynomial or rational functions. Consider the real-valued polynomial function F ( x ) of degree 2 d , x R n . Testing for non-negativity is NP-hard [22]. However, a sufficient condition for F ( x ) to be nonnegative is that it can be decomposed into a sum of squares (SOS) [23]:
F ( x ) = i f i 2 ( x ) 0 ,
where f i are polynomial functions. Now, F ( x ) is a SOS if and only if there exists a positive semidefinite matrix R and
F ( x ) = i f i 2 ( x ) = χ T R χ , χ = [ 1 , x ( 1 ) , x ( 2 ) , , x ( n ) , x ( 1 ) x ( 2 ) , , x ( n ) d ] .
The entries of vector χ consist of all monomial combinations of the elements of vector x up to degree d (including x ( i ) 0 = 1 ) and, thus, its length is = n + d d . Note that R is not necessarily unique. However, Equation (8) poses certain constraints on R of the form trace ( A j R ) = c j , where A j and c j are appropriate matrices and constants respectively. As an illustration, in the following, we reproduce Example 3.5 of [23].
Consider the following polynomial function, where z 1 = x 1 2 , z 2 = x 2 2 , z 3 = x 1 x 2 :
F ( x ) = 2 x 1 4 + 2 x 1 3 x 2 x 1 2 x 2 2 + 5 x 2 4 = x 1 2 x 2 2 x 1 x 2 T q 11 q 12 q 13 q 12 q 22 q 23 q 13 q 23 q 33 x 1 2 x 2 2 x 1 x 2 = q 11 x 1 4 + q 22 x 2 4 + ( q 33 + 2 q 12 ) x 1 2 x 2 2 + 2 q 13 x 1 3 x 2 + 2 q 23 x 1 x 2 3 .
This leads to the following set of linear equalities:
q 11 = 2 , q 22 = 5 , q 33 + 2 q 12 = 1 , 2 q 13 = 2 , 2 q 23 = 0 .
Thus, j = 1 , , 5 , and, for example, c 1 = 2 , c 3 = 1 ,
A 1 = 1 0 0 0 0 0 0 0 0 , and A 3 = 0 1 0 1 0 0 0 0 1 .
In general, in order to find R, we solve the optimisation problem associated with the following semidefinite program:
minimise trace ( A 0 R ) subject to trace ( A j R ) = c j , j = 1 , , m R = R T 0 .
In this paper, to solve SOS programmes, we use SOSTOOLS [24], a free, third-party MATLAB toolbox. For instance, the problem of determining Lyapunov function V ˜ ( x 1 , X i ) , given by Equation (2), such that (4) holds can be cast as follows
given f 1 ( x 1 ) , f i ( x i ) , D , N , k search for V ˜ ( x 1 , X i ) subject to ε ( x 1 T x 1 + x i T x i ) + V ˜ ( x 1 , X i ) is SOS , ε ( x 1 T x 1 + x i T x i ) V ˜ ˙ ( x 1 , X i ) is SOS , 0 < ε 1 .
Here, some additional remarks:
  • Consider a rational function F ( x ) , F ( x ) = f ( x ) g ( x ) , where f ( x ) and g ( x ) are polynomial functions. Then, F ( x ) 0 if (9) is feasible with χ T R χ = F ( x ) g 2 ( x ) or with χ T R χ = F ( x ) g ( x ) if g ( x ) > 0 .
  • If
    F ( x ) + p ( x ) h ( x ) = i g i 2 ( x ) 0 , p ( x ) 0 , h ( x ) = 0 if a i x ( i ) b i i > 0 otherwise ,
    then F ( x ) 0 if a i x ( i ) b i for all i, where a i , b i are constants. This can be used to show that F ( x ) is nonnegative in a specific region of the state and/or parameter space.
  • Finally, the computational time necessary to solve SOS decomposition problems scales badly with the size of the problem, as the length of vector is = n + d d .

2.4. Different Coupling Configurations

The presented results can be extended to any symmetric linear coupling topology with corresponding Laplacian matrix C by letting V ¯ ( X i ) = X i T P X i , for all i, in Equation (2), where P = P T 0 . To show this, first, we define the Laplacian matrix of the completely connected graph given by U R N × N through
U = N I N E 0 , E = e e T , e = [ 1 1 ] T .
Note that the eigenvalues of U are λ 1 = 0 and λ i = N for all i > 1 . Then, the following lemma establishes a relation between C and U.
Lemma 1.
Let U , C R N × N , where U is the Laplacian matrix of the completely connected graph defined in Equation (11) and C is a symmetric Laplacian matrix of rank ( N 1 ) . Then, for all y R N ,
y T C y λ min ( C ) N y T U y ,
where λ min ( C ) denotes the smallest positive eigenvalue of the matrix C .
Proof. 
First, note that C ( i , j i ) 0 . Now, U E = 0 and C E = 0 . Hence, U U = N U and C U = N C . By the Spectral Theorem, there exists a unitary matrix S such that S 1 = S T and S C S T = Λ C , where Λ C is the diagonal matrix containing the eigenvalues of C and, without loss of generality, its first entry is 0. Moreover,
S C U S T = S C S T S U S T = S C S T N = S C S T Λ U
implies that S U S T = Λ U , where Λ U is the diagonal matrix containing the eigenvalues of U, which are N but for the first one, which is 0. Now, for all y R N and Λ C being the diagonal matrix containing the eigenvalues of C,
y T C y = y T S T Λ C S y = v T Λ C v λ min ( C ) N v T Λ U v = λ min ( C ) N y T S T Λ U S y = λ min ( C ) N y T U y ,
which completes the proof. □
Consider, again, difference X i x i x 1 , where X i R n and X 1 = 0 [25]. Let matrix F R N × N be defined as follows,
F = 0 0 0 0 1 1 0 0 1 0 1 0 1 0 0 1 .
Then,
X = ( F I n ) x , X , x R N n ,
and
X ˙ = ( F I n ) x ˙ .
Moreover, note that for symmetric coupling determined by Laplacian matrix C , C = C F and F U = N F . Let f ˜ ( x ) = [ f 1 ( x 1 ) f N ( x N ) ] T . Then,
j = 2 N V ¯ ˙ ( X i ) = X T ( I N P ) X ˙ = X T ( F P ) x ˙ = X T ( F P ) ( f ˜ ( x ) + k ( C D ) x ) = X T ( F P ) ( f ˜ ( x ) + k ( C F D ) x ) = X T ( F P ) ( f ˜ ( x ) + k ( C D ) X ) X T ( F P ) ( f ˜ ( x ) + α k ( U D ) X ) = i = 2 N X i T P ( f i ( x i ) f 1 ( x 1 ) ) + α k X T ( F U P D ) X = i = 2 N X i T P ( f i ( x i ) f 1 ( x 1 ) ) α N k X T ( F P D ) X = i = 2 N X i T P ( f i ( x i ) f 1 ( x 1 ) α k N D X i ) ,
where the inequality follows from Lemma 1, α N = λ min ( ( F C ) gL ) , and ( F C ) gL is the matrix obtained by deleting the first row and the first column of F C . It is the so-called grounded Laplacian [26,27], where the deleted row and deleted column correspond to the grounded variable—in our case it is x 1 . A grounded Laplacian matrix is a principal sub-matrix of the Laplacian matrix. The fact that, for square matrices A , B , eig ( A B ) = eig ( B A ) , has been proven in [28] and, thus, eig ( F C ) = eig ( C F ) = eig ( C ) . Since the first row of matrix F C is zero, it follows that λ min ( ( F C ) gL ) = λ min ( C ) . For illustration, we provide examples of λ min ( C ) of different unweighted coupling configurations.
(a)
All-to-all coupling ( C = U ): λ min ( C ) = N .
(b)
Star-configuration: λ min ( C ) = 1 .
(c)
Ring of diffusively coupled systems: λ min ( C ) = 4 sin 2 π N .
(d)
Ring of 2 k -nearest neighbour coupled systems [29]: λ min ( C ) 2 π 2 k ( k + 1 ) ( 2 k + 1 ) / 3 N 2 if k N .
The different coupling configurations are illustrated in Figure 1.

3. Coupled Continuous Stirred Tank Reactors

Consider a continuous stirred tank reactor for biomass production held at a constant temperature [30]. Concentrations of feed or substrate are denoted by s and those of the product by x (Figure 2). The reaction describing the inflow and outflow of substrate and outflow/harvesting of product is given by x s . The reaction describing the growth of product by consuming substrate is given by x s k 2 x , where reaction constant k is described below.
The corresponding mathematical model based on ordinary differential equations is given by
x ˙ = 2 x s ( s + 1 ) 2 0.4 x , s ˙ = u 2 x s ( s + 1 ) 2 0.4 s ,
where we denote the input rate of substrate by u. Moreover, we assume an outflow rate of 0.4 (unitless) and a reaction rate of k x s , where k = 2 ( s + 1 ) 2 implies inhibition for high concentrations of substrate (poisoning). We denote equilibrium points by x * and s * . Figure 3 shows the phase plot with the vector field for u = 2 , which depicts the behaviour of the system, determined through the direction set by the vector arrows, in dependence of its state. For u = 2 , there are three equilibrium points, which are clearly visible. The equilibrium point that corresponds to higher concentrations of product is given by x * = 4.618 and s * = 0.382 . Note that the system described by Equation (16) is a nonnegative system such that for nonnegative initial conditions, solutions remain in the nonnegative orthant [31]. To obtain a set that bounds solutions, we let x ¯ = x + x * and s ¯ = s + s * , first, solve the following problem
given Equation ( 16 ) , u = 2 , β R search for M subject to p 1 ( x ¯ , s ¯ ) , p 2 ( x ¯ , s ¯ ) , p 3 ( x ¯ , s ¯ ) are SOS functions of degree 4 , 0.001 ( x ¯ 2 + s ¯ 2 ) + V 1 is SOS , V 1 = 1 2 x ¯ s ¯ T M x ¯ s ¯ , 0.001 ( x ¯ 2 + s ¯ 2 ) V ˙ 1 · ( s + 1 ) 2 + p 1 ( x ¯ , s ¯ ) ( x ¯ 2 + s ¯ 2 β ) p 2 ( x ¯ , s ¯ ) x p 3 ( x ¯ , s ¯ ) s is SOS ,
for different values of β , and, then, solve
given V 1 minimise γ , γ R subject to p 1 ( x ¯ , s ¯ ) , p 2 ( x ¯ , s ¯ ) , p 3 ( x ¯ , s ¯ ) are SOS functions of degree 2 , V 1 γ p 1 ( x ¯ , s ¯ ) ( x ¯ 2 + s ¯ 2 β ) p 2 ( x ¯ , s ¯ ) x p 3 ( x ¯ , s ¯ ) s is SOS .
The boundary of the set is given by V 1 = γ . Figure 3 depicts the sets defined by β = 7 and β = 9.95 . Their boundaries are given by the dashed red lines. For β > 9.95 , programme (17) is infeasible and remains so even after increasing the degree of polynomials p 1 , p 2 , and p 3 . Note that outside the regions, whose boundaries are given by the dashed red line in Figure 3, one might risk washout or depletion of product by being attracted to the equilibrium point given by x * = 0 and s * = 5 .
Now, by interconnecting reactors as in Equation (1), we can mitigate the danger of washout that might result, for example, from disturbances. Note that in our modelling, we assume full control over all flow rates, that is, in all inlets, where we let u = 2 , and between all reactors, where k = 5 . For N > 1 , we solve the following problem, which is analogous to problem (17),
given Equation ( 16 ) , Equation ( 1 ) , u = 2 , k = 5 , D = I 2 , N > 1 , β R search for V ˜ ( x ¯ 1 , s ¯ 1 , x ¯ i , s ¯ i ) subject to p j ( x ¯ 1 , s ¯ 1 , x ¯ i , s ¯ i ) are SOS functions of degree 4 , j { 1 , , 5 } , V ˜ = 1 2 1 N 1 x 1 T M x 1 + V ¯ ( X i ) is a SOS function of degree 2 , 0.001 ( x ¯ 1 2 + s ¯ 1 2 + x ¯ i 2 + s ¯ i 2 ) + V ˜ is SOS , 0.001 ( x ¯ 1 2 + s ¯ 1 2 + x ¯ i 2 + s ¯ i 2 ) V ˜ ˙ · ( s ¯ 1 + 1 ) 2 ( s ¯ i + 1 ) 2 + p 1 · ( x ¯ 1 2 + s ¯ 1 2 + x ¯ i 2 + s ¯ i 2 β ) p 2 · x 1 p 3 · s 1 p 4 · x i p 5 · s i is SOS .
By, first, solving a problem analogous to problem (18), with V 1 = 1 2 1 N 1 x ¯ 1 s ¯ 1 T M x ¯ 1 s ¯ 1 , for β = 12 and β = 17.5 , for different N > 1 , and, then, assuming that all systems are at equilibrium but the first one, we find that solutions remain bounded and return to equilibrium if x 1 and s 1 are initially within the regions bounded by the black solid lines in Figure 3. Finally, for investigating the robustness of our previous results, we repeat above computation with additional free variables that account for uncertainty in model parameters and corresponding constraints. For instance, for an uncertain input, such as 2 u 2.153 , which translates to 4.618 x * 5 , x * becomes an additional free variable. We add the additional constraints to the last line of programme (19). They have the following form
+ p 6 · ( x 1 * 4.168 ) ( x 1 * 5 ) + p 7 · ( x i * 4.168 ) ( x i * 5 ) ,
where p 6 ( x ¯ 1 , s ¯ 1 , x ¯ i , s ¯ i ) and p 7 ( x ¯ 1 , s ¯ 1 , x ¯ i , s ¯ i ) are SOS functions of degree 4. We, again, observe that the regions that bound solutions are given by the black solid lines.

4. Discussion and Conclusions

In this paper, we presented a novel approach for determining asymptotic stability and robustness of a network consisting of coupled dynamical systems, where individual system dynamics are represented through polynomial or rational functions. The analysis framework relies on local analysis; thus, making it computationally implementable. With respect to the latter, we presented an efficient method that relies on semidefinite programming. Importantly, for cases, where multiple equilibrium points exist, we show how to determine regions around an asymptotically stable equilibrium point that bounds solutions. The novelty of the presented research lies in the analysis of complex networks that consist of realistic nonlinear systems, which are of relevance to science and industry, by employing barrier certificates to provide guarantees for asymptotic stability and robustness, as opposed to previous research, which does rather not consider realistic dynamical systems as network nodes.
For example, for biomass production in continuous stirred tank reactors, on the one hand, one wishes to increase yield by providing large concentrations of substrate; on the other hand, too much of it can lead to a decrease or even washout of product. Intuitively, coupling a few reactors should allow for greater deviations of individual reactors from the optimal operating point without risking washout. Using a model of such a networked system, we showed how to determine just how much deviation is (surely) safe. Importantly, assuming all reactors at the optimal operating point but one, we showed that deviations deemed unsafe for a single reactor are now still within the safe zone. Finally, wind farms, the industrial internet of things, and the spread of epidemics, have all in common complex dynamical behaviour because of the interconnectedness of an increasingly large number of individual systems or individuals. However, the fact that they consist of complex networks makes achieving the goal of robustly controlling such systems particularly challenging. The research presented in this paper provides an approach to overcoming this challenge.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Notation

R , R m × n real numbers, m × n real matrices
A ( i , j ) ( i , j ) th entry of matrix A R m × n
I n the identity matrix, I n R n × n
A T transpose of matrix A R m × n
x ˙ derivative of x with respect to time variable t
A 0 , B 0 matrix A is positive definite, matrix B is positive semidefinite
A 0 , B 0 matrix A is negative definite, matrix B is negative semidefinite
A B , A R m × n , B R p × q The Kronecker product: A ( 1 , 1 ) B A ( 1 , n ) B A ( 2 , 1 ) B A ( 2 , n ) B A ( m , 1 ) B A ( m , n ) B

References

  1. Bulatov, Y.N.; Kryukov, A.V. A multi-agent control system of distributed generation plants. In Proceedings of the International Conference on Industrial Engineering, Applications and Manufacturing (ICIEAM), St. Petersburg, Russia, 16–19 May 2017. [Google Scholar]
  2. Olfati-Saber, R.; Fax, J.A.; Murray, R.M. Consensus and Cooperation in Networked Multi-Agent Systems. Proc. IEEE 2007, 95, 215–233. [Google Scholar] [CrossRef] [Green Version]
  3. Ge, X.; Yang, F.; Han, Q.L. Distributed networked control systems: A brief overview. Inf. Sci. 2017, 380, 117–131. [Google Scholar] [CrossRef]
  4. Pajic, M.; Sundaram, S.; Pappas, G.J.; Mangharam, R. The wireless control network: A new approach for control over networks. IEEE Trans. Autom. Control 2011, 56, 2305–2318. [Google Scholar] [CrossRef] [Green Version]
  5. dos Santos, A.M.; Lopes, S.R.; Viana, R.L. Rhythm synchronization and chaotic modulation of coupled Van der Pol oscillators in a model for the heartbeat. Physica A 2004, 338, 335–355. [Google Scholar] [CrossRef]
  6. Vespignani, A. Twenty years of network science. Nature 2018, 558, 528–529. [Google Scholar] [CrossRef]
  7. Liu, Y.Y.; Slotine, J.J.; Barabási, A.L. Controllability of complex networks. Nature 2011, 473, 167–173. [Google Scholar] [CrossRef]
  8. Liu, Y.Y.; Barabási, A.L. Control principles of complex systems. Rev. Mod. Phys. 2016, 88, 035006. [Google Scholar] [CrossRef] [Green Version]
  9. Zanudo, J.G.T.; Yang, G.; Albert, R. Structure-based control of complex networks with nonlinear dynamics. Proc. Natl. Acad. Sci. USA 2017, 114, 7234–7239. [Google Scholar] [CrossRef] [Green Version]
  10. Leitold, D.; Vathy-Fogarassy, Á.; Abonyi, J. Controllability and observability in complex networks—The effect of connection types. Sci. Rep. 2017, 7, 151. [Google Scholar] [CrossRef]
  11. Prajna, S.; Jadbabaie, A.; Pappas, G. Stochastic Safety Verification Using Barrier Certificates. In Proceedings of the IEEE Conference on Decision and Control (CDC) (IEEE Cat. No.04CH37601), Nassau, Bahamas, 14–17 December 2005; pp. 929–934. [Google Scholar]
  12. Vandenberghe, L.; Boyd, S. Semidefinite programming. SIAM Rev. 1996, 38, 49–95. [Google Scholar] [CrossRef] [Green Version]
  13. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  14. August, E. Using Noise for Model-Testing. J. Comput. Biol. A J. Comput. Mol. Cell Biol. 2012, 19, 968–977. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, L.; Ames, A.; Egerstedt, M. Safety Barrier Certificates for Collisions-Free Multirobot Systems. IEEE Trans. Robot. 2017, 33, 661–674. [Google Scholar] [CrossRef]
  16. Li, M.; Shuai, Z. Global-stability problem for coupled systems of differential equations on networks. J. Differ. Equ. 2010, 248, 1–20. [Google Scholar] [CrossRef] [Green Version]
  17. Khalil, H.K. Nonlinear Systems, 3rd ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 2000. [Google Scholar]
  18. Feng, J.; Luan, X. Stability analysis for coupled systems with time-varying coupling structure. Adv. Differ. Equ. 2017, 2017, 191. [Google Scholar] [CrossRef] [Green Version]
  19. Boyd, S.; El Ghaoui, L.; Feron, E.; Balakrishnan, V. Studies in Applied Mathematics. In Linear Matrix Inequalities in System and Control Theory; SIAM: Philadelphia, PA, USA, 1994; Volume 15. [Google Scholar]
  20. Löfberg, J. YALMIP: A Toolbox for Modeling and Optimization in MATLAB. In Proceedings of the CACSD Conference, New Orleans, LA, USA, 2–4 September 2004. [Google Scholar]
  21. Löfberg, J. Automatic robust convex programming. Optim. Methods Softw. 2012, 27, 115–129. [Google Scholar] [CrossRef] [Green Version]
  22. Murty, K.G.; Kabadi, S.N. Some NP-complete problems in quadratic and nonlinear programming. Math. Program. 1987, 39, 117–129. [Google Scholar] [CrossRef]
  23. Parrilo, P.A. Semidefinite programming relaxations for semialgebraic problems. Math. Program. Ser. B 2003, 96, 293–320. [Google Scholar] [CrossRef]
  24. Papachristodoulou, A.; Anderson, J.; Valmorbida, G.; Prajna, S.; Seiler, P.; Parrilo, P.A. SOSTOOLS—Sum of Squares Optimization Toolbox, User’s Guide. 2018. Available online: http://sysos.eng.ox.ac.uk/sostools/ (accessed on 29 October 2020).
  25. Belykh, V.N.; Belykh, I.V.; Hasler, M. Connection graph stability method for synchronized coupled chaotic systems. Physica D 2004, 195, 159–187. [Google Scholar] [CrossRef]
  26. Miekkala, U. Graph properties for splitting with grounded laplacian matrices. BIT 1993, 33, 485–495. [Google Scholar] [CrossRef]
  27. Jadbabaie, A.; Motee, N.; Barahona, M. On the Stability of the Kuramoto Model of Coupled Nonlinear Oscillators. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004; pp. 4296–4301. [Google Scholar]
  28. Thurston, H.S. On the Characteristic Equations of Products of Square Matrices. Am. Math. Mon. 1931, 38, 322–324. [Google Scholar] [CrossRef]
  29. Barahona, M.; Pecora, L.M. Synchronization in Small-World Systems. Phys. Rev. Lett. 2002, 89, 054101. [Google Scholar] [CrossRef] [Green Version]
  30. Chmiel, H. Bioprozesstechnik, 3rd ed.; Springer: Heidelberg, Germany, 2001. [Google Scholar]
  31. Haddad, W.M.; Chellaboina, V.; August, E. Stability and Dissipativity Theory for Nonnegative Dynamical Systems: A Thermodynamic Framework for Biological and Physiological Systems. In Proceedings of the 40th IEEE Conference on Decision and Control (Cat. No.01CH37228), Orlando, FL, USA, 4–7 December 2001; pp. 442–458. [Google Scholar]
Figure 1. Four different coupling configurations. Here, N = 10 in (a,c,d), N = 11 in (b), and k = 2 in (d).
Figure 1. Four different coupling configurations. Here, N = 10 in (a,c,d), N = 11 in (b), and k = 2 in (d).
Systems 08 00039 g001
Figure 2. Schematic view of a continuous stirred tank reactor.
Figure 2. Schematic view of a continuous stirred tank reactor.
Systems 08 00039 g002
Figure 3. Phase plot with vector field and corresponding contour lines, for the system described by Equation (16), and regions that bound solutions (red dashed lines). For coupling of a few such systems as in Equation (1) and for N > 1 , solutions are bounded by the black solid lines.
Figure 3. Phase plot with vector field and corresponding contour lines, for the system described by Equation (16), and regions that bound solutions (red dashed lines). For coupling of a few such systems as in Equation (1) and for N > 1 , solutions are bounded by the black solid lines.
Systems 08 00039 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

August, E. Determining Asymptotic Stability and Robustness of Networked Systems. Systems 2020, 8, 39. https://0-doi-org.brum.beds.ac.uk/10.3390/systems8040039

AMA Style

August E. Determining Asymptotic Stability and Robustness of Networked Systems. Systems. 2020; 8(4):39. https://0-doi-org.brum.beds.ac.uk/10.3390/systems8040039

Chicago/Turabian Style

August, Elias. 2020. "Determining Asymptotic Stability and Robustness of Networked Systems" Systems 8, no. 4: 39. https://0-doi-org.brum.beds.ac.uk/10.3390/systems8040039

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