Download PDF
Research Article  |  Open Access  |  13 Dec 2024

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Views: 43 |  Downloads: 4 |  Cited:  0
Intell Robot 2024;4(4):439-56.
10.20517/ir.2024.26 |  © The Author(s) 2024.
Author Information
Article Notes
Cite This Article

Abstract

This paper presents a novel trajectory tracking controller for an underactuated unmanned surface vehicle (USV). The controller incorporates an event-triggered extended state observer (ETESO), a minimum learning parameter neural network, an integral non-singular terminal sliding mode (INTSM) control strategy, and a dynamic event-triggered mechanism (DETM). Firstly, an ETESO is developed to estimate unmeasurable velocities and lumped disturbances, differentiating it from most existing extended state observers without the necessity for real-time output measurements. To further alleviate the communication burden and minimize actuator wear, a DETM with an adjustable threshold is introduced. In contrast to traditional event-triggered methods, which employ fixed threshold parameters, this mechanism allows for online adaptive updates of the triggering thresholds, thereby enhancing resource efficiency. Additionally, an INTSM is designed to ensure rapid convergence of the position and velocity errors of the USV. To effectively counteract external disturbances and internal modeling uncertainties, a minimum learning parameter (MLP) neural network algorithm is implemented to approximate and compensate for these uncertainties. Finally, using Lyapunov's theory, it is demonstrated that all signals within the closed-loop tracking control system remain bounded. Simulation results are given to illustrate the effectiveness of theoretical results.

Keywords

Event-triggered extended state observer, dynamic event-triggered mechanism, USV, integral non-singular terminal sliding mode control

1. INTRODUCTION

The ocean is the largest and most complex ecosystem on Earth, covering approximately 71% of the planet's surface and harboring rich biodiversity and resources. In recent years, the rapid growth of the global population has led to a sharp increase in energy consumption, which in turn has intensified environmental pollution, posing a continuous threat to the living conditions of humans on land[1]. Ocean exploration is not only an important area of scientific research but also a key factor in achieving sustainable development for humanity. At the same time, with increasing attention on maritime rights by various countries, the ocean has become a focal point of international competition[2]. Unmanned surface vehicles (USVs) have garnered significant interest from the control science community[3]. This is mainly due to their advantages in a simple control structure and low manufacturing cost[4]. In terms of application, USVs have been employed in several contexts, including seabed resource extraction, bathymetry, marine search and rescue, oceanographic surveys, and sampling and patrolling[5]. However, the complexity of the marine environment, the unmodeled dynamics in the mathematical models of the USV, and the inability to measure the speed of USV in real-world ocean conditions can lead to a degradation of the maneuvering performance of the USV, potentially causing instability in system performance. Moreover, the real-time transmission of command signals from the controller to the actuators can lead to unnecessary communication burdens and excessive wear on the actuators[6]. Given these challenges, further research into USVs and their integration into offshore engineering is crucial.

Typically, the data exchange of the USV is facilitated through wireless communication networks. This involves connecting system components, such as sensors, controllers, filters, and actuators, through the communication network, with controllers typically situated in remote motherships or land-based stations[7]. Therefore, it is paramount to reduce the computational cost and conserve communication resources to achieve adequate control of the USV. In light of the circumstances above, event-triggered control has garnered increasing attention in recent years, as evidenced by the work of[8]. Jiang et al. presented a method for addressing the adverse effects of denial of service (DoS) attacks by incorporating the phenomenon of signal non-transmission into the event-triggered interval[9]. Zhou et al. proposed an event-triggered control mechanism scheme based on dynamic surface control and adaptive dynamic programming[10]. This approach ensures that tracking the reference trajectory effectively saves computation and reduces the number of controller executions. However, the event-triggered thresholds in the above studies are static, which imposes rigid communication at the expense of performance. In order to further reduce the computational cost and save communication resources, Girard and Antoine (2014) presented a dynamic triggering mechanism for event-triggered control[11]. Subsequently, dynamic event-triggered mechanisms (DETM) have been widely applied in the field of control. For example, Cao et al. combined network-induced errors and relative threshold strategies to establish two new DETM and dynamic rules for threshold parameters, reducing communication between the controller and the actuator[12]. He et al. proposed a distributed dynamic event-triggered strategy, introducing an auxiliary parameter for each agent to dynamically adjust its threshold[13]. Wang and Chang proposed a new DETM[14]. By analyzing data in real time, the event-triggering parameters are adjusted to ensure better triggering performance. Compared with the traditional static event-triggered control, dynamic event triggering has a longer time interval. However, there is currently no trajectory tracking control algorithm for USVs based on DETM. Therefore, researching trajectory tracking control for USVs based on DETM is of significant importance.

The theory of sliding mode control was mainly developed in the early 1990s[15]. In recent years, sliding mode variable structure control has been widely used in the USV control field because of its strong robustness, fast response and simple implementation[16,17]. In order to reduce the jitter of the system and improve the exemplary reading of the control system, Liu et al. utilized super-twisting sliding mode control to design the system[18]. Yu et al. employed an integral form of the sliding mode surface for the control of the USV[19]. This approach effectively reduces the system's steady-state error and enhances its steady-state performance. In addition, Zhang et al. used a terminal sliding mode surface to achieve the convergence of the control error of the USV system in a finite time[20]. Although this terminal sliding surface can make the system state converge in finite time, it may make the system exist at a singular point of zero. Based on those mentioned above, this paper employs an integral form of non-singular terminal sliding mode surface for trajectory tracking control of the USV; That is, it can make the system error converge in finite time without a singular point and eliminate the system's steady-state error.

Many methods have been used to estimate disturbances and interferences in the control of USVs. Sun et al. studied the optimal coverage control problem for multiple USVs under time-varying disturbances[21]. A disturbance vector observer was designed to approximate the unknown time-varying disturbances to address this issue. In the study of Han et al., a nonlinear version of the Kalman filter-based active modeling method was proposed to provide online estimates of the unstructured model to eliminate the errors due to the structural inaccuracies between the quasi-linear parameter-varying (qLPV)-structured model and the natural system[22]. Chen et al. proposed a disturbance-observer-based sliding mode control design to achieve good tracking performance, where the observer estimated and compensated for the modeling uncertainties and external disturbance[23]. In the study of Chen et al., an adaptive sliding mode control design for nonlinear USVs is proposed, incorporating a radial basis function neural network to approximate system modeling uncertainty and a disturbance observer to estimate external unknown disturbances[24]. However, radial basis function neural network algorithms require online adjustment of the ownership vector of the network, which increases the computational effort. Therefore, the minimum learning parameter (MLP) neural network algorithm can solve the problem effectively. Most of the above studies assume that the speed of unmanned vessels is measurable and controllers are directly designed based on this assumption. However, the actual speed of marine vessels during navigation is not measurable. Therefore, it is necessary to estimate the state variables using available information (inputs and outputs) through models. Using observer-based feedback control is a very important strategy. Typically, an observer is added to establish a simulation system with the same dynamic equations as the actual system, thereby enabling state estimation, system performance analysis, and control design[25].

Based on the information provided, this paper puts forward an integral non-singular terminal sliding mode (INTSM) trajectory tracking control algorithm. An event-triggered extended state observer (ETESO) is developed to estimate the velocity and utilizes an MLP neural network algorithm to estimate the internal and external disturbances of the USV. Furthermore, dynamic event triggers are incorporated into the control channel to conserve communication resources while tracking a specified trajectory. The primary contributions of this paper compared to existing research results are as follows:

(1) firstly, this paper introduces a DETM with adjustable thresholds into the USV's control input. In contrast to the existing event-triggered mechanisms widely used in USV trajectory tracking control[26,27], which have fixed threshold parameters for triggering results, the mechanism employed in this study can adaptively update the triggering thresholds online, thereby conserving communication resources more effectively. Additionally, to further save communication resources, an event-triggered mechanism has also been incorporated into the extended state observer;

(2) compared to the traditional sliding mode surfaces used for trajectory tracking control of the USV in ref.[28,29], the INTSM surface employed in this paper not only ensures that the tracking error of the USV converges within a finite time but also reduces the steady-state error of the system;

(3) in relation to ref.[30,31], which uses radial basis function neural networks to directly approximate USV's modeling uncertainties and external disturbances, this paper further employs the MLP technique to compress the weights of the neural network. Replacing the online learning of all weight vectors with single-parameter online learning reduces the amount of computation required to estimate and compensate for external disturbances.

2. PROBLEM DESCRIPTION AND PRELIMINARY PREPARATION

2.1. Model of USV

The mathematical modeling of USVs is essential for analyzing their motion dynamics. In recent studies addressing trajectory tracking control for underactuated vessels, the Fossen model is predominantly employed, where the longitudinal thrust and steering torque of the USV serve as control outputs to regulate its movement[32]. Following this, the kinematic and dynamic models of the USV can be represented mathematically as follows:

$$ \dot{\eta}=R\left(\psi\right)\nu, M\ddot{\eta}+C_{0}\dot{\eta}+D_{0}\dot{\eta}=R\tau+\Delta f+\tau_{w} . $$

