Preprint

Article

Research on Practical Path Tracking Control of Autonomous Underwater Vehicle Based on Constructive Dynamic Gain Controller

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

Fang Wang,Yudong Peng,

Longchuan Guo^*

Fang Wang,Yudong Peng,

Longchuan Guo^*

This version is not peer-reviewed

Submitted:

04 August 2023

Posted:

04 August 2023

You are already at the latest version

Alerts

Abstract

This paper studies the trajectory tracking control of an autonomous underwater vehicle (AUV) based on a stochastic uncertain nonlinear system. We investigate the time-varying gain adaptive control method looking for possible approaches to alleviate the heavy computational burden. Under nonlinear growth, conditions satisfy polynomial growth conditions. These two problems are resolved the fast response time and good path tracking, respectively. Novel adaptive algorithms are developed exploiting the dynamic properties of the AUV motion. By appropriately turning the original controller design problems into parameters choosing problems and then solving them in a functional time-varying observer technical theorem. In order to deal with the station error of system converges to arbitrarily small domains with stochastic uncertain disturbance. A coordinate transformation is proposed for all system states to meet boundedness. We show that the convergence of the AUV trajectory errors can be guaranteed by the proposed contraction constraints in the stability analysis, closed-loop stability is proved, and the system is asymptotically probabilistic in the global scope. They are taking advantage of the guaranteed stability. Extensive simulation studies on the AUV model demonstrate the effectiveness and robustness of the proposed approach. A real-time time-varying gain constructive control strategy is further developed for the Hardware-in-the-loop simulation; the controller design results are imported into the AUV actuator model to verify the effectiveness of the controller design.

Keywords:

Subject: Engineering - Control and Systems Engineering

1. Introduction

An Autonomous Underwater Vehicle (AUV) is widely used in underwater exploration and scientific research. In AUV research, the control system of underwater thrusters can be said to be its neural center and the most complex technology, so the control of underwater thrusters is one of the most challenging problems in the field of control. Due to the crucial role played by AUVs in deep-sea environmental resource exploration. In recent decades, tracking control, the state maintenance, and path planning of AUVs have received attention from researchers, leading to an urgent need to develop new control theories in the theoretical and engineering fields of AUV output feedback control, making them fit for complex systems affected by various factors such as randomness, nonlinearity, and time delay. In addition, the control theory system built under the framework of certainly systems is no longer applicable to stochastic nonlinear systems [1,2,3,4,5,6,7,8,9,10,11]. Therefore, it is a great challenge to solve the output feedback control problem of stochastic nonlinear systems by using stochastic control theory and related mathematical tools. In engineering practice, there is a high probability that the input and output of a controlled system are random signals. Therefore, the stochastic nonlinear system that combine the nonlinear and stochastic characteristics are a research area of considerable interest in the control theory community.

The application of trajectory tracking control is essential for autonomous underwater robots. However, due to the complex system, uncertain working environment, and highly coupled nonlinear characteristics of underwater robots [1,2], the controller design of AUV becomes very challenging. In the past decades, the design of trajectory tracking controller [3] for underwater vehicles has been a problem worthy of attention. Traditional PID, LQR, and Kalman filtering [4] can also control the trajectory of the underwater robot. In addition, nonlinear PID also has good performance in control. However, when the target trajectory being tracked is nonlinearity, the nonlinearity of the curved trajectory makes nonlinear PID no longer applicable in vehicle driving and class linear control techniques. Feedback linearization [5] presents a powerful tool for dealing with nonlinear characteristics. However, applying feedback linearization to the AUV requires a system model with highly accurate hydraulic dynamic coefficients [6]. In this context, the adaptive [7,8,9] Lyapunov method has become the mainstream method for AUV trajectory control. In Article [10], the dynamic controller is used for the first time in the backstepping control technique. AUV output backstep control can be found in Article [11,12]; another popular AUV trajectory tracking control method is sliding mode control. It is well known that this method has an excellent robust effect when the parameters are uncertain, but the sliding effect of sliding mode control tends to appear discrete. In order to reduce this phenomenon, sliding mode control is often used in combination with control methods such as robust control, adaptive control, and PID [13,14,15].

Nonetheless, the above control methods have a general deficiency the inability to deal with the stochastic characteristics and uncertain disturbances of the system. For AUVs, stochastic and uncertain are very common [16,17], which the actuator must suffer in the working environment or the system’s characteristics. Adaptive control is featured in dealing with this randomness and uncertainties [18], which proposes a robust control structure in the face of this wide range of control problems. In addition, adaptive control can solve complex nonlinear problems for solving the control problems of AUV dynamic systems. So far, we do not have a very effective control method for trajectory tracking control problems. In article [19], a self-correcting adaptive method is proposed, a complete method to solve the problem that its dynamic parameters and even the model structure often change. The trajectory tracking control algorithm for AUV stochastic uncertain nonlinear system is proposed. In paper [20], where the concept of stochastic and uncertain is creatively introduced. In article [21,22], a comprehensive trajectory tracking control and path planning problem is researched. A unified optimization framework [23] is developed for the problem of AUVs combined motion control.

Although an adaptive control based on an integral high-order sliding mode concept provides a good algorithm for nonlinear systems with random and uncertain disturbances [13], when solving, this method heavily burdens the computational bottleneck. In theoretical research, the computational time is often ignored, but in practice, the computational difficulty increases exponentially with the increase of stochastic and uncertain interference problems. Due to the short sampling period, many strategies such as iterative methods, precomputation, and numerical continuation [24] have been proposed to decrease the computational complexity and reduce the computational time. By exploring the motion characteristics of underwater robots, a dynamic-static combined with a high-gain observer is successfully applied to the output feedback adaptive control algorithm. However, the stability proof of the closed-loop system in paper [25] is not provided. Since the dynamic gain control, algorithm solves the trajectory tracking control problem of the underwater robots. The implicit coupling of the system state and the control signal complicates the stability analysis of the closed-loop system. It is urgent to find a better stability analysis method.

