Fractional Order Calculus : Basic Concepts and Engineering Applications

The fractional order calculus FOC is as old as the integer one although up to recently its application was exclusively in mathematics. Many real systems are better described with FOC differential equations as it is a well-suited tool to analyze problems of fractal dimension, with longterm “memory” and chaotic behavior. Those characteristics have attracted the engineers’ interest in the latter years, and now it is a tool used in almost every area of science. This paper introduces the fundamentals of the FOC and some applications in systems’ identification, control, mechatronics, and robotics, where it is a promissory research field.


Introduction
The fractional order calculus FOC was unexplored in engineering, because of its inherent complexity, the apparent self-sufficiency of the integer order calculus IOC , and the fact that it does not have a fully acceptable geometrical or physical interpretation 1, 2 .Notwithstanding it represents more accurately some natural behavior related to different areas of engineering, and now it is used as a promissory tool in bioengineering 3, 4 , viscoelasticity 5, 6 , electronics 7, 8 , robotics 9-11 , control theory 12, 13 , and signal processing 14, 15 between others.
In the latter years FOC attracted engineers' attention, because it can describe the behavior of real dynamical systems in compact expressions, taking into account nonlocal characteristics like "infinite memory" [16][17][18] .Some instance are thermal diffusion phenomenon 19 , botanical electrical impedances 20 , model of love between humans 21 , the relaxation of water on a porous dyke whose damping ratio is independent of the mass of moving water 22 , and so forth.On the other hand, direction the behavior of a process with fractional order controllers would be an advantage, because the responses are not restricted to a sum of exponential functions, therefore a wide range of responses neglected by integer order calculus would be approached 23 .Bearing these ideas in mind, this paper is organized as follows.Section 2 presents the fundamentals and analytical definitions.Section 3 introduces several approaches to the solution of the 300-years-old problem of the geometrical interpretation of the FOC.Then in Section 4 some applications in systems' identification, control, and robotics are presented.Finally Section 5 introduces the main conclusions and future applications of FOC.

Fractional Order Calculus (FOC)
The intuitive idea of FOC is as old as IOC, it can be observed from a letter written by Leibniz to L'Hopital in 1695 24 .It is a generalization of the IOC to a real or complex order 25 .Formally the real order generalization is introduced as follows: Its applications in engineering were delayed because FOC has multiple definitions 18, 26 , there is not a simple geometrical interpretation and the IOC seems, at first sight, to be enough to solve engineering problems.However, many natural phenomena may be better described by a FOC formulation, because it takes into account the past behavior and it is compact when expressing high-order dynamics 27, 28 .Some common definitions of FOC are listed as follows 18, 26 : i Riemann-Liouville: Integral: Derivative: ii Gr ünwald-Letnikov: Integral: Derivative: iii Caputo: iv Cauchy: where the function Γ α is the generalization of factorial function 29 and it is defined as: or without a restriction for x We can choose one definition or another, depending on the application and the preference of the designer.In 26 the authors compare these definitions in applications of control and signals processing, finding that the Cauchy definition preserves some important frequency properties, that also exist in IOC simplifying the data's interpretation.Some other tools of interest for engineers are the classical transforms of Laplace and Fourier, that are valid and used in order to simplify operations like convolution and can be applied to solve FOC differential equations.In FOC the Laplace transform is defined as 30

Mathematical Problems in Engineering
As shown, this transform takes into account all initial conditions from the first to the nth − 1 derivative.In practice, the Fourier transform can be achieved by replacing s by jw in 2.10 .
In addition to the problem for which definition must be chosen based on its properties or implementation complexity, the engineers may know the implications of using a mathematical tool.An easy way to understand it, is by plotting it in a figure and seeing what is happening when it is applied.Pitifully for FOC it is a lack, but some approaches were proposed in the last decade, as will be presented in Section 3.

