Robust Tracking Control of Wheeled Mobile Robot Based on Differential Flatness and Sliding Active Disturbance Rejection Control: Simulations and Experiments

This paper proposes a robust tracking control method for wheeled mobile robot (WMR) against uncertainties, including wind disturbances and slipping. Through the application of the differential flatness methodology, the under-actuated WMR model is transformed into a linear canonical form, simplifying the design of a stabilizing feedback controller. To handle uncertainties from wheel slip and wind disturbances, the proposed feedback controller uses sliding mode control (SMC). However, increased uncertainties lead to chattering in the SMC approach due to higher control inputs. To mitigate this, a boundary layer around the switching surface is introduced, implementing a continuous control law to reduce chattering. Although increasing the boundary layer thickness reduces chattering, it may compromise the robustness achieved by SMC. To address this challenge, an active disturbance rejection control (ADRC) is integrated with boundary layer sliding mode control. ADRC estimates lumped uncertainties via an extended state observer and eliminates them within the feedback loop. This combined feedback control method aims to achieve practical control and robust tracking performance. Stability properties of the closed-loop system are established using the Lyapunov theory. Finally, simulations and experimental results are conducted to compare and evaluate the efficiency of the proposed robust tracking controller against other existing control methods.


Introduction
The domain of robotics finds mobile robots to be particularly intriguing, attracting considerable fascination and study.Designed to operate in dynamic settings, be it indoors or outdoors, these robots demonstrate the capacity to navigate autonomously or with minimal human input.Central to their functionality is their mobility, achieved through diverse locomotion methods, such as wheels, tracks, or legs.This mobility empowers them to traverse diverse terrains, overcoming obstacles encountered during their journeys.Recently, mobile robots have been used in various domains, including civilian, industrial, and military, to carry out diverse tasks such as surveillance [1], transportation [2], agricultural operations [3], and exploration [4].Given the broad application spectrum and critical nature of tasks involving mobile robots, there exists an urgent need to develop performance tracking controllers to execute proposed missions with exceptional accuracy.However, achieving this objective remains a significant challenge due to the inherent under-actuation and nonlinearity in WMRs, constrained by nonholonomic limitations.Consequently, researchers have directed their efforts towards investigating the control of mobile robotic systems.
In the past few decades, substantial progress has been made in the field of tracking control for wheeled mobile robots (WMR) through the application of nonlinear control theory [5][6][7].Among these control methodologies, linearization controllers, such as the flatness controller [8], have risen as a popular approach that can significantly simplify the controller design process.The flatness property is a technique used to define the dynamic behavior of nonlinear underactuated models by identifying a set of core system variables known as flat outputs.This perspective has significant implications for control systems, as will be demonstrated.The first step in flatness control involves generating a desired realizable trajectory that implicitly incorporates the system model.Following that, the nonlinear WMR model can be linearized, resulting in the canonical Brunovsky form [9,10].This special form simplifies the concept of a feedback controller capable of achieving exact trajectory tracking.In fact, controlling a linear system is easier than controlling an underactuated nonlinear system, and this feature has encouraged researchers to use the properties of flatness in several application domains, such as the control of hydraulic systems [11], exoskeleton robots [12], microgrid [13], underwater robot [14], and quadrotor [15,16].
Numerous research studies on WMR have utilized the concept of flatness control.Abadi [17] introduced an approach for optimal path planning for WMR using the collocation method, flatness control, and spline curves.This method effectively reduces the time needed to compute optimal robot trajectories during navigation, which is crucial for real-world applications.Kaniche [18] proposed a flatness visual servoing control for WMR subjected to disturbances.Salah [19] developed an approach to generate the upper coverage trajectory of a mobile robot by leveraging flatness.Yakovlev [20] combined flatness control with predictive control to enable safe navigation of a WMR among static and dynamic obstacles.
There is always is a difference between the mathematical model describing the movement of WMR and reality.This difference is due to environmental phenomena neglected during modeling, such as wind, slipping, etc.The question that arises is how flatness control applied to WMR can ensure the accurate tracking of a desired trajectory despite the presence of uncertainties.To resolve this problem, a robust feedback controller must be combined with flatness, taking into account the impact of uncertainties to the model.Up to the present, there have been limited methods in the literature concerning the robustness issues of flatness systems.Among these approaches, the sliding mode control (SMC) has been successfully utilized in a variety of systems [21][22][23].
SMC is a robust control technique used to manage dynamic systems in the presence of uncertainties and disturbances.At its core, SMC aims to drive the system state onto a designated sliding surface within the state space.Once on this surface, the system's behavior is constrained, allowing for effective regulation.SMC achieves this through discontinuous control actions, known as switching control, which dynamically alternate between different control laws.This switching mechanism ensures that the system remains on the sliding surface, enhancing robustness against external influences.Despite its effectiveness, SMC is associated with a phenomenon called chattering [24], characterized by rapid switching between control actions near the sliding surface.While chattering can theoretically improve tracking accuracy, it can lead to practical issues such as mechanical wear and high-frequency oscillations.To resolve this problem, numerous approaches have been suggested in the existing literature, such as high-order SMC [25], boundary layer [26], and active adaptive continuous nonsingular terminal sliding mode algorithm [27].A frequently utilized approach for mitigating the chattering phenomenon involves incorporating the boundary layer technique within SMC.This entails replacing the sign function with a smooth function.However, this strategy presents its own set of challenges.Firstly, there exists a trade-off between the size of the boundary layer and the performance of SMC, which impacts the effectiveness of chattering reduction.Secondly, the robustness and accuracy of the system may not always be guaranteed within the boundary layer.Additionally, beyond addressing the chattering issue, achieving precise control of a robotic system necessitates knowledge about its state, typically obtained through real instruments, which can incur high costs and complicate the system's structure.Moreover, in many instances, directly measuring certain system parameters may be impractical.To overcome these limitations, one potential solution involves implementing software sensors or observers, commonly referred to as virtual sensors.Therefore, to tackle both the reduced robustness resulting from the boundary layer approach and the challenge of state estimation, we propose a novel robust feedback controller that integrates the boundary layer sliding method with a disturbance observer.
In recent times, disturbance observers have emerged as potent tools for handling consolidated uncertainties, closely tied to disturbance-observer-based control.Within the domain of nonlinear disturbance observers, two notable approaches stand out: the uncertainty and disturbance estimator (UDE) [28] and the active disturbance rejection control [29] based on the extended state observer (ESO).In the UDE, only the disturbance is estimated, though in general, the observer equations depend on system states and inputs.Thus, a state observer is necessary unless all states are measurable.The idea of the ESO is to extend the original state vector by the disturbance vector and possibly, some of its time-derivatives, and then design a state observer for the extended system.ESO distinguishes itself by incorporating a dynamic model of disturbances or uncertainties into its estimation methodology, enabling it to identify and mitigate uncertainties not explicitly accounted for in the system model.The design of ESO is distinguished by its minimal dependence on system data and its freedom from the traditional system model, which simplifies its implementation process.Furthermore, several other types of disturbance observers are available, such as the nonlinear extended state observer (NLESO) [30], the adaptive extended state observer (AESO) [31], and the extended high-gain observer (EHGO) [32].EHGO, part of the ESO family, stands out in two key aspects: it does not necessitate slow variations in disturbances, and it estimates a matched disturbance term originating from model uncertainty and external disturbances.Given the advantages offered by ESO, considerable research effort has been devoted to developing advanced controls for robotic systems.
In Ref. [33], Xie introduced a controller that integrates the backstepping technique with ESO to improve tracking performance for underwater robots.Additionally, Qi [34] improved the bandwidth of ESO to achieve more accurate disturbance estimation.Subsequently, they utilized a simple feedback controller to ensure attitude stabilization over a 3D hovering quadrotor system.In the work by Aole [35], an improved ADRC methodology for controlling lower limb exoskeletons is presented.The proposed approach integrates Linear ESO with a tracking differentiator, nonlinear state error feedback, and a proportional controller.Simulation results demonstrate the effectiveness of the suggested ADRC in efficiently regulating the hip and knee movements of the robot in the presence of disturbances.Hu [36] integrated a predictive control technique with ESO for unmanned underwater vehicles, offering a solution to concurrently handle external disturbances and system measurement noises.Based on this observation, the main contributions of our research can be summarized as follows: 1.
The kinematic model for WMR is structured in a standard format that systematically tackles underactuation and transforms nonmatching disturbances into matching ones through a flatness-based approach; 2.
The designated trajectory is feasible in practice because of the concept of differential flatness, which equates differential flatness with controllability, ensuring its physical achievability; 3.
Continuous sliding mode control (SMC) is employed to eliminate chattering, an essential necessity for the efficient application of control in real-world scenarios; 4.
SMC is integrated with ESO for the uncertain kinematic WMR model.This strategy seeks to improve the practicality and resilience of the tracking controller by reducing chattering through boundary layer SMC and estimating the lumped disturbance affecting the WMRs via ESO, which is then employed as a feedforward compensation; 5.
The proposed control method was compared with several other control methods, including traditional flatness control, backstepping tracking control flatness-based sliding control, and flatness active disturbance rejection control and backstepping sliding active disturbance rejection control.These comparisons were validated through simulations conducted in Matlab/Simulink and experiments carried out on the Turtle-Bot WMR.
The structure of the remaining sections of this article is outlined as follows.Section 2 provides a thorough overview of the flatness control technique for WMR.Section 3 elaborates on the concept of flatness-based sliding mode tracking control of WMRs.The proposed robust tracking controller is delineated in Section 4. Sections 5 and 6 present and discuss the results of simulations and experiments.Finally, Section 7 concludes the paper by summarizing the key findings and suggesting potential future directions.

Flatness-Based Tracking Control
In our study, we analyzed a differential two-wheeled mobile robot (see Figure 1) that consists of two independent active wheels and a third passive wheel (a standard freewheel).This robotic system is widely regarded as an effective trade-off between control ease and the degrees of freedom that enable the robot to meet mobility requirements.The configuration of the mobile robot with wheels can be described by the vector q r = [x, y, θ].In this notation, x and y represent the coordinates of the robot's center position in the stationary frame (O, X, Y), while θ represents the orientation angle of the robot.The state equation of the WMR kinematic model, neglecting uncertainties, is represented as follows: The robot's translational and rotational velocities are denoted by v and w, respectively.The angular velocities of the right and left wheels (w r and w l ) can be defined as functions of the robot's translational and rotational velocities as follows: where the variables r and 2b represent the radius and distance between the wheels, respectively.The nonholonomic limitation is defined as follows, based on the nonslip requirement: The accuracy of the tracking will be guaranteed through the flatness property, which involves describing all system states and inputs, as well as their finite time derivatives, within the framework of a flat output.Considering the following nonlinear system: where x ∈ R n and u ∈ R m represent the state and the input vector.
The nonlinear system ( 5) is differentially flat if there exists an output λ in the following form: such that the state and the input can be expressed as follows: where a and c are finite multi-indices, and ξ, κ 1 , and κ 2 are smooth vector functions of the output vector λ and its derivatives.By introducing the functions κ 1 and κ 2 , this flat output is composed of a set of variables that enable the parameterization of all other system variables: the state, the command, and also the output λ.Indeed, if the output of the system is defined by a relation of the form λ = Ξ(x, u, u, . . ., u (p) ), then necessarily, the quantities described in Equations ( 7) and ( 8) make it possible to affirm that there exists an integer c such that: The flat output combines all unconstrained variables of the system since the components of λ are differentially independent.Alternatively, based on Equation ( 9), we can argue that the flat output λ solely relies on the state and the command.This would make it an endogenous variable of the system, in contrast to the state of an observer, which would be an example of an exogenous variable of the observed system.In addition, Lie-Bäcklund's notion of differential equivalence [8] shows that the number of components of λ is the same as the number of components of the control: This fundamental characteristic allows us to determine the requisite number of independent variables needed in a model to establish its flatness.A key benefit of the flatness property lies in its facilitation of various transformations, such as diffeomorphism and feedback linearization.These transformations enable the conversion of a nonlinear system into a controllable linear system, where the flat outputs represent the state vector.
Several studies in the literature, including Ref. [37], have shown that the WMR kinematic modeling can be defined as a differentially flat model, where the positional coordinates denoted as λ = [λ 11 , λ 21 ] T = [x, y] T serve as the flat outputs.Therefore, the entire set of state and control components pertaining to the WMR system are expressed using the flat variable λ and its derivatives, as demonstrated below: The differentially flat nature of the WMR's kinematic model has been demonstrated in the literature by various researchers [37].This implies that all the states and controls of the kinematic WMR model can be expressed as functions of λ and its derivatives.However, the noninvertible relationship between the control input vectors w and v and the highest derivatives of the flat output limits the development of static feedback linearization for the nonlinear WMR.To address this constraint, we incorporate the control input v into the kinematic model defined by Equation ( 1) by treating it as an additional state.As a result, we obtain a revised system that can be defined as follows: The state and control inputs of the modified system defined by Equation ( 14) are represented by X r = [x, y, v, θ] T and u r1 = v and u r2 = w.In order to establish a bijective relationship between the inputs u r1 , u r2 , and higher-order derivatives of λ 11 = x, λ 21 = y, we apply successive differentiations to the flat outputs until at least one of the input variables appears in the resulting expressions, as illustrated below: where B rob is described as follows: The matrix B rob is not singular if v ̸ = 0.In this case, we can define the control as follows: To arrive at the linearized system, referred to as the Burnovsky Form (BF), we can substitute the control input (17) into Equation (15).This substitution yields the following modified expression: where v 1 and v 2 represent a suitable feedback controller defined as follows: where λ xd and λ yd denote the desired trajectories for the flat output λ 11 and λ 21 , respectively.Meanwhile, the controller gains are represented by σ x1 , σ x2 , σ y1 , and σ y2 .The polynomial of the Burnovsky system (18) can be defined as follows: where the parameters m x and m y are the damping coefficients, and ϵ xc and ϵ yc are the frequencies in Equations ( 21) and (22).We can calculate the controller gain as follows: By integrating the feedback law, as described in Equations ( 19) and ( 20), into the system (17), we can express the flatness-based tracking control (FBTC) utilized for the mobile robot in the following manner: where e 1 = λ 11 − λ xd and e 2 = λ 21 − λ yd .
In ideal conditions where uncertainties such as wind and wheel slip are negligible in the kinematic model of the WMR, the control input defined by Equation ( 24) can achieve satisfactory tracking performance for the desired trajectory.However, it is practically impossible to have a model that accurately represents the real-world movement of the robot in all environmental conditions.As a result, the following section will focus on developing a robust tracking control for a WMR kinematic model that is subject to uncertainties.

Flatness-Based Sliding Tracking Control
In order to account for real-world conditions, we consider uncertainties such as slippage and external environmental disturbances when describing the kinematic model of WMR (Figure 2).As a result, the model is defined differently, as shown below.(Uncertain Kinematic Model) where the variables p x and p y represent the external environmental disturbances, indicating the potential influences from the surrounding conditions.On the other hand, v t and v s represent the slip velocities, where v t denotes the slip velocity along the forward direction and v s represents the slip velocity normal to it.Additionally, w s denotes the angular slip velocity.According to [37], it is assumed that the slippage phenomenon can be defined and bounded as follows: where κ 1 , ε 1 , ε 2 and ε 3 are positive constants.
Assuming that λ xd and λ yd are the reference trajectories for λ 11 and λ 21 , respectively, we can define the error dynamics as e i = λ i1 − λ id for i = 1, 2. To achieve convergence of the tracking error e i to zero in the presence of uncertainties, we employ a sliding mode control approach that relies on the principles of the flatness law.By incorporating this control strategy, we aim to ensure robust and accurate tracking performance even in the face of system uncertainties.The design of the sliding mode control involves two essential stages: the choice of the sliding surface and the development of the control law.These steps play a crucial role in establishing an effective and stable sliding mode control strategy.The selection of the sliding surface determines the desired system behavior and convergence properties, while the design of the control law focuses on generating control signals that guide the system towards the desired sliding surface and ensure its maintenance on that surface.In the context of the tracking example for the WMR, we make use of the sliding variable σ r = [s x , s y ] T to represent the tracking error.To define the sliding surface, we consider the desired tracking behavior and express it as follows, taking into account the specific requirements of the system: where the gains β 1 and β 2 can be selected using pole-placement techniques to ensure the asymptotic convergence of the tracking errors e 1 = λ 11 − λ xd and e 2 = λ 21 − λ yd to zero.In this tracking example, the sliding variable σ r = [s x , s y ] T is chosen as the tracking error.Therefore, the sliding surface for the WMR can be defined as follows: As suggested by Mauledoux [38], to guarantee that the sliding surface σ r = 0 is attractive, we can enforce the dynamics of σ r as follows: where the standard signum function is denoted by sgn, and k i (i = 1, 2) is a constant.One approach to proving the error dynamics stability is to analyze the following Lyapunov function: The derivative of V s is defined as follows: We can conclude that V s is a positive function and its derivative Vs is negative or zero.Hence, the system exhibits asymptotic Lyapunov stability.Using Equations ( 28), (29) and (32) we obtain: As a result, by using Equations ( 35) and (36), we can obtain: Substituting λ11 and λ21 with their new expressions defined by Equations ( 35) and (36) in the control defined by ( 17), the flatness-based sliding mode tracking controller (FSMC) applied to WMR is defined as follows: The FSMC defined by Equation ( 39) contains a discontinuous control term due to the function sgn(σ).Although selecting sufficiently large values for k 1 and k 2 can achieve convergence to sliding variable in limited time and provide robustness against perturbations, it also causes the phenomenon of chattering.Thus, to avoid this problem, the function sgn(σ) can be replaced by the function Sat defined as follows: where a s is the width of the threshold of the saturation function.
The thickness of the boundary layer, denoted as a s , within the saturation function stands as a pivotal parameter influencing the efficacy of the sliding mode controller.As the value of a s increases, the approximation diverges more from the ideal sgn function, resulting in enhanced reduction of chattering.However, this improvement comes at the cost of diminished robustness.Conversely, if the value of the parameter a s is reduced, the change of the control signal will be too frequent, which leads to inevitable chatter of the control signal.Therefore, a variable-thickness boundary layer a s is tailored to strike a balance between mitigating chattering and upholding system robustness amid uncertainties.In the upcoming section, the FSMC described by Equation ( 39) will be integrated with active disturbance rejection control to enhance the robustness lost by the Sat function and maintain the advantage of reducing chattering.

Proposed Robust Tracking Controller
In this section, we introduce a novel cascade control strategy that utilizes a combination of flatness property, active disturbance rejection control (ADRC), and boundary layer sliding mode control to solve the problem of reduced robustness obtained when replacing the function sgn by the function sat in the FSMC defined by Equation (39).Given the uncertain kinematic model (25), we can obtain the following relationship by differentiating λ 11 and λ 21 until the input terms u 1 and u 2 become evident: where C rob and D rob are defined as follows: By utilizing the control input described in Equation (17) on system (41), we achieve: where (43) in terms of two linear integrator systems subject to perturbation yields the following expressions: Consider ∆ 1 and ∆ 2 as the differentials of δ 1 and δ 2 with respect to time t, respectively.We assume that both δ i and ∆ i (i = 1, 2) are bounded.In practical applications, determining the actual values of the lumped disturbances δ 1 and δ 2 that affect the system is considered a challenging problem.Hence, an observer is required to estimate these values.

ESO Design
The extended state observer (ESO) plays a vital role in system control by simultaneously estimating the system states and uncertainties.This capability enables the ESO to effectively reject or compensate for disturbances, enhancing the system's robustness and performance.The ESO takes into account all factors that affect the system and treats parameter uncertainties and external perturbations as a single observed disturbance.The ESO is named as such because it estimates uncertainties as an extended state.Its benefits include not being reliant on the mathematical model of the system, as well as having a straightforward implementation and demonstrating good performance.Consider λ 13 = δ 1 , α 23 = δ 2 as an extended state for system (44).The latter can be expressed as follows: We can express systems (45) in matrix form as follows: where

The expression for the Extended State
Observer (ESO) corresponding to each extended system (46) and (47) can be given as follows: where α gx = [α 11 , α 12 , α 13 ] T , α gy = [α 21 , α 22 , α 23 ] T .To determine the observer gains α ij (i = 1, 2, 3), (j = 1, 2, 3), we can adopt the methodology proposed by Gao [39] outlined in the following manner: The choice of γ xo and γ yo is made to ensure that Equations (50) and (51) form Hurwitz polynomials with respect to the complex variable.The observer gain can be formulated as a function of the ESO bandwidth by utilizing Equations ( 50) and (51), as demonstrated below: The observer error associated with each ESO can be defined by employing Equations ( 46)-(49) as follows: It is possible to express Equations ( 53) and (54) in matrix form as shown below: where ê = [ê x , ˙ê x , ëx , êy , ˙ê y , ëy ] T , Ĥ Lemma 1.In Equation ( 55), the boundedness of lim t→∞ ê(t) can be guaranteed if at least one of the following two conditions is satisfied: Asymptotic stability of the estimated error dynamics can be achieved when the values of δ i , i = 1, 2, are either directly obtained or assumed to be constant, leading to ∆ i = 0, i = 1, 2. In this scenario, the positive constants n 1 and n 2 play a vital role in ensuring the system's stability.Lemma 1, as stated in Zhang et al. [40] , establishes that the roots of the matrix Ĥ in Equation (55) reside in the left half plane.This result is ensured by the nonnegativity of the bandwidths γ xo and γ yo .Consequently, it can be deduced that the estimated error dynamics described by Equations ( 53) and (54) are asymptotically stable.

New Robust Feedback Controller
The feedback controller presented in Equations (19) and (20) relies on state measurements, but except for λ 11 and λ 21 , the remaining states cannot be accurately measured.To solve this problem, the state estimation obtained through the two ESOs defined in Equations ( 48) and (49) are used instead.Furthermore, in order to simplify the compensation of the lumped disturbances δ 1 and δ 2 , they are replaced by their approximations, δ1 and δ2 .By incorporating the results of the extended state observers (ESOs), a robust feedback controller can be developed in the following manner: according to the sliding mode active disturbance rejection control feedback given in Equations ( 56) and (57), we can obtain the new robust tracking controller named Flatness-Sliding-Active-Disturbance-Rejection Control (FSADRC), defined as follows: where êr1 = λ11 − λ xd and êr2 = λ21 − λ yd .The schematic diagram presented in Figure 3 illustrates the principle of trajectory tracking control for a mobile robot.

Stability Analysis of the Closed-Loop System
This section will address the stability analysis of the tracking error systems for x and y, utilizing the estimation error defined by Equations ( 53) and (54).In order to prove the stability of the error dynamics of position x, Lyapunov's function is chosen as follows: where We can define the derivative of the Lyapunov function V sx as follows: When replacing λ11 by its Equation (44) defined by λ11 = v 1 + δ 1 , we obtain: When v 1 represents the feedback controller, substituting it with the proposed robust feedback tracking control, denoted as v SADRCx defined by Equation (56), yields: where s x is defined as follows: Concerning the stability and boundedness of the ESO defined by Equation ( 53), it can be achieved by choosing α gx in such a way that the eigenvalues of A x − α gx C x are negative, indicating poles in the left-hand plane, and ensuring that uncertainty is bounded.As a result, the error ˙êx → 0. This implies that λ11 → λ 11 , δ1 → δ 1 , and ŝx → s x .In this scenario, the Lyapunov function defined by Equation ( 62) is formulated as follows: Since Sat(s x ), defined by Equation ( 63), is divided into two segments, the proof process will be analyzed in two cases.In the first scenario, when the saturation function is defined as described by: Moreover, the Lyapunov function is defined as follows: Alternatively, when the saturation function is given by: the Lyapunov function takes the form: Thus, based on Equations ( 66) and ( 68), it can be concluded that the Lyapunov function Vsx is negative regardless of the definition of the function Sat(s x ).As a result, the tracking error of the position x is stable.Similarly, the same conclusions about the stability of the closed-loop system y can be drawn.

Simulation Results
This section presents simulation tests to validate the efficacy and superiority of the suggested controller, flatness sliding active disturbance rejection control (FSADRC), as defined by Equation (58).The proposed control is evaluated against flatness sliding mode control (FSMC), represented by Equation ( 39), and flatness-based tracking control (FBTC), as defined in Equation ( 24), using computer simulation results.The parameters of the WMR are r = 0.1 m, b = 0.15 m.To enhance the observation and comparison of the simulation results, we have chosen two types of reference trajectories: a circular path and a Bézier curve.Additionally, we also consider two different scenarios of perturbation.The controller design parameters of FBTC, FSMC, and FSADRC are chosen as m x = m y = 1, ϵ xc = ϵ yc = 2, β 1 = β 2 = 5, and k 1 = k 2 = 10.As suggested by Gao [39], it is advisable to select the observer bandwidth to be sufficiently higher than the controller bandwidth.This ensures that the observer dynamics remain faster than the system dynamics, enabling effective disturbance estimation and compensation.In our case, we have chosen observer bandwidths of γ xo = γ yo = 6 rad/s to fulfill this requirement and ensure robust performance of the control system.To ensure that the sliding mode control system achieves both satisfactory dynamic and steady-state performance, and to prevent chatter in the control signal, the cut and dry method is frequently employed to establish the thickness of the boundary layer.Specifically, in this case, a sx = a sy = 0.3 is chosen.

First Scenario
In this simulation, we consider that slip velocities v t and v s can be up to 30% of the forward speed.Thus, κ 1 = 0.3.In addition, the WMR is subjected to constant wind perturbation defined as follows: The reference trajectory considered in this scenario is a circle, which is defined by the following equation: x r = cos(t), y r = sin(t) (70) The performance of the uncertain WMR systems under different control strategies, namely FBTC, FSMC, and FSADRC, is depicted in Figure 4. Figure 5 shows the results of the estimated lumped disturbance affecting the x and y channels obtained using the extended state observer (ESO).Figure 6 illustrates the control input applied to the wheeled mobile robot under the conditions of the first scenario.The simulation results indicate that the uncertainty caused by slow wind perturbation and slip decreases the tracking performance in trajectory following, rendering FBTC ineffective as a controller.On the other hand, both FSMC and FSADRC demonstrate robustness in handling the overall disturbance affecting the WMR model.These controllers exhibit the ability to mitigate disturbances and successfully maintain the desired trajectory of the WMR system.Consequently, it can be inferred that controllers that disregard uncertain models, despite being feedback controllers, may exhibit unsatisfactory performance.The fundamental distinction between the FSMC and FSADRC controllers lies in their design methodologies and approaches.FSMC relies on finely-tuned gains to achieve disturbance rejection, which can lead to chattering due to the relatively high gain values.In contrast, FSADRC combines the advantages of the boundary layer method to minimize chattering and an ESO to estimate and eliminate lumped disturbance.

Second Scenario
The objective of this simulation is to create and follow a trajectory for a robot, starting from an initial state where x(0) = y(0) = 0, and reaching a final state specified by x(10) = 3.5 and y(10) = 5.This trajectory must navigate through a room containing obstacles, while also considering time-varying wind disturbances and slipping.The desired trajectory should meet the following criteria: minimizing energy consumption, maneuvering around static obstacles, and adhering to the specified state constraints as follows: The optimal trajectory generation method proposed in [17] offers a solution to obtain the desired trajectory by solving a nonlinear optimization problem.By integrating the principles of flatness, the collocation method, and B-spline functions, this method efficiently generates trajectories while guaranteeing constraint satisfaction.To ensure consistency in the simulation results, the parameters for all three controllers remain unchanged from the previous simulations.Considering an uncertain initial condition of x(0) = 1 and ŷ(0) = 1 for the wheeled mobile robot (WMR), we further specify that the slip velocities v t and v s can potentially reach up to 50% to 70%.In addition, we take into account the influence of sinusoidal wind disturbances.In contrast to the initial scenario, the disturbance signals consist of combinations of multi-frequency sinusoidal signals representing time-varying disturbances, particularly wind, defined as follows: The simulation results regarding trajectory tracking performance of the second scenario are depicted in Figure 7. Based on these figures, it can be observed that the WMR system experiences significant divergence from the desired trajectory when affected by slippage and external disturbances, rendering FBTC ineffective as a controller.The FSMC controller's intervention through the sliding mode's discontinuous term eliminates uncertainty effects and maintains the stability of the closed-loop control.However, as shown in Figure 8, the presence of chattering in the FSMC control signals negatively impacts the system's behavior.Hence, it can be inferred that while FSMC is a robust control approach, its practical applicability is quite restricted.Therefore, developing a control approach capable of mitigating the chattering effect while maintaining the robustness advantage provided by SMC is necessary.The results of the lumped disturbance estimation for this simulation are illustrated in Figure 9.According to the simulation findings, the mobile robot satisfactory trajectory tracking performance when confronted with model disturbances and uncertain initial conditions while employing the FSADRC controller.Of greater significance, the proposed control methodology achieves superior tracking of the desired trajectory, devoid of the chattering phenomenon, and enhances tracking performance against aggressive disturbances.

Tracking the Experimental Results of a Wheeled Mobile Robot
This section outlines experiments conducted with the TurtleBot3, a Wheeled Mobile Robot (WMR), to evaluate a proposed methodology.The TurtleBot provides a cost-effective platform for researchers to explore and validate control algorithms without requiring expensive robotic systems.Its compatibility with the Robot Operating System (ROS) enhances its functionalities, offering resources for algorithm development and experimentation.With Li-DAR, IMU, and wheel encoders onboard, the TurtleBot3 provides precise environmental feedback, facilitating algorithm optimization.Researchers can augment the system with additional sensors or hardware components to evaluate various control algorithms across diverse scenarios.To facilitate the observation and comparison of experimental results, we have selected two types of reference trajectories: an eight-shaped path and a Bézier curve.Additionally, we have considered two different scenarios of perturbation: the first involves slowly time-varying disturbances, while the second entails aggressive time-varying disturbances.For further validation, the performance of the proposed control method is compared with other state-of-the-art control techniques such as backstepping tracking control (BTC) [41], flatness active disturbance rejection control (FADRC) introduced in [42], flatness-based tracking control (FBTC) as defined by Equation (39), and backstepping sliding active disturbance rejection control (BSADRC) [43].The controller design parameters selected for the experimental results are identical to those chosen for the simulation results.

First Experiment with Slowly Time-Varying Disturbances
In this experiment, eight shapes were chosen for the reference trajectory, as outlined below: x r = 2cos(t), y r = −2sin(t) (73) To replicate real-world navigation conditions for the WMR, high-speed fans are utilized in the laboratory to simulate windy environments.Additionally, a stick is employed to disturb the castors of the WMR, creating slipping incidents, thus adding further realism to the testing environment.Figure 10 illustrates the real-time tracking of the eight-shaped reference trajectory of the WMR using the proposed control method described in this paper.Simulation of the experiment under identical conditions reveal tracking trajectories in Figure 11. Figure 12 illustrates lumped disturbance estimation, while Figure 13 displays control torques.Based on the experimental results shown in Figure 11, it is evident that FADRC, FSADRC, and BSADRC methods excel at tracking trajectories even in the face of genuine uncertainty.Conversely, the FBTC and BTC methods demonstrate significant shortcomings when it comes to handling uncertainties.To assess the superiority of the proposed control, we will conduct a thorough study in the subsequent section.This study will include a quantitative analysis of the controllers under more severe disturbance conditions.

Second Experiment with Aggressive Time-Varying Disturbances
In this experiment, we intensify the frequency of disturbance variation generated by the industrial ventilator and subject the robot to aggressive impacts with a stick to assess the effectiveness of the proposed controller.Additionally, we adopt the eighth-order Bézier curve as a reference trajectory for both the x and y positions, defined as follows: where P xj , P yj , and j = 0 . . .8 represent the control parameters of the reference trajectory.These parameters may vary depending on several factors, including the robot's initial position, the desired final position, and constraints such as obstacle avoidance.As an example, we select control parameters that allow the WMR to transition from its initial state q r (0) = [0, 0, 0] T to the desired final state q r (20) = [2, 2, 0] T .The tracking experiment results of the WMR under aggressive time-varying disturbances are depicted in Figure 14.In Figure 15, the lumped disturbance affecting the WMR within the context of the second experimental scenario is displayed.Similarly, Figure 16 illustrates the proposed control input applied to the wheeled mobile robot within the same context.To quantitatively assess the tracking performance of the WMR, we employed the integral absolute error (IAE) and the control effort performance index as comparison metrics.The IAE is computed for each of the control strategies in the following manner: where t f is the total simulation duration and i = 1, 2, represents the position in the x and y direction, respectively.The control effort is given as follows: where N indicates the total count of samples.The associated key performance indicators IAE and P avg for both strategies are provided in Table 1.Examining the data in Table 1, it is evident that the FSADRC controller outperforms the BTC, FBTC, FADRC, and FSMC methods in terms of tracking performance.Although its tracking performance is nearly comparable to that of the BSADRC, the FSADRC requires minimal effort to accomplish its task compared to the BSADRC.This characteristic is particularly crucial in contexts where energy resources are limited, such as in mobile or autonomous applications.The enhanced efficiency of the FSADRC over the BSADRC is explained by the advantage of flatness control, which simplifies controller design by transforming the nonlinear system into a linear one.This feature makes all control based on the concept of flatness less complex than control based on backstepping.Ultimately, the experiment and table findings show that the disturbance rejection function simplifies the system model by addressing real-time modeling uncertainties.Consequently, the FSADRC method relies less on an exact analytical model description, treating unknown dynamics as internal disturbances compensated for by the rejection function.This enhances the robustness of FSADRC, which also incorporates the boundary layer technique to alleviate chattering effects.

Conclusions
This paper aims to introduce a robust control methodology for uncertain wheeled mobile robots (WMR).By employing flatness-based control, the nonlinear kinematic model of the WMR undergoes transformation into a canonical form, enabling the implementation of a robust feedback controller that incorporates boundary layer sliding mode control and extended state observer techniques.Simulation results conducted under various scenarios of uncertainties illustrate the effectiveness of FSADRC in enhancing the trajectory tracking performance of the WMR when compared to BTC, FBTC, and FADRC, even amid variations in slipping and external wind disturbances.Furthermore, within the same context, FSADRC demonstrates comparable efficiency to BSADRC in terms of trajectory tracking, while exhibiting an advantage in effort usage due to its flatness property.The smooth operation of FSADRC, coupled with its resilience against parameter variations and external disturbances, renders it a practical choice for real-world applications.Moreover, experimental findings using the TurtleBot3 validate the efficacy of the proposed FSADRC in real-world navigational tasks.In future studies, the application of FSADRC will extend to other robotic systems, such as quadrotors and arm manipulators, to assess its effectiveness and explore its potential for broader deployment.
FSADRC Flatness sliding active disturbance rejection control BSADRC Backstepping sliding active disturbance rejection control

Figure 2 .
Figure 2. Two-wheeled mobile robot subject to uncertainties.

Figure 3 .
Figure 3. Mobile robot trajectory tracking control principle scheme.

Figure 4 .
Figure 4. Simulation tracking results of the wheeled mobile robot under the conditions of the first scenario.

Figure 5 .
Figure 5. Lumped disturbance affecting the x and y position channels in the context of the first scenario.

Figure 6 .
Figure 6.Control input applied to the wheeled mobile robot under the conditions of the first scenario.

Figure 7 .
Figure 7. Simulation tracking results of the wheeled mobile robot in the conditions of the second scenario.

Figure 8 .
Figure 8.Control input applied to the wheeled mobile robot under the conditions of the second scenario.

Figure 9 .
Figure 9. Lumped disturbance affecting the x and y channels in the context of the second scenario.

Figure 11 .Figure 12 .
Figure 11.Results of the wheeled mobile robot's tracking under the conditions of the first experiment scenario.

FSADRCFigure 13 .
Figure 13.Control torques applied to the right and left wheels to track the eight-shaped reference trajectory.

Figure 14 .
Figure 14.Results of the wheeled mobile robot's tracking under the conditions of the second experiment scenario.

Figure 15 .FSADRCFigure 16 .
Figure 15.Estimated values of the lumped disturbances under the conditions of the second experiment scenario.

Table 1 .
Performance indexes IAE and P avg .