This paper hope to provide a method to eliminate the influence of random and uncertain disturbances on AUV trajectory tracking and simultaneously ensure the closed-loop system stability. Here, the reference augmentation technique is applied to modeling the AUV system so that the coupled motion between the systems is weakened [26]. Then the dynamic gain method is studied to reduce the computational stress of stochastic uncertain nonlinear systems. Since the novel dynamic gain, adaptive control algorithm was proposed, the complex control problem has been changed into a parameter selection and construction problem through dynamic gain. The computational complexity is significantly reduced and the computation time is significantly reduced. The contributions of this article are as follows:

1. Aiming at the problem of AUV trajectory tracking control, a well-known stochastic uncertain nonlinear dynamic adaptive control algorithm is provided to study the AUV motion’s dynamic characteristics.

2. A novel dynamic-static combination of high-gain observers is proposed, which dramatically reduces the computational complexity of the control algorithm and transforms the controller design problem into a parameter calculation and selection problem.

3. Investigate fundamental properties of closed-loop systems. The control time and control error of AUV trajectory tracking is remarkably reduced, the control effect is incredibly optimized, and the sensitivity is improved.

4. Numerical simulations and additional experiments reveal that the proposed dynamic-static high-gain adaptive control algorithm has excellent robust performance against stochastic and uncertain disturbances.

The rest of this paper is organized as follows: The second section describes the AUV motion model. The third section is the dynamic and static high-gain adaptive control algorithm. Section 4 presents simulation studies and additional experiments. The fifth section is the conclusion and prospects.

2. Description of Dynamic Modelling of Robot

In this section, the dynamics analysis of the AUV is carried out. According to the analysis of the power, resistance, and inertial force when moving underwater, we can obtain the general dynamic equation. Considering the complexity and uncertainty of the underwater environment [2,27], a stochastic process was added to the model. Based on the established underactuated underwater robot model, as shown in Figure 1.

The 6-DOF dynamic model of the AUV can be expressed as follows:

M \dot{v} + C (v) v + D (v) v + g (ϵ) = F

(1)

where

M = M_{A} + M_{R B}

is the inertia mass matrix, which contains AUV additional quality.

C (v)

are Coriolis force matrix of the underwater vehicles.

D (v)

is for underwater vehicle fluid damping matrix.

g (ϵ)

stands for the gravity of the underwater vehicle in operation and the restoring force (moment) matrix generated by the buoyancy.

F = {[F_{x}, F_{y}, F_{z}, F_{k}, F_{m}, F_{n}]}^{T}

stands for the resultant forces and moments.

\dot{v} = {[\dot{u}, \dot{v}, \dot{ω}, \dot{p}, \dot{q}, \dot{r}]}^{T}

are the linear and angular accelerations of the body (moving) frame in the direction of pitch, roll and heave,

v = {[u, v, ω, p, q, r]}^{T}

stands for the linear and angular velocity with respect to body (moving) frame,

u

– surge velocity,

v

– sway velocity,

ω

– heave velocity,

p

– roll rate,

q

– pitch rate,

r

– yaw rate.

ϵ = {[x, y, z, α, β, γ]}^{T}

are the position and orientation in inertial(fixed) frame,

x

– surge position,

y

– sway position,

z

– heave position,

α

– roll angle,

β

– pitch angle,

γ

– yaw angle.

According to the above, the six-degree-of-freedom (6-DOF) model of the underwater robot has complex nonlinearity and state coupling. If the paper wants to design a controller with six degrees of freedom, it will pose a great challenge to the controller design and physical characteristics. Here, we decompose the 6-DOF motion model into two kinematic models, which is linear velocity variables

v_{1} = {[u, v, ω]}^{T}

and angular velocity variables

v_{2} = {[p, q, r]}^{T}

. In the earth-fixed frame, the velocity can generally be decomposed into two kinematic models, underwater robot position

ϵ_{1} = {[x, y, z]}^{T}

and underwater robot orientation

ϵ_{2} = {[α, β, γ]}^{T}

, which can highly simplify the AUV model. In this paper, we only consider the velocity variables

v_{1} = {[u, v, ω]}^{T}

and linear motion

ϵ_{1} = {[x, y, z]}^{T}

of the AUV underwater.

The relationship between the fixed frame and body frame in linear velocity is as follows:

{\dot{ϵ}}_{1} = J (ϵ) v_{1}

(2)

Where J (ϵ) J (ϵ) = (\begin{matrix} c γ c β & - s γ c α + c γ s β s α & s γ s α + c γ c β s α \\ s γ c β & c γ c α + s γ s β s α & - c γ s α + s γ s β c α \\ - s β & c β s α & c β c α \end{matrix})

is the kinematic transformation matrix of the following form:

Where,

s = s i n

c = c o s

and

ϵ_{2} = {[α, β, γ]}^{T}

is the angle between the surge, sway and heave direction and the earth frame, respectively. When

β = \pm 90^{°}

, this transformation is undefined and the quaternion method has to be considered. While, most robots are designed to operate at pitch angles below

\pm 90^{°}

; Thus, this restriction is of no great significance here. In order to better understand and analyze the motion state, we will study the AUV system based on the earth reference frame. In order to unify the signal states, using Eq. (2) for coordinate transformation

(ϵ_{1}, v_{1}) \overset{μ}{\to} (ϵ_{1}, {\dot{ϵ}}_{1})

, we obtain:

(\begin{matrix} ϵ_{1} \\ {\dot{ϵ}}_{1} \end{matrix}) = [\begin{matrix} I & 0 \\ 0 & J (ϵ) \end{matrix}] (\begin{matrix} ϵ_{1} \\ v_{1} \end{matrix})

M_{ϵ_{1}} {\ddot{ϵ}}_{1} + C_{ϵ_{1}} {\dot{ϵ}}_{1} + D_{ϵ_{1}} {\dot{ϵ}}_{1} + g_{ϵ_{1}} = F_{ϵ_{1}}

(3)

The coordinate transformation is a global diffeomorphism analogous to a similarity transformation in the linear system. The dynamic underwater vehicle model with the earth fixed reference frame is as follows:

where

