Abstract

In this paper, a trajectory tracking control algorithm is proposed based on the fractional-order PD (FOPD) controller for a Wheeled Mobile Robot (WMR). Firstly, an improved flat phase property is put forward as a robust controller tuning specification. This specification is capable of guaranteeing the flatness of the phase curve in a frequency interval, so the controlled system robustness can be improved. Then, the stabilization process is discussed with respect to the parameters of the FOPD controller through a visualized 3-dimensional surface, so both the stability and robustness of the controlled system can be guaranteed under the proposed controller. Furthermore, the implementation of the proposed robust FOPD controller is presented, which makes the control algorithm easy to be realized. At last, the effectiveness of the proposed trajectory tracking control algorithm is verified by the simulation and experiment results.

1. Introduction

Wheeled Mobile Robots (WMRs) are capable of working in different situations, including the inclement, dangerous, or even harmful ones. In the past few decades, WMRs have been widely applied in a great variety of civilian and military tasks, for example, space exploration, materials transportation, supplies delivery, mine clearance, and search and rescue [16]. Currently, trajectory tracking and regulation is one of the most concerned problems in WMRs related studies [7]. Without an effective control strategy, a predefined tracking strategy is hard to follow, especially in long distance or complicated environment tasks. The nonholomic properties, internal dynamics, feedback sensors of WMRs, and external load disturbance may bring in different kinds of immeasurable uncertainties [8]. Therefore, more precise and robust trajectory tracking strategy will certainly help in improving the operation efficiency of WMRs.

Typically, Proportional-Integral-Derivative (PID) controller is always used in the control process of industrial robotics including WMRs. But its control effect has been suspicious at times when better robustness and transient performance are required. However, the combination of fractional calculus and traditional PID controller, namely, fractional-order PID (FOPID) controller, provides novel potential and opportunity of solving this kind of problem. Many studies have investigated the design and application of FOPID type controllers [911]. One of the most representative works is the controller proposed by Podlubny, which is an extension of the traditional PID controller with two extra order parameters [10]. Another kind of FOPID controller designed based on phase and magnitude margin frequency specifications is presented by Vinagre et al. [12]. Except for phase and magnitude margin specifications, Monje et.al. added the extrasensitivity and complementary sensitivity functions in the design specifications of the FOPID controller [13]. A similar FOPI controller whose parameters were tuned by symmetrical optimization function was studied by Maione and Lino [14]. In [15, 16], the original and modified Ziegler–Nichols methods were used in the design process of FOPID controllers. The Newton–Raphson recursive algorithm was applied by Feliu et al. in the parameter optimization process of an FOPI controller [17]. Moreover, some intelligent control algorithm-based FOPID tuning methods have also been studied [1822]. An autotuning variable-order fuzzy FOPID controller is investigated in [18] by Liu et al. in which all the parameters can be tuned online to deal with random time delay or system parameter uncertainties. The detailed advantages of the FOPID type controller can be concluded as following. Firstly, the FOPID type controller inherits the characteristics of the traditional PID controller, such as simple structure, clear physical, and meaning. Secondly, through the extradifferential and integral order parameters, the design flexibility of the controller is increased so better dynamic performance and robustness may be achieved with this type of controller. Subsequently, the FOPID type controller has the memory property and can adjust the output of the controller properly following the error history information in order to obtain better performance. At last, the added fractional-order terms can help adjust the high-frequency and low-frequency characteristics of the closed-loop easily [23].

In addition, the practical applicability of the FO controller has been verified in some practical applications. Several FO control examples, including industrial plants and electromechanic actuator, are presented in [24]. A brief summary and introduction of FOPID controllers applied in various kinds of industrial situations is given in [25]. The CRONE controller which can compensate the disturbance caused by parameter uncertainties and load changes has been successfully used in vehicle suspension [26] and path tracking [27]. A series of robust FOPID controllers tuned based on the flat phase property have been used in motion control [28, 29]. Moreover, in [30], an FO controller is used in the trajectory tracking task of a differential drive WMR. The experimental results show that the control performance of the FO controller is more satisfactory compared with the traditional PD controller. Another speed and direction control task accomplished by the FOPI controller for skid-steered WMR is presented in [31], which also achieves superior control performance. Nevertheless, the related studies are still quite limited at present. The stability of the controlled systems and the advantages of FO controllers applied in speed regulation, trajectory tracking, and other practical applications of WMRs need more exploration.

