Advances in Mathematical Physics

Volume 2013 (2013), Article ID 186037, 6 pages

http://dx.doi.org/10.1155/2013/186037

Research Article

The Proposed Modified Liu System with Fractional Order

Alireza K. Golmankhaneh,¹ Roohiyeh Arefi,¹ and Dumitru Baleanu2,3,4

¹Department of Physics, Urmia Branch, Islamic Azad University, P.O. Box 969, Oromiyeh, Iran

²Department of Mathematics and Computer Science, Çankaya University, 06530 Ankara, Turkey

³Institute of Space Sciences, P.O. Box MG-23, 76900 Magurele-Bucharest, Romania

⁴Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, P.O. Box 80204, Jeddah 21589, Saudi Arabia

Received 27 February 2013; Accepted 15 March 2013

Academic Editor: J. A. Tenreiro Machado

Copyright © 2013 Alireza K. Golmankhaneh et al. 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.

The chaos in a new system with order 3 is studied. We have shown that this chaotic system again will be chaotic when the order of system is less than 3. Generalized Adams-Bashforth algorithm has been used for investigating in stability of fixed points and existence of chaos.

1. Introduction

It is well known that the nonlinear equations of dynamical systems with special condition have chaotic behavior [1]. Subsequently, additional studies were performed on the chaos and chaotic systems. Therefore solutions of different systems display their chaotic behavior such as Chen's system, Chua's dynamical system, the motion of double pendulum, and Rossler system amongst others. At first, it was thought that the chaos exists only when the order of system of differential equation is exactly 3. When the system of differential equation is composed of three first order differential equations, the order of the system is the sum of orders. But later on, a very interesting thing was realized; that is, it is also possible to observe chaotic behavior in a fractional order system. The system is composed of differential equations with fractional order derivatives [1–28]. For example, Sheu et al. reviewed the chaotic behavior of the Newton-Leipnik system with fractional order [10]. The important thing in the study of fractional-order systems is the minimum effective dimension of the system for which the system remains chaotic. This quantity has been numerically calculated for different systems including fractional order Lorenz system [11], fractional order Chua's system [12], and fractional order Rössler system [13]. Recently, the chaos has been studied in fractional ordered Liu system, where the numerical investigations on the dynamics of this system have been carried out, and properties of the system have been analyzed by means of Lyapunov exponent [14]. In this paper we study the chaotic behavior of a generalization of the Liu system with fractional order.

The framework of the paper is as follows.

In Section 2, we study the behavior of a new fractional order system (modification of Liu system), and we study commensurate and incommensurate ordered systems and find lowest order at which chaos exists by numerical experiments. We have investigated the instability of fixed points and used Lyapunov exponent for the existence of chaos. In Section 3, we state the main conclusions.

2. The Proposed Modified Liu System

In this section, we review the condition for asymptotic stability of the commensurate and incommensurate fractional ordered systems. We suggest the readers to see [15–21] for the following section.

2.1. Asymptotic Stability of the Commensurate Fractional Ordered System

Let us consider the commensurate fractional ordered dynamical system equation

An equilibrium point of (1) is which is , , and a small disturbance from a fixed point is given by . Thus, we have

Write the system in the matrix form

where

where is Jacobian matrix of the system, and if does not have purely imaginary eigenvalues, therefore the trajectories of the nonlinear system in the neighborhood of the equilibrium point have the same form as the trajectories of the linear system [18]. So we arrive at the following linear autonomous system:

where is matrix and . The system (5) is asymptotically stable if and only if for all eigenvalues of . So for this condition the solutions of (5) tend to as . Therefore, the equilibrium point of the system is asymptotically stable if , for all eigenvalues of . For example,

2.2. Asymptotic Stability of the Incommensurate Fractional Ordered System

Consider the following incommensurate fractional ordered dynamical system [19, 20]. Now suppose the commensurate fractional ordered dynamical system equation

where . One can write it as , , so , are positive integers. The definition of common multiple of ’s is . Equilibrium point and small disturbance of the system and respectively as above. So we get

where it can be written as

where is Jacobian matrix evaluated at point . The definition is

that is, if all the roots of equation satisfy the condition [21], the solution of linear system is asymptotically stable as follows:

The term in the left side of (11) is an instability measure for equilibrium point in fractional ordered system (IMFOS). Then, fractional order equation (7) exhibits chaotic attractor if the condition is [19, 20]

2.3. Modified Liu System

In this section we introduce the following system and show that the system is chaotic:

where , , , , and with the initial conditions lead to the chaotic trajectories. Meanwhile, we want to show chaotic behavior of (13) involving fractional order. Also, we will calculate the minimum effective dimension by which the system remains chaotic. The corresponding fractional order system is

