Next Article in Journal
Machine Tool Wear Prediction Technology Based on Multi-Sensor Information Fusion
Next Article in Special Issue
Anomaly Detection for Asynchronous Multivariate Time Series of Nuclear Power Plants Using a Temporal-Spatial Transformer
Previous Article in Journal
Vehicle-Type Recognition Method for Images Based on Improved Faster R-CNN Model
Previous Article in Special Issue
Tool Condition Monitoring Using Machine Tool Spindle Current and Long Short-Term Memory Neural Network Model Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults

College of Electrical Engineering and Control Science, Nanjing Tech University, Nanjing 211816, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(8), 2651; https://doi.org/10.3390/s24082651
Submission received: 11 March 2024 / Revised: 6 April 2024 / Accepted: 18 April 2024 / Published: 21 April 2024
(This article belongs to the Special Issue AI-Assisted Condition Monitoring and Fault Diagnosis)

Abstract

:
This article investigates the robust cooperative fault-tolerant control problem of multi-agent systems subject to mismatched uncertainties and actuator faults. During the design process of the intermediate variable estimator, there is no need to satisfy fault estimation matching conditions, and this overcomes a crucial constraint of traditional observers and estimators. The feedback term of the designed estimator contains the centralized estimation errors and the distributed estimation errors of the agent, and this further improves the design freedom of the proposed estimator. A novel fault-tolerant control protocol is designed based on the fault estimation information. In this work, the bounds of the fault and its derivatives are unknown, and the considered method is applicable to both directed and undirected multi-agent systems. Furthermore, the parameters of the estimator are determined through the resolution of a linear matrix inequality (LMI), which is decoupled by employing coordinate transformation and Schur decomposition. Lastly, a numerical simulation result is used to demonstrate the effectiveness of the proposed method.

1. Introduction

Spawned by the rapid advances in networked systems and distributed cooperation, there has been a flurry of activity on the topic of multi-agent systems (MASs). During the past few years, MASs have found extensive application in both industrial and military domains, such as multi-robot systems [1], satellite formation control [2], sensor networks [3], smart grids [4], flight control systems [5], etc.
As a type of interconnected system, MASs complete control tasks by exchanging information among neighboring agents. However, due to the system configuration and communication topology, a fault occurring in one agent can potentially propagate through the network, affect other agents, or even cause network disruption [6,7,8]. As a result, one specific agent fault could significantly impact MASs, potentially causing agent instability or even a system crash. Indeed, it is essential to highlight that faults may occur more frequently as the number of agents increases and as each agent expands in size and complexity. Furthermore, the security and reliability of MASs become increasingly crucial when they operate in complex and adverse environments.
To ensure the reliable operation of control systems, fault diagnosis has gained widespread attention in recent years and become a key focus of research. Fault diagnosis can be categorized into three main themes: fault detection (FD), fault isolation (FI), and fault estimation (FE) [9,10]. FD and FI primarily serve to detect and determine the location of faults, but they cannot obtain accurate fault information. In contrast, FE can provide detailed information about the shape and size of faults [11]. In general, the more fault information that can be obtained, the higher the fault-tolerant control performance [12]. Recently, a lot of fault estimation methods have been developed, including robust estimators, adaptive observers, sliding-mode controllers, and others [13,14]. However, there has been relatively limited research on fault estimation for MASs. For instance, ref. [15] addressed distributed fault estimation for interconnected systems with actuator faults. A novel distributed FD and FI method was proposed in [16] for heterogeneous MASs. To estimate sensor faults, a category of unknown input observers (UIO) was introduced in [17] by incorporating the system state and fault into a new extended vector. A robust fault estimation method using sliding-mode observers was derived in [18] for linear MASs subject to actuator faults. For the case of nonlinear MASs, an intermediate observer was introduced in [19] for estimating system states and actuator faults. Similar approaches were applied in [20] to construct nonlinear fault estimation observers for nonlinear MASs with sensor faults. In [21], timing FE observers were developed for interconnected systems with both multiplicative and additive actuator faults.
It is noteworthy that the majority of fault estimation methods in the aforementioned studies necessitate knowledge of the upper bound of the fault signal. However, since the bounds of fault information are mostly unknown for most control systems, these methods tend to be conservative in terms of fault estimation performance and applicability. Furthermore, in some existing works, the fault signal is designed as a constant, which is unrealistic since faults occurring in dynamic systems are inevitably influenced by environmental conditions and exhibit dynamic characteristics. Hence, these issues must be taken into consideration in the context of fault estimation for MASs.
Additionally, the fault estimation matching condition, which states that the rank of the product of the system output matrix and the fault distribution matrix should equal the rank of the system output matrix, is necessary for most existing fault estimation protocols. The aforementioned condition exhibits a higher level of conservatism when compared to the strict positivity assumption. In response to these challenges, a pioneering fault estimation method grounded in an intermediate observer was introduced in [22]. Notably, this approach dispenses with the necessity for both the fault estimation matching condition and the strictly positive real assumption.
Fault-tolerance control is the follow-up work of fault estimation. In the study of fault-tolerant control for MASs, there has been significant attention placed on observer-based active fault-tolerant control algorithms. In the process of designing fault-tolerant controllers, understanding the magnitude of the fault is crucial, and this information must be estimated by the constructed observers. For a class of MASs under a directed fixed topology, ref. [23] sidestepped the restriction of the requirement of zero initial conditions in most existing works and pioneered a novel delay-dependent fault-tolerant controller that solves the fault-tolerant constrained consensus problem of MASs in the presence of communication delays and actuator faults. Similarly, in order to sidestep the limitation of zero initial conditions in the H control method, by using disturbance estimation information and a new performance function, ref. [24] designed a disturbance rejection adaptive fault-tolerant constrained consensus algorithm and solved the constraint consensus problem of perturbation-resilient adaptive fault-tolerant multi-agent systems under actuator failures. For linear and Lipschitz non-linear systems, an adaptive fault-tolerant controller utilizing the compensated actuator fault estimation information was proposed in [25] to mitigate adverse effects on consensus tracking. The problem of distributed fault-tolerant control for linear systems was explored in [26]. This approach utilized fuzzy logic systems to approximate unknown non-linear functions and implemented local observers to estimate the system states. In [27], a form of distributed adaptive fault-tolerant control strategy was proposed for linear MASs. This strategy aimed to alleviate the adverse impacts of actuator failures and losses in actuator effectiveness.
However, precise system models are necessary in the aforementioned works. Actually, model uncertainties are inevitable for practical systems and can stem from unknown internal or external noise, environmental influences, poor plant knowledge, and uncertain/slowly varying parameters. For multi-area power systems with actuator failures, ref. [28] considered the influence of uncertainty factors in the fault-tolerant controller and designed a fault-tolerant control scheme considering both multiplicative perturbation and additive perturbations. In [29], radial basis function neural networks (RBFNNs) were used to estimate dynamic uncertainties in the system model to ensure the accuracy of formation control. Similarly, the fuzzy logic method was used in [30] to address the problem of dynamic uncertainty estimation and further improve the accuracy of the proposed formation control algorithm. In [31], a distributed fault diagnosis method was created for a category of uncertain MASs with actuator faults. There was a need to design a state feedback controller ahead of the observer, suggesting the assumption that the full-dimensional state of the system can be measured. This does not align with the conditions observed in practical systems. A fault diagnosis algorithm based on an unknown input observer (UIO) was proposed in [32] for linear uncertain MASs. This work is only applicable to undirected MASs.
Unfortunately, most existing works have paid limited attention to fault estimation for MASs with uncertainties. Moreover, in the majority of these works, especially in the context of SMOs and UIOs, it is essential to fulfill the fault estimation matching condition and determine the bounds of the fault signal. Undoubtedly, these constraints greatly reduce the practicality and generality of traditional fault estimation and fault-tolerant control approaches.
Inspired by the above discussions, in this paper, the fault estimation and fault-tolerant control problem is addressed for a class of MASs with uncertainties. The innovations of this paper are as follows:
(1)
A type of intermediate variable observer has been designed for each agent to estimate the fault and system state information. The fault estimation matching condition is not necessary for this approach, which is also applicable to time-varying faults.
(2)
In the observer design process, both the centralized and distributed estimation errors of the agents are considered, which has the advantages of centralized and distributed structures and can enhance design flexibility and improve estimation performance.
(3)
In the current work, it is not necessary to obtain the bounds of faults and their derivatives. Consequently, the distributed fault estimation method proposed in this paper has great generality and practicality.
(4)
Compared with most existing results, this FE and FTC scheme is suitable for both directed and undirected MASs.
The structure of the paper is as follows. Section 2 gives some basic assumptions and describes the problem considered in this paper. Section 3 and Section 4 present the main theoretical results of this paper, including the construction of the intermediate variable observers and the convergence analysis. Section 5 presents a simulation example to demonstrate the effectiveness of the proposed method. Finally, the concluding remarks are given in Section 6.
Notations: In this article, R n is the n-dimensional Euclidean spaces and I N is an identity matrix of size N × N . For a vector s R n , s is the Euclidean norm. For a matrix D , H e [ D ] = D + D T , λ m a x ( D ) and λ m i n ( D ) represent the maximum and minimum eigenvalues of a matrix D , respectively.

