Next Article in Journal
Intelligent Reflecting Surface-Assisted Wireless Communication Using RNNs: Comprehensive Insights
Previous Article in Journal
A Zassenhaus Lemma for Digroups
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation and Control of Positive Complex Networks Using Linear Programming

School of Information and Communication Engineering, Hainan University, Haikou 570228, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(19), 2971; https://doi.org/10.3390/math12192971
Submission received: 27 July 2024 / Revised: 11 September 2024 / Accepted: 15 September 2024 / Published: 25 September 2024
(This article belongs to the Section Network Science)

Abstract

:
This paper focuses on event-triggered state estimation and control of positive complex networks. An event-triggered condition is provided for discrete-time complex networks by which an event-based state estimator and an estimator-based controller are designed through matrix decomposition technology. Thus, the system is converted to an interval uncertain system. The positivity and the L 1 -gain stability of complex networks are ensured by resorting to a co-positive Lyapunov function. All conditions are solvable in terms of linear programming. Finally, the effectiveness of the proposed state estimator and controller are verified by a numerical example. The main contributions of this paper are as follows: (i) A positive complex network framework is constructed based on an event-triggered strategy, (ii) a new state estimator and an estimator-based controller are proposed, and (iii) a simple analysis and design approach consisting of a co-positive Lyapunov function and linear programming is presented for positive complex networks.

1. Introduction

A complex network is an interconnected structure consisting of numerous nodes and connections, characterized by intricate relationships and interactions. With the continuous advancement of modern technology, the research on complex networks has received considerable attention due to their pervasive applications in biology [1], neural networks [2], physiological processes [3], and social sciences [4].
As expected, the work [5] was concerned with the synchronization problem of coupled complex-valued neural networks, thereby ensuring the stability of the state estimation system. The work [6] achieved fixed synchronization of complex networks by extensively borrowing network interconnections and interactions between nodes. On the basis of the pulse control theory, the work [7] provided sufficient conditions for the exponential stability of impulse perturbations. With the deepening of research, it has been gradually recognized that in some specific contexts, some networks present a special property-positivity, where all variables in the network remain in a non-negative state. This positive condition not only simplifies the network analysis, but also provides a new perspective for the study of positive systems. Positive systems are a class of complex dynamic systems whose state values and outputs are non-negative when the initial conditions of the system are non-negative [8,9]. In recent years, positive systems have gained significant traction in diverse fields such as multi-agent systems, biology, network communication, environmental science, and many others. For more examples of positive systems and complex networks, the reader can refer to [10,11]. Currently, few results are reported for positive complex networks, for the following reasons: (i) Determining how to deal with the positivity of complex network systems is a difficult task, (ii) determining how to design state observers of positive complex networks is a considerable challenge, and (iii) determining how to design the synchronization control gain needs to introduce a novel approach.
Unsurprisingly, notable advancements have been achieved in the analysis of dynamic processes in complex networks [12]. Alongside the extensive research on synchronization control in complex networks, the equally important problem of state estimation should be emphasized due to its potential applications in understanding the internal dynamics of complex networks [13,14]. The work [15] introduced a new observer based on unknown disturbance models to reduce network complexity. In work [16], a set of event-based state estimators was constructed to reduce unnecessary data transmission in communication. While numerous studies have focused on state estimators and controllers in complex networks, a significant portion of them rely on quadratic forms to optimize expected values, leaving a gap in utilizing linear programming methods. Therefore, exploring research avenues that employ linear programming techniques can provide valuable insights and contribute to the advancement of this field, where a broad overview of the topic can be found in [17,18] and the references therein. With the rapid advancement of network communication technologies, numerous transmission strategies have emerged to effectively utilize limited network resources and prevent unnecessary waste. Another possible factor to consider is the limitation of network resources, which results in failures or interruptions in continuous communication between the sender and receiver in networks. As research progressed, it was discovered that time-triggering methods [19,20] may not respond immediately to events, potentially missing crucial state changes within fixed time intervals. Consequently, some new event-triggered mechanisms [21,22] were proposed and widely embraced in contemporary research. For general systems, the event-triggered conditions are usually presented in linear matrix inequalities and the quadratic Lyapunov function is commonly used. Since the states of positive systems are non-negative, it is unnecessary to apply the quadratic approach directly to positive systems. It had been shown that a linear approach is more powerful for handling the issues of positive systems [23,24,25]. Therefore, it is of significance to employ linear approach to describe event-triggered conditions of positive complex network systems.
With the motivations as described above, in this paper, the event-triggered problem of discrete-time complex networks using linear programming methods is investigated. Firstly, the research delves into the problem of event-triggered state estimation for discrete-time complex networks. Subsequently, by utilizing matrix decomposition technology, an event-triggered state estimator and controller are designed to address the problem at hand. Finally, the Lyapunov stability theory is employed to guarantee the positivity and stability of the system. The rest of the paper is organized as follows: Section 2 introduces a discrete-time complex network model, Section 3 provides the main results, simulation examples are included in Section 4, and Section 5 concludes the paper.