The main contribution of this paper includes the design, implementation, and validation of an FOPD controller for precise trajectory tracking of WMRs. For robust FO controller regulation, the flat phase property has been widely used as a design specification [29, 32]. This is because it can guarantee the flatness of the phase curve around the interested crossover frequency , so the system phase margin can remain constant and the system is robust to load variations. However, there are mainly two drawbacks of the original flat phase property. The first one is that it can only guarantee the flatness of the phase curve at one frequency, namely, . The other one is there are several nonlinear equations to be solved in the controller design process, so it is hard to find a solution. The improved flat phase property proposed in this paper can guarantee the robust flat phase frequency interval width instead of only one flat phase frequency and at least one parameter set which satisfies all the specifications can be found. In addition, the system stabilization and implementation problems are also discussed. The simulation and experiments results are shown to verify the effectiveness of the proposed trajectory tracking control algorithm.

The rest of this paper is organized as follows: In Section 2, the modelling of XQ WMR is given; then, the fractional calculus and FOPD controller are presented in Section 3; Section 4 shows the proposed controller design specifications; Section 5 discusses the stabilizing process for FO delay systems with the FOPD controller; furthermore, the experiment results and discussion are illustrated in Section 6 to verify the effectiveness and flexibility of the proposed controller; finally, the conclusion is drawn in Section 7.

2. Modelling of WMR

The XQ WMR studied in this paper is a three-wheeled differential drive robot, as shown in Figure 1.

2.1. Kinematic Model of XQ WMR

There are two front drive wheels and a universal follower wheel on the chassis of XQ WMR. The two front wheels driven by DC motors provide the forward power, so the movement and orientation of XQ WMR can be accomplished by the velocity difference between the two drive wheels; the universal follower wheel is mounted on the rear of the chassis. The kinematic model of XQ WMR is illustrated in Figure 2. As it is shown in Figure 2, denote as the chassis linear, right wheel, and left wheel velocities, respectively, θ as the angular displacement and ω as the angular velocity, where , and as the distance between the two drive wheels. Denote the state vector S of XQ WMR as follows:where are coordinates of the center mass of XQ WMR. Therefore, the forward kinematic model of XQ WMR can be obtained as follows:

Then, in order to achieve the input desirable right and left wheel velocities, the inverse kinematic model of XQ WMR is achieved as

The control scheme proposed in this paper is shown in Figure 3. Note that the centre of the two drive wheels O is assumed to be the mass centre of XQ WMR.

2.2. Motor Dynamic Model of XQ WMR

Before controller design procedure, the motor dynamic model should be achieved first. Figure 4 is the schematic diagram of the DC motor, where and are the total resistance and total inductance of the armature winding; is the armature voltage; is the electromagnetic torque; is the rotational angular velocity of the motor; and is the rotational inertia on the deceleration shaft.

The transfer function of XQ WMR DC motor can be achieved as follows:where and are constants related to the magnetic flux. Ignoring some of the negligible values and the parameters which have less impact on the overall mathematical model, a simplified DC mathematical model can be obtained as a second-order transfer function as follows:where the time constant is and .

From the experimental identification, the parameters of the DC motor are achieved as .

3. Fractional Calculus and FOPD Controller

3.1. Fractional-Order Derivative

Fractional calculus, which is an extension of traditional calculus, has not got a unified definition so far. There are three definitions which have been extensively used [10, 33], namely, Grunwald–Letnikov definition, Riemann–Liouville definition, and Caputo definition. Each of the definitions has its own properties. The definitions should be applied appropriately in different research fields, such as engineering, applied mathematics, and computer science. The initial values of Caputo definition are the same with that of the integer-order system and own practical physical meaning [34, 35]. Therefore, the Caputo definition is used in this paper.

An FO integral-differential operator can be depicted as follows:where is the integral or differential order and are lower and upper limits of the FO operator, respectively.

The Caputo derivative of order δ for a function is defined as follows[33]:where , and m is a positive integer.

The Laplace transformation corresponding to Caputo definition can be obtained as follows:where s is the Laplace transformation operator and represents the transformation.

3.2. The Considered Fractional-Order PD Controller

Different from the traditional PD controller, there are three parameters in the proposed FOPD controller which are formulated as follows:where . This is a specific form of the FOPID controller which involves an integrator order λ ( in this paper) and a differentiator order γ [36, 37].

4. FOPD Controller Design Specifications with Guaranteed Flat Phase Interval

The open-loop transfer function of WMR can be achieved from the above sections as follows:

Based on the definitions of gain and phase margins, the following design specifications are introduced:

4.1. Gain Crossover Frequency Specification