\begin{array}{l} M_{ϵ_{1}} = J (ϵ)^{- T} M J (ϵ)^{- 1} \\ C_{ϵ_{1}} = J (ϵ)^{- T} (C (ν) - M J (ϵ)^{- 1} \dot{J} (ϵ)) J (ϵ)^{- 1} \\ D_{ϵ_{1}} = J (ϵ)^{- T} D (ν) J (ϵ)^{- 1} \\ g_{ϵ_{1}} = J (ϵ)^{- T} g (ϵ) \\ F_{ϵ_{1}} = J (ϵ)^{- T} F \end{array}

(4)

We need some assumptions as follows.

Assumption 1. In this paper, only the velocity between the self-position and the fixed frame is considered, the velocity is

\dot{l} = {({\dot{x}}^{2} + {\dot{y}}^{2} + {\dot{z}}^{2})}^{\frac{1}{2}}

. Suppose the origin state

ϵ_{1} (0) = [0 0 0]

and

{\dot{ϵ}}_{1} (0) = [0 0 0]

Assumption 2. In this paper, the velocity

ϵ_{1} = [x y z]

and the angle

ϵ_{2} = [α β γ]

are known using the sensors.

In this paper, we focus on the problem of trajectory tracking control for a class of stochastic uncertain nonlinear systems. For the trajectory tracking control problem of the underwater robot under nonlinear dynamics, to better design the control algorithm, we convert the dynamic model of the system into a broad numerical model as follows.

\begin{matrix} D l_{1} = l_{2} d t + φ_{1} (t, l (t), u) d ω \\ d l_{2} = u d t + φ_{2} (t, l (t), u) d ω \\ y = l_{1} - y_{r} \end{matrix}

(5)

In the stochastic nonlinear systems

l = {(l_{1}, l_{2})}^{T} \in R^{2}

u \in R

y \in R

are the states, input and output of system, respectively;

y_{r}

is the target trajectory to be tracked. Where state

l_{2}

and

y_{r}

are unmeasurable. The stochastic process of the system is introduced in this paper:

ω

is an m – dimensional standard Wiener process defined on the complete probability space

(Ω, Γ, P)

, where

Ω

is the sample space,

Γ

is the filter, and

P

is the probability measure. Nonlinear term

φ_{i} : R^{+} \times R \times R^{2}, i = 1,2

is satisfied continuous polynomial growth conditions and local Lipschitz.

In this paper, we assume that system (5) satisfies the following assumptions and implements output feedback tracking control.

Assumption 3. There exists a positive constant

c_{0} \geq 0

and a known integer

q \geq 1

such that the following inequality holds

|φ_{i} (t, l, u)| \leq c_{0} (1 + {|l_{1}|}^{q}) (|l_{1}| + |l_{i}|) + c_{0}, i = 1,2

(6)

This assumption shows that the output polynomial growth rate governed this system (5).

Assumption 4. The target trajectory

m a x (| y_{r} (t) |, | {\dot{y}}_{r} (t) |) \leq c_{1}

c_{1}

is known constants. According to the description of the above system, this section’s goal is to design a time-varying gain constructive control algorithm: for any given in advance tolerance

δ > 0

. In the closed-loop system, all the states are well-defined and globally bounded on

[0, + \infty]

in finite time

T_{δ} > 0

, for instance

| y (t) | = | l_{1} (t) - y_{r} (t) | \leq δ, \forall t \geq T_{δ}

Remark 1. It should be pointed out that the output feedback tracking control problems studied in previous literatures [28,29] are mainly aimed at nonlinear systems or nonlinear systems with parameter uncertainty and unknown control direction. The system (5) studied in this paper is a stochastic uncertain nonlinear system after introducing random factors. At the same time, system (5) is observable by Assumption 3, which depends on the unmeasurable state and nonlinear growth conditions satisfying the polynomial growth condition, according to Assumption 4. It can be seen that for the reference trajectory

y_{r}

, its upper bound value and the upper bound value of the derivative are given, which means that it is not necessary to give a specific description function or give more information to the reference trajectory

y_{r}

Based on this assumption 3, 4 and (5), for the

l_{1} = y + y_{r}, i = 1, 2

. We have:

\begin{matrix} |φ_{i} (t, l, u)| \leq c_{0} (1 + {|y + y_{r}|}^{q}) (|l_{1}| + |l_{2}| + \dots + |l_{i}|) + c_{0} \\ \leq c (1 + {|y|}^{q}) (|l_{1}| + |l_{2}| + \dots + |l_{i}|) + c \end{matrix}

(7)

where

c = c_{0} m a x {1 + 2^{q - 1} c_{1}^{q}, 2^{q - 1}}

is a known constant.

3. Dynamic and Static High Gain Adaptive Control Algorithm

In this section, we propose a dynamic gain observer strategy to achieve dynamic adaptive tracking control, so that the complexity can be greatly reduced, which will study the stochastic uncertain nonlinear underwater robot system in two steps. The first step is the dynamic gain observer and controller design of the AUV system; the second step is stability analysis and implementation of the AUV system.

3.1. Dynamic Gain Observer and Controller Design

First step, this section is to introduce the following state transformation:

\begin{array}{l} z_{1} = y = l_{1} - y_{r} \\ z_{2} = l_{2} \end{array}

(8)

So, we can get the system (9).

\begin{array}{l} {\dot{z}}_{1} = z_{2} + f_{1} (t, z, u) \dot{ω} \\ {\dot{z}}_{2} = u + f_{2} (t, z, u) \dot{ω} \\ y = z_{1} \end{array}

(9)

where

\begin{array}{l} f_{1} (t, z, u) = φ_{1} (t, z_{1} + y_{r}, z_{2}, u) - \frac{d y_{r}}{d ω} \\ f_{2} (t, z, u) = φ_{2} (t, z_{1} + y_{r}, z_{2}, u) \end{array}

(10)

Since the states of the stochastic system are unmeasurable except for

l_{1}

, which is a measurable state, the time-varying observers are first constructed:

\begin{matrix} {\dot{\hat{z}}}_{1} = z_{2} + K^{1} h_{1} (z_{1} - {\hat{z}}_{1}), \\ {\dot{\hat{z}}}_{2} = u + K^{2} h_{2} (z_{1} - {\hat{z}}_{1}), \end{matrix}

(11)

Meantime, we choose design parameters

h_{I}, i = 1, 2

, are constants and satisfying Hurwitz condition and

s^{2} + h_{1} s + h_{2}

. The

K = A B (t)

is a dynamic gain parameter, which

A > 1

is a constant and

B (t)

is a time-varying function updated by

\begin{array}{l} \dot{B} (t) = - μ_{1} B^{2} + μ_{2} {(1 + {|y|}^{q})}^{2} B \\ B (0) = 1 \end{array}

(12)

Theorem 1. So far, we can conclude that all closed-loop system states are bounded on

[0, T]

, then

T

must be infinite. Otherwise, at least one state will escape at

T

, contradicts the continuity of the closed-loop system. For a class stochastic uncertain nonlinear system (5),

y = l_{1} - y_{r}

is the output, and the nonlinear term of the system meets Assumptions 3 and 4. According to choose the suitable parameters

K_{I}, g_{I}, i = 1, 2, A, μ_{1}

and

μ_{2}

and based on observer (11), the output feedback controller is designed as follows:

u = - (K^{2} a_{1} {\hat{z}}_{1} + K a_{2} {\hat{z}}_{2})

(13)

where the parameter

a = {[a_{1}, a_{2}]}^{T}

is a constant vector and satisfying Hurwitz condition and

h_{i} = a_{3 - i}

. Next, we will to prove that practical global tracking can be achieved when

A

is large enough. Furthermore, practical path tracking can be achieved. For any

δ > 0

, there exists a finite time

T_{δ}

, such that

| y (t) | \leq δ, \forall t \geq T_{δ}

Remark 2. Through equation (8), we can turn complex systems into simple ones. According to the unmeasurable state of the system, we introduce a dynamic gain observer (11) to construct the observation time of the system state. Unlike the previous observers, the observer gain in this paper is an observer combined with dynamic

B (t)

and static

A

gain parameters. The observed effect is closer to the actual value of the system, which provides a good condition for the following controller to control the error variable. On account of the insufficient information on the system and the tracking signal and the system’s instability after introducing random factors, the common tracking control methods, such as the MPC method, can no longer solve the problems in this chapter. For these challenging problems, this article will pay close attention to stochastic nonlinear systems satisfying assumptions 3 and 4 by constructing an output feedback dynamic adaptive tracking controller ground on the combination of the static and dynamic gain observer to achieve the original system (5) the actual tracking control target.

3.2. Stability Analysis and Implementation

In this subsection, the main results of the paper are presented and harshly proved. For the stability analysis, we need the following scaling transformation for the closed-loop system.

The state estimation error is defined as

e_{i} = z_{i} - {\hat{z}}_{I}, (i = 1, 2)

. Then the dynamics of

e_{I}

satisfy:

\begin{array}{l} {\dot{e}}_{1} = e_{2} - K^{1} h_{1} e_{1} + f_{1} (t, z, u) \dot{ω} \\ {\dot{e}}_{2} = - K^{2} h_{2} e_{1} + f_{2} (t, z, u) \dot{ω} \end{array}

(14)

To facilitate the design of the controller, we need to transform the estimated state

{\hat{z}}_{I}

and the error state

e_{I}

as follows:

\begin{matrix} ε_{i} = \frac{e_{i}}{K^{b + i - 1}} (i = 1,2) \\ τ_{i} = \frac{{\hat{z}}_{i}}{K^{b + i - 1}} (i = 1,2) \end{matrix}

(15)

\begin{array}{l} \dot{ε} = K H ε - \frac{\dot{K}}{K} C_{b} ε + G (a, K) \dot{ω} \\ \dot{τ} = K H_{b} τ - K h ε_{1} + \frac{\dot{h}}{h} C_{b} τ \end{array}

(16)

where

0 < b < \frac{1}{4 q}

is a known constant,

ε = {(ε_{1}, ε_{2})}^{T}, τ = {(τ_{1}, τ_{2})}^{T}, h = {(h_{1}, h_{2})}^{T}

. According to (15), stochastic nonlinear systems (9) and (11) can be converted as follows:

Meantime,

G (z, K) = {[\frac{f_{1}}{K^{b}}, \frac{f_{2}}{K^{b + 1}}]}^{T}, H = (\begin{matrix} - h_{1} & 1 \\ - h_{2} & 0 \end{matrix}), H_{b} = (\begin{matrix} 0 & 1 \\ - h_{1} & - h_{2} \end{matrix}), C_{b} = (\begin{matrix} b & 0 \\ 0 & b + 1 \end{matrix})

Then, we discuss the gain

B

and the boundedness of the states

ε

and

τ

. Above all, choose the proper parameters

r_{1}, r_{2}, r_{3}, r_{4}

, so that the relationship between the positive definite matrices

P, Q

and the matrix

H, H_{b}, C_{b}

satisfies:

\begin{matrix} H^{T} P + P A H \leq - I_{2} & r_{1} I_{2} \leq C_{b} P + P C_{b} \leq r_{2} I_{2} \\ H_{b}^{T} Q + Q H_{b}^{T} \leq - 2 I_{2} & r_{3} I_{2} \leq C_{b} Q + Q C_{b} \leq r_{4} I_{2} \end{matrix}

(17)

Let

V (ε (t), τ (t)) = ε^{T} (t) P ε (t) + τ^{T} Q τ (t)

be the Lyapunov function candidate. Then along the trajectories of (16), the time derivative of

V

[0, T]

satisfies the following inequality:

L V = + 2 ε^{T} G P - {|ε|}^{2} K - 2 ε^{T} \frac{\dot{K}}{K} (C_{b} P + P C_{b}) ε + 2 K τ^{T} ε_{1} h - 2 K {|τ|}^{2} - T r (G^{T} (P + Q) G) - \frac{\dot{K}}{K} τ^{T} (C_{b} P + P C_{b}) τ

(18)

From Remark 2, we can know

h > 1

, and from (12),

\frac{\dot{h}}{h} = \frac{\dot{D}}{D} = - μ_{1} D + μ_{2} {(1 + {| y |}^{I})}^{2}

\begin{matrix} - \frac{\dot{K}}{K} ε^{T} (C_{b} P + P C_{b}) ε \leq r_{2} μ_{1} B {|ε|}^{2} - r_{1} μ_{2} {(1 + {|y|}^{q})}^{2} {|ε|}^{2} \\ - \frac{\dot{K}}{K} τ^{T} (C_{b} Q + Q C_{b}) τ \leq r_{4} μ_{1} B {|τ|}^{2} - r_{3} μ_{2} {(1 + {|y|}^{q})}^{2} {|τ|}^{2} \end{matrix}

(19)

Then, from this and (17), we have

By Assumption 4, (7), (8), (15) and the fact

h > 1,

we have

|G| = {({|G_{1}|}^{2} + {|G_{2}|}^{2})}^{\frac{1}{2}} = {({(\frac{f_{1}}{K^{b}})}^{2} + {(\frac{f_{2}}{K^{b + 1}})}^{2})}^{\frac{1}{2}}

(20)

Where

\begin{array}{l} |G_{1}| \leq c [(1 + {|y|}^{q}) (|ε| + |τ| + c_{1}) + 1] + c_{1} \\ |G_{2}| \leq c [(1 + {|y|}^{q}) (\sqrt{2} (|ε| + |τ|) + c_{1}) + 1] + c_{1} \end{array}

(21)

Split and enlarge the terms in formula (18), and we can get

\begin{matrix} 2 {G ε}^{T} P \leq 2 ‖ε‖ ‖P‖ (c (1 + {|y|}^{q}) (c_{1} + (\sqrt{2} (|ε| + |τ|))) + c_{1} + c) \\ \leq ‖P‖ c (1 + c_{1}) {({|ε|}^{2} + {|τ|}^{2}) (1 + {|y|}^{q})}^{2} + ‖P‖ (5 c + c_{1}) {|ε|}^{2} + ‖P‖ (c c_{1} + c + c_{1}) \\ 2 K ε_{1} h τ^{T} \leq K ({{|τ|}^{2} + {|ε|}^{2} |h|}^{2}) \end{matrix}

(22)

Therefore, combining equations (18)-(22), the system Ito differential equation can be obtained as:

\begin{matrix} \begin{matrix} L V \leq - K {|ε|}^{2} + c ‖P‖ (1 + c_{1}) {(1 + {|y|}^{q})}^{2} ({|ε|}^{2} + {|τ|}^{2}) + ‖P‖ (5 c + c_{1}) {|ε|}^{2} + ‖P‖ (c c_{1} + c + c_{1}) \\ + α_{1} B ({r_{2} |ε|}^{2} + r_{4} {|τ|}^{2}) - α_{2} {(1 + {|y|}^{q})}^{2} ({r_{1} |ε|}^{2} + r_{3} {|τ|}^{2}) - 2 K ({|τ|}^{2} + τ^{T} ε_{1} h) + T r (G^{T} (P + Q) G) \\ \leq - K {|ε|}^{2} + c ‖P‖ (1 + c_{1}) {(1 + {|y|}^{q})}^{2} ({|ε|}^{2} + {|τ|}^{2}) + ‖P‖ (5 c + c_{1}) {|ε|}^{2} + ‖P‖ (c c_{1} + c + c_{1}) \end{matrix} \\ + r_{2} α_{1} B {|ε|}^{2} - α_{2} {(1 + {|y|}^{q})}^{2} ({r_{1} |ε|}^{2} + {|τ|}^{2}) - 2 K {|τ|}^{2} + r_{4} α_{1} B {|τ|}^{2} + K ({|τ|}^{2} + {|h|}^{2} {|ε|}^{2}) \\ - (‖P‖ + ‖G‖) {(c (1 + {|y|}^{q}) (\sqrt{2} (|ε| + |τ|) + c_{1}) + c + c_{1})}^{2} \end{matrix}

(23)

Where the last term of formula (23) can be enlarged as:

\begin{matrix} \begin{matrix} (‖G‖ + ‖P‖) {(c (1 + {|y|}^{q}) (\sqrt{2} (|ε| + |τ|) + c_{1}) + c + c_{1})}^{2} \\ \leq 2 (‖P‖ + ‖G‖) (c^{2} {(1 + {|y|}^{q})}^{2} {(\sqrt{2} (|ε| + |τ|) + c_{1})}^{2} + {(c + c_{1})}^{2}) \\ \leq 2 (‖P‖ + ‖G‖) (2 c^{2} {(1 + {|y|}^{q})}^{2} (2 {(|ε| + |τ|)}^{2} + {c_{1}}^{2}) + {(c + c_{1})}^{2}) \end{matrix} \\ \leq 2 (‖P‖ + ‖G‖) (2 c^{2} {(1 + {|y|}^{q})}^{2} (4 ({|ε|}^{2} + {|τ|}^{2}) + {c_{1}}^{2}) + {(c + c_{1})}^{2}) \\ = (‖P‖ + ‖G‖) (16 c^{2} {(1 + {|y|}^{q})}^{2} ({|ε|}^{2} + {|τ|}^{2}) + 4 c^{2} {(1 + {|y|}^{q})}^{2} c_{1}^{2} + 2 {(c + c_{1})}^{2}) \end{matrix}

(24)

After finishing, we obtain

\begin{matrix} L V \leq - (K (1 - {|h|}^{2}) - (5 c + c_{1}) ‖P‖ - r_{2} B μ_{2}) {|ε|}^{2} - {({|y|}^{q} + 1)}^{2} (r_{1} μ_{2} - ‖P‖ c (c_{1} + 1) - 16 (‖Q‖ + ‖P‖)) {|ε|}^{2} \\ + (r_{4} μ_{1} B - K) {|τ|}^{2} - {({|y|}^{q} + 1)}^{2} (r_{3} μ_{2} - c ‖P‖ (1 + c_{1}) - 16 c^{2} (‖Q‖ + ‖P‖)) {|ε|}^{2} \\ + ‖P‖ (c c_{1} + c + c_{1}) + (4 c^{2} {(1 + {|y|}^{q})}^{2} c_{1}^{2} + 2 {(c + c_{1})}^{2}) (‖P‖ + ‖Q‖) \end{matrix}

(25)

