Next Article in Journal
Acceptability, Feasibility, and Effectiveness of Immersive Virtual Technologies to Promote Exercise in Older Adults: A Systematic Review and Meta-Analysis
Next Article in Special Issue
Parametric Dynamic Distributed Containment Control of Continuous-Time Linear Multi-Agent Systems with Specified Convergence Speed
Previous Article in Journal
BlockEdge: A Privacy-Aware Secured Edge Computing Framework Using Blockchain for Industry 4.0
Previous Article in Special Issue
Indoor Location Technology with High Accuracy Using Simple Visual Tags
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Model Predictive Control for Mobile Robots with Localization Fluctuation Estimation

1
Intelligent Transportation Systems Research Center, Wuhan University of Technology, 1178 Heping Avenue, Wuhan 430000, China
2
Chongqing Research Institute, Wuhan University of Technology, 598 Liangjiang Avenue, Chongqing 400000, China
3
Hubei Institute of Measurement and Testing Technology, 2 Maodianshan Middle Road, Wuhan 430000, China
4
The State Key Laboratory of New Textile Materials and Advanced Processing Technologies, Wuhan Textile University, 1 Yangguang Avenue, Wuhan 430000, China
5
School of Mechanical Science and Engineering, Huazhong University of Science and Technology, 1037 Luoyu Road, Wuhan 430000, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(5), 2501; https://doi.org/10.3390/s23052501
Submission received: 31 December 2022 / Revised: 12 February 2023 / Accepted: 16 February 2023 / Published: 23 February 2023
(This article belongs to the Special Issue Advanced Sensing and Control Technologies for Autonomous Robots)

Abstract

:
Mobile robots are widely employed in various fields to perform autonomous tasks. In dynamic scenarios, localization fluctuations are unavoidable and obvious. However, common controllers do not consider the impact of localization fluctuations, resulting in violent jittering or poor trajectory tracking of the mobile robot. For this reason, this paper proposes an adaptive model predictive control (MPC) with an accurate localization fluctuation assessment for mobile robots, which balances the contradiction between precision and calculation efficiency of mobile robot control. The distinctive features of the proposed MPC are three-fold: (1) Integrating variance and entropy—a localization fluctuation estimation relying on fuzzy logic rules is proposed to enhance the accuracy of the fluctuation assessment. (2) By using the Taylor expansion-based linearization method—a modified kinematics model that considers that the external disturbance of localization fluctuation is established to satisfy the iterative solution of the MPC method and reduce the computational burden. (3) An improved MPC with an adaptive adjustment of predictive step size according to localization fluctuation is proposed, which alleviates the disadvantage of a large amount of the MPC calculation and improves the stability of the control system in dynamic scenes. Finally, verification experiments of the real-life mobile robot are offered to verify the effectiveness of the presented MPC method. Additionally, compared with PID, the tracking distance and angle error of the proposed method decrease by 74.3% and 95.3%, respectively.

1. Introduction

Mobile robots are being progressively used in numerous scenarios such as unmanned factories, logistics centers, and exhibition halls, thanks to their superior flexibility and maneuverability [1,2,3]. For unmanned operations, autonomous navigation technology is intuitively important for robots [4]. To follow a given trajectory, mobile robots have to be able to control their pose precisely and robustly based on the localization results [5]. However, localization and control issues are often studied independently, leaving robot control performers to be improved.
The control system needs accurate localization results as a reference to maintain good trajectory-tracking accuracy [6]. In traditional control methods, kinematic or dynamic modelling or the control theory have been given more attention, and localization results are always seen as an absolute truth value [7,8]. However, in practice, mobile robots, whether using vision-based or LiDAR-based localization solutions, are subject to noise interference from their external sensors, which may lead to fluctuations in localization [9,10]. If the controller still treats the localization result as absolute truth, this can lead to severe jittering or large tracking deviations. To this end, this paper solves the robust control problem of mobile robots under the localization fluctuations, including the accurate localization fluctuation estimation and the robust controller design.
Through the above analysis, to improve the operation accuracy of the four-wheel differential mobile robot, this paper analyses the localization fluctuation state in the dynamic scene and then realizes the adaptive adjustment of the predictive step size of the model predictive control. Therefore, the operation adaptability of the four-wheel differential mobile robot is improved.
(1)
Integrating variance and information entropy—an enhanced localization fluctuation estimation method based on fuzzy logic rules is proposed to improve the accuracy of the fluctuation assessment.
(2)
A modified kinematics model with external disturbance using the Taylor expansion-based linearization is established, which is convenient for controller design under localization fluctuations.
(3)
An improved MPC with an adaptive adjustment of predictive step size related to localization fluctuation is proposed, which ensures the stability of the control system in dynamic scenes.
(4)
The proposed method has been tested in real dynamic scenarios and compared with mainstream methods, and its effectiveness has been demonstrated.

2. Related Works

2.1. Localization Fluctuation Estimation

The problem of localization can be divided into simultaneous localization and map building (SLAM) and localization based on an a priori map, depending on the presence or absence of an a priori map [11]. Mainstream SLAM methods such as ORB-SLAM [12] and LEGO-LOAM [13] can generate a map of the environment while obtaining localization results. However, this method is subject to cumulative errors and has poor real-time performance. In contrast, a priori map-based localization provides accurate and efficient positional information and is often used as a reference for control [11]. Bayesian filter-based localization frameworks are currently the dominant approach using a priori map, such as Kalman filtering or Monte Carlo localization (MCL) [14,15,16]. In particular, MCL is widely used due to its ability to adapt to non-Gaussian non-linear scenarios [14,17]. Although there is much research around the robustness of MCL, existing algorithms are not immune to localization fluctuations in highly dynamic scenarios. The accurate description of the localization fluctuations is of great significance for the design of subsequent navigation systems. For this reason, Zapata et al. propose to use the maximum particle weight in the MCL as a benchmark to determine the current localization reliability [18]. When the maximum particle weight is less than the weight threshold, this is an indication that the currently estimated pose is not reliable. Nevertheless, it is not robust to use the weight of only one particle to measure localization reliability. In turn, variance and entropy values are common in addition to valid metrics for estimating localization fluctuations, which take into account the set of particles with weights in an integrated manner [19,20,21]. The variance is well understood mathematically and physically—the larger the variance, the greater the localization fluctuation—but is poorly described for localization data with multi-peaked distributions. Higher information entropy indicates smaller differences in particle weights, demonstrating greater uncertainty in localization. In particular, entropy is more accurate for non-convex data evaluations [22,23]. However, relying solely on a numerical metric is prone to misclassification. In this regard, the main objective of this paper is to design a robust assessment method that integrates more localization fluctuation metrics into consideration.

