Pseudo-State Sliding Mode Control of Fractional SISO Nonlinear Systems

This paper deals with the problem of pseudo-state sliding mode control of fractional SISO nonlinear systems with model inaccuracies. Firstly, a stable fractional sliding mode surface is constructed based on the Routh-Hurwitz conditions for fractional differential equations. Secondly, a sliding mode control law is designed using the theory of Mittag-Leffler stability. Further, we utilize the control methodology to synchronize two fractional chaotic systems, which serves as an example of verifying the viability and effectiveness of the proposed technique.


Introduction
Fractional calculus has a long history of three hundred years, over which a firm theoretical foundation has been established.In the past few decades, with deep understanding of the power of fractional calculus and rapid development of computer technology, an enormous number of interesting and novel applications have emerged in physics, chemistry, engineering, finance, and other sciences [1][2][3].In particular, engineers and scientists from various fields have developed plenty of fractional dynamical systems, that is, systems which are better characterized by noninteger order mathematics models.As has been stated in [3], a suitable way to more efficient control of fractional dynamical systems is to design fractional controllers.
Several pioneering attempts to develop fractional control methodologies have been made, such as TID controller [4], CRONE controller [5], fractional PID controller [6], and fractional lead-lag compensator [7].Basic ideas and technical formulations of the above four fractional control schemes with comparative comments have been presented in [8].
Very recently, by applying fractional calculus to advanced nonlinear control theory, several fractional nonlinear control schemes have been proposed, such as fractional sliding mode control, fractional adaptive control, and fractional optimal control.Exactly, to design fractional sliding mode controls, various fractional sliding surfaces have been constructed in [9][10][11][12][13][14][15][16][17][18][19][20][21].In particular, adaptive sliding mode controls have been proposed in [12,[19][20][21], fractional terminal sliding mode controls in [15,16], and sliding mode controls for linear fractional systems with input and state delays in [17].In [22], the authors have presented two ideas to extend the conventional model reference adaptive control (MRAC) by using fractional parameter adjustment rule and fractional reference model.In [23], a fractional model reference adaptive control algorithm for SISO plants has been proposed, which can guarantee the stability and ability to reject disturbances.In [24][25][26], adaptive controllers have been designed to control and synchronize fractional chaotic systems.In [27], Agrawal has proposed a general formulation for a class of fractional optimal control problems, which is specialized for a system with quadratic performance index subject to a fractional system dynamic constraint.In [28], the authors have generalized the optimality conditions of [27].Besides the above fractional control methodologies, fractional optimal synergetic control has been presented in [29] and active disturbance rejection control in [30].
Motivated by the above contributions, this paper proposes a sliding mode control design for fractional SISO nonlinear systems in the presence of model inaccuracies.By constructing a stable fractional sliding mode surface on the basis of Routh-Hurwitz conditions, a sliding mode control law is designed.Further, stability analysis is performed using Mittag-Leffler stability theory.Comparing this with methods Advances in Mathematical Physics in the previous papers, we utilize the fractional derivative of the sliding mode surface instead of first-order derivative, to obtain the equivalent control law.Moreover, to carry out the stability analysis of the closed-loop fractional nonlinear system, we use the fractional derivative of the Lyapunov function candidate in terms of Theorem 2 in [18].
The rest of the paper is organized as follows.Section 2 reviews some basic definitions for fractional calculus.Section 3 proposes the sliding control design for fractional SISO nonlinear systems.Numerical simulations of synchronization of the fractional Genesio-Tesi system and the fractional Arneodo system are presented in Section 4. Finally, Section 5 concludes this paper with some remarks on future study.

Basic Definitions for Fractional Calculus
Fractional calculus is a generalization of integration and differentiation to noninteger order fundamental operator     , where  and  are the bounds of the operation and  ∈ .The continuous integrodifferential operator is defined as [1] The three most frequently used definitions for fractional calculus are the Grünwald-Letnikov definition, the Riemann-Liouville definition, and the Caputo definition [1][2][3].
Definition 1.The Grünwald-Letnikov derivative definition of order  is described as Definition 2. The Riemann-Liouville derivative definition of order  is described as However, applied problems require definitions of fractional derivatives allowing the utilization of physically interpretable initial conditions, which contains (),   (), and so forth.Unfortunately, the Riemann-Liouville approach fails to meet this practical need.It is M. Caputo who solved this conflict.

Definition 3. The Caputo definition of fractional derivative can be written as
In the rest of the paper, we use the Caputo approach to describe the fractional systems and the Grünwald-Letnikov approach to propose numerical simulations.To simplify the notation, we denote the fractional-order derivative of order  as   instead of 0    in this paper.

Sliding Control of Fractional SISO Nonlinear Systems
Consider fractional SISO nonlinear systems where  ∈ (0, 1] is the order of the dynamic system, x = ( 1 ,  2 . . .,   )  ∈ R  is denoted as state vector, and the scalar ,  ∈  are system's input and output, respectively.The dynamic (x) (possibly nonlinear or time varying) is not exactly known, but estimated as f(x).The control gain (x) (possibly time varying or state dependent) is an unknown function.
Assumption 4. The estimation error on (x) is assumed to be bounded by some known function (x): Assumption 5. (x) is assumed to be bounded by where both  min (x) and  max (x) are known functions.
In this paper, our goal is to design a suitable fractional sliding mode controller to make the output  of the system (5) track the desired signal   ; that is, Now we are ready to give the design steps.Firstly, the tracking error is defined as Then, the tracking error vector is defined as A fractional sliding surface is proposed as where the polynomial () =  1 +  2  + ⋅ ⋅ ⋅ +  −1  −2 +  −1 satisfies Routh-Hurwitz conditions for fractional differential equations [31].
Advances in Mathematical Physics 3 The fractional sliding surface (11) can be rewritten as Taking its th order derivative with respect to time yields Substituting the second equation of system ( 5) into ( 13), one has In terms of Assumption 5, one of the estimated values of (x) can be chosen as Bounds of () are written as where  = √ max / min .Design the control law as where and  is a positive constant.
To ensure the stability of the fractional system (5), we have the following theorem.Theorem 7.Under Assumptions 4-6, the fractional SISO nonlinear system (5) can be controlled using the sliding mode control law (17) with the fractional sliding surface (11).
Proof.Consider the following candidate Lyapunov function: where () is the fractional sliding surface (11) constructed previously.
Taking its th order derivative with respect to time along with the fractional sliding surface (11), one has where  is assumed to be an arbitrarily large positive constant which is a bound on the series of (19).Substituting ( 14) into (19) yields Substituting the control law ( 17) into (20) gives Substituting the estimation error ( 6), the equivalent control law û , and control gain  of ( 17) into (21), one derives which implies that the sliding mode dynamic is globally stable and the tracking error vector (10) converges to zero according to Theorem 2 in [18].This proves that the fractional SISO nonlinear system (5) can be controlled using the sliding mode control law (17) with the fractional sliding surface (11).

Numerical Simulations
In this section, we apply the fractional sliding mode control method proposed in Section 3 to deal with the problem of   synchronization between the fractional Arneodo system and the fractional Genesio-Tesi system.To carry out numerical simulations, we utilize the algorithm for numerical calculation of fractional derivatives introduced in [1].This method is derived from the Grnwald-Letnikov Definition 1 based on the fact that the three Definitions 1, 2, and 3 are equivalent for a wide class of functions.
The fractional Arneodo system is represented as with the desired tracking signal   =  1 .
The fractional Genesio-Tesi system is described as Figure 5: Synchronization of the fractional Arneodo system and the fractional Genesio-Tesi system with the control input (28) (red line represents the trajectories of the drive system, while blue line represents the trajectories of the response system).
where (), (), (), and () are assumed to be unknown parameters and satisfy |() + 6| ≤ 0.1, |() + 2.92| ≤ 0.1, |() + 1.2| ≤ 0.1, and |() − 1| ≤ 0.1, from which one has Both of the fractional Arneodo system (23) and the fractional Genesio-Tesi system (24) exhibit chaos with certain values of fractional order .In the sequel we investigate synchronization between the two fractional chaotic systems using the technique proposed previously.The former system is taken as drive system, while the latter one is considered as response system.
The control law is designed as Initial conditions for the above two systems are, respectively, chosen as  0 = (−0.1,0.5, 0.2) and  0 = (−0.2,0.5, 0.2), while parameters in the control law (27) are selected as  = 0.1 and  = 0.5.Numerical simulations of synchronization between the fractional chaotic systems (23) and (24) are presented in Figures 1, 2, 3, and 4, with the simulation time  sim = 10 and time step ℎ = 0.0005.For interpretations of the references to color in these figure legends, the reader is referred to the online version of this paper.
From the simulation results, we see that synchronization performance is excellent but is obtained at the price of high control chattering.It can be eliminated by replacing the discontinuous switching control law sign (⋅) by the smooth function tanh (⋅) in (27); that is, Numerical simulations with the modified control law (28) are presented in Figures 5, 6, 7, and 8.For interpretations of the references to color in these figure legends, the reader is referred to the online version of this paper.
From the above simulation results, one can easily see that the fractional Arneodo system and the fractional Genesio-Tesi system can be effectively synchronized via the proposed sliding mode control technique.Furthermore, the control   chattering caused by the discontinuous control law (17) is successfully eliminated by the modification of (28).

Concluding Remarks
In this paper, we have investigated the pseudo-state sliding control design for fractional SISO nonlinear systems with model inaccuracies.A stable fractional sliding mode surface has been constructed based on the Routh-Hurwitz conditions for fractional differential equations.Then, a sliding mode control law is designed using the Mittag-Leffler stability theorem.Finally, numerical simulations of synchronization of the fractional Arneodo system and the fractional Genesio-Tesi system have been performed to demonstrate the effectiveness of the proposed control technique.
As for the future perspectives, our research activities will be on, (i) designing adaptive sliding control to deal with parametric uncertainties in (⋅), (ii) generalizing the method to fractional MIMO nonlinear systems, (iii) generalizing the method to incommensurate nonlinear systems.

Figure 1 :
Figure1: Synchronization of the fractional Arneodo system and the fractional Genesio-Tesi system with the control input(27) (red line represents the trajectories of the drive system, while blue line represents the trajectories of the response system).

Figure 2 :Figure 3 :
Figure 2: Synchronization errors with the control input (27) (red line represents the first error state  1 , blue line represents the second error state  2 , and black line represents the third error state  3 ).

Figure 6 :
Figure 6: Synchronization errors with the control input (28) (red line represents the first error state  1 , blue line represents the second error state  2 , and black line represents the third error state  3 ).

Figure 7 :
Figure 7: Time history of fractional sliding mode with the control input (28).