Next Article in Journal
A Collaborative System Design for Avoiding and Removing the Unexpected Obstacles Encountered during Wheelchair Navigation
Next Article in Special Issue
Optimization of Mode-Switching Quality of Hybrid Tractor Equipped with HMCVT
Previous Article in Journal
Compilation Method of Flight Service Environment Spectrum Based on Altitude Features and Fuzzy Clustering
Previous Article in Special Issue
Optimization and Analysis of Clutch Switching Timing for Hybrid Tractors Equipped with Hydraulic Mechanical Combined Transmission
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compensation Observer-Based Adaptive Output Feedback Control for Multi-Agent Systems

School of Electronic and Electrical Engineering, Shanghai University of Engineering Science, Shanghai 261201, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(13), 5406; https://doi.org/10.3390/app14135406
Submission received: 2 March 2024 / Revised: 14 June 2024 / Accepted: 18 June 2024 / Published: 21 June 2024

Abstract

:
This study presents a modified compensating observer control strategy for nonlinear multi-agent systems affected by unknown hysteresis signal loops. Compared to conventional high-gain observers, this approach introduces a novel compensation signal, effectively reducing the tracking error of traditional observers. Then, by utilizing a backstepping method, an adaptive output feedback controller is designed, such that the tracking error converges to the small neighborhood around the origin. Simulation experiments with and without the compensation term demonstrate that this control strategy can effectively reduce error, but it increases input chattering to some extent.

1. Introduction

The coordinated control of multi-agent systems has emerged as a pivotal area of contemporary research, driven by the proliferation of distributed network resources and the imperative for collaboration among various mechanical devices [1,2,3]. Key aspects of such systems encompass group interaction and communication, coordination and cooperation, and conflict resolution. A prominent example of these concepts is the multi-flight simulator system [4]. Additionally, multi-agent systems have notable applications in the healthcare sector [5].
Exploration of multi-agent systems has yielded numerous methodologies and control strategies aimed at ensuring consistent stability. For instance, dynamic surface control combined with a first-order low-pass filter has been utilized to address the inherent “differential explosion” issue in traditional recursive methods [6,7,8,9]. While this method simplifies computational complexity, it also introduces errors associated with filter computation and increases the number of control parameters, potentially compromising transient performance management. Recent research [10] has attempted to resolve the algebraic loop problem intrinsic to the backstepping method, thereby ensuring transient tracking performance, albeit at the cost of increased control complexity due to the introduction of new variables.
In addressing the problem of unknown states in high-order multi-agent systems, observers are often employed. Among these, the use of high-gain observers is a common research approach [11,12]. Through the application of such observers, disturbances within the system, including external input disturbances, system coupling, and other factors, can be effectively managed [13,14,15,16,17]. These mechanisms ensure the stability of the entire closed-loop system while optimizing the tracking performance of each subsystem. Various types of observers employ different methodologies; for example, interval observers utilize neural networks, resilient control techniques, and fuzzy equations to estimate unmeasurable states of the system [18]. By narrowing the interval width, this method ensures that the closed-loop system remains semi-globally ultimately bounded, allowing for a convergence of tracking and observation errors within a small neighborhood near the origin [19]. Similar developments have been documented in discrete systems [20]. Optimizations of high-gain observers involve dimensionality reduction, as seen in the current literature [21,22] where the observer model’s complexity is optimized within manageable limits while ensuring that all closed-loop signals are ultimately bounded under the influence of a continuously differentiable switching function. Therefore, building upon the research foundation of high-gain observers, this study considers the addition of compensatory terms to achieve error reduction and structural simplification.
Furthermore, with the expansion of material science applications, particularly in multi-agent systems equipped with intelligent materials, significant challenges arise. For example, piezoelectric ceramics and ionic polymer metal composites have gained attention due to their exceptional performance and desirable physical characteristics [23]. However, these materials also exhibit performance deficiencies, such as non-smoothness, nonlinearity, and hysteresis, which can significantly impact the precision and stability of control systems, particularly due to the effects of hysteresis signals and coupling characteristics [24]. To address nonlinear hysteresis inputs, two primary strategies are typically employed. The first strategy involves mitigating hysteresis effects through the development of an adaptive hysteresis inverse controller [25]. The second strategy, similar to the methods outlined in [26], employs an adaptive algorithm to reduce the impacts of hysteresis by modeling it with both linear and nonlinear components. Commonly referenced models include the dead zone hysteresis model, the Prandtl–Ishlinskii (PI) hysteresis model, and the Bouc–Wen hysteresis model. This study discusses the processing of hysteresis signals using the Prandtl–Ishlinskii model. The specific work content is as follows:
(1)
This study introduces a modified compensating observer scheme. By incorporating this compensatory mechanism, the proposed method effectively reduces tracking errors and confines them to a small region near zero while ensuring stability, compared to traditional control schemes.
(2)
Unlike most output feedback systems, the system investigated in this study involves unknown parameters, nonlinear coupling, and hysteresis interference, and the proposed compensating observer has more practical significance.
The proposed scheme in this study demonstrates superior performance in handling nonlinear multi-agent systems through numerical simulations and case studies while maintaining low computational complexity. This provides new insights and methods for controlling nonlinear multi-agent systems and lays a foundation for future research directions. By comparing with existing mainstream control strategies, this study not only highlights the advantages of the new method but also reveals its potential limitations in practical applications, thereby providing valuable references for further optimization and improvement.
The remainder of this paper is organized as follows. Section 2 outlines the relevant background information, problem models, assumptions, and lemmas necessary for understanding the subsequent analysis. In Section 3, we propose a decentralized adaptive output feedback backstepping approach that incorporates a high-gain observer. Section 4 provides a comprehensive stability analysis, which substantiates the stability and tracking capabilities of the system. Finally, Section 5 corroborates the effectiveness of our proposed design through simulation results based on an example presented earlier in the paper.

2. Establishing the Problem Model

2.1. Preliminaries