The variables are defined as follows: $$ \eta=[x, y, \psi]^\mathrm{T} $$ represents the surge, sway, and yaw angle of the ship in the inertial frame; $$ \nu =[u, v, r]^\mathrm{T} $$ represents the velocity vector of the ship in the body-fixed frame; $$ R(\psi) $$ is the transformation matrix between the inertial frame and the body-fixed frame, and it satisfies $$ R^{- 1}( \psi ) = R^{\mathrm{T} }( \psi ) $$; $$ M $$ is the inertia matrix of the ship and satisfies $$ M=M^T>0 $$; $$ C(v) $$ is the Coriolis and centripetal matrix and satisfies $$ C(v)=-C(v)^T $$; $$ D $$ is the damping matrix; $$ \tau=[\tau_u, 0, \tau_r]^{\mathrm{T}} $$ is the control force and moment vector for the underactuated ship, where $$ \tau_u $$ and $$ \tau_r $$ represent the longitudinal thrust and the steering moment, respectively. For underactuated ships, there is no lateral thruster, so $$ \tau_v=0 $$; $$ \tau_w=[\tau_{wu}, 0, \tau_{wr}]^\mathrm{T} $$ denotes the unknown environmental time-varying interference in the body-fixed frame, where $$ \tau_{wu} $$ is the environmental interference forces of the ship in the longitudinal directions, and $$ \tau_{wr} $$ is the interference force of the ship in the yaw direction. The expressions for $$ R(\psi) $$, $$ M $$, $$ C(v) $$, and $$ D $$ are given by:

$$ \begin{array}{c} R(\psi)=\begin{bmatrix}\cos\psi&-\sin\psi&0\\\sin\psi&\cos\psi&0\\0&0&1\end{bmatrix}, M=\begin{bmatrix}m_{11}&0&0\\0&m_{22}&0\\0&0&m_{33}\end{bmatrix}\\C(v)=\begin{bmatrix}0&0&-m_{22}v\\0&0&m_{11}u\\m_{22}v&-m_{11}u&0\end{bmatrix}, \quad\boldsymbol{D}=\begin{bmatrix}d_{11}&0&0\\0&d_{22}&0\\0&0&d_{33}\end{bmatrix} \end{array} $$

Moreover, the parameters are defined as follows: $$ m_{11}=m-X_{\dot{u}}, m_{22}=m-Y_{\dot{v}} $$, $$ m_{33}=I_{z}-N_{\dot{t}} $$, $$ d_{11}=-X_{u} $$, $$ d_{22}=-Y_{v} $$, $$ d_{33}=-N_{\tau} $$. $$ m $$ represents the mass of the underactuated ship, $$ I_\mathrm{z} $$ is the moment of inertia, and other parameters such as $$ X_\mathrm{i}=\partial X/\partial\dot{u} $$ are hydrodynamic derivatives. These symbols were standardized by the Society of Naval Architects and Marine Engineers (SNAME) in 1950, and will not be elaborated here. The mathematical model of the USV can also be written in the following form:

$$ \left\{ \begin{array}{l} \dot{x}=u\cos\left(\psi\right)-v\sin\left(\psi\right)\\ \dot{y}=u\sin\left(\psi\right)+v\cos\left(\psi\right)\\ \dot{\psi}=r\\ \dot{u}=\frac{1}{m_{11}}\big(m_{22}v r-d_{11}u+\Delta f_{u}+\tau_{u}+\tau_{wu}\big)\\ \dot{v}=\frac{1}{m_{22}}\big(-m_{11}ur-d_{22}v+\Delta f_{v}+\tau_{wv}\big)\\ \dot{r}=\frac{1}{m_{33}}\Big(\big(m_{11}-m_{22}\big)uv-d_{33}r+\Delta f_{r}+\tau_{r}+\tau_{wr}\Big) \end{array} \right. $$

2.2. Design of ETESO

In this section, ETESO is designed to obtain accurate estimates of velocities with fewer communication resources. According to the USV mathematical model in Equation (1), the composite dynamics model of the system is given as:

$$ \ddot{\boldsymbol{\eta}}={\mathbf R}(\boldsymbol{\psi})\dot{\boldsymbol{\upsilon}}+\dot{{\mathbf R}}(\boldsymbol{\psi}){\mathbf \nu}={\mathbf R}(\boldsymbol{\psi}){\mathbf M}^{-1}[-{\mathbf C}({\mathbf \nu})-{\mathbf D}\boldsymbol{v}-\boldsymbol{\tau}_w-\Delta{\mathbf f}]+\dot{{\mathbf R}}(\boldsymbol{\psi}){\mathbf v}+{\mathbf R}(\boldsymbol{\psi}){\mathbf M}^{-1}\boldsymbol{\tau}=\boldsymbol{\chi}+{\mathbf R}(\boldsymbol{\psi}){\mathbf M}^{-1}\boldsymbol{\tau}. $$

We define $$ {\mathbf x}_{1}={\mathbf \eta} $$, $$ {\mathbf x}_2=\dot{\eta} $$. Consequently, the dynamics model of the USV can be expressed as:

$$ \begin{cases}\dot{{\mathbf x}}_1={\mathbf x}_2\\\dot{{\mathbf x}}_2={\mathbf \chi}+{\mathbf R}(\psi){\mathbf M}^{-1}{\mathbf \tau}\end{cases}, $$

where $$ \chi={\mathbf R}(\psi){\mathbf M}^{-1}[-{\mathbf C}({\mathbf v})-{\mathbf D}v-{\mathbf \tau}_w-\Delta{\mathbf f}]+\dot{{\mathbf R}}(\psi)\nu $$ represents the lumped uncertainty, which includes unmodeled dynamics and external disturbances.

To design ETESO, we introduce Assumption 1, which is often used in the control of the USV[33].

Assumption 1[34] The composite disturbance vector $$ \tau_{w}=[\tau_{wu}, 0, \tau_{wr}] $$ is time-varying and unknown, with its Euclidean norm bounded by a constant, i.e., $$ \parallel \tau_w\parallel < \bar{\tau}_w $$, where $$ \bar{\tau}_w $$ is a constant. Additionally, the lumped uncertainty $$ \chi $$ satisfies $$ \parallel\dot{\chi}\parallel\leq H $$, where $$ H > 0 $$ is an unknown bounded constant.

The following ETESO is designed

$$ \begin{cases}\dot{\hat{{\mathbf x}}}_1=\hat{{\mathbf x}}_2+{\mathbf \varepsilon}_e g_1\left ( \frac{{\mathbf{\xi(t)}}-\hat{{\mathbf x}}_1(t)}{{\mathbf \varepsilon}_e^2}\right ) \\ \dot{\hat{ {\mathbf x}}}_2=\hat{{\mathbf x}}_3+g_2\left ( \frac{{\mathbf{\xi(t)}}-\hat{{\mathbf x}}_1(t)}{{\mathbf \varepsilon}_e^2}\right)+R(\psi)M^{-1}\tau \\ \dot{{\mathbf \chi}}= {\mathbf \varepsilon}_e^{-1} g_3\left ( \frac{{\mathbf{\xi(t)}}-\hat{{\mathbf x}}_1(t)}{{\mathbf \varepsilon}_e^2}\right ) \end{cases}, $$

where $$ \widehat{{\mathbf x}}_1 $$, $$ \widehat{{\mathbf x}}_2 $$ and $$ {\mathbf \chi} $$ are estimate of $$ {\mathbf x}_1 $$, $$ {\mathbf x}_1 $$ and $$ {\mathbf \chi} $$, respectively. The parameter $$ {\mathbf \varepsilon}_e $$ is a high-gain parameter to be determined. The functions $$ g_{i}(*), i=1, 2, 3 $$ are defined as follows: $$ g_{1}(\vartheta)=a_{1}\vartheta+\varphi(\vartheta) $$, $$ g_{j}(\vartheta)=a_{j}\vartheta, j=2, 3 $$, are constants to be chosen to ensure that the matrix $$ N=\begin{bmatrix}-a_{1}&1&0\\-a_{2}&0&1&\\-a_{3}&0&0\end{bmatrix} $$ is Hurwitz, and $$ \varphi(\vartheta) $$ is defined as $$ \left.\varphi(\vartheta)=\left\{\begin{array}{ll}\frac{1}{4}sign(\vartheta), &if |\vartheta|>\frac{\pi}{2}\\\frac{1}{4}sin(\vartheta), &if |\vartheta|\leq\frac{\pi}{2}\end{array}\right.\right. $$. where $$ \zeta(t) $$ is the measurement output, which is defined as:

$$ \zeta(t)=\begin{cases}\eta(t_{o, k}), if \Gamma_1\leq0, \\ \eta(t), \quad otherwise, \end{cases} $$

with an event-triggered condition given by $$ \Gamma_1=\sum_{i=1}^{3}L_i|\sigma(t)|-\overline{B}\varepsilon_e $$, where $$ \overline{B} $$ is a positive constant to be determined, and $$ L_i>0 $$ is a Lipschitz constant for $$ g_{i}(*) $$. The $$ \sigma(t) $$ is an sampled error defined as $$ \sigma(t)=\frac{\eta(t_{o, k})-\eta(t)}{\varepsilon^{2}_e} $$, for $$ t\in[t_{o, k}, t_{o, k+1}) $$. Here, $$ t_{o, k} $$ denotes the triggered instant, which depends on condition 6.

Lemma 1[35] For any $$ \ell > 0 $$ and $$ a\in\mathbb{R} $$, the following inequality holds:

$$ 0\le\mid a\mid-atanh(\frac{a}{\ell})\le\aleph a . $$

Given that $$ \aleph =e^{-(1+\aleph)} $$, we take $$ \aleph=0.2785 $$.

Lemma 2[36] For a given unknown continuous function $$ f({\mathbf x}){:}\Re^m\to\Re $$, it can be approximated over a compact set $$ \Omega_{\mathrm{x}}\subset\Re^m $$ with the following radial basis function neural network:

$$ f=W^{*T}h(z)+\varepsilon $$

where $$ W^{*T}=[W_1^{*}, W_1^{*}, …, W_k^{*}]^T $$ is the ideal weight vector with $$ k $$ being the number of NN nodes and matrix $$ z=\left[u, v, r\right]^{\tau} $$ represents the input signals of radial basis function neural network. $$ \varepsilon $$ represents the estimated error of neural network with the upper bound $$ |\varepsilon|\le\varepsilon^* $$, where $$ \varepsilon^* $$ is an unknown positive constant. The basis function $$ h(x)=\left[h_{1}(x), h_{2}(x), ..., h_{n}(x)\right]^{\mathrm{r}} $$ is chosen as the Gaussian function, that is

$$ h_{j}\left(x\right)=\exp\Biggl[-\frac{\left\|x-c_{j}\right\|^{2}}{2{b_{j}}^{2}}\Biggr], \bigl(j=1, \cdots, n\bigr) $$

where $$ b_{j}>0 $$ and $$ c_{j}=\left[c_{1}, c_{2}, \cdots, c_{m}\right]^{\mathrm{T}}\in{\mathbf R}^{m} $$ are the center of the receptive field and the width of the Gaussian function. Formally, the optimal weight vector $$ W^* $$ is computed as:

$$ W^*=arg\min\limits_{\hat{\boldsymbol{\gamma}}\in\mathbb{R}^k}\left\{sup\left|f(\boldsymbol{x})-\hat{\boldsymbol{W}}^T\boldsymbol{S}(\boldsymbol{x})\right|\right\} $$

where $$ \hat{W} $$ is the estimate of $$ W^* $$.

Define $$ \tilde{W}=W^*-\hat{W} $$, then

$$ \tilde{f}=f-\hat{f}=W^{*T}h(z)+\varepsilon-\hat{W}^Th(z)=\tilde{W}^Th(z)+\varepsilon $$

Lemma 3[37] Consider a nonlinearsystem $$ \dot{{\mathbf x}}=f({\mathbf x}(t)) $$, $$ {\mathbf x} \in \Re $$, if there exists a continuous and positive definite Lyapunov function $$ V({\mathbf x}) $$ satisfying $$ k_1(\|{\mathbf x}\|) \leq V({\mathbf x}) \leq k_2(\|{\mathbf x}\|) $$ such that $$ \dot{V} \leq-l V+ $$$$ Q $$, where $$ k_1, k_2 $$ are class $$ \kappa $$ functions and $$ l $$, $$ Q $$ are positive constants, and then the solution $$ {\mathbf x}(t) $$ is ultimately uniformly bounded.

3. TRAJECTORY TRACKING CONTROLLER DESIGN

In this section, we will design an INTSM controller to effectively track the trajectory of the USV. The complete design procedure will be elaborated upon in the following subsections.

3.1. Virtual control design

For future use, we will define the following tracking errors:

$$ \left. \left[\begin{array}{c}x_{e}\\y_{e}\\\psi_{e}\end{array}\right.\right]=J^{\mathrm{T}}(\psi)\left[\begin{array}{c}x_{d}-x\\y_{d}-y\\\psi_{d}-\psi\end{array}\right], $$

where $$ (x, y) $$ denotes the position of the USV in the earth - fixed frame, , and $$ \psi $$ represents the vessel's yaw angle. The terms $$ x_\mathrm{d}, y_\mathrm{d}, \psi_{e} $$ denote the longitudinal position error, transverse position error, and heading angle error, respectively. Additionally, $$ z_e=\sqrt{x_{e}^{2}+y_{e}^{2}} $$, where $$ x_\mathrm{e} $$ and $$ y_\mathrm{e} $$ represent the vertical and horizontal position errors. These errors correspond to the horizontal and vertical coordinates of the reference trajectory. The expected course can be given as:

$$ \left.\psi_{d}=\begin{cases}&\arctan\left(\frac{\dot{y}_{d}}{\dot{x}_{d}}\right), z_{e}=0\\&0.5\left[1-\operatorname{sgn}\left(x_{e}\right)\right]\operatorname{sgn}\left(y_{e}\right)\pi+\arctan\left(\frac{y_{e}}{x_{e}}\right), \operatorname{other}\end{cases}\right.. $$

In Figure 1, the coordinate system $$ OXY $$ represents the inertial reference frame, with point $$ O $$ as the initial position. The $$ OX $$ axis points due north, while the $$ OY $$ is due east. Point $$ Ob $$ is defined as the midpoint between the bow and stern of the USV. The vector $$ o_bx_b $$ extends toward the bow along the USV's midline, and the vector $$ obyb $$ is directed along the port side of the USV. Additionally, line segment $$ AB $$ represents the reference trajectory.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 1. Basic schematic diagram of USV trajectory tracking. USV: Unmanned surface vehicle.

According to the relationships depicted in Figure 1, the variables $$ x_e $$, $$ y_e $$ and $$ z_e $$ can be expressed as:

$$ x_{e}=z_{e}\cos\psi_{d}\quad y_{e}=z_{e}\sin\psi_{d}. $$

Taking the time derivative of the variables $$ z_e $$ and $$ \psi_e $$, along with utilizing Equations (2) and (14), we can derive

$$ \begin{array}{l} \dot{z}_{e}=\dot{x}_{d}\cos\psi_{d}+\dot{y}_{d}\sin\psi_{d}-u\cos\psi_{e}-v\sin\psi_{e}\\ \dot{\psi}_{e}=\dot{\psi}_{d}-r . \end{array} $$

Design the virtual control functions $$ \alpha_{u} $$ and $$ \alpha_{r} $$ as follows:

$$ \alpha_{u}= \cos^{-1}\left(\psi_{e}\right)\left [k_{ze1}z_{e}+k_{ze2}\mathrm{sgn}^{\delta}\left(z_{e}\right)\right] + \dot{x}_{d}\cos\left(\psi_{d}\right)+\dot{y}_{d}\sin\left(\psi_{d}\right)-\upsilon\sin\left(\psi_{e}\right), $$

$$ \alpha_{r}=k_{\psi e1}\psi_{e}+k_{\psi e2}\mathrm{sgn}^{\delta}\left(\psi_{e}\right)+\dot{\psi}_{d}, $$

where $$ k_{ze1}>0 $$, $$ k_{ze2}>0 $$, $$ k_{\psi e1}>0 $$, $$ k_{\psi e2}>0 $$, $$ 0.5\leq\delta<1 $$ are the design parameters.

Note that when $$ \psi_{e} = \pm\frac{\pi}{2} $$, $$ \alpha_{u} $$ becomes undefined. Therefore, in actual engineering, the first assumption is that condition $$ |\psi_{e}| < \frac{\pi}{2} $$ holds, and the assumption is guaranteed by the transformation of[38]

$$ \left.\psi_e=\left\{\begin{array}{lc}\psi_e-\pi, &\quad\psi_e\geq0.5\pi\\\psi_e, &\quad-0.5\pi<\psi_e<0.5\pi\\\psi_e+\pi, &\quad\psi_e\leq-0.5\pi\end{array}\right.\right. $$

3.2. Design of INTSM controller

To prevent dimensional explosion, we introduce the following first-order filter.

$$ \mu\dot{\gamma}+\gamma=\alpha\quad\gamma(0)=\alpha(0), $$

where $$ \gamma $$ is first-order filter and $$ \mu $$ denotes the time constant.

Define the longitudinal error and the yaw velocity error as:

$$ e_{u}=u-\gamma_{u}, e_{r}=r-\gamma_{r}. $$

Introduction of an INTSM surface for longitudinal error and yaw velocity error:

$$ \begin{eqnarray} S_i=\lambda_i\int_{0}^{t} sig^{r_{i}}(e_i)d\tau+\beta_i e_i, i=u, r \end{eqnarray} $$

given that $$ \lambda_i $$, $$ \beta_i $$ and $$ r_i>1 $$ are designed constants, we proceed to differentiate $$ S_i $$ to obtain

$$ \dot{S}_i = \lambda_isig^{r_i}(e_i)+\beta_{i}\dot{e}_i, i=u, r . $$

The surge motion control law $$ \tau_u $$ and yaw motion control law $$ \tau_r $$ can be designed as:

$$ \tau_u=m_{11}\dot{\gamma}_u-\hat{\tau}^{*}_{wu}tanh(S_u/\hat{\varepsilon}_{u})-\bigtriangleup\hat{f}_u+d_{11}\hat{u}-m_{22}\hat{v}\hat{r}-\frac{m_{11}\lambda_usig^{r_u}e_u}{\beta_u}-k_utanh(S_u), $$

$$ \tau_r=m_{33}\dot{\gamma}_r-\hat{\tau}^{*}_{wr}tanh(S_r/\varepsilon_{r})-\bigtriangleup\hat{f}_r+d_{33}\hat{r}-(m_{11}-m_{22})\hat{u}\hat{v}-\frac{m_{33}\lambda_rsig^{r_r}e_r}{\beta_r}-k_rtanh(S_r) , $$

where $$ \Delta\widehat{f}_u=\frac12s_u\widehat{\phi}_uhh^T $$ and $$ \Delta\widehat{f}_r=\frac12s_r\widehat{\phi}_rhh^T $$ represent the estimate of $$ \bigtriangleup f_u $$ and $$ \bigtriangleup f_r $$, respectively. The corresponding adaptive laws are given in

$$ \dot{\hat{\phi}}_i=\hat{\lambda}_i\left(\frac12s_i^2h^\mathrm{T}h-\kappa_u\hat{\phi}_i\right), i=u, v $$

The adaptive law for model uncertainty in the USV is defined as follows:

$$ \dot{\widehat{\tau}}_{Wi}^{*}=\gamma_i\left[S_i\psi(S_i)-\mu_i\left(\widehat{\tau}_{wi}^{*}-\tau_{wi}^0\right)\right] , i=u, v $$

3.3. Dynamic event-triggered input scheme with ETESO

In this subsection, we present a dynamic event-triggered control scheme in the last step. First, we define $$ \hat{\tau}=\begin{bmatrix} \hat{\tau}_u&\hat{\tau}_r \end{bmatrix}^T $$, the dynamic event-triggered control law $$ \hat{\tau} $$ is given as:

$$ \hat{\tau}(t) = \tau\left(t_{j}\right), \forall t\in[t_{j}, t_{j+1}) , j\in\mathbb{N}. $$

It is evident that $$ \hat{\tau}(t) $$ is updated only at triggering instants $$ t_j $$, where $$ j\in\mathbb{N} $$. We define the measurement error $$ z(t) $$ as

$$ z\left(t\right) = \hat{\tau}\left(t\right)-\tau\left(t\right), \forall t\in\left[t_{j}, t_{j+1}\right), j\in\mathbb{N}. $$

The DETM is designed as:

$$ t_{j+1}=\inf\limits_{t>t_j}\{t\in\mathbb{R}_{\geq0}|\eta_i(t)+\left(k_i-1\right)S_i^{2}-Z^2_i\leq0\}, t_{0}=0, $$

where $$ \eta\left(t\right) $$ is an internal dynamic variable satisfying

$$ \dot{\eta}=-\chi\eta+\left(k_i-1\right)S_i^{2}-Z_i^{2}, \eta(0)=\eta_0. $$

It is evident that the inequality $$ -\eta_k(t)\leq\left(k_i-1\right)S_i^{2}-Z^2_i $$ holds all the time. We will now demonstrate that $$ \eta_k(t)>0 $$ is indeed true. Based on the definition of $$ Z(t) $$ in Equation (28), we can conclude that $$ ||Z(t)|| \leq ||Z(t^-)|| $$ for all $$ t\ge 0 $$. Then, the DETM (29) ensures that for all $$ t>0 $$

$$ \eta_{i}(t)+(k_{i}-1) S_{i}^{2}-Z_{i}^{2}\ge0 $$

from Equation (31), we have

$$ (k_{i}-1) S_{i}^{2}-Z_{i}^{2}\ge-\eta_{i}(t) $$

Combining Equations (30) and (32), it can be deduced that for all $$ t \geq 0 $$

$$ \dot{\eta}(t)\geq-\chi\eta(t)-\eta(t), \eta(0)=\eta_0>0. $$

Then by the Comparison Lemma, we can obtain that for all $$ t \geq 0 $$

$$ \eta(t)\geq\eta_0e^{-(\chi+1)t}\geq0. $$

4. STABILITY ANALYSIS

To prove the stability of ETESO, we introduce Assumption 2, which is often used to ensure the stability of the nominal system considered[39].

Assumption 2[40] There exist two nonnegative definite functions $$ P(*) $$ and $$ W(*):\mathbb{R}^{n+1}\to\mathbb{R} $$ satisfying

$$ \lambda_{1}||\vartheta||^{2}\leq P(\vartheta)\leq\lambda_{2}||\vartheta||^{2}, \lambda_{3}||\vartheta||^{2}\leq W(\vartheta)\leq\lambda_{4}||\vartheta||^{2}, $$

$$ \sum\limits_{j=1}^{n}\frac{\partial P}{\partial\vartheta_{j}}(\vartheta_{j+1}-g_{j}(\vartheta_{1}))-\frac{\partial P}{\partial\vartheta_{n+1}}g_{n+1}(\vartheta_{1})\leq-W(\vartheta), $$

$$ \left|\left|\frac{\partial P}{\partial\vartheta}\right|\right|\leq\beta||\vartheta||, $$

where $$ \beta $$ and $$ \lambda_i $$ are positive constants, $$ i=1, \ldots, 4 $$.

Theorem 1 Consider the closed-loop system comprising the underactuated USV dynamics given in Equation (2), subject to modeling uncertainty and external disturbances. This system satisfies Assumptions 2, utilizing the intermediate control laws specified in Equations (23) and (24), the triggering instants defined in Equation (29), the ETESO in Equation (5), the MLP neural network update laws in Equation (25), and the adaptive laws given in Equation (26). By selecting appropriate parameters, all error signals in the system converge quickly to an arbitrarily small vicinity of the origin, while also ensuring that the Zeno phenomenon is excluded.

Proof 1 From Equations (4) and (5), we have for $$ \tilde{x}_i $$ and $$ t\in[t_{k}, t_{k+1}) $$

$$ \begin{cases}\dot{{\mathbf e}}_i=\frac{1}{\varepsilon }(e_{i+1}(t)-g_i\left (e_1(t)\right )-\alpha_{i}(e_1(t)-\sigma(t))), i=1, 2, \\ \dot{{\mathbf e}}_3=\dot{{\mathbf x}}_3 -\frac{1}{\varepsilon }(g_3\left (e_1(t)\right )+\alpha_{3}(e_1(t)-\sigma(t))), \end{cases}. $$

Considering a positive semidefinite function $$ V_0=P(e) $$, and $$ P(e) $$ is given by $$ P(\vartheta)=\langle\tilde{P} \vartheta, \vartheta\rangle $$$$ P(*): \mathbb{R}^{n+1} \rightarrow \mathbb{R} $$, the matrix $$ \tilde{P} $$ is a positive definite one satisfying $$ \tilde{P} N+N^T \tilde{P}=-I $$, and $$ I $$ is the identity matrix. Taking the derivative of $$ V_0 $$, we can compute its derivative, which yields:

$$ \frac{d}{dt}V_{0}(t)=\sum\limits_{i=1}^{3}\frac{\partial P}{\partial e_i}\dot{e}_i(t)=\frac{1}{\varepsilon}\sum\limits_{i=1}^{2}\frac{\partial P}{\partial e_{i}}(e_{i+1}(t)-g_{i}(e_{1}))-\frac{1}{\varepsilon}\frac{\partial P}{\partial e_{3}}g_{3}(e_{1}(t))-\frac1\varepsilon\sum\limits_{i=1}^{3}\frac{\partial P}{\partial e_i}\alpha_i(e_1(t), \sigma(t))+\frac{\partial P_0}{\partial e_{3}}\dot{{\mathbf x}}_3 . $$

From Equations (35) and (36) in Assumptions 2, we have

$$ \dot{V}_{0}(t)\leq-\frac{\lambda_{3}}{\varepsilon}||e||^{2}+\beta\dot{x}_{3}||e||-\frac{1}{\varepsilon}\sum\limits_{i=1}^{3}\frac{\partial P}{\partial e_{i}}\alpha_{i}(e_{1}(t), \sigma(t)). $$

considering the fact that the functions $$ g_i $$ are assumed to be Lipschitz with $$ Li>0 $$, it follows

$$ |\alpha_i(e_1, \sigma)|=|g_i(e_1+\sigma)-g_i(e_1)|\leq L_i|e_1+\sigma-e_1|=L_i|\sigma|. $$

By combining Equation (6), it yields

$$ -\frac{1}{\varepsilon}\sum\limits_{i=1}^{n+1}\frac{\partial P}{\partial e_{i}}\alpha_{i}(e_{1}(t), \sigma(t))\leq\left(\frac{\varepsilon^{-\frac{1}{3}}||e||^{2}}{2}+\frac{\varepsilon^{\frac{1}{3}}}{2}\right)\overline{B}\beta. $$

According to Young's inequality, it follows that $$ \beta\dot{x}_{n+1}||e||\leq\frac{\beta}{2\varepsilon}\left(\varepsilon^{\frac{2}{3}}||e||^{2}+\varepsilon^{\frac{4}{3}}|\dot{x}_{n+1}|^{2}\right) $$. Then, from Assumption 2, we can get

$$ \beta\dot{x}_{n+1}||e||\leq\frac{\beta}{2\varepsilon}\left(\varepsilon^{\frac{2}{3}}||e||^{2}+\varepsilon^{\frac{4}{3}}M^{2}\right). $$

Substituting Equations (42) and (43) into Equation (40), from Assumption 2, it yields:

$$ \dot{V}_{0}(t)\leq-\frac{2\lambda_3-\beta\varepsilon^{\frac{2}{3}}-\beta\varepsilon^{\frac{2}{3}}\overline{B}}{\varepsilon\lambda_{2}}V_0+\frac{\varepsilon^{\frac{1}{3}}}{2}\beta(M^{2}+\overline{B}), $$

where we choose parameters such that $$ \frac{2\lambda_3-\beta\varepsilon^{\frac23}-\beta\varepsilon^{\frac23}\overline{B}}{\varepsilon\lambda_2}>0 $$. Next, consider

$$ V_1=\frac{m_{11}}{2\beta_u}S_{u}^{2}+\frac{m_{33}}{2\beta_r}S_{r}^{2}+\frac{1}{2\gamma_{u}}\tilde{\tau}_{\mathrm{wu}}^{*}{}^{2}+\frac{1}{2\gamma_{r}}\tilde{\tau}_{\mathrm{wr}}^{*}{}^{2}+\frac{1}{2\hat{\lambda}}_u\tilde{\phi}_{u}^{2} +\frac{1}{2\hat{\lambda}}_r\tilde{\phi}_{r}^{2}+\eta . $$

Taking the derivative of $$ V_1 $$ and substituting Equations (23)-(26) into it have:

$$ \begin{array}{l} \dot{V}_1 =-k_uS^2_u-k_rS^2_r+S_{u}\boldsymbol{\phi}_{\mathrm{u}}^{*\mathrm{T}}\boldsymbol{h}(\boldsymbol{z})+\xi _{u}S_{u}+\dot{\eta}-\frac{1}{2}S_{u}^{2}\boldsymbol{\hat{\phi}}_{u}\boldsymbol{h}^{\mathrm{T}}\boldsymbol{h}+S_{r}\boldsymbol{\phi}_{\mathrm{r}}^{*\mathrm{T}}\boldsymbol{h}(\boldsymbol{z})+\xi_{r}S_{r}-\frac{1}{2}S_{r}^{2}\boldsymbol{\hat{\phi}}_{r}\boldsymbol{h}^{\mathrm{T}}\boldsymbol{h}\\ \;\;\;\;+\tau_{\mathrm{wu}}^{*}\left(\left|S_{u}\right|-S_{u}\phi(S_{u})\right)+\tau_{\mathrm{wr}}^{*}\left(\left|S_{r}\right|-S_{r}\phi(S_{r})\right)-\frac{1}{\hat{\lambda}_{u}}\tilde{\phi}_{u}\dot{\hat{\phi}}_{u}-\frac{1}{\hat{\lambda}_{r}}\tilde{\phi}_{r}\dot{\hat{\phi}}_{r}+\sigma_{u}\Big(\tau_{\mathrm{wu}}^{*}-\hat{\tau}_{\mathrm{wu}}^{*}\Big)\Big(\hat{\tau}_{\mathrm{wu}}^{*}-\tau_{\mathrm{wu}}^{0}\Big)\\ \;\;\;\; +\sigma_r\Big(\tau_\mathrm{wr}^*-\hat{\tau}_\mathrm{wr}^*\Big)\Big(\hat{\tau}_\mathrm{wr}^*-\tau_\mathrm{wr}^0\Big)+Z_uS_u+Z_rS_r+\dot{\eta} . \end{array} $$

From Young's inequality, we have the following inequalities:

$$ \begin{array}{l} 2S_{u}\phi_{\mathrm{u}}^{*\mathrm{T}}h\leq S_{u}^{2}\left\|\phi_{\mathrm{u}}^{*}\right\|^{2}\left\|h\right\|^{2}+1=S_{u}^{2}\phi_{u}h^{\mathrm{T}}h+1, 2S_{r}\phi_{\mathrm{r}}^{*\mathrm{T}}h\leq S_{r}^{2}\left\|\phi_{\mathrm{r}}^{*}\right\|^{2}\left\|h\right\|^{2}+1=S_{r}^{2}\phi_{r}h^{\mathrm{T}}h+1\\ -\sigma_u\widetilde{\tau}_{wu}^{*}\left(\widehat{\tau}_{wu}^{*}-\tau_{wu}^{0}\right)\leq-\frac{\sigma_{u}}{2}\widetilde{\tau}_{wu}^{*}+\frac{\sigma_{u}}{2}\Big(\tau_{wu}^{*}-\tau_{wu}^{0}\Big)^{2}-\sigma_r\widetilde{\tau}_{wr}^{*}\left(\widehat{\tau}_{wr}^{*}-\tau_{wr}^{0}\right)\leq-\frac{\sigma_{r}}{2}\widetilde{\tau}_{wr}^{*}+\frac{\sigma_{r}}{2}\Big(\tau_{wr}^{*}-\tau_{wr}^{0}\Big)^{2} . \end{array} $$

Then, we have the following inequalities

$$ \begin{array}{l} \dot{V} \leq-\eta_{1}s_{1}^{2}-\eta_{2}s_{2}^{2}+\frac{1}{2}s_{1}^{2}\tilde{\phi}_{1}h^{\mathrm{T}}h+\frac{1}{2}s_{2}^{2}\tilde{\phi}_{2}h^{\mathrm{T}}h+\frac{\zeta_{\mathrm{U}}^{2}}{2}+\frac{s_{1}^{2}}{2}+\frac{\zeta_{\mathrm{R}}^{2}}{2}+\frac{s_{2}^{2}}{2}-\frac1{{\mathbf {\hat{\lambda}}}_1}\tilde{\phi}_1\dot{\hat{\phi}}_1-\frac1{{\mathbf {\hat{\lambda}}}_2}\tilde{\phi}_2\dot{\hat{\phi}}_2+0.2785\varepsilon_1\tau_{{\mathbf wu}}^*+0.2785\varepsilon_2\tau_{{\mathbf wr}}^*\\ \;\;\;\;-\frac{\sigma_1}2\left(\hat{\tau}_{{\mathbf wu}}^*-\tau_{{\mathbf wu}}^*\right)^2+\frac{\sigma_1}2\Big(\tau_{\mathrm{wu}}^*-\tau_{\mathrm{wu}}^0\Big)^2-\frac{\sigma_2}2\Big(\hat{\tau}_{\mathrm{wr}}^*-\tau_{\mathrm{wr}}^*\Big)^2+\frac{\sigma_2}2\Big(\tau_{\mathrm{wr}}^*-\tau_{\mathrm{wr}}^0\Big)^2+1 \\ \;\;\;\;=-\frac{2\eta_{1}-1}{2}s_{1}^{2}-\frac{2\eta_{2}-1}{2}s_{2}^{2}-\frac{\sigma_{1}}{2}\left(\hat{\tau}_{\mathrm{wu}}^{*}-\tau_{\mathrm{wu}}^{*}\right)^{2}-\frac{\sigma_{2}}{2}\left(\hat{\tau}_{\mathrm{wr}}^{*}-\tau_{\mathrm{wr}}^{*}\right)^{2}+\tilde{\boldsymbol{\phi}}_1\left(\frac12s_1^2\boldsymbol{h}^\mathrm{T}\boldsymbol{h}-\frac1{\hat{\boldsymbol{\lambda}}_1}\dot{\hat{\boldsymbol{\phi}}}_1\right)+\tilde{\boldsymbol{\phi}}_2 \left(\frac12s_2^2\boldsymbol{h}^\mathrm{T}\boldsymbol{h}-\frac1{\hat{\boldsymbol{\lambda}}_2}\dot{\hat{\boldsymbol{\phi}}}_2\right)\\ \;\;\;\;+0.2785\boldsymbol{\varepsilon}_1\dot{\boldsymbol{\tau}}_\mathrm{wu}+0.2785\boldsymbol{\varepsilon}_2\dot{\boldsymbol{\tau}}_\mathrm{wr}+\frac{{\zeta_{\mathrm{U}}}^{2}}{2}+\frac{{\zeta_{\mathrm{R}}}^{2}}{2}+\frac{\sigma_{1}}{2}\Big(\tau_{\mathrm{wu}}^{*}-\tau_{\mathrm{wu}}^{0}\Big)^{2}+\frac{\sigma_{2}}{2}\Big(\tau_{\mathrm{wr}}^{*}-\tau_{\mathrm{wr}}^{0}\Big)^{2}+\frac{N_{\mathrm{u}}^{2}}{4\alpha_{\mathrm{l}}}+1 \\ \;\;\;\;=-\frac{2\eta_{1}-1}{2}s_{1}^{2}-\frac{2\eta_{2}-1}{2}s_{2}^{2}-\left(\frac{1}{T}-\alpha_{1}\right)y_{1}^{2}-\frac{\sigma_{1}}{2}\left(\hat{\tau}_{\mathrm{wu}}^{*}-\tau_{\mathrm{wu}}^{*}\right)^{2}-\frac{\sigma_{2}}{2}\left(\hat{\tau}_{\mathrm{wr}}^{*}-\tau_{\mathrm{wr}}^{*}\right)^{2}+\kappa_1\tilde{\phi}_1\hat{\phi}_1+\kappa_2\tilde{\phi}_2\hat{\phi}_2+0.2785\varepsilon_1\tau_\mathrm{wu}^*\\ \;\;\;\;+0.2785\varepsilon_2\tau_\mathrm{wr}^*+\frac{\zeta_\mathrm{U}^2}2+\frac{\zeta_\mathrm{R}^2}2+\frac{\sigma_1}2\left(\tau_\mathrm{wu}^*-\tau_\mathrm{wu}^0\right)^2+\sigma_2\left(\tau_{\mathrm{wr}}^*-\hat{\tau}_{\mathrm{wr}}^*\right)\left(\hat{\tau}_{\mathrm{wr}}^*-\tau_{\mathrm{wr}}^0\right)-\frac1{\hat{\lambda}_1}\tilde{\phi}_1\dot{\hat{\phi}}_1-\frac1{\hat{\lambda}_2}\tilde{\phi}_2\dot{\hat{\phi}}_2 \end{array} $$

$$ \tilde{\phi}_1\hat{\phi}_1 $$ and $$ \tilde{\phi}_2\hat{\phi}_2 $$ satisfy the following inequality

$$ 2\tilde{\phi}_1\hat{\phi}_1= \phi_1^2-\hat{\phi}_1^2-\tilde{\phi}_1^2\leq\phi_1^2-\tilde{\phi}_1^2, 2\tilde{\phi}_{2}\hat{\phi}_{2}=\phi_2^2-\tilde{\phi}_2^2-\tilde{\phi}_2^2\leq\phi_2^2-\tilde{\phi}_2^2 . $$

Substituting Equation (49) into Equation (48) yields

$$ \begin{array}{l} \dot{V} =-\frac{2\eta_{1}-1}{2}s_{1}^{2}-\frac{2\eta_{2}-1}{2}s_{2}^{2}-\frac{\sigma_{1}}{2}\left(\hat{\tau}_{\mathrm{wu}}^{*}-\hat{\tau}_{\mathrm{wu}}^{*}\right)^{2}-\frac{\sigma_{2}}{2}\left(\hat{\tau}_{\mathrm{wr}}^{*}-\hat{\tau}_{\mathrm{wr}}^{*}\right)^{2}-\frac{\kappa_1}2\tilde{\phi}_1^2-\frac{\kappa_2}2\tilde{\phi}_2^2+0.2785\varepsilon_1\tau_\mathrm{wu}^*\\ \;\;\;\; +0.2785\varepsilon_2\tau_\mathrm{wr}^*+\frac{\xi_\mathrm{U}^2}2+\frac{\xi_\mathrm{R}^2}2+\frac{\sigma_1}2\left(\tau_\mathrm{wu}^*-\tau_\mathrm{wu}^0\right)^2+\frac{\sigma_{2}}{2}\Big(\tau_{\mathrm{wr}}^{*}-\tau_{\mathrm{wr}}^{0}\Big)^{2}+\frac{N_{\mathrm{u}}^{2}}{4\alpha_{1}}+1+\frac{\kappa_{1}}{2}\phi_{1}^{2}+\frac{\kappa_{2}}{2}\phi_{2}^{2}\end{array} . $$

Letting $$ V=V_0+V_1 $$ and taking the derivative of $$ V $$, we can get:

$$ \dot{V}\le-C_1V+C_2, $$

where $$ C_1=min\left \{\frac{2\eta_{1}-1}{2} , \frac{2\eta_{2}-1}{2}, \frac{\sigma_{1}}{2}, \frac{\sigma_{2}}{2}, \frac{2\lambda_3-\beta\varepsilon^{\frac{2}{3}}-\beta\varepsilon^{\frac{2}{3}}\overline{B}}{\varepsilon\lambda_{2}}\right \} $$ and $$ C_2=0.2785\varepsilon_1\tau_\mathrm{wu}^*+0.2785\varepsilon_2\tau_\mathrm{wr}^*+\frac{\xi_\mathrm{U}^2}2+\frac{\xi_\mathrm{R}^2}2+\frac{\sigma_1}2\left(\tau_\mathrm{wu}^*-\tau_\mathrm{wu}^0\right)^2+\frac{\sigma_{2}}{2}\Big(\tau_{\mathrm{wr}}^{*}-\tau_{\mathrm{wr}}^{0}\Big)^{2}+\frac{N_{\mathrm{u}}^{2}}{4\alpha_{1}}+1+\frac{\kappa_{1}}{2}\phi_{1}^{2}+\frac{\kappa_{2}}{2}\phi_{2}^{2}+\frac{\varepsilon^{\frac{1}{3}}}{2}\beta(M^{2}+\overline{B}) $$. Solving Equation (51) gives

$$ 0<V(t)\le(V(0)-\frac{C_2}{C_1})e^{-C_1t}+\frac{C_2}{C_1}. $$

Based on Lemma 3, Equations (51) and (52), we conclude that $$ V(t) $$ is uniformly ultimately bounded, which implies that all error signals $$ S_u $$, $$ S_r $$, $$ \widetilde{\phi}_u $$, $$ \widetilde{\phi}_r $$, $$ \widetilde{\tau}^*_u $$ and $$ \widetilde{\tau}^*_r $$ are ultimately uniformly bounded. Furthermore, since the error signals $$ u_e $$ and $$ r_e $$ are ultimately bounded, it follows that $$ x_e $$ and $$ y_e $$ are also bounded. By appropriately choosing the parameter $$ \frac{C_2}{C_1} $$, it can be made arbitrarily small, ensuring that the tracking error becomes negligible, allowing the USV to track the trajectory with high accuracy.

Next, we first prove that there exists a Zeno-free phenomenon in the ETESO. By the definition of ETESO in Equation (6) and the definition of Equation ($$ \sigma(t) $$), we can get

$$ \left | y(t(o, k)-y(t))\right |\le \frac{\bar{B}\varepsilon^3}{ { \sum\nolimits_{3}^{i=1}}L_i } . $$

Then, for the error $$ \sigma(t) $$ in Equation (6), as $$ t\in[t_{o, k}, t_{o, k+1}) $$, the following result $$ \left | y(t(o, k)-y(t))\right |=\int_{t_{o, k}}^{t}|\dot{x}_1| d\tau=\int_{t_{o, k}}^{t}|x_2| d\tau $$ holds. Since the system states are bounded, there exists $$ B>0 $$ such that the inequality $$ \left | y(t(o, k)-y(t))\right |=\int_{t_{o, k}}^{t}|x_2| d\tau \le (t-t(o, k))B $$ holds. Therefore, for $$ t=t(o, k)+\breve{\tau} $$, there exists $$ \breve{\tau}=\frac{\bar{B}\varepsilon^3}{B { \sum_{3}^{i=1}}L_i} $$ such that $$ y(t)-y(to, k)\le\frac{\bar{B}\varepsilon^3}{ { \sum_{3}^{i=1}} } $$. Accordingly, it comes to the result that the interval between $$ t_{0, k} $$ and $$ t_{0, k+1} $$ satisfies $$ \{t_{o, k+1}-t_{o, k}\}\geq\breve{\tau}>0 $$.

Finally, we prove that no Zeno phenomenon occurs by contradiction. Suppose that there exists Zeno phenomenon, then $$ \lim_{j \to\infty} t_j=T_0<\infty $$, where $$ T_0 $$ is a positive constant, and $$ \lim_{j \to\infty} \bigtriangleup j=0 $$ with $$ \bigtriangleup j=t_{j+1}-t_j $$. By the existence of the limit, for any constant $$ \epsilon _0>0 $$, there exists a positive integer $$ N(\epsilon_0) $$, such that

$$ t_j\in \left [ T_0-\epsilon_0, T_0\right ], \forall \sigma >N(\epsilon_0 ). $$

Define the triggering instant as $$ \tau=t_{j+1} $$. Note that once one event is triggered, then from Equation (29), the measurement error $$ Z(t) $$ is reset to zero, that is, $$ ||Z(\tau^+)||=0 $$. Just before the triggering, the following inequality can be obtained from Equation (30)

$$ \eta_{i}(\tau^-)+(k_{i}-1) S_{i}^{2}(\tau^-)-Z^{2}_{i}(\tau^-)\le0, $$

it is easy to see that

$$ \eta_{i}(\tau^-)\le Z^{2}_{i}(\tau^-), $$

there exists a constant $$ L_i>0 $$, such that

$$ Z^2_n(t) =\left (||Z(t_j)-Z(t)|| \right)^2\le L_i(t-t_j) . $$

From Equations (56) and (57), we can have

$$ L_i \Delta_j \geq Z^2_n(t^-_{j+1}) \geq \eta(t^-_{j+1})\geq\frac{\eta_{0}e^{-(\chi+\frac{1}{\lambda})t_{j+1}}}{\lambda}\geq\frac{\eta_{0}e^{-(\chi+\frac{1}{\lambda})T_{0}}}{\lambda}, $$

It is easy to see by Equation (58) that when $$ \operatorname*{lim}_{j\to\infty} t_{j}=T_{0}<\infty $$, the assumption about $$ \lim_{j\to\infty} \Delta_j=0 $$ cannot be obtained, which results in a contradiction. Therefore, the Zeno phenomenon does not occur. Based on the above calculation, Theorem 1 is proved.

5. SIMULATION ANALYSIS

To provide a more convincing demonstration of the effectiveness of the proposed trajectory tracking controller, the algorithm presented in this paper is compared with a commonly used event-triggered sliding mode algorithm based on LOS. The benchmark algorithm integrates the LOS guidance rate[38], event-triggered mechanism[26], linear sliding surface[28], and state observer extension[41]. In order to verify the effectiveness of the proposed unmanned boat tracking control strategy and formation control strategy, a model boat Cybership Ⅱ USV developed by the Norwegian University of Science and Technology was selected for the simulation experimental study. The model mass $$ m $$ = 23.8 kg, length $$ L_i $$ = 1.225 m width $$ B_i $$ = 0.29 m, the center of gravity is at a distance $$ x $$ axis from the origin $$ x_g $$ = 0.046 m, and the inertia moment $$ I_z $$ = 0.046 m. The relevant model parameters are in ref.[42].

In the ETESO, $$ \alpha_i(\vartheta) $$ is selected as $$ \alpha_1=3 $$, $$ \alpha_2(\vartheta)=3 $$, $$ \alpha_{3}=1 $$, $$ L_1=3.25 $$, $$ L_2=3 $$, $$ L_3=1 $$, $$ \overline{B}=10 $$, $$ \varepsilon_e=0.1 $$ We set the reference trajectory as $$ x_d=300sin(0.03t) $$ and $$ y_d=300cos(0.03t) $$. The initial position of the USV is $$ x(0) $$=1 m, $$ y(0) $$=29 m, $$ \psi(0) $$=0 rad, $$ u(0) $$=0 m/s, $$ v(0) $$=0 m/s and $$ r(0) $$=0 rad/s. The control parameters are $$ k_u=10 $$, $$ k_r=10 $$, $$ \lambda_u=10 $$, $$ \lambda_r=10 $$, $$ \beta_u=1.5 $$, $$ \beta_r=1.5 $$, $$ r_u=1.1 $$, $$ r_r=1.2 $$. Parameters of adaptive laws are $$ \tau_{\mathrm{wu}}^{0}=0.1 $$, $$ \quad\tau_{\mathrm{wr}}^{0}=0.1 $$, $$ \gamma_{u}=12 $$, $$ \gamma_{r}=12 $$, Parameters of update laws are $$ \hat{\lambda}_{u}=10 $$ and $$ \hat{\lambda}_{r}=10 $$. To testify to the robustness and effectiveness of the proposed controller, the disturbances are assumed to be as follows: $$ \Delta f_{\mathrm{u}}=0.5+0.25\sin(0.01t) $$, $$ \Delta f_{\mathrm{r}}=-40+2.5\cos(0.10t) $$. $$ \tau_{wu}=\sin\bigl(0.2t\bigr)+\cos\bigl(0.5t\bigr) $$, $$ \tau_{wr}=\sin\bigl(0.5t\bigr)+\cos\bigl(0.3t\bigr) $$. The time-varying target trajectory is chosen as follows:

$$ \begin{bmatrix} x_{d}\\y_{d} \end{bmatrix}=\begin{bmatrix}30sin(0.03t)\\-40+2.5cos(0.03t)\end{bmatrix}. $$

The numerical results are shown in Figures 2-8. Figure 2 shows the trajectory tracking comparison simulation results between the proposed algorithm and the control algorithm under model uncertainty and external disturbances. The experimental results indicate that both controllers can accurately track the desired trajectory. However, the proposed solution demonstrates better performance in terms of convergence speed during the initial control phase. Meanwhile, the corresponding tracking error is shown in Figure 3, which indicates that the USV can track the reference position within approximately six seconds and achieve faster error stabilization compared to the control scheme. Figure 4 compares the velocity tracking results between the proposed algorithm and the control algorithm. From the zoomed-in view of the components, it can be seen that both the proposed algorithm and the control algorithm can track the velocity well, but the proposed algorithm requires fewer trigger occurrences. The control inputs $$ \tau_u $$ and $$ \tau_r $$ are shown in Figure 5. It can be observed that, compared to the control scheme, the proposed scheme exhibits smaller fluctuations in the control signals. As shown in Figure 6, the model uncertainty of USVs can be accurately estimated and approximated using MLP neural networks. According to the above simulation results, it can be concluded that the proposed controller has faster convergence and better robustness. Figure 7 shows the estimation of the lumped disturbance. Therefore, both the proposed and contrast algorithms can observe and estimate the unmeasurable velocity and lumped disturbance with minimal error, demonstrating excellent estimation performance. Figure 8 shows that the adaptive law can estimate the upper bound of the external disturbance by selecting the appropriate parameters. Figure 9 shows triggering instants and triggering time of $$ \tau_u $$, $$ \tau_r $$ and ETESO. It is clear that, compared to the static event-triggering mechanism in the comparative experiments, the DETM used in the proposed algorithm results in longer intervals between controller-triggered events.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 2. The reference and actual trajectories in the x-y plane.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 3. The curves of reference and actual positions.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 4. The comparison of actual and estimated velocities in the ETESO. ETESO: Event-triggered extended state observer.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 5. Controller inputs $$ \tau_u $$ and $$ \tau_r $$.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 6. The model uncertainty curves approximated by MLP. MLP: Minimum learning parameter.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 7. The lumped disturbances and observer estimations.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 8. The external disturbances and their boundary estimations.

Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network

Figure 9. Triggering instants and triggering time of $$ \tau_u $$, $$ \tau_r $$ and ETESO. ETESO: Event-triggered extended state observer.

6. CONCLUSION

Based on the nonlinear mathematical model with model uncertainties and external disturbances, a INTSM control method based on DETM is proposed. This approach integrates ETESO, NTSM, DETM, MLP neural networks, and adaptive techniques. To address the issue of difficult-to-measure velocity in practical applications, an ETESO is designed to accurately estimate the unmeasurable velocity and lumped disturbances. To improve the system's convergence speed and reduce steady-state error, INTSM is used to design the control input. Additionally, a DETM is incorporated into the controller, considering communication and actuator wear issues in the USV. To enhance the robustness of the control system, an MLP technique is employed to approximate and compensate for model uncertainties. An adaptive law is also designed to compensate for the neural network approximation errors and disturbances. Furthermore, through rigorous theoretical analysis, the overall stability of the closed-loop control system is proved, and it is demonstrated that the tracking error converges to a small neighborhood of the origin. Finally, comparative numerical simulations are conducted to analyze the tracking performance. The simulation results fully demonstrate the effectiveness and superiority of the proposed approach. In the future, to further conserve communication resources, we will consider replacing the event-triggered mechanism in ETESO with a DETM and verifying the proposed algorithm in practical applications.

DECLARATIONS

Authors' contributions

Conceptualization, methodology, manuscript writing and technical support: Wang X

Data acquisition: Liu R

Reviewing and editing: Li Q, Aleksander S

Supervision: Jiang R

Availability of data and materials

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Financial support and sponsorship

This work was supported by the Shandong Provincial Natural Science Foundation (Grant No. ZR2024QF114) and the major innovation project for the science education industry integration pilot project of Qilu University of Technology (Shandong Academy of Sciences) (Grant No. 2023JBZ03).

Conflicts of interest

All authors declared that there are no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2024.

REFERENCES

1. Liu Z, Zhang Y, Yu X, Yuan C. Unmanned surface vehicles: an overview of developments and challenges. Annu Rev Control 2016;41:71-93.

2. Su M, Pu R, Wang Y, Yu M. A collaborative siege method of multiple unmanned vehicles based on reinforcement learning. Intell Robot 2024;4:39-60.

3. Barrera C, Padrón Armas I, Luis F, Llinas O, Marichal Plasencia GN. Trends and challenges in unmanned surface vehicles (USV): from survey to shipping. TransNav 2021;15:135-42.

4. Xue K, Liu J, Xiao N, Ji X, Qian H. A bio-inspired simultaneous surface and underwater risk assessment method based on stereo vision for USVs in nearshore clean waters. IEEE Robot Autom Lett 2023;8:360-7.

5. Zhao L, Bai Y. Unlocking the ocean 6G: a review of path-planning techniques for maritime data harvesting assisted by autonomous marine vehicles. J Mar Sci Eng 2024;12:126.

6. Zhang G, Liu S, Huang J, Zhang W. Dynamic event-triggered path-following control of underactuated surface vehicle with the experiment verification. IEEE T Veh Technol 2022;71:10415-25.

7. Wang YL, Han QL. Network-based fault detection filter and controller coordinated design for unmanned surface vehicles in network environments. IEEE Trans Ind Informat 2016;12:1753-65.

8. Wang X, Fei Z, Gao H, Yu J. Integral-based event-triggered fault detection filter design for unmanned surface vehicles. IEEE Trans Ind Informat 2019;15:5626-36.

9. Jiang C, Zhang G, Huang C, Zhang W. Memory-based event-triggered path-following control for a USV in the presence of DoS attack. Ocean Eng 2024;310:118627.

10. Zhou L, Shen Z, Nie Y, Yu H. Event-triggered adaptive dynamic programming for optimal tracking control of unmanned surface vessel with input constraints. Trans Inst Meas Control 2023;45:2835-47.

11. Girard A. Dynamic triggering mechanisms for event-triggered control. IEEE Trans Autom Control 2015;60:1992-7.

12. Cao L, Yao D, Li H, Meng W, Lu R. Fuzzy-based dynamic event triggering formation control for nonstrict-feedback nonlinear MASs. Fuzzy Set Syst 2023;452:1-22.

13. He W, Xu B, Han QL, Qian F. Adaptive consensus control of linear multiagent systems with dynamic event-triggered strategies. IEEE Trans Cybern 2020;50:2996-3008.

14. Wang XY, Chang XH. Nonlinear continuous-time system H∞ control based on dynamic quantization and event-triggered mechanism. Neural Proc Lett 2023;55:12223-38.

15. Alfaro-Cid E, McGookin EW, Murray-Smith DJ, Fossen TI. Genetic algorithms optimisation of decoupled Sliding Mode controllers: simulated and real results. Control Eng Pract 2005;13:739-48.

16. Zhang Y, McEwen RS, Ryan JP, Bellingham JG. Design and tests of an adaptive triggering method for capturing peak samples in a thin phytoplankton layer by an autonomous underwater vehicle. IEEE J Oceanic Eng 2010;35:785-96.

17. Rodriguez J, Castañeda H, Gonzalez-Garcia A, Gordillo JL. Finite-time control for an unmanned surface vehicle based on adaptive sliding mode strategy. Ocean Eng 2022;15:111255.

18. Liu W, Ye H, Yang X. Super-twisting sliding mode control for the trajectory tracking of underactuated USVs with disturbances. J Mar Sci Eng 2023;11:636.

19. Yu XN, Hao LY. Integral sliding mode fault tolerant control for unmanned surface vessels with quantization: less iterations. Ocean Eng 2022;260:111820.

20. Zhang Q, Zhang S, Liu Y, Zhang Y, Hu Y. Adaptive terminal sliding mode control for USV-ROVs formation under deceptive attacks. Front Mar Sci 2024;11:1320361.

21. Sun Q, Liu ZW, Chi M, Ge MF, He D. Robust coverage control of multiple USVs with time- varying disturbances. Intell Robot 2023;3:242-56.

22. Han J, Xiong J, He Y, Gu F, Li D. Nonlinear modeling for a water-jet propulsion USV: an experimental study. IEEE Trans Ind Electron 2017;64:3348-58.

23. Chen Z, Zhang Y, Zhang Y, Nie Y, Tang J, Zhu S. Disturbance-observer-based sliding mode control design for nonlinear unmanned surface vessel with uncertainties. IEEE Access 2019;7:148522-30.

24. Chen Z, Zhang Y, Nie Y, Tang J, Zhu S. Adaptive sliding mode control design for nonlinear unmanned surface vessel using RBFNN and disturbance-observer. IEEE Access 2020;8:45457-67.

25. Chang XH, Han X. Observer-based fuzzy l2-l control for discrete-time nonlinear systems. IEEE Trans Fuzzy Syst 2024;32:2523-8.

26. Ning J, Ma Y, Li T, Chen CLP, Tong S. Event-triggered based trajectory tracking control of under-actuated unmanned surface vehicle with state and input quantization. IEEE Trans Intell Veh 2024;9:3127-39.

27. Wang H, Zhang S. Event-triggered reset trajectory tracking control for unmanned surface vessel system. In: Proceedings of the Institution of Mechanical Engineers, Part Ⅰ: Journal of Systems and Control Engineering. 2021;235:633–45.

28. Zhao Y, Sun X, Wang G, Fan h. Adaptive backstepping sliding mode tracking control for underactuated unmanned surface vehicle with disturbances and input saturation. IEEE Access 2021;9:1304-12.

29. Zhu F, Peng Y, Cheng M, Luo J, Wang Y. Finite-time observer-based trajectory tracking control of underactuated USVs using hierarchical non-singular terminal sliding mode. Cyber Phys Syst 2022;8:263-85.

30. Wan L, Su Y, Zhang H, Tang Y, Shi B. Global fast terminal sliding mode control based on radial basis function neural network for course keeping of unmanned surface vehicle. Int J Adv Robot Syst 2019;16:1729881419829961.

31. Zhao C, Yan H, Gao D, Wang R, Li Q. Adaptive neural network iterative sliding mode course tracking control for unmanned surface vessels. J Math 2022;2022:1417704.

32. Bateman A, Hull J, Lin Z. A backstepping-based low-and-high gain design for marine vehicles. Int J Robust Nonlinear Control 2009;19:480-93.

33. Zhang J, Liu X, Wang X, Wang Y, Wang Y. Adaptive prescribed performance tracking control for underactuated unmanned surface ships with input quantization. Intell Robot 2024;4:146-63.

34. Huang T, Xue Y, Xue Z, Zhang Z, Miao Z, Liu Y. USV-tracker: a novel USV tracking system for surface investigation with limited resources. Ocean Eng 2024;312:119196.

35. Qiu B, Wang G, Fan Y, Mu D, Sun X. Adaptive sliding mode trajectory tracking control for unmanned surface vehicle with modeling uncertainties and input saturation. Appl Sci 2019;9:1240.

36. Li J, Zhang G, Li B. Robust adaptive neural cooperative control for the USV-UAV based on the LVS-LVA guidance principle. J Mar Sci Eng 2022;10:51.

37. Wu GX, Ding Y, Tahsin T, Atilla I. Adaptive neural network and extended state observer-based non-singular terminal sliding modetracking control for an underactuated USV with unknown uncertainties. Appl Ocean Res 2023;135:103560.

38. Yu S, Lu J, Zhu G, Yang S. Event-triggered finite-time tracking control of underactuated MSVs based on neural network disturbance observer. Ocean Eng 2022;253:111169.

39. Guo BZ, Zhao Zl. On the convergence of an extended state observer for nonlinear systems with uncertainty. Syst Control Lett 2011;60:420-30.

40. Huang Y, Wang J, Shi D, Shi L. Toward event-triggered extended state observer. IEEE Trans Autom Control 2018;63:1842-49.

41. Xu Z, Xie N, Shen H, Hu X, Liu Q. Extended state observer-based adaptive prescribed performance control for a class of nonlinear systems with full-state constraints and uncertainties. Nonlinear Dynam 2021;105:345-58.

42. Skjetne R, Fossen TI, Kokotović PV. Adaptive maneuvering, with experiments, for a model ship in a marine control laboratory. Automatica 2005;41:289-98.

Cite This Article

Research Article
Open Access
Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network
Xingmin Wang, ... Ru Jiang

How to Cite

Wang, X.; Liu, R.; Sladkowski, A.; Li, Q.; Jiang, R. Dynamic event-triggered integral non-singular terminal sliding mode tracking of unmanned surface vehicle via an event-triggered extended state observer and adaptive neural network. Intell. Robot. 2024, 4, 439-56. http://dx.doi.org/10.20517/ir.2024.26

Download Citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click on download.

Export Citation File:

Type of Import

Tips on Downloading Citation

This feature enables you to download the bibliographic information (also called citation data, header data, or metadata) for the articles on our site.

Citation Manager File Format

Use the radio buttons to choose how to format the bibliographic data you're harvesting. Several citation manager formats are available, including EndNote and BibTex.

Type of Import

If you have citation management software installed on your computer your Web browser should be able to import metadata directly into your reference database.

Direct Import: When the Direct Import option is selected (the default state), a dialogue box will give you the option to Save or Open the downloaded citation data. Choosing Open will either launch your citation manager or give you a choice of applications with which to use the metadata. The Save option saves the file locally for later use.

Indirect Import: When the Indirect Import option is selected, the metadata is displayed and may be copied and pasted as needed.

About This Article

© The Author(s) 2024. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Data & Comments

Data

Views
43
Downloads
4
Citations
0
Comments
0
2

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

0
Download PDF
Share This Article
Scan the QR code for reading!
See Updates
Contents
Figures
Related
Intelligence & Robotics
ISSN 2770-3541 (Online)
Follow Us

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/