where . In (14) if we choose , the system is called commensurate, and otherwise it is incommensurate. Now, we have four real equilibrium points for (13) which are shown in Table 1. In Table 1, we see the equilibrium points and the eigenvalues of the corresponding Jacobian matrix

A saddle point is stable, if the Jacobian matrix has at least one eigenvalue with a negative real part. Otherwise, one eigenvalue with a nonnegative real part is called unstable. And saddle points have index one or two if there is exactly one or two unstable eigenvalues, respectively. It is established in the literature [22–26] that scrolls are generated only around the saddle points of an index one or two. Saddle points of index one are to connect scrolls. Table 1 shows that the equilibrium points and are saddle points of index two; therefor we have two-scroll attractor [22], in the fractional system given by (14).

Table 1: Equilibrium points and corresponding eigenvalues.

Equilibrium point
Eigenvalue
Nature

()
Saddle point

()
Saddle point

()
Saddle point

()
Saddle point

2.4. Commensurate Ordered System

Consider the system equation (14), and let , so it is called commensurate order for this case. In this case a system shows regular behavior if it satisfies , then we have [15–21]

From Figure 1 we can see that the Lyapunov exponent for the case of commensurate ordered equation (14) is positive if [27, 28]. Thus we can realize that the system does not indicate chaotic behavior for the value . This corollary has been figured out by numerical outcome. Moreover, Figure 2 illustrates the phase portrait in -plane for the . Numerical experiments and Figures 3 and 4 demonstrate that the system has chaotic behavior for . In addition, Figures 5 and 6 represent solutions and for , respectively. Adams-Bashforth predictor-corrector algorithm is used for numerical result with step size .

Figure 1: Lyapunov exponent.

Figure 2: Phase portrait .

Figure 3: Phase portrait x-y  .

Figure 4: Phase portrait y-z .

Figure 5: Solution for .

Figure 6: Solution for .

2.5. Incommensurate Ordered System

In this section we indicate that the condition for being chaotic system in the case of commensurate is not sufficient for the incommensurate case. So let us consider the fractional order system equation (14). Figures 7 and 8 display that the Lyapunov exponent is positive for , , and for the case , respectively. Now, let us consider the following cases for (14).

Figure 7: Lyapunov exponent-.

Figure 8: Lyapunov exponent-.

(1) In the first condition we choose . Therefore, one can acquire = LCM. Since we have , then

And the instability measure for equilibrium point in fractional ordered system (IMFOS) will be for the system

In Figure 9, there is no chaos, whereas IMFOS [15–21]. This consequence leads to that the IMFOS is not sufficient for existence of chaos.

Figure 9: Phase portrait x-y .

(2) As second case, suppose , so that in this case the = LCM. Also we have , as well as we obtain

Therefore, the IMFOS of the system for this state is

If we look at Figure 10 we conclude the chaotic behavior for the system with the mentioned condition. Here we remark the lowest dimension of the system which has chaos (Figure 11).

Figure 10: Phase portrait x-z .

Figure 11: Phase portrait y-z .

(3) As third case, let . So with some manipulation we get = LCM,

If we compute the IMFOS for the system, we will arrive at

Figure 11 and (22) point that the system does not have chaotic condition.

(4) For the fourth case, we consider . Therefore, we will obtain = LCM,

IMFOS for the fourth case is

Numerical results and Figure 12 indicate for the this case that the system is chaotic.

Figure 12: Phase portrait y-z .

3. Conclusions

A fractional order of a new system is investigated. Numerical calculations are performed for different values of fractional order. Lyapunov exponents and analytical conditions given in the literature have been used to check the existence of chaos. A minimum effective dimension is calculated for commensurate fractional order. Mathematica 7 has been used for computations in this paper.

References

E. N. Lorenz, “Deterministic nonperiodic flow,” Journal of The Atmospheric Sciences, vol. 20, pp. 130–141, 1963.

J. Lü, T. Zhou, G. Chen, and S. Zhang, “Local bifurcations of the Chen system,” International Journal of Bifurcation and Chaos, vol. 12, no. 10, pp. 2257–2270, 2002.

L.J.M. Kocić, S. Gegovka-Zajkovka, and S. Kostadinova, “On Chua dynamical system,” Applied Mathematics, Informatics & Mechanics, Series A, vol. 2, pp. 53–60, 2010.

T. Stachowiak and T. Okada, “A numerical analysis of chaos in the double pendulum,” Chaos, Solitons and Fractals, vol. 29, no. 2, pp. 417–422, 2006.

W. Xuedi and W. Chen, “Bifurcation analysis and control of the Rossler,” in Proceedings of the 7th International Conference on System Natural Computation (ICNC '11), vol. 3, pp. 1484–1488, IEEE, Shanghai, China, 2011.

K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, John Wiley & Sons, New York, NY, USA, 1993.

C. M. Ionescu and R. De Keyser, “Relations between fractional-order model parameters and lung pathology in chronic obstructive pulmonary disease,” IEEE Transactions on Biomedical Engineering, vol. 56, no. 4, pp. 978–987, 2009.

J. A. T. Machado, “Entropy analysis of integer and fractional dynamical systems,” Nonlinear Dynamics, vol. 62, no. 1-2, pp. 371–378, 2010.

A. M. Lopes, J. A. T. Machado, C. M. A. Pinto, and A. M. S. F. Galhano, “Fractional dynamics and MDS visualization of earthquake phenomena,” Computers and Mathematics with Applications, 2013.

L.-J. Sheu, H.-K. Chen, J.-H. Chen et al., “Chaos in the Newton-Leipnik system with fractional order,” Chaos, Solitons & Fractals, vol. 36, no. 1, pp. 98–103, 2008.

I. Grigorenko and E. Grigorenko, “Chaotic dynamics of the fractional Lorenz system,” Physical Review Letters, vol. 91, no. 3, pp. 034101/1–034101/4, 2003.

T. T. Hartley, C. F. Lorenzo, and H. K. Qammer, “Chaos in a fractional order Chua's system,” IEEE Transactions on Circuits and Systems I, vol. 42, no. 8, pp. 485–490, 1995.

C. Li and G. Chen, “Chaos and hyperchaos in the fractional-order Rössler equations,” Physica A, vol. 341, no. 1–4, pp. 55–61, 2004.

V. Daftardar-Gejji and S. Bhalekar, “Chaos in fractional ordered Liu system,” Computers & Mathematics with Applications, vol. 59, no. 3, pp. 1117–1127, 2010.

I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.

S. G. Samko, A. A. Kilbas, and O. I. Marichev, Fractional Integrals and Derivatives, Gordon and Breach, New York, NY, USA, 1993.

K. Diethelm, N. J. Ford, and A. D. Freed, “A predictor-corrector approach for the numerical solution of fractional differential equations,” Nonlinear Dynamics, vol. 29, no. 1–4, pp. 3–22, 2002.

Z. Vukic, Lj. Kuljaca, D. Donlagic, and S. Tesnjak, Non-linear Control Systems, Marcel Dekker, New York, NY, USA, 2003.

A. Razminia, V.J. Majd, and D. Baleanu, “Chaotic incommensurate fractional order Rössler system: active control and synchronization,” Advances in Difference Equations, article 15, 2011.

M. S. Tavazoei and M. Haeri, “Chaotic attractors in incommensurate fractional order systems,” Physica D, vol. 237, no. 20, pp. 2628–2637, 2008.

W. Deng, C. Li, and J. Lü, “Stability analysis of linear fractional differential system with multiple time delays,” Nonlinear Dynamics, vol. 48, no. 4, pp. 409–416, 2007.

M. S. Tavazoei and M. Haeri, “A necessary condition for double scroll attractor existence in fractional-order systems,” Physics Letters A, vol. 367, no. 1-2, pp. 102–113, 2007.

L. O. Chua, M. Komuro, and T. Matsumoto, “The double scroll family. I. Rigorous proof of chaos,” IEEE Transactions on Circuits and Systems, vol. 33, no. 11, pp. 1072–1097, 1986.

C. P. Silva, “Shilnikov theorem—a tutorial,” IEEE Transactions on Circuits and Systems I, vol. 40, no. 10, pp. 675–682, 1993.

D. Cafagna and G. Grassi, “New 3D-scroll attractors in hyperchaotic Chua's circuits forming a ring,” International Journal of Bifurcation and Chaos, vol. 13, no. 10, pp. 2889–2903, 2003.

J. Lü, G. Chen, X. Yu, and H. Leung, “Design and analysis of multiscroll chaotic attractors from saturated function series,” IEEE Transactions on Circuits and Systems I, vol. 51, no. 12, pp. 2476–2490, 2004.

A. Wolf, J. B. Swift, H. L. Swinney, and J. A. Vastano, “Determining Lyapunov exponents from a time series,” Physica D, vol. 16, no. 3, pp. 285–317, 1985.

M. T. Rosenstein, J. J. Collins, and C. J. De Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,” Physica D, vol. 65, no. 1-2, pp. 117–134, 1993.