According to the principles of the graph theory presented in [10,11,12], the relationship among followers is represented here by a directed topological graph G = { U , T , M } , where U = { 1 , , N }  denotes the set of nodes, T = U × U  defines the set of edges, and M = { a i j } R i  is the adjacency matrix. If node j  can send data to node i , then a i j = 1  ; otherwise, a i j = 0 . To facilitate system analysis, we employ the following lemmas and inequality theorems.
Lemma 1
[5]. The following formula is valid as long as there exist normal numbers Δ , σ 1  , and σ 2  and a symmetrical positive definite matrix P  , where E  is the identity matrix. D = diag [ 0 1 n i 1 ]  :
σ 1 E D P + P D + 2 Δ P σ 2 E .
Lemma 2
[17]. There is a positive definite matrix P  that meets the following requirements if matrix A i q i c i T  is a Hurwitz matrix:
A i q i c i T T P + P A i q i c i T = E .
Lemma 3
[9]. The following inequality is met for any d > 0  and x R i  (where R i  represents the set of real numbers):
0 x x 2 x 2 + d 2 < d .
Theorem 1.
Young’s inequality:
Θ i Γ e i j ρ 2 θ   Γ 2 e i j 2 + 1 2 ρ e i j ϵ 2 j ρ l i 1 + 2 Δ e i j 2 + 1 4 ρ l i ϵ 2 j 2
α i 1 y i ϵ 2 j e i j ρ l i 1 + 2 Δ α i 1 y i 2 e i j 2 + 1 4 ρ l i ϵ 2 j 2 α i 1 y i Θ i Γ e i j ρ 2 α i 1 y i 2 θ   Γ 2 e i j 2 + 1 2 ρ .
The parameters that are already in place are known positive design parameters. The other parameters will be determined later.

2.2. Basic Model and Problem Hypothesis

We consider a set of nonlinear systems defined as follows. First, Equation (6) delineates the state equation:
x ˙ i j = A i x i j + φ 0 , i , j y i + k = 1 r i φ k , i , j y i a k , i , j + F i , j y i x i j + b i H i , j ( y i ) w i y i = c i T x i .
Let x i j = x i 1 x i 2 x i n i T , x i j R i  denote the j -th state variable of the i -th subsystem, where j = 1 , 2 , , n i  and i = 1 , 2 , , n . The actual output of the i -th subsystem is y i R i . Here, the smooth nonlinear function vectors φ 0 , i , j ( y i ) R i  and φ k , i , j ( y i ) R i  are known, as is the lower triangular nonlinear function matrix F i , j ( y i ) . The matrix H i , j ( y i )  is a known function matrix, while a k , i , j  and b i  are unknown, non-zero constant parameters. The structure of the matrix is as follows:
A i = 0 E i 1 0 0
c i = 1 0 0 0 0 = 1 0 n i 1 × 1
F i , j y i = 0 0 0 0 f 2 , 2 0 0 f 3 , 2 f 3 , 3 0 0 0 0 f n , 2 f n , 3 f n , n i
H i , j = 0 z 2 , j 0 H z , j H n , j T
φ 0 , i , j y i = φ 0 , 1 , j y i φ 0 , 2 , j y i φ 0 , n , j y i T
φ k , i , j y i = φ k , 1 , j y i φ k , 2 , j y i φ k , n , j y i T
where w i  is an input for hysteresis. The Prandtl–Ishlinskii model was selected as the hysteresis model in this study. The relevant formula and mathematical model are as follows:
w i ( t ) = σ i u i ( t ) + F [ u i ] ( t ) .
In the expression F [ u i ] ( t ) = 0 R k ( r ) F r Ψ [ u i ] ( t ) dr , the function k ( r )  is a continuous, non-negative known function, while σ i  is an unknown, non-zero constant. This equation represents an area-based hysteresis model.
The design parameter R  serves as the upper limit of integration. In practical analysis, the influence of k ( r )  can be neglected when R  approaches infinity ( R = ). Consequently, both k ( r )  and F r Ψ [ u i ] ( t ) dr  are bounded, and k ( r )  is determined by experimental data. As a result, the functions F [ u i ] ( t )  are also bounded. Thus, the range of F [ u i ] ( t )  is constrained to F , or F r k [ u i ] ( t ) dr F .
This formulation describes the saturation play operator F r Ψ [ u i ] . The expression uses u i ( t )  and w i , respectively, as the input and output of symmetrical hysteresis, which is constrained by the curve Ψ  and represented by a well-known, continuous, non-decreasing function. Here, R  indicates the hysteresis threshold:
F r Ψ [ u i ] ( 0 ) = f r Ψ u i 0
F r Ψ [ u i ] ( t ) = f r Ψ u i t , F r Ψ [ u i ] ( t k )
f r Ψ u i t , w i = max Ψ u i r , min Ψ u i + r , w i
t k < t < t k + 1 , 0 k M 1 .
Next, we provide a monotone input signal, which is monotone in the interval [ t k , t k + 1 ] , and select the parameter number using the formula above, as follows:
σ i = 0.6 , Ψ u i = 6 tanh u i , R = 70 , k r = 0.9 e 0.04 r 0.5 2 u i = 4 sin 10 π t + cos 0.5 π t
In this way, the hysteresis model diagram shown in Figure 1 can be obtained.
The hysteresis effect set in this study has an upper and lower limit, as shown in Figure 1. Here, we use the parameter T  to represent the range of the hysteresis effect, which is T = H i , j ( y i ) F [ u i ] ( t ) . The parameter H  represents the last term of the hysteresis effect, and the parameter h  represents the minimum value of the hysteresis effect. Here, H = 20, h = −20.
Remark 1.
The hysteresis phenomenon represents one of the most prevalent nonlinear influences in real systems. However, the existing literature on hysteresis research, such as the research in [5,7,9,13,17,23], has not yet proposed efficient solutions. Therefore, this study continues to focus on hysteresis phenomena as the subject of analysis.
We suggest the following presumptions for the current system to simplify system analysis.
Proposition 1.
The output of the theory leader y r  and its derivatives ( y ˙ r  and y ¨ r  ) are known, bounded, and smooth functions. The simulation will provide this function’s information.
Proposition 2.
Let the parameter λ i  in the equation of state be a normal integer that can satisfy the following expression to maintain generality using the Hurwitz Equation (14):
χ i s = λ i , m s m + + λ i , 1 s 1 + λ i , 0 .
Remark 2
[27]. As the polynomial χ i s = λ i , m s m + + λ i , 1 s 1 + λ i , 0  is a Hurwitz polynomial, the positivity of coefficients ensures that the real parts of the roots (zeroes) of the system are negative, thereby guaranteeing system stability.
Proposition 3.
The leader transmits at least one directed spanning tree to the follower in the directed graph G  composed of system (7).
Proposition 4.
If a smooth and integrable bounded positive function σ i t  exists, the following expression (20) is fulfilled:
lim t + 0 t σ i ρ i d ρ i σ + .