The amplitude of the open-loop system should be zero at its gain crossover frequency in logarithmic frequency domain. In other words, the amplitude should equal to 1 at its gain crossover frequency, which can be expressed as follows:where is the interested crossover frequency.

4.2. Phase Margin Specification

Denote as the required phase margin, and this specification can be depicted as follows:

4.3. Flat Phase Interval Guaranteed Robustness Specification

To realize the robustness for the controlled system, the flat phase robust tuning specification can be described by the following condition at crossover frequency :

However, there are a few drawbacks of the flat phase specification in equation (13), which have already been discussed in Section 1. Here, we give the following condition to improve the existing flat phase specification:where is a small enough scalar and smaller ε means better flatness of the phase curve; given ε, , and are the minimum and maximum frequency near () between which the flat phase interval can be guaranteed. Define as the width of the flat phase with regards to ε. Longer l means better robust performance under the same ε. Next, we present the steps to obtain l.Step 1: given ε, compute whether the following inequality holds or not:

If yes, go to Step 2. If not, it obtains .Step 2: if inequality (3) holds, solve the next equation and gain its solutions:Let and be the solutions which are the nearest to in negative and positive directions, respectively.Step 3: give the width .

Remark 1. All the above steps need to be processed under a given . When is fixed, longer l could offer better robustness for the controlled system. However, an excessively small may lead to . Thus, should be well chosen according to the practical requirements.

5. Design Procedure

Stability is the primary concern in control systems. So, the proposed controller design procedure is finding out the complete parameter sets which can stabilize the controlled system first, and then pick up the parameters which satisfy the proposed design specifications from the complete sets.

5.1. Control Process Stabilizing

The stabilizing process is investigated with respect to parameters of the proposed FOPD controller, namely, .

One single closed-loop transfer function depicted in Figure 3 can be yielded according to equation (10) as follows:

So, the characteristic equation of the XQ WMR motor system is

The controlled system stability is determined by the root locations of its characteristic equation in (11). If all the roots locate in the left half of s-domain, the system is supposed to be bounded-input bounded-output stable. Therefore, the stability region of the XQ WMR motor system can be obtained by seeking for the parameter sets which can ensure the roots of equation (11) lie in the left half of s-domain in . The boundaries of the stability region derived from can be achieved by finding the corresponding IRB (infinite root boundary), CRB (complex root boundary), and RRB (real root boundary) [38]. Note that if a system is strictly proper, it will not have IRB [39]. So, only RRB and CRB are taken into consideration in this paper.(i)RRB is defined by

So that  = 0.(ii)CRB can be expressed as follows:where

Here, both the real part and imaginary part of equation (20) should be equal to 0:where

Then, one can get from equation (21) that

Therefore, with from 0, the stability region of a fixed fractional order γ composed by the corresponding RRB and CRB can be achieved. All parameter sets () in the region can guarantee the stability of the controlled XQ WMR motor system. Then, a three-dimensional surface of () can be determined by sweeping the differential order . is the biggest γ which can ensure that the control system is strictly proper. The achieved surface is the maximum stability surface of the XQ WMR system under the FOPD controller because the phase margin is assumed to be 0 at this moment. Parameter sets () on this surface can ensure the stability of the controlled system; however, no design specification is satisfied at this stage.

Substitute s in equations (10) and (11) by as follows:where

Then, one can get

Multiply both sides of equation (16) with and yields,

Hence, equation (13) is obtained as follows:

Denote

From equation (28), one can obtain

From equations (26) and (30), E can be obtained as follows:

Substituting E in equation (26) by equation (30), are achieved as follows:

When is predefined, with from 0 and from 0, an available stability surface can be obtained from equation (32). Then, if an interested is chosen, an available stability curve can be achieved. At this stage, the design specifications in equations (11) and (12) are fulfilled by parameters () on the available stability curve. Remark that cannot be bigger than , where is determined by the intersection points of RRB and CRB. It is the maximum ω value which can ensure the stability of the control system.

From equations (14) and (28), one can obtainwhere

At last, the proposed FOPD parameters () are determined by checking all the parameter sets on the obtained available stability curve and finding out the one which satisfy the robust design specification in equation (14). So, the obtained FOPD controller satisfies all the three design specifications in Section 3 and can also ensure the stability of the closed-loop system.

On a whole, the design procedure can be summarized as follows:

Step 1:. find the maximum stability surface of XQ WMR under the FOPD controller with from 0 by sweeping . The parameter sets on maximum stability surface can only guarantee the stability of the controlled system.

