Adaptive Disturbance-Observer-Based Continuous Sliding Mode Control for Small Autonomous Underwater Vehicles in the Trans-Atlantic Geotraverse Hydrothermal Field with Trajectory Modeling Based on the Path

<jats:p>Considering intense hydrothermal activities and rugged topography in a near-bottom environment of the trans-Atlantic geotraverse (TAG) hydrothermal mound, a small autonomous underwater vehicle (S-AUV) will suffer from time-varying disturbances, model uncertainties, actuator faults, and input saturations. To handle these issues, a fault-tolerant adaptive robust sliding mode control method is presented in this paper. Firstly, unknown disturbances, model uncertainties, and actuator faults of the S-AUV are synthesized into a lumped uncertain vector. Without requiring the upper bound and gradient of the uncertainties, a continuous adaptive finite-time extended state observer is designed to estimate the lumped uncertain vector. Then, an auxiliary dynamic system composed of continuous functions is introduced to deal with input saturations, thereby contributing to achieving fixed-time trajectory tracking control of S-AUVs. Based on a designed continuous fixed-time nonsingular fast sliding mode surface, the proposed continuous adaptive controller is chattering free. Simulated topography is built according to topographic data of the TAG mound, and a smooth trajectory model is constructed by cubic spline interpolation. Comprehensive simulations performed on an actual S-AUV model are given to validate the effectiveness and superiority of the presented algorithm.</jats:p>

section_content: The deep-sea hydrothermal field is regarded as one of the most significant scientific discoveries in marine science in recent years [1] [2] [3]. To meet the urgent needs of hydrothermal research, the near-bottom exploration of hydrothermal fields is fast becoming a key instrument to conduct acoustic and optical data acquisition [4, 5]. Considering the rugged topography and concentrated regions of hydrothermal deposits in the near-bottom environment [6, 7], small autonomous underwater vehicles (S-AUVs) are applied to track the trajectory at 3 m height, identifying and observing hydrothermal deposits. Note that the battery capacity carried by S-AUVs is limited. Therefore, the combination of remotely operated vehicles (ROVs) and S-AUVs is chosen to save the battery energy of S-AUVs, that is, an ROV is applied to carry an S-AUV to sail from sea level to the vicinity of hydrothermal fields, then release the vehicle to complete exploration tasks. In active hydrothermal fields (such as the trans-Atlantic geotraverse (TAG) active mound), S-AUVs would inevitably suffer from unknown disturbances caused by the deep-sea environment and hydrothermal activities. Moreover, the highly non-linear strong coupling and uncertain characteristics of the S-AUV dynamics influence the controller design [8, 9]. Consequently, all these issues make the fast and accurate trajectory tracking control of S-AUVs in hydrothermal fields a challenging and complicated task. 

As one of the most efficient approaches to deal with non-linear systems under unknowns, perturbations, and external time-varying disturbances, the sliding mode control (SMC) method based on the sign function has been widely applied in trajectory tracking controls of multiple types of AUVs [9, 10]. It is worth noting the fact that finite-time (FT) control methods based on sliding mode technology have drawbacks [11]. Recently, fixedtime control methods have been developed in multi-agent systems, manipulators, and unmanned surface vehicles, which makes the settling time function completely independent of the initial states of control systems. There are few results in the tracking control system of AUVs. In Ref. [12], an adaptive fixed-time trajectory tracking control strategy of an S-AUV was applied, and it could realize a fixed-time convergence of tracking errors. For a special S-AUV, a novel adaptive continuous fixed-time sliding mode controller has been presented in Ref. [13]. An interval type-2 fuzzy neural network (NN) approximator has been utilized to approximate time-varying uncertainties. However, self-learning algorithms, including adaptive NN techniques [14] and fuzzy methods [15], give rise to some drawbacks in hardware [16]. Uncertainties and disturbances have not been adequately addressed in the above control algorithms. 

In order to further handle external disturbances and uncertainties, FT disturbance observers have been widely incorporated in various manners [17, 18]. Compared with other FT observers, the FT observer based on a sliding mode technology has great advantages of simple structure, fast convergence, and strong robustness. However, traditional sliding mode observers may cause the chattering phenomenon due to the existence of sign functions. In this context, the FT second-order sliding mode observer based on the extended state observer (ESO) structure becomes much more attractive, since flexible frameworks can not only be developed for eliminating the chattering effect and lumped uncertainties but also observe state variables quickly. In Ref. [17], a finite-time ESO (FTESO) has been innovatively designed to estimate lumped uncertainties. For the S-AUV under unknown parameters and external disturbances, an FT sliding mode disturbance observer has been designed in Ref. [18], whereby FT stable convergence of the whole FTESO-based control system can be guaranteed. 

On the other hand, actuator faults and input saturations are inevitable and widely exist in trajectory tracking systems of S-AUVs [19] [20] [21] [22]. In practice, the output power of thrusters is limited, especially underwater thrusters. Given the special working environment of underwater thrusters, whereby the density of seawater differs from that of air by about 900 times, underwater thrusters rarely exceed 3000 rpm. Therefore, the saturation of control input must be considered for the control system of S-AUVs [23] [24] [25]. In Ref. [23], a robust model predictive control with input and state constraints has been designed. In a constrained workspace, an online optimal control scheme has been applied. For tracking systems of S-AUVs with input saturations, an adaptive robust SMC method [24] and a robust tracking control method [25] have been presented, respectively. A novel auxiliary dynamic system (ADS) has been introduced into control laws to handle the impact of input saturations. As another issue that must be considered, the actuator fault is the most typical actuator problem appearing in S-AUV tracking systems [26, 27]. Therefore, it is optional and essential to solving the fault-tolerant control (FTC) problem of S-AUV systems under actuator faults. In the literature, faults are classified into lumped uncertainty terms for convenience of analysis and calculation. In Ref. [26], an adaptive FTC has been designed to achieve the fault reconstruction of S-AUVs, where FTESOs are applied to compensate for the effect of faults on control systems. With FT second-order sliding mode observers, an adaptive sliding mode FTC has been presented to achieve superior robustness under actuator faults [27]. 

In this paper, for the trajectory tracking control problem of an S-AUV in the TAG mound subject to unknown time-varying disturbances under uncertain model parameters, actuator faults, and input saturations, an adaptive fixed-time fault-tolerant method is constructed by combining the adaptive ESO, ADS, and nonsingular continuous fixedtime sliding mode (NCFTSM) surface. The three main contributions are described in the following parts:

(1) An adaptive observer-based trajectory tracking control method is proposed. It is proved that tracking errors of S-AUVs can converge to the residual set in a fixed time. (2) We design an adaptive FTESO to estimate the lumped uncertainties and realize the self-adaption of observer parameters. Compared with Ref. [28], the upper bound of lumped uncertainties and its gradient need not to be known in the observer design. (3) Simulated topography is built according to topographic data of the TAG hydrothermal mound. Based on the obtained path, a trajectory model is constructed by cubic spline interpolation. 

Our paper is organized in the following five main sections. The tracking problem and nonlinear model of an S-AUV are given in Section 2. Section 3 gives the design of fixed-time trajectory tracking control methods and the stability analysis. Three types of control methods are compared in Section 4. Section 5 draws the conclusions.
section_name: Introduction

section_name: Nonlinear Model and Problem Formulation

section_content: The motion of AUVs in three-dimensional space can be described in the body-fixed motion frame (G-XYZ) and earth-fixed frame (E-X 0 Y 0 Z 0 ), as shown in Figure 1. Similarly, the general kinematics and dynamics of the vehicle in three-dimensional space are modeled by these coordinate frames. The position and attitude of AUVs in the earth-fixed frame are denoted by η = [x, y, z, ϕ, θ, ψ] T, and υ = [u, v, w, p, q, r] T is used to represent the velocity of AUVs with respect to the body-fixed motion frame. Thus, the general kinematics of AUVs may be expressed as [29]. η = J(η)υ (1) where J(η) = diag(J 1 (η), J 2 (η)) ∈ R 6×6 is the spatial transformation matrix between the body-fixed motion frame and the earth-fixed frame.
section_name: Kinematics and Dynamics of the AUV

section_content: ). For the transformation matrix J(η), there exists an unknown positive constant J r > 0 such that J(η) ≤ J r.
section_name: Assumption 1 ([8]

section_content: ). The pitch cannot enter the neighborhood of θ = ±π/2 in the actual operation of an AUV due to the existence of static restoring forces. 

The dynamic equation of the AUV can be written as follows [29] :

where M ∈ R 6×6 represents the bounded matrix of the system inertia, including the added mass, C(υ) ∈ R 6×6 is the term of Coriolis-centripetal forces composed of Coriolis forces and centripetal forces, D(υ) ∈ R 6×6 is the hydrodynamic damping matrix, including damping forces and damping moments, g(η) ∈ R 6×1 denotes the vector of restoring forces and moments caused by gravity and buoyancy, τ dη = τ dη,1, τ dη,2, τ dη,3, τ dη,4, τ dη,5, τ dη,6

T denotes unknown external disturbances of AUVs, and

] T denotes the actual control input vector. 

With the actuator faults and saturations, in practice, the control input can be mathematically constructed as follows [30] :

where τ maxi (i = 1, 2,. . ., 6) represents the known maximum of the ith control input. In Equation ( 4 ), τ ci is the control action exerted on AUVs considering the actuator fault, τ ni is the torque command of the ith actuator control input from trajectory tracking controllers, τ i denotes the uncertain fault entering the AUV in an additive manner, e ii is an actuator fault effectiveness with 0 < e ii ≤ 1, and b i (tt 0i ) describes the time profile of a fault action on the ith actuator that occurs at time t 0i, which can be formulated as

where t 0i denotes the fault occurrence time, and a i > 0 denotes the evolution rate of unknown faults [30]. The dynamic equation of the AUV can be written as follows [29] :
section_name: Remark 1 ([8]

section_content: ,,,,,   τ denotes the actual control input vector. 

With the actuator faults and saturations, in practice, the control input can be mathematically constructed as follows [30] :

Earth-fixed frame Assumption 2. For the ∆τ = [∆τ 1,. . ., ∆τ 6 ] T, it always satisfies the following condition:

Assumption 3. The additive fault of the actuator τ i in ( 4 ) is bounded by a constant f 0 > 0

Taking the actuator faults and input saturations into account, the actual control input vector of AUVs can be described by

where

) is a function describing the time profile of a fault, and E = diag(e 11,. . ., e 66 ) is the effectiveness of the actuating component.
section_name: AUVs, and

section_content: According to the characteristics of the TAG active hydrothermal mound, the effects of the hydrothermal field on S-AUVs are discussed as follows. 

There are great topographic changes in the TAG active hydrothermal mound due to the intense hydrothermal activity and crustal movement. As shown in Figure 2, the active TAG mound consists of two circular terraces. There are escarpments of different heights with the slope of 25 In the active hydrothermal field shown in Figure 3, the plume is produced when hydrothermal fluids are expelled from the vent. Due to the great friction, the hydrothermal fluid continues to entrain the surrounding cold seawater and rise together [32, 33]. Correspondingly, the seawater near the vents will also be sucked into hydrothermal fluids below 0. 6 max Z ( max Z denotes the maximum rising height of hydrothermal plumes). 

When the S-AUV sails near the vent, it will be affected by the rising seawater being sucked into plumes (i. e., d f τ ). Meanwhile, the S-AUV also suffers from unknown disturbances from the deep-sea environment (i. e., dη τ ). In the active hydrothermal field shown in Figure 3, the plume is produced when hydrothermal fluids are expelled from the vent. Due to the great friction, the hydrothermal fluid continues to entrain the surrounding cold seawater and rise together [32, 33]. Correspondingly, the seawater near the vents will also be sucked into hydrothermal fluids below 0. 6 Z max (Z max denotes the maximum rising height of hydrothermal plumes). When the S-AUV sails near the vent, it will be affected by the rising seawater being sucked into plumes (i. e., τ d f ). Meanwhile, the S-AUV also suffers from unknown disturbances from the deep-sea environment (i. e., τ dη ). 

When the high-temperature hydrothermal activity occurs, the rising fluid continuously dissipates heat to the surrounding seawater and thereby leads to the change of the surrounding seawater in density. Distinctly, it has a great impact on the parameters g(η) of the S-AUV model in Equation (2). Considering that the S-AUV needs to be equipped with different equipment when carrying out different tasks, the variation of loads will also cause the perturbation of the dynamic model of S-AUVs. According to the above observations, the dynamic parameters M, C(υ), D(υ), and g(η) have great uncertainties in practical application. Thus, the parameters of the dynamic model in Equation ( 2 ) are regarded as the sum of the nominal dynamic terms and the uncertain terms

where M, Ĉ(υ), D(υ), and ĝ(η) denote the nominal dynamic terms, M, C(υ), D(υ), and g(η) are the uncertain terms. In the active hydrothermal field shown in Figure 3, the plume is produced when hydrothermal fluids are expelled from the vent. Due to the great friction, the hydrothermal fluid continues to entrain the surrounding cold seawater and rise together [32, 33]. Correspondingly, the seawater near the vents will also be sucked into hydrothermal fluids below 0. 6 max Z ( max Z denotes the maximum rising height of hydrothermal plumes). 

When the S-AUV sails near the vent, it will be affected by the rising seawater being sucked into plumes (i. e., d f τ ). Meanwhile, the S-AUV also suffers from unknown disturbances from the deep-sea environment (i. e., dη τ ). Accordingly, the S-AUV dynamic model in the active hydrothermal field can be rewritten as

where τ d f ∈ R 6×1 represents the effect from hydrothermal plumes and is defined as unknown disturbances. ∆N(η, υ) denotes the model uncertainties, which are defined as

Substituting Equation (7) into Equation ( 9 ), the complete dynamics of S-AUVs with actuator faults and saturations can be given by

where the lumped system uncertainties are defined as N(η, υ), which can be expressed as
section_name: Effect of the Hydrothermal Field on AUVs

section_content: Considering the S-AUV system described by Equations ( 1 ) and ( 11 ), the control objective is to design a fixed-time fault-tolerant adaptive robust sliding mode control (FTF-TARSMC) algorithm for trajectory tracking tasks in the TAG mound. The objective of this work is met under the conditions of unknown disturbances, model uncertainties, actuator faults, and input saturations.
section_name: Problem Objective

section_content: In this section, a new adaptive continuous disturbance observer based on the ESO technique is designed, firstly to estimate lumped system uncertainties in finite time. Then, the FTFTARSMC algorithm with FT disturbance observers is developed to achieve the control object, which can guarantee that the tracking errors of S-AUVs are bounded and can converge to a small residual set within a fixed time.
section_name: Control Design and Stability Analysis

section_content: The following definitions and lemmas for the fixed-time controller design and the stability analysis are presented here. Definition 1 ([35] ). Consider the nonlinear dynamical system x = f (x), x(0) = x 0 (13) where x ∈ R n is the system state, and f (x) : Ω → R n is continuous on an open neighborhood Ω a of the origin. The system ( 13 ) is said to be fixed-time stable if it is FT stable and the setting time T(x 0 ) is uniformly bounded for ∀x 0 ∈ R n, that is, ∃T m > 0 such that T(x 0 ) < T m. 

Lemma 1 ([35] ). For the nonlinear system (13), suppose that there exists a Lyapunov function V(x) : R n → R, such that the following inequality holds. 

V(x) ≤ -aV r (x) -bV l (x) + ς (14) where a > 0, b > 0, r > 1, 0 < l < 1, and ς are constants. Then, the system (13) is practical fixed-time stable. In addition, the residual set of the solution of system (13) can be given by

with 0 < υ < 1 and the settling time T(x 0 ) being bounded by T(x 0 ) ≤ 1 aυ(r-1) + 1 bυ(1-l). 

Lemma 2 ([36] ). For any δ ∈ R and h ∈ R, the following relationship holds

where m > 0, q > 1, and p > 1 are constants satisfying (q -1) (p -1) = 1.
section_name: Definitions and Lemmas

section_content: ). For any scalars ∏ i ∈ R, i = 1,. . ., n, 0 < r 1 < 1, and r 2 > 1, the following relationships hold
section_name: Lemma 3 ([37]

section_content: To eliminate the effect of uncertainties on the trajectory tracking system of S-AUVs, the adaptive FTESO is presented. In the analysis that follows, transformation variables of state space will be introduced to facilitate the design of disturbance observers, which can be defined as

Then, the mathematic model ( 1 ) and ( 11 ) can be converted into

. 

where

where

Assumption 4. The variable D including the lumped system uncertainties N(x 1, x 2 ) satisfies D ≤ H and.
section_name: Finite-Time Disturbance Observer Design

section_content: The first step of the disturbance observer design is to define an auxiliary variable as

where z 1 := [z 1, z 2,. . ., z 6 ] T and A = diag(A 1, A 2,. . ., A 6 ) are positive definite matrices. From Equation ( 21 ), the time derivative of z 1 is calculated as

where

x 1 + J M-1 ∆τ. Then, we introduce a new variable z 2 = D with considering the variable D as an extended state, the system (22) will be expanded as follows

where g(t) is the first-order derivative of the variable D. Finally, after defining the estimation error e 1 of auxiliary variable z 1 and the estimation error e 2 of disturbances z 2, the FT disturbance observer can be constructed as follows:

where ẑ1 and ẑ2 represent the estimates of z 1 and z 2, respectively. Each element of the vectors φ 1 = [φ 11, φ 12,. . ., φ 16 ] T and φ 2 = [φ 21, φ 22,. . ., φ 26 ] T is expressed as

where i = 1, 2,. . ., 6, 2/3 < β 1 < 1, and β 2 = 2β 1 -1. The adaptive gain parameters ρ j (t) (j = 1, 2,. . ., 5) are selected to satisfy

where the initial value R(0) > 0, ε 0 is a small positive constant, and k > 0 is the increased rate of adaptive gains. The gain C k (k = 1, 2, 3, 4) is the positive parameter and satisfies the following conditions:

Theorem 1. Under assumption 4, consider the expanded system (23). The FT disturbance observer described by Equations ( 25 )-( 27 ) can ensure that the estimation errors converge to a bounded neighborhood of the origin within a finite time T 0 under the condition (28). 

Proof of Theorem 1. Subtracting Equations ( 25 ) and ( 26 ) into ( 24 ) yields

. 

In order to facilitate the Lyapunov analysis, the following three variables are introduced:

It can be known from Equation (30) that if ξ l → 0 (l = 1, 2, 3) in a finite time, then e 1 → 0 and e 2 → 0 in a finite time. Substituting Equation ( 30 ) into (29), we have. 

Taking the derivative of Equation ( 30 ) and combining it with Equation (31), we can obtain

From Equation (32), we further have

where ξ = [ξ 1, ξ 2, ξ 3 ] T is a column vector. The matrix ∏ i (i = 1, 2, 3) is designed as follows:

Then, we need to prove the FT convergence of the matrix ξ. The following Lyapunov function for S-AUV tracking systems is constructed as

where P ∈ R 3×3 denotes a positive definite matrix and is chosen as follows

Substituting Equation ( 36 ) into ( 35 ), V can be written as

Differentiating V with respect to time and applying Equation ( 33 ) yields

From Equation ( 38 ), we further have

where

and A 3 = ξ T P∏ 3. In addition, matrices A 1 and A 2 are calculated as

Defining ∏ 31 =. 

R/2R ξ 1 2ξ 2 0 T and ∏ 32 = -0 0 g(t) i T, we have the following condition:

Combining with Equations ( 39 ) and ( 41 ), A 3 can be expressed as

with

From the above equations, A 31 and A 32 satisfy the following conditions:

, where λ 1

For the convenience of analysis and calculation,. 

V satisfies the following inequality:

. 

with

It can be easily verified that A 1 and A 2 are positive definite matrices satisfying the conditions (28). The following inequalities hold due to the properties of positive definite matrices

Similarly, for the matrix P defined in Equation ( 36 ), we can obtain

As a result, we have the following inequalities:

Substituting Equations ( 45 ), (46), and (51) into (47) yields

where the obtained known parameters υ 1 = λ max ∏ 31 /(λ min (P)),

and

)/λ min (P) are all constants. Since. R = k > 0, it is easy to find that Rυ 2 -. R/R υ 1 ≥ 0 will hold within a finite time T 0. After that, we have

According to the proof of lemma 4 in Ref. [38], the vector ξ will converge with the origin within a finite time T 0, and

/2 with a small positive constant θ ∈ (0, 1 ). Thus, we can obtain the result of Theorem 1. 

Remark 2. Based on the work in Ref. [39], although greater constraints are applied to the parameters, a more detailed stability analysis is performed, and the term (1/2)ρ 4i sig β 1 (e 1i ) is added to obtain a faster convergence rate of unknown terms. Compared with the observer in Ref. [37], the proposed observer needs not to know the upper bound of lumped uncertainties and its gradient. Therefore, it gives rise to superior performance in practical engineering.
section_name: D ≤ G. H and G exist but are not necessarily known.

section_content: In this subsection, the FTFTARSMC law is presented. Firstly, a fast NCFTSM manifold consisting of trajectory tracking errors is defined. Next, an ADS is introduced to deal with the input saturations. At last, we present the adaptive robust controller, which ensures that the tracking errors of S-AUVs are bounded and can converge to a small residual set Ω within the fixed time. 

Define the trajectory tracking errors as

where e 3 = [e 31, e 32, e 33, e 34, e 35, e 36 ] T, and η d ∈ R 6×1 denotes the desired trajectory. 

In order to obtain the fast fixed-time convergence performance of trajectory tracking errors e 3 and e 4, an NCFTSM manifold is constructed as

where β 3 = diag(β 31, β 32,. . ., β 36 ) and β 4 = diag(β 41, β 42,. . ., β 46 ) are all positive definite design matrices. The matrix F(e 3 ) ∈ R 6×1 is designed as

where s = e 4 + β 3 |e 3 | r 1 sign(e 3 ) + β 4 |e 3 | r 2 sign(e 3 ) + 2β 4 e 3. The designed constant terms should satisfy

Theorem 2. Consider the trajectory tracking error systems (55) and (56). For the sliding mode manifold satisfying s i = s i = 0, then errors e 3 and e 4 can converge with the residual set within a fixed time T 1, which can be further expressed as

Proof of Theorem 2. If s = s = 0 is reached, we can have the following condition:

The Lyapunov function is introduced as V 0 = 0. 5e T 3 e 3. Taking the time derivative of V 0, one can obtain. 

According to the definition |e 3 | r 1 sign(e 3 ) and the parameter r m (m = 1, 2), the following inequalities for. V 0 hold:. 

Based on Lemma 3, we obtain. 

To facilitate the proof, we define a new variable ζ = V 0 (1-r 1 )/2. Combining with Equation (62), the time derivative of ζ is calculated as

with the constant n 1 = (r 2 -1)/(1r 1 ). 

For |ζ| ≥ 1, by integrating and solving the differential Equation (63), one obtains

Further, we have

Subsequently, the convergence time for |ζ| ≥ 1 is calculated as

For |ζ| ≤ 1, by integrating and solving the differential Equation (63), we obtain

Next, the convergence time for |ζ| ≤ 1 is calculated as

Thus, the total convergence time T 1 is estimated as

The proof is finished.
section_name: Fixed-Time Tracking Control Design

section_content: Considering the NFSM surface in Ref [38], the terms 2β 4 e 3 and β 4 |e 3 | r 2 sign(e 3 ) are applied to improve the designed NCFTSM surface. After that, it guarantees that faster fixed-time convergence performance is realized in the region near and far away from the origin. 

Remark 4. In order to further improve the convergence rate and reliability,r 2 is replaced by the term κ 1i = 0. 5(r 2 + 1) + 0. 5(r 2 -1)sign(|e 3i | -1). Then, Equation ( 57 ) is equivalent to the following form:

Although the above-mentioned one has the advantage in convergence rate, its control law tends to be discontinuous. For that matter, it will give rise to the chattering problem. In the TAG mound, the accuracy of controllers will be inevitably affected if the control algorithm with chattering is applied to the tracking task of S-AUVs. 

In order to handle input saturations, an ADS is introduced as follows [40] :

where ξ = [ξ 1, ξ 2,. . ., ξ 6 ] T denotes the state of the auxiliary system, k 1 = diag(k 11, k 12,. . ., k 16 ) and k 2 = diag(k 21, k 22,. . ., k 26 ) are all positive definite matrices. The design parameters λ 1 > 1 and 0 < λ 2 < 1 are constants. The matrix h(ξ) is designed as

where ĥ(ξ) = 0. 5 -0. 5 sin(L(ξ)) and L(ξ) = (-π/2) (2 ξµ 1µ 2 )/(µ 2µ 1 ). The constants, µ 1, and µ 2 respectively, satisfy > 1/(λ 2 + 1) and 0 < µ 1 < µ 2. 

Based on the sliding mode manifold (57), ADS (70), and adaptive FTESO (25), we propose an FTFTARSMC algorithm for tracking errors (55) and (56) of S-AUVs as

where the matrix

) is positive definite, ε 2 denotes a small positive constant, and k A > 1. With â(0) > 0, â is the estimate of a and will be updated by the adaptive law (76). According to Lemma 1 of Ref [41], the gain â is bounded, i. e., there exists a positive constant a * satisfying a * ≥ â. 

In the following, the proof of the FTFTARSMC algorithm is given in Theorem 3.
section_name: Remark 3.

section_content: Considering the S-AUV model based on Equations ( 1 ) and ( 11 ) subject to unknown disturbances, model uncertainties, actuator faults, and input saturations under Assumptions 1-5, if the conditions λ min (k 4 ) > (λ max (k 5 ))

and λ min (k 3 ) > 6 (λ 1 -1)/2 (λ max (k 5 )) λ 1 +1 + 1 /(λ 1 + 1) are satisfied, the proposed controller (72)-( 75 ) with the adaptive law (76), ADS (70), adaptive FTESO (25), and NCFTSM manifold (57) can guarantee that the tracking errors (55) and (56) of S-AUVs are bounded and can converge to a residual set Ω E within a fixed time T 2. 

Proof of Theorem 3. Define the error a = a *â and consider the following Lyapunov function

1+λ 1 2

