Journal of Control Science and Engineering
Volume 2016 (2016), Article ID 4873083, 12 pages
http://dx.doi.org/10.1155/2016/4873083

Research Article

Optimal Control Problem Investigation for Linear Time-Invariant Systems of Fractional Order with Lumped Parameters Described by Equations with Riemann-Liouville Derivative

V. A. Kubyshkin and S. S. Postnov

V. A. Trapeznikov Institute of Control Sciences of Russian Academy of Sciences, Profsoyuznaya Street 65, Moscow 117997, Russia

Received 29 November 2015; Accepted 3 May 2016

Academic Editor: Francisco Gordillo

Copyright © 2016 V. A. Kubyshkin and S. S. Postnov. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper studies two optimal control problems for linear time-invariant systems of fractional order with lumped parameters whose dynamics is described by equations which contain Riemann-Liouville derivative. The first problem is to find control with minimal norm and the second one is to find control with minimal control time at given restriction for control norm. The problem setting with nonlocal initial conditions is considered which differs from other known settings for integer-order systems and fractional-order systems described in terms of equations with Caputo derivative. Admissible controls are allowed to belong to the class of functions which are -integrable on half segment. The basic investigation approach is the moment method. The correctness and solvability of moment problem are validated for considered problem setting for the system of arbitrary dimension. It is shown that corresponding conditions are analogous to those derived for systems which are described in terms of equations with Caputo derivative. For several particular cases of one- and two-dimensional systems the posed problems are solved explicitly. The dependencies of basic values from derivative index and control time are analyzed. The comparison is performed of obtained results with known results for analogous integer-order systems and fractional-order systems which are described by equations with Caputo derivative.

1. Introduction

Affairs of dynamics and control for fractional-order systems attract sufficient attention of modern research community. This field develops impetuously and is characterized by both of significant theoretical and very actual applied results. And this field is much wider than integer-order dynamics and contains some open problems concerning the foundations of fractional calculus (e.g., the problem of unified definition of fractional derivative and its interpretation). It is known that Caputo derivative and Riemann-Liouville derivative are the most popular definitions of fractional-order derivative. First of them is recognized by many researches as more “physical”: more realistic and similar to ordinary derivative. But this definition imposes rather essential requirements (differentiability) on function from which the fractional derivative is calculated. The Riemann-Liouville definition is used frequently in theoretical investigations although it has a physical sense but less similar to ordinary derivative. Particularly, this derivative is nonzero for constant function and initial and boundary problems for equations with derivative of this kind require posing nonlocal conditions. On the other hand, for existence of Riemann-Liouville derivative from any function only summability of the function is required.

Today optimal control problems for fractional-order dynamical systems actively develop. Many interesting results are already obtained in this area for different types of fractional derivative [1]. So, up to recent times there are no theorems analogous to Pontryagin maximum principle. In 2014 an analogue of this principle was formulated and proved in [2] for dynamical system described by equations with Riemann-Liouville derivative. Later the formulation and proof of Pontryagin-like maximum principle were proposed in [3] for linear systems, described by equations with Caputo derivative. Another effective and quite universal approach to the search of optimal control is moment method [4]. Based on these methods in [5–8] the approach was developed to investigation of optimal control problems for linear dynamic systems of fractional order with lumped parameters, described by equations with Caputo derivative. In [6, 8] this approach was generalized on systems with distributed parameters, described by diffusion-like equation with Caputo time-derivative.

In this paper the moment method is applied to investigation of optimal control problems for dynamical systems of fractional order with lumped parameters, described by equations with Riemann-Liouville derivative. The distinctive feature of problems considered is the posing of nonlocal initial conditions. The optimal control problem for multidimensional linear time-invariant system is considered in general form. Possibility of this problem reduction to moment problem is demonstrated. Then, correctness and solvability are analyzed for the last problem. Further, one-dimensional system of general view and single and double integrators is investigated in detail. Explicit solutions of optimal control problem are obtained and analyzed, including comparison with analogous integer-order systems and systems described by equations with Caputo derivative.

