Design of a Model Predictive Trajectory Tracking Controller for Mobile Robot Based on the Event-Triggering Mechanism

A novel model predictive control(MPC-) based trajectory tracking controller for mobile robot is proposed using the eventtriggering mechanism, and the aim is to solve the problem that the MPC optimization problem requires a large amount of online computation and communication resources. )is method includes two different event-triggering strategies, namely, the eventtriggering based on threshold curve and the event-triggering based on threshold band. )e selection of the triggering threshold is achieved by applying the statistical method to the historical data of the trajectory tracking of the mobile robot under the classic MPC method. Simulation and experimental tests illustrate that the proposed approach is able to significantly reduce the computation and communication burdens without affecting the control performance. Furthermore, the experimental results show that compared with the classic MPC-based tracking method, the proposed two event-triggering strategies can reduce 28.1% and 75.7% of the computation load and 0.886 s and 2.385 s communication time.


Introduction
With the rapid development of computer engineering, electronic engineering, network communication engineering, and other disciplines, mobile robot technology has made great progress and has been widely used in scientific exploration, national defense and military, risk relief, life services, and other fields. erefore, the trajectory tracking method of mobile robot using a conventional control strategy is difficult to achieve the desired control performance, so it is urgent to study the trajectory tracking strategy of a robot based on advanced control algorithms.
In recent years, many research achievements have been made on the trajectory tracking control of mobile robots. e mainstream tracking control algorithms at this stage included sliding mode control, model predictive control (MPC), and optimal control [1][2][3][4]. Due to the advantages in handling system constraints and balancing multiple control objectives, MPC has been applied in more and more practical mobile robots, such as wheeled robots [5][6][7][8][9] and underwater autonomous navigation robots [10,11]. In [8], the authors designed a MPC-based trajectory tracking control method, which could ensure that the unmanned vehicle track the reference trajectory quickly and stably; the distance error and heading error are in a reasonable range, and the real-time performance meet the requirements. [9] proposed a real-time optimization scheme to reduce the control horizon and the control update frequency for MPCbased robot path tracking control, which could balance the realtime requirements and control accuracy. An MPC-based path tracking controller for the constrained under-driven autonomous underwater vehicle (AUV) system is designed in [10], and it is pointed out that a long prediction horizon could be chosen to ensure the final convergence of the control system. In addition, [11] proposed a new MPC method for the trajectory tracking of underwater robots, which used the linearized error model to formulate the MPC optimization problem such that the speed jump problem could be avoided without violating the control constraints.
As MPC controllers solve a constrained optimization problem at each sampling instant, the online computation burden is quite large and needs to be further optimized. In [12], the authors proposed a MPC strategy based on the event-triggering mechanism, in which the optimization is only triggered when the difference between the predicted and the actual trajectories exceeded a certain threshold, which effectively reduced the online computation. [13] proposed an event-based MPC method and further proved the closed-loop stability property by constructing terminal constraint and terminal cost function. In [14], the authors proposed a novel aperiodic adaptive event-triggered communication mechanism for master-slave synchronization problem with aperiodic sampled data, which can effectively reduce the transmission load. For more details regarding communication links in tracking control, the reader is referred to [15] and references therein. [16] proposed an event-based MPC approach for linear discrete systems subject to external perturbation, in which the optimization problem included a time-vary tightened state constraint. It should be emphasized that compared with the standard MPC applied in practical robots, the above eventtriggered MPC algorithms including additional items such as terminal and tightened constraints, which cannot be directly adopted, and therefore, the event-triggering mechanism should be further studied for MPC-based tracking control of robot systems.
Aiming at solving the abovementioned problems, this study proposes a more practical event-triggered MPC-based trajectory tracking method for mobile robot subject to external disturbances. Two novel event-triggering strategies, i.e., event-triggering based on the threshold curve and eventtriggering based on the threshold band are developed. Furthermore, the selection of the threshold curve and threshold band is achieved via applying the statistical approach to the recorded historical trajectory data of the robot system with the classic MPC method. Finally, simulation and experimental tests show that the developed event-triggered MPC method can significantly reduce the computation and communication resources occupied by the MPC controller.
It is worth mentioning that the main motivation of this manuscript is that in order to facilitate the theoretical analysis, the majority of the existing event-triggered MPC methods have extra assumptions about the robot system and require additional items (e.g., terminal cost) in the MPC setting, which makes the corresponding MPC controller too computational expensive for many real robot systems. Since robot systems with relatively low hardware configuration are an important portion of industrial robots, it is of great importance to develop the event-triggered MPC method suitable for those robot systems. Furthermore, the main contribution of the proposed method can be illustrated in two aspects. On the one hand, compared with the event-triggered MPC with terminal cost/ constraint, the proposed method employed the standard robot kinematic model with no extra assumption and has no additional terms in the MPC cost function, which can be directly applied to real robot systems. On the other hand, compared with event-triggered MPC with standard MPC framework, the proposed technique incorporated the external disturbances into the development of the triggering condition and adopted statistical methods to construct two different kinds of triggering strategies for performance improvement (Section 5.2).