2. Problem Statement

In this section, the problem formulation is presented. Consider the following linear MASs with unknown mismatched uncertainties
x ˙ i ( t ) = ( A + Δ A ( t ) ) x i ( t ) + B u i ( t ) + E f i ( t ) y i ( t ) = C x i ( t )
where x i ( t ) R n is the state information of agent i, and u i ( t ) R m and y i ( t ) R q are the control input and measured output of the system. f i ( t ) R r is the fault signal and represents actuator faults when E = B . Δ A ( t ) is the perturbed matrix that satisfies Δ A ( t ) = M F ( t ) N , where F T ( t ) F ( t ) I . A, B, E, and C are the known constant matrices of the system, which have appropriate dimensions. Without loss of generality, in this paper, we assume that ( A , C ) is observable and ( A , B ) is stabilizable.
Assumption 1.
In this paper, the fault f i ( t ) satisfies f i ˙ ( t ) θ with θ 0 .
Assumption 2.
r a n k E = r .
Assumption 3.
The following equation holds with respect to every complex number s
R a n k A s I E C 0 = n + r a n k E
Assumption 4.
R a n k B , E = r a n k B .
Remark 1.
Assumption 1 gives the L-2 norm bounds of the fault and its derivative, which implies that the fault and its derivative are energy-bounded. This assumption is common in the field of FD and FTC [22,33] for estimating time-varying signals. In this paper, it is not necessary to know the specific information about the fault and its derivative bounds, i.e., θ is unknown, making the proposed approach more general than most traditional observers [13,14,17], in which the bounds of faults and their first derivatives must be known.
Remark 2.
It is important to note that Assumption 2 is more common than the fault estimation matching conditions proposed in many existing works. Assumption 3 is natural and pervasive in the majority of published articles that investigate fault detection, isolation, and estimation. This assumption implies that the system ( A , C , E ) has a constant amount of zeros in the left half-plane. Such a condition finds frequent application in the realms of system control and fault diagnosis.
Remark 3.
Assumption 4 is very common in the existing results of fault-tolerant control, implying that the fault is situated within the channel responsible for controlling the system’s input, and there is a likelihood that it can be mitigated through compensation via the control input.
The following lemma is used for the subsequent work:
Lemma 1
((Young’s Inequality) [19]). For any two vectors, a R n and b R n , we have
a T b 1 p α p a p + 1 q α q a q
where α > 0 , p > 0 , q > 0 , and p q = p + q .
In this paper, we design a fault estimator to acquire state information and fault information in real time and use the acquired information to design a fault-tolerant control protocol to compensate for the adverse effects of fault signals on MASs for the purpose of robust fault-tolerant cooperative control for MASs.