2. Preliminaries

We will model the observation and control of discrete-time complex networks. To improve readability, the notations used in this paper are standardized, as shown in Table 1.
Consider the dynamical networks with N coupled nodes in the following form:
x i ( k + 1 ) = A i x i ( k ) + j = 1 N ϖ i j T x j ( k ) + E i ω i ( k ) + B i u i ( k ) , y i ( k ) = C i x i ( k ) ,
where for the ith node, x i ( k ) R n , u i ( k ) R n and y i ( k ) R m are the state vector, control input, and measurement output, respectively. ω i ( k ) R n is the external disturbance, and the 1 -norm of the vector ω ( k ) is ω ( k ) 1 = k = 1 ω ( k ) 1 . The notation T = d i a g { t 1 , t 2 , , t n } is the inner coupling matrix between two connected nodes for all 1 i and j N . The notation W = [ ϖ i j ] N × N is the coupled configuration matrix representing the coupling structure of the dynamical network. If there is a connection between node i and j for i j , then ϖ i j 0 , otherwise ϖ i j = 0 . The notations A i R n × n , B i R n × n , C i R m × n and E i R n × n are constant matrices with appropriate dimensions. The notation M i represents the gain of the controller to be designed.
Lemma 1 
([8,26]). System (1) is considered positive if and only if the matrices A i 0 , B i 0 , C i 0 , and E i 0 . Based on the assumptions made for systems (1), it can be concluded that the system is positive.
Lemma 2 
([8,26]). The equivalence between the following conditions holds for a matrix A 0 : (i) The matrix A is a Schur matrix. (ii) There exists some vector v 0 such that ( A I ) v 0 .
Definition 1 
([8,26]). The system (1) is said to be 1 -gain stable if the following statements hold: (i) When ω ( k ) = 0 , the system (1) is asymptotically stable. (ii) Under zero-initial conditions, the following inequality holds for ω ( k ) 0 ,
k = 1 e ( k ) 1 k = 1 γ ω ( k ) 1 ,
where γ > 0 is the 1 -gain value.

3. Main Results

This section will design the event-triggered observer and observer-based controller for positive complex networks.

3.1. Event-Triggered Observer