2.2. Mobile Robot Control

As one of the key modules, control technology has a great impact on the stability and accuracy of mobile robots [24]. In the process of motion, it is usually subject to external disturbances such as model uncertainty and parameter perturbation, resulting in motion oscillation and deviation, and even skidding and rollover [25]. To suppress disturbances, many scholars have made efforts to improve the performance of mobile bots. For example, the literature [26] designs a sliding mode control (SMC) scheme based on a reduced-order-extended-state observer, which realizes the active compensation of friction under the condition of uncertain parameters and ensures the stability of an operation. Through the smooth fitting of the path and design of MPC, the literature [27] realizes the high-speed movement of the four-wheeled independently steering robot with action delay. For the tracking problem of wheeled mobile robots with bounded disturbances and various practical constraints, a robust MPC method is proposed to ensure the safety and comfort during an operation [28]. At present, SMC and MPC have been the focus of research because of their excellent characteristics in high-speed and high-precision motion control methods [6,29,30]. Although SMC has the advantage of being insensitive to disturbances, the oscillation produced by itself is difficult to be eliminated, which makes it difficult to be widely used [31]. The gradual iterative optimization is brought by MPC, which can well handle the model constraints caused by structure, dynamic system, etc., and is conducive to achieving smooth motion [32,33]. This advantage improves the robustness of the control. In the design process of existing MPC methods, the observed localization data are usually treated as accurate values, which can ensure the stability of motion in static- or high-localization accuracy scenes [34,35]. Therefore, in the dynamic scene, the design of the controller needs to take into account the localization fluctuations, so as to avoid causing motion oscillation. For this purpose, how to improve the robustness and accuracy of the controller in the localization fluctuation scenario has become a key issue to be studied in our work.

3. System Modelling and Problem Formulation

3.1. System Modelling

Figure 1 shows the four-wheel differential platform model. The four-wheel differential platform has good motion performance and can achieve zero radius turning by adjusting the speed of the left and right wheels, which improves its adaptability to complex scenes. As the special case of mobile robots, the general modelling method can improve the general adaptability of the model [36,37]. Therefore, to further analyze the four-wheel differential platform and improve its motion controllability, the following general kinematics model of a mobile robot is given:
χ = f ( χ , u )
where χ = [ x , y , θ ] T is the state variable; u = [ v , ω ] T represents the control variable, x and y are the position of the mobile robot center point in the global fixed coordinate system, and θ denotes the robot heading angle; v and ω are the linear velocity and angular velocity of the mobile robot, respectively. The Taylor formula is used to expand the nonlinear mobile robot model at the reference point ( χ r , u r ) to obtain the linear model of the mobile robot, so that the modelling accuracy can be guaranteed.
At the same time, the linear model reduces the amount of computation, convenient for the design of the controller and is conducive to the actual implementation. Furthermore, at the Taylor expansion for Equation (1) at the point ( χ r , u r ) , we have:
χ ˙ = f ( χ r , u r ) + f χ r ( χ χ r ) + f u r ( u u r ) + O r
where χ r is the reference control input; u r denotes the calculated reference input, f χ r ( ) and are coefficient matrices of function f ( ) expanded at ( χ r , u r ) ; O ( r ) is the higher order remainder of the Taylor expansion. By defining e = χ χ r , we have
e ˙ = f χ r e + f u r u ˜ + O r
where e represents the following error of the mobile robot; u ˜ = u u r is the change of input control law. In the actual motion process, it is difficult to realize the continuous control of the mobile robot. Therefore, obtaining a discrete-time model of the four-wheel mobile robot is necessary. Setting the sampling period as T , where
e ( k + 1 ) = e ( k ) + T e ˙ ( k )
where k is the sampling time. The above model is rewritten as
e ( k + 1 ) = A ( k ) e ( k ) + B ( k ) u ˜ ( k ) + O r ( k )
where,
A ( k ) = 1 0 T v r sin θ r ( k ) 0 1 T v r cos θ r ( k ) 0 0 1 , B ( k ) = T cos θ r ( k ) 0 T sin θ r ( k ) 0 0 T
Therefore, in order to improve the controllability of the four-wheel mobile robot, we realized the linear modelling of the mobile robot kinematics model. This facilitates the design of the controller and makes the control of the mobile robot simpler.

3.2. Localization Problem Formulation

The mobile robot localization problem is often regarded as a typical Bayesian estimation problem. Bayesian filter-based localization algorithms solve robot localization problems by estimating the probability distribution of robot poses in the pose space and assigning a probability to each possible hypothetical pose using a confidence level (Belief), which is expressed as
B ( s k ) = p ( s k | o 1 : k , u 1 : k , M )
where s k is the robot’s pose in the two-dimensional plane at time k , which can be expressed as x k , y k , θ k ; x k and y k are the robot’s position, and θ k is the robot’s heading; o 1 : k and u 1 : k represent the sensor observations and motion control from the initial time to time k , respectively; M is the a priori map.
According to the Markov hypothesis and Bayes rule, the Bayesian filter-based localization method can recursively estimate the robot’s poses, but it involves a large number of integral operations and nonlinear non-Gaussian features of the observation and motion models, which make the efficient and accurate solution a concern. To this end, this paper uses the MCL approach to compute Equation (7), which uses a particle set with weights to represent B ( s k ) , i.e.,
B ( s k )   n p = 1 N p ω k [ n p ] δ ( s k s k [ n p ] )
where the set of particles with weights is denoted as s k [ n p ] , ω k [ n p ] n p = 1 N p ; ω k [ n p ] is the weight of the n p -th particle s k [ n p ] and N p means the total number of particles; s k [ n p ] can be written as x k [ n p ] , y k [ n p ] , θ k [ n p ] ; x k [ n p ] , y k [ n p ] and θ k [ n p ] are the position and orientation of the particle; δ ( · ) is the Dirichlet function. Usually, the particle with the highest weight is selected as the current localization result.

4. Localization Fluctuations Estimation