3. Intermediate Observer Design

In this subsection, an observer is designed for each agent i. To start with, the following intermediate variable ξ i ( t ) is denoted:
ξ i ( t ) = f i ( t ) ω E T x i ( t )
where ω is an intermediate constant selected based on experience.
From Equations (1) and (2), we have the dynamics of the intermediate variable:
ξ ˙ i ( t ) = f ˙ i ( t ) ω E T ( ( A + Δ A ) x i ( t ) + B u i ( t ) + E ξ i ( t ) + ω E E T x i ( t ) )
Based on Equations (1)–(3), the following intermediate variable estimator is designed:
x ^ ˙ i ( t ) = A x ^ i ( t ) + B u i ( t ) + E f ^ i ( t ) + ρ 1 L 1 Δ 1 i + ρ 2 L 2 Δ 2 i ξ ^ ˙ i ( t ) = ω E T E ξ ^ i ( t ) ω E T ( A x ^ i ( t ) + B u i ( t ) + ω E E T x ^ i ( t ) ) f ^ i ( t ) = ξ ^ i ( t ) + ω E T x ^ i ( t ) y ^ i ( t ) = C x ^ i ( t )
where x ^ i ( t ) , ξ ^ i ( t ) , f ^ i ( t ) , and y ^ i ( t ) are the estimations of x i ( t ) , ξ i ( t ) , f i ( t ) , and y i ( t ) , respectively.
In addition, Δ 1 i and Δ 2 i are defined as
Δ 1 i = y i ( t ) y ^ i ( t )
Δ 2 i = j = 1 N a i j [ ( y i ( t ) y ^ i ( t ) ) ( y j ( t ) y ^ j ( t ) ) ]
which represent the centralized and distributed output estimation errors. L 1 R n and L 2 R n are the observer gain matrices, which are designed later. The non-negative constants ρ 1 and ρ 2 are the corresponding weight values that satisfy 0 ρ 1 1 , 0 ρ 2 1 , and ρ 1 + ρ 2 = 0 .
Remark 4.
According to the designed intermediate variable estimator, it can be observed that both the centralized output estimation errors, Δ 1 i , and distributed output estimation errors, Δ 2 i , are taken into account. Here, ρ 1 and ρ 2 represent the weights assigned to the centralized and distributed output estimation errors, respectively. The magnitude of these values signifies the extent to which neighboring nodes influence the observer: a smaller ρ 1 (resulting in a larger ρ 2 ) amplifies the impact of neighboring nodes, whereas a larger ρ 1 (resulting in a smaller ρ 2 ) reduces the influence of neighboring nodes. In practice, the choice of ρ 1 and ρ 2 could increase the flexibility in designing the observers, and their specific selection should be tailored to the actual operating conditions.

4. Estimation Error Analysis