Construct the event-based state estimator for the ith node without input as
x ^ i ( k + 1 ) = A i x ^ i ( k ) + j = 1 N ϖ i j T x ^ j ( k ) L i C i x ^ i ( k ) + K i [ y i ( k t i ) C i x ^ i ( k ) ] ,
where K i and L i denote the gains of the state estimator to be designed.
Remark 1. 
It is clear that the state estimator gain K i ( k ) given in (2) is designed based on the event-triggered mechanism. However, in general complex network systems the state estimators are designed in the non-event-triggered form [27,28]. In addition, an additional state estimator gain L i ( k ) is introduced in (2) to ensure the positivity of the system. It implies that existing estimation design cannot be applied for positive complex networks. This is also one of novelties of (2).
Let ϱ i ( k ) = y i ( k t i ) y i ( k ) be defined as the difference between the ith sensor’s measurements at the latest triggering instant and the current sampling instant. Then the estimator (2) can be rewritten as
x ^ i ( k + 1 ) = A i x ^ i ( k ) + j = 1 N ϖ i j T x ^ j ( k ) + K i [ C i x i ( k ) C i x ^ i ( k ) + ϱ i ( k ) ] L i C i x ^ i ( k ) .
By denoting e i ( k ) = x i ( k ) x ^ i ( k ) , the estimator (3) and the estimation error dynamics without input can be rewritten as
x ^ i ( k + 1 ) = ( A i L i C i ) x ^ i ( k ) + j = 1 N ϖ i j T x ^ j ( k ) + K i C i e i ( k ) + K i ϱ i ( k ) , e i ( k + 1 ) = ( A i K i C i ) e i ( k ) + j = 1 N ϖ i j T e j ( k ) + E i ω i ( k ) K i ϱ i ( k ) + L i C i x ^ i ( k ) .
For convenience, denote
A = d i a g N { A i } , C = d i a g N { C i } , E = d i a g N { E i } , K N = d i a g N { K i } , L = d i a g N { L i } , e ( k ) = { e 1 ( k ) , e 2 ( k ) , , e N ( k ) } , ω ( k ) = { ω 1 ( k ) , ω 2 ( k ) , , ω N ( k ) } , x ^ ( k ) = { x ^ 1 ( k ) , x ^ 2 ( k ) , , x ^ N ( k ) } , ϱ ( k ) = { ϱ 1 ( k ) , ϱ 2 ( k ) , , ϱ N ( k ) } .
By selecting x ˜ ( k + 1 ) = ( x ^ ( k + 1 ) , e ( k + 1 ) ) and applying the Kronecker product, the system without input can be expressed in the following form:
x ˜ ( k + 1 ) = L C K C L C K C x ^ ( k ) e ( k ) + K K ϱ ( k ) + 0 E ω ( k ) ,
where = A + W T .
Theorem 1. 
If there exist constants 0 < λ 1 < λ 2 < 1 , 0 < μ ̲ < μ < 1 , β > 0 , ϖ i i < 0 , R n vectors υ 2 i υ 1 i 0 , δ i ( ι ) 0 , ζ i ( ι ) 0 , δ i ̲ 0 , δ i ¯ 0 , ζ i 0 such that
1 n υ 2 i A i ι = 1 n 1 n ( ι ) ζ i ( ι ) C i + 1 n υ 2 i ϖ i i T ι = 1 n 1 n ( ι ) δ i ( ι ) β 1 m × m C i 0 ,
ι = 1 n 1 n ( ι ) δ i ( ι ) C i ι = 1 n 1 n ( ι ) δ i ( ι ) β 1 m × m C i 0 ,
ι = 1 n 1 n ( ι ) ζ i ( ι ) C i ι = 1 n 1 n ( ι ) δ i ( ι ) β 1 m × m C i 0 ,
1 n υ 2 i A i ι = 1 n 1 n ( ι ) δ i ( ι ) C i + 1 n υ 2 i ϖ i i T ι = 1 n 1 n ( ι ) δ i ( ι ) β 1 m × m C i 0 ,
A υ 1 + ( 1 λ 1 ) C ζ + ( W T ) υ 1 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ μ υ 1 0 ,
A υ 2 ( 1 λ 2 ) C δ ̲ + ( W T ) υ 2 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ + I N 1 n μ υ 2 0 ,
E υ 2 I N γ 1 m 0 ,
λ 1 υ 2 υ 1 λ 2 υ 2 ,
δ i ̲ δ i ( ι ) δ i ¯ , ζ i ( ι ) ζ i ,
hold for ι = 1 , 2 , , N , then under the event-based state estimator (2) with gain matrices satisfying
K i = ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i , L i = ι = 1 n 1 n ( ι ) ζ i ( ι ) 1 n υ 2 i ,
the resulting closed-loop system (4) is positive and stable.
Remark 2. 
From the studies [15,16,27,28], it can be seen that the state estimator gain is designed based on linear matrix inequalities and the corresponding variables are matrices. However, in this paper, a new gain matrix design is introduced in (14) using a matrix decomposition approach and its gain variables are vectors. This matrix decomposition approach has been applied to the design of observers and controllers of positive systems [10,11,23,24]. Compared with the linear matrix inequality method, the linear programming method is more convenient and more suitable for solving positively complex networks. It should be noticed that the presented design approach can also be extended to other designs of positive complex networks such as compensators, synchronization controllers, and filters.
Proof. 
Using the event-triggering condition, it follows that
ϱ ( k ) 1 β y ( k ) 1 ,
where x i ( k ) 0 implies y i ( k ) 0 through the equation y i ( k ) = C i x i ( k ) . Then, it holds that
β 1 m × m y i ( k ) ϱ i ( k ) β 1 m × m y i ( k ) .
From (4), the matrix form can be readily derived that
x ˜ ( k + 1 ) L C Λ K C Λ L C Λ K C Λ x ^ ( k ) e ( k ) + 0 E ω ( k ) ,
where
L C Λ = Θ 1 ϖ 12 T ϖ 1 N T ϖ 21 T Θ 2 ϖ 2 N T ϖ N 1 T ϖ N 2 T Θ N , K C Λ = d i a g N { i } , K C Λ = P 1 ϖ 12 T ϖ 1 N T ϖ 21 T P 2 ϖ 2 N T ϖ N 1 T ϖ N 2 T P N , L C Λ = d i a g N { Υ i } ,
and Λ = K ( I N β 1 m × m ) C , Θ i = A i L i C i + ϖ i i T K i β 1 m × m C i , i = K i C i K i β 1 m × m C i , Υ i = L i C i K i β 1 m × m C i , P i = A i K i C i + ϖ i i T K i β 1 m × m C i .
From (14), it follows that
Θ i = A i ι = 1 n 1 n ( ι ) ζ i ( ι ) 1 n υ 2 i C i + ϖ i i T ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i β 1 m × m C i , i = ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i C i ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i β 1 m × m C i , Υ i = ι = 1 n 1 n ( ι ) ζ i ( ι ) 1 n υ 2 i C i ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i β 1 m × m C i , P i = A i ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i C i + ϖ i i T ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i β 1 m × m C i .
With the fact that ϖ 1 i T 0 and utilize Equations (5)–(8), it follows that L C Λ 0 , K C Λ 0 , L C Λ 0 and K C Λ 0 . Thus, it can be inferred that system (4) is positive.
Next, the stability of system (4) is addressed. From (4) and (16), the system can be rewritten as
x ˜ ( k + 1 ) Z + Λ K C + Λ L C + Λ Ξ + Λ x ^ ( k ) e ( k ) + 0 E ω ( k ) ,
where Z = A L C + W T , Ξ = A K C + W T .
Consider the choice of a positive Lyapunov function V ( x ( k ) ) = x ˜ ( k ) υ , where υ = ( υ 1 , υ 2 ) , λ 1 υ 2 υ 1 λ 2 υ 2 , υ 1 = ( υ 11 , υ 12 , υ 1 N ) , υ 2 = ( υ 21 , υ 22 , υ 2 N ) and υ 1 i R n , υ 2 i R n . Considering the case ω ( k ) = 0 , it is not hard to derive that
V ( k ) = V ( x ( k + 1 ) ) V ( x ( k ) ) x ^ ( k ) e ( k ) A υ 1 + C L ( υ 2 υ 1 ) + ( W T ) υ 1 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) υ 1 A υ 2 + C K ( υ 1 υ 2 ) + ( W T ) υ 2 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) υ 2 .
From (13) and (14), it holds that 1 n δ i ̲ 1 n υ 2 i K i 1 n δ i ¯ 1 n υ 2 i , L i 1 n ζ i 1 n υ 2 i . From (19), it follows that
V ( k ) + ( 1 μ ) V ( k ) x ^ ( k ) e ( k ) A υ 1 + ( 1 λ 1 ) C ζ + ( W T ) υ 1 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ μ υ 1 A υ 2 ( 1 λ 2 ) C δ ̲ + ( W T ) υ 2 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ μ υ 2 .
With the facts that x ^ ( k ) 0 , e ( k ) 0 in mind, using (9) and (10) leads to V ( K ) + ( 1 μ ) V ( K ) 0 . Consequently, it can be obtained that V ( k + 1 ) μ V ( k ) . As a result, the system (4) is exponentially stable.
Considering the case ω ( k ) 0 , the Equation (19) can be rewritten as
V ( k ) x ^ ( k ) e ( k ) A υ 1 + C L ( υ 2 υ 1 ) + ( W T ) υ 1 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) υ 1 A υ 2 + C K ( υ 1 υ 2 ) + ( W T ) υ 2 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) υ 2 + ω ( k ) E υ 2 .
Define ψ ( k ) = γ ω ( k ) 1 e ( k ) 1 . Then the inequality holds that
V ( k ) + ( 1 μ ) V ( k ) x ^ ( k ) e ( k ) A υ 1 + ( 1 λ 1 ) C ζ + ( W T ) υ 1 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ μ υ 1 A υ 2 ( 1 λ 2 ) C δ ̲ + ( W T ) υ 2 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ μ υ 2 + ω ( k ) ( E υ 2 I N γ 1 m ) .
Applying (9)–(11), it yields that V ( k ) + ( 1 μ ) V ( k ) + ψ ( k ) 0 . Therefore, it can be guaranteed that V ( K + 1 ) μ V ( K ) + ψ ( K ) . Then, it yields that
V ( k + 1 ) V ( 0 ) ( μ 1 ) ( V ( k ) + V ( k 1 ) + + V ( 0 ) ) + k = 1 γ ω ( k ) 1 k = 1 e ( k ) 1 .
With the fact that V ( k + 1 ) 0 , V ( 0 ) = 0 , μ 1 0 , the inequality (23) can be rewritten as k = 1 e ( k ) 1 k = 1 γ ω ( k ) 1 . This implies that the system exhibits 1 -gain stability with a performance level of γ . □
Remark 3. 
Different from the state estimation of general complex networks [5], the positivity of state observer and the stability of the error dynamics are required for the state estimation of positive complex networks. In Theorem 1, a novel observer is proposed for positive complex networks by virtue of matrix decomposition approach. Due to the property that all states of positive systems are non-negative, quadratic Lyapunov functions and linear matrix inequalities may no longer be the optimal choice [13,14,15,16]. A co-positive Lyapunov function V ( x ( k ) ) = x ˜ ( k ) υ used in Theorem 1 is presented to analyze the stability of the corresponding systems. In contrast to the linear matrix inequalities, the conditions (5)–(11) and (26)–(28) are given in the form of linear programming. It is less difficult to calculate and more convenient to deal with the problems related to large-scale systems. It should be pointed out that the presented state estimation framework in Theorem 1 can be developed for other issues of positive complex networks.