Both variance and information entropy have their own unique advantages in expressing localization fluctuations and can reflect the characteristics of the particle set. However, a single performance metric still has a large randomness that affects the accuracy of the fluctuation assessment. To solve the above problems, a fuzzy logic rule incorporating variance and entropy is proposed to evaluate the localization fluctuations. Localization fluctuations are represented as follows:
L f = [ L f x , L f y , L f θ ] T = [ f V x ( V x ) + f E ( E ) , f V y ( V y ) + f E ( E ) , f V θ ( V θ ) + f E ( E ) ] T
f V i ( V i ) = α V 0 V i < η 1 V T i β V η 1 V T i V i < η 2 V T i λ V η 2 V T i V i ,   f E ( E ) = α E 0 E < η 3 E T β E η 3 E T E < η 4 E T λ E η 4 E T E
where V i and E are the variance and entropy, respectively; i = x , y , θ ; f V i ( V i ) and f E i ( E i ) are the localization fluctuation factors based on V i and E ; α V , β V , λ V , α E , β E and λ E are the fluctuation parameters; 0 < α V < β V < λ V and 0 < α E < β E < λ E ; η 1 , η 2 , η 3 and η 4 are the weight coefficient; 0 < η 1 < η 2 and 0 < η 3 < η 4 ; V T i and E T are the variance and entropy threshold value, severally. L f is a three-dimensional vector [ L f x , L f y , L f θ ] T that represents the fluctuations of x, y, and θ . Meanwhile, L f is determined by the variance and entropy values, where the entropy values are only related to the particle weights, so all three dimensions are set uniformly.
It is worth stating that the calculation of V i and E usually requires a series of localization results over a period of time. However, it is not practical to perform a large number of localization experiments in situ to determine the current localization fluctuation state, and we would prefer to conduct the evaluation depending on the localization data at a certain time. Fortunately, the MCL result is expressed as a particle set s k [ n p ] , ω k [ n p ] n p = 1 N p , so that V i and E can be represented as
V p x = i = 1 N p ω k [ n p ] ( x k n p x ¯ k ) 2 ,   V p y = i = 1 N p ω k [ n p ] ( y k n p y ¯ k ) 2 ,   V p θ = i = 1 n ω t i ( θ t i θ ¯ t ) 2
E p = n p = 1 N p ω k [ n p ] log ω k [ n p ]
where V p i and E p are the variance and entropy from s k [ n p ] , ω k [ n p ] n p = 1 N p ; ¯ denotes the mean value of the related particle’s pose.
Although V i and E can be reflected by the localization results at a certain time based on s k [ n p ] , ω k [ n p ] n p = 1 N p , the exact relationship between V p i , V i , E p and E are still difficult to be indicated analytically. Fuzzy logic rules are an effective means to infer an output based on input variables. Hence, we integrate V p i and V i into a fuzzy formula using the following fuzzy logic rules:
rule   1 :   i f 0 V p i V p i 1 t h e n f V i 1 = f V i ( V i ) s . t . 0 V i < η 1 V T i rule   2 :   i f V p i 2 V p i V p i 3 t h e n f V i 2 = f V i ( V i ) s . t . η 1 V T i V i < η 2 V T i rule   3 :   i f V p i 4 V p i t h e n f V i 3 = f V i ( V i ) s . t . η 2 V T i V i
where V p i 1 , V p i 2 , V p i 3 and V p i 4 are fuzzy demarcation boundaries for V p i ; f V i 1 , f V i 3 and f V i 3 are the fluctuation values based on variance.
Similarly, the fuzzy formula for describing the mapping relation between E p and E can be written as
rule   1 :   i f 0 E p E p 1 t h e n f E 1 = f E ( E ) s . t . 0 E < η 3 E T rule   2 :   i f E p 2 E p E p 3 t h e n f E 2 = f E ( E ) s . t . η 3 E T E < η 4 E T rule   3 :   i f E p 4 E p t h e n f E 3 = f E ( E ) s . t . η 4 E T E
where E p 1 ,   E p 2 ,   E p 3 and E p 4 are fuzzy demarcation boundaries for E p ; f E 1 ,   f E 3 and f E 3 are the fluctuation values based on entropy.
For implementation, the boundaries V p i 1 to V p i 4 and E p 1 to E p 4 can be learned from a mass of variance and entropy of localization results in different dynamic environments. As the next step of a standard procedure of constructing a fuzzy logic system, defuzzification is achieved by the weighted average method. Additionally, the fluctuation factor L f is obtained as follows:
L f = [ L f x , L f y , L f θ ] T = m = 1 r χ m f V x m m = 1 r χ m + m = 1 r χ m f E m m = 1 r χ m , m = 1 r χ m f V y m m = 1 r χ m + m = 1 r χ m f E m m = 1 r χ m , m = 1 r χ m f V θ m m = 1 r χ m + m = 1 r χ m f E m m = 1 r χ m T
where r is the number of fuzzy rule bases; χ m denotes the trigger strength of m-th rule. If the V p i or E p satisfies the fuzzy rule, then χ m is 1. Otherwise, χ m is set as 0.

5. Adaptive MPC Considering Localization Fluctuation