Denote e x i ( t ) = x i ( t ) x ^ i ( t ) , e ξ i ( t ) = ξ i ( t ) ξ ^ i ( t ) , and e f i ( t ) = f i ( t ) f ^ i ( t ) . Since e f i ( t ) = e ξ i ( t ) + ω E T e x i ( t ) , we can obtain the dynamics of the estimation error system as follows:
e ˙ x i ( t ) = ( A ρ 1 L C ) e x i ( t ) + E e ξ i ( t ) + ω E E T e x i ( t ) + Δ A x i ( t ) ρ 2 L 2 Δ 2 i
e ˙ ξ i ( t ) = ω E T E e ξ i ( t ) ω E T A e x i ( t ) ω 2 E E T E e x i ( t ) ω E T Δ A x i ( t ) + f ˙ i ( t )
According to Assumption 4, there exists a matrix B * satisfying ( I B B * ) E = 0 . Based on the fault estimation information, the control protocol proposed in this paper is presented as follows:
u i t = K x i ^ t B * E f ^ i t
where K is designed to guarantee that A B K is Hurwitz. Substituting Equation (9) into Equation (1), we can obtain
x ˙ i ( t ) = ( A B K ) x i ( t ) + Δ A x i ( t ) + B K e x i ( t ) + E e ξ i ( t ) + ω E E T e x i ( t )
From Equations (7), (8) and (10), the global error dynamic is
x ˙ ( t ) = ( I N ( A B K ) ) x ( t ) + ( I N Δ A ) x ( t ) + ( I N B K ) e x ( t ) + ( I N E ) e ξ ( t ) + ( I N ω E E T ) e x ( t )
e ˙ x ( t ) = ( I N ( A ρ 1 L 1 C ) ) e x ( t ) + ( I N E ) e ξ ( t ) + ( I N ω E E T ) e x ( t ) + ( I N Δ A ) x ( t ) ( L ρ 2 L 2 C ) e x ( t )
e ˙ ξ ( t ) = ( I N ω E T E ) e ξ ( t ) ( I N ω E T A ) e x ( t ) ( I N ω 2 E E T E ) e x ( t ) ( I N ω E T Δ A ) x ( t ) + f ˙ ( t )
where x ( t ) = [ x 1 ( t ) , x 2 ( t ) , , x N ( t ) ] T , e x ( t ) = [ e x 1 ( t ) , e x 2 ( t ) , , e x N ( t ) ] T , e ξ ( t ) = [ e ξ 1 ( t ) , e ξ 2 ( t ) , , e ξ N ( t ) ] T , f ( t ) = [ f 1 ( t ) , f 2 ( t ) , , f N ( t ) ] T .
Theorem 1.
Under Assumptions 1–3, the intermediate variable estimator in Equation (4) guarantees that the global error dynamic in Equations (11)–(13) is uniformly ultimately bounded for the given intermediate constant ω > 0 , ε i > 0 ( i = 1 , 2 , 3 , 4 ), and there exists matrix P j > 0 ( j = 1 , 2 , 3 ), Q, and a constant δ > 0 satisfying the following inequality:
Π i = Π ˜ 11 i Π ˜ 12 i * Π ˜ 22 i < 0 , i = 1 , 2 , , N
where Π ˜ 11 i = Π 11 i Π 12 i Π 13 i * Π 22 i Π 23 i * * Π 33 i , with Π 11 i = H e P 1 A P 1 B K + 1 ε 1 N T N + 1 ε 2 N T N + 1 ε 3 N T N , Π 12 i = P 1 B K + ω P 1 E E T , Π 13 i = P 1 E , Π 22 i = H e [ ( P 2 A ρ 1 Q C ) + ( ω P 2 E E T ) ] + δ ρ 2 λ ¯ i C T C , Π 23 i = P 2 E ω A T E P 3 ω 2 E E T E P 3 , Π 33 i = H e [ ω P 3 E T E ] , Π ˜ 22 i = diag 1 ε 1 , 1 ε 2 , 1 ε 3 , 1 ε 4 , Π ˜ 12 i = P 1 M 0 0 0 0 P 2 M 0 0 0 0 ω P 3 E T M P 3 . The observer gain matrices can be designed as L 1 = P 2 1 Q , L 2 = δ P 2 1 C T .
Proof of Theorem 1.
Choose a Lyapunov function as follows:
V ( t ) = x T ( t ) ( I N P 1 ) x ( t ) + e x T ( t ) ( I N P 2 ) e x ( t ) + e ξ T ( t ) ( I N P 3 ) e ξ ( t )
According to Equations (11)–(13) and Δ A ( t ) = M F ( t ) N , one has
V ˙ ( t ) = x T ( t ) H e I N ( P 1 A P 1 B K ) x ( t ) + 2 x T ( I N P 1 M F ( t ) N ) x ( t ) + 2 x T ( I N P 1 B K ) e x ( t ) + 2 x T ( t ) ( I N P 1 E ) e ξ ( t ) + 2 x T ( t ) ( I N ω P 1 E E T ) e x ( t ) 2 e ξ T ( t ) ( I N ω P 3 E T A ) e x ( t ) + e x T ( t ) H e I N ( P 2 A ρ 1 P 2 L 1 C ) e x ( t ) + e x T ( t ) H e L ( P 2 A ρ 2 P 2 L 2 C ) e x ( t ) + 2 e x T ( t ) ( I N P 2 E ) e ξ ( t ) + 2 e x T ( t ) ( I N ω P 2 E E T ) e x ( t ) 2 ω e ξ T ( t ) ( I N P 3 E T M F ( t ) N ) x ( t ) + 2 e x T ( t ) ( I N P 2 M F ( t ) N ) x ( t ) 2 e ξ T ( t ) ( I N ω P 3 E T E ) e ξ ( t ) 2 e ξ T ( t ) ( I N ω 2 P 3 E T E E T ) e x ( t ) + 2 e ξ T ( t ) ( I N P 3 ) f ˙ ( t )
From Lemma 1 and F ( t ) F T ( t ) I , the following inequalities hold for the positive constants ε 1 > 0 , ε 2 > 0 , ε 3 > 0 , and ε 4 > 0 :
2 x T ( t ) ( I N P 1 M F ( t ) N ) x ( t ) ε 1 x T ( I N P 1 M ) ( I N P 1 M ) T x ( t ) + 1 ε 1 x T ( I N N T N ) x ( t ) 2 e x T ( t ) ( I N P 2 M 1 F 1 ( t ) N 1 ) x ( t ) ε 2 e x T ( t ) ( I N P 2 M ) ( I N M T P 2 ) e x ( t ) + 1 ε 2 x T ( t ) ( I N N T N ) x ( t ) 2 e ξ T ( t ) ( I N ω P 3 E T M F ( t ) N ) x ( t ) ε 3 e ξ T ( t ) ( I N ω P 3 E T M ) ( I N ω M T E P 3 ) e ξ ( t ) + 1 ε 3 x T ( t ) ( I N N T N ) x ( t )
Based on Assumption 1, it can be inferred that there exists a scalar θ N > 0 such that the following inequality consistently holds:
2 e ξ T ( t ) ( I N P 3 ) f ˙ ( t ) 1 ε 4 e ξ T ( t ) ( I N P 3 ) ( I N P 3 ) e ξ ( t ) + ε 4 θ N
Furthermore, based on Equations (17) and (18), we have
V ˙ ( t ) x T ( t ) H e [ I N ( P 1 A P 1 B K ) ] x ( t ) + ε 1 x T ( t ) ( I N P 1 M M T P 1 ) x ( t ) + 1 ε 1 x T ( t ) ( I N N T N ) x ( t ) + 2 x T ( t ) ( I N ω P 1 E E T ) e x ( t ) + 2 x T ( I N P 1 B K ) e x ( t ) + 2 x T ( t ) ( I N P 1 E ) e ξ ( t ) + e x T ( t ) H e [ I N ( P 2 A ρ 1 Q C ) ] e x ( t ) + e x T ( t ) ( L + L T ) ( δ ρ 2 C T C ) e x ( t ) + 2 e x T ( t ) ( I N P 2 E ) e ξ ( t ) + 2 e x T ( t ) ( I N ω P 2 E E T ) e x ( t ) + ε 2 e x T ( t ) ( I N P 2 M M T P 2 ) e x ( t ) + 1 ε 2 x T ( t ) ( I N N T N ) x ( t ) 2 e ξ T ( t ) ( I N ω P 3 E T E ) e ξ ( t ) 2 e ξ T ( t ) ( I N ω P 3 E T A ) e x ( t ) 2 e ξ T ( t ) ( I N ω 2 P 3 E T E E T ) e x ( t ) + ε 3 e ξ T ( t ) ( I N ω 2 P 3 E T M M T E P 3 ) e ξ ( t ) + 1 ε 3 x T ( t ) ( I N N T N ) x ( t ) + 1 ε 4 e ξ T ( t ) ( I N P 3 P 3 ) e ξ ( t ) + ε 4 θ N
where Q = P 2 L 1 , L 2 = δ P 2 1 C T .
Obviously, Equation (19) is equivalent to
V ˙ ( t ) x ˜ T ( t ) Φ x ˜ ( t ) + ε 4 θ N
where
x ˜ ( t ) = x ( t ) e x ( t ) e ξ ( t ) , Φ = Φ 11 Φ 12 Φ 13 * Φ 22 Φ 23 * * Φ 33
and Φ 11 = H e [ I N ( P 1 A P 1 B K ) ] + ε 1 ( I N P 1 M M T P 1 ) + 1 ε 1 ( I N N T N ) + 1 ε 2 ( I N N T N ) + 1 ε 3 ( I N N T N ) , Φ 12 = ( I N P 1 B K ) + ( I N ω P 1 E E T ) , Φ 13 = ( I N P 1 E ) , Φ 22 = H e [ ( I N ( P 2 A ρ 1 Q C ) ) + ( I N ω P 2 E E T ) ] + ( L + L T ) ( δ ρ 2 C T C ) + ( I N P 2 M M T P 2 ) , Φ 23 = ( I N P 2 E ) ( I N ω A T E P 3 ) ( I N ω 2 E E T E P 3 ) , Φ 33 = H e [ ( I N ω P 3 E T E ) ] + ε 3 ( I N ω 2 P 3 E T M M T E P 3 ) + ε 4 ( I N P 3 P 3 ) .
From Equation (15), we can be obtain
V ( t ) max λ max P 1 , λ max P 2 , λ max P 3 ( x ( t ) 2 + e x ( t ) 2 + e ξ ( t ) 2 ) = max λ max P 1 , λ max P 2 , λ max P 3 e x ˜ ( t ) 2
It follows that
e x ˜ ( t ) 2 V ( t ) max λ max P 1 , λ max P 2 , λ max P 3
From Equation (20), we can obtain
V ˙ ( t ) λ max Φ e x ˜ ( t ) 2 + ε 4 θ N
It is obvious that λ max Φ < 0 , and then we have
V ˙ ( t ) κ V ( t ) + α
where κ = λ max Φ max λ max P 1 λ max P 2 , λ max P 3 < 0 , α = ε 4 θ N .
Denote a set Ω satisfying the following condition:
Ω = { x ( t ) , e x ( t ) , e ξ ( t ) λ min ( P 1 ) x ( t ) 2 + λ min P 2 e x ( t ) 2 + λ min P 3 e ξ ( t ) 2 α κ }
Let Ω ¯ be the supplementary set of Ω , and then the following inequality holds:
V ( t ) λ min P 1 ) x ( t ) 2 + λ min P 2 e x ( t ) 2 + λ min P 3 e ξ ( t ) 2 α κ
if x ( t ) , e x ( t ) , e ξ ( t ) Ω ¯ . From Equations (25) and (27), it is obvious that for x ( t ) , e x ( t ) , e ξ ( t ) Ω ¯ , we have
V ˙ ( t ) 0
which means that x ( t ) , e x ( t ) , e ξ ( t ) is uniformly ultimately bounded and converges to Ω exponentially with a rate greater than e α t from the Lyapunov stability theory.
Obviously, the inequality Φ < 0 is a preliminary condition to ensure the stability of the global error dynamic in Equations (11)–(13). However, it is not hard to find that Φ < 0 is still a high-dimensional and nonlinear characteristic. Therefore, in order to further ensure the solvability of the inequality Φ < 0 , this condition needs to be further transformed. Note that the term L + L T in Equation (21) is real symmetric, which means that it must have N real eigenvalues. By spectral decomposition of the real symmetric matrix L + L T , we obtain
L + L T = V ¯ Λ ¯ V ¯ T
where V ¯ is constructed from the eigenvectors of L + L T , and Λ ¯ = diag λ ¯ 1 , λ ¯ 2 , , λ ¯ N and λ ¯ i ( i = 1 , 2 , , N ) are the corresponding eigenvalues of L + L T . Certainly, the matrix V ¯ is orthogonal. Then, an orthogonal transformation matrix is defined as follows:
H = V ¯ T I n 0 0 0 V ¯ T I n 0 0 0 V ¯ T I r
By pre-multiplying and post-multiplying Φ < 0 with H and its transpose, we can obtain
Φ ˜ = Φ 11 Φ 12 Φ 13 * Φ ˜ 22 Φ 23 * * Φ 33 < 0
where Φ ˜ 22 = H e [ ( I N ( P 2 A ρ 1 Q 2 C ) ) + ( I N ω P 2 E E T ) ] + Λ ¯ ( δ ρ 2 C T C ) + I N P 2 M 1 M 1 T P 2 , and it can be found that the other terms are the same as in Equation (21). Finally, based on the Schur complement lemma, it can be found that Φ ˜ < 0 is equivalent to Equation (14) for i = 1 , 2 , , N . The proof of the theorem is complete. □
Remark 5.
In fact, excessively high dimensionality can negatively affect the accuracy of the LMI solution and may even result in no feasible solution for the LMI. Therefore, decoupling and dimension reduction for the preliminary condition Φ < 0 is very necessary.
Remark 6.
According to Equations (11)–(13), it can be found that there is a coupling between x ( t ) , e x ( t ) , and e ξ ( t ) . In order to prove the convergence of the system state and the observation error at the same time, it is necessary to choose a Lyapunov function containing x ( t ) , e x ( t ) , and e ξ ( t ) .
Remark 7.
In the proof of Theorem 1, by considering the term e x T ( t ) ( L + L T ) ( δ ρ 2 C T C ) e x ( t ) in Equation (20), it is found that L + L T is a real symmetric matrix, that is, the Laplacian matrix L of MASs can be asymmetric. In other words, the topology of MASs can be either undirected or directed in this paper. Therefore, the method proposed in this paper has better universality and practicability compared with most existing results.
Remark 8.
In this paper, the system state and the estimation error are analyzed strictly for convergence, and the explicit boundary, α κ , is obtained. Obviously, by adjusting the gain matrix K, parameters ω, and ε j ( j = 1 , 2 , 3 , 4 ) , parameter α is sufficiently large, and thus a relatively small bound is obtained. In addition, it is not difficult to find the convergence rate of the system state, and the global error dynamic in Equations (11)–(13) can be quantified by directly selecting the proper matrix K and parameter ω. Once the matrix K and parameter ω are determined, the observer gain matrices, L 1 and L 2 , can be obtained by solving the LMI in Equation (14). At the same time, the feasibility of the LMI in Equation (14) can be enhanced by adjusting parameters ε j ( j = 1 , 2 , 3 , 4 ) . On the other hand, if the LMI in Equation (14) has a feasible solution, the estimation and control performance can be enhanced by adjusting the matrix K and parameter ω. Therefore, this paper makes full use of the design freedom of parameters to ensure feasibility and effectiveness.