3.2. Observer-Based Control

Design an observer-based controller: u i ( k ) = M i x ^ i ( k ) , where M i is the control gain. Conduct an event-based state estimator for the ith node:
x ^ i ( k + 1 ) = A i x ^ i ( k ) + j = 1 N ϖ i j T x ^ j ( k ) L i C i x ^ i ( k ) + B i M i x ^ ( k ) + K i [ y i ( k t i ) C i x ^ i ( k ) ] .
Then the system with input can be expressed in the following form:
x ˜ ( k + 1 ) = L C + B M K C L C K C x ^ ( k ) e ( k ) + K K ϱ ( k ) + 0 E ω ( k ) .
Theorem 2. 
If there exist constants 0 < λ 1 < λ 2 < 1 , 0 < μ ̲ < μ < 1 , β > 0 , ϖ i i < 0 , φ > 0 , R n vectors υ 2 i υ 1 i 0 , δ i ( ι ) 0 , ζ i ( ι ) 0 , δ i ̲ 0 , δ i ¯ 0 , ζ i 0 , ε i ( ι ) 0 , ε i 0 such that
1 n B i υ 2 i A i φ ι = 1 n 1 n ( ι ) ζ i ( ι ) C i + B i ι = 1 n 1 n ( ι ) ε i ( ι ) + 1 n B i υ 2 i ϖ i i T φ ι = 1 n 1 n ( ι ) δ i ( ι ) β 1 m × m C i 0 ,
A υ 1 + ( 1 λ 1 ) C ζ + ( W T ) υ 1 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ + λ 1 ε μ υ 1 0 ,
A υ 2 ( 1 λ 2 ) C δ ̲ + ( W T ) υ 2 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ + I N 1 n μ υ 2 0 ,
1 n B i υ 2 i φ 1 n υ 2 i ,
ε i ( ι ) ε i ,
hold for ι = 1 , 2 , , N , then the event-based state estimator with gain matrix the same as (14) and the event-triggered control law with gain matrix satisfy
M i = ι = 1 n 1 n ( ι ) ε i ( ι ) 1 n B υ 2 i ,
the resulting closed-loop system (25) is positive and stable.
Proof. 
From (25), the matrix is straightforward to obtain
x ˜ ( k + 1 ) Ξ Λ K C Λ L C Λ K C Λ x ^ ( k ) e ( k ) + 0 E ω ( k ) ,
where Ξ = A L C + B M + W T , and Λ , , K C Λ , L C Λ , K C Λ are the same as (17).
Ξ Λ = T 1 ϖ 12 T ϖ 1 N T ϖ 21 T T 2 ϖ 2 N T ϖ N 1 T ϖ N 2 T T N ,
with T i = A i L i C i + ϖ i i T + B i M i K i β 1 m × m C i .
From (14) and (31), it follows that T i = A i ι = 1 n 1 n ( ι ) ζ i ( ι ) 1 n υ 2 i C i + ϖ i i T + B i ι = 1 n 1 n ( ι ) ε i ( ι ) 1 n B υ 2 i ι = 1 n 1 n ( ι ) δ i ( ι ) 1 n υ 2 i β 1 m × m C i .
With the fact that ϖ 1 i T 0 , by utilizing Equation (26), it follows that Ξ Λ 0 and from (6)–(8), it obtains that K C Λ 0 , L C Λ 0 and K C Λ 0 . Thus, it is not hard to conclude that system (25) is positive.
Next, it comes time to prove the stability of system (25). By using (25) and (16), the system can be rewritten as
x ˜ ( k + 1 ) Ξ + Λ K C + Λ L C + Λ K C + Λ x ^ ( k ) e ( k ) + 0 E ω ( k ) ,
Considering the case ω ( k ) = 0 , it derives that
V ( k ) x ^ ( k ) e ( k ) Π 1 Π 2 ,
where Π 1 = A υ 1 + C L ( υ 2 υ 1 ) + ( W T ) υ 1 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) + M B υ 1 υ 1 , and Π 2 = A υ 2 + C K ( υ 1 υ 2 ) + ( W T ) υ 2 + C ( I N β 1 m × m ) K ( υ 1 + υ 2 ) υ 2 . In addition, from (30) it drives that M i 1 n ε i 1 n B i υ 2 i . Thus Π 1 A υ 1 + ( 1 λ 1 ) C ζ + ( W T ) υ 1 + ( 1 + λ 2 ) C ( I N β 1 m × m ) δ ¯ + λ 1 ε μ υ 1 .
With the facts that x ^ ( k ) 0 , e ( k ) 0 in mind, it can be inferred that V ( k + 1 ) μ V ( k ) . As a result, the system (25) is exponentially stable.
Considering the case ω 0 , and define ψ ( k ) = γ ω ( k ) 1 e ( k ) 1 . Then the following inequality holds as
V ( k ) + ( 1 μ ) V ( k ) ψ ( k ) x ^ ( k ) e ( k ) Π 1 Π 2 + I N 1 n + ω ( k ) ( E υ 2 I N γ 1 m ) .
With the fact that V ( k + 1 ) 0 , V ( 0 ) = 0 , μ 1 0 , the inequality (35) can be rewritten as k = 1 e ( k ) 1 k = 1 γ ω ( k ) 1 . As a result, the system is 1 -gain stable with performance γ . □
Remark 4. 
Theorem 1 provides a design scheme for the gain matrices K i and L i of the state estimator, while Theorem 2 outlines the design method for the gain matrix M i of the controller, building upon the foundation established in Theorem 1. It is worth noting that, in many cases, complex networks are commonly addressed using quadratic forms, as evident in references [5,6,7]. However, such an approach is not practical for positive complex networks. Therefore, it becomes crucial to carefully select appropriate methods for designing the state estimator and controller in order to ensure their effectiveness in positive complex networks.