In order to realize the robust control of the four-wheel mobile robot, based on the state space discrete model of Equation (5), we construct the cost function required in the optimal control process
J e ( k ) , u ˜ ( k ) = i = 1 N p e ( k + i | k ) T P e ( k + i | k ) + i = 0 N c 1 u ˜ ( k + i | k ) T R [ u ˜ ( k + i | k ) ]
where J ( ) is the cost function; e ( k + i | k ) and u ˜ ( k + i | k ) represent the error value and control fluctuation value of time k + 1 predicted at time k , respectively; N p and N c are prediction and control horizon, respectively, and the value of the control horizon is not greater than the prediction horizon; P and R are the weight matrix.
To ensure the feasibility of optimal prediction, the control variables of the prediction process are as follows:
e min e ( k ) e max u min u ( k ) u max u ˜ min u ˜ ( k ) u ˜ max u min u ( k ) + T u ˜ ( k ) u max O r min O r ( k ) O r max
where e min and e max are the minimum and maximum errors, respectively. u min and u max are the minimum and maximum values of the control law increment, respectively. O r min and O r max are the minimum and maximum perturbations, respectively.
Considering the existence of localization fluctuation, to ensure the stability of operation, N p and N c need to be further adjusted as follows:
N p = [ k 1 max ( L f ) ] + k p N c = [ k 2 max ( L f ) ] + k c
where L f is the estimated localization fluctuation value obtained from Equation (9), max ( ) is the maximum value function of a vector, k 1 , 2 is the adjustment coefficient. k c , k p N + is the minimum adjustment coefficient. [ X ] is the maximum integer value not greater than X . Considering the constraint state of the predictive control method, there are N c N p .
Furthermore, using the cost function of Equation (16), the following state equation of the prediction horizon is obtained:
e ¯ ( k ) = F ( k ) u ˜ ¯ ( k ) + L ( k ) e ( k ) + G ( k ) u ˜ ( k 1 ) + O ¯ r ( k )
with
e ¯ ( k ) = [ e ( k + 1 | k ) , , e ( k + N p | k ) ] T u ˜ ¯ ( k ) = [ u ˜ ( k + 1 | k ) , , u ˜ ( k + N c 1 | k ) ] T O ¯ r ( k ) = [ O r ( k + 1 | k ) , , O r ( k + N p 1 | k ) ] T
The coefficient matrix is expressed as:
F ( k ) = B ( k ) 0 A ( k + 1 ) B ( k ) 0 B ( k + N p 1 ) + i = k + 1 i = k + N p 1 A ( i ) B ( k ) B ( k + N p 1 ) + i = k + 1 i = k + N p 1 A ( i ) B ( k + N c )
L ( k ) = A ( k ) , A ( k + 1 ) A ( k ) , , i = k i = k + N p 1 A ( i ) T
G ( k ) = B ( k ) B ( k + 1 ) + A ( k + 1 ) B ( k ) B ( k + N p 1 ) + + i = k + 1 i = k + N p 1 A ( i ) B ( k )
Therefore, in order to ensure the optimal operation process, the following optimization objectives are obtained by combining Formulas (16), (17) and (19):
min e ¯ T ( k ) P e ¯ ( k ) + u ˜ ¯ T ( k ) R u ˜ ¯ ( k )
Considering modelling error and disturbance, Equation (19) is introduced into the optimization objective (24), thus:
min O ¯ r T ( k ) P O ¯ r ( k ) + u ˜ ¯ T ( k ) R u ˜ ¯ ( k )
subject to
min ( e max M ( k ) , O ¯ r max ) O ¯ r ( k ) max ( e min M ( k ) , O ¯ r min ) u min u ( k ) u max u ˜ min u ˜ ( k ) u ˜ max   u min u ( k ) + T u ˜ ( k ) u max
where M ( k ) = F ( k ) u ˜ ¯ ( k ) L ( k ) e ( k ) G ( k ) u ˜ ( k 1 ) is the intermediate variable.
By adjusting the above equation, the control law optimization under the scenario of disturbance fluctuation is realized.
Furthermore, the stability proof of the designed controller is given as follows:
Theorem 1.
When the optimal control strategy of (19)–(25) is adopted and the following conditions are satisfied:
e ( k + N p | k ) = 0
u ˜ ( k + N c ) = 0
Then, the system will approach stability.
Proof. 
The optimization function of the system can be rewritten as:
J e ( k + 1 ) , u ˜ ( k + 1 ) = i = 1 N p e ( k + i + 1 | k + 1 ) T P e ( k + i + 1 | k + 1 ) + i = 0 N c 1 u ˜ ( k + i + 1 | k + 1 ) T R u ˜ ( k + i + 1 | k + 1 ) = i = 1 N p e * ( k + i + 1 | k ) T P e * ( k + i + 1 | k ) + i = 0 N c 1 u ˜ * ( k + i + 1 | k ) T R u ˜ ( k + i + 1 | k )
Combined with Formulas (27) and (28), there is
J e ( k + 1 ) , u ˜ ( k + 1 ) = i = 1 N p e * ( k + i | k ) T P e * ( k + i | k ) + i = 0 N c 1 u ˜ * ( k + i | k ) T R u ˜ * ( k + i | k ) e * ( k + 1 | k ) T P e * ( k + 1 | k ) u ˜ * ( k | k ) T R u ˜ * ( k | k ) = J * e ( k ) , u ˜ ( k ) e * ( k + 1 | k ) T P e * ( k + 1 | k ) u ˜ * ( k | k ) T R u ˜ * ( k | k )
with J e ( k + 1 ) , u ˜ ( k + 1 ) , which is the optimal control value at time k . Then, e ( k + i | k ) and u ˜ ( k + i | k ) are the predicted values obtained by solving the optimization function at time k to k + i . Hence, combining Equations (27) and (28), we have:
e ( k + 1 ) = e ( k + 1 | k ) ,   u ˜ ( k ) = u ˜ ( k | k )
Then,
J e ( k + 1 ) , u ˜ ( k + 1 ) J e ( k + 1 ) , u ˜ ( k + 1 ) = J * e ( k ) , u ˜ ( k ) [ e * ( k + 1 ) ] T P [ e * ( k + 1 ) ] [ u ˜ ( k ) ] T R [ u ˜ * ( k ) ] J e ( k ) , u ˜ ( k )
Therefore, the optimal value of the cost function will gradually decrease in the process of iterative updating. It shows that the optimization method can guarantee the asymptotic convergence of the system. □

6. Experimental Validations

6.1. Experimental Implementation

The experimental scenes and four-wheel differentially driven mobile robot are shown in Figure 2. The mobile robot basically consists of an industrial computer (Intel(R) Core (TM) i7-6500U CPU @2.50 GHz, 8 GB of RAM, 64-bit operating system), two LiDARs, four motor encoders and some related sensors, such as an ultrasonic transducer and anti-collision strip. More specifically, two UTM-30LX 2-D LiDARs with a range of 30 m and scanning rate of 40 Hz are used, which guarantee that the robot has enough field of view to ensure safety and real-time pose tracking. To avoid the blind spot of LiDAR, we installed four MaxBotix MB7360 ultrasonic sensors around the robot with a high resolution of 1 mm and a measuring distance of 5 m, which is sufficient to detect possible dynamic obstacles during the robot’s movement. As demonstrated in Figure 2b, the experimental environment is a dynamic and wide square surrounded by overgrown grass and a large number of pedestrians. Three sites are selected for the localization experiments to determine the fuzzy rule parameters. To comprehensively demonstrate the performance of the proposed control system, we first conduct localization experiments in different dynamic environments. Then, through real control experiments, the effectiveness of the proposed MPC method is verified.

6.2. Experimental Results and Discussions