Step 2:. define phase margin and find the available stability surface. Then, choose the interested crossover frequency , so the available stability curve can be obtained on the available stability surface. The parameter sets on available stability curve satisfy both the gain crossover frequency and phase margin specifications.

Step 3:. check all the parameter sets on the achieved available stability curve and find the one fulfill the robust design specification mostly.

6. Experiment and Discussion

6.1. Implementation of FOPD Controller

The implementation of the FOPD controller is one of the critical problems to achieve desirable trajectory tracking performance. The key point of the FOPD controller implementation concentrates on the approximation of the fractional operator in equation (9). In this section, the widely used Oustaloup et al. approximation method is applied without loss of generality [40].

A fractional operator can be expressed as follows:

Consider the interested approximation frequency interval as , so the operator can be substituted bywhere .

Then, equation (35) can be updated as follows:where .

The transfer function above can be transformed into a zero-pole form as follows:where

The 7th Oustaloup approximation with expected frequency range is used in this section.

6.2. Controller Parameters Regulation

According to the controller design procedures in Section 4, the maximum stability surface is obtained for XQ WMR in Figure 5. So, all the parameter sets on the surface can guarantee the stability of the control system.

Set the phase margin as , the available stability surface can be obtained, as shown in Figure 6. It can be seen that the area of the available stability surface is inversely proportional to the chosen phase margin. Then, is set as the interested crossover frequency, so the available stability curve can be obtained on the available stability surface. The parameters on the achieved curve satisfy the gain crossover frequency and phase margin specifications.

Finally, the parameter sets on the available stability curve are checked according to the robustness specification and the proposed FOPD controller parameters are obtained as . Correspondingly, for fair comparison, a PID controller which is tuned by the optimal tuning method in [41] with similar specifications is also designed here.

Figure 7 illustrates the speed regulation comparison under the FOPD controller and PID controller with step input. It is shown that the speed regulation process under the FOPD controller has smaller rising time, accommodation time, steady state error, and almost no overshoot, which outperforms the control performance under the PID controller.

The robustness and precision comparisons are demonstrated in Figure 8 with different gain variations and external disturbance at . Note that the gain variations here can be regarded as load disturbance of XQ WMR. It is shown that the speed regulation process under the FOPD controller changes little with respect to different kinds of disturbance. However, there are obvious changes in the overshoot and accommodation time of the speed regulation process under the PID controller, which means the XQ WMR system controlled by the PID controller may not achieve precise and robust control performance in practice.

6.3. Experiment

In this section, the trajectory tracking experiment results of XQ WMR are shown to verify the effectiveness of the proposed FOPD controller. For fair comparison, a traditional PD controller, whose proportional and derivative parameters are the same with that of the proposed FOPD controller, is also implemented. The reference signal is depicted as follows:where are the x-axis and y-axis reference trajectories. The sampling time is .

The comparison of XQ WMR trajectory tracking performance is shown in Figure 9. It can be seen that the tracking performances under both controllers are acceptable. However, from the zoomed in areas of Figure 9, it is observed that the tracking performance under the proposed FOPD controller is more precise and the settling time is much shorter, which is also verified by the absolute error and x-axis error shown in Figures 10 and 11 illustrates the left and right speed comparisons of XQ WMR, and the control signal comparisons are demonstrated in Figure 12. These comparisons show that there is smaller oscillation in the speed control process of the proposed controller compared with the traditional PD controller, though the two controllers have the same proportional and derivative parameters.

7. Conclusion

In this paper, a synthesis of the FOPD controller for XQ WMR trajectory tracking control is proposed. Different from other robust FO controller tuning methods, the flat phase property can be guaranteed in a frequency interval instead of one frequency point, so the robustness of the control system is enhanced. The stabilizing process is also studied, so that the proposed FOPD controller can guarantee both the robustness and stability of the controlled system In addition, there is no complicated nonlinear equations which should be solved in this paper, so at least one parameter set which satisfies all the specifications can be found. The implementation and experiment results are presented to show the advantages of the proposed trajectory tracking algorithm. The proposed control algorithm is also capable of fulfilling different control requirements.

Data Availability

The data used to support the findings of this study are currently under embargo while the research findings are commercialized. Requests for data, 6 months after publication of this article, will be considered by the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Nos. 11902252 and 51979229), the Fundamental Research Funds for the Central Universities of China (Nos. G2018KY0305 and G2018KY0302), China Postdoctoral Science Foundation (No. 2019M663811), the Natural Science Foundation of Shaanxi Province (No. 2019JQ-164), and the Opening Foundation of Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University (No. 1817).