A Bounded-Function-Based Scheme for Finite-Time Stabilization of a NWMR with Input Constraints

Article Open Access

A Bounded-Function-Based Scheme for Finite-Time Stabilization of a NWMR with Input Constraints

Author Information
The Seventh Research Division, Beihang University (BUAA), Beijing 100191, China
*
Authors to whom correspondence should be addressed.
Views:217
Downloads:25
Drones and Autonomous Vehicles 2025, 2 (1), 10001;  https://doi.org/10.70322/dav.2025.10001

Received: 14 October 2024 Accepted: 02 January 2025 Published: 06 January 2025

Creative Commons

© 2025 The authors. This is an open access article under the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/).

ABSTRACT: This paper addresses the finite-time stabilization problem for a nonholonomic wheeled mobile robot (NWMR) with input constraints. By utilizing the hyperbolic tangent function tanh(·), bounded finite-time stabilization controllers are developed. In addition, an explicit upper-bound estimate for the closed-loop settling time is given, and the level of input constraints is characterized by parameters that depend on the actuator’s capacity. A thorough finite-time stability analysis is carried out using appropriate Lyapunov functions. For a compact set contained in the domain of attraction, a guideline is presented to clarify how to construct it. Finally, simulation results show the effectiveness of the developed controllers.
Keywords: Finite-time stability/Stabilization; Bounded input; Input constraints

1. Introduction

Nonholonomic wheeled mobile robot (NWMR) is a typical nonlinear control system subject to nonholonomic constraints. It has been widely applied in daily life, scientific research and engineering [1]. Since the stabilization of a NWMR is an important topic, then there have been fruitful results [2,3,4,5]. It should be pointed out that many of the existing results stabilize the NWMR asymptotically, which implies that the closed-loop system has an infinite settling time. However, for accomplishing specific tasks within a limited timeframe, finite settling time is preferred. Therefore, finite-time control is more suitable for practical applications. Lots of methods have been used to address the finite-time control problems, such as the sliding mode control method [6,7,8], the homogeneity method [9,10] and the backstepping method [11]. In addition, these methods have been applied in various systems, including quadrotors [12], spacecraft [6], brushless DC (BLDC) motors [13], and industrial robots [14]. Due to the nonholonomic constraints of the mobile robot, only limited results have been achieved in finite-time stabilization for the NWMR [15,16,17]. It is noteworthy that many of these results address finite-time stabilization without considering input constraints. Input constraints are critical factors in controller design. Ignoring these constraints during the controller design can result in degraded or even deteriorated control performance. Therefore, it is imperative to design controllers considering input constraints. Most importantly, addressing the finite-time control problem while considering input constraints is challenging. For NWMR stabilization, a switching strategy has been employed in [18], introducing state and input transformations to decouple the original system into a first-order integrator system and a second-order integrator system. However, when input constraints are considered, the input transformation makes it difficult to determine whether the original inputs exceed the allowable limits. To address this issue, this paper proposes a bounded finite-time controller for NWMR stabilization without introducing input transformation. Based on the above clarifications, the main contributions of this paper are summarized as follows: (1) Bounded finite-time controllers are proposed, using the hyperbolic tangent function tanh(·), to stabilize the mobile robot and the upper-bound estimate for the closed-loop settling time is given explicitly. In addition, the level of the input constraints is given by parameters selected according to the actuator’s capacity. (2) For the proposed bounded finite-time controllers, a compact set contained in the domain of attraction is explicitly characterized, and a step-by-step guideline is provided for constructing this compact set. The remainder of this paper is organized as follows. In Section 2, some useful lemmas are recalled and the problem formulation is conducted. Section 3 presents the main results. Section 4 presents the simulation results to validate the effectiveness of the proposed controllers. Section 5 concludes the paper.

2. Preliminaries

Consider a general system defined by
```latex\dot{x}=f(x),x(0)=x_0,```
where $$x \in \mathbb{R}^{n}$$ is the system state, f:DRn represents a continuous nonlinear function on an open neighborhood $$D \subset \mathbb{R}^{n}$$ of the origin and $$f \left(0\right) = 0$$. In this section, some important lemmas are recalled and the problem formulation is conducted, which will facilitate the derivation of the main results. 2.1. Some Useful Lemmas Lemma 1. [19] For a system (1), if there exists a continuous positive definite function V:DR on an open neighborhood $$U\subset D$$ of the origin and real numbers $$c>0,\rho\in\begin{pmatrix}0,1\end{pmatrix}$$ such that $$\dot{V}(x)+cV^{\rho}(x)\leq0,x\in U\backslash\{0\}$$, then the origin of the system (1) is finite-time stable and the settling time can be estimated as $$T\left(x_0\right)\le V^{1-\rho}\left(x_0\right)/c\left(1-\rho\right)$$ for a given initial condition $$x_0\in U$$, where $$T(x_0)$$ is the settling time function. Lemma 2. [20] If $$0<\lambda=\lambda_1/\lambda_2\le1$$ and both $$\lambda_1,\lambda_2$$ are positive odd integers, then the following inequality holds:
```latex\left|x^\lambda-y^\lambda\right|\leq2^{1-\lambda}\left|x-y\right|^\lambda.```
Lemma 3. [21] For any $$x\in\mathbb{R},y\in\mathbb{R},c,d$$ are positive real numbers and $$\Lambda(x,y)>0$$ is a real-valued function, the following inequality can be obtained
```latex|x|^c|y|^d\leq\frac{c}{c+d}\Lambda(x,y)|x|^{c+d}+\frac{d}{c+d}\Lambda^{-\frac{c}{d}}(x,y)|y|^{c+d}.```
Lemma 4. [22,23] For any number $$r>0$$, when $$\left|x\right|\le r$$, the following inequality holds:
```latex\frac{\tanh(k_2r)}{r}|x|\leq|\tanh(k_2x)|,```
where $$k_{2}>0$$ is a constant. Lemma 5. [24] If $$x_i \in R, i =1,2,…n$$ and $$0<\lambda\leq1$$ is a real number, then the following inequality holds:
```latex(\sum\nolimits_{i=1}^n\lvert x_i\rvert)^\lambda\leq\sum\nolimits_{i=1}^n\lvert x_i\rvert^\lambda.```
2.2. Problem Formulation The kinematics of a mobile robot can be represented as
```latex\begin{cases}\dot{x}=v\mathrm{cos}\theta,\\\dot{y}=v\mathrm{sin}\theta,\\\dot{\theta}=\omega,&&\end{cases}```
where $$\begin{pmatrix}x,y\end{pmatrix}$$ represents the coordinate of the center of mass, $$\theta$$ is the robot’s orientation in the global frame, v and $$\omega$$ are the linear and angular velocities, respectively. It is noteworthy that the kinematic model (2) is quite common in real-world applications and can often be found in the control of Automated Guided Vehicles (AGVs) and robot vacuum cleaners. It can be seen from (2) that there exists the following nonholonomic constraint:
```latex\dot{x}\sin\theta-\dot{y}\cos\theta=0.```
This nonholonomic constraint reflects the limitation of the robot’s freedom of motion: the robot can only move along its heading direction and cannot slide laterally. In this paper, considering control input constraints, the control objective is to design bounded controllers v and $$\omega$$ to stabilize (2) within finite time. In addition, provide an explicit upper-bound estimate for the closed-loop settling time and ensure that the saturation level of the control inputs can be predefined based on the actuator’s capacity. Finally, offers a guideline for constructing a compact set contained in the domain of attraction.