6.2.1. Experimental Results of Localization Fluctuation Estimation

To construct fuzzy logic rules for localization fluctuation estimation, 100 frames of LiDAR measurement are collected at each of the three sites in different dynamic scenarios using a mobile robot. In total, 900 frames of LiDAR data were provided for the fuzzy estimation. The LiDAR data from low to high dynamic scenes at Site 2 in the square are shown in Figure 3a. The LiDAR data from low to high dynamic scenes at Site 2 in the square are shown in Figure 1. Specifically, as shown in Figure 1, the low dynamic scenes have less than 10 percent of the observed noise, which is ideal for localization. The medium dynamic scenario is depicted in Figure 3b, where the observed noise accounts for 20 to 30 percent of the overall observed data, a common and realistic scenario. Extremely dynamic scenarios, such as the one shown in Figure 3c, will have noise levels close to 50 percent, which can easily lead to localization fluctuations or even localization failure.
Next, we determine the relationship between V p and E , and the mapping between E p and E . We put a different number of obstacles around the robot and record 100 localization results. V i and E correspond to each localization result at three sites. In this way, as indicated in Table 1, we get the V i and E of 100 localization results at three sites in different dynamic environments. And, as the dynamic level of the scene increases, V i and E become correspondingly larger.
According to the V i and E in different dynamic scenes, we can know the status of the localization fluctuations to which V p and E belong. Meanwhile, for example, when i = x, Equation (10) can be set to
f V x ( V x ) = 0.00 0.00 V x < 1.93 × 10 4 0.25 1.93 × 10 4 V x < 3.23 × 10 4 0.5 3.23 × 10 4 V x ,   f E ( E ) = 0 0.00 E < 6.37 0.25 6.37 E < 6.47 0.5 6.47 E
Specifically, the relationship between V p i and V i , and the mapping between E p and E , are shown in Table 2. We have grouped V p i and E p , which fit into different fluctuation ranges, into one category. The 300 V p i and V i ranges are linked. Similarly, the mapping of the ranges of E and E p is found.
The distribution of V p i and E p in the different ranges of V i and E are shown in Figure 4. We have used box plots to illustrate the distribution characteristics. It is easy to see that there are overlapping parts of V p i and E p in different ranges, which provide for the construction of fuzzy rules. According to Figure 4, the boxes in different ranges have a distinct overlap. Hence, we can easily choose the fuzzy partition boundary V p i 1 V p i 4 and E p 1 E p 4 for V p i and E p , respectively. Finally, the “premise” of the fuzzy rule base defines V p i and E p , and the fluctuation values f V i 1 f V i 3 and f E 1 f E 3 are the “conclusion”. Equations (13) and (14), when i = x, can be rewritten as
rule   1 :   i f 0.00 V p x 5 × 10 4 t h e n f V x 1 = f V x ( V x ) s . t . 0.00 V x < 1.93 × 10 4 rule   2 :   i f 1.62 × 10 4 V p x 6.07 × 10 4 t h e n f V i 2 = f V x ( V x ) s . t . 1.93 × 10 4 V x < 3.23 × 10 4 rule   3 :   i f 4.94 × 10 4 V p x t h e n f V i 3 = f V i ( V x ) s . t . 3.23 × 10 4 V x
rule   1 :   i f 0.00 E p 5.21 t h e n f E 1 = f E ( E ) s . t . 0.00 E < 6.37 rule   2 :   i f 4.46 E p 6.27 t h e n f E 2 = f E ( E ) s . t . 6.37 E < 6.47 rule   3 :   i f 6.04 E p t h e n f E 3 = f E ( E ) s . t . 6.47 E

6.2.2. Experimental Results of Adaptive MPC

Furthermore, to verify the superiority of the control method considering location uncertainty, we choose the following comparison methods: (1) The traditional PID control method with optimizing parameters, where K p = 0.5 , K i = 1.5 and K d = 0.1 . (2) The proposed control method does not consider the localization fluctuation (NMPC) with N c = N p = 5 . In the proposed control method, the control parameters are set as follows: k p = k c = 5 , k 1 = k 2 = 10 . The specific tracking process of the comparison method is described as follows:
The trajectory tracking of various comparison methods is shown in Figure 5. The initial point of the four-wheel differential mobile robot is (10, 0). It can be seen from Figure 5 that there is an error between the attitude angle of the initial point and the slope of the curve, resulting in a large tracking error. With time adjustment, all control methods can achieve fast convergence. It can be seen from the enlarged drawing of Figure 5 that compared with the proposed method and NMPC method, PID has a large oscillation. Further analysis shows that the proposed method has a better tracking effect.
Figure 6 and Figure 7 are the distance error and angle error, respectively. Specifically, Figure 7 shows that the overshoots of the PID and NMPC are 0.5903 m and 0.5502 m, while the proposed method is 0.4603 m. Furthermore, by comparing the proposed method with NMPC, it is found that the tracking process is optimized through adaptive step adjustment. On closer inspection, the average errors in the stable stage (10 s~60 s) of PID, NMPC and proposed methods are 0.0374 m, 0.0177 m and 0.0096 m, which show that the errors are reduced by 74.3% and 45.8%, respectively. It is found in 0 that the angle error of the PID method fluctuates greatly, while the NMPC and the proposed method will be limited to a small range. The average angle errors of PID and NMPC are 0.1001 rad and 0.0062 rad during a 10~60 s stabilization period, while the proposed method with a small error is 0.0047 rad. The stable errors are reduced up to 95.3% and 31.9%, respectively. By comparison, the proposed method achieves better error suppression.
From Figure 8, the adjustment of the prediction horizon is given through the localization fluctuation output by the localization module to ensure the control’s stability and accuracy under the high localization fluctuation scenario. Therefore, compared with the traditional PID and NMPC methods, the proposed method can realize the dynamic adjustment of the step size, improve tracking accuracy and reduce motion fluctuation.

7. Conclusions and Outlook