\begin{array}{l} μ_{2} \geq \frac{c (1 + c_{1}) ‖P‖ - 16 c^{2} (‖Q‖ + ‖P‖)}{r_{1}} \\ μ_{2} \geq \frac{c (1 + c_{1}) ‖P‖ - 16 c^{2} (‖Q‖ + ‖P‖)}{r_{3}} \\ μ_{2} \geq μ_{1} > 0 \end{array}

(26)

According to formula (25), selection

μ_{2}

satisfies:

According to

B > 1

, choose

A

to satisfy:

\begin{array}{l} A \geq \frac{‖P‖ (5 c + c_{1}) + r_{2} μ_{1}}{1 - {|h|}^{2}} \\ A \geq r_{4} μ_{1} \\ A \geq 1 \end{array}

(27)

According to the above selection of parameters, the formula (25) can be transformed into:

\begin{matrix} L V \leq - K ({|τ|}^{2} {+ |ε|}^{2}) + ‖P‖ (c c_{1} + c_{1} + c) + (‖P‖ + ‖Q‖) (4 c_{1}^{2} c^{2} {(1 + {|y|}^{q})}^{2} + 2 {(c + c_{1})}^{2}) \\ \leq - \frac{A V}{m a x (λ_{m a x} (P) + λ_{m a x} (Q))} + (c c_{1} + c_{1} + c) ‖P‖ + (‖P‖ + ‖Q‖) (4 c^{2} c_{1}^{2} {(1 + {|y|}^{q})}^{2} + 2 {(c + c_{1})}^{2}) \end{matrix}

(28)

The relevant parameters are defined as follows:

\begin{array}{l} E_{1} = \frac{A}{m a x (λ_{m a x} (P) + λ_{m a x} (Q))} \\ E_{2} = ‖P‖ (c c_{1} + c + c_{1}) + (4 c^{2} {(1 + {|y|}^{q})}^{2} c_{1}^{2} + 2 {(c + c_{1})}^{2}) (‖P‖ + ‖Q‖) \end{array}

(29)

As previously mentioned, the closed-loop system has a unique solution on the maximal time interval

[0, T)

, where

0 < T \leq + \infty

V (ε (t), τ (t)) \leq V (ε (0), τ (0)) e^{- \frac{E_{1}}{A} t} + \frac{E_{2}}{E_{1}}, t \in [0, T)

(30)

which shows that

τ

are bounded on

[0, T)

Then prove the boundedness of

B

[0, T)

because

ε a n d τ

[0, T)

on the bounded, and

ε_{1} + τ_{1} = \frac{z_{1}}{K^{b}} = \frac{y}{K^{b}}

, we get

|y| \leq {(A B)}^{b} E_{3}

. Then from this, (12) and

0 < b < \frac{1}{4 q}

, it is easily to get that

\begin{matrix} \dot{B} = - μ_{1} B^{2} + μ_{2} {(1 + {|y|}^{q})}^{2} B \\ \begin{matrix} \leq - μ_{1} B^{2} + μ_{2} {(1 + {({(A B)}^{b} E_{3})}^{q})}^{2} B \\ \leq - μ_{1} B^{2} + 2 μ_{2} (1 + {({(A B)}^{b} E_{3})}^{2 q}) B \\ \leq - μ_{1} B^{2} + 2 μ_{2} B (1 + {(A B)}^{\frac{1}{2}} E_{3}^{2 q}) \end{matrix} \end{matrix}

(31)

It means that

B (t) \leq \frac{4 μ_{2}^{2} {(1 + E_{3}^{2 q} \sqrt{A})}^{2}}{μ_{1}^{2}}, \forall t \in [0, T)

(32)

And hence

B

is ounded on

[0, T)

Proof 1 It is easy to verify that the resulting closed-loop system is locally Lipschitz in the open neighborhood of the initial conditions of

(x, \hat{x}, A)

. Therefore, the closed-loop system has a unique solution on a small interval

[0, t)

. Let

[0, T)

be the largest interval for which a unique solution exists, where

0 < T \leq + \infty

. As stated in Lemma 1, where

T = + \infty

, the closed-loop system states are defined on

[0, + \infty)

V (ε (0), τ (0)) e^{- \frac{E_{1}}{A} t} \leq \frac{E_{2}}{E_{1}}

(33)

Next step we can get

V (ε (t), τ (t)) \leq \frac{E_{2}}{E_{1}}, \forall t > T

(34)

It is important to point out here that

E_{1}

contains parameter

A

. From this, it is evident that

λ_{m i n} (P) {|ε|}^{2} + λ_{m i n} (Q) {|τ|}^{2} \leq V (ε (t), τ (t))

(35)

Then we have

ε_{1}^{2} (t) + τ_{1}^{2} (t) \leq \frac{E_{4}}{A}

(36)

Combining (34), (36) and

|y (t)| \leq {(A B)}^{b} E_{3}

, we can get:

|y (t)| \leq A^{b} B^{b} (|ε_{1} (t)|, |τ_{1} (t)|) \leq \sqrt{2 (ε_{1}^{2} (t) + τ_{1}^{2} (t))} \leq \sqrt{\frac{{2 E}_{4}}{A}}, \forall t > T

(37)

B (t) \leq \frac{4 μ_{2}^{2}}{μ_{1}^{2}} {(1 + {(\frac{{2 E}_{4}}{A})}^{q} \sqrt{A})}^{2}

(38)

Then, we know that for all

\forall t \in [T, + \infty)

Combining (37) and in turn

y^{2} (t) = {(A^{b} B^{b})}^{2} {(|ε_{1} (t)|, |τ_{1} (t)|)}^{2} \leq 2 A^{2 b} B^{2 b} (ε_{1}^{2} (t) + τ_{1}^{2} (t)) \leq 2 \frac{16^{b} μ_{2}^{4 b} E_{4}}{μ_{1}^{4 b} A^{1 - 2 b}} {(1 + \frac{2^{q} E_{4}^{q}}{A^{q - \frac{1}{2}}})}^{2 b}

(39)

From this and

1 - 2 b > 1 - \frac{1}{2 q} > 0 (s i n c e 0 < b < \frac{1}{4 q})

, when

t > T

, by choosing

A

big enough to make any small

|y (t)|

, so as to realize the actual path tracking.

Remark 3. By analyzing the dynamic gain parameter

B (t)