4. Numerical Example

In this section, a numerical simulation example is given to demonstrate the validity of the proposed design. Consider system (1), where
A 1 = 0.34 0.36 0.33 0.34 0.32 0.34 0.32 0.34 0.33 , A 2 = 0.33 0.36 0.32 0.34 0.32 0.36 0.32 0.35 0.35 , A 3 = 0.32 0.36 0.36 0.34 0.33 0.34 0.33 0.35 0.33 , B 1 = 0.02 0.01 0.03 0.01 0.03 0.02 0.03 0.01 0.01 , B 2 = 0.01 0.01 0.01 0.03 0.03 0.02 0.02 0.02 0.01 , B 3 = 0.03 0.02 0.01 0.02 0.03 0.02 0.01 0.02 0.03 , C 1 = 0.01 0.03 0.02 0.02 0.03 0.01 , C 2 = 0.02 0.02 0.01 0.01 0.03 0.01 , C 3 = 0.03 0.01 0.02 0.02 0.02 0.03 , E 1 = 0.4 0.2 0.2 0.1 0.3 0.3 , E 2 = 0.3 0.3 0.2 0.1 0.1 0.2 , E 3 = 0.2 0.1 0.2 0.2 0.3 0.3 , W = 0.02 0.01 0.01 0.01 0.02 0.01 0.01 0.01 0.02 , T = 1 0 0 0 1 0 0 0 1 .
Given the values β = 0 , 13 , φ = 0.08 , λ 1 = 0.4 , λ 2 = 0.5 , μ = 0.98 , γ = 354.9149 and the disturbance signal ω ( k ) = ( 20 × 1 ( 1 + k ) 6 / 7 ) , the event-based state estimator with gain matrices satisfy:
K 1 = 0.8744 4.2970 0.8740 4.2967 0.8743 4.2970 , K 2 = 5.5099 0.0669 5.5098 0.0668 5.5099 0.0699 , K 3 = 0.0767 6.1071 0.0764 6.1068 0.0769 6.1073 , L 1 = 0.6965 0.6970 0.6882 0.6890 0.6950 0.6953 , L 2 = 0.9133 0.6093 0.9105 0.6104 0.9134 0.6093 , L 3 = 0.7842 0.8658 0.7915 0.8502 0.7947 0.8635 ,
and the event-triggered control law with gain matrices:
M 1 = 3.7188 3.1974 3.5826 3.8179 3.1867 3.6875 3.6253 3.1845 3.5325 , M 2 = 3.2985 2.5194 3.0847 3.2985 2.5194 3.0847 1.5763 1.6570 1.5947 , M 3 = 3.2569 3.5010 3.3832 6.9703 5.3239 6.5186 3.3310 3.5015 3.3699 ,
Figure 1 illustrates the event-triggering signal. Unlike the time-triggering strategy, the number of the event-triggering strategy is smaller than the number of the time-triggering strategy at the same instant, which reduces unnecessary communication and avoids a waste of resources. Additionally, Figure 2, Figure 3 and Figure 4 display the initial states of x and the corresponding state estimates. It can be found that the state of the estimator converges with the state of the system node. This shows that the state estimator designed in this paper has a good effect on the system. Figure 5 presents the states of the error system and all three subsystems eventually converge to 0. It reaches the positivity and stability of the systems.