In this paper, an adaptive MPC for mobile robots with the fuzzy estimation of localization fluctuations was discussed, which accurately evaluates the localization state and effectively improves control accuracy and robustness. First of all, a localization fluctuation estimation based on fuzzy logic rules was proposed, which takes into account both the effects of variance and information entropy, thereby obtaining accurate estimates of localization fluctuations. Then, a modified kinematics model with external disturbance using the Taylor expansion-based linearization was constructed. In addition, according to localization fluctuation, an enhanced MPC with an adaptive adjustment of predictive step size was proposed to maintain the stability of the control system in dynamic scenes. The experimental results show that the step size of MPC can be effectively and adaptively adjusted with localization fluctuation in different dynamic scenes, avoiding the generation of oscillation. In the stable tracking phase, the proposed method reduced the tracking distance error by 74.3% and 45.8% compared with PID and NMPC. In addition, in comparison to PID and NMPC, the tracking angle errors of the proposed method decreased by 95.3% and 31.9%, respectively. In future work, we will explore the potential of deep learning for applications in localization fluctuation assessments and improve the control framework to better integrate localization fluctuation information.

Author Contributions

Conceptualization, J.M. and N.J.; methodology, J.M., N.J. and L.J. (Liquan Jiang); software, H.X. and L.J. (Liyu Jiang); validation, L.J. (Liyu Jiang) and Z.H.; formal analysis, L.J. (Liyu Jiang); investigation, H.X.; resources, J.M. and N.J.; data curation, L.J. (Liquan Jiang); writing—original draft preparation, J.M. and N.J.; writing—review and editing, J.M.; visualization, L.J. (Liquan Jiang) and Z.H.; supervision, N.J.; project administration, N.J.; funding acquisition, N.J. and J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by the Natural Science Foundation of Hubei Province, China under Grant no. 20221j0060, in part by the Fund of National Engineering Research Center for Water Transport Safety under Grant no. A202303, in part by the Natural Science Foundation of Chongqing, China under Grant no. CSTB2022NSCQ-MSX1566 and in part by Hubei Provincial Engineering Research Center for Intelligent Textile and Fashion (Wuhan Textile University) under Grant no. 2022HBITF01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Notations and Abbreviations

The following key notations and abbreviations are used in this paper.
Notations/AbbreviationsDescriptions
SMCsliding mode control
MPCmodel predictive control
J ( ) cost function
N p , N c Prediction/control horizon
P , R weight matrix
v linear velocity
ω angular velocity
f χ r ( ) , f u r ( ) coefficient matrices
O ( r ) higher order remainder of the Taylor expansion
P l estimated localization fluctuation value
k 1 , 2 adjustment coefficient
k c , k p N + minimum adjustment coefficient