According to the characteristics of the ADS, Equations ( 70 ) and ( 71 ), (78) can be arranged as

where

The value of the matrix χ is not required to be known; it is only used for theoretical analysis. 

According to Lemma 1, it is known that the S-AUV trajectory tracking control system in hydrothermal fields can realize practical fixed-time stability. From Equation ( 77 ),

T settles within the residual set

with 0 < δ < 1, the fixed time T 2 is expressed by T 2 ≤ 2/(aδ(λ 1 -1)) + 2/(bδ(1λ 2 )). 

Further, according to Equation (79), it is worth noting that s, ξ, and a in Equation (77) will converge to Ω E and be stabilized in Ω E if the initial variables are designed as

with Ω E 1 being a compact set and a(0) ∈ Ω a. Thus, it is concluded that the S-AUV trajectory tracking control system in hydrothermal fields is semi-globally uniformly ultimately bounded. For the detailed derivation processes of (77) and (78), please refer to Appendix A. 

Remark 5. Compared with Refs [42, 43], the main improvements of the proposed trajectory tracking algorithm include: (1) Based on the analysis of the hydrothermal environment, the S-AUV control algorithm is modified with considering unknown disturbances, parameter perturbations, input saturations, and actuator faults; (2) A multi-parameter adaptive FTESO is constructed without knowing the upper bound of the lumped uncertainties and its gradient. 

Remark 6. Due to the finite-time convergent observer errors, -ẑ 2 + J M-1 N(x 1, x 2 ) may not converge to zero at T 2, and s T J M-1 N(x 1, x 2 ) -ẑ2 ≤ κ holds with κ ≥ 0 being a positive constant. The constant κ only affects the parameter χ in Equation (79). When χ is large, the tracking errors may increase and lead to loss of control accuracy. The complex TAG mound and the poor performance of the forward-looking obstacle avoidance sonar at 3 m height make the fixed-time convergence ability of S-AUV tracking systems a top priority. Actually, the control performance is accepted based on the simulation results.
section_name: Theorem 3.

section_content: In this section, we will build the TAG active hydrothermal mound to construct the actual rugged topography in the simulink. According to the chosen path points, the desired trajectory is modeled by cubic spline interpolation. 

The topography in simulations is built based on topographic data of the TAG hydrothermal mound [44, 45], as shown in Figure 4. According to task requirements, the S-AUV starts from the release point, sails along the trajectory at 3 m height, passes through dense areas of hydrothermal deposits, and takes circular observations for the large active hydrothermal vents at the northwest of mound, as shown in Figure 5. 

avoidance sonar at 3 m height make the fixed-time convergence ability of S-AUV tracking systems a top priority. Actually, the control performance is accepted based on the simulation results.
section_name: Topography Building and Trajectory Modeling

section_content: In this section, we will build the TAG active hydrothermal mound to construct the actual rugged topography in the simulink. According to the chosen path points, the desired trajectory is modeled by cubic spline interpolation. 

The topography in simulations is built based on topographic data of the TAG hydrothermal mound [44, 45], as shown in Figure 4. According to task requirements, the S-AUV starts from the release point, sails along the trajectory at 3 m height, passes through dense areas of hydrothermal deposits, and takes circular observations for the large active hydrothermal vents at the northwest of mound, as shown in Figure 5. avoidance sonar at 3 m height make the fixed-time convergence ability of S-AUV tracking systems a top priority. Actually, the control performance is accepted based on the simulation results.
section_name: Topography Building and Trajectory Modeling

section_content: In this section, we will build the TAG active hydrothermal mound to construct the actual rugged topography in the simulink. According to the chosen path points, the desired trajectory is modeled by cubic spline interpolation. 

The topography in simulations is built based on topographic data of the TAG hydrothermal mound [44, 45], as shown in Figure 4. According to task requirements, the S-AUV starts from the release point, sails along the trajectory at 3 m height, passes through dense areas of hydrothermal deposits, and takes circular observations for the large active hydrothermal vents at the northwest of mound, as shown in Figure 5. To construct a three-dimensional trajectory from the starting point P 1 to the final point P f in the time interval of T 1, T f, points on the path are collected according to the path characteristics at first. Then, the whole trajectory is divided into an approaching part and an encircling part. In addition, appropriate boundary conditions are derived to guarantee the continuity of velocities. Next, cubic spline interpolation based on collected points is applied to construct the approaching trajectory to ensure that the interpolated trajectory is as close to the corresponding path as possible. 

The approaching trajectory describes the spatial curve of S-AUVs from the starting point P 1 to the point near the vent P 38 in [T 1, T 38 ]. To create the approaching trajectory, 38 accuracy points on this path are considered as

where T j is the time instant at which the S-AUV plans to point P j. Based on cubic spline interpolation methods, the mathematical models X j, Y j, and Z j with respect to time are constructed, and boundary conditions are guaranteed at the point (X 38, Y 38, Z 38 ). Thus, the approaching trajectory can be defined as follows:

x 1d = spline(X, T, t 1 )

where

, and T a = [T 1. . . T 38 ] ∈ R 1×38 are constant vectors. t 1 denotes the total time of the approaching trajectory. According to cubic spline interpolation, x 1d, y 1d, and z 1d in the approaching trajectory will be composed of 37 cubic polynomials for 38 points in [T 1, T 38 ] and can be written as follows:

x 1d (i, 1) = x 1d. coefs(i, 1) (tx 1d. breaks(i)) 3 +x 1d. coefs(i, 2) (tx 1d. breaks(i)) 2 +x 1d. coefs(i, 3) (tx 1d. breaks(i)) + x 1d. coefs(i, 4) y 1d (i, 1) = y 1d. coefs(i, 1) (x 1d (i, 1)y 1d. breaks(i)) 3 +y 1d. coefs(i, 2) (x 1d (i, 1)y 1d. breaks(i)) 2 i f T i ≤ t ≤ T i+1 +y 1d. coefs(i, 3) (x 1d (i, 1)y 1d. breaks(i)) + y 1d. coefs(i, 4) z 1d (i, 1) = z 1d. coefs(i, 1) (x 1d (i, 1)z 1d. breaks(i)) 3 +z 1d. coefs(i, 2) (x 1d (i, 1)z 1d. breaks(i)) 2 +z 1d. coefs(i, 3) (x 1d (i, 1)z 1d. breaks(i)) + z 1d. coefs(i, 4)

where i = 1, 2,. . ., 37, x 1d. coefs ∈ R 37×4, y 1d. coefs ∈ R 37×4, and z 1d. coefs ∈ R 37×4 are the polynomial coefficients obtained by applying the cubic spline interpolation. 

x 1d. breaks ∈ R 1×38, y 1d. breaks ∈ R 1×38, and z 1d. breaks ∈ R 1×38 represent the known interpolated points along the curve. For the attitude in the approaching trajectory, we choose roll 1d = 0rad and pitch 1d = 0rad based on task requirements. Similar to x 1d, yaw 1d can be designed as follows:

where Yaw = [Yaw 1. . . Yaw 38 ] ∈ R 1×38, yaw 1d. coefs ∈ R 37×4 is the polynomial coefficients, and yaw 1d. breaks ∈ R 1×38 is the interpolated points. 

Remark 7. For cubic spline interpolation, it is common to apply as many points as possible to achieve superior performance in terms of interpolation and community. However, the disadvantages are as follows: (1) More points means more precise division and interpolation for the time interval, and more workload is required. (2) The generation of more cubic polynomials gives rise to the reduction in the original interval width, which may cause an increase and a large change in the control input. 

The encircling trajectory denotes the observation trajectory of S-AUVs for the hydrothermal vents and the surrounding rocks. It is required that the high-definition camera at the head of S-AUVs always aims at hydrothermal deposits, and the 360 • video acquisition of surrounding rocks and vents needs to be guaranteed. According to the location of the hydrothermal vent x c = 210, 547. 24 mm and y c = 19, 000. 57 mm, the circumvent trajectory is designed as follows:

x od = 10 sin(πt/20) + 210. 54724 m y od = -10 cos(πt/20) + 19. 00057 m z od = 2t/20 + 6. 03737 m roll od = 0 rad pitch od = 0 rad

Since the heading direction of S-AUVs always points to hydrothermal deposits, the yaw should be classified based on the desired position of S-AUVs as follows:

where L 2 = (x odx c ) 2 + (y ody c ) 2. Combined with Equations ( 83 ) and ( 85 )-(87), the whole trajectory of S-AUVs can be constructed. For evaluating the interpolation accuracy of the created approaching trajectory, 335 path points (including the 38 points used) are selected to construct the original path. Then, the difference between the trajectory model and the original path is represented by horizontal and vertical planes. Figure 6 shows the slight difference between the trajectory model and the original path in the horizontal plane. In the vertical plane, the trajectory in the rising region (126, 412 mm ≤ x ≤ 145, 153 mm) is largely different from the original path. In the area of vertical descent and horizontal motion, the original path exhibits vertical change at x = 162747mm. Thus, a large difference is made. where ( ) ( ) In order to further demonstrate a clear comparison, the single error ( )

and the accumulated error errorall between the approaching trajec- tory and the original path are calculated by the quantitative analysis: In order to further demonstrate a clear comparison, the single error error(i) (i = 1, 2,. . ., 335) and the accumulated error errorall between the approaching trajectory and the original path are calculated by the quantitative analysis:

where x 1di, y 1di, and z 1di denote the position value of the ith point calculated in the approaching trajectory model. After that, the total accumulated error errorall is acceptable, since the total accumulated error is less than 30 m. However, it is necessary to ensure that the total accumulated error is not excessively concentrated at some points. 

According to the results
section_name: Topography Building and Trajectory Modeling

section_content: 5 : 38 errorall and ( )
section_name: ( )

section_content: , the accumulated error of the trajectory in the rising area, as well as the area of vertical descent and horizontal motion, accounts for about two-thirds of the total accumulated error. It can be seen from Figures 6 and 7 that the accumulated error mainly comes from the vertical position, i. e.,

In the rising area, due to the small change in hydrothermal terrain, the accumulated error with an average single error of less than 0. 7 m has little effect on the tracking motion of S-AUVs at 3 m height. As shown in Figure 7b, the intersection angle of the vertical descent curve and the horizontal motioning curve at 162747 mm x = is equal to 90 °. However, it is impossible for a smooth trajectory. The accuracy of the tra- jectory needs to be reduced. Nevertheless, it is acceptable for the difference in the area of vertical descent and horizontal motion, since the open terrain is nearby, the curve formed by the trajectory is far away from the fault, and the accumulated error mainly concentrates A three-dimensional approach trajectory with the starting point P 1 (125. 623, 0, 35. 729) m to the point P 335 (210. 54724, 9. 00057, 60. 3737) m is acceptable, since the total accumulated error is less than 30 m. However, it is necessary to ensure that the total accumulated error is not excessively concentrated at some points. According to the results errorall(5 : 38) and errorall(84 : 97), the accumulated error of the trajectory in the rising area, as well as the area of vertical descent and horizontal motion, accounts for about two-thirds of the total accumulated error. It can be seen from Figures 6 and 7 that the accumulated error mainly comes from the vertical position, i. e., ∑|Z iz 1di |. In the rising area, due to the small change in hydrothermal terrain, the accumulated error with an average single error of less than 0. 7 m has little effect on the tracking motion of S-AUVs at 3 m height. As shown in Figure 7b, the intersection angle of the vertical descent curve and the horizontal motioning curve at x = 162,747 mm. is equal to 90

•. However, it is impossible for a smooth trajectory. 

The accuracy of the trajectory needs to be reduced. Nevertheless, it is acceptable for the difference in the area of vertical descent and horizontal motion, since the open terrain is nearby, the curve formed by the trajectory is far away from the fault, and the accumulated error mainly concentrates on the vertical position in this area. Based on the above analysis, the whole desired trajectory in the constructed TAG mound can be obtained by applying models (83) and ( 85 )-(87), as shown in Figure 8. 

Remark 8. Different from most AUV trajectory tracking control works, the desired trajectory in our paper is constructed by interpolating the task path. The S-AUV can complete the near-bottom exploration task of TAG active mounds by tracking this trajectory. This means that engineering practicability is obtained.
section_name: : 97 errorall

section_content: In this section, the effect of hydrothermal fluids on S-AUVs d f τ is constructed. The tracking performance of FTFTARSMC algorithms in the TAG mound is verified by numerical simulations. 

As shown in Figure 9, the torpedo-type S-AUV has a length of more than 1700 mm, a diameter of 200 mm, and a weight of 70 kg in air. The propulsion system of the S-AUV consists of eight independent thrusters (including two vertical channel thrusters, two vertical conduit thrusters, two lateral channel thrusters, and two main thrusters) to achieve a fully actuated system given in Figure 10. Then, according to the thrust curves, the saturation limits for the main thrusters, channel thrusters, and vertical conduit thruster are roughly taken as ±600 N, ±500 N, and ±150 N, respectively. In the simulation, the thruster outputs from eight thrusters can be obtained in the form of total control forces and moments v τ at each degree of freedom. Remark 8. Different from most AUV trajectory tracking control works, the desired trajectory in our paper is constructed by interpolating the task path. The S-AUV can complete the near-bottom exploration task of TAG active mounds by tracking this trajectory. This means that engineering practicability is obtained.
section_name: Numerical Simulations and Discussions

section_content: In this section, the effect of hydrothermal fluids on S-AUVs τ d f is constructed. The tracking performance of FTFTARSMC algorithms in the TAG mound is verified by numerical simulations. 

As shown in Figure 9, the torpedo-type S-AUV has a length of more than 1700 mm, a diameter of 200 mm, and a weight of 70 kg in air. The propulsion system of the S-AUV consists of eight independent thrusters (including two vertical channel thrusters, two vertical conduit thrusters, two lateral channel thrusters, and two main thrusters) to achieve a fully actuated system given in Figure 10. Then, according to the thrust curves, the saturation limits for the main thrusters, channel thrusters, and vertical conduit thruster are roughly taken as ±600 N, ±500 N, and ±150 N, respectively. In the simulation, the thruster outputs from eight thrusters can be obtained in the form of total control forces and moments τ v at each degree of freedom. 

When the S-AUV tracks the designed trajectory near the vents (below 0. 6 Z max ), it suffers from the disturbance τ d f caused by the hydrothermal plume. Based on the results of the flow field structure for hydrothermal plumes in Refs. [46, 47], the direction of cold seawater near the active vents is roughly horizontal to the hydrothermal vent. In addition, the velocity is inversely proportional to the distance from the current position of S-AUVs to the vent in the range with the vent as the center and D (D ≤ 30 m) as the diameter. However, according to the velocity vector fields of hydrothermal plumes shown in Figure 11, the exact value of the velocity varies depending on the model parameters of plumes. When the S-AUV tracks the designed trajectory near the vents (below 0. 6 max Z ), it suffers from the disturbance d f τ caused by the hydrothermal plume. Based on the results of the flow field structure for hydrothermal plumes in Refs. [46, 47], the direction of cold seawater near the active vents is roughly horizontal to the hydrothermal vent. In addition, the velocity is inversely proportional to the distance from the current position of S-AUVs to the vent in the range with the vent as the center and D (
section_name: Numerical Simulations and Discussions

section_content: ) as the diameter. However, according to the velocity vector fields of hydrothermal plumes shown in Figure 11, the exact value of the velocity varies depending on the model parameters of plumes. 

In addition, there are also some small-scale hydrothermal vents in the TAG active mound. Actually, the effect of small vents can be ignored due to less heat flux and buoyancy flux of the vents. It is worth noting that the S-AUV only suffers from the disturbance d f τ when sailing above the vent, since the entrainments of hydrothermal plumes have a weak influence on the cold seawater below the vent [45]. ) as the diameter. However, according to the velocity vector fields of hydrothermal plumes shown in Figure 11, the exact value of the velocity varies depending on the model parameters of plumes. 

In addition, there are also some small-scale hydrothermal vents in the TAG active mound. Actually, the effect of small vents can be ignored due to less heat flux and buoyancy flux of the vents. It is worth noting that the S-AUV only suffers from the disturbance d f τ when sailing above the vent, since the entrainments of hydrothermal plumes have a weak influence on the cold seawater below the vent [45]. In addition, there are also some small-scale hydrothermal vents in the TAG active mound. Actually, the effect of small vents can be ignored due to less heat flux and buoyancy flux of the vents. It is worth noting that the S-AUV only suffers from the disturbance τ d f when sailing above the vent, since the entrainments of hydrothermal plumes have a weak influence on the cold seawater below the vent [45]. 

According to the above analysis, the disturbance τ d f is modeled on the earth-fixed frame (E-X 0 Y 0 Z 0 ), which can be defined as

where 2 is the distance from S-AUVs to the vent. The resultant force coefficient F d f denotes a small positive constant, since the maximum flow outside plumes is no more than 10 cm/s. For verifying the robustness of FTFTARSMC algorithms, the disturbance τ d f is appropriately amplified. The proposed fixed-time fault-tolerant adaptive robust sliding mode controller (FTF-TARSMCer) of S-AUVs is applied to the designed trajectory and topography in the TAG hydrothermal mound. In order to evaluate the ability of control algorithms in estimating uncertainties and eliminating saturated constraints, two modified control algorithms based on control laws (72)-(75) will be designed and compared with the proposed FTFTARSMCer, i. e., the proposed controller without the ADS (CWA) and the proposed controller without observers (CWO). The initial variables η(0) and υ(0) are selected as η 1 (0) = [126, 0, 35] T m, η 2 (0) = [0, 0, 0] T rad, υ 1 (0) = [0, 0, 0] T m/s, and υ 2 (0) = [0, 0, 0] T rad/s. To guarantee the stable operation of sonars, the speed of S-AUVs needs to be controlled within a certain range. The main control parameters are selected in Table 1. The proposed fixed-time fault-tolerant adaptive robust sliding mode controller (FTF-TARSMCer) of S-AUVs is applied to the designed trajectory and topography in the TAG hydrothermal mound. In order to evaluate the ability of control algorithms in estimating uncertainties and eliminating saturated constraints, two modified control algorithms based on control laws (72)-(75) will be designed and compared with the proposed FTF-TARSMCer, i. e., the proposed controller without the ADS (CWA) and the proposed controller without observers (CWO). The initial variables ( )

. To guarantee the stable operation of sonars, the speed of S-AUVs needs to be controlled within a certain range. The main control parameters are selected in Table 1. In order to reflect the lumped uncertainties (10) in the dynamic model, an additive variance of 40% of the nominal values is considered, i. e., ∆N(η, υ) = 0. 4 M. υ + Ĉ(υ)υ + D(υ)υ + ĝ(η). The disturbances τ dη are chosen as τ dη,1 = 4 + 3 sin(0. 3t)N, τ dη,2 = 5 sin(0. 2t)N, τ dη,3 = 2 + 4 sin(0. 1t)N, τ dη,4 = 3 + 2 cos(0. 1t) Nm, τ dη,5 = 4 cos(0. 2t)Nm, τ dη,6 = 4 + cos(0. 3t)Nm. 

(94)

The conditions of actuator faults are set in Table 2. The trajectory tracking simulation of S-AUVs is created under actuator faults, and the results are presented in Figures 12 13 14 15 16. Figure 12 displays the trajectory tracking results under the CWA, CWO, and FTFTARSMCer, respectively. Further, there is little difference between tracking results due to the long-distance trajectory. It can be obviously seen from Figures 13 and 14 that the S-AUV based on the FTFTARSMCer can quickly and stably track the desired trajectory. The CWO can achieve similar control performance on tracking but shows poor performance in terms of control accuracy and stability. Compared with the CWA, the FTFTARSMCer has remarkable performance in the control stability, the extreme value of control input, and the convergence of tracking errors. When the fault occurs at t 0 = [15; 15; 15; 15; 15; 15]s, the disturbance observer will accurately and quickly estimate the effect of the fault on S-AUVs to improve the robustness of S-AUV systems as shown in Figure 15. From Figure 16, the ADS plays a significant role in the whole control system. From the simulation results, the tracking performance of the fixed-time stable system of S-AUVs is accepted. In addition, in order to thoroughly compare the performance differences between the CWA, CWO, and FTFTARSMCer, the integral of time multiplied by the absolute error (ITAE), integral time absolute square error (ITASE), and integral of the absolute square error (IASE) are introduced as follows:

The error results for the CWA, CWO, and FTFTARSMCer in the actuator fault case are shown in Table 3. From this figure, the FTFTARSMCer provides the minimum value of error results and, thereby, it gives rise to a better performance than CWO and CWA. In addition, in order to thoroughly compare the performance differences between the CWA, CWO, and FTFTARSMCer, the integral of time multiplied by the absolute error (ITAE), integral time absolute square error (ITASE), and integral of the absolute square error (IASE) are introduced as follows:

The error results for the CWA, CWO, and FTFTARSMCer in the actuator fault case are shown in Table 3. From this figure, the FTFTARSMCer provides the minimum value of error results and, thereby, it gives rise to a better performance than CWO and CWA.
section_name: m D ≤

section_content: In the near-bottom explorations of the TAG mound, the S-AUV will suffer from timevarying disturbances, model uncertainties, actuator faults, and input saturations. It is caused by the S-AUV dynamic characteristics, intense hydrothermal activities, and complex terrain. An FTFTARSMC method is proposed to handle these issues. By utilizing the ESO technique, an adaptive FT disturbance observer is constructed to eliminate the effect of lumped uncertainties. To obtain a remarkable fixed-time tracking performance under actuator faults, a continuous fixed-time fast sliding mode surface is designed. It is theoretically verified that fixed-time stabilization for the trajectory tracking system of S-AUVs can be realized with continuous fixed-time adaptive robust controllers. In addition, the simulated terrain is built according to the actual TAG hydrothermal mound, and a smooth trajectory model is constructed by cubic spline interpolation. Finally, comprehensive simulations performed on an actual S-AUV model are given to validate the effectiveness and superiority of the FTFTARSMCer. 

In the future, much work needs to be performed, according to Refs [7, 23]. We will devote to constructing the semi-physical simulation platform in Unity with C#. The state constraints, model mismatches, and uncertainties of the underactuated AUVs will be considered. Further, we will develop the controller with less computational complexity and computational resources. 

After substituting Equations ( 70 ), (A6), and (A7) into (A5), we can obtain. 

Similarly, from the Lemma 3, one obtains

Applying the above inequalities to (A8) yields. 

According to the designed error a = a *â ≥ 0 and a =| a|, we obtain

(A11) 1/2 < (λ 2 + 1)/2 holds due to 0 < λ 2 < 1. It is apparent that with the existence of the constant 0 < h 1 < 1, then the following inequality condition is satisfied: V 1 can be written as. 

V ≤ s T -ẑ 2 + J M-1 N(x 1, x 2 ) + s T J M-1 ∆τλ min (k 3 )6 (A14) 1/2 < (λ 1 + 1)/2 holds due to λ 1 > 1. With the existence of the constant 0 < h 2 < 1, then the following inequality condition holds:

Combining with the inequality condition (A15), we have. 

V ≤ s T -ẑ 2 + J M-1 N(x 1, x 2 ) + s T J M-1 ∆τλ min (k 3 )6

According to the characteristics of the ADS (70) and (71), the detailed discussion about the stability analysis of S-AUV trajectory tracking control systems in the TAG mound will be shown separately in the following cases. 

( (λ 2 + 1)/2 < (λ 2 + 1)/(2λ 2 ) holds due to 0 < λ 2 < 1. It is apparent that with the existence of the constant 0 < h 3 < 1, then the following inequality condition is satisfied:

Synthesizing (A19) and (A20) yields

Applying the observer output to the above inequality and using the Lemma 1, we obtain. 

where (A23)

and χ 1 = â. In addition, the matrix k 4 satisfies Equation (2) when 1 < ξ < µ 2, ĥ(ξ) ≤ 1 and ξ are bounded. According to Lemma 4, we obtain

Substituting (A17), (A18), (A20), and (A24) into (A16), one obtains. 

V 1 ≤λ min (k 3 )6 λ min (k 1 )6

where a 2 = 2 1+λ 1 2 min λ min (k 3 )6

1-λ 1 2

-1 λ 1 +1, 1 h 2 (1+λ 1 ) s, λ min (k 1 )6

and χ 2 = â + λ 1 λ 1 +1 J M-1 ∆τ λ 1 +1 λ 1. In addition, the matrix k 3 satisfies

Equation ( 3 ) when ξ ≤ min{µ 2, 1}, ĥ(ξ) ≤ 1 and ξ are bounded. Based on Lemma 4, the following inequality holds:

(λ 1 + 1)/2 > (λ 1 + 1)/(2λ 1 ) holds due to λ 1 > 1. Similarly, with the existence of the constant 0 < h 4 < 1, then the following inequality condition holds:

Substituting (A27)-(A29) and (A24) into (A16), we obtain. 

V 1 ≤λ min (k 3 )6

1-λ 1 2

+ χ 3 (A30) where and

λ 2. The matrices k 2 and k 3 satisfy

Then, with synthesizing the inequality conditions (A26) and (A32), it is concluded that inequalities (A22), (A25), and (A30) hold if the appropriate matrices k 2, k 3, k 4, k 5 and parameters λ 1, λ 2 are designed to guarantee that conditions (A23), (A31), and (A32) are satisfied. Thus, synthesizing the inequalities (A22), (A25), and (A30), we can conclude that. 

V 1 ≤ -aV 1 (1+λ 1 )/2 -bV 1 (1+λ 2 )/2 + χ (A33)
section_name: Conclusions

section_name: Institutional Review Board Statement: Not applicable.

section_content: Funding: This work is supported in part by the National Natural Science Foundation of China under Grant ( 51979058, 62173104, 51979057, 51779058 ).

section_content: Informed Consent Statement: Not applicable.
section_name: Institutional Review Board Statement: Not applicable.

section_content: Data Availability Statement: Not applicable.

section_name: Conflicts of Interest:

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest:

section_content: The detailed derivation of Equation (77). We consider the following Lyapunov function:

Substituting Equations ( 56 ), (57), and (20) into the time derivative of

Applying the control law ( 72 )-( 75 ) into (A2) results in

Substituting the adaptive law into (A3), we have s T s λ 1 = -
section_name: Appendix A