Description of the Robot System
2.1. Kinematic Model of the Robot System. First, the kinematic model of the mobile robot is constructed. It is assumed that the robot system conforms to nonholonomic constraints, and the robot body has no lateral sliding. Considering that the velocity of the robot is generally low and the lateral acceleration such as centrifugal acceleration has little influence during steering [17], the kinematic model for the mobile robot can be established as where the state is denoted as x y φ T , consisting of the position of the mobile robot x y , as well as its orientation φ. e input signal is represented as v ω T , including linear velocity and angular velocity of the mobile robot. e above kinematic model can then be expressed in a more general state space form as where ξ is the state, and u is the input; f(·) represents the corresponding mapping relationship. Since the trajectory status of the reference system is known, the relationship of the state vector and control vector of the reference system is expressed as Expand equation (2) by Taylor series at any reference point (ξ r , u r ) and retain only the first-order terms and ignore the higher-order terms. e following equation can be obtained: By making a difference between equations (3) and (4), the linearized error model is obtained as e linearized error model of the mobile robot is further discretized by Euler method as where is the sampling time.

Augmented State Space Model.
To ease the development of the model predictive trajectory tracking controller, an augmented state space model is constructed. Given the discrete time linear model of the mobile robot in (6), an augmented state is set as us, the desired discrete state space equation is obtained: where A k,t � A k B k ; 0 m×n I m is the system matrix, B k,t � B k ; I m is the control matrix, and C k,t � I n 0 is the output matrix, and m, n are the dimensions of the state and control vectors. To simplify the representation, the following expressions are employed: where N h is the horizon (including prediction horizon and control horizon).

Design of a Model Predictive Trajectory
Tracking Controller e objective function of the predictive controller contains information such as system state error and the change of control, based on which the optimal control problem can be formulated to ensure the mobile robot to track the reference trajectory quickly and smoothly. In this work, the cost function is formulated as where Q and R denote the state-weighting matrix and inputweighting matrix, respectively, ε is the relaxing factor, and α is the associated weighting coefficient. In this objective function, the first element shows the capability of the system to track the set point, and the second element illustrates its ability to constrain the change of the control input. e cost function (10) determines the closed-loop performance of the MPC system. More specifically, the dynamic and steady-state performances as well as the energy consumption of the robot system can be tuned by the design parameters (such as prediction horizon N p , control horizon N c , and weighting matrices Q and R) in (10). e prediction horizon N p mainly affects the stability and rapidity of the system. When N p is small, the rapidity of the system response is guaranteed, but the robustness of the system is relatively poor. erefore, the choice of N p mainly involves the trade-off between robust stability and rapidity of the system. e control horizon N c mainly affects the dynamic characteristics of the system. e smaller the N c is selected, the worse the tracking performance of the system will be. A larger N c makes the control signal more flexible and improves the overall dynamic characteristics of the system, but it is also prone to robust stability problems. e weighting matrices Q and R restrict each other. e rapid response performance of the system can be improved by reducing R (or increasing Q), and the steady-state performance and the energy consumption can be improved via increasing R (or reducing Q), since the change of the control signal can be effectively suppressed and the control cost can be reduced. e relaxing factor and its weighting coefficient can be tuned to guarantee that MPC admits an optimal input trajectory at each sampling instant, which indirectly guarantees the feasibility of the optimal problem. e prediction of the system output Y(t) during the prediction horizon is obtained: where Y(t) � By substituting equation (11) into the cost function in (10), the output deviation in the prediction horizon can be denoted as where rough some standard matrix operations, the cost function can be adjusted as Mathematical Problems in Engineering en, it can be shown that the constrained optimization problem of the predictive controller at each step is equivalent to the quadratic programming (QP) problem as where ΔU min and ΔU max are the system control incremental sequence constraints, U min and U max are the system control sequence constraints, and Y min and Y max are the system output sequence constraints.
By solving problem (14) at a given sampling instant, T can be obtained. e first item in the sequence is applied as the actual control signal increment, i.e., u(t) � u(t − 1) + Δu * t . en, the newly obtained state is utilized to update the optimization problem for the next instant. e aforementioned procedure is recursively performed, until the control process is completed. Besides, the convergence speed of the proposed MPC control algorithm can be controlled by the tunable parameters in the MPC cost function, namely, prediction horizon N p , control horizon N c , and weighting matrices Q and R. In general, when convergence speed of the system response needs to be improved, a larger N c and a smaller N p are usually selected, and R is reduced (or Q is increased). It should be noted that in the physical sense, N c <� N p needs to be satisfied, and the functions of the weighting matrices Q and R are mutually restricted.

Design of Event-Triggering Strategies
Compared with classic time-triggered MPC, event-triggered MPC has the characteristic of only performing actions at the time when a predetermined event occurs, such as the error exceeds a certain threshold, and thus, it can effectively reduce the frequency of system sampling and update, which could effectively reduce computation and communication load.
Due to the inevitable existence of various disturbances in actual mobile robot systems, such as model mismatch and parameter perturbation, the design of triggering conditions should further consider the external disturbances of the considered system. Denote the disturbance signal as w ∈ W with W being a compact set. en, on the basis of the nominal model in equation (1), a perturbed model of the robot system can be described as It is worth nothing that although the external disturbances can be handled via incorporating terminal cost function, terminal constraint, and tightening set to the MPC algorithm in (10); these additional items will significantly increase the amount of online computation and seriously affect the real-time property and performance of the controller. For the robot system with limited computation ability, such modified MPC optimization problem may not be solved in time to update the control signals, which will lead to system performance deterioration and even instability. To solve the abovementioned problems, this study proposes an event-triggered MPC-based trajectory tracking method for the robot systems, which can effectively deal with the influence of external disturbances on the system without increasing the amount of online computation.

Event-Triggering Strategy
where t k is the current triggering time, ξ 1 , ξ 2 , and ξ 3 are the values of position components of the systems x, y, and φ.
And σ x , σ y , and σ φ denote the event-triggering threshold corresponding to the state components, respectively.
Note that the threshold curve parameters σ x , σ y , and σ φ can be selected via taking the mean values of the position coordinates of the robot from N groups of historical data collected during the similar working condition under the classic MPC, and thus, the three thresholds at each sampling time can be obtained through statistical processing via where x i,k , y i,k , and φ i,k are the position and the orientation information of the mobile robot in i th group historical data, k is the k th sampling time; σ x (k), σ y (k), and σ φ (k) are the obtained thresholds corresponding to the position and orientation of the robot at sampling time k, respectively.

Event-Triggering Strategy Based on
reshold Band. To further reduce the frequency of solving the optimization problem and save computation resources, the second eventtriggering strategy, which is based on threshold band, is proposed in this section. More specifically, at some sampling time, when the position coordinate of any state component of the mobile robot exceeds the upper or lower bound of the threshold band, it is considered to be the case where the triggering condition is satisfied. To formulate the eventtriggering strategy, the triggering instant t k+1 is defined as where σ x_u , σ x_d , σ y_u , σ y_d , σ φ_u , and σ φ_d denote the upper and lower threshold bounds of the position and orientation. e selection method of the threshold band is that under the condition that the control effect does not decline; the threshold band can be formed by using the relationship between the threshold curve and the maximum disturbance. In actual working conditions, due to the existence of uncertain factors such as measurement noise, load changes, and external disturbances of the system, it is difficult to obtain accurate values of the parameters of the mobile robot model. Based on this, the model mismatch and parameter perturbation of the system are usually treated as a total disturbance.
According to a large number of historical data, the disturbance upper bound of the usual working scene can be measured, and a threshold band can be obtained by subtracting the disturbance upper bound from the threshold curve. Using the threshold band as the triggering condition can improve the robustness of the system, and the threshold band can be obtained as where σ x_u (k), σ y_u (k), σ φ_u (k), σ x_d (k), σ y_d (k), and σ φ_d (k) are the upper and lower threshold limits corresponding to the position and orientation of the robot at sampling time k. Remark 1. e method proposed in this study relies on the robot state information mainly in two aspects, namely, calculating the triggering threshold and determining the triggering instant; the former needs to process a large amount of historical complete state information, and the latter needs to know the current complete state information of the robot. However, in certain scenarios (such as linear trajectory tracking), the complete state information can be reduced to only state information containing the horizontal and vertical coordinates. In addition, we can also reduce the limit of requiring the complete state information by changing the triggering condition, e.g., developing a new condition with only the horizontal and vertical coordinates, but it may affect the control performance to some extent and needs further investigation.

Design of the Event-Triggered MPC Algorithm.
Based on the event-triggering mechanism proposed in Sections 4.1 and 4.2, the event-triggered MPC algorithm is provided. Note that N p � N c is utilized based on the practical experience. In the developed algorithm, the online optimization is executed only when the triggering condition is satisfied. To further lessen the times to solve the optimization problem, the last element in the control sequence is continuously used after the control horizon expires. us, the control law is constructed as e pseudocode of the proposed event-triggered MPC algorithm is given in Algorithm 1.

Remark 2.
Note that the computational complexity is minor while applying the proposed event-triggering mechanism. e reason is that only a simple-triggering condition is checked at each sampling instant, and the associated computation is quite simple in the machine level implementation, as operands in the operation register of the instruction can be used directly. More specifically, the (classic) MPC controller needs to solve a nonlinear optimization problem at each sampling instant, and if the associated computational complexity is O(MIN_ITER) with MIN_ITER being the minimum number of iterations to find a solution to the nonlinear optimization, the computational complexity for the developed triggering mechanism is only O(1), and therefore, the complexity of the triggering mechanism is minor and would not add much computational burden to the robot system. Besides, introducing the event-triggering mechanism to the MPC can effectively reduce the number of the sampling and updating instants, and thus, the entire computation resources will actually be significantly saved compared with its time-triggered counterpart.
Interpretations of Algorithm 1: (1) e initialization parameters of the algorithm include prediction horizon N p , initial state ξ 0 � x 0 y 0 φ 0 T , triggering condition level σ (including threshold curve and threshold band), and weighting matrices Q and R Mathematical Problems in Engineering 5 (2) When t k � 0, the optimization problem in (14) is solved to obtain the control sequence ΔU * t k � [Δu * t k , Δu * t k +1 , ... , Δu * t k +N c −1 ] T , which uses the first element to form the control input u(t k ) � u(t k − 1) + ΔU * t k (1), where u(t k − 1) � 0, and u(t k ) is applied to the system to obtain the new state vector ξ(i + 1|t k ) (3) e triggering conditions are checked via substituting ξ(i + 1|t k ) into equations (15) and (17) to see whether the predetermined thresholds are met (4) When the triggering conditions are not valid, the elements in the input sequence are continuously applied without computing the optimization; otherwise, the optimization problem (14) is solved for the new control sequence, and the new input is utilized for getting the new state vector; then, the algorithm is back to step (3).
In order to quantify the computation and communication resources saved by using the event-triggering mechanism, S 1 and S 2 are defined as where T e and T t are the event-triggering times and timetriggering times, and d is the measured mean network delay in a standard experimental environment. erefore, equations (20) and (21) can quantificationally calculate the saved computation and communication resources. It is worth nothing that a large amount of historical data is normally collected when the robot is controlled by the classical MPC before using the above strategy. erefore, in view of different working conditions, the control strategy can also be used after statistical analysis of the previous historical data. erefore, the control strategies proposed in this study also have a strong adaptability.

Simulation and Experimental Verifications
To verify the effectiveness of the proposed threshold curve/ band-based MPC method, the simulation test and experimental verification are carried out, respectively.

Simulation Test.
e simulation is performed on the MATLAB platform. Given the differential motion model of the mobile robot in equation (1) e threshold curve and band are then obtained by using the threshold selection method (equations (16) and (18)) proposed in Section 4, and the obtained values are shown in Figure 1.
en, the control effect is compared with the standard time-triggering strategy, and the results are shown in Figure 2, in which TT denotes the time-triggered method, while ET1 and ET2 denote the event-triggering strategies using the threshold curve and threshold band methods, respectively. Figure 2 illustrates that the proposed two event-triggering strategies have similar tracking performance compared with the standard time-triggered method. In fact, we can see that the event-triggering strategies both have faster response speed, indicating the advantage of the event-triggered controller with some practical and meaningful condition. In addition, the proposed method has been compared with an existing event-triggered MPC approach developed in literature [18], and the results show that event-triggering strategy proposed in this study can achieve better tracking performance and response speed (as illustrated in Figure 2).
Besides, Table 1 provides the effect of changing the design parameters on the closed-loop performance of the robot system for the proposed method. From Table 1, it is observed that weighting matrices Q and R are mutually restricted, and by reducing R (or increasing Q), the rapid response performance of the system can be improved; by increasing R (or reducing Q), the control cost can be reduced and the steady-state performance of the system can be improved. As for the prediction horizon N p and control horizon N c , due to physical limitations, N c <� N p must be satisfied, and in order to maintain the degree of freedom for the MPC optimization, N c � N p is normally considered in real applications. Given such a setting, it is observed that a larger prediction horizon N p /control horizon N c can provide a faster response speed but a certain compromise in steady-state performance will also be resulted. Figure 3 shows the tracking error of state components x, y, and φ from the corresponding reference. Specifically, Figure 3 shows that the event-triggering strategies have faster response speeds than their time-triggering counterparts, and the threshold band method also has a faster response speed than the threshold curve method. Note that due to the existing of the external disturbance, the state components are changing slightly around the origin, rather than remain constants.
Furthermore, the triggering instants are shown in Figure 4, in which the triggered time steps of ET1 and ET2 are marked with red and blue circles with value 1. Note that the time-triggered MPC needs to update at every single sampling time (250 times in total), and the curve-based event-triggering strategy needs 184 times, and the band-based strategy only

Experimental Verification.
en, the experiment test is carried out. Figure 5 shows the utilized crawler mobile robot and the experimental environment. e robot body is driven by two DC servo motors equipped with rotary encoders. e   reference trajectory is a diagonal line within an 8 m × 8 m rectangle. e operating system on the robot is the Indigo version of the robot operating system (ROS). In order to establish communication with the same version of the ROS laptop, the 54 M bandwidth 802.11 g wireless communication protocol is used to realize remote wireless connection and control the movement of the robot. After network testing, the average network communication delay is about 17.04 ms, as shown in Figure 6.
Due to the error between the expected velocity command transmitted by the MPC node to the motion controller and the actual velocity measured by the encoder, it is necessary to calibrate the velocity of the mobile robot before the experiment. Let the robot travel a distance of 3 m at velocities of 0.1, 0.2, and 0.3 m/s, and each velocity was measured three times, and the average velocity is recorded, based on which the error can be calculated, and the detail data are shown in Table 2.
e initial position is [0, 0, 0], the target position is [8, 8, π/4], the prediction and control horizons are N p � N c � 20, the controller sampling frequency is 20 Hz, the input constraint is −0.3 ≤ u ≤ 0.3, the disturbance limit is ρ � 0.05, the group number for constructing the triggering condition is N � 10, and other parameters are the same as the simulation. e upper bound of disturbance is selected as the maximum steady-state error in the experiment. e real-time position and speed information of the mobile robot is obtained through wireless communication, and the ROS visualization (Rviz) tool and stageros tool in the host computer are used to record the real-time trajectory, as shown in Figure 7.
In addition, coordinate values recorded by a mobile robot in the global coordinate system are imported into MATLAB. According to equations (16) and (18), the threshold curve and threshold band can be calculated, as shown in Figure 8. rough the recorded data, we can get the tracking error of each state component under timetriggered and event-triggered methods, as shown in Figure 9.
It is observed that compared with the time-triggered method, the maximum error of each state component is  0.0268, 0.0141, and 0.0157 m for the curve-based eventtriggering strategy, and the maximum error is 0.0258, 0.0333, and 0.0275 m for the band-based strategy, which illustrates that the developed event-triggering mechanism has a similar control effect as the standard method. Besides, from the 3rd subplot of Figure 9, the robot mainly completes the orientation adjustment process at 1.9 s and then gradually reduces the horizontal and vertical coordinate errors. Figure 10 illustrates the corresponding triggering instants. Note that the time-triggered method needs to

Conclusion
is study proposes a novel MPC-based trajectory tracking control method based on the event-triggering mechanism for mobile robot. e method includes two event-triggering strategies, the triggering strategy based on the threshold curve and the threshold band, respectively. e construction of the threshold curve and the threshold band is to employ the statistical method to process the recorded historical trajectory data. e simulation result and experimental verification show that compared with the classic timetriggered MPC, the proposed method can significantly reduce the update times of the optimization problem such that the computation and communication burdens can be effectively reduced without affecting the control performance of the system. In the future, the method to optimally select the design parameters in the event-triggering conditions and the way to handle more complex working conditions, e.g., other types of model uncertainties, time delays, and cyber attacks, are considered as potential research directions, and the target is to enhance the practicability of the eventtriggered control technology. Data Availability e main contribution of the work is the proposed two new event-triggering strategies, and thus, the data are not the main supporting material.

Conflicts of Interest
e authors declare that they have no conflicts of interest.