3. Controller Design and Derivation

3.1. Design of the Gain Compensation Observer

We devised a number of compensation gain filters in the manner described below to estimate the unknown state of the nonlinear multi-agent system (7):
ε ˙ 0 i j = ( A i l i L i q i c i T ) ε 0 i j + F i , j y i ε 0 i j + φ 0 , i , j y i + l i L i q i y i ε 0 i j = ( A i l i L i q i c i T ) ε k i j + F i , j y i ε k i j + φ k , i , j y i ε ˙ s i j = ( A i l i L i q i c i T ) ε s i j + F i , j y i ε s i j + H i , j ( y i ) u i ( t ) ε ˙ z i j = ( A i l i L i q i c i T ) ε c i j + F i , j y i ε c i j + T H h 2 + σ 2
The factors in this case have values of i [ 1 n ] , q i = [ q 1 q n ] T R i , and q i , which correspond to the observer’s design parameters. The dynamic gain of the filter is l i , and its expression is l ˙ i = Λ i l i 2 + Λ i l i l 0 σ 2 + l i ζ ( y i ) , while the dynamic gain at position Λ i  denotes a favorable design feature. A function with slick positive elements is represented by ζ ( y i ) . A positive design parameter, represented by σ , will be explained subsequently. The gain begins with a number of l 0 .
We estimate the status as
x ^ i j = ε 0 i j + i = 1 r a k , i , j ε k i j + b i σ i ε s i j + b i ε z i j .
The definition error is ϵ ˙ i j = x i j x ^ i j , and its derivative is expressed as
ϵ ˙ i j = ( A i l i L i q i c i T ) ϵ i j + F i , j y i ϵ i j + b i T b i T H h 2 + σ 2 .
Transformation of this error is achieved as follows:
ϵ ¯ i j = l i Δ L i 1 ϵ i j .
We next define the Lyapunov equation as an expression of V ϵ ^ i j = ϵ ¯ i j T P i ϵ ¯ i j  and A i l i L i q i c i T = l i L i A i q i c i T L i 1  and combine Formulas (14)–(18):
V ˙ ϵ ¯ i j = l ˙ i l i 2 Δ P + D P + P D ϵ ¯ i j 2 l i ϵ ¯ i j 2 + ϵ ¯ i j T F T P + P F ϵ ¯ i j + 2 l i Δ L i 1 ϵ ¯ i j b i T T H h 2 + σ 2 .
Remark 3.
Compared to the methods described in [28,29], the quadratic form employed in this study not only exhibits fundamental quadratic characteristics but also offers a simpler form and more favorable mathematical properties. The convex nature of quadratic functions facilitates easier analysis and the derivation of system stability and convergence properties, as well as controller design.
Proposition 5.
There exists a smooth non-negative function ζ ( y i )  that fulfills F i y i ζ ( y i ) . The dynamic gain is subjected to the flattening function whose word is 1.6 y i 2 + 0.1 y i + 0.15  according to us.
The compensation component and hysteresis effect can be reduced as follows by combining Lemma 1, Lemma 2, Proposition 5, and the aforementioned theorems:
T T H h 2 + σ 2 = T h 2 + T σ 2 T H h 2 + σ 2 = T h 2 H + T σ 2 h 2 + σ 2 H + H σ 2 H h
l ˙ i l i 2 Δ P + D P + P D ϵ ¯ i j 2 σ 1 l ˙ i l i ϵ ¯ i j 2
F T P + P F 2 ζ ( y i ) P .
Remark 4.
The primary purpose of introducing T is to compensate for the self-interference term in the traditional high-gain observer, thereby reducing the impact of interference terms and consequently decreasing observational errors.
On the other hand, M = H + H σ 2 H h 2 l i Δ L i 1 ϵ ¯ i j , where M  represents a finite parameter. Finally, we can obtain the following:
V ˙ ϵ ¯ i j l ˙ i l i σ 1 + l i 2 P ζ ( y i ) ϵ ¯ i j 2 + M .
Remark 5.
The value of M is very small due to the influence of observer compensation, dynamic gains, and related parameters, which are consistently greater than 1. This fact corresponds with the stability proof content, ensuring the validity of the stability analysis.
Remark 6.
The selection criterion for σ  is to ensure that the value of this expression is minimized, ideally approaching the value of H . The purpose of this design is to ensure that the value of M  is sufficiently small, which will be further explained in the stability analysis section.

3.2. Controller Design