3. Main Results

In this section, a bounded finite-time control scheme is proposed to stabilize the mobile robot. Before designing the controllers v and $$\omega$$, a classical transformation is introduced as follows:
```latex\begin{cases}\,\,\,\,\,\,\,\,\,\,\,\,\,\,z_1=\theta,\\z_2=x\mathrm{sin}\theta-y\mathrm{cos}\theta,\\z_3=x\mathrm{cos}\theta+y\mathrm{sin}\theta.&\end{cases}```
Under this transformation (3), the mobile robot system (2) can be written as
```latexS_1:\dot{z}_1=\omega,S_2:\dot{z}_2=\omega z_3,\dot{z}_3=-\omega z_2+\nu,```
where $$z_1,z_2,z_3$$ are the system states, $$\omega$$ and v are the control inputs. It can be seen that through (3) the mobile robot (2) is decoupled into a first-order subsystem S1 and a second-order subsystem S2. Inspired by [18], a switching strategy is used here. Specifically, to achieve the finite-time stabilization of the mobile robot, one needs to finish the following two steps: (1) Set $$\omega=c$$, where c is a non-zero constant. Then (4) can be represented as
```latexS_1:\dot{z}_1=c,\overline{S}_2:\dot{z}_2=\overline{z}_3,\dot{\overline{z}}_3=-c^2z_2+c\nu,```
where $$\overline{z}_3=cz_3$$. One needs to design a bounded controller v to stabilize $$\bar{S}_2$$ within a finite time $$T_1(z_2(0),\bar{z}_3(0))$$. (2) When $$t>T_1\left(z_2(0),\bar{z}_3(0)\right)$$, the second-order subsystem $$\bar{S}_2$$ has been stabilized at the origin, i.e., $$z_2=\overline{z}_3=0\Rightarrow z_2=z_3=0$$. Then, at this instant, one needs to use a bounded finite-time controller $$\omega$$ to stabilize the first-order subsystem S1. From the analysis of the above two steps, design switching controllers as
```latex\omega=\begin{cases}\quad c,t\leq T_1,\\-k_{\omega1}\mathrm{tanh}^{\frac{m}{n}}(k_{\omega2}z_1),t>T_1,&\end{cases}```
where $$k_{\omega1},k_{\nu1}$$ are positive constants, and c is a non-zero constant, and m,n are positive odd constants satisfying m < n. 1 < p = p1/p2 < 2 and both $$p_{1}$$ and $$p_{2}$$ are positive odd integers satisfying $$p_1>p_2$$, and $$k_{\omega2},k_{\nu2},k_{\nu3}$$ are appropriate positive constants. Without any ambiguity, we use $$\left(\cdot\right)^p$$ for simplicity instead of $$\left|\cdot\right|^p\operatorname{sign}(\cdot)$$, where $$p=p_1/p_2$$. Before presenting the main results, Lemma 6 and Lemma 7 are given as follows. Lemma 6. For the subsystem $$\bar{S}_2$$ in (5), if the initial states $$\left(z_2(0),\bar{z}_3(0)\right)\in\Omega_\beta$$, where $$\Omega_{\beta}$$ is a given compact set, then $$\bar{S}_2$$ can be finite-time stabilized by the bounded controller (7) at the origin. The controller parameters $$k_{v1}$$, $$k_{v2}$$ and $$k_{v3}$$ satisfy.
```latexk_{v3}>\frac{2^{1-\frac{1}{p}}p+1}{1+p}+\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^p}\frac{2-p}{1+p},```
```latex\frac{k_{v1}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^p}\left(\frac{\tanh(k_{v2}r)}{r}\right)^{\frac{2}{p}-1}>\frac{2^{1-\frac{1}{p}}}{1+p}+2^{1-\frac{1}{p}}+\frac{pk_{v3}^{1+\frac{1}{p}}}{1+p}+\frac{(2p-1)c^{\frac{2(1+p)}{2p-1}}E^{\frac{2\left(p^2-1\right)}{{(2p-1)p}}}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^p(1+p)},```
where E is a given positive constant, and r is a positive constant satisfying $$r\leq k_{\nu3}^pE$$, and $$k_{\nu1}>c^2E$$, and the control input is bounded by
```latex|v|\leq\frac{k_{v1}}{|c|}.```
Proof. The proof is divided into three steps. In the first two steps, the time derivative of the Lyapunov function is given. In Step 3, the finite-time stability of the closed-loop system $$\bar{S}_2$$ and (7) is proved in the given domain $$\Omega_{\beta}$$. Step 1: Substituting controller (7) into $$\bar{S}_2$$ in (5) derives that
```latex\begin{cases}\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\dot{z}_2=\bar{z}_3,\\\dot{\bar{z}}_3=-c^2z_2-k_{v1}\mathrm{tanh}^{\frac{2}{p}-1}\left(k_{v2}\left(\bar{z}_3^p+k_{v3}^pz_2\right)\right).&\end{cases}```
Choose an appropriate Lyapunov function V as
Taking the time derivative of V0 and using $$\bar{z}_3^∗=-k_{\nu3}z_2^\frac{1}{p}$$ derive
```latex\begin{aligned}&\dot{V}_{0}=z_{2}(\bar{z}_{3}-\bar{z}_{3}^{∗})+z_{2}\bar{z}_{3}^{∗}\\&\leq|z_{2}||\bar{z}_{3}-\bar{z}_{3}^{∗}|-k_{v3}|z_{2}|^{1+\frac{1}{p}}.\end{aligned}```
From Lemma 2, for the term |$$\bar{z}_3-\bar{z}_3^∗$$|, it can be derived that
```latex|\bar{z}_3-\bar{z}_3^∗|\leq2^{1-\frac{1}{p}}|\bar{z}_3^p-\bar{z}_3^∗p|^{\frac{1}{p}}.```
Denote $$\overline{z}_{3}^{p}-\overline{z}_{3}^{∗p}=\overline{z}_{3}^{p}+k_{\nu3}^{p}z_{2}=\zeta$$, and then (13) can be written as
```latex\begin{aligned}&\dot{V}_{0}\leq-k_{v3}|z_{2}|^{1+\frac{1}{p}}+2^{1-\frac{1}{p}}|z_{2}||\bar{z}_{3}^{p}-\bar{z}_{3}^{*p}|^{\frac{1}{p}}\\&=-k_{v3}|z_{2}|^{1+\frac{1}{p}}+2^{1-\frac{1}{p}}|z_{2}||\zeta|^{\frac{1}{p}}.\end{aligned}```
For the term $$\left|z_2\right|\left|\zeta\right|^\frac{1}{p}$$, from Lemma 3, it can be obtained
```latex|z_2||\zeta|^{\frac{1}{p}}\leq\frac{p}{1+p}|z_2|^{1+\frac{1}{p}}+\frac{1}{1+p}|\zeta|^{1+\frac{1}{p}}.```
Then $$\dot{V}_0$$ can be derived as
```latex\dot{V}_0\leq-k_{v3}|z_2|^{1+\frac{1}{p}}+2^{1-\frac{1}{p}}\left(\frac{p}{1+p}|z_2|^{1+\frac{1}{p}}+\frac{1}{1+p}|\zeta|^{1+\frac{1}{p}}\right).```
Step 2: Taking the time derivative of V1 derives
```latex\dot{V}_{1}=\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\left(\bar{z}_{3}^{p}-\bar{z}_{3}^{∗p}\right)^{2-\frac{1}{p}}\dot{\bar{z}}_{3}+\frac{\bar{z}_{3}}{2^{1-\frac{1}{p}}}\int_{\bar{z}_{3}^{∗}}^{\bar{z}_{3}}\left(s^{p}-\bar{z}_{3}^{∗p}\right)^{1-\frac{1}{p}}ds.```
Invoking (14), it is noted that
```latex\frac{\bar{z}_3}{2^{1-\frac{1}{p}}}\int_{\bar{z}_3^∗}^{\bar{z}_3}\left(s^p-\bar{z}_3^{∗p}\right)^{1-\frac{1}{p}}ds\leq|\bar{z}_3||\bar{z}_3^p-\bar{z}_3^{∗p}|.```
Then, substituting $$\dot{\bar{z}}_3$$ in (11) into (16) derives
```latex\begin{aligned}&\dot{V}_{1}\leq|\bar{z}_{3}||\zeta|+\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\zeta^{2-\frac{1}{p}}\dot{\bar{z}}_{3}\\&=|\bar{z}_{3}||\zeta|+\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\zeta^{2-\frac{1}{p}}\left(-c^{2}z_{2}-k_{v1}\mathrm{tanh}^{\frac{2}{p}-1}(k_{v2}\zeta)\right),\end{aligned}```
where $$\overline{z}_3^p+k_{\nu3}^pz_2=\overline{z}_3^p-\overline{z}_3^{∗p}=\zeta$$ has been used. For the term $$|\bar{z}_3||\zeta|$$, from (14) and Lemma 3, it can be derived that
```latex\begin{aligned}&|\bar{z}_3||\zeta|\leq|\zeta||\bar{z}_3-\bar{z}_3^∗|+|\zeta||\bar{z}_3^∗|\\&\leq2^{1-\frac{1}{p}}|\zeta|^{1+\frac{1}{p}}+|k_{v3}\zeta||z_{2}|^{\frac{1}{p}}\\&\leq2^{1-\frac{1}{p}}|\zeta|^{1+\frac{1}{p}}+\frac{pk_{v3}^{1+\frac{1}{p}}}{1+p}|\zeta|^{1+\frac{1}{p}}+\frac{1}{1+p}|z_{2}|^{1+\frac{1}{p}}.\end{aligned}```
For the term $$-c^2\zeta^{2-\frac{1}{p}}z_2$$, from Lemma 3, it can be derived that
```latex\begin{aligned}&-c^2\zeta^{2-\frac{1}{p}}z_2\leq\left(c^{\frac{2p}{2p-1}}|z_2|^{\frac{2p-2}{2p-1}}|\zeta|\right)^{2-\frac{1}{p}}|z_2|^{\frac{2}{p}-1}\\&\leq\frac{2p-1}{1+p}\left(c^{\frac{2p}{2p-1}}|z_2|^{\frac{2p-2}{2p-1}}\right)^{1+\frac{1}{p}}|\zeta|^{1+\frac{1}{p}}+\frac{2-p}{1+p}|z_2|^{1+\frac{1}{p}}.\end{aligned}```
From (17) and (18) and the fact $$|\zeta|^{1+\frac{1}{p}}=|\zeta|^{\frac{p_{1}+p_{2}}{p_{1}}}=\zeta^{\frac{p_{1}+p_{2}}{p_{1}}}=\zeta^{1+\frac{1}{p}}$$, then $$\dot{V}_{1}$$ can be written as
```latex\begin{aligned}&\dot{V}_{1}\leq-\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\zeta^{2-\frac{1}{p}}\left(k_{v1}\mathrm{tanh}^{\frac{2}{p}-1}(k_{v2}\zeta)\right)\\&+\left(2^{1-\frac{1}{p}}+\frac{pk_{v3}^{1+\frac{1}{p}}}{1+p}+\frac{(2p-1)c^{\frac{2(1+p)}{2p-1}}|z_{2}|^{\frac{2\left(p^{2}-1\right)}{(2p-1)p}}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}(1+p)}\right)\zeta^{1+\frac{1}{p}}\\&+\left(\frac{1}{1+p}+\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^p}\frac{2-p}{1+p}\right)\left|z_2\right|^{1+\frac{1}{p}}.\end{aligned}```
Step 3: In view of Lemma 4, we consider
```latex|\zeta|=\left|\bar{z}_3^p+k_{v3}^pz_2\right|\leq r,```
where r is selected to satisfy $$r\leq k_{\nu3}^pE$$. In terms of (20), denote a domain containing the origin by
```latexD=\left\{(z_2,\bar{z}_3)\in\mathbb{R}^2\left|\left|\bar{z}_3^p+k_{v3}^pz_2\right|\leq r\right\}\subset\mathbb{R}^2.\right.```
Choose a set as
```latexB_r=\left\{(z_2,\bar{z}_3)\in\mathbb{R}^2||\bar{z}_3|^p+k_{v3}^p|z_2|\leq r\right\}.```
$$B_r\subset D$$ is obvious from the fact $$\left|\bar{z}_3^p+k_{v3}^pz_2\right|\leq|\bar{z}_3|^p+k_{v3}^p|z_2|\leq r$$. Let $$\alpha = \underset{\left|\overset{-}{z}_{3}\right|^{p} + k_{v 3}^{p} \left|z_{2}\right| = r}{\mathrm{min}} V > 0$$, where $$r\leq k_{\nu3}^pE$$ to ensure that $$|z_2|\leq E$$. Take $$\beta\in\begin{pmatrix}0,\alpha\end{pmatrix}$$ and let
```latex\Omega_\beta=\{(z_2,\bar{z}_3)\in B_r|V\leq\beta\}.```
Thus, it can be inferred that the set $$\Omega_{\beta}$$ is a compact set and $$\Omega_{\beta} \subset B_{r} \subset D$$. By Lemma 4, for the states in the domain $$D$$, the following inequality holds:
```latex\frac{\tanh(k_{v2}r)}{r}|\zeta|\leq|\tanh(k_{v2}\zeta)|.```
In addition, for the states in the domain $$B_{r}$$, it can be known that $$\left|z_{2}\right| \leq E$$. From (15), (19), (21), and $$\left|z_{2}\right| \leq E$$, the time derivative $$\overset{\cdot}{V}$$ can be written as
```latex\begin{aligned}\dot{V}&=\dot{V}_{0}+\dot{V}_{1}\\&\leq-k_{v3}|z_{2}|^{1+\frac{1}{p}}+2^{1-\frac{1}{p}}\left(\frac{p}{1+p}|z_{2}|^{1+\frac{1}{p}}+\frac{1}{1+p}|\zeta|^{1+\frac{1}{p}}\right)\\&-\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\zeta^{2-\frac{1}{p}}\left(k_{v1}\tanh^{\frac{2}{p}-1}(k_{v2}\zeta)\right)+\left(2^{1-\frac{1}{p}}+\frac{pk_{v3}^{1+\frac{1}{p}}}{1+p}+\frac{(2p-1)c^{\frac{2(1+p)}{2p-1}}|z_{2}|^{\frac{2(p^{2}-1)}{(2p-1)p}}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}(1+p)}\right)\zeta^{1+\frac{1}{p}}\\&+\left(\frac{1}{1+p}+\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\frac{2-p}{1+p}\right)\left|z_{2}\right|^{1+\frac{1}{p}}\\&-\left(\frac{k_{\nu1}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{\nu3}^p}\left(\frac{\tanh(k_{\nu2}r)}{r}\right)^{\frac{2}{p}-1}-\frac{2^{1-\frac{1}{p}}}{1+p}-2^{1-\frac{1}{p}}-\frac{pk_{\nu3}^{1+\frac{1}{p}}}{1+p}-\frac{(2p-1)c^{\frac{2(1+p)}{2p-1}}E^{\frac{2(p^2-1)}{(2p-1)p}}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{\nu3}^p(1+p)}\right)\zeta^{1+\frac{1}{p}}.\end{aligned}```
Denote
```latex\begin{aligned}&\varphi_{0}=k_{v3}-\frac{2^{1-\frac{1}{p}}p+1}{1+p}-\frac{1}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\frac{2-p}{1+p},\\&\varphi_{1}=\frac{k_{v1}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}}\left(\frac{\tanh(k_{v2}r)}{r}\right)^{\frac{2}{p}-1}-\frac{2^{1-\frac{1}{p}}}{1+p}-2^{1-\frac{1}{p}}-\frac{pk_{v3}^{1+\frac{1}{p}}}{1+p}\\&-\frac{(2p-1)c^{\frac{2(1+p)}{2p-1}}E^{\frac{2(p^2-1)}{(2p-1)p}}}{\left(2-\frac{1}{p}\right)2^{1-\frac{1}{p}}k_{v3}^{p}(1+p)}.\end{aligned}```
Invoking conditions (8) and (9), it can be known that $$\varphi_{0} , \varphi_{1} > 0$$ always holds. Recalling (12), it can be derived that
```latexV_1\leq\frac{1}{\left(2-\frac{1}{p}\right)k_{v3}^p}\zeta^2.```
In view of (23) and Lemma 5, the inequality (22) can be further represented as
```latex\begin{aligned}&\dot{V}\leq-2^{\frac{1+p}{2p}}\varphi_0V_0^{\frac{1+p}{2p}}-\left(\left(2-\frac{1}{p}\right)k_{v3}^p\right)^{\frac{1+p}{2p}}\varphi_1V_1^{\frac{1+p}{2p}}\\&\leq-\varphi\left(V_0^{\frac{1+p}{2p}}+V_1^{\frac{1+p}{2p}}\right)\\&\leq-\varphi V^{\frac{1+p}{2p}},\end{aligned}```
where
```latex\varphi=\min\left\{2^{\frac{1+p}{2p}}\varphi_0,\left(\left(2-\frac{1}{p}\right)k_{\nu3}^p\right)^{\frac{1+p}{2p}}\varphi_1\right\}.```
Invoking Lemma 1, it can be derived that the settling time of the subsystem $$\overline{S}_{2}$$ can be estimated as
```latexT_{\bar{s}_2}\leq\frac{V^{1-\frac{1+p}{2p}}\left(z_2(0),\bar{z}_3(0)\right)}{\varphi\left(1-\frac{1+p}{2p}\right)}.```
From the definition of $$\Omega_{\beta}$$, the relation $$\Omega_{\beta} \subset B_{r} \subset D$$ and (24), it can be known that $$\Omega_{\beta}$$ is a positively invariant set. Choose $$\Omega_{\beta}$$ as the compact set where the initial states are included. Then, from Lemma 1, it can be proved that $$\overline{S}_{2}$$ under (7) is finite-time stable at the origin. Due to the boundedness of tanh(·), the control input v satisfies (10). □ Remark 1. A systematic procedure for selecting the controller parameters in (7) and constructing the compact set $$\Omega_{\beta}$$ is presented as follows: (1) Set the constant value $$c$$. Choose a constant $$E$$ in $$\left|z_{2}\right| \leq E$$. It is obvious that $$E$$ represents a range in which the user desires the state $$z_{2}$$ locates. (2) Based on the actuator’s capacity, select a value for $$k_{v 1}$$ satisfying $$k_{v 1} > c^{2} E$$ such that the control action could be large enough to counteract the effect of $$- c^{2} z_{2}$$ in (5). (3) Choose appropriate parameters $$k_{v 3}$$ and $$p$$ satisfying (8). (4) Select a large enough $$k_{v 2}$$ and an appropriate $$r$$ satisfying (9) and $$r \leq k_{v 3}^{p} E$$. (5) Once $$k_{v 3} , p$$ and $$r$$ are selected, the set $$D$$ and $$B_{r}$$, in the proof of Lemma 6, are determined. Then, calculate $$\alpha$$ and choose $$\beta \in \left(0 , \alpha\right)$$. Finally, one can construct the compact set $$\Omega_{\beta}$$. After the subsystem $$\overline{S}_{2}$$ has been stabilized at the origin, select the switching time as
```latexT_1\left(z_2(0),\bar{z}_3(0)\right)=\frac{V^{1-\frac{1+p}{2p}}\left(z_2(0),\bar{z}_3(0)\right)}{\varphi\left(1-\frac{1+p}{2p}\right)}.```
When $$t>T_1(z_2(0),\bar{z}_3(0))$$, it is time to prove that the first-order subsystem S1, under the bounded controller $$\omega$$, is finite-time stable at the origin. Actually, this problem, i.e., the first-order finite-time stabilization problem considering input constraints, has been addressed in our previous result [23]. Thus, the following lemma are presented here and the detailed proof can be found in [23]. Lemma 7. For the first-order subsystem S1 in (4), design the control input as
```latex\omega=-k_{\omega1}\mathrm{tanh}^{\frac{m}{n}}(k_{\omega2}z_1).```
Then the subsystem $$S_{1}$$ with (26) is globally finite-time stable at the origin and the settling time function $$T_{s 1}$$ is bounded by
```latexT_{s1}(z_{10})\leq\frac{n\bar{r}^{\frac{m}{n}}V_{z1}^{\frac{n-m}{2n}}(z_{10})}{k_{\omega1}(n-m)\mathrm{tanh}^{\frac{m}{n}}(k_{\omega2}\bar{r})},```
where $$z_{10}$$ represents the initial value of $$z_{1}$$, and $$V_{z 1} = z_{1}^{2}$$, and $$k_{\omega 2}$$ and $$\overline{r}$$ are defined in Lemma 4. In addition, the control input is bounded by
```latex|\omega|\leq k_{\omega1}.```
Based on Lemma 6 and Lemma 7, the main results are presented in the following. Theorem 1. For the mobile robot system (2), under the transformation (3), if the initial states satisfy $$z_1\left(0\right)\in\mathbb{R}$$, $$\left(z_{2}(0),\bar{z}_{3}(0)\right)\in\Omega_{\beta}$$, where $$\overline{z}_3\begin{pmatrix}0\end{pmatrix}=cz_3\begin{pmatrix}0\end{pmatrix}$$, c is a non-zero constant and $$\Omega_{\beta}$$ is a given compact set, then (2) can be finite-time stabilized by the switching controllers (6) and (7) at the origin and the switching time is given as (25). The settling time function $$T_{2} \left(z_{0}\right)$$ of the closed-loop system (2), (6) and (7) is bounded by
```latexT_2(z_0)\leq T_1\left(z_2(0),\bar{z}_3(0)\right)+\frac{n\bar{r}^{\frac{m}{n}}V_{z1}^{\frac{n-m}{2n}}\left(z_1(T_1)\right)}{k_{\omega1}(n-m)\mathrm{tanh}^{\frac{m}{n}}(k_{\omega2}\bar{r})},```
where $$z_0=[z_1(T_1),z_2(0),\bar{z}_3(0)]^T$$. In addition, the control inputs are bounded by
```latex|\omega|\leq\max\{|c|,k_{\omega1}\},|v|\leq\frac{k_{v1}}{|c|}.```
Proof. First, following the same line as the proof of Lemma 6, it is proved that the subsystem $$\bar{S}_{2}$$ can be finite-time stabilized by the bounded controller (7) at the origin and the upper-bound estimate for the settling time of the subsystem $$\bar{S}_{2}$$, denoted as $$T_{1}$$, can be provided explicitly. Then, select $$T_{1}$$ as the switching time. For the subsystem $$S_{1}$$, since $$\overset{\cdot}{z}_{1} = c$$, when $$t \leq T_{1}$$, there is no finite-time escape phenomenon. Next, invoking Lemma 7 originating from [23], it can be proved that the subsystem $$S_{1}$$ can be globally finite-time stabilized by (6) at the origin. It is noted that $$z_{10} = z_{1} \left(T_{1}\right) = z_{1} \left(0\right) + c T_{1}$$ and thus the settling time function $$T_{2}$$ is bounded by (27). Under the transformation (3), it is obvious that $$z_{1} = z_{2} = z_{3} = 0 \Rightarrow x = y = \theta = 0$$. In addition, it can be easily checked that the control inputs are bounded by (28). □ To help understand how to implement the control algorithm in a NWMR, the pseudocode and schematic are provided below in Algorithm 1 and Figure 1, respectively:

Algorithm 1: Finite-Time Stabilization of a NWMR With Input Constraints

1

Input the initial states $$x \left(0\right) , y \left(0\right) , \theta \left(0\right)$$;

2

Based on the guideline in Remark 1 and the actuator’s capacity, select a set of controller parameters related to controller (7): $$c,E,k_{\nu1},k_{\nu3},p,k_{\nu2},r$$ and construct the given compact set $$\Omega_{\beta}$$;

3

From $$x (0) , y (0) , \theta (0) , c ,$$ Equation (3) and $$\overset{-}{z}_{3} = c z_{3}$$, compute $$z_{2} (0)$$ and $$\overset{-}{z}_{3} (0)$$;

4

If $$\left(z_{2} (0) , \overset{-}{z}_{3} (0)\right) \in \Omega_{\beta}$$, go to 5. Otherwise, go back to 2 to adjust the controller parameters;

5

From (25), compute the switching time $$T_{1}$$, and then select the controller parameters $$k_{\omega1},k_{\omega2},m,n,\overset{-}{r}$$;

6

From (27), compute the explicit upper-bound estimate for the closed-loop settling time;

7

Output the controller $$v , \omega$$ and the closed-loop settling time $$T_{2}$$.

Figure 1. The schematic of the proposed algorithm.
From Lemma 6, Lemma 7, and Theorem 1, it is evident that the finite-time stabilization of a NWMR with input constraints involves two steps. In the first step, based on a constant controller $$\omega = c$$, one needs to design a bounded finite-time controller $$v$$ to stabilize the second-order subsystem $$\bar{S}_{2}$$ in (5), while also computing the upper-bound estimate of the settling time for subsystem $$\bar{S}_{2}$$ under controller (7). In the second step, the upper-bound estimate obtained in the first step is selected as the switching time. Subsequently, a bounded finite-time controller $$\omega$$ is designed to stabilize the first-order system $$S_{1}$$ in (4). Through these two steps, an explicit upper-bound estimate of the closed-loop settling time is derived, as shown in (27). Furthermore, due to the boundedness of tanh(·), the controllers $$\omega , v$$ can be bounded by the controller parameters, as described in (28), which are selected in advance based on the actuator’s capacity. It is important to note that the upper-bound estimate of the settling time for the subsystem $$\bar{S}_{2}$$ under (7) is relatively conservative. As a result, the switching time is also conservative. Therefore, reducing this conservatism will be a key focus of future research.

4. Simulation Results

In this section, for the proposed bounded finite-time switching controllers (6) and (7), the simulation results are presented to validate their effectiveness. In addition, our method is compared with the switching method in [18], which has been wildly used in recent studies on nonholonomic systems [25,26,27,28]. For a fair comparison, the initial conditions are set as $$\left[x \left(0\right) , y \left(0\right) , \theta \left(0\right)\right]^{T} = \left[3 , - 1 , - 0 . 5\right]^{T}$$. Here, for our controllers (6) and (7), we choose the controller parameters as $$c=0.5,E=1,k_{\omega1}=3,k_{\omega2}=1,m=1,n=3,\bar{r}=620$$, $$k_{v1}=7,k_{v2}=10,k_{v3}=1.245,p=27/17,r=1$$. From (25), the switching time is calculated as $$T_1\left(z_2(0),\bar{z}_3(0)\right)=1236.1s$$. In addition, from Step 3 in the proof of Lemma 6, it can be obtained that $$\alpha = 0.0343$$. Then we choose $$\beta = 0.034$$. It can be checked that the initial conditions are included in the given compact set $$\Omega_{\beta}$$. The responses of the states $$x, y$$ and $$\theta$$ are shown in Figure 2a,b, respectively. The control inputs $$\omega$$ and $$v$$ are shown in Figure 3.
Figure 2. State trajectories of the mobile robot using the proposed algorithm. (<b>a</b>) State trajectories <i>x</i> and <i>y</i>; (<b>b</b>) State trajectories $$\theta$$.
From Figure 2a,b, it can be seen that the state trajectories $$x , y$$ converge to zero at about $$2 s$$ and $$\theta$$ converges to zero at about $$1440 s$$. It is evident that all states converge to zero in finite time. In addition, from (27), it can be computed that the upper-bound estimate of the closed-loop settling time is $$1545 . 2828 s$$. It should be noted that, before the switching time $$T_{1} \left(z_{0}\right)$$, since $$\omega$$ is a non-zero constant $$c$$, the value of $$\left|\theta\right|$$ increases continuously. When $$t > T_{1} \left(z_{0}\right)$$, the angle $$\theta$$ is stabilized by the controller (6) in a finite time. From Figure 3, invoking (28), it can be seen that the control inputs are bounded by the prescribed parameters, i.e., $$\left|\omega\right| \leq 3$$ rad/s and $$\left|v\right| \leq 14$$ m/s.
Figure 3. Stabilization inputs $$\omega$$ and <i>v</i>.
Using the switching method in [18], for the nonholonomic system (2), introduce the following state transformation and input transformation:
```latex\left.\left\{\begin{array}{c}z_{1D}=\theta_D,\\z_{2D}=x_D\sin\theta_D-y_D\cos\theta_D,\\z_{3D}=x_D\cos\theta_D+y_D\sin\theta_D,\\v_{1D}=\omega_D,\\v_{2D}=v_D-z_{2D}\omega_D,\end{array}\right.\right.```
where the subscript D denotes the corresponding variables from Defoort’s paper [18], distinguishing them from the related variables used in this paper. Under this transformation (29), the mobile robot system (2) can be written as
```latexS_{1D}{:}\dot{z}_{1D}=v_{1D},S_{2D}{:}\dot{z}_{2D}=v_{1D}z_{3D},\dot{z}_{3D}=v_{2D}.```
Design a controller as
```latex\begin{aligned}&v_{1D}=\begin{cases}\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,c,t\leq T_{1D},\\-\alpha_3|z_{1D}|^2\mathrm{sign}(z_{1D})-\beta_3\mathrm{sign}(z_{1D}),t>T_{1D},&&\end{cases}\\&v_{2D}=\begin{cases}-\frac{\alpha_1+3\beta_1z_{2D}^2}{2c}\operatorname{sign}(s)-\frac{|\alpha_2s+\beta_2|s|^3\operatorname{sign}(s)|^\frac{1}{2}\operatorname{sign}(\alpha_2s+\beta_2|s|^3\operatorname{sign}(s))}{c},&t\leq T_{1D},\\\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,\,0,t> T_{1D},&&\end{cases}\end{aligned}```
where $$s=cz_{3D}+||cz_{3D}|^2\mathrm{sign}(cz_{3D})+\alpha_1z_{2D}+\beta_1|z_{2D}|^3\mathrm{sign}(z_{2D})|^\frac{1}{2}\mathrm{sign}\left(|cz_{3D}|^2\mathrm{sign}(cz_{3D})+\alpha_1z_{2D}+\beta_1|z_{2D}|^3\mathrm{sign}(z_{2D})\right)$$. We choose the controller parameters as $$c = 0 . 5 , \alpha_{1} = \beta_{1} = 0 . 1 , \alpha_{2} = \beta_{2} = 1 . 5 , \alpha_{3} = \beta_{3} = 0 . 1$$, and then from $$T_{1 D} = \frac{2}{\sqrt{\alpha_{2}}} + \frac{2}{\sqrt{\beta_{2}}} + \frac{2 \sqrt{2}}{\sqrt{\alpha_{1}}} + \frac{2 \sqrt{2}}{\sqrt{\beta_{1}}}$$, the switching time is computed as $$T_{1 D} = 21 . 1545 s$$. The responses of the states $$x_{D} , y_{D}$$ and $$\theta_{D}$$ are shown in Figure 4a,b, respectively. The control inputs $$\omega_{D}$$ and $$v_{D}$$ are shown in Figure 5.
Figure 4. State trajectories of the mobile robot using the method in [18]. (<b>a</b>) State trajectories $$x_D$$ and $$y_D$$; (<b>b</b>) State trajectories $$\theta_D$$.
Figure 5. Stabilization inputs $$\omega_D$$ and $$v_D$$.
By comparing Figure 3 and Figure 5, it is evident that the control inputs $$\omega$$ and $$v$$, in this paper, are bounded by (28), with computed values $$\left|\omega\right| \leq 3$$ rad/s and $$\left|v\right| \leq 14$$ m/s. However, for the same initial states, since the method in [18] does not account for input constraints, the input signals $$\omega_{D}$$ and $$v_{D}$$ exceed $$3$$ rad/s and $$14$$ m/s, respectively, and even $$\left|\omega_{D}\right| > 10$$ rad/s. In practice, excessive input values that exceed the actuator’s limits can potentially damage the actuator. Therefore, it is crucial to consider input constraints in control design. It is important to note that this comparison focuses solely on whether the constraint level of the control input can be pre-defined by parameters, and whether the designed controllers can be bounded by these predefined values. Thus, the comparison above demonstrates the validity of the proposed algorithm in this paper. In addition, from Figure 2b and Figure 4b, it is clear that the switching time in this paper is significantly longer than that in [18]. Consequently, achieving finite-time stabilization of a NWMR with input constraints while reducing the switching time remains a key direction for our future research.

5. Conclusions

In this paper, the finite-time stabilization control problem for a nonholonomic wheeled mobile robot with input constraints is addressed. Specifically, based on the bounded odd function tanh(·), bounded finite-time controllers are proposed to stabilize the mobile robot. In addition, the upper-bound estimate for the closed-loop settling time is provided explicitly and the level of input constraint is presented by parameters dependent on actuator’s capacity. Using appropriate Lyapunov functions, a comprehensive finite-time stability analysis is conducted. Moreover, a guideline is provided to clarify how to choose the controller parameters and construct a compact set contained in the domain of attraction. Finally, the effectiveness of the proposed controllers is validated through simulation results.

Acknowledgments

The authors wish to thank the associate editor and reviewers who made some valuable suggestions to improve the paper significantly.

Author Contributions

Conceptualization and Methodology, X.L. and Z.Z.; Software and Validation, X.L.; Writing—Original Draft Preparation, X.L.; Writing—Review & Editing, Z.Z.

Ethics Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 62473020).

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

1.
Zheng Y, Zheng J, Shao K, Zhao H, Man Z, Sun Z. Adaptive fuzzy sliding mode control of uncertain nonholonomic wheeled mobile robot with external disturbance and actuator saturation. Inf. Sci. 2024, 663, 120303. [Google Scholar]
2.
Li J-W. Adaptive tracking and stabilization of nonholonomic mobile robots with input saturation. IEEE Trans. Autom. Control. 2022, 67, 6173–6179. [Google Scholar]
3.
Wang Z, Li G, Chen X, Zhang H, Chen Q. Simultaneous stabilization and tracking of nonholonomic WMRs with input constraints: Controller design and experimental validation. IEEE Trans. Ind. Electron. 2019, 66, 5343–5352. [Google Scholar]
4.
Wang Y, Miao Z, Zhong H, Pan Q. Simultaneous stabilization and tracking of nonholonomic mobile robots: A lyapunov-based approach. IEEE Trans. Control. Syst. Technol. 2015, 23, 1440–1450. [Google Scholar]
5.
Urakubo T. Feedback stabilization of a nonholonomic system with potential fields: application to a two-wheeled mobile robot among obstacles. Nonlinear Dyn. 2015, 81, 1475–1487. [Google Scholar]
6.
Guo Z, Wang Z, Ding S, Li S, Wang X. Global finite-time stabilization of spacecraft attitude with disturbances via continuous nonlinear controller. IEEE Trans. Aerosp. Electron. Syst. 2023, 59, 2608–2620. [Google Scholar]
7.
Feng Y, Yu X, Han F. On nonsingular terminal sliding-mode control of nonlinear systems. Automatica 2013, 49, 1715–1722. [Google Scholar]
8.
Feng Y, Yu X, Man Z. Non-singular terminal sliding mode control of rigid manipulators. Automatica 2002, 38, 2159–2167. [Google Scholar]
9.
Zhou Y, Ríos H, Mera M, Polyakov A, Zheng G, Dzul A. Homogeneity-based control strategy for trajectory tracking in perturbed unicycle mobile robots. IEEE Trans. Control. Syst. Technol. 2024, 32, 274–281. [Google Scholar]
10.
Bhat SP, Bernstein DS. Geometric homogeneity with applications to finite-time stability. Math.Control. Signals Syst. 2005, 17, 101–127. [Google Scholar]
11.
Li Y-X. Finite time command filtered adaptive fault tolerant control for a class of uncertain nonlinear systems. Automatica 2019, 106, 117–123. [Google Scholar]
12.
Chen Q, Ye Y, Hu Z, Na J, Wang S. Finite-time approximation-free attitude control of quadrotors: Theory and experiments. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 1780–1792. [Google Scholar]
13.
Li X, Tang J, Lin K, Zuo Z. Semi-global finite-time stabilization of second-order systems with bounded control input. IEEE Trans. Ind. Electron. 2024, in press.
14.
Cong Y, Du H, Zhu W, Chen C-C. Position regulation of industrial robots via bounded integral terminal sliding mode control algorithm.  IEEE Trans. Ind. Electron. 2024, 71, 7403–7412. [Google Scholar]
15.
Kimura S, Nakai T, Nakamura H, Ibuki T, Sampei M. Finite-time control of two-wheeled mobile robot via generalized homogeneous locally semiconcave control lyapunov function. In Proceedings of the 55th Annual Conference of the Society of Instrument and Control Engineers of Japan, Tsukuba, Japan, 20–23 September 2016.
16.
Defoort M, Djemaȉ M. Finite-time controller design for nonholonomic mobile robot using the heisenberg form. In Proceedings of the IEEE International Conference on Control Applications, Yokohama, Japan, 8–10 September 2010.
17.
Zhu X, Dong G, Hu D. A finite-time controller for a class of wheeled mobile robots. In Proceedings of the 6th World Congress on Intelligent Control and Automation, Dalian, China, 21–23 June 2006.
18.
Defoort M, Demesure G, Zuo Z, Polyakov A, Djemaȉ M. Fixed-time stabilisation and consensus of non-holonomic systems. IET Control. Theory Appl. 2016, 10, 2497–2505. [Google Scholar]
19.
Bhat SP, Bernstein DS. Finite-time stability of continuous autonomous systems. SIAM J. Control. Optim. 2000, 38, 751–766. [Google Scholar]
20.
Huang X, Lin W, Yang B. Global finite-time stabilization of a class of uncertain nonlinear systems. Automatica 2005, 41, 881–888. [Google Scholar]
21.
Qian C, Lin W. Non-lipschitz continuous stabilizers for nonlinear systems with uncontrollable unstable linearization. Syst. Control. Lett. 2001, 42, 185–200. [Google Scholar]
22.
Deng Z, Xu Y, Sun H, Shen X. Distributed, bounded and finite-time convergence secondary frequency control in an autonomous microgrid. IEEE Trans. Smart Grid 2019, 10, 2776– 2788. [Google Scholar]
23.
Zuo Z, Li X, Ning B, Han Q-L. Global finite-time stabilization of first-order systems with bounded controls. IEEE Trans. Circuits Syst. II Express Briefs 2023, 70, 2440–2444. [Google Scholar]
24.
Hardy G, Littlewood J, Polya G. Inequalities; Cambridge Univ. Press.: Cambridge, UK, 1952.
25.
Sanchez-Torres J, Defoort M, Munoz-Vazquez A. Predefined-time stabilisation of a class of nonholonomic systems. Int. J. Control. 2020, 93, 2941–2948. [Google Scholar]
26.
Boubaker S, Labbadi M, Alsubaei F. A Generalized Supertwisting Algorithm Based on Terminal Sliding Manifold Control of Perturbed Unicycle Mobile Robots. IEEE Access 2024, 12, 26324–26333. [Google Scholar]
27.
Labbadi M, Boubaker S, Djemai M, Mekni S, Bekrar A. Fixed-Time Fractional-Order Global Sliding Mode Control for Nonholonomic Mobile Robot Systems under External Disturbances. Fractal Fract. 2022, 6, 177. [Google Scholar]
28.
Cai Y, Shang Y, Gao F, Huang J, Wu Y. Prescribed-Time Stabilization for a Class of Switched Nonholonomic Systems. In Proceedings of 2nd Conference on Fully Actuated System Theory and Applications, Qingdao, China, 14–16 July 2023.
TOP