5. Numerical Simulation

In this section, a numerical simulation is used to demonstrate the effectiveness of the method proposed in this paper.
In this numerical simulation, a MAS with four agents is considered. The topology of the MAS is shown in Figure 1.
According to Figure 1, the adjacency matrix A and the Laplacian matrix L of the MAS are as follows.
A = 0 0 0 1 1 0 0 0 1 1 0 0 0 0 1 0 , L = 1 0 0 1 1 1 0 0 1 1 2 0 0 0 1 1
Obviously, since we consider a directed graph, the Laplacian matrix L is not symmetric.
In this paper, the actuator fault is considered. The system parameters of the ith ( i = 1 , 2 , 3 , 4 ) agent are given as
A = 1 2 3 2 2 1 3 1 1 , B = E = 0 0 2 , C = 1 0 0 0 1 0
M = 1 0 0 0 1 0 0 0 1 , N = 0.5 0 0 0 0.2 0 0 0 0.1 , F t = sin t
Obviously, it is not hard to find that r a n k ( C E ) r a n k ( E ) , that is, the fault estimation matching condition is not satisfied. In order to make the matrices ( A B K ) Hurwitz, the gain matrix K is chosen as K = 7.1153 4.6538 7.50 .
Without loss of generality, we make the following assumptions about the fault of each of the four agents, respectively:
f 1 t = 0 0 s t 10 s sin 0.5 t + cos 0.2 t 10 s < t 100 s
f 2 t = 0 0 s t 20 s 0.1 t 2 20 s < t 60 s 0.1 t 5 60 s < t 100 s
f 3 t = 0 0 s t 30 s 2 1 e t 30 30 s < t 60 s 1 2 1 e t 60 60 s < t 100 s
f 4 t = 0 0 s t 100 s
Based on the fault estimator given in the previous section, we choose the weight values as ρ 1 = 0.3 and ρ 2 = 0.7 , respectively. The intermediate constant ω in Equation (2) is selected as ω = 0.2 . Based on Theorem 1, the gain matrices of the intermediate variable estimator in Equation (4) can be expressed as follows:
L 1 = 17.5310 3.9315 32.1705 2.0079 69.3836 66.3622 , L 2 = 0.0261 0.3609 0.0088 0.4167 0.9831 2.3536
The simulation results are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. Figure 2, Figure 3, Figure 4 and Figure 5 show the effectiveness of fault estimation for the actuators of the four agents, respectively, where the red dashed line represents the result of the proposed method and the blue dashed line represents the result of the method proposed in [32]. As can be seen in Figure 2, Figure 3, Figure 4 and Figure 5, the fault estimation results of all agents are more satisfactory. Figure 6, Figure 7 and Figure 8 validate the efficacy of the fault-tolerant control protocol presented in this paper, and the results are better than those of the method proposed in [32]. Obviously, the state component of each agent converges to a small enough range under the designed fault-tolerant control protocol, and the influence of the actuator fault on the convergence process is greatly reduced. In addition, it is found that the fault estimation matching condition of fault estimation is not satisfied, and the bounds of the fault and its derivatives are not required. Therefore, the distributed fault estimation method based on the intermediate variable observer and the fault-tolerant control protocol based on fault estimation proposed in this paper are validated on a MAS with uncertainty.