Geometrical Interpretation
In the case of integral order calculus, there is a well-accepted geometrical explanation which clearly relates some physical quantities, for example, instant rate of change of a function completely explains the relationship between concepts like position and speed of an object.Unfortunately, until the last decade there was no geometrical interpretation of the fractional order derivatives.One of them was proposed in 1 explaining FOC from a probabilistic point of view, using the Gr ünwald-Letnikov definition 2.4 and 2.5 .If α is a value between 0 and 1, and γ is defined as Then for m 0 we obtain γ 1, that is, the value of the function at evaluation time present appears with probability of 1.
If m > 0, For values of m / 0 the γ value vanishes when the analysis point is far from the evaluation one.Therefore, the author suggests that the expression − ∞ m 1 γ α, m x t − mh is the expected value of a random variable X, where P X x mh γ α, m , m 1, 2, . . ., 0 < α < 1.

3.3
Therefore the values near to the evaluation time present have more influence over the result than those that are far from it.This interpretation is shown in Figure 1.
A geometric interpretation based on Riemann-Louville-definition 2.2 and 2.3 was presented in 31 .This definition can be written as  With this information a tridimensional graph is drawn with axes g τ , f τ , and τ as shown in Figure 2. The projection of the area below f t , over the plane τ, f τ , is t 0 f τ dτ, the same as the integer integral definition.The projection of the area below the curve, over the plane f τ , g τ , is t 0 f τ dg τ .Note that it is the same definition as 3.4 , that is, the integral of the function with a non homogeneous time scale that depends on the parameter α.
Another geometrical interpretation, this time based fractal dimension was proposed in 32 .Here the author argues that the Riemman-Lioville 2.2 is the convolution of the function f t with kernel: For α 0, the function h ∞ is undetermined.By increasing α fractional integral case , the kernel h ∞ takes into account the effect of the past values, weighting them Figure 3 .If α 1 integer integral case , then J α c has perfect memory and all the past is equally weighted.In the derivative case −1 < α < 0 , the interpretation cannot be explicit obtained from 2.2 .Using the Leibniz rule on 2.3 we obtain Note that the kernel of 3.6 naturally appears when α is replaced by 1 − β and it regulates the effect of the past in a β-proportional rate.The derivative value is the sum of the effect of the initial condition and the value of the integer derivative, both regulated by the kernel h ∞ .The kernel behavior is similar to the Cantor set 33 , that is, an iterative function that removes the middle section of a line Figure 4 .The fractal dimension of the Cantor set is defined by In this case the Dth dimension represents the density of the remaining bars and it is analogous to the kernel h ∞ , but in discrete time.

Engineering Applications
Fractional order can represent systems with high-order dynamics and complex nonlinear phenomena using few coefficients 6, 34, 35 , since the arbitrary order of the derivatives gives an additional degree of freedom to fit an specific behavior.Another important characteristic is that the fractional order derivatives depend not only on local conditions of the evaluated time, but also on all the history of the function.This fact is often useful when the system has long-term"memory" and any evaluation point depends on the past values of the function.However, it is also a problem when fractional derivative functions are implemented in logical circuits, because they require a huge quantity of physical memory.The strategies to simulate fractional order systems are classified in three groups 36, 37 .
i Computational methods based on the analytic equation.These methods present multiple parameters and are complicated to analyze, as it is necessary to evaluate every single point in the function and its history; moreover, the explicit equation if often difficult to obtain.
ii Approximation through a rational system in discrete time.The analytical system is replaced by its discrete equivalent in frequency space.Those methods result in irrational coefficients, that are approximated again by truncating the polynomial series, which is equivalent to truncate the model in the time domain; therefore, it requires as minimum the same number of coefficients as samples, losing the characteristic of "infinite memory".On the other hand, if the series has a lot of coefficients, it limits simulation in real time, as it requires more processing cycles.
iii Approximation of the fractional system using rational function in continuous time.This is approximated by rational continuous approach, but the series must be truncated; therefore, it must be limited to a specific frequency range of operation.

Electronics Applications
Another way to obtain the response of a fractional order system is by using analogical circuits with fractional order behavior as shown in Figure 5 or systems with fractal configuration as shown in Figure 6 a .Here three methods are introduced.
i Component by component implementation 29, 32 .The approximation of the transfer function is done by the recursive circuit shown in Figure 5.The gain between V o and V i in Laplace transform is the continuous fraction approximation to the original system 38 , that is where This circuit has two principal disadvantages: 1 it has a limited frequency band of work, and 2 this is an approximation, therefore it requires a lot of low tolerance components, depending on the accuracy required by the designer.ii Field Programmable Analog Array (FPAA) 39 .The designer implements the circuit component by component into a FPAA.It allows changing of the dynamical behavior of the fractional order system with a few simple modifications and each element has custom tolerance.
iii Fractional order impedance component.It is a capacitor with fractional order behavior introduced in 40 .In general it consists in a capacitor of parallel plates, where one of them presents fractal dimension Figure 6 a .Each branch could be modeled as a low pass resistor/capacitor RC circuit filter and it is linked to the principal branch, as shown in Figures 6 b and 6 c .
Anyone of these approaches could be used in engineering applications, in this paper we introduce its use in systems' identification, control theory and robotics.

Fractional Order Identification of Dynamical Systems
Fractional order dynamical systems can be modeled using the Laplace transform-like transfer functions 41 as Some high-order systems would be approximated with a compact fractional order expression, it is useful in cases where an approach between holistic and detailed description of the process is required.As an instance the model of the 5th order 7 G s s 4 36s 3 126s 2 84s 9 9s 4 84s 3 126s 2 36s 1 .

4.3
This 8-parameter system would be well approximated by G s ≈ 1/s 0.5 , a compact fractional order system with just a parameter, valid in the frequency range from 100 to 10000 Hz, as shown in Figure 7.
Many real systems are better identified as fractional order equations 16, 42 than integer ones.In fact, some responses cannot be approximated just as a linear combination of exponential functions 43 , and the arbitrary order is an additional degree of freedom that a Fractal tree yields a better approximation to the real system while describing it in a compact way 44 .In 45 it was used this fact to identify a fractal system, typically modeled in frequency as:  Adjusting the model is accomplished by finding the parameters {K, a, α} that minimize the mean error with the real data.
Another instance of the fractional order formulation is presented in 46 , the authors approximated a complex system, a flexible structure with five vibration modes, modeling it with few parameters, being still valid for a wide range of frequencies.They propose the transfer function: with α 1, α 2, and α 0.5.A real value of α models the damper behavior without increasing the order of the system, and maintaining a compact expression too, valid for the frequency range 0.1 Hz-200 Hz .Another example of identification of a biological system was presented in 20 , the authors note that their frequency response does not decay/increase in multiples of 20 dB/dec in the Bode's plot.It may occurs because fruit and vegetable's electrical properties depends on several parameters as type of fruit/vegetable, size, temperature, and pressure between others.As the author demonstrate by the experiment shown in Figure 8, by applying a sine voltage and analyzing the current over the object.They found that the response in frequency has a fractional order behavior with a constant slope, depending on the fruit/vegetable.
A nonparametric method introduced in 47 uses a continuous neural network CNN in order to identify nonlinear systems.This type of networks uses integral blocks instead of time delays.This fact makes the model continuous and its behavior is not a "black box" anymore.From this kind of network is possible to separate the static nonlinear system neural network from the dynamical one integral blocks .If the integral blocks are fractional order blocks, then the CNN captures the fractional behavior too.In order to train the network, the authors used the square mean error between the system output y r and the neural network output y c Figure 9 .
Just as an example, we propose an experiment with synthetic data, simulating the vibration present in a gearbox.These kinds of systems are highly complex as several frequencies and their harmonics are exited by the rotation of the axes, unbalanced pieces, meshing between gears, bearing balls interaction, backslash between pieces among others.
When the system has a failure, harmonics and side-bands are added to the frequency Magnitude of the Bode's plot of two complex systems, one represents the vibration signal of a rotational system without a failure a and the other is b a system with a teeth broken on the transmission box.Note that when approximating by a fractional order equations the order changes from a system to another.
spectrum and the dynamical model of the system may change.If these models were known a predictive maintenance strategy would be proposed based on comparison between them.
Unfortunately as there are many components interacting and some have nonlinear behavior, a dynamical model of integer order is frequently difficult to obtain and involve several parameters that are hardly comparable.Notwithstanding, as shown in Figure 10 the signal on the Bode's plot does not decay by 20 dB/dec, hence the systems would be approximate by a fractional order equations.When a failure is introduced, the model of the system change.In this case the failure was identified with just one parameter, the order of the equation.

Fractional Order Control
Dynamic systems are typically fractional order, but often just the controller is designed as that, as the plant is modeled with integer order differintegral operators.A robust fractional order controller requires less coefficients than the integer one 48 .Grouping by type of plant and controller, the systems are classified in four sets 49 , as shown in Table 1.In 49 it is proved that fractional order controllers are more robust than integer order.The authors proposed two dynamic systems with three coefficients, 1 an integer system of second order and 2 a system of fractional order with three coefficient.They optimized those controllers and found that fractional algorithms were more stable taking into account stationary error and the overshoot percentage.
The typical fractional controller in literature are 27 as follows.
i Tilted Proportional and Integral (TID).It is a controller similar to the PID of integer order in its architecture, but replacing the proportional component by a function like s −α , with α ∈ R. It gives an additional degree of freedom to the system and allows a better behavior than that of the integer order controller.A block diagram of TID controllers is shown in Figure 11.
ii Acronym in French of Crontrôle Robuste d'Ordre Non Entier (CRONE).These type of controllers are based on "fractal robustness" a damping behavior that is independent of the mass observed in water dykes 50 in which the conjugated roots of the characteristic equation of the system can move over a fixed angle in the complex plane.When analyzed in feedback, the system has a constant phase 4.8 .This result is identical to the phase of the proposed system in open loop for high frequencies.Therefore, it implies that the controller is robust in this characteristic, which is directly related with the overshoot and the damper factor.
The function approximation to the dyke behavior was Therefore, in feedback with a negative gain G s 1 e α ln τω G jω e −α ln τω , 4.7 iii Algorithm PI α D μ .This is the generalization of the integer PID.The general structure of this kind of controllers is There is not a rigorous formula to design this type of controller, some techniques to adjust it are artificial intelligence, as swarm intelligence 51 , genetic algorithms 52 or other where the parameter space has five variables K p , K i , K d , α, μ .A block diagram of PI α D μ is shown in Figure 12.
iv Fractional lead-lag controller.It is the generalization of the lead-lag controller of integer order.It can be written as where 0 < ω b < ω h , C 0 > 0 and r ∈ 0, 1 In 53 the author proposes a general optimization architecture, based on Caputo formula, where the system equation is optimized in Lagrange terms as follow: where λ is the Lagrange multiplier and the initial conditions are known.

Applications in Robotics
In industrial environments the robots have to execute their task quickly and precisely, minimizing production time.It requires flexible robots working in large workspaces; therefore, they are influenced by nonlinear and fractional order dynamic effects 10 .For instance in 54, 55 the authors analyze the behavior of two links in a redundant robot a robot that has more degree of freedom than required to carry out its task following a circular trajectory in the Cartesian space.By calculating the inverse kinematics, the pseudoinverse matrix does not converge into an optimal solution either for repeatability or manipulability.
In fact, the configuration of those links has a chaotic behavior that can be approximated by fractional order equations, since it is a phenomena that depends on the long-term history, as introduced in 56 .Another case-fractional order behavior in robotics was presented in 10 , where a robot of three degrees of freedom was analyzed by following a circular trajectory, controlled with a predictive control algorithm on each joint.Despite it has an integer order model, the current of all motors at the joints presents clearly a fractional order behavior.
In 57 , the authors analyze the effect of a hybrid force and position fractional controller applied to two robotic arms holding the same object, as shown in Figure 13.The load of the object varied and some disturbances are applied as reference of force and position.A PI α D μ controller was tuned by trial and error.The resulting controller was demonstrated to be robust to variable loads and small disturbances at the reference.
Another interesting problem in robotics which can be treated, with FOC is the control of flexible robots, as this kind of light robots use low power actuators, without self-destruction effects when high impacts occurs.Nevertheless significant vibrations over flexible links make a position control difficult to design, because it reveals a complex behavior difficult to approximate by linear differential equations 58 .However in 59 , the authors propose a PD α for a flexible robot of one degree of freedom with variable load, resulting in a system with static phase and constant overshoot, independent of the applied load.
Another case was analyzed in 60 , simulating a robot with two degrees of freedom, and some different physical characteristics, as an ideal robot, a robot with backslash and a robot with flexible joints.In each one of these configurations they applied PID and PI α D μ controllers and their behavior was compared.These controllers were tuned by trial and error in order to achieve a behavior close to the ideal and tested 10000 trajectories with different type of accelerations 61 .Over the ideal robot, the PID controller had a smaller response time and smaller overshoot peak than the fractional order PID.When any kind of nonlinearity is added to the model, the fractional controller has a smaller overshoot and a smaller stationary error, demonstrating that these type of controllers are more robust than classical PID to nonlinear effects.
A small overshoot in fractional order controllers is an important characteristic when accuracy and speed are desired in small spaces.In 62 the authors used CRONE controllers in order to reduce the overshoot on small displacement over a XY robot.The workspace is of 1 mm 2 and the overshoot obtained was lower than 1%.
An application in a robot with legs was presented in 63, 64 , designing a set of PD α algorithms in order to control position and force, applied to an hexapod robot with 12 degrees of freedom.The authors defined two performance metrics, one for quantity of energy and the other for position error.The controllers with α 0.5 had the best performance in this robot.

Conclusions
In this paper some basic concepts of FOC and some applications in engineering were presented.However, its inherent complexity, the lack of a clear geometrical interpretation and the apparent sufficiency of the integer calculus have delayed its use outside the area of mathematics.Nowadays, some applications have begun to appear but they are still at the initial stage of development.In the near future, with a deep understanding of FOC's implications, its use in systems' identification will increase, as it captures very complex behavior neglected by IOC, and in control of systems this tool open a wide range of desired behavior, where the integer one is just a special case.

Figure 2 :
Figure 2: Podlubny fractional order derivative interpretation.f t t and α 0.3 evaluated in the interval 0, 3 .The fractional derivative is the projection of the area below the functions over a nonlinear time scale g τ , with a deformation parameter α order of the derivative .

6 MathematicalFigure 3 :
Figure 3: Function kernel h ∞ evaluated for different α values.When convoluted with a function, f t drives the weigh of the "memory" in the output.

Figure 4 :
Figure 4: Cantor set with α 1/3.Note that the density of the remaining segments is similar to the decay of the kernel function h ∞ .

Figure 6 :
Figure 6: Fractor is a parallel capacitance with fractional order behavior.It uses fractal geometry when fabricated. a Introducing a type of fractal tree.b Presenting the link diagram, and c The circuit equivalence.

Figure 7 :
Figure 7:Comparison between a high-order integer system and its approximation by a fractional one.

Figure 8 :Figure 9 :
Figure 8: Circuit used for identification of fractance in fruits and vegetables.

Figure 10 :
Figure10: Magnitude of the Bode's plot of two complex systems, one represents the vibration signal of a rotational system without a failure a and the other is b a system with a teeth broken on the transmission box.Note that when approximating by a fractional order equations the order changes from a system to another.

L o a d Robot 1 Robot 2 Figure 13 :
Figure 13: A cooperative cell of robots achieving a desired task.

Table 1 :
Classification of dynamical system grouping by order of the plant and the controller.