5. Conclusions

This paper designs an event-triggered state estimator and an estimator-based controller for discrete-time complex network systems with disturbances. The event-triggered conditions for discrete complex networks are formulated using the 1-norm. Matrix factorization techniques are utilized to design the gain matrices for the state estimator and the estimator-based controller. The positivity and stability of the lower bound and upper bound of the system are rigorously proven through the use of linear positive Lyapunov functions and linear programming, thereby ensuring the positivity and stability of the L 1 -gain of the system. Future research could further explore the identification of key nodes and critical paths in the networks to improve the precision of control. It could also include how nonlinear interactions in the network affect system stability and how these interactions can be exploited to achieve precise guidance of network behavior. Additionally, some possible research can be devoted to understanding how small interventions guide the network from one state to another.

Author Contributions

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

Funding

This work was supported by Hainan Provincial Natural Science Foundation of China (622RC618), Science Research Funding of Hainan University (KYQD(ZR)21007) and Postgraduate Innovative Research Funding of Hainan Province (Qhys2023-279, Qhys2023-278 and Qhys2023-280).

Data Availability Statement

No new data were created or analyzed in this study.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Gupta, O.P.; Deshmukh, R.; Kumar, A.; Singh, S.K.; Sharma, P.; Ram, S.; Singh, G.P. From gene to biomolecular networks: A review of evidences for understanding complex biological function in plants. Curr. Opin. Biotechnol. 2022, 74, 66–74. [Google Scholar] [CrossRef] [PubMed]
  2. Long, H.; Ci, J.; Guo, Z.; Wen, S.; Huang, T. Synchronization of coupled switched neural networks subject to hybrid stochastic disturbances. Neural Netw. 2023, 166, 459–470. [Google Scholar] [CrossRef] [PubMed]
  3. Sethi, A.; Melamud, E. Joint inference of physiological network and survival analysis identifies factors associated with aging rate. Cell Rep. Methods 2022, 2, 100356. [Google Scholar] [CrossRef] [PubMed]
  4. Jia, C.; Zhang, R.; Wang, D. Evolutionary game of cooperative behavior among social capitals in PPP projects: A complex network perspective. Ain Shams Eng. J. 2023, 14, 102006. [Google Scholar] [CrossRef]
  5. Liu, Y.; Shen, B.; Zhang, P. Synchronization and state estimation for discrete-time coupled delayed complex-valued neural networks with random system parameters. Neural Netw. 2022, 150, 181–193. [Google Scholar] [CrossRef]
  6. Lu, R.; Yu, W.; Lü, J.; Xue, A. Synchronization on complex networks of networks. IEEE Trans. Neural Netw. Learn. Syst. 2014, 25, 2110–2118. [Google Scholar] [CrossRef]
  7. Liu, Y.; Wang, Z.; Liang, J.; Liu, X. Synchronization and state estimation for discrete-time complex networks with distributed delays. IEEE Trans. Syst. Man Cybern. Part B 2008, 38, 1314–1325. [Google Scholar]
  8. Farina, L.; Rinaldi, S. Positive Linear Systems: Theory and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2000. [Google Scholar]
  9. Zhao, P.; Kang, Y.; Niu, B.; Zhao, Y. Input-to-state stability and stabilization for switched nonlinear positive systems. Nonlinear Anal. Hybrid Syst. 2023, 47, 101298. [Google Scholar] [CrossRef]
  10. Zhang, J.; Zheng, G.; Feng, Y.; Chen, Y. Event-triggered state-feedback and dynamic output-feedback control of PMJSs with intermittent faults. IEEE Trans. Autom. Control 2022, 68, 1039–1046. [Google Scholar] [CrossRef]
  11. Su, H.; Wu, H.; Chen, X.; Chen, M.Z. Positive edge consensus of complex networks. IEEE Trans. Syst. Man Cybern. Syst. 2017, 48, 2242–2250. [Google Scholar] [CrossRef]
  12. Li, Q.; Wang, Z.; Sheng, W.; Alsaadi, F.E.; Alsaadi, F.E. Dynamic event-triggered mechanism for H non-fragile state estimation of complex networks under randomly occurring sensor saturations. Inf. Sci. 2020, 509, 304–316. [Google Scholar] [CrossRef]
  13. Wang, L.; Wang, Z.; Huang, T.; Wei, G. An event-triggered approach to state estimation for a class of complex networks with mixed time delays and nonlinearities. IEEE Trans. Cybern. 2015, 46, 2497–2508. [Google Scholar] [CrossRef] [PubMed]
  14. Liu, Y.; Wang, Z.; Ma, L.; Alsaadi, F.E. A partial-nodes-based information fusion approach to state estimation for discrete-time delayed stochastic complex networks. Inf. Fusion 2019, 49, 240–248. [Google Scholar] [CrossRef]
  15. Wu, Z.; Chen, J.; Zhang, X.; Xiao, Z.; Tao, J.; Wang, X. Dynamic event-triggered synchronization of complex networks with switching topologies: Asynchronous observer-based case. Appl. Math. Comput. 2022, 435, 127413. [Google Scholar] [CrossRef]
  16. Zhao, Y.; Li, X.; Rao, R. Synchronization of nonidentical complex dynamical networks with unknown disturbances via observer-based sliding mode control. Neurocomputing 2021, 454, 441–447. [Google Scholar] [CrossRef]
  17. Zhou, X.; Zhang, J.; Lin, P.; Shi, X. Event-triggered PI control of positive switched systems. In Proceedings of the 2022 41st Chinese Control Conference (CCC), Hefei, China, 25–27 July 2022; IEEE: Piscataway, NJ, USA, 2022; pp. 1593–1598. [Google Scholar]
  18. Liu, L.J.; Karimi, H.R.; Zhao, X. New approaches to positive observer design for discrete-time positive linear systems. J. Frankl. Inst. 2018, 355, 4336–4350. [Google Scholar] [CrossRef]
  19. Kirner, R.; Puschner, P. A qualitative cybersecurity analysis of time-triggered communication networks in automotive systems. J. Syst. Archit. 2023, 136, 102835. [Google Scholar] [CrossRef]
  20. Zhang, X.; Wang, Q.; Fu, B. Further stabilization criteria of continuous systems with aperiodic time-triggered intermittent control. Commun. Nonlinear Sci. Numer. Simul. 2023, 125, 107387. [Google Scholar] [CrossRef]
  21. Liu, Y.; Shen, B.; Shu, H. Finite-time resilient H state estimation for discrete-time delayed neural networks under dynamic event-triggered mechanism. Neural Netw. 2020, 121, 356–365. [Google Scholar] [CrossRef]
  22. Liu, C.; Li, H.; Shi, Y.; Xu, D. Codesign of event trigger and feedback policy in robust model predictive control. IEEE Trans. Autom. Control 2019, 65, 302–309. [Google Scholar] [CrossRef]
  23. Zhang, J.F.; Liu, L.Y.; Fu, S.Z.; Li, S. Event-triggered control of positive switched systems with actuator saturation and time-delay. Int. J. Autom. Comput. 2021, 18, 141–154. [Google Scholar] [CrossRef]
  24. Zhang, S.; Zhang, J.; Li, S.; Yu, S. Event-triggered control of positive semi-markovian jump systems without/with input saturation. J. Shanghai Jiaotong Univ. 2022, 27, 723–736. [Google Scholar] [CrossRef]
  25. Xiao, S.; Zhang, Y.; Zhang, B. Event-triggered networkbased state observer design of positive systems. Inf. Sci. 2018, 469, 30–43. [Google Scholar] [CrossRef]
  26. Kaczorek, T. Positive 1D and 2D Systems; Springer: London, UK, 2001. [Google Scholar]
  27. Shen, B.; Wang, Z.; Liu, X. Bounded H synchronization and state estimation for discrete time-varying stochastic complex networks over a finite horizon. IEEE Trans. Neural Netw. 2010, 22, 145–157. [Google Scholar] [CrossRef]
  28. Shen, B.; Wang, Z.; Ding, D.; Shu, H. H state estimation for complex networks with uncertain inner coupling and incomplete measurements. IEEE Trans. Neural Netw. Learn. Syst. 2013, 24, 2027–2037. [Google Scholar] [CrossRef]