In this section, a controller is constructed by utilizing the adaptive backstepping method with a gain compensation observer. Subsystems are recursively designed.
Step 1: i = 1
The tracking error is defined as e 1 j = y 1 y r . The origin of the phrase is
e ˙ 1 j = x 2 j + φ 0 , 1 , j y i + k = 1 r i φ k , 1 , j y i a k , 1 , j + F 1 , j y i x 1 j y ˙ r .
The following changes were made:
Θ i = a 1 , 1 , j a 2 , 1 , j a r 1 , 1 , j a 1 , 2 , j a 2 , 2 , j a r 2 , 2 , j b i
Γ = φ 1 , 1 , j y i φ 2 , 1 , j y i φ r 1 , 1 , j y i ε 12 j ε 22 j ε r 2 2 j ε z 2 j T .
Substituting Equations (31) and (32) into (30) yields
e ˙ 1 j = ε 02 j + b i σ i ε s 2 j + ϵ 2 j + φ 0 , 1 , j y i + Θ i Γ + F 1 , j y i x 1 j y ˙ r
where γ θ  and γ p  are positive design parameters. Then, we synchronously set θ = max   Θ i 2 , σ i 2 . Here, θ  is an unknown parameter and ρ  is a positive design parameter.
Next, we apply Formula (34):
ε s i j = e i j + α i 1 .
At the same time, we set the auxiliary controller as
α 1 = p ^ α ¯ 1
α ¯ 1 = c 1 j e 1 j ε 02 j φ 0 , 1 , j y i + y ˙ r ρ 2 θ ^   Γ 2 e 1 j ρ l i 1 + 2 Δ e 1 j
where c 1 j  is the known normal number; θ ^  is the estimated value; and θ ˜ = θ θ ^ , θ ^  is the estimated value of the unknown parameter θ . Additionally, p ˜ = p p ^ , p ^  are estimates of unknown parameter p .
According to inequality Theorem 1, the two parameters are scaled as follows:
Θ i Γ e 1 j ρ 2 θ   Γ 2 e 1 j 2 + 1 2 ρ e 1 j ϵ 2 j ρ l i 1 + 2 Δ e 1 j 2 + 1 4 ρ l i ϵ ¯ 2 j 2 .
The Lyapunov equation can be changed into the following shape by defining the new parameter Ξ 1 . The expression of this function is defined below:
V ˙ 1 c 1 j e 1 j 2 + θ ^ e 1 j e 2 j + 1 2 ρ + 1 4 ρ l i ϵ ¯ 2 j 2 1 γ θ θ ˜ θ ^ ˙ Ξ 1 b i σ i γ p p ˜ p ^ ˙ + α ¯ 1 e 1 j + V ˙ ϵ ¯ i j
Step 2: i = 2
As above, we set the Lyapunov function under the condition of i = 2  and obtain the following expression according to Formula (34):
V ˙ 2 j = V ˙ 1 j + e 2 j e ˙ 2 j = V ˙ 1 j + e 2 j ε ˙ s 2 j α ˙ 1 V ˙ 1 j + e 2 j ε s 3 j l i L i q i ε s 1 j + F 2 , j y i ε s 2 j α ˙ 1 .
One can then derive the following equation from the auxiliary controller:
α ˙ 1 = α 1 θ ^ θ ^ ˙ + α 1 p ^ p ^ ˙ + α 1 y r y ˙ r + α 1 y ˙ r y ¨ r + α 1 l i l ˙ i + α 1 φ 0 , 1 , j φ ˙ 0 , 1 , j + α 1 ε 02 j ε ˙ 02 j + α 1 ε k 2 j ε ˙ k 2 j + α 1   Γ Γ + α 1 y φ 0 , 1 , j y i + ε 02 j + ϵ 2 j + Θ i Γ
In addition, we define the expression of the following parameters:
B 2 = l i L i q i ε s 1 j + F 2 , j y i ε s 2 j α 1 y r y ˙ r α 1 y ˙ r y ¨ r α 1 l i l ˙ i α 1 φ 0 , 1 , j φ ˙ 0 , 1 , j α 1 ε 02 j ε ˙ 02 j α 1   Γ ε ˙ k 2 j + ε ˙ z 2 j + φ ˙ k 2 j α 1 y φ 0 , 1 , j y i + ε 02 j α 1 p ^ p ^ ˙
Then, we can obtain the following expression:
V ˙ 2 j = V ˙ 1 j + e 2 j e 3 j + α 2 + B 2 α 1 θ ^ θ ^ ˙ α 1 y i ϵ 2 j + Θ i Γ .
At this time, we use inequality Theorem 1 again:
α 1 y i ϵ 2 j e 2 j ρ l i 1 + 2 Δ α 1 y i 2 e 2 j 2 + 1 4 ρ l i ϵ ¯ 2 j 2 α 1 y i Θ i Γ e 2 j ρ 2 α 1 y i 2 θ   Γ 2 e 2 j 2 + 1 2 ρ .
This theorem is then substituted into the Lyapunov equation, and the following parameters are defined:
  Ξ 1 = γ θ ρ 2   Γ 2 e 1 j 2   Ξ 2 = Ξ 1 + γ θ ρ 2 α 1 y i 2   Γ 2 e 2 j 2 + γ θ e 1 j e 2 j .
The expression of the auxiliary controller can be obtained with Formula (44):
α 2 = c 2 j e 2 j B 2 ρ l i 1 + 2 Δ α 1 y i 2 e 2 j θ ^ e 1 j ρ 2 α 1 y i 2 θ ^ Γ 2 e 2 j + α 1 θ ^ Ξ 2 γ θ γ θ ˜ θ ^
where γ θ ˜  is positive design parameter.
Substituting Formula (44) into (39), we can obtain the expression of the Lyapunov equation as
V ˙ 2 j c 2 j e 2 j 2 c 1 j e 1 j 2 + e 2 j e 3 j + 2 2 ρ + 2 4 ρ l i ϵ ¯ 2 j 2 1 γ θ θ ˜ θ ^ ˙ Ξ 2 b i σ i γ p p ˜ p ^ ˙ + α ¯ 1 e 1 j e 2 j α 1 θ ^ θ ^ ˙ Ξ 2 + γ θ γ θ ˜ θ ^ + V ˙ ϵ ¯ i j
Step 3: i [ 3 , z ) .
The following formula can be obtained using the same inversion calculation method:
α i = c i j e i j e i 1 j B i ρ l i 1 + 2 Δ α i 1 y i 2 e z j ρ 2 α i 1 y i 2 θ ^   Γ 2 e i j 2 + α i 1 θ ^ Ξ i γ θ γ θ ˜ θ ^ + γ θ q = 2 i 1 e q α q 1 y i ρ 2 α i 1 y i 2   Γ 2 e i j
Ξ i = Ξ i 1 + γ θ ρ 2 α 1 y i 2   Γ 2 e i j 2
V ˙ i j a = 1 i c a 2 e a 2 + e i 1 j e i j + i 2 ρ + i 4 ρ l i ϵ ¯ 2 j 2 1 γ θ θ ˜ θ ^ ˙ Ξ i b i σ i γ p p ˜ p ^ ˙ + α ¯ 1 e 1 j b = 2 i e b j α b 1 θ ^ θ ^ ˙ Ξ b + γ θ γ θ ˜ θ ^ + V ˙ ϵ ¯ i j
Step 4: i [ z , n ] .
This section adds the H i , j  function and hysteresis input environment. The input is defined as follows:
u i ( t ) = 1 H i , j ( y i ) α i ε s i + 1 j .
The expression of α i  is the same as that in Formula (46). Next, we set the Lyapunov function to (50)
V i j = 1 2 e i j 2 + V i 1 j .
The expression can be obtained by substituting the above parameters into the following calculation:
V ˙ n j i = 1 n c i 2 e i 2 + n 2 ρ + n 4 ρ l i ϵ ¯ 2 j 2 b i σ i γ p p ˜ p ^ ˙ + α ¯ 1 e 1 j 1 γ θ θ ˜ θ ^ ˙ Ξ n i = 2 n e i j α i θ ^ θ ^ ˙ Ξ n + γ θ γ θ ˜ θ ^ + V ˙ ϵ ¯ i j
Let the parameter value be
θ ^ ˙ = Ξ n γ θ γ θ ˜ θ ^ p ^ ˙ = α ¯ 1 e 1 j .
In combination with Equation (51), the final expression of the Lyapunov equation can be obtained as follows:
V ˙ n j i = 1 n c i 2 e i 2 + n 2 ρ + n 4 ρ l i ϵ ¯ 2 j 2 + γ θ ˜ θ ˜ θ ^ l ˙ i l i σ 1 + l i 2 P ζ ( y i ) ϵ ¯ i j 2 + M .

4. Stability Analysis

This section is dedicated to analyzing the stability and tracking performance of the proposed control scheme. The globally stable operation and precise tracking capabilities of the closed-loop, decentralized, adaptive control system are ensured by incorporating a specially designed high-gain compensation observer.
Remark 7.
Consider a closed-loop system composed of the controlled system (7), the hysteresis signal (11), the compensation filter (14), the update law (24), and Equations (31) and (33). Next, select appropriate values for the gain’s positive design parameters c i , γ θ ˜ , Λ i , l 0 , and ζ ( y i ) . The tracking error of each subsystem within the closed-loop decentralized control system converges to an arbitrarily small residual set, while all signals remain globally uniformly bounded.
Remark 8.
Notably, the proposed approach can effectively handle a wider range of interference signals, and system (7) encompasses more interference factors than systems in the existing literature. Thus, this system represents a broader category of multi-agent systems. The following are the key certification elements.
To summarize, let the Lyapunov function of the whole be V n j = V n 1 j + 1 2 e n j 2 , where the value is 2 , , n .
Next, use Young’s inequality, inequality Theorem 1, to extend and reduce Formula (55), as shown below:
γ θ ˜ θ ˜ θ ^ γ θ ˜ 2 θ 2 γ θ ˜ 2 θ ˜ 2 .
Because of the universal existence of n 4 ρ l i ϵ ¯ 2 j 2 n 4 ρ l i ϵ ¯ i j 2 , parameter ρ = n  and the following equations are designed:
L R = l ˙ i l i σ 1 + l i 1 4 l i 2 P ζ ( y i )
M ° = M + 1 2 + γ θ ˜ 2 θ 2 .
Combining (55)–(58), we can obtain Formula (59):
V ˙ n j i = 1 n c i 2 e i 2 γ θ ˜ 2 θ ˜ 2 L R ϵ ¯ i j 2 + M ° .
Then, set the parameters in dynamic gain l ˙ i = Λ i l i 2 + Λ i l i l 0 σ 2 + l i ζ ( y i )  as follows:
σ 1 = 2 P , σ 2 = σ 1 + 0.2 , Λ i = 1 2 σ 2 , l 0 = 1
where P  is set to 0.5 in the simulation. The following outcomes can be obtained using this setting technique, which effectively reduces the interference of parameter variables:
L R = 3 4 σ 1 2 σ 1 + 0.2 l i + P l 0 .
For Formula (59), we selected the design parameters to meet the standard requirements:
c i Ω γ θ ˜ 2 Ω γ θ L R Ω
where Ω  is the positive term parameter designed in this study. Combining this formula with Formula (59), we can obtain
V ˙ n j 2 Ω V n j + M ° .
Next, solve the integral inequality:
0 V n j t M ° 2 Ω + V n j 0 M ° 2 Ω e 2 Ω T .
In addition, when t , the limit value can be obtained:
lim t V n j t = M ° 2 Ω .
Here, the parameters ϵ i j , e i j , θ ˜ , and p ˜  are bounded. Through the design process of Formulas (22)–(34), we can obtain expressions of the following main variables:
ϵ ˙ i j = x i j x ^ i j e 1 j = y 1 y r e 2 j = ε s 2 j α 1 e i j = ε s i j α i 1 .
The adjusted observer parameter ϵ ¯ i j  is also bounded, as both the expected value y r  and dynamic gain l i  are bounded according to our theoretical framework. Once the parameters of the compensatory observer, such as ε 0 i j , ε k i j , ε s i j , and ε z i j , are confirmed to be finite, it can be demonstrated that the state variable x i j  and the auxiliary controller α i  are also bounded. As a result, the integral discrepancy shown in Equation (62) can be resolved as follows:
lim t e 1 j = lim t 2 V n j t = M ° Ω .
Ultimately, variables M ° and Ω  are unrelated to one another. Through properly increasing c i  , γ θ ˜  , γ θ  , L R  , and l 0  , the larger value of Ω  in Formula (59) can be selected. This step provides all remaining evidence.
Remark 9.
The compensation-based control strategy effectively reduces errors by lowering the key system parameters M ° . With appropriately chosen design parameters, the system errors can be reduced to zero under specific conditions.

5. Simulation Examples and a System Analysis of Control Strategies

This section presents two cases to illustrate the feasibility of our design strategy using numerical examples and gear system models. The objective of this control strategy is to progressively align each subsystem’s actual output y i  with its desired output y r  by utilizing an adaptive control law u i ( t ) , developed through a compensatory observer control scheme. In this context, the differential equation y ¨ r + 2 y ˙ r + y r = sin t  defines the optimal output y r .

5.1. Numerical Example

We consider the following second-order numerical system:
x ˙ 11 = x 12 x ˙ 12 = 0.2 sin 2 y 1 + f 1 , 2 x 2 + b 1 H 12 w 1 y 1 = x 11
where f 1 , 2 = 2 , b 1 = 1 , H 12 = 0.1 . The remaining parameters are all set to zero.
To validate the effectiveness of the control strategy, we set the initial state of the control system to zero and compare the simulation results between the compensation observer control scheme and the traditional observer control scheme.
The compensatory observer control scheme proposed in this design introduces a compensation term T H h 2 + σ 2  into the traditional observer. Consequently, the traditional observer control scheme follows the observer structure shown in Equation (67), but without the compensation term. To ensure a fair comparison, all other parameters remain identical:
ε ˙ 0 i j = ( A i l i L i q i c i T ) ε 0 i j + F i , j y i ε 0 i j + φ 0 , i , j y i + l i L i q i y i ε ˙ k i j = ( A i l i L i q i c i T ) ε k i j + F i , j y i ε k i j + φ k , i , j y i ε ˙ s i j = ( A i l i L i q i c i T ) ε s i j + F i , j y i ε s i j + H i , j ( y i ) u i ( t ) ε ˙ z i j = ( A i l i L i q i c i T ) ε c i j + F i , j y i ε c i j
Figure 2 presents the compensation observer control strategy, while Figure 3 illustrates the traditional controller control strategy. The error range in Figure 2b is smaller, indicating better tracking performance. This result further validates the superiority of the compensation observer control strategy.
Comparing the input waveforms between Figure 4 and Figure 5 indicates that under both control strategies, the input waveform of the overall system exhibits periodic variations. However, under the compensation strategy, the system input presents localized oscillations. This result represents a drawback of the control scheme, indicating a limitation in its ability to mitigate localized oscillations in the system input.
The dynamic behavior depicted in Figure 6, showing periodic states of the system output in three-dimensional state space, indicates the stable periodic responses of the system to certain inputs or initial conditions, indirectly suggesting a degree of stability in the system.
As shown in Figure 7a, the derivative of the dynamic gain exhibits periodic oscillations. This phenomenon is similar to observations in Figure 7b. Additionally, the waveforms of θ ^ ˙ i  and p ^ ˙  trend towards zero, indicating that the integrated values of these parameters are bounded. This characteristic may reflect the stability of the system dynamics, as the integrated values of the parameters are constrained, leading to limited variations in the system state within a finite range.
Based on the results shown in Figure 8, the waveforms of x i  and x ^ i  are nearly identical, indicating that the observer can accurately track and predict the system’s state variables. This result demonstrates that the observer effectively estimates the system state in dynamic environments, thereby providing a necessary foundation for the performance of the stringent controller.

5.2. Pratical Example

In this section, we employ a dual-motor coupled control system as an example of the simulation. The instance model of this control system is illustrated in Figure 9. Here, we treat the two motors as independent agents within a multi-agent system. Each agent is independently controlled, and overall stability is achieved by coordinating the two agents via controllers or information transmission.
The physical meanings of the parameters are as follows: θ 2  is the drive rotation angle, and θ 3  is the load rotation angle, which are represented by x 21  and x 31 , respectively. J 1  is the moment of inertia of the load, J d  is the moment of inertia of the driver, g r  is the gear tooth ratio, k 1  is the damping of the driver, k 2  is the load damping, and k  is the torsional elastic constant.
The state equation of the system is
x ˙ 21 = x 22 + φ 0 , 2 , 1 x ˙ 22 = φ 0 , 2 , 2 + f 2 , 2 x 2 + a k , 2 , 2 φ k , 2 , 2 + b 2 H 2 , 2 w 2 y 2 = x 22
x ˙ 31 = x 32 + φ 0 , 3 , 1 x ˙ 32 = φ 0 , 3 , 2 + f 3 , 2 x 2 + a k , 3 , 2 φ k , 3 , 2 + b 3 H 3 , 2 w 3 y 3 = x 32 .
In Formula (68) of subsystem 2,
f 2 , 2 = k 1 J d , b 2 = 1 , H 22 = 0.1 , a k , 2 , 2 = 1 φ 0 , 2 , 1 = 0.05 cos ( 3 y 2 ) 0.02 sin ( y 2 0.5 ) φ 0 , 2 , 2 = 0.5 sin 2 y 2 0.5 , φ k , 2 , 2 = k A J d A y 2 y 3 .
In Formula (69) in subsystem 3,
f 3 , 2 = k 1 J d , b 3 = 1 , H 32 = 0.1 , a k , 3 , 2 = 1 φ 0 , 3 , 1 = 0.02 cos ( y 3 ) 0.02 sin ( y 3 0.5 ) φ 0 , 3 , 2 = 0.5 tan 2 y 3 , φ k , 3 , 2 = k J 1 A y 2 B y 3 .
Next, we set the parameter values as follows:
A = 1 g r , k = 1 , k 1 = 0.8 , k 2 = 0.6 , J 1 = 1.2 , J d = 1 , g r = 0.8 , B = 1 , Δ = 1 , γ θ ˜ = 2 , γ θ = 4 c 1 = c 2 = 20 , q 1 = 2 , q 2 = 1 , ρ = 2 , .
To summarize, the simulation operation effect of the system can be obtained as follows:
Comparing Figure 10 with Figure 11 yields a consistent conclusion, indicating that employing the compensation observer control strategy can further reduce tracking errors. This finding underscores the efficacy of the compensation observer control strategy in enhancing the system’s tracking performance.
A comparative analysis of Figure 12 and Figure 13 similarly reveals an inherent flaw in the compensation observer control strategy: the exacerbation of input oscillations. This observation highlights a limitation of the compensational observer control strategy, indicating its tendency to amplify input oscillations. However, the degree of oscillation amplification is relatively mild and remains within an acceptable range.
Figure 14 graphically depicts the trajectories of Systems 2 and 3 within the state space, indicating a nearly complete overlap between them. This significant overlap not only underscores the high effectiveness of the control strategy in tracking the desired trajectory but also confirms the existence of a stable periodic response in both systems. The results from the comparative experiments further support these findings.
Figure 15 presents the actual and estimated values of the first state variable. Here, the overlap between x 21   and x ^ 21 , as well as that between x 31   and x ^ 31 , indicates that the observer can closely match the true state of the system. The superior fitting accuracy of the second state variable can be observed in Figure 16, as indicated by the minimal deviation between x 22   and x ^ 22 , as well as x 32   and x ^ 32 . This level of precision in estimating the second state variable highlights the enhanced performance of the observer in multi-dimensional settings, which is crucial to effectively control multi-agent systems.

6. Conclusions and Further Research

This study introduced a modified compensatory observation scheme that markedly differs from traditional observation methods and was specifically designed for a class of nonlinear multi-agent systems influenced by an unknown hysteresis signal loop. This innovative approach not only significantly reduces tracking errors but also substantially enhances overall system tracking performance. However, we observed that under certain conditions, the compensatory observer control strategy may inadvertently increase input oscillations. Although these oscillations are relatively mild, it is essential to consider their potential impacts on overall system control and stability.
Future research will focus on further optimizing the system’s design. This includes detailed studies on parameter adjustments, structural enhancements, sensor layout optimization, and other aspects. These efforts aim to refine the compensatory observation scheme, minimize any adverse effects, and improve the robustness and reliability of the control strategy. Ultimately, such advancements will contribute to significant technological progress and development in related fields.

Author Contributions

Algorithm derivation, Z.L.; simulation model, Y.L.; paper writing and editing, Z.L.; paper conception and ideas, Z.L. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China No. 61703269.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Acknowledgments

The authors would like to thank all of the authors cited in this article and the anonymous referees for their helpful suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Montagna, S.; Castro Silva, D.; Henriques Abreu, P.; Ito, M.; Schumacher, M.I.; Vargiu, E. Autonomous Agents and Multi-Agent Systems Applied in Healthcare. Artif. Intell. Med. 2019, 96, 142–144. [Google Scholar] [CrossRef] [PubMed]
  2. Cong, Z. Distributed ESO Based Cooperative Tracking Control for High-Order Nonlinear Multiagent Systems with Lumped Disturbance and Application in Multi Flight Simulators Systems. ISA Trans. 2018, 74, 217–228. [Google Scholar] [CrossRef] [PubMed]
  3. Wang, L.; Chen, M.Z.Q.; Gao, Z.; Wu, C.W.; Lin, Z. Special Issue on “Distributed Coordination Control for Multi-Agent Systems in Engineering Applications”. ISA Trans. 2017, 71, 1–2. [Google Scholar] [CrossRef] [PubMed]
  4. De Wilde, P.; Briscoe, G. Stability of Evolving Multiagent Systems. IEEE Trans. Syst. Man Cybern. Part B (Cybern.) 2011, 41, 1149–1157. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, Y.; Fang, Z.; Li, Y. Global Decentralized Output Feedback Control for Interconnected Time Delay Systems with Saturated PI Hysteresis. J. Frankl. Inst. 2020, 357, 5460–5484. [Google Scholar] [CrossRef]
  6. Nguyen, T.T.; Nguyen, N.D.; Nahavandi, S. Deep Reinforcement Learning for Multiagent Systems: A Review of Challenges, Solutions, and Applications. IEEE Trans. Cybern. 2020, 50, 3826–3839. [Google Scholar] [CrossRef] [PubMed]
  7. Cui, G.; Xu, S.; Ma, Q.; Li, Y.; Zhang, Z. Prescribed Performance Distributed Consensus Control for Nonlinear Multi-Agent Systems with Unknown Dead-Zone Input. Int. J. Control 2017, 91, 1053–1065. [Google Scholar] [CrossRef]
  8. Gong, A.; Xie, X.; Yue, D.; Xia, J. Multi-instant observer design of discrete-time fuzzy systems via an enhanced gain-scheduling mechanism. IEEE Trans. Cybern. 2022, 53, 2876–2885. [Google Scholar] [CrossRef] [PubMed]
  9. Hua, C.; Liu, S.; Li, Y.; Guan, X. Distributed adaptive output feedback leader-following consensus control for nonlinear multiagent systems. IEEE Trans. Syst. Man Cybern. Syst. 2018, 50, 4309–4317. [Google Scholar] [CrossRef]
  10. Huang, D.; Yang, C.; Pan, Y.; Cheng, L. Composite learning enhanced neural control for robot manipulator with output error constraints. IEEE Trans. Ind. Inform. 2019, 17, 209–218. [Google Scholar] [CrossRef]
  11. Li, W.; Zhang, Z.; Ge, S.S. Dynamic Gain Reduced-Order Observer-Based Global Adaptive Neural-Network Tracking Control for Nonlinear Time-Delay Systems. IEEE Trans. Cybern. 2022, 53, 7105–7114. [Google Scholar] [CrossRef] [PubMed]
  12. Li, Z.; Liu, J.; Huang, Z.; Peng, Y.; Pu, H.; Ding, L. Adaptive impedance control of human–robot cooperation using reinforcement learning. IEEE Trans. Ind. Electron. 2017, 64, 8013–8022. [Google Scholar] [CrossRef]
  13. Liu, G.; Zhao, L.; Yu, J. Adaptive finite-time consensus tracking for nonstrict feedback nonlinear multi-agent systems with unknown control directions. IEEE Access 2019, 7, 155262–155269. [Google Scholar] [CrossRef]
  14. Liu, Y.H. Adaptive dynamic surface asymptotic tracking for a class of uncertain nonlinear systems. Int. J. Robust Nonlinear Control 2018, 28, 1233–1245. [Google Scholar] [CrossRef]
  15. Liu, Y.J.; Tong, S. Barrier Lyapunov functions-based adaptive control for a class of nonlinear pure-feedback systems with full state constraints. Automatica 2016, 64, 70–75. [Google Scholar] [CrossRef]
  16. Liu, Y.J.; Tong, S. Barrier Lyapunov functions for Nussbaum gain adaptive control of full state constrained nonlinear systems. Automatica 2017, 76, 143–152. [Google Scholar] [CrossRef]
  17. Liu, Y.; Lin, Y. Global adaptive output feedback tracking for a class of non-linear systems with unknown backlash-like hysteresis. IET Control Theory Appl. 2014, 8, 927–936. [Google Scholar] [CrossRef]
  18. Liu, Y.; Lin, Y.; Huang, R. Decentralised adaptive output feedback control for interconnected nonlinear systems preceded by unknown hysteresis. Int. J. Control 2015, 88, 1712–1725. [Google Scholar] [CrossRef]
  19. Schwartz, M.; Krebs, S.; Hohmann, S. Guaranteed State Estimation Using a Bundle of Interval Observers with Adaptive Gains Applied to the Induction Machine. Sensors 2021, 21, 2584. [Google Scholar] [CrossRef] [PubMed]
  20. Shi, X.; Lu, J.; Li, Z.; Xu, S. Robust Adaptive Distributed Dynamic Surface Consensus Tracking Control for Nonlinear Multi-Agent Systems with Dynamic Uncertainties. J. Frankl. Inst. 2016, 353, 4785–4802. [Google Scholar] [CrossRef]
  21. Tao, G.; Kokotovic, P.V. Adaptive control of plants with unknown hystereses. IEEE Trans. Autom. Control. 1995, 40, 200–212. [Google Scholar] [CrossRef]
  22. Tong, S.; Li, Y.; Liu, Y. Observer-based adaptive neural networks control for large-scale interconnected systems with nonconstant control gains. IEEE Trans. Neural Netw. Learn. Syst. 2020, 32, 1575–1585. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, H.; Shen, H.; Xie, X.; Hayat, T.; Alsaadi, F.E. Robust Adaptive Neural Control for Pure-Feedback Stochastic Nonlinear Systems with Prandtl–Ishlinskii Hysteresis. Neurocomputing 2018, 314, 169–176. [Google Scholar] [CrossRef]
  24. Ye, W.; Aghili, F.; Wu, J.; Su, C.-Y. Robust Adaptive Control of Uncertain Systems Preceded with Unknown Hysteresis Actuators. In Proceedings of the 2021 IEEE 7th International Conference on Control Science and Systems Engineering (ICCSSE), Qingdao, China, 30 July–1 August 2021; pp. 23–28. [Google Scholar]
  25. Zhao, W.; Liu, Y.; Liu, L. Observer-based adaptive fuzzy tracking control using integral barrier Lyapunov functionals for a nonlinear system with full state constraints. IEEE/CAA J. Autom. Sin. 2021, 8, 617–627. [Google Scholar] [CrossRef]
  26. Zhou, M.; Zhang, Y.; Ji, K.; Zhu, D. Model Reference Adaptive Control Based on KP Model for Magnetically Controlled Shape Memory Alloy Actuators. J. Appl. Biomater. Funct. Mater. 2017, 15, 31–37. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, Q.; Liu, Y.; Zhang, G.; Li, H.; Liu, H. Novel Dynamic Surface Funnel Control for Input Delay Nonlinear Systems. Control. Theory Appl. 2022, 39, 1426–1432. [Google Scholar]
  28. Hosseinzadeh, M.; Yazdanpanah, M.J. Performance enhanced model reference adaptive control through switching non-quadratic Lyapunov functions. Syst. Control Lett. 2015, 76, 47–55. [Google Scholar] [CrossRef]
  29. Tao, G. Model reference adaptive control with L 1+ α tracking. Int. J. Control 1996, 64, 859–870. [Google Scholar] [CrossRef]