6. Conclusions

In this paper, a novel robust distributed cooperative fault-tolerant control protocol is designed for a class of MASs with uncertainty and actuator faults. Unlike most existing approaches, in our method, the fault estimation matching condition is not necessary, and the bounds of the fault and its derivative are unknown. By introducing an intermediate variable and using both the centralized estimation errors and distributed estimation errors, the actuator fault is estimated, providing the basis for the fault-tolerance control scheme. Coordinate transformation and Schur decomposition are used to further reduce and decouple the LMI with high dimension and interference coupling, based on which the gain matrix of the estimator can be guaranteed. In future work, distributed cooperative fault-tolerant control for MASs with uncertainties in output and control input channels will be further considered.

Author Contributions

Methodology, J.S.; formal analysis, X.C.; investigation, S.X.; writing—original draft preparation, A.L.; project administration, C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by National Natural Science Foundation (NNSF) of China under Grants 62373184 and 62303217, and is also supported by Natural Science Foundation of Jiangsu Province under Grant BK20211502, in part by the Natural Science Foundation of the Jiangsu Higher Education Institutions of China under Grant 23KJB510006.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Dunbabin, M.; Marques, L. Robots for environmental monitoring: Significant advancements and applications. IEEE Robot. Autom. Mag. 2012, 19, 24–39. [Google Scholar] [CrossRef]
  2. Kamel, M.A.; Yu, X.; Zhang, Y. Formation control and coordination of multiple unmanned ground vehicles in normal and faulty situations: A review. Annu. Rev. Control 2020, 49, 128–144. [Google Scholar] [CrossRef]
  3. Chen, F.; Ren, W. On the control of multi-agent systems: A survey. Found. Trends Syst. Control 2019, 6, 339–499. [Google Scholar] [CrossRef]
  4. Loia, V.; Vaccaro, A. Decentralized economic dispatch in smart grids by self-organizing dynamic agents. IEEE Trans. Syst. Man Cybern. Syst. 2013, 44, 397–408. [Google Scholar] [CrossRef]
  5. Chang, J.; Guo, Z.; Cieslak, J.; Henry, D. Robust fault accommodation strategy of the reentry vehicle: A disturbance estimate-triggered approach. Nonlinear Dyn. 2021, 103, 2605–2625. [Google Scholar] [CrossRef]
  6. Qin, J.; Ma, Q.; Shi, Y.; Wang, L. Recent advances in consensus of multi-agent systems: A brief survey. IEEE Trans. Ind. Electron. 2016, 64, 4972–4983. [Google Scholar] [CrossRef]
  7. Qin, J.; Zhang, G.; Zheng, W.X.; Kang, Y. Adaptive sliding mode consensus tracking for second-order nonlinear multiagent systems with actuator faults. IEEE Trans. Cybern. 2018, 49, 1605–1615. [Google Scholar] [CrossRef] [PubMed]
  8. Yang, H.; Han, Q.L.; Ge, X.; Ding, L.; Xu, Y.; Jiang, B.; Zhou, D. Fault-tolerant cooperative control of multiagent systems: A survey of trends and methodologies. IEEE Trans. Ind. Inform. 2019, 16, 4–17. [Google Scholar] [CrossRef]
  9. Ding, S.X. Advanced Methods for Fault Diagnosis and Fault-Tolerant Control; Springer: Berlin/Heidelberg, Germany, 2021. [Google Scholar]
  10. Varga, A. Solving fault diagnosis problems. In Studies in Systems, Decision and Control, 1st ed.; Springer International Publishing: Berlin/Heidelberg, Germany, 2017; Volume 84, pp. 8–9. [Google Scholar]
  11. Gao, Z. Fault estimation and fault-tolerant control for discrete-time dynamic systems. IEEE Trans. Ind. Electron. 2015, 62, 3874–3884. [Google Scholar] [CrossRef]
  12. Liu, C.; Jiang, B.; Patton, R.J.; Zhang, K. Hierarchical-structure-based fault estimation and fault-tolerant control for multiagent systems. IEEE Trans. Control Netw. Syst. 2018, 6, 586–597. [Google Scholar] [CrossRef]
  13. Zhao, S.; Huang, B.; Zhao, C. Online probabilistic estimation of sensor faulty signal in industrial processes and its applications. IEEE Trans. Ind. Electron. 2020, 68, 8853–8862. [Google Scholar] [CrossRef]
  14. Shi, Y.; Hua, Y.; Yu, J.; Dong, X.; Ren, Z. Fully data-driven robust output formation tracking control for heterogeneous multiagent system with multiple leaders and actuator faults. IEEE Trans. Cybern. 2022, 54, 3183–3196. [Google Scholar] [CrossRef] [PubMed]
  15. Zhang, K.; Jiang, B.; Chen, M.; Yan, X.G. Distributed fault estimation and fault-tolerant control of interconnected systems. IEEE Trans. Cybern. 2019, 51, 1230–1240. [Google Scholar] [CrossRef] [PubMed]
  16. Chadli, M.; Davoodi, M.; Meskin, N. Distributed state estimation, fault detection and isolation filter design for heterogeneous multi-agent linear parameter-varying systems. IET Control Theory Appl. 2017, 11, 254–262. [Google Scholar] [CrossRef]
  17. Gao, M.; Yang, S.; Sheng, L. Distributed fault estimation for time-varying multi-agent systems with sensor faults and partially decoupled disturbances. IEEE Access 2019, 7, 147905–147913. [Google Scholar] [CrossRef]
  18. Menon, P.P.; Edwards, C. Robust fault estimation using relative information in linear multi-agent networks. IEEE Trans. Autom. Control 2013, 59, 477–482. [Google Scholar] [CrossRef]
  19. Han, J.; Liu, X.; Gao, X.; Wei, X. Intermediate observer-based robust distributed fault estimation for nonlinear multiagent systems with directed graphs. IEEE Trans. Ind. Inform. 2019, 16, 7426–7436. [Google Scholar] [CrossRef]
  20. Liu, X.; Gao, X.; Han, J. Distributed fault estimation for a class of nonlinear multiagent systems. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 3382–3390. [Google Scholar] [CrossRef]
  21. Liu, Q.; Zhang, K.; Jiang, B. Fixed-time fault estimation and prescribed performance fault-tolerant control for interconnected systems. IEEE Trans. Cybern. 2022, 54, 1084–1095. [Google Scholar] [CrossRef]
  22. Zhu, J.W.; Yang, G.H.; Wang, H.; Wang, F. Fault estimation for a class of nonlinear systems based on intermediate estimator. IEEE Trans. Autom. Control 2015, 61, 2518–2524. [Google Scholar] [CrossRef]
  23. Li, J.N.; Ren, W. Finite-horizon H fault-tolerant constrained consensus for multiagent systems with communication delays. IEEE Trans. Cybern. 2019, 51, 416–426. [Google Scholar] [CrossRef] [PubMed]
  24. Li, J.N.; Liu, X.; Ru, X.F.; Xu, X. Disturbance rejection adaptive fault-tolerant constrained consensus for multi-agent systems with failures. IEEE Trans. Circuits Syst. II Express Briefs 2020, 67, 3302–3306. [Google Scholar] [CrossRef]
  25. Zuo, Z.; Zhang, J.; Wang, Y. Adaptive fault-tolerant tracking control for linear and Lipschitz nonlinear multi-agent systems. IEEE Trans. Ind. Electron. 2014, 62, 3923–3931. [Google Scholar] [CrossRef]
  26. Deng, C.; Yang, G.H. Adaptive fault-tolerant control for a class of nonlinear multi-agent systems with actuator faults. J. Frankl. Inst. 2017, 354, 4784–4800. [Google Scholar] [CrossRef]
  27. Deng, C.; Yang, G.H. Distributed adaptive fault-tolerant control approach to cooperative output regulation for linear multi-agent systems. Automatica 2019, 103, 62–68. [Google Scholar] [CrossRef]
  28. Li, J.N.; Feng, H.; Wang, Y.; Liu, G.Y. A novel failure-distribution-dependent non-fragile H fault-tolerant load frequency control for faulty multi-area power systems. IEEE Trans. Power Syst. 2023, 39, 2936–2946. [Google Scholar] [CrossRef]
  29. Aryankia, K.; Selmic, R.R. Formation control and target tracking for a class of nonlinear multi-agent systems using neural networks. In Proceedings of the 2020 European Control Conference (ECC), St. Petersburg, Russia, 12–15 May 2020; pp. 160–165. [Google Scholar]
  30. Cui, Y.; Liu, X.; Deng, X.; Wang, Q. Observer-based adaptive fuzzy formation control of nonlinear multi-agent systems with nonstrict-feedback form. Int. J. Fuzzy Syst. 2021, 23, 680–691. [Google Scholar] [CrossRef]
  31. Liu, C.; Jiang, B.; Zhang, K. Distributed adaptive observers-based fault estimation for leader-following multi-agent linear uncertain systems with actuator faults. In Proceedings of the 2016 IEEE Chinese Guidance, Navigation and Control Conference (CGNCC), Nanjing, China, 12–14 August 2016; pp. 1162–1167. [Google Scholar]
  32. Liu, Y.; Yang, G.H. Integrated design of fault estimation and fault-tolerant control for linear multi-agent systems using relative outputs. Neurocomputing 2019, 329, 468–475. [Google Scholar] [CrossRef]
  33. Huang, S.J.; Yang, G.H. Fault tolerant controller design for T–S fuzzy systems with time-varying delay and actuator faults: A K-step fault-estimation approach. IEEE Trans. Fuzzy Syst. 2014, 22, 1526–1540. [Google Scholar] [CrossRef]
Figure 1. Topology of the MAS.
Figure 1. Topology of the MAS.
Sensors 24 02651 g001
Figure 2. f 1 ( t ) and its estimations in agent 1.
Figure 2. f 1 ( t ) and its estimations in agent 1.
Sensors 24 02651 g002
Figure 3. f 2 ( t ) and its estimations in agent 2.
Figure 3. f 2 ( t ) and its estimations in agent 2.
Sensors 24 02651 g003
Figure 4. f 3 ( t ) and its estimations in agent 3.
Figure 4. f 3 ( t ) and its estimations in agent 3.
Sensors 24 02651 g004
Figure 5. f 4 ( t ) and its estimations in agent 4.
Figure 5. f 4 ( t ) and its estimations in agent 4.
Sensors 24 02651 g005
Figure 6. The first component of x i ( t ) .
Figure 6. The first component of x i ( t ) .
Sensors 24 02651 g006
Figure 7. The second component of x i ( t ) .
Figure 7. The second component of x i ( t ) .
Sensors 24 02651 g007
Figure 8. The third component of x i ( t ) .
Figure 8. The third component of x i ( t ) .
Sensors 24 02651 g008
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

Shi, J.; Chen, X.; Xing, S.; Liu, A.; Chen, C. Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults. Sensors 2024, 24, 2651. https://doi.org/10.3390/s24082651

AMA Style

Shi J, Chen X, Xing S, Liu A, Chen C. Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults. Sensors. 2024; 24(8):2651. https://doi.org/10.3390/s24082651

Chicago/Turabian Style

Shi, Jiantao, Xiang Chen, Shuangqing Xing, Anning Liu, and Chuang Chen. 2024. "Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults" Sensors 24, no. 8: 2651. https://doi.org/10.3390/s24082651

APA Style

Shi, J., Chen, X., Xing, S., Liu, A., & Chen, C. (2024). Robust Cooperative Fault-Tolerant Control for Uncertain Multi-Agent Systems Subject to Actuator Faults. Sensors, 24(8), 2651. https://doi.org/10.3390/s24082651

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