Figure 1. Event-triggering signal.
Figure 1. Event-triggering signal.
Mathematics 12 02971 g001
Figure 2. Observation signal and status signal of x 1 .
Figure 2. Observation signal and status signal of x 1 .
Mathematics 12 02971 g002
Figure 3. Observation signal and status signal of x 2 .
Figure 3. Observation signal and status signal of x 2 .
Mathematics 12 02971 g003
Figure 4. Observation signal and status signal of x 3 .
Figure 4. Observation signal and status signal of x 3 .
Mathematics 12 02971 g004
Figure 5. The status of the error system.
Figure 5. The status of the error system.
Mathematics 12 02971 g005
Table 1. Notations.
Table 1. Notations.
NotationsExpression
R n n-dimensional Euclidean space
R n × m n × m real matrices
a i j The element in the ith row and jth column of matrix A
A 0 ( 0 ) All elements in A are non-negative (or positive)
A The transposition of matrix A
d i a g N { A i } The block-diagonal matrix consisting of A 1 , A 2 , A N
1 m × m R m × m An m × m -dimensional matrix with all elements being 1
1 n An n-dimensional matrix with all elements being 1
I N The N × N identity matrix
1 n ( ι ) 1 n ( ι ) = ( 0 , , 0 ι 1 ) , 1 , ( 0 , , 0 n ι )
The Kronecker product
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

Zhang, Y.; Wu, Y.; Sun, Y.; Zhang, P. Estimation and Control of Positive Complex Networks Using Linear Programming. Mathematics 2024, 12, 2971. https://doi.org/10.3390/math12192971

AMA Style

Zhang Y, Wu Y, Sun Y, Zhang P. Estimation and Control of Positive Complex Networks Using Linear Programming. Mathematics. 2024; 12(19):2971. https://doi.org/10.3390/math12192971

Chicago/Turabian Style

Zhang, Yan, Yuanyuan Wu, Yishuang Sun, and Pei Zhang. 2024. "Estimation and Control of Positive Complex Networks Using Linear Programming" Mathematics 12, no. 19: 2971. https://doi.org/10.3390/math12192971

APA Style

Zhang, Y., Wu, Y., Sun, Y., & Zhang, P. (2024). Estimation and Control of Positive Complex Networks Using Linear Programming. Mathematics, 12(19), 2971. https://doi.org/10.3390/math12192971

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