Figure 1. PI hysteresis input relation model diagram.
Figure 1. PI hysteresis input relation model diagram.
Applsci 14 05406 g001
Figure 2. Results under the compensation control strategy. (a) Trajectory tracking diagram of System 1. (b) System 1 output error.
Figure 2. Results under the compensation control strategy. (a) Trajectory tracking diagram of System 1. (b) System 1 output error.
Applsci 14 05406 g002
Figure 3. Results under the traditional control strategy. (a) Trajectory tracking diagram of System 1. (b) System 1 output error.
Figure 3. Results under the traditional control strategy. (a) Trajectory tracking diagram of System 1. (b) System 1 output error.
Applsci 14 05406 g003
Figure 4. Input to System 1 under the compensation control strategy.
Figure 4. Input to System 1 under the compensation control strategy.
Applsci 14 05406 g004
Figure 5. Input to System 1 under the traditional control strategy.
Figure 5. Input to System 1 under the traditional control strategy.
Applsci 14 05406 g005
Figure 6. State space tracking performance of System 1.
Figure 6. State space tracking performance of System 1.
Applsci 14 05406 g006
Figure 7. Results under the compensation control strategy. (a) Derivation of three parameters of System 1. (b) Waveform diagram of parameter l i . (c) Waveform diagram of parameter θ ^ ˙ i . (d) Waveform diagram of parameter p ^ ˙ .
Figure 7. Results under the compensation control strategy. (a) Derivation of three parameters of System 1. (b) Waveform diagram of parameter l i . (c) Waveform diagram of parameter θ ^ ˙ i . (d) Waveform diagram of parameter p ^ ˙ .
Applsci 14 05406 g007aApplsci 14 05406 g007b
Figure 8. Comparison of observer performance. (a) x 1 and x ^ 1  in System 1; (b) x 2  and x ^ 2  in System 1.
Figure 8. Comparison of observer performance. (a) x 1 and x ^ 1  in System 1; (b) x 2  and x ^ 2  in System 1.
Applsci 14 05406 g008
Figure 9. A dual-motor coupled control system.
Figure 9. A dual-motor coupled control system.
Applsci 14 05406 g009
Figure 10. Results under the compensation control strategy. (a) Trajectory tracking diagram of System 2 and 3. (b) System 2 and 3 output error.
Figure 10. Results under the compensation control strategy. (a) Trajectory tracking diagram of System 2 and 3. (b) System 2 and 3 output error.
Applsci 14 05406 g010
Figure 11. Results under the traditional control strategy. (a) Trajectory tracking diagram of System 2 and 3. (b) System 2 and 3 output error.
Figure 11. Results under the traditional control strategy. (a) Trajectory tracking diagram of System 2 and 3. (b) System 2 and 3 output error.
Applsci 14 05406 g011
Figure 12. Input to System 2 and 3 under the compensation control strategy.
Figure 12. Input to System 2 and 3 under the compensation control strategy.
Applsci 14 05406 g012
Figure 13. Input to System 2 and 3 under the traditional control strategy.
Figure 13. Input to System 2 and 3 under the traditional control strategy.
Applsci 14 05406 g013
Figure 14. State space tracking performance of system 2 and 3.
Figure 14. State space tracking performance of system 2 and 3.
Applsci 14 05406 g014
Figure 15. The comparison of observer performance.
Figure 15. The comparison of observer performance.
Applsci 14 05406 g015
Figure 16. The comparison of observer performance.
Figure 16. The comparison of observer performance.
Applsci 14 05406 g016
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

Liu, Z.; Liu, Y. Compensation Observer-Based Adaptive Output Feedback Control for Multi-Agent Systems. Appl. Sci. 2024, 14, 5406. https://doi.org/10.3390/app14135406

AMA Style

Liu Z, Liu Y. Compensation Observer-Based Adaptive Output Feedback Control for Multi-Agent Systems. Applied Sciences. 2024; 14(13):5406. https://doi.org/10.3390/app14135406

Chicago/Turabian Style

Liu, Zhaoyuan, and Ye Liu. 2024. "Compensation Observer-Based Adaptive Output Feedback Control for Multi-Agent Systems" Applied Sciences 14, no. 13: 5406. https://doi.org/10.3390/app14135406

APA Style

Liu, Z., & Liu, Y. (2024). Compensation Observer-Based Adaptive Output Feedback Control for Multi-Agent Systems. Applied Sciences, 14(13), 5406. https://doi.org/10.3390/app14135406

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