of the dynamic observer, we can see that this method turns the control problem into a parameter selection problem through analysis, which can be found that as long as the more prominent the parameter

B (t)

selection area is, the better the control effect will be. A dynamic high-gain state observer is constructed. With the help of Ito stochastic calculus theory, the output feedback actual tracking controller is obtained. By selecting appropriate design parameters, the state and high gain parameters of the closed-loop stochastic nonlinear system are guaranteed to be bounded, and the system tracking error can be converge to zero within a small neighborhood. The study of output feedback tracking control for nonlinear systems is extended to stochastic systems. The practical tracking of output feedback for a class of stochastic nonlinear systems satisfying the growth condition of the output polynomial function is studied for the first time.

4. Underwater Robot Transportation Model Example

In this literature, we will simulate and verify the proposed stochastic uncertain dynamic gain adaptive control algorithm through specific numerical simulation. The dynamic model is converted into a specific mathematical model according to the dynamic characteristics of the underwater vehicle. Consider the following stochastic nonlinear system:

\begin{array}{l} {\dot{l}}_{1} = l_{2} + s i n l_{1} \dot{ω} \\ {\dot{l}}_{2} = u + l_{2} {l n}^{(1 + l_{1}^{2})} \dot{ω} \\ y = l_{1} - y_{r} \end{array}

(40)

In this section, we apply Theorem 1, Assumptions 3 and 4, and Lemma 1 to specify the underwater robot path tracking control system described in Figure 1 with

c_{0} = 0.21, I = 2

, and

c_{1} = 1

. Therefore, an adaptive output feedback controller can be designed using Remark. The related control laws are implemented in the following from

u = - K^{2} a_{1} {\hat{z}}_{1} + K a_{2} {\hat{z}}_{2}

. As well-known in Section 2,

{\hat{z}}_{1}

and

{\hat{z}}_{2}

are the estimations of

y

and

l_{2}

, respectively. These parameters are selected as

b = 0.1, h = {[0.5,1]}^{T}, h = {[1,0.5]}^{T}, μ_{1} = 7.2, μ_{2} = 9.8

and

N = 117

According the initial stations

l_{0} = {[1, - 0.5]}^{T}

and

{\hat{z}}_{0} = 0

the simulation results are shown in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8. These figure show

l, \hat{z}, D

and

\dot{D}

are all bounded, which have to verify the effectiveness of the path tracking controller. Besides, the tracking error

|l_{1} - y_{r}| \leq 0.001

after 0.04s.

5. Concluding and Future Prospects

In this paper, we have studied the application of this dynamic adaptive control to AUV trajectory tracking. The nonlinear characteristics of underwater robot motion are studied, and a novel dynamic-static gain observer strategy is proposed to reduce the computational burden. Numerical simulations and additional experiments validate the efficiency and robustness of the proposed strategy and highlight the advantages of dynamic-static high-gain adaptive control algorithms. It can reduce not only the reaction time of the controlled system and the error of the control state but also has an excellent inhibitory effect on stochastic and uncertain interference.

We discovered a class of nonlinear systems with piecewise linearization in practical engineering and research. Piecewise linearization decomposes a nonlinear system into a combination of finite or infinite linear subsystems. This combination can approximate nonlinear systems. In the control process, piecewise linearization has the same characteristics as nonlinear and simplifies the control process and method. According to the known situation, this method has been applied in many practical projects and has achieved a perfect control effect. However, the current application scope of this method is still minimal and has special conditions for nonlinear systems. The following work will expand the application range of piecewise linearization of nonlinear systems, which will be the next hotspot of nonlinear control research.

Funding

This work was supported by the Stable Supporting Fund of Science and Technology on Underwater Vehicle Technology (JCKYS2022SXJQR-01).

Data Availability Statement

All data generated or analyzed during this research are included in this paper.

6. Acknowledgments

This work was supported by the Stable Supporting Fund of Science and Technology on Underwater Vehicle Technology (JCKYS2022SXJQR-01).

Conflicts of Interest

The author declares that there is no conflict of interest regarding the publication of this paper.

References

