Microsoft Word - 4_Vol17_No1_2016_(CompiledVersion_4).docx IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 115 CHAOS CONTROL AND SYNCHRONIZATION USING SYNERGETIC CONTROLLER WITH FRACTIONAL AND LINEAR EXTENDED MANIFOLD MORTEZA POURMEHDI, ABOLFAZL RANJBAR NOEI*, AND JALIL SADATI Department of Electrical & Computer Engineering, Babol University of Technology, Babol, Iran. morteza.pourmehdi1366@gmail.com; *a.ranjbar@nit.ac.ir; and j.sadati@nit.ac.ir (Received: 02 Jun. 2015; Accepted: 20 Sept. 2015; Published on-line: 30 Apr. 2016) ABSTRACT: In this manuscript, for the first time, a fractional-order manifold in a synergetic approach using a fractional order controller is introduced. Furtheremore, in the synergetic theory a macro variable is expended into a linear combination of state variables. An aim is to increase the convergence rate as well as time response of the whole closed loop system. Quality of the proposed controller is investigated to control and synchronize a nonlinear chaotic Coullet system in comparison with an integer order manifold synergetic controller. The stability of the proposed controller is proven using the Lyapunov method. In this regard stabilizing control effort is yielded. Simulation result confirm convergence of states towards zero. This is achieved through a control effort with fewer oscillations and lower amplitude of signls which confirm feasibility of the control effort in practice. ABSTRAK: Manuskrip ini memperkenalkan, buat kali pertama, manifold pecahan pesanan dalam pendekatan sinergi yang menggunakan pengawal pecahan pesanan. Dalam teori sinergi pembolehubah makro diperluaskan menjadi satu kombinasi linear pembolehubah keadaan. Tujuannya adalah untuk meningkatkan kadar penumpuan serta masa tindak balas keseluruhan sistem lingkaran tertutup. Kualiti pengawal yang dicadangkan dikaji dari segi mengawal dan menyelaras sistem Coullet tak linear berbanding dengan pengawal sinergi pesanan integer manifold. Kestabilan pengawal yang dicadangkan terbukti dengan penggunaan kaedah Lyapunov. Keadaan ini menghasilkan kawalan yang stabil. Simulasi menghasilkan penumpuan keadaan ke arah sifar. Ini dicapai melalui usaha kawalan yang mempunyai kurang pengayunan dan amplitud isyarat lebih rendah dan ini mengesahkan bahawa usaha kawalan dalam amalan boleh dilaksanakan. KEYWORDS: Synergetic control theory; Fractional order system; Synchronization; Nonlinear chaotic Coullet system; Chaos control 1. INTRODUCTION Synergetic control theory is primarily reported by Russian researchers [1]. In essence, a synergetic controller is of an analytical topic with high flexibility in defining dynamical manifold. Fractional order controllers are shown providing quality and robust performance in the presence of uncertainties and disturbances for nonlinear systems [2]. During the design procedure, undesired dynamics can be eliminated by introducing dynamical constraints, defined in a manifold. In this regard specifications of desired performance can be applied to the system according to this manifold. Elimination of undesired dynamics and reduction of the system’s order are some issues of desired control specifications. The aim of the synergetic controller is to force states of the system to reach and remain on a manifold IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 116 through macro variables. These variables are selected according to the control goals. When system reaches to the desired manifold, its behavior is usually controlled by a first order differential equation. Innovatively a combination of fractional controller and the synergetic controller is used to improve the tuning process and the convergence of states. Although synergetic controller is recently proposed, several applications are successfully reported in different fields of engineering. This controller is successfully applied in nonlinear power electronic and industrial processes. Nonlinear power system stabilizers [3],power converters for pulse current charging [4], DC-DC boost converters [5] and control of chaotic oscillation in power systems [6] are such of these applications. Fractional calculus by over 300 years’ history is a generalization of integer order derivative and integral into a real order one. In recent decades, applications of fractional calculus in modeling and control are widely reported. Several works such as [7-16] deals with fractional calculus. Fractional order dynamics are involved in deferent systems such as viscoelasticity materials, electrochemical processes, electrical machines and etc. It is shown in many fields of science and engineering that the fractional operators can provide more accurate modeling[17,18]. During modeling of physical systems which exhibit hereditary, diffusion and viscoelasticity properties [19-22], heat transfer [23], [24], [25], [12] and fractional PID [26] controllers can be mentioned as such examples of applying fractional order controllers in different field. Recently [27] incorporates a combination of synergetic controller with fuzzy theory. In this manuscript, fractional derivative in the Caputo sense [15] will be utilized. The current research proposes an idea to combine both synergetic controller and fractional operators to control a fractional order chaotic system. This proposed approach shows improving time response and the convergence rate. The rest of the paper is organized as follows: In section 2, basics of synergetic theory are introduced. Section 3 briefly describes the Caputo fractional calculus. Synergetic controller with fractional order and together with combined linear extension of the manifold is proposed for the first time in section 4. The stability of the control system is analyzed in this section. Section 5 investigates with quality of the proposed method through a chaos control of a fractional order system. Synchronization is made possible in section 6 using the proposed fractional-order synergetic structure. Finally, section 7 closes the work by a conclusion. 2. BASICS OF SYNERGETIC THEORY Synergetic method generates a control effort through defining a macro variable introducing required dynamic as a manifold. Assume that the system is represented in the following nonlinear state space equations: (1) ̇ = ( , , ) Where denotes vector of the state variable, u refers to the vector of the control input, defines the time variable and finally ( , , ) indicates a nonlinear relation between states and input. Synergetic design uses macro variable ψ( ) in either linear or nonlinear function of the state variable such as: (2) = ( ) IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 117 Appropriate control law drives states towards and stays on manifold of ( ) = 0. The designer may choose macro variable according to the control specifications and requirements [28]. Macro variables are usually considered as a liner combination of the state variables. Number of macro variables may be defined as much as the control input channels (multi input). A dynamical equation of manifold may be defined as follows [2]: (3) ̇ + = 0, > 0 where denotes the design parameter. This eventually adjusts rate of the convergence of the states in the manifold (2). Equation (3) immediately yields: (4) ( ) = , ≥ 0 The convergence of ( ) towards the manifold ( ) = 0 is guaranteed with any bounded initial condition where here the rate of convergence will be tuned by appropriate choice of the design parameter . Roll of parameter can be found using the following chain rule of Eqn.(4): (5) ̇ = = ̇ Substituting equation Eqn.(5) into Eqn.(3), involving the system Eqn.(1) achieves: (6) T (x, u, t) + = 0 Finally, an appropriate control law will be achieved using Eqn.(6) considering other specifications. Each manifold imposes a new constraint to the system which decreases the order of the system. 3. THE CAPUTO FRACTIONAL DERIVATIVE In fractional calculus several methods of Grünwald–Letnikov, Riemann–Liouville and Caputo are commonly used. In this manuscript, the Caputo fractional derivative [15] is used. Benefits of the Caputo technique are studied in [29-32]. This is because initial conditions in the Caputo are of an integer order which is physically realizable. In the Riemann–Liouville definition the initial condition is of a fractional order that is less realizable in practice as well as increases the complexity [15]. The Caputo fractional derivative is defined as follows: (7) ( ) = 1 Γ( − ) ( )( ) ( − ) where Γ(. ) is the Euler Gamma function. For − 1 < < , initial condition of the fractional order differential equation is the same as the integer order one [15]. A unified formula for fractional order integral of ( ) i.e. ( ), with ∈ (0,1) in sense of the Caputo in the time interval [0,t] is defined as follows [15]: (8) ( ) = 1 Γ( ) ( − ) ( ) IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 118 Likewise ( ), − 1 < < is a Caputo fractional order derivative of as defines as follows [15]: (9) ( ( )) = ( ) − ( )( ) ! ( − ) For = 0, = 0 and 0 ≤ < 1 [15] It is shown that he following relation is yielded: (10) ( ( )) = ( ) − (0) Furthermore, the following relations are also held for the Caputo method [15]: (11) ( ( )) = ( ), ( = 0,1,2, … ; − 1 < < ) Meanwhile: (12) ( ( )) = ( ( )) = ( ) ( )(0) = 0 , ∈ ℕ , = , + 1, … , ( = 0,1,2, … ; − 1 < < ) In the following, the above fractional calculus is combined with the synergetic control theory to obtain appropriate control law. 4. METHODOLOGY: ANALYTICAL INVESTIGATION OF THE STABILITY OF FRACTIONAL-ORDER MANIFOLD SYNERGETIC CONTROLLER In this section, the synergetic controller in Eqn.(2) to Eqn.(6) gains a fractional manifold. The combination gives advantages of both the synergetic theory and the fractional calculus. The latter provides more flexibility and degree of freedom to design the fractional controller. When a higher integer order is needed a small size fractional controller provides better results [17, 33-35].The idea is also developed to consider a linear combination of state variables. In this regard, a fractional order form of the synergetic manifold of Eqn.(3) is proposed as follows: (13) ( ) + ( ) = 0, > 0 where is a real number in the interval ∈ (0, 1). In the following, performance of the proposed controller is compared with conventional synergetic controller of integer order manifold when fractional order chaotic Coullet system [36] is controlled. Chaotic systems are found very sensitive to initial conditions. This means two identical systems but with a minor deviation in their initial conditions may produce a completely different result. Hence, it is difficult predict the behavior of these systems. The term Chaos arises from a deterministic dynamic with nearly stochastic (almost unpredictable) behavior. The chaos effect and control is a more challenging topic. Therefore the following section addressess control of chaotic Coullet dynamic. 4.1 Stabilizing the Chaos The Coullet system is used in a synchronization approach, which is as follows: IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 119 (14) = = = + + + Parameters of the Coullet system are chosen as a = 0.8, b = -1.1, c = -0.45, d = -1 and = 0.95 together with the state initial conditions: x (0) = −0.8, x (0) = 1.2 and x (0) = 0.2. State response of the system (4.2) to initial conditions is simulated in 50 seconds where depicted in Fig. 1 for − . Fourth plot in Fig. 1 represents a phase portrait of the system. Fig. 1: States (x1-x3) responses of together with a phase portrait (fourth). As can be seen from the first three graphs in Fig. 1, state x1-x3 are of chaotic (neither periodic nor converging). The fourth plot also confirms states fail to converge to any point. To deeply verify the chaos, an FFT of x1 is assessed and plotted in Fig. 2, using power GUI facility in the MATLABTM simulink. As can be seen from Fig. 2, all of the frequency components occur in the frequency spectrum. This is in confirmation that the signal is not pure oscillatory. This is main characteristics of chaos (similar to frequency characteristics of noise) where all frequencies can be seen in the chaos [5-6, 37-38]. Although this is for x1, similar results are achieved for x2 and x3. In order to control the chaos, a control effort ( )is applied to the third equation of Eqn.(4.2): (15) = = = + + + + 0 20 40 60 -1 -0.5 0 0.5 1 1.5 Time (s) A m pl itu de X1 State 0 10 20 30 40 50 -1 -0.5 0 0.5 1 1.5 Time (s) A m pl itu de X2 State 0 20 40 60 -1 -0.5 0 0.5 1 Time (s) A m pl itu de X3 State -2 0 2 -2 0 2 -1 0 1 x1x2 x3 IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 120 Fig. 2: Time (the above) and frequency (bottom) response of state x1 of the Coullet system. The control takes places using the synergetic technique. The duty is to provide zero convergence of the state variables. The macro variable is defined as follows [2]: (16) = + + In order to make more degree of freedom, arbitrary gains of , and are innovatively added to the macro variable in Eqn. (16). This generalization expands the macro as in the following: (17) = + + It will be shown that the proposed macro in Eqn. (17) increases the performance. The required control effort can be obtained when the conventional manifold in Eqn. (16) i.e. ψ = x + x + x = 0 is substituted into Eqn. (13), which is as follows: (18) = − 1 ( ) + ( ) + ( + + + ) + ( ) + ( ) + ( ) Likewise, anew control law is derived for the proposed manifold in Eqn. (17) using the following fractional dynamical equation of macro variables: (19) ( ) + ( ) = 0, > 0 Substitution of Eqn. (19) into Eqn. (13) yields: (20) ( ) + ( ) + ( ) + ( ) + ( ) + ( ) = 0 0 10 20 30 40 50 -0.5 0 0.5 1 Selected signal: 50 cycles. FFT window (in red): 1 cycles Time (s) 0 10 20 30 40 50 0 10 20 30 40 50 60 Frequency (Hz) Fundamental (1Hz) = 0.3419 , THD= 78.96% M ag ( % o f F un da m en ta l) IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 121 where ( ) is -order derivative of state variable in the Caputo sense as: (21) ( ) = = / By substituting dynamical equation Eqn. (15) into Eqn. (20), the following equation is obtained: (22) ( ) + ( ) + ( + + + + ) + ( ) + ( ) + ( ) = 0 The control input u is finally extracted form (22) which is as follows (23) = − 1 ( ) + ( ) + ( + + + ) + ( ) + ( ) + ( ) The stability of the proposed controller will be analyzed through the following theorem: Theorem: Consider the Coullet fractional order system in Eqn. (15). If the control law in Eqn. (23) is applied, system states asymptotically converge to zero. Proof: The Lyapunov function is candidate as in the following: (24) = 1 2 (ψ( )) ≥ 0 The conventional derivative yields: (25) V̇ = (ψ( ))( )ψ( ) = ψ( )ψ( ) A fractional form of the macro-variable i.e. ψ( ), is achieved as follows: (26) ( ) = ( ) + ( ) + ( ) = ( ) + ( ) + ( + + + + ) Replacement the control effort ( ) from Eqn. (23) into Eqn. (26) provides: (27) ( ) = ( ) + ( ) + ( + + + − 1 ( ) + ( ) + ( + + + ) + ( ) + ( ) + ( ) = − 1 ( ) + ( ) + ( ) = − 1 ( ) Substituting ψ( ) = − ψ( )from Eqn. (27) into Eqn. (25) achieves: IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 122 (28) V̇ = − 1 ( ) ( ) = − 1 ( ( ) ) ≤ 0 Inequality Eqn. (28) guarantees strictly decreasing nature of the Lyapunov function in Eqn. (22). This confirms that the control input Eqn. (21) is a stabilizing control effort for system Eqn. (15). 5. RESULTS: APPLICATION OF THE PROPOSED SYNERGETIC CONTROLLER TO CONTROL A CHAOS A simulation is carried out for the fractional order Coullet system gaining both the manifold in Eqn. (16) [2] and the proposed manifold in Eqn. (17). Time responses of three states are shown in Fig. 3, Fig. 4 and Fig. 5 when initial conditions are assumed as (0) = −0.8, (0) = 1.2, (0) = 0.2 and T = 0.1 whilst gains are chosen = 3, = 5 and = 1 using Eqn. (3) to Eqn. (6). From Fig. 3 it can be seen that by considering extra gains into Eqn. (16) the macro variable in Eqn. (17) causes reduction in the overshoot. Meanwhile better performance with respect to the conventional macro variable in Eqn. (17) is also achieved. However, errors in two methods approach zero. The time behavior of the second state variable in Fig. 4 is similar to the first state variable as in Fig. 3. In contrast to those two state variables, the third state variable in Fig. 5 provides more negative peak with respect to that of the conventional macro in Eqn. (16). However, the proposed macro in Eqn. (17) provides a faster convergence rate (with the price of more negative peak). This necessitates to tune those gains in an optimal way rather than simple selection of = 3, = 5 and = 1. Fortunately, this probable draw back requires no more attention to provide the control effort which is depicted in Fig. 6 which makes the lack of negative peak negligible. Fig. 3: Time response of state variable x in fractional order Coullet system with conventional macro Eqn. (16) as = + + = 0 and the proposed macro Eqn. (17) as = + + = 0. Fig. 4: Time response of state variable x in fractional order Coullet system with conventional macro Eqn. (16) as = + + = 0 and the proposed macro Eqn. (17) as = + + = 0. 0 5 10 15 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 Time(s) x 1 The proposed method method in [2] 0 5 10 15 -0.2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 Time(s) x 2 The proposed method method in [2] IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 123 Fig. 5: Time response of state variable x in fractional order Coullet system with conventional macro Eqn. (16) as = + + = 0 and the proposed macro Eqn. (17) as = + + = 0. Fig. 6: The required control input for using the conventional macro in Eqn. (16) as = + + = 0 and the proposed macro Eqn. (17) as = + + = 0. 6. SYNERGETIC SYNCHRONIZATION OF FRACTIONAL-ORDER COULLET SYSTEM Basically, chaos synchronization means forcing two systems work in a same way in a master-slave structure. The designed nonlinear control system obtains signals from the master to control the slave dynamics. Block diagram of the master and slave synchronization using the synergetic controller is shown in Fig. 7. Fig. 7: Block diagram of the master and slave synchronization using the synergetic controller. The goal is to synchronize a fractional-order Coullet system assuming the master as in Eqn. (29) to be followed by the slave as in Eqn. (5.2) using the proposed synergetic controller Eqn. (23). (29) = = = + + + Supposing (0) = −0.8, (0) = 1.2 and (0) = 0.2 where the makes the system chaotic. The slave dynamics is similarly shown by: 0 5 10 15 -1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 Time(s) x 3 The proposed method method in [2] 0 5 10 15 -6 -5 -4 -3 -2 -1 0 1 Time(s) co nt ro l e ff or t The proposed method method in [2] IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 124 (30) = = = + + + + Initial conditions of the slave are assumed different as (0) = 1, (0) = −1.2 and (0) = −0.4. The error is defined as discrepancy of the corresponding states, i.e. = − for i = 1,2,3. Deducing the master dynamic from the slave leads to: (31) = = = + + + + The same procedure as in Eqn. (19) to Eqn. (23) using the proposed linear gains leads to generate the following control effort: (32) = − 1 ( ) + ( ) + ( + + + ) + ( ) + ( ) + ( ) A proper selection of control parameters ki, i.e. k1 = 3, k2 = 5.75 and k3 = 1 provides lim → = 0. Outcomes of the synchronization are shown in Fig. 108, Fig. 9 and Fig. 10. The proof is similar to that of in Eqn. (24) to Eqn. (28). This immediately means that output of the slave asymptotically follows the master state i.e. = . This confirms that the synchronization between master and slave is made possible. Fig. 8: Synchronization of states and together with the = − . Fig. 9: Synchronization of states and together with the = − . Fig. 10: Synchronization of states and together with the = − . 0 5 10 15 20 -4 -2 0 2 Time(s) x 1 ,y 1 x 1 y 1 0 5 10 15 20 -4 -2 0 2 Time(s) e 1 =x 1- y 1 0 5 10 15 20 -2 0 2 Time(s) x 2 ,y 2 x 2 y 2 0 5 10 15 20 0 1 2 3 Time(s) e 2 =x 2- y 2 0 5 10 15 20 -2 -1 0 1 Time(s) x3 ,y 3 x3 y3 0 5 10 15 20 -2 -1 0 Time(s) E rr or = x3 -y 3 IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 125 7. DISCUSSION AND CONCLUSION In this paper, a new controller is proposed by considering a fractional derivative in linear weighted combination of state variables in the synergetic controller. The degree of freedom is also increased by adding control coefficients , , and to the designed synergetic manifold. By tuning parameters , and frequency of oscillations is decreased when convergence is finally achieved. The stability of the proposed configuration is proven using the Lyapunov theory. The achievement confirms performance of the proposed structure. The parameters , and are determined by a trial and error process, while they can be determined using intelligent algorithm such as PSO (Particle swarm optimization) or other algorithms to obtain better results. REFERENCES [1] Kolesnikov A, et al. (2000) Modern applied control theory: Synergetic approach in control theory. TRTU, Moscow, Taganrog. [2] Djennoune S, Bettayeb M. (2013) Optimal synergetic control for fractional-order systems. Automatica, 49(7):2243-2249. [3] Jiang Z. (2009) Design of a nonlinear power system stabilizer using synergetic control theory. Electric Power Sys. Res., 79(6):855-862. [4] Jiang Z, Dougal RA. (2004) Synergetic control of power converters for pulse current charging of advanced batteries from a fuel cell power source. IEEE Trans. Power Electronic, 19(4):1140-1150. [5] Santi E, et al. (2003) Synergetic control for DC-DC boost converter: implementation options.IEEETrans.Ind.App., 39(6):1803-1813. http://dx.doi.org/10.1109/TIA.2003.818967 [6] Ni J, et al. (2014) Variable speed synergetic control for chaotic oscillation in power system. Nonlinear Dynamics, 78(1):681-690. [7] Delavari H, Lanusse P, Sabatier J. (2013) Fractional order controller design for a flexible link manipulator robot. Asian J. Control, 15(3):783-795. [8] El-Khazali R. (2013) Fractional-order controller design. Computers & Mathematics with Applications, 66(5):639-646. [9] Faieghi MR, Delavari H, Baleanu D. (2012) Control of an uncertain fractional-order Liu system via fuzzy fractional-order sliding mode control. J. Vib. Control, 18(9), 1366-1374. [10] Kilbas AAA, Srivastava HM, Trujillo JJ. (2006) Theory and applications of fractional differential equations, vol. 204, Elsevier Science Limited. [11] Magin RL. (2006) Fractional calculus in bioengineering. Begell House, Rodding. [12] Mainardi F. (1997) Fractional calculus. Springer. [13] Monje, CA, et al. (2010) Fractional-order systems and controls: fundamentals and applications. Springer. [14] Odibat Z, Momani S. (2006) Application of variational iteration method to nonlinear differential equations of fractional order. Int. J. Nonlinear Sci. Num. Simulation, 7(1):27-34. [15] Podlubny I. (1998) Fractional differential equations: an introduction to fractional derivatives, fractional differential equations, to methods of their solution and some of their applications. vol. 198. Academic press. [16] Sabatier J, Agrawal OP, Machado JT. (2007) Advances in fractional calculus. Springer. [17] Luo Y, et al. (2010) Tuning fractional order proportional integral controllers for fractional order systems. J. Process Control, 20(7):823-831. [18] Meng L, Wang DF, Han P. (2012) Identification of fractional order system using particle swarm optimization. IEEE International Conference on Machine Learning and Cybernetics (ICMLC). [19] Saxena R, Mathai A, Haubold H. (2004) Unified fractional kinetic equation and a fractional diffusion equation. Astrophysics and Space Science, 290(3-4):299-310. [20] Bagley RL, Torvik P. (1983) A theoretical basis for the application of fractional calculus to viscoelasticity. J. Rheology, 27(3):201-210. IIUM Engineering Journal, Vol. 17, No. 1, 2016 Pourmehdi et al. 126 [21] Di Paola M, Zingales M. (2012) Exact mechanical models of fractional hereditary materials. J. Rheology, 56(5):983-1004. [22] Hilfer R, et al. (2000) Applications of fractional calculus in physics. vol. 128. World Scientific. [23] Poinot T, Trigeassou J. (2004) Identification of fractional systems using an output-error technique. Nonlinear Dynamics, 38,133-154. [24] Cois O, Oustaloup A, Battaglia E. (2000) Non-integer model from modal decomposition for time domain system identification. 12th IFAC Symposium on System Identification, Santa Barbara, CA, USA. [25] Sabatier J, et al. (2002) CRONE control: Principles and extension to time-variant plants with asymptotically constant coefficients. Nonlinear Dynamics, 29(1-4):363-385. [26] Podlubny I. (1999) Fractional-order systems and PI/sup/spl lambda//D/sup/spl mu//- controllers. IEEE Trans. Autom. Control, 44(1):208-214. [27] Nechadi E, et al. (2012) Type-2 fuzzy based adaptive synergetic power system control. Electric Power Sys. Res., 88:9-15. [28] Santi E, et al. (2004) Synergetic control for power electronics applications: a comparison with the sliding mode approach. J. Circ. Sys. Comp., 13(4):737-760. [29] Caputo M. (1969) Elasticita´ e dissipazione. Zanichelli. Bologna. [30] Caputo M. (1989) The rheology of an anelastic medium studied by means of the observation of the splitting of its eigenfrequencies. J. Acoustical Soc. America, 86(5):1984-1987. [31] El-Sayed AM. (1995) Fractional order evolution equations. J. Fract. Calc, 7:89-100. [32] El-Sayed AM, Ibrahim A. (1995) Multivalued fractional differential equations. App. Math. Comp., 68(1):15-25. [33] Vinagre B, et al. (2002) Using fractional order adjustment rules and fractional order reference models in model-reference adaptive control. Nonlinear Dynamics, 29(1-4):269-279. [34] Xue D, Zhao C, Chen Y. (2006) Fractional order PID control of a DC-motor with elastic shaft: a case study. Proceedings of American Control Conference. [35] Jesus IS, Machado JT. (2008) Fractional control of heat diffusion systems. Nonlinear Dynamics, 54(3):263-282. [36] Arneodo A, Coullet P, Tresser C. (1981) Possible new strange attractors with spiral structure. Comm. Math. Phys., 79(4):573-579. [37] Liu Y, Chen L. (2013) Chaos in attitude dynamics of spacecraft. Springer. [38] Petržela J, Kolka Z, Hanus S. (2011) Simple chaotic oscillator: from mathematical model to practical experiment. Models and Applications of Chaos Theory in Modern Sciences, p. 317.