2. Problem Statement

Let the system state and control be defined on a half segment , by vector-functions and correspondingly.

We will consider dynamical systems, described by the following equation:

where are perturbations (known) and are known coefficients, . The repeated indices suppose summation. The fractional derivative of arbitrary order from function is comprehended in our investigation as left-side Riemann-Liouville derivative [9]:

Also we will compare obtained results with results described in [5, 6, 8] for system (1) with Caputo left-side derivative [9]:

Definition 1. System (1) with fractional derivative operator in sense of Riemann-Liouville (Caputo) will be named as Riemann-Liouville (Caputo) system or RL-system (C-system).

The initial conditions for RL-system are determined in nonlocal form [9]:

where

is a left-side Riemann-Liouville integral of order [9]. The substitution is comprehended in sense of limit of the expression in square brackets at . Such initial conditions differ from ordinary (local) form used in optimal control problems. Inherently, in this case no initial state is determined (defined by phase coordinates values in initial time point) but the value of some integral functional from phase coordinates. One of the first papers in which the optimal control problem with nonlocal initial conditions is considered was [10]. Later such statement was studied by many authors, particularly, in [2]. Initial conditions of type (4) can have a quite definite physical interpretation [11, 12]. For example, in some problems of viscoelasticity the stress and strain appear to be coupled by integrodifferential operator of fractional order. Then the initial condition of type (4) for one of these functions can be rewritten in local form for other functions [11].

The final state of system (1) is determined in ordinary form:

Let the control belong to the space , , with norm [4, 13]:

We will consider also the limit case when and the control norm is determined as follows [4, 13]:

Let us believe that functions ,   possessing all properties required for existence of solutions of equations studied further on are, in particular, summable.

We will study the two following statements of optimal control problem.

Problem 2 (OCP A). Find a control , , such that system (1) transfers in given final state (6) from the state determined by (4) and norm of control will be minimal with the assigned control time .

Problem 3 (OCP B). Find a control , , such that system (1) transfers in given final state (6) from the state determined by (4) and control time will be minimal provided , , where is the assigned constant.

3. The Moment Problem

3.1. Preliminaries

The classical -problem of moments can be formulated in the following way.

Problem 4. Let us have a system of functions , ,   . Let us also have the assigned numbers (called moments), and . We should find function , , that satisfies the following conditions:

In order for the moment problem (9) and (10) to be solvable, the existence of number and numbers that give a solution to the following equivalent conditional minimum problem is necessary and sufficient [4, 14, 15].

Problem 5. Find

with additional condition

In case problem (12) and (13) is solvable the optimal control for OCP A will be given by the following expression [4]:

OCP B may be resolved by the following formula:

where is the minimal nonnegative real root of equation

For solvability of Problem 4 it is necessary and sufficient to satisfy one of two equivalent conditions [4]:    and functions are linearly independent. The single question is the possibility of correct posing of Problem 4, which depends on the existence of norm of functions in space and the existence of at least one nonzero component in number set .

Definition 6. The moment problem (9) and (10) is called correct if the norm of functions is determined in space , , and there exists at least one nonzero component in number set .

It is known that optimal control problems in form of Problems 2 and 3 can be reduced to Problem 4 in case of integer-order systems with lumped parameters [4]. It was shown in [5–8] that the same is valid for Caputo systems. Below we will demonstrate that for Riemann-Liouville systems Problems 2 and 3 also can be reduced to Problem 4.

3.2. The Correctness and Solvability of Moment Problem for Multidimensional Fractional-Order System

The general solution of (1) in case of RL-system at can be represented by the following expression [9, .4.2]:

where is matrix -exponent [9], is a matrix of coefficients , and is the vector of values , defined by initial conditions (4). It is also known that general solution of (1) in case of Caputo system can be written as follows [9, 7.4.2]:

As in case of integer-order systems, solutions (17) and (18) at can be written in the form of expression (9) and, consequently, can be represented in the form of moment problem (in case of Caputo system it is demonstrated in detail in [5–8]). And in both of these solutions the components of matrix -exponent act as functions . Other terms in expressions form the moments. Consequently, the moment problem for Riemann-Liouville system differs from the problem for analogous Caputo system only by expressions (and values) for the moments. So, the theorems concerning correctness and solvability of the moment problem, proved for Caputo systems [5–8], will be valid also for analogous Riemann-Liouville systems. Thus, if there exists at least one nonzero moment in the set ,  , we have the following conditions for multidimensional Riemann-Liouville system:

(1)In case of equal differentiation indices the moment problem of type (9) derived from (17) at will be correct and solvable for all which satisfy the following condition:

(2)In case of system (1) at , , and , , the moment problem of type (9) derived from (17) at will be correct and solvable for every and which satisfy the following condition:

3.3. The Moment Problem for One-Dimensional System of Fractional-Order: Problem Setting and Study

In case of the solution of (1) with initial condition (4) can be written as follows [9, 4.1.1] (the subscripts are omitted):

where is two-parameter Mittag-Leffler function [9, 1.8]. By direct calculation one can obtain that at expression (21) with regard to (6) may be written over as (9) with the following symbols:

Expression (22) is identical to analogous expression for one-dimensional Caputo system (formula (17) at in [5]). So, expressions that result from (12) and (13) will match with analogous expressions obtained for the Caputo system [5, 6, 8] (the difference will appear only after substitution of expression (23) for the moment). Consequently, if we keep in mind the noted difference, we can use the solutions for Problem 5 obtained for Caputo system. As shown in [5, 6, 8], the analytical solution of Problem 4 for one-dimensional system can be obtained for arbitrary (if condition (19) is satisfied).

Let . Using (13) we can reduce Problem 5 to simple integral calculation which can be carried out similarly to case of Caputo system.

Consider the case . Direct calculation by formula (12) with regard to (13) and (22) leads us to the following expression (which matches with analogous result for Caputo system [5, 6, 8]):

where . Using (24) we can obtain from (14) and (15) the solutions of OCP A and OCP B correspondingly:

where may be received from (16) with (24). It is seen that controls (25) have not one switching point. It is similar to the system of order behaviour in accordance with Feldbaum’s theorem about intervals [16].

In case of , , we can obtain the following from (12) and (13) subject to (22):

where

Taking into account (14) one can obtain further for OCP A the following:

For OCP B the solution can be expressed by the following formula:

where can be calculated from (16) using (26).

Let us now consider a single integrator, the special case of one-dimensional system at . Instead of (22) and (23) we will have more simple expressions written in terms of elementary functions:

Taking into account formula (30) one can obtain the following from (12) and (13):

From formula (32) using (14) we can derive the following expression for control:

Analogously, in case of OCP B we can derive the following using (15):

where can be calculated from (16) using (32).

Note that expressions (32), (33), and (34) provide an explicit solution of optimal control problem for Riemann-Liouville single integrator at arbitrary .

3.4. The Moment Problem for Double Integrator of Fractional Order: Problem Setting and Study

Two-dimensional system (1) at , , , , and represents itself as a double integrator. Solution (17) subject to initial conditions (4) for this system can be written in the following form:

It is clear that solutions (35) at can be written in form of expression (9) with the following functions and moments:

The functions defined by (36) are identical to the analogous functions obtained in [5] for Caputo double integrator. On the contrary, the moments defined by (37)-(38) differ, in general, from the moments obtained in [5] for Caputo double integrator. Consequently, as for one-dimensional system, the form of solution for moment problem will match with that for Caputo double integrator. So, we can use the solutions of OCP A and OCP B obtained in [5] taking into account the moments (37).

Firstly, consider the case when . It is shown [5] that in this case the minimization problem (12) and (13) reduces to the following algebraic equation:

There is no explicit solution of (39) that can be found at arbitrary and . In some particular cases the explicit solution exists [5], for example, in case of zero-second moment, . Then the solutions of (12) and (13) lead us to the following expression:

The solution of OCP A will be written as

The OCP B explicit solution can be found at additional assumption :

The minimal control time can be calculated based on (16) and (40):

In case of calculations by (12) and (13) give an expression:

where , , , and .

According to (14) the OCP A solution will be written as

The OCP B solution can be derived from (15) and (16) and will be represented by the following expression:

where can be found as least positive real root of the equation

Note also that for Caputo double integrator the case when corresponds to the case when [5]. From formula (38) one can see that the condition is not so clear.

4. Qualitative Analysis of Results

In this chapter the obtained solutions of OCP are analyzed. It will be shown that at integer (equal to 1) values of differentiation indices these solutions reduce to the known results for corresponding integer-order systems. For double integrator analytical expressions derived for system phase trajectories in different modes.

4.1. Optimal Control Behaviour at Integer-Order Differentiation Indices

For one-dimensional system of first order one can obtain the solution of moment problem and explicit expressions for optimal controls. For example, in case of single integrator at the OCP A solution is expressed by the following formula:

The OCP B solution gives the following [4]:

It is easy to see that solution of OCP A defined by (33) (taking into account (31)) at is identical to expression (48). The same is true for OCP B solution: expressions (34) and (16) subject to (32) at give formulas (49).

For double integrator of first order at arbitrary initial and final conditions the solutions of OCP A and OCP B lead to quadratic equation for [4]. Equation (39) reduces to that equation at . The solution of OCP A for double integrator of first order at for is given by the following formula [4]:

By direct calculation one can obtain that formula (41) (i.e., OCP A solution for Riemann-Liouville double integrator) subject to (37) at and reduces to (50).

Consider OCP B in case of . Then the solution of the problem for double integrator of first order will be given by the following formulas:

Expressions (42) and (43) at transform into formulas (51), which can be proved by corresponding substitution.

Thus, as for Caputo systems [5, 6, 8], all the results obtained for single and double integrators of fractional order reduce to corresponding formulas for integer-order systems when differentiation indices are equated to 1.

4.2. Investigation of Qualitative Dynamics for Double Integrator

For two-dimensional systems the analysis of qualitative dynamics is interesting itself. We will calculate below the boundary trajectories for double integrator and its trajectories corresponding to optimal control mode. Consider .

Definition 7. The boundary trajectories of some system are the phase trajectories corresponding to boundary values of control .

In case of differentiation indices equal to 1 the boundary trajectories of some system represent the boundaries of integral vortex for differential inclusion corresponding to this system [17]. This manifold bounds the phase space region which contains all admissible trajectories of the system. In case of fractional-order systems (Caputo and Riemann-Liouville) this property is also valid caused by comparison theorem [18, 19].

Substituting the boundary values of control to solutions (35) one can obtain the following explicit expressions for boundary trajectories:

It is seen that expressions (52) at arbitrary and do not allow eliminating time and obtaining the explicit dependency as opposed to double integrator of integer order [17, 35] and Caputo double integrator [5]. Nevertheless, in some particular cases the time can be eliminated from expressions (52). For example, in case of one can calculate directly that these formulas reduce to the following equation:

Analogously, in case of it is possible to eliminate time from (52) and to derive the following equation:

Note that expressions (52) at reduce to the quadratic equation for integral vortex boundaries of differential inclusion corresponding to the double integrator of integer order [17, 35].

In order to calculate the phase trajectories of Riemann-Liouville double integrator in optimal control mode we will substitute control (41) into solutions (35). Then the following motion laws can be obtained:

where is the switching point of control (41) and is the Heaviside’s function. Formulas (55) are obtained for control (41) which is valid in case of . Consequently, these expressions allow some simplification.

5. Computational Results and Analysis

In this chapter the computational results are represented and analyzed (including comparison with analogous Caputo systems). Calculations and its visualization were performed in MatLab 7.9. Mittag-Leffler function values were calculated using special procedure [20]. Integrals were calculated by Gauss-Kronrod method [21].

5.1. Single Integrator

We will consider the task of system transfer from any state with in the final state .

Figures 1–3 show the dependencies of control norm from differentiation index calculated by formula (32) at and different values for , , and correspondingly (dash-dot lines). Also the analogous dependencies are shown for Caputo single integrator calculated by formulas from [5, 8] at different for the same values of (solid lines). It is seen that curves for Riemann-Liouville and Caputo single integrator differ from each other qualitatively but converge to the same point at . In addition, the control norm increases with or growing for both of integrator types.

Figure 1: The dependencies for norm of control from at (logarithmic scale used for ordinates).

Figure 2: The dependencies for norm of control from at (logarithmic scale used for ordinates).

Figure 3: The dependencies for norm of control from at (logarithmic scale used for ordinates).

In Figure 4 the dependencies of control norm from control time are shown at different values of differentiation index and (for Riemann-Liouville integrator) and (for Caputo integrator). The solid lines correspond to , the dotted lines correspond to , and the dash-dotted lines correspond to . It is clear that all of these curves decrease monotonically which correspond to the case of Caputo integrator [8] and integrator of integer order [4]. And the norm of control increases with growing of and .

Figure 4: The dependencies for norm of control from . Logarithmic scale used for ordinates.

Consider now the minimal control time dependency from differentiation index . In case of and (16) subject to formula (32) can be solved explicitly and have only one real root:

In Figure 5 the computational results are shown for the obtained formula used at different at and (dash-dotted lines). Also the similar dependencies (calculated by formulas from [8]) are represented for Caputo integrator at and (solid lines). As can be seen from Figure 5 the dependency at is linear for both of integrator types. In case of the dependency for Caputo integrator becomes nonlinear.

Figure 5: The dependencies for minimal control time from .

5.2. One-Dimensional System of General Type

Consider (as above) the task of system transfer from any state with to the state .

Figures 6 and 7 represent the dependencies of control norm from differentiation index at and calculated by formulas (24) and (26) correspondingly (dash-doted lines). The analogical dependencies for Caputo one-dimensional system are also shown (solid lines) calculated by formulas from [5, 8]. The curves were calculated for different values of parameter and the following values of other parameters: , , (RL-system), and (C-system). It is seen from Figures 6 and 7 that difference in derivative type has qualitatively effect on behaviour of dependencies. This behaviour changes at and . In case of the curves for C-system decrease and curves for RL-system increase at (analogically to single integrator). At the curves for both of systems increase. For the dependencies increase for both of system types.

Figure 6: The dependencies for control norm from at and different values of parameter .

Figure 7: The dependencies for control norm from at and different values of parameter .

The dependencies of control norm from control time are similar, in general, for single integrator and are not shown here.

The dependencies of minimal control time from at different values of parameter are shown in Figures 8 and 9 for and correspondingly. The following values of other parameters were chosen: , (RL-system), and (C-system). Dash-dotted lines correspond to RL-system and solid lines correspond to C-system. It is seen that in case of RL-system the curves have the quasi-linear trend and differ qualitatively from the curves for C-system.

Figure 8: The dependencies for minimal control time from at and different values of parameter .

Figure 9: The dependencies for minimal control time from at and different values of parameter .

5.3. Double Integrator

As in case of one-dimensional systems, we will consider the task of system transfer from any state with nonzero values of phase coordinates into the origin at phase plane: we will suppose that and .

The dependencies of control norm from one of differentiation indices for Riemann-Liouville double integrator at several fixed values of other indexes are shown in Figures 10 and 11 (dash-dotted lines). The curves were calculated by expression (40) for , , and . Also the similar dependencies for Caputo double integrator are shown at these figures calculated for the same values of parameters and and for (solid lines). The dependencies for are similar and not demonstrated here. It is seen from Figure 10 that curves corresponding to different systems (RL- and C-systems) at identical values of converge into the same point at . In Figure 11 the distance between different curves for C- and RL-systems decreases when increases.

Figure 10: The dependencies for control norm from at and several values of .

Figure 11: The dependencies for control norm from at and several values of .

The dependencies of control norm from control time for double integrator are analogous, generally, to that in case of one-dimensional systems and are not demonstrated here.

Let us consider the dependency of minimal control time from differentiation indices and . In case of this dependency is given by formula (43). In case of with additional assumption the explicit formula also can be found for such dependency (as solution of (16) subject to (40)):

In Figures 12–15 the dependencies are shown for minimal control time from one of differentiation indices at fixed value of other indices. The curves are calculated at for (Figures 12 and 13) and (Figures 14 and 15). Dash-dotted lines correspond to Riemann-Liouville double integrator and are calculated for . Solid lines correspond to Caputo double integrator and are calculated for . As seen from these figures, the dependencies for RL- and C-systems differ qualitatively from each other.

Figure 12: The dependencies for minimal control time from at different values of for . Logarithmic scale used for ordinates.

Figure 13: The dependencies for minimal control time from at different values of for . Logarithmic scale used for ordinates.

Figure 14: The dependencies for minimal control time from at different values of for . Logarithmic scale used for ordinates.

Figure 15: The dependencies for minimal control time from at different values of for . Logarithmic scale used for ordinates.

Let us now analyze the qualitative dynamics of Riemann-Liouville double integrator, its boundary trajectories, and phase trajectories in optimal control mode. Note that last type of trajectories in case of Riemann-Liouville double integrator will not have an obvious initial point at phase plane since the initial condition (4) defines not the system state but the value of some integral functional at initial time. The system state in this case can be undefined as follows from solutions (35) and expressions (55).

The boundary trajectories for studied system are represented in Figure 16 for Riemann-Liouville double integrator (dash-dotted lines) and for analogous Caputo system (solid lines). For the last system calculations were performed using formulas from [5]. It is clear that trajectories for both of systems are very close to each other. In Figure 17 the fragment of trajectories is shown at enlarged scale.

Figure 16: The boundary trajectories for double integrator of fractional order at and .

Figure 17: The boundary trajectories for double integrator of fractional order at and .

The phase trajectories of Riemann-Liouville double integrator in optimal control mode were calculated using expressions (55). In Figure 18 the results are shown for different values of differentiation indices, and . In Figure 19 the phase trajectories of Caputo double integrator in optimal control mode are represented which were calculated by formulas from [5] at and . It is clear from Figures 18 and 19 that both of systems demonstrate the overcontrolling effect. In case of this effect does not appear and phase trajectories coincide for both of systems studied (Figure 20).

Figure 18: The phase trajectories in optimal control mode for Riemann-Liouville double integrator.

Figure 19: The phase trajectories in optimal control mode for Caputo double integrator.

Figure 20: The phase trajectories in optimal control mode for Riemann-Liouville and Caputo double integrator in case of .

6. Summary

In this paper the optimal control problem is investigated for linear dynamical systems of fractional order described by equations with Riemann-Liouville derivative. The problem reduced to the moment problem. For the last problem the correctness and solvability conditions derived. Some examples of one- and two-dimensional systems were considered and explicit solutions of corresponding optimal control problems were found. The properties of these solutions were analyzed and compared with analogous integer-order systems and fractional-order systems, described by equations with Caputo derivative.

Some computational results are demonstrated that show the dependency of control norm from differentiation indices ant control time and the dependency of minimal control time from differentiation indices. Qualitative dynamics of double integrator are analyzed and show that boundary trajectories of Riemann-Liouville system of this type are very close to the same trajectories for corresponding Caputo system. It is shown also that the overcontrolling effect appears in optimal control mode of system studied.

Obtained results can be useful in investigation of optimal control problems for fractional-order systems and in design of particular control systems.

Competing Interests

The authors declare that they have no competing interests.

References

A. G. Butkovskii, S. S. Postnov, and E. A. Postnova, “Fractional integro-differential calculus and its control-theoretical applications. II. Fractional dynamic systems: modeling and hardware implementation,” Automation and Remote Control, vol. 74, no. 5, pp. 725–749, 2013.

R. Kamocki, “Pontryagin maximum principle for fractional ordinary optimal control problems,” Mathematical Methods in the Applied Sciences, vol. 37, no. 11, pp. 1668–1686, 2014.

R. Kamocki and M. Majewski, “Fractional linear control systems with Caputo derivative and their optimization,” Optimal Control Applications and Methods, vol. 36, no. 6, pp. 953–967, 2015.

A. G. Butkovskiy, Distributed Control Systems, American Elsevier, New York, NY, USA, 1969.

V. A. Kubyshkin and S. S. Postnov, “Optimal control problem for a linear stationary fractional order system in the form of a problem of moments: problem setting and a study,” Automation and Remote Control, vol. 75, no. 5, pp. 805–817, 2014.

V. A. Kubyshkin and S. S. Postnov, “Optimal control problem for linear fractional-order systems,” in Proceedings of the International Conference on Fractional Differentiation and Its Applications (ICFDA '14), p. 6, IEEE, Catania, Italy, June 2014, Paper ID 3189701.

V. A. Kubyshkin and S. S. Postnov, “Analysis of two optimal control problems for a fractional-order pendulum by the method of moments,” Automation and Remote Control, vol. 76, no. 7, pp. 1302–1314, 2015.

V. A. Kubyshkin and S. S. Postnov, “The optimal control problem for linear systems of non-integer order with lumped and distributed parameters, Discontinuity,” Nonlinearity and Complexity, vol. 4, no. 4, pp. 429–443, 2015.

A. A. Kilbas, H. M. Srivastava, and J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204, Elsevier, Amsterdam, The Netherlands, 2006.

R. Almeida and D. F. M. Torres, “Calculus of variations with fractional derivatives and fractional integrals,” Applied Mathematics Letters, vol. 22, no. 12, pp. 1816–1820, 2009.

N. Heymans and I. Podlubny, “Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives,” Rheologica Acta, vol. 45, no. 5, pp. 765–771, 2006.

V. E. Tarasov, Fractional Dynamics, Springer, Berlin, Germany, 2010.

A. N. Kolmogorov and S. V. Fomin, Elements of the Theory of Functions and Functional Analysis, Dover Publications, Mineola, NY, USA, 1999.

N. I. Akhiezer, The Classical Moment Problem and Some Related Questions in Analysis, Hafner, New York, NY, USA, 1965.

M. G. Krein and A. A. Nudelman, “The Markov moment problem and extremal problems,” in Ideas and problems of P. L. Chebyshev and A. A. Markov and their further development, vol. 50 of Translations of Mathematical Monographs, American Mathematical Society, Providence, RI, USA, 1977.

A. A. Feldbaum, Optimal Control Systems, Academic Press, New York, NY, USA, 1965.

A. G. Butkovskiy, Phase Portraits of Control Dynamical Systems, Kluwer Academic, Dordrecht, The Netherlands, 1991.

C. Li, D. Qian, and Y. Chen, “On Riemann-Liouville and caputo derivatives,” Discrete Dynamics in Nature and Society, vol. 2011, Article ID 562494, 15 pages, 2011.

Y. Li, Y. Q. Chen, and I. Podlubny, “Mittag-Leffler stability of fractional order nonlinear dynamic systems,” Automatica, vol. 45, no. 8, pp. 1965–1969, 2009.

http://www.mathworks.com/matlabcentral/fileexchange/8738-mittag-leffler-function

L. F. Shampine, “Vectorized adaptive quadrature in Matlab,” Journal of Computational and Applied Mathematics, vol. 211, no. 2, pp. 131–140, 2008.