References

  1. Rubio, F.; Valero, F.; Llopis-Albert, C. A review of mobile robots: Concepts, methods, theoretical framework, and applications. Int. J. Adv. Robot. Syst. 2019, 16, 1729881419839596. [Google Scholar] [CrossRef] [Green Version]
  2. Skoczeń, M.; Ochman, M.; Spyra, K.; Nikodem, M.; Krata, D.; Panek, M.; Pawłowski, A. Obstacle detection system for agricultural mobile robot application using RGB-D cameras. Sensors 2021, 21, 5292. [Google Scholar] [CrossRef] [PubMed]
  3. Meng, J.; Wang, S.; Jiang, L.; Xie, Y.; Zheng, S.; Wu, H. Robust lateral stabilization control of in-wheel-motor-driven mobile robots via active disturbance suppression approach. Sensors 2020, 20, 5238. [Google Scholar] [CrossRef]
  4. Higgins, J.; Bezzo, N. Negotiating visibility for safe autonomous navigation in occluding and uncertain environments. IEEE Robot. Autom. Lett. 2021, 6, 4409–4416. [Google Scholar] [CrossRef]
  5. Zhang, X.; Xie, Y.; Jiang, L.; Li, G.; Meng, J.; Huang, Y. Fault-tolerant dynamic control of a four-wheel redundantly-actuated mobile robot. IEEE Access 2019, 7, 157909–157921. [Google Scholar] [CrossRef]
  6. Jiang, L.; Wang, S.; Xie, Y.; Xie, S.; Zheng, S.; Meng, J.; Ding, H. Decoupled Fractional Supertwisting Stabilization of Interconnected Mobile Robot Under Harsh Terrain Conditions. IEEE Trans. Ind. Electron. 2021, 69, 8178–8189. [Google Scholar] [CrossRef]
  7. Tzafestas, S.G. Mobile robot control and navigation: A global overview. J. Intell. Robot. Syst. 2018, 91, 35–58. [Google Scholar] [CrossRef]
  8. Jiang, L.; Wang, S.; Xie, Y.; Xie, S.Q.; Zheng, S.; Meng, J. Fractional robust finite time control of four-wheel-steering mobile robots subject to serious time-varying perturbations. Mech. Mach. Theory 2022, 169, 104634. [Google Scholar] [CrossRef]
  9. Meng, J.; Wan, L.; Wang, S.; Jiang, L.; Li, G.; Wu, L.; Xie, Y. Efficient and reliable LiDAR-based global localization of mobile robots using multiscale/resolution maps. IEEE Trans. Instrum. Meas. 2021, 70, 1–15. [Google Scholar] [CrossRef]
  10. Zhang, M.; Liu, X.; Xu, D.; Cao, Z.; Yu, J. Vision-based target-following guider for mobile robot. IEEE Trans. Ind. Electron. 2019, 66, 9360–9371. [Google Scholar] [CrossRef]
  11. Meng, J.; Wang, S.; Xie, Y.; Li, G.; Zhang, X.; Jiang, L.; Liu, C. A safe and efficient LIDAR-based navigation system for 4WS4WD mobile manipulators in manufacturing plants. Meas. Sci. Technol. 2021, 32, 045203. [Google Scholar] [CrossRef]
  12. Mur-Artal, R.; Montiel, J.M.M.; Tardos, J.D. ORB-SLAM: A versatile and accurate monocular SLAM system. IEEE Trans. Robot. 2015, 31, 1147–1163. [Google Scholar] [CrossRef] [Green Version]
  13. Shan, T.; Englot, B. Lego-loam: Lightweight and ground-optimized lidar odometry and mapping on variable terrain. In Proceedings of the 2018 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Madrid, Spain, 1–5 October 2018; pp. 4758–4765. [Google Scholar]
  14. Chung, M.-A.; Lin, C.-W. An improved localization of mobile robotic system based on AMCL algorithm. IEEE Sens. J. 2021, 22, 900–908. [Google Scholar] [CrossRef]
  15. Li, G.; Meng, J.; Xie, Y.; Zhang, X.; Huang, Y.; Jiang, L.; Liu, C. Reliable and fast localization in ambiguous environments using ambiguity grid map. Sensors 2019, 19, 3331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Meng, J.; Wang, S.; Jiang, L.; Hu, Z.; Xie, Y. Accurate and Efficient Self-localization of AGV Relying on Trusted Area Information in Dynamic Industrial Scene. In IEEE Transactions on Vehicular Technology; IEEE: Manhattan, NY, USA, 2023. [Google Scholar]
  17. Ge, G.; Zhang, Y.; Wang, W.; Jiang, Q.; Hu, L.; Wang, Y. Text-MCL: Autonomous mobile robot localization in similar environment using text-level semantic information. Machines 2022, 10, 169. [Google Scholar] [CrossRef]
  18. Zhang, L.; Zapata, R.; Lepinay, P. Self-adaptive Monte Carlo localization for mobile robots using range finders. Robotica 2012, 30, 229–244. [Google Scholar] [CrossRef] [Green Version]
  19. Sun, H.; Wang, S.; Meng, J.; Liu, Y.; Xie, Y. Accurate Pose Tracking of Mobile Robot Using Entropy-based TrimICP in Dynamic Environment. In Proceedings of the IECON 2022–48th Annual Conference of the IEEE Industrial Electronics Society, Brussels, Belgium, 17–20 October 2022; pp. 1–6. [Google Scholar]
  20. Bukhori, I.; Ismail, Z.H. Detection of kidnapped robot problem in monte carlo localization based on the natural displacement of the robot. Int. J. Adv. Robot. Syst. 2017, 14, 1729881417717469. [Google Scholar] [CrossRef] [Green Version]
  21. Meng, J.; Wang, S.; Li, G.; Jiang, L.; Zhang, X.; Liu, C.; Xie, Y. Iterative-learning error compensation for autonomous parking of mobile manipulator in harsh industrial environment. Robot. Comput.-Integr. Manuf. 2021, 68, 102077. [Google Scholar] [CrossRef]
  22. Toomaj, A.; Di Crescenzo, A. Connections between weighted generalized cumulative residual entropy and variance. Mathematics 2020, 8, 1072. [Google Scholar] [CrossRef]
  23. Zidek, J.V.; Van Eeden, C. Uncertainty, entropy, variance and the effect of partial information. Lect. Notes-Monogr. Ser. 2003, 42, 155–167. [Google Scholar]
  24. Chen, Y.; Li, Z.; Kong, H.; Ke, F. Model predictive tracking control of nonholonomic mobile robots with coupled input constraints and unknown dynamics. IEEE Trans. Ind. Inform. 2018, 15, 3196–3205. [Google Scholar] [CrossRef]
  25. Singh, P.; Nandanwar, A.; Behera, L.; Verma, N.K.; Nahavandi, S. Uncertainty compensator and fault estimator-based exponential supertwisting sliding-mode controller for a mobile robot. In IEEE Transactions on Cybernetics; IEEE: Manhattan, NY, USA, 2021. [Google Scholar]
  26. Ren, C.; Li, X.; Yang, X.; Ma, S. Extended state observer-based sliding mode control of an omnidirectional mobile robot with friction compensation. IEEE Trans. Ind. Electron. 2019, 66, 9480–9489. [Google Scholar] [CrossRef]
  27. Liu, X.; Wang, W.; Li, X.; Liu, F.; He, Z.; Yao, Y.; Ruan, H.; Zhang, T. MPC-based high-speed trajectory tracking for 4WIS robot. ISA Trans. 2022, 123, 413–424. [Google Scholar] [CrossRef] [PubMed]
  28. Dai, L.; Lu, Y.; Xie, H.; Sun, Z.; Xia, Y. Robust tracking model predictive control with quadratic robustness constraint for mobile robots with incremental input constraints. IEEE Trans. Ind. Electron. 2020, 68, 9789–9799. [Google Scholar] [CrossRef]
  29. Sun, Z.; Xia, Y.; Dai, L.; Liu, K.; Ma, D. Disturbance rejection MPC for tracking of wheeled mobile robot. IEEE/ASME Trans. Mechatron. 2017, 22, 2576–2587. [Google Scholar] [CrossRef]
  30. Jiang, L.; Xie, Y.; Jiang, Z.; Meng, J.; Li, W. Adaptive model predictive control of mobile robot with local path refitting. In Proceedings of the 2022 IEEE 17th Conference on Industrial Electronics and Applications (ICIEA), Chengdu, China, 16–19 December 2022; pp. 596–601. [Google Scholar]
  31. Li, J.; Wang, J.; Peng, H.; Hu, Y.; Su, H. Fuzzy-torque approximation-enhanced sliding mode control for lateral stability of mobile robot. IEEE Trans. Syst. Man Cybern. Syst. 2021, 52, 2491–2500. [Google Scholar] [CrossRef]
  32. Ding, T.; Zhang, Y.; Ma, G.; Cao, Z.; Zhao, X.; Tao, B. Trajectory tracking of redundantly actuated mobile robot by MPC velocity control under steering strategy constraint. Mechatronics 2022, 84, 102779. [Google Scholar] [CrossRef]
  33. Zhang, Y.; Zhao, X.; Tao, B.; Ding, H. Point stabilization of nonholonomic mobile robot by Bézier smooth subline constraint nonlinear model predictive control. IEEE/ASME Trans. Mechatron. 2020, 26, 990–1001. [Google Scholar] [CrossRef]
  34. Hu, Y.; Su, H.; Fu, J.; Karimi, H.R.; Ferrigno, G.; De Momi, E.; Knoll, A. Nonlinear model predictive control for mobile medical robot using neural optimization. IEEE Trans. Ind. Electron. 2020, 68, 12636–12645. [Google Scholar] [CrossRef]
  35. Taheri, H.; Zhao, C.X. Omnidirectional mobile robots, mechanisms and navigation approaches. Mech. Mach. Theory 2020, 153, 103958. [Google Scholar] [CrossRef]
  36. Wang, J.; Luo, Z.; Wang, Y.; Yang, B.; Assadian, F. Coordination control of differential drive assist steering and vehicle stability control for four-wheel-independent-drive EV. IEEE Trans. Veh. Technol. 2018, 67, 11453–11467. [Google Scholar] [CrossRef]
  37. Cui, M.; Liu, W.; Liu, H.; Jiang, H.; Wang, Z. Extended state observer-based adaptive sliding mode control of differential-driving mobile robot with uncertainties. Nonlinear Dyn. 2016, 83, 667–683. [Google Scholar] [CrossRef]
