Robust Output Feedback Stabilization and Tracking for an Uncertain Nonholonomic Systems with Application to a Mobile Robot

This paper presents a novel robust output feedback control that simultaneously performs both stabilization and trajectory tracking for a class of underactuated nonholonomic systems despite model uncertainties, external disturbance, and the absence of velocity measurement. To solve this challenging problem, a generalized normal form has been successfully created by employing an input–output feedback linearization approach and a change in coordinates (diffeomorphism). This research mainly focuses on the stabilization problem of nonholonomic systems that can be transformed to a normal form and pose several challenges, including (i) a nontriangular normal form, (ii) the internal dynamics of the system are non-affine in control, and (iii) the zero dynamics of the system are not in minimum phase. The proposed scheme utilizes combined backstepping and sliding mode control (SMC) techniques. Furthermore, the full-order high gain observer (HGO) has been developed to estimate the derivative of output functions and internal dynamics. Then, full-order HGO and the backstepping SMC have been integrated to synthesize a robust output feedback controller. A differential-drive type (2,0) the wheeled mobile robot has been considered as an example to support the theoretical results. The simulation results demonstrate that the backstepping SMC exhibits robustness against bounded uncertainties.


Introduction
In recent years, nonholonomic systems have been considered a challenging benchmark problem due to their restricted mobility and underactuated behavior.These systems satisfy non-integrable velocity constraint, which originates as a result of physical kinematic constraints (such as pure rolling conditions) imposed on the system.This classical subject has gained significant attention due to the prevalence of nonholonomic constraints in numerous advanced robotic structures, including space manipulators, legged robots, wheeled mobile robots (WMRs), and multi-fingered robot hands.WMRs are a prime example of nonholonomic dynamic systems and are commonly used as a benchmark problem for this class of systems due to their limited mobility and underactuated behavior.The two main control tasks of WMRs are tracking to a reference trajectory and posture stabilization at the desired pose [1].From a control perspective, posture stabilization of a mobile robot is a more challenging task than trajectory tracking.This is because posture stabilization is considered an input-state problem, with two control inputs responsible for zeroing three independent configuration errors (e x , e y , e θ ).On the other hand, trajectory tracking can be recognized as an input-output problem, with two control inputs available to zero two-dimension errors (e x , e y ).In fact, it is understood that these systems cannot achieve asymptotic stabilization via smooth time-invariant feedback control [2].
In the literature, many researchers have proposed solutions for trajectory tracking and posture stabilization using kinematic models [3][4][5][6][7].However, few research papers have focused on integrating mobile robot kinematic and dynamic models without considering the parameter uncertainties and external disturbances [8][9][10].Input-output feedback linearization has been widely utilized to develop feedback controllers for WMRs due to its challenging nonlinear model [1,5,7,[11][12][13].These papers use the exact model of WMR for the cancellation of nonlinear terms that may cause inaccuracy during trajectory tracking and posture stabilization.However, in practical scenarios, the exact knowledge of kinematic and dynamic model parameters is not known due to uncertainties.Therefore, robust/adaptive-based tracking and stabilization controllers have been proposed to mitigate the effect of external disturbances and model uncertainties [14][15][16][17][18].In particular, SMC is preferred because of its robustness against parameter uncertainties and external disturbance, satisfactory transient behavior, and fast response [19][20][21][22][23][24].It is also notable that most of these papers have studied stabilization and tracking problems separately.
Another major issue found in a mobile robot is the measurement of velocity components and other states through sensors.However, some states are difficult to measure due to the cost, weight limitation, and sensor unavailability.Furthermore, in practical situations, velocity measurement using a sensor may degrade the desired performance due to the measurement noise.For that reason, output feedback control is the most effective solution to avoid sensor-based measurement of velocity components.In [3,15,[25][26][27][28][29][30], trajectory tracking and stabilization [31][32][33] based on output feedback control has been developed to estimate the velocities of a mobile robot.Among these approaches, HGO is the most popular approach used to estimate the derivative of the output signal while eliminating the disturbance effects [34,35].Bocker and Khalil in [36] and other researchers in [37] designed full-order observers containing an extended Kalman filter observer and an extended HGO for the estimation of internal and external dynamics, respectively.These techniques were applied to the underactuated systems that can be transformed into strict feedback normal form structure.However, the key difficulties in developing an observer-based output feedback control for nonholonomic systems are that (i) the Coriolis matrix for Lagrangebased systems possesses quadratic cross terms of unmeasured velocities, (ii) nonholonomic constraints are imposed on the system, and (iii) the inability to transform into a standard strict feedback normal form.Therefore, these existing techniques cannot be directly applied to nonholonomic WMRs.
This paper presents a novel robust output feedback control approach for a class of underactuated nonholonomic systems with multi-inputs and multi-outputs (MIMO).The control technique is based on the globally defined normal form representation of the kinematic and dynamic models.Moreover, the proposed method takes into account bounded model uncertainties and external disturbances, and the absence of velocity measurements in the nonholonomic system.The controller developed in this study is applicable to both stabilization and trajectory tracking problems simultaneously.
To address the challenging task of stabilizing nonholonomic systems, a generalized normal form is introduced.This normal form representation aims to overcome the inherent difficulties associated with stabilizing such systems.By using the principles of differential geometry, it becomes feasible to change the system dynamics into a normal form through a skillfully chosen change in coordinates, known as a diffeomorphism, along with the implementation of an input-output linearization method.This transformation allows for the separation of the nonholonomic system dynamics into external and internal dynamics.The proposed normal form for uncertain nonholonomic systems integrates both kinematic and dynamic models uncertainties and external disturbance.The effect of these uncertainties is added as a perturbation term.This perturbation term enters the system through the same channel as the control input, i.e., matched uncertainty.In comparison with previous research, [1,5,7,[11][12][13][14]38], the nominal values of system parameters were used for feedback linearization, and the systems were not designed to accommodate such uncertainties.
The primary effort of this research is on the stabilization problem of the normal form for the class of nonholonomic dynamic systems.This presents several challenges that need to be overcome.Firstly, the normal form representation for nonholonomic systems exhibits a nontriangular structure.Secondly, the internal dynamics of the system are non-affine in control, meaning they are highly nonlinear.Lastly, the zero dynamics of the system are not minimum phase, which further complicates the stabilization problem.
The primary contributions of this paper are as follows: 1.

Introduction of a Generalized Normal Form:
We propose a novel normal form representation for underactuated nonholonomic systems that can handle both kinematic and dynamic model uncertainties as well as external disturbances.

2.
Novel Recursive Backstepping Sliding Mode Control (BSMC): To tackle the challenges of stabilizing nonholonomic systems, we develop a recursive backstepping SMC technique that ensures asymptotic stabilization of both external and internal dynamics.This method provides a unified solution for both trajectory tracking and stabilization problems while being robust to bounded parameter uncertainties and external disturbances.

3.
Output Feedback Control with High Gain Observer (HGO): We design an output feedback control that employs a full-order HGO to estimate the derivative of the output function and velocity components without sensors, even in the presence of bounded uncertainties.This method achieves comparable performance to full-state feedback control.

4.
Application to Wheeled Mobile Robots: The theoretical development is applied to the case of wheeled mobile robots, demonstrating the broad applicability of the proposed control approach to a typical benchmark problem in nonholonomic dynamic systems.

Kinematic and Dynamic Modeling for an Uncertain Nonholonomic System
Consider a class of nonholonomic dynamic system with n generalized coordinates q = [q 1 . . .q n ] T , q ∈ R n , subject to m nonholonomic constraints (m < n) is described as: let J(q) ∈ R m * n be a full rank matrix linked with kinematic constraints.Suppose [S 1 (q) . . .S n−m (q)] are linearly independent vector fields in the null space of J(q): J(q)S(q) = 0 (2) using Equations ( 1) and (2) generated by the nonholonomic constraints, there exists a velocity vector ϑ(t) ∈ R n−m q = S(q)ϑ(t) Moreover, the velocity vector ϑ(t) in the kinematic model (3) can be formulated in terms of systems dynamics, such as where Q ∈ R (n−m) * (n−m) is a matrix of constant parameters, and ϖ(t) ∈ R n−m is a vector of angular wheel velocities.The kinematic model of the nonholonomic system under parameter uncertainties is derived by substituting Equation (4) into Equation (3): where C(q) = S(q)Q.The dynamic model of a nonholonomic system under parameter variations and external disturbance can be expressed using the Euler-Lagrange equation: where M(q) ∈ R n * n is a positive definite and symmetric inertia matrix, B(q) ∈ R n * (n−m) defines the input transformation matrix, V(q, q) ∈ R n * n denotes the Coriolis and centripetal forces, τ ∈ R n−m is the input torque vector, τ d ∈ R n−m represents the bounded external disturbance, and λ ∈ R m is constrain forces vector.The parameter uncertainties in both kinematic and dynamic models arise due to the incorrect measurement of parameters, such as moment of inertia and mass.The actual matrix Q, M(q), and V(q, q) can be defined as where Q 0 , M 0 (q), and V 0 (q, q) represent nominal functions, and ∆Q, ∆M(q), and ∆V(q, q) denote internal uncertainties because of parametric and nonparametric uncertainties.
The kinematic (5) and dynamic models (6) can be integrated to develop the state-space model for an uncertain nonholonomic system to enhance the performance of stabilizing and tracking controllers.Taking the time derivative of (5) yields q = Ċ(q)ϖ(t) + C(q) π(t) (7) substituting Equations ( 5) and ( 7) into Equation ( 6), then multiplying the resultant equation by C T and considering C T J T λ = 0 because of Equation ( 2), one obtains after simplification of Equation ( 8) for π, which gives where N(q) = C T MC, F(q, q) = C T M Ċ + C T VC, and B = C T B. The state-space model is determined by combining the kinematic model (5) and dynamic Equation (9), yields (10) where x ∈ R 2n−m is the state vector; τ ∈ R n−m is the control input; and f (x), p(x), g(x), and δ(x, t) are sufficiently smooth functions in domain D ∈ R 2n−m .In Equation (10), p(x) represents the parameter uncertainties in the dynamic model, while δ(x, t) accounts for the effects of external disturbances.Table 1 describes the symbols and nomenclature used in the research paper, providing a clear reference for all variables and parameters.Assumption 2. Assume the input transformation matrix g(x) has full rank (n − m), with g(0) ̸ = 0 to hold the controllability of the system.Furthermore, we assume that f (x) is known, while p(x), g(x), and δ(x, t) could be uncertain.
Remark 2. The system (10) is not input-state feedback linearizable due to nonholonomic constraints.However, input-output feedback linearization can be achieved by selecting the suitable output functions [13].substituting x = 0 in input transformation matrix g(x)

Input-Output Feedback Linearization
The input-output feedback linearization technique has been used to change the dynamics of an uncertain nonholonomic system (10) into a generalized normal form.For this purpose, a suitable set of output functions and internal dynamics variables are chosen skillfully.Then, trajectory tracking and stabilization problems can be solved using the same output functions under model uncertainties and external disturbance.Furthermore, it also provides ease in the development of HGO design for the estimation of internal dynamics and derivative of output functions.The proposed output functions can be described as Proposition 1.Consider that the multi-input multi-output nonholonomic system (10) and ( 11) is partially input-output feedback linearizable.Furthermore, the largest feedback linearizable subsystem of the system (10) has 2(n − m) dimension in x ∈ R 2n−m and the relative degree ρ = {ρ 1 , ρ 2 , . . ., ρ n−m } of each subsystem is ρ i = 2 for each output.
Proof.Taking the Lie derivative of Equation (11), yields where Φ i (x) = (∂h i (q)/∂q)C(q).Equation ( 12) is not explicitly dependent on the input function.Therefore, another derivative of output function (12) can be computed until input appears ÿi substituting π from Equation (10) into Equation ( 13), gives where The decoupling matrix β(x) is described as The decoupling matrix will be nonsingular if the regular conditions are satisfied in Equation (11).As a result, system (10) is partially input-output feedback linearizable by the proposed output.
The relative degree of a nonholonomic system (10) is 2(n − m) in x ∈ R 2n−m , with each output having a relative degree of two.We require m more components to achieve a diffeomorphic transformation.Therefore, we define the following change in coordinates . . .
The variables η 1 to η m and ξ 1,1 to ξ n−m,2 define the internal and external dynamics of a nonholonomic system (10), respectively.The internal dynamics η is chosen to ensure that T(x) in ( 15) is a legitimate diffeomorphism on a domain D 0 ⊂ D in R (2n−m) and satisfy the following requirements: The suitable choice of internal dynamics ψ(q) that not only satisfies the conditions in ( 16) but also provides an ease to attain asymptotic stabilization of η dynamics, would be If the Jacobian matrix ∂T ∂x of the map T(x) is invertible at a point x 0 ∈ D, then Equation ( 15) will be a valid diffeomorphism.
Assuming that (∂T/∂x) has full rank for all x ∈ R (2n−m) , the map T(x) is globally diffeomorphism.Therefore, one can apply the coordinates transformation (15) to change the dynamics of the system (10) into a normal form representation: where ᾱ(η, where ξ ∈ R 2(n−m) ; η ∈ R m ; and A, B, and C can be defined as The following nonlinear feedback control is applied to the external dynamics (19b) of the system for obtaining input-output feedback linearization and decoupling where u is the auxiliary control input, and β0 (η, ξ) and Γ0 (η, ξ) are the nominal values of β(η, ξ) and Γ(η, ξ), respectively.Due to the parameter uncertainties, these values are replaced by their estimated values in the feedback control law (20).Through this process, the control action will be divided into continuous and switching components, to minimize the amplitude of the switching component to mitigate the chattering effect.After, substituting control law ( 20) into (19b), we obtained where perturbation term ∆ d (η, ξ, u, t) can be defined as We assume that perturbation term ∆ d (η, ξ, u, t) satisfies the bound Thus, to satisfy Equation ( 22), we need the inequalities to hold over a domain that contains the origin for some continuous function σ(η, ξ).Inequality (23a) is restrictive because it puts a definite limit on the perturbation ( β β−1 0 − I) to achieve asymptotic stabilization of equilibrium point at the origin.Inequality (23b), conversely, is not restrictive, because we do not require σ to be small.Substituting Equation (22) into Equation (21) gives using the results of Equation ( 24), we consider a generalized nontriangular normal form (19a)-(19c) for the class of nonholonomic system with (n − m) = 2 inputs that can be transformed to the following form: The generalization of obtained results to the case of (n − m) inputs is trivial.If we set ξ 1,2 = ξ 2,1 = 0, the origin of η = 0 is not asymptotically stable for the zero dynamics of η = F (η, 0)0 Thus, the zero dynamics in (26) is not minimum phase.Remark 3. In contrast to nonlinear systems in a strict feedback form, this paper primarily focuses on the stabilization and tracking issues of nonholonomic systems represented in normal form (25) and (26), which pose the following difficulties: 1.
The proposed normal form has a nontriangular structure.

2.
Internal dynamics of the system is non-affine in control, i.e., highly nonlinear.

3.
The zero dynamics of the system is not minimum phase.
where tracking errors are described as and error dynamics of trajectory tracking can be determined by computing the time derivative of (28), which gives Step 1: Backstepping SMC for Subsystem Σ 1 Assumption 4. Suppose there exists a virtual state feedback control is asymptotically stable, i.e., ξ 1,1e → zero as t → ∞.Furthermore, there is a Lyapunov function ) that is positive definite and satisfies the equation To design a backstepping SMC, we start by designing the sliding manifold, s 1 = ξ 1,2e − φ 1 (ξ 1,1e ) = 0 such that, when the motion is restricted to a manifold, Assumption 4 applies.We then apply change in variables s 1 = ξ 1,2e − φ 1 (ξ 1,1e ) to transform the dynamics of system (29a) and (29b) into an equivalent representation where φ1 (ξ 1,1e ) = ∂φ 1 ∂ξ 1,1e (φ 1 + s 1 ).Later, u 1 will be formulated to bring s 1 to zero in finite time and remain there for all future time.Theorem 1.Consider the subsystem (31a).Suppose Assumption 4 holds and the uncertainty satisfies the inequality (22).Let the following backstepping SMC be defined by u 1 where µ 1 is a positive constant and Suppose µ 1 is chosen such that the error dynamics of subsystem (31a) with arbitrary initial conditions are bounded and exponentially converge to zero for all t ≥ 0 and reach the set Proof.Consider the Lyapunov function Taking V2 (ξ 1,1e , s 1 ) along the trajectories of (31a), yields To reduce the chattering effect, a high-slope saturation function has been used in place of the signum function as substituting control law u 1 (32) and (33) in the above equation, gives Therefore, every trajectory reaches the manifold (surface) One also observes that whenever |s 1 (0)| > µ 1 , |s 1 (t)| will reach the set {|s 1 | ≤ µ 1 } in finite time and maintain inside for all future time.
Step 2: Backstepping SMC for Subsystem Σ 2 Assumption 5.The internal dynamics η e in (31b) is asymptotically stable using the virtual state feedback control Furthermore, the internal dynamics in (31b) is input-to-state stable with respect to ξ 1,1e and s 1 , a positive definite Lyapunov candidate V 3 (η e ) exists that satisfies the inequality where b 1 , b 2 , and k 3 are class K functions.
Proof.Consider the Lyapunov function Taking V4 along the trajectories of subsystem ( 35) and substituting virtual feedback control (36) in the above equation Thus, V4 (η e , z) is negative definite.Therefore, η e and z will remain bounded and converge to the origin as t → ∞.
To design a robust backstepping SMC for the entire closed-loop subsystem (35), we start by designing the sliding manifold (surface) s 2 = ξ 2,2e − φ 3 (η e , z) = 0 such that, when the motion is restricted to manifold, Proposition 2 applies.We then apply change in variables s 2 = ξ 2,2e − φ 3 (η e , z) that transforms the dynamics of subsystem (35) into an equivalent representation where φ3 (η e , z) = ∂φ 3 ∂η e ηe + ∂φ 3 ∂z ż.Now, u 2 will be developed to bring s 2 to zero in finite time and remain there for all future time.Theorem 2. Consider the closed-loop subsystem (37).Suppose Proposition 2 is satisfied and the uncertainty meets the inequality (22).Let the following backstepping SMC be defined by u 2 where µ 2 is a positive constant and Furthermore, suppose µ 2 is chosen such that the error dynamics of subsystem (37) with arbitrary initial conditions are bounded and exponentially converge to the origin ∀ t ≥ 0 and reach the set {|s 2 | ≤ µ 2 } in finite time.

Backstepping SMC for Stabilization
Consider the nonholonomic system (25) and (26).After applying the same virtual state feedback control developed in (30), (34), and (36) and the change in coordinates defined in previous Section 4 (with all reference signals set to zero), we can rewrite subsystems (31a) and (37) for the stabilization of the nonholonomic system as Theorem 3. Consider the system (41a) and (41b), where the uncertainty term satisfies the inequality (22).Then, the robust backstepping SMC u 1 (32) and ( 33) and u 2 (38) and (39) (with all reference signals set to zero), yields asymptotic stabilization to the origin and reach the set {|s i | ≤ µ i } in finite time.

Output Feedback Control
Suppose the uncertain nonholonomic system (18a)-(18d).Consider an internal observer comprises of internal dynamics η and a virtual output ζ for the auxiliary system Remark 4. η = ψ(q) dynamics is a function of generalized coordinates q = [q 1 . . .q n ] T as described in (15).Assuming that output functions y i = h i (q) in Equation ( 11) are selected to ensure that all generalized coordinates q can be measured through the output functions.
Thus, the internal observer can be formulated as where L η is the observer gain.The high gain observer for the estimation of external dynamics of the nonholonomic system (18b) and (18c) is given by The observer gains T and the gain parameters h i,1 , h i,2 are selected such that the polynomial s 2 + h i,1 s + h i,2 is Hurwitz.Furthermore, the parameter ϵ i is a small positive constant.
Theorem 4. Consider the closed-loop system (46a) and (46b) under output feedback controller (47a) and (47b), designed using the full order HGO (48a) and (48b).Theorems 1-3 are satisfied for the asymptotic tracking and stabilization of the closed-loop system (19), respectively, and assume R is its region of attraction.Let Q be any compact subset of R ρ(n−m)+m , and S be any compact set in the interior of R.Then, µ * i > 0 and ϵ * i > 0 exist such that for every 0 < µ i < µ * i and 0 < ϵ i < ϵ * i , the solutions of the closed-loop systems (46a)-(46b) and (48a)-(48b), starting in S × Q, are bounded and the estimation error approaches zero as t → ∞.

Proof. Consider the Lyapunov function for the estimation error of internal dynamics (51), such as
V η (t, η) = ηT P −1 η satisfying the inequality 0 < c 1 I m ≤ P −1 ≤ c 2 I m , where c 1 , c 2 are positive constants, and P −1 is a positive definite matrix.It also verifies that where c 3 , c 4 > 0 and are independent of ϵ.The compact set Q 1 ⊂ R m will ensure that inequality (53) and (54) are valid.Now, applying the change in coordinate Υ = t/ϵ in (52) to obtain the boundary layer model, then substituting ϵ = 0, yields The Lyapunov function candidate for the boundary layer system (55) can be defined as where W is the positive definite solution of Lyapunov function WΛ + Λ T W = −I, satisfies Additionally, it can be shown that time derivative of V η and V χ also satisfy where k 5 and k 6 are positive constants independent of ϵ.Therefore, for all 0 < ϵ i < ϵ * i , t ≥ t 0 , and trajectories ( η(0) and ξ(0)) starting in Q 1 × Q 2 , the dynamics of estimation errors ( η, χ) converge to origin as t → ∞.

Design Example: Differential Drive WMR 7.1. Problem Formulation
This paper presents the differential drive type (2,0) WMR as an example of performing robust output feedback stabilization and trajectory tracking, as illustrated in Figure 1.The posture of a mobile robot in an inertial Cartesian frame O, X, Y can be described by q = [x 0 , y 0 , θ] T .The parameter r denotes the radius of each driving wheel, separated by a distance 2L.The constraint vector J(q), S(q), and velocity vector ϑ(t) as defined in Equations ( 1)-( 3) are described as where ϑ 1 (t) and ϑ 2 (t) are the linear and angular velocities of WMR, respectively.Similarly, transformation matrix Q and angular wheel velocity vector ϖ from Equation ( 4) are defined as The kinematic model of WMR by Equation ( 5) where p 1 = r 2 and p 2 = r 2L .The following matrices define the dynamic model as introduced in Equation ( 6) where m = mc + 2 mw and . The parameters mw and mc are the mass of driving wheels, including rotors of the DC motor, and the mass of the mobile robot platform, respectively.The parameter I w represents the moment of inertia of each wheel with a motor rotor, and I c is the moment of inertia of the WMR platform.The output vector (11) for (n − m = 2) outputs is selected as and the suitable choice of internal dynamic (17) would be where the decoupling matrix Φ(η, ξ) (12) can be defined as The proposed nontriangular normal form for an uncertain wheeled mobile robot, according to the notations introduced in Equations ( 15)-( 19), we obtained where P 3 = L rI , P 4 = 1 r m , and P 5 = md I , are calculated from the physical parameters of WMR.The mobile robot parameters such as mass ( m), wheel radius (r), distance between two wheels (2L), and moment of inertia (I) are assumed to be uncertain.Therefore, P 1 , P 2 , P 3 , P 4 , and P 5 are also uncertain.

Trajectory Tracking of a WMR
Applying the nonlinear feedback control (20) to obtain input-output feedback linearization and decoupling of the WMR described in Equations (60a)-(60c), and then using the results of Equations ( 21)-( 28), the error dynamical model (29) for trajectory tracking can be described as follows: The virtual control inputs ( 30), (34), and ( 36) are chosen as Then, apply the change in variables to transform the error dynamics of a mobile robot from Equations (61a) and (61b) into a similar form as described in Equations (31a) and (37), respectively, The internal dynamics η e in (64b) will remain bounded and exponentially converges to zero if the following requirement is fulfilled to evade singularity: |ξ 1,2r (t) + φ 1 (ξ 1,1e )| = 0, as stated in Equation (62b), To meet the above requirement labeled C1, it is essential to initialize the reference trajectory ξ 1,2r (0) and the error trajectory ξ 1,1e (0) appropriately.Meeting the condition of C1 means ensuring that |ξ 1,2r (0)| > |φ 1 (ξ 1,1e (0))|, will protect against singularity, when error dynamics ξ 1,1e approaches zero in the initial transient phase.
We can apply the backstepping sliding mode control law u 1 (32) and ( 33) and u 2 (38) and (39) to the error dynamical model Equations (64a) and (64b) to perform trajectory tracking of a WMR.The values of u 1eq and u 2eq in (33) and (39), respectively, can be formulated as

Posture Stabilization of a WMR
Consider the system (61a) and (61b) with all reference signals set to zero, we obtained Then, employing the same virtual state feedback control developed in Equations (62a)-(62c), with all reference signals set to zero Later, using the results of Equations (63a) and (64b) (setting all reference signals to zero) to transform the dynamics of the system (67a) and (67b) into equivalent representation Proposition 3. The internal dynamics η in (69b) will remain bounded and exponentially converge to zero if the given requirements are fulfilled to achieve the bounded solution of Equation (68b) Proof.In this Proposition, we will prove that Equation (68b) will remain bounded as the solution of ξ 1,1 (t) converges to the origin, ∀ t ⩾ 0. For this purpose, let us derive the solution of subsystem (69a) by substituting the control law u 1 as defined in the previous section (setting all reference signals to zero) Integrating the closed-loop system (70), we obtain where where the eigenvalue is calculated as Considering that condition C2. in Proposition 3 emphasizes that the gain parameters of the proposed control law u 1 are taken to ensure that the eigenvalues of subsystem (70) lie in the open left-half plane.Therefore, both ξ 1,1 and ξ 1,2 states of WMR will never cross the zero.In addition, C3. restricts the initial condition of ξ 1,1 (0) ̸ = 0 to avoid singularity during the initial transient in (68b).Hence the solution of Equation (68b) will remain bounded ∀ t ⩾ 0.
Posture stabilization of a nonholonomic mobile robot (69a) and (69b) can be attained by employing the backstepping SMC u 1 from Equations ( 32) and (65) and u 2 from Equations ( 38) and (39) (by setting all reference signals to zero).

Output Feedback Control
The closed-loop system (60a)-(60c) is further designed under output feedback control using the notation described in Equations (47a) and (47b) and Equations (48a) and (48b).The proposed full-order high gain observer for differential drive WMR are

Simulation Results
This section presents numerical simulations of robust trajectory tracking and posture stabilization of differential drive WMR using state feedback and output feedback backstepping sliding mode control (BSMC) under bounded parameter uncertainties and external disturbance.The nominal values of the physical parameters of the mobile robot are as follows: r = 0.05 m, m = 4 kg, 2L = 0.27 m, d = 0.05 m, and I = 2.5 kg.m 2 .During the simulation, parameter uncertainty on all physical parameters are assumed up to 20% of their nominal values, and an external disturbance of τ d = [5 sin(3t), 5 sin(3t)] T is applied.The control parameters for trajectory tracking are chosen as follows: 0) = 57 • , and θr = 5 (rad/s).
The control parameters a and γ in backstepping SMC u 1 (32) and u 2 (38) can be calculated using the inequalities defined by Equations (23a) and (23b).By choosing the small value of µ in backstepping SMC, we can recover the performance of discontinuous SMC.However, if µ is excessively small, it may lead to chattering.Hence, a value of µ = 0.1 is selected to mitigate the occurrence of chattering effects.The controller gain parameters k 1 , k 3 , and k 4 are chosen to ensure asymptotic stabilization of WMR towards the origin.The observer gain parameters h i,1 and h i,2 are selected such that the polynomial s 2 + h i,1 s + h i,2 is Hurwitz.Furthermore, it is desirable to have a small value of ϵ i in the high gain observer, to ensure that estimation errors quickly converge to zero.
The first simulation displays the robust trajectory tracking of a WMR using the state feedback control, as illustrated in Figure 2. The performance of the proposed backstepping SMC is compared with the approach presented in [1] under bounded parameter uncertainties and external disturbances.In this simulation, a lemniscate curve reference trajectory is selected as follows: x r (t) = sin(0.04t),y r (t) = sin(0.08t) In particular, Figure 2a,b display the trajectory tracking of a mobile robot to a lemniscate curve reference trajectory.The convergence of tracking errors (x e , y e , θ e ) to zero is presented in Figure 2c.
It can be observed that the proposed robust backstepping sliding mode control minimizes tracking errors compared to previous backstepping control [1], which cannot efficiently deal with these uncertainties.(c) (d) Figure 3a-d display the performance of the tracking controller in the presence of white noise.The presented results demonstrate that the suggested backstepping SMC control accurately tracked the reference trajectory in the x − y plane compared to the previous control approach described in [1].Furthermore, Figure 3c shows that the tracking errors (x e , y e , and θ e ) asymptotically converge to zero even in the presence of white noise using the proposed controller.The second simulation (Figure 4) compares the performance of trajectory tracking under state feedback (SFB) and output feedback (OFB) controllers.For simulation purposes, a circular trajectory is defined as: x r (t) = cos(0.05t),y r (t) = sin(0.05t) Figure 4a shows that trajectory tracking under output feedback with estimated states approaches the response obtained under the state feedback control, despite parameter uncertainties and time-varying external disturbance.Asymptotic convergence of the tracking error, i.e., x e = x − x r to zero is presented for both scenarios in Figure 4b.The performance of linear and angular velocities using both control techniques is depicted in Figure 4c.The estimation errors e x = x − x, e y = y − ŷ, and e θ = θ − θ, under the output feedback design are displayed in Figure 4d, showing that the estimation errors converge to zero quickly.Furthermore, Table 2, compares the performance of the proposed controller in comparison with the previous controller presented in [25].Simulation results illustrated in Figure 4 are compared with simulation results (Figures 2-7) presented in [25].The control parameters selected for robust posture stabilization are as follows: 0) = ŷc (0) = 0, and θ(0) = 45 • .Figure 5 displays the robust posture stabilization of a mobile robot using state feedback control.The performance of posture stabilization using the proposed backstepping sliding mode control (BSMC) is compared with the backstepping controller [1] under bounded parameter uncertainties and external disturbance.As demonstrated in these graphs, the proposed control law successfully achieves posture stabilization compared to the previous approach [1].Finally, Figure 6a displays the posture stabilization performance analysis under output feedback and state feedback control schemes.Figure 6b,c show that the output feedback controller based on HGO can recover the state feedback performance by converging the estimated state to the true state rapidly even in the presence of parameter uncertainties and time-varying external disturbance.Figure 6d displays the estimation errors of posture variables, i.e., e x = x − x, as well as the linear and angular velocities under output feedback control.It illustrates that estimation errors quickly converge to zero due to a small value to ϵ = 0.01 in the high gain observer.

Conclusions
This paper focuses on the design of a robust output feedback control that is simultaneously valid for both trajectory tracking and stabilization for a class of MIMO underactuated nonholonomic systems.The control design is based on kinematic and dynamic models in a globally defined normal form.We also consider that the nonholonomic system is under the influence of bounded model uncertainties and disturbances and in the absence of velocity measurements.The control and stabilization of these systems are widely regarded as the most difficult benchmark problems because of various inherent difficulties.These include: (i) it is not exactly feedback linearizable (input-state linearizable) because of nonholonomic constraints, (ii) obtaining a cascade normal form structure in the strict feedback form is not always possible, (iii) internal dynamics of nontriangular normal form are non-affine in control, i.e., highly nonlinear, and (iv) zero dynamics of the system are not minimum phase.The significant contribution of this research is the formulation of a generalized normal form using an input-output feedback linearization and change in coordinates approach that provides ease in designing (1) unified backstepping sliding mode control solutions for both stabilization and trajectory tracking subject to external disturbances and parameter uncertainties and (2) full-order HGO for the estimation of the derivative of output functions and internal dynamics.Moreover, a full-order high-gain observer and the backstepping sliding mode control are integrated to synthesize a robust output feedback controller.The proposed backstepping SMC reduces or eliminates the chattering phenomenon by dividing the control law into continuous and switching components, i.e., the high-slope saturation function.A differential drive type (2,0) WMR is used as an example to show the performance of the proposed controller.The simulation results illustrate that the suggested control law is robust against the bounded uncertainties, and the output feedback controller with estimated states can retrieve state feedback control performance (trajectory tracking and posture stabilization) even in the presence of such uncertainties.Therefore, both asymptotic trajectory tracking and posture stabilization are achieved in semi-global regions, with nonzero initial condition of the heading angle compared with nonzero initial condition of desired velocities in previous research.

Future Work
One promising area of potential future research in the field of nonholonomic systems is the design of robust control methodologies specifically tailored for systems affected by unmatched uncertainties.Another possible direction for future work is the study of formation/consensus control of wheeled mobile robots, as cooperative control of these vehicles and its real-world implementation continue to be active areas of investigation.
Angular Velocity via the BSMC Angular Velocity via the Ref.[1]

Figure 2 .
Figure 2. Lemniscate curve Trajectory tracking performance of WMR for BSMC (dashed line) and backstepping controller in reference [1] (solid line) under parameter uncertainties and external disturbance (a-d).
Linear Velocity via the BSMC Linear Velocity via the Ref.[1]

Figure 3 .
Figure 3. Circular trajectory tracking performance of WMR for BSMC (dashed line) and backstepping controller in reference [1] (solid line) under parameter uncertainties and external disturbance (a-d).

Figure 4 .
Figure 4. Trajectory tracking performance analysis of WMR using SFB and OFB controllers under parameter uncertainties and external disturbance (a-d).

Figure 5 .
Figure 5. Posture stabilization performance of WMR for BSMC (dashed line) and backstepping controller in reference [1] (solid line) under parameter uncertainties and external disturbance (a,b).

Figure 6 .
Figure 6.Posture stabilization performance analysis of WMR using SFB and OFB controllers under parameter uncertainties and external disturbance (a-d).

Table 1 .
Table of nomenclature.

Table 2 .
Output feedback trajectory tracking performance analysis.≈5 Estimation error e v 1 = v 1 − v1 zero not equal to zero Estimation error e v 2 = v 2 − v2 zero not equal to zero