Zhang, Y.; Liu, X.; Luo, M.; Yang, C. MPC-based 3-D trajectory tracking for an autonomous underwater vehicle with constraints in complex ocean environments. Ocean Eng 2019, 189, 106309. [Google Scholar] [CrossRef]
Bi, F.; Wei, Y.; Zhang Z, J.; Cao, W. Position-tracking control of underactuated autonomous underwater vehicles in the presence of unknown ocean currents. Iet Control Theory and Applications 2010, 4, 2369–2380. [Google Scholar] [CrossRef]
Yuh, J. Design and Control of Autonomous Underwater Robots: A Survey. Autonomous Robots 2000, 8(1), 7–24. [Google Scholar] [CrossRef]
Haj-Ali, A.; Ying, H. Structural analysis of fuzzy controllers with nonlinear input fuzzy sets in relation to nonlinear PID control with variable gains. Automatica 2004, 40(9), 1551–1559. [Google Scholar] [CrossRef]
Krstic, M. Feedback Linearizability and Explicit Integrator Forwarding Controllers for Classes of Feedforward Systems. IEEE Transactions on Automatic Control 2004, 49(10), 1668–1682. [Google Scholar] [CrossRef]
Tsai, WC.; Wu, CH.; Cheng, MY. Tracking Accuracy Improvement Based on Adaptive Nonlinear Sliding Mode Control. IEEE/ASME Transactions on Mechatronics 2021; 26(1): 179–190.
Zhu, C.; Huang, B.; Zhou, B.; Su, Y.; Zhang, E. Adaptive model-parameter-free fault-tolerant trajectory tracking control for autonomous underwater vehicles. ISA Transactions 2021, 114, 57–71. [Google Scholar] [CrossRef]
Wadi, A.; Mukhopadhyay, S.; Lee, JH. Adaptive observer design for wave PDEs with nonlinear dynamics and parameter uncertainty. Automatica 2021, 123, 109295. [Google Scholar]
Wen, L.; Tao, G.; Song, G. Higher-order tracking properties of nonlinear adaptive control systems. Systems & Control Letters 2020, 145, 104781. [Google Scholar]
Zhang, X.; Baron, L.; Liu, Q.; Boukas, EK. Design of Stabilizing Controllers With a Dynamic Gain for Feedforward Nonlinear Time-Delay Systems. IEEE Transactions on Automatic Control 2011, 56(3), 692–697. [Google Scholar] [CrossRef]
Elmokadem, T.; Zribi, M.; Youcef-Toumi, K. Terminal sliding mode control for the trajectory tracking of underactuated Autonomous Underwater Vehicles. Ocean Eng 2017, 129, 613–625. [Google Scholar] [CrossRef]
Wang, Y.; Pu, H.; Shi, P.; Ahn, CK.; Luo, J. Sliding mode control for singularly perturbed Markov jump descriptor systems with nonlinear perturbation. Automatica 2021, 127, 109515. [Google Scholar] [CrossRef]
Taleb, M.; Plestan, F.; Bououlid, B. An adaptive solution for robust control based on integral high-order sliding mode concept. International Journal of Robust and Nonlinear Control 2015, 25(8), 1201–1213. [Google Scholar] [CrossRef]
Liu, K.; Liu, K.; Wang, Y.; Ji, H.; Wang, S. Adaptive saturated tracking control for spacecraft proximity operations via integral terminal sliding mode technique. International Journal of Robust and Nonlinear Control 2021, 31(18), 9372–9396. [Google Scholar] [CrossRef]
Khan, S.; Guivant, J.; Li, X. Design and experimental validation of a robust model predictive control for the optimal trajectory tracking of a small-scale autonomous bulldozer. Robotics and Autonomous Systems 2022, 147, 103903. [Google Scholar] [CrossRef]
Hyun-Wook. ; Jo, Jong-Tae.; Lim.; Ho-Lim.; Choi. Observer based output feedback regulation of a class of feedforward nonlinear systems with uncertain input and state delays using adaptive gain. Systems and Control Letters 2014, 71, 45–53. [Google Scholar]
Chen, X.; Zhang, X.; Liu, Q. Prescribed-time decentralized regulation of uncertain nonlinear multi-agent systems via output feedback. Systems & Control Letters 2020, 137, 104640. [Google Scholar]
Krishnamurthy, P.; Khorrami, F. Feedforward Systems with ISS Appended Dynamics: Adaptive Output-Feedback Stabilization and Disturbance Attenuation. IEEE Transactions on Automatic Control 2008, 53(1), 405–412. [Google Scholar] [CrossRef]
Machado, JE.; Ortega, R.; Astolfi, A. ; Jos. An Adaptive Observer-Based Controller Design for Active Damping of a DC Network with a Constant Power Load. IEEE Transactions on Control Systems Technology 2021, 29, 2312–2324. [Google Scholar] [CrossRef]
Chaudhari, S.; Shendge, PD.; Phadke, SB. Disturbance Observer Based Controller Under Noisy Measurement for Tracking of n DOF Uncertain Mismatched Nonlinear Interconnected Systems. IEEE/ASME Transactions on Mechatronics 2020, 25(3), 1600–1611. [Google Scholar] [CrossRef]
Liu, S.; Wang, D.; Poh, E. Output feedback control design for station keeping of AUVs under shallow water wave disturbances. International Journal of Robust & Nonlinear Control 2010, 19(13), 1447–1470. [Google Scholar]
Negahdaripour, S.; Xu, X.; Khamene, A.; Awan, Z. 3-D motion and depth estimation from sea-floor images for mosaic-based station-keeping and navigation of ROVs/AUVs and high-resolution sea-floor mapping. IEEE Proceedings of the 1998 Workshop on Autonomous Underwater Vehicles 1998, PP:191-200. [Google Scholar]
Zhuan, X.; Xia, X. Optimal Scheduling and Control of Heavy Haul Trains Equipped With Electronically Controlled Pneumatic Braking Systems. IEEE Transactions on Control Systems Technology 2007, 15(6), 1159–1166. [Google Scholar] [CrossRef]
Du, T.; Hughes, J.; Wah, S.; Matusik, W.; Rus, D. Underwater Soft Robot Modeling and Control with Differentiable Simulation. IEEE Robotics and Automation Letters 2021, PP(99), 1–1. [Google Scholar] [CrossRef]
Chen, YL.; Ma, XW.; Bai, GQ.; Sha, Y.; Liu, J. Multi-autonomous underwater vehicle formation control and cluster search using a fusion control strategy at complex underwater environment. Ocean Eng 2020, 216(7), 108048. [Google Scholar] [CrossRef]
Hu, C.; Fu, L.; Yang, Y. Cooperative navigation and control for surface-underwater autonomous marine vehicles. IEEE 2018. [Google Scholar]
Liang, Z.; Liu, Q. Design of stabilizing controllers of upper triangular nonlinear time-delay systems. Systems & Control Letters 2015, 75, 1–7. [Google Scholar]
Li, DJ. Adaptive output feedback control of uncertain nonlinear chaotic systems based on dynamic surface control technique. Nonlinear Dynamics 2012, 68(1-2), 235–243. [Google Scholar] [CrossRef]
Huang, Y.; Meng, Z. Global finite-time distributed attitude synchronization and tracking control of multiple rigid bodies without velocity measurements. Automatica 2021, 132, 109796. [Google Scholar] [CrossRef]

Figure 1. Hardware-in-the-loop simulation device with underwater robot.

Figure 2. The trajectory of Output.

Figure 3. The trajectory of Station

l_{1}

Figure 3. The trajectory of Station

l_{1}

Figure 4. The trajectory of Station

l_{2}

Figure 4. The trajectory of Station

l_{2}

Figure 5. The trajectory of Observers

{\hat{l}}_{1}

and

{\hat{l}}_{2}

Figure 5. The trajectory of Observers

{\hat{l}}_{1}

and

{\hat{l}}_{2}

Figure 6. The trajectory of Errors

l_{1}

and

l_{2}

Figure 6. The trajectory of Errors

l_{1}

and

l_{2}

Figure 7. The trajectory of Controller

u

Figure 7. The trajectory of Controller

u

Figure 8. The trajectory of

D

and

\dot{D}

Figure 8. The trajectory of

D

and

\dot{D}

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.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.