Figure 1. The four-wheel model and single-track model of the considered IWMD-MR.
Figure 1. The four-wheel model and single-track model of the considered IWMD-MR.
Sensors 23 02501 g001
Figure 2. Experimental scene and platform. (a) Platform prototype. (b) Dynamic scene. (c) Grid map of the scene.
Figure 2. Experimental scene and platform. (a) Platform prototype. (b) Dynamic scene. (c) Grid map of the scene.
Sensors 23 02501 g002
Figure 3. LiDAR measurement in different dynamic scenes at site 2. (a) Low dynamic environment. (b) Medium dynamic environment. (c) High dynamic environment.
Figure 3. LiDAR measurement in different dynamic scenes at site 2. (a) Low dynamic environment. (b) Medium dynamic environment. (c) High dynamic environment.
Sensors 23 02501 g003
Figure 4. The distribution of V p and E p in the different ranges of V and E . (a) V p x in the different ranges of V x . (b) V p y in the different ranges of V y . (c) V p θ in the different ranges of V θ . (d) E p in the different ranges of E .
Figure 4. The distribution of V p and E p in the different ranges of V and E . (a) V p x in the different ranges of V x . (b) V p y in the different ranges of V y . (c) V p θ in the different ranges of V θ . (d) E p in the different ranges of E .
Sensors 23 02501 g004
Figure 5. Tracking response of the comparison method.
Figure 5. Tracking response of the comparison method.
Sensors 23 02501 g005
Figure 6. Distance error of the comparison method.
Figure 6. Distance error of the comparison method.
Sensors 23 02501 g006
Figure 7. Orientation error of the comparison method.
Figure 7. Orientation error of the comparison method.
Sensors 23 02501 g007
Figure 8. Adjustment of the prediction horizon.
Figure 8. Adjustment of the prediction horizon.
Sensors 23 02501 g008
Table 1. The results of the V i and E of multiple localization results in different dynamic environments.
Table 1. The results of the V i and E of multiple localization results in different dynamic environments.
  Low Dynamic SceneMedium Dynamic SceneHigh Dynamic Scene
  Vx/(m2)Vy/(m2)Vθ/(rad2)EVx/(m2)Vy(m2)Vθ/(rad2)EVx/(m2)Vy(m2)Vθ/(rad2)E
Site 11.24 × 10−43.23 × 10−41.53 × 10−66.332.19 × 10−48.55 × 10−41.93 × 10−66.414.67 × 10−46.28 × 10−37.48 × 10−56.56
Site 21.62 × 10−45.71 × 10−45.57 × 10−76.322.09 × 10−47.45 × 10−41.04 × 10−66.434.33 × 10−41.26 × 10−32.70 × 10−66.51
Site 31.80 × 10−42.63 × 10−41.48 × 10−66.372.66 × 10−49.09 × 10−44.76 × 10−66.393.45 × 10−41.83 × 10−38.97 × 10−66.54
Table 2. The results of V p i and E p in different ranges.
Table 2. The results of V p i and E p in different ranges.
RangeRange1:
0.00 V x < 1.93 × 10 4
Range2:
1.93 × 10 4 V x < 3.23 × 10 4
Range3:
3.23 × 10 4 V x
V p x /(m2)7.03 × 10−43.06 × 10−45.05 × 10−41.17 × 10−23.54 × 10−43.92 × 10−41.18 × 10−19.10 × 10−45.84 × 10−4
98 lines of V p x
4.76 × 10−43.37 × 10−44.57 × 10−43.96 × 10−42.88 × 10−41.49 × 10−41.39 × 10−37.85 × 10−48.83 × 10−4
Rangerange1:
0.00 V y < 6.11 × 10 4
range2:
6.11 × 10 4 V y < 1.98 × 10 3
range3:
1.98 × 10 3 V y
V p y /(m2)8.01 × 10−44.08 × 10−41.45 × 10−42.61 × 10−23.24 × 10−49.01 × 10−42.82 × 10−11.20 × 10−21.37 × 10−3
98 lines of V p x
2.75 × 10−41.86 × 10−42.69 × 10−42.18 × 10−43.21 × 10−41.40 × 10−41.12 × 10−38.20 × 10−41.14 × 10−3
Rangerange1:
0.00 V θ < 1.88 × 10 6
range2:
1.88 × 10 6 V θ < 1.57 × 10 5
range3:
1.57 × 10 5 V θ
V p θ /(rad2)4.59 × 10−61.03 × 10−63.50 × 10−73.14 × 10−48.30 × 10−79.34 × 10−75.13 × 10−31.20 × 10−41.99 × 10−5
98 lines of V p θ
7.66 × 10−62.45 × 10−63.01 × 10−69.88 × 10−72.12 × 10−61.69 × 10−61.09 × 10−42.11 × 10−55.95 × 10−5
Rangerange1:
0.00 ≤ E <6.37
range2:
6.37 ≤ E <6.47
range3:
6.47 ≤ E
E p 6.222.515.185.056.305.001.118.857.79
98 lines of E p
4.885.034.075.735.395.197.297.987.94
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Meng, J.; Xiao, H.; Jiang, L.; Hu, Z.; Jiang, L.; Jiang, N. Adaptive Model Predictive Control for Mobile Robots with Localization Fluctuation Estimation. Sensors 2023, 23, 2501. https://doi.org/10.3390/s23052501

AMA Style

Meng J, Xiao H, Jiang L, Hu Z, Jiang L, Jiang N. Adaptive Model Predictive Control for Mobile Robots with Localization Fluctuation Estimation. Sensors. 2023; 23(5):2501. https://doi.org/10.3390/s23052501

Chicago/Turabian Style

Meng, Jie, Hanbiao Xiao, Liyu Jiang, Zhaozheng Hu, Liquan Jiang, and Ning Jiang. 2023. "Adaptive Model Predictive Control for Mobile Robots with Localization Fluctuation Estimation" Sensors 23, no. 5: 2501. https://doi.org/10.3390/s23052501

APA Style

Meng, J., Xiao, H., Jiang, L., Hu, Z., Jiang, L., & Jiang, N. (2023). Adaptive Model Predictive Control for Mobile Robots with Localization Fluctuation Estimation. Sensors, 23(5), 2501. https://doi.org/10.3390/s23052501

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop