The Scientific World Journal

Volume 2013 (2013), Article ID 645080, 13 pages

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

## Stability Analysis of Distributed Order Fractional Chen System

^{1}Department of Applied Mathematics, School of Mathematical Sciences, University of Guilan, P.O. Box 1914, Rasht, Iran^{2}Department of Applied Mathematics, Faculty of Mathematical Sciences, Islamic Azad University, Lahijan Branch, P.O. Box 1616, Lahijan, Iran

Received 31 August 2013; Accepted 3 October 2013

Academic Editors: C. Li, F. Liu, and R. Magin

Copyright © 2013 H. Aminikhah 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

We first investigate sufficient and necessary conditions of stability of nonlinear distributed order fractional system and then we generalize the integer-order Chen system into the distributed order fractional domain. Based on the asymptotic stability theory of nonlinear distributed order fractional systems, the stability of distributed order fractional Chen system is discussed. In addition, we have found that chaos exists in the double fractional order Chen system. Numerical solutions are used to verify the analytical results.

#### 1. Introduction

The history of fractional calculus is more than three centuries old, yet only in the past 20 years has the field received much attention and interest. The reader may refer to [1, 2]. The generalization of dynamical equations using fractional derivatives proved to be useful and more accurate in mathematical modeling related to many interdisciplinary areas. Applications of fractional calculus and fractional-order differential equations include dielectric relaxation phenomena in polymeric materials [3], transport of passive tracers carried by fluid flow in a porous medium in groundwater hydrology [4], transport dynamics in systems governed by anomalous diffusion [5, 6], long-time memory in financial time series [7], and so on [8, 9]. Stability analysis and control systems are two of the most important problems such that, in 1996, Matignon [10] firstly studied stability of -dimensional linear fractional systems from a point of view of control. Since then, many researchers have completed further studies on the stability of linear fractional differential systems [11–14]. For the nonlinear fractional differential systems, the stability analysis is much more difficult and only a few are available. Some authors [15, 16] studied the nonlinear fractional differential system.

Nonlinear systems are very interesting to engineers, physicists, and mathematicians because most real physical systems are inherently nonlinear in nature. Nonlinear dynamical systems can exhibit completely an unpredictable behavior, the so-called deterministic chaos. Chaos is an important dynamical phenomenon that has been extensively studied and developed by scientists and engineers in the last decades in various fields such as physical [17], chemical [18], and ecological systems [19]. Since chaos is useful and has potential applications in many technological disciplines, the discovery and the creation of chaos are important. In 1963, Lorenz found the first chaotic attractor in a three-dimensional autonomous system [20]. Later, several dynamical systems exhibiting chaos have been presented in various branches of science [21]. For example, in 1999, Chen and Ueta found another simple three-dimensional autonomous system, which is not topologically equivalent to Lorenz’s system and which has a chaotic attractor [22, 23]. Similar to integer-order chaotic systems, the fractional order chaotic systems have also interested several researchers. I. Grigorenko and E. Grigorenko extended the study of this prototypical system to equations of fractional order but, unfortunately, the results presented in this paper are not correct [24]. In [25, 26], the chaotic behavior and its control in the fractional order Chen system are investigated. Also, chaotic behaviors have been found in the fractional order systems of Chuna [27], Rossler [28], Coullet [29], modified Van der Pol-Duffing [30], and Liu [31]. In addition, the fractional order Chen system has been studied with time-delay in [32].

The idea of fractional derivative of distributed order is stated by Caputo [33] and later developed by Caputo himself [34, 35], and Bagley and Torvik [36, 37]. Other researchers used this idea, and interesting reviews appeared to describe the related mathematical models of partial fractional differential equation of distributed order. For example, Diethelm and Ford [38] used a numerical technique along with its error analysis to solve the distributed order differential equation and analyze the physical phenomena and engineering problems; see [38] and references therein. Recently Saberi Najafi et al. [39, 40] studied stability analysis of distributed order fractional differential equations with respect to the nonnegative density function. The aim of the present work is twofold. First, we consider the stability of -dimensional nonlinear distributed order fractional differential system with respect to the nonnegative density function and then we study the stability of distributed order fractional Chen system.

This paper is organized as follows. In the next section, we recall some basic definitions of the Caputo fractional derivative operator, systems with fractional derivatives of distributed order. Section 3 contains the main definitions and theorems for checking the stability analysis of distributed order fractional system. In Section 4, we present the distributed order fractional Chen system and afterwards based on the stability theorem of distributed order fractional systems, the stability of the distributed order fractional Chen system is discussed. Finally, the numerical solution to illustrate the validity of the results is presented in Section 5.

#### 2. Elementary Definitions

In this section, we consider the main definitions and properties of fractional derivative operators of single and distributed order.

The differential and integral operator in fractional calculus is denoted by , where and are the bounds of the operation and is the fractional order, which can be rational, irrational, or even complex. For simplicity and without loss of generality, in the following, we assume that and . The continuous integrodifferential operator is defined as follows: There are many definitions of fractional derivatives of order [1, 2], such as Grunwald-Letnikov’s definition (GL), Riemann-Liouville’s definition (RL), and Caputo’s fractional derivative. The RL definition is given as where is the first integer which is not less than ; that is, and is a Gamma function. The Caputo fractional derivative of is defined as Finally, Grunwald-Letnikov definition is given by Fortunately, the Laplace transform of the Caputo fractional derivative satisfies where and is the Laplace variable. The Laplace transform of Caputo fractional derivative requires the knowledge of the initial values of the function and its integer derivatives of order . When , is given by When the fractional calculus operators act on , and we integrate with respect to the order, then distributed order fractional differential equations can be obtained. In this brief, the following distributed order fractional differential operator notation is adopted: where the derivative is taken to be a fractional derivative of Caputo type of order with respect to the nonnegative density function .

The idea of distributed order is stated by Caputo [33, 34]. Further the Laplace transform of the Caputo distributed order satisfies where is the Laplace transform of and A distributed order fractional equations can be defined by the following model: where and .

In [38] the following results about the existence and uniqueness of solutions for (10) are further presented.

Theorem 1. *Let the function be absolutely integrable on the interval and satisfy for and and is such that for for every ; then initial value problem (10) has a unique solution.*

Furthermore, the above definition in one dimension can naturally be generalized to the case of multiple dimensions; that is, let and , . The -dimension distributed order fractional system is described as follows: where The result of Theorem 1 can be easily generalized to the initial value problem of (11).

#### 3. Stability Analysis of Distributed Order Fractional Systems

In this section, we generalize the main stability properties for systems of distributed order fractional. The linear distributed order fractional systems are expressed as where , the matrix , and , , . Then Saberi Najafi et al. [39] have obtained the general solution of the distributed order fractional systems (13), which is written by where , , , and .

Now, we will recall some theorems and definitions about linear distributed order fractional equations and then we will show this theorem for nonlinear distributed order fractional equations as well.

Theorem 2 (see [39]). *The distributed order fractional system of (13) is asymptotically stable if and only if all roots of have negative real parts.*

*Definition 3 (see [39]). *The value of is the characteristic function of the matrix with respect to the distributed function , where is the distributed function with respect to the density function .

The inertia of a matrix is the triplet of the numbers of eigenvalues of with positive, negative, and zero real parts. As pointed out in [39], authors have generalized the inertia concept for analyzing the stability of linear distributed order fractional systems.

*Definition 4. *The inertia of the system (13) is the triple
where , , and are, respectively, the number of roots of with positive, negative, and zero real parts, where is the distributed function with respect to the density function .

Theorem 5 (see [39]). *The linear distributed order fractional system (13) is asymptotically stable if and only if any of the following equivalent conditions holds:*(1)*,
*(2)*all roots of the characteristic function of with respect to satisfy .*

Next, we will mainly discuss the stability of a nonlinear autonomous distributed order fractional system, which can be described by with the initial value , where , and , .

Theorem 6. *Let be the equilibrium of system (16); that is, and is the Jacobian matrix at the point ; then the point is asymptotically stable if and only if all roots of the characteristic function of with respect to satisfy .*

*Proof. *Let , where is a small disturbance from a fixed point. Therefore
since and ; thus, we have
System (18) can be written as
with the initial value .

The analytical procedure of linearization is based on the fact that if the matrix has no purely imaginary eigenvalues, then the trajectories of the nonlinear system in the neighborhood of the equilibrium point have the same form as the trajectories of the linear system [41]. Hence, by applying Theorem 2, the linear system (20) is asymptotically stable if and only if all roots of the characteristic function of with respect to satisfy , which implies that the equilibrium of the nonlinear distributed order fractional system (16) is as asymptotically stable.

*Remark 7. *The nonlinear distributed order fractional system (16) in the point is asymptotically stable if and only if .

#### 4. Distributed Order Fractional Chen System

The Chen system is described by the following nonlinear differential equations on [22, 23]: where , , and are the state variables and , , and are three system parameters. The above system has a chaotic attractor when , , and as shown in Figure 1. The corresponding distributed order fractional Chen system (21) can be written in the form where for denote the nonnegative density function of order . As a generalization of nonlinear fractional order differential equation into nonlinear distributed order fractional differential equation, the linearized form of the system (22) at the equilibrium points , , and , that is, , can be written in the form where , , , and for . The Jacobian matrix of distributed order fractional Chen system (22) at the equilibrium point is given by

*Remark 8. *If , where for and is the Dirac delta function, then, we have the following fractional incommensurate-order Chen system [25]:

Based upon Theorem 6, the stability of the distributed order fractional Chen system can be reached with ease. For analyzing the stability of the distributed order fractional Chen system, we compute in the case that the density function varies. The results are shown in Table 1 for some parameters , , and .

#### 5. Numerical Methods

As pointed out in [38], distributed order fractional differential equations may be regarded as a generalization of single-term fractional differential equations or multiterm fractional differential equations Therefore, in this section, numerical method to solving of (26) and (27) is presented. The approximate Grunwald-Letnikov’s definition is given below, where the step size of is assumed to be very small [2, 42, 43] where and are binomial coefficients, which can be computed as [43] Then, general numerical solution of (26) and (27) can be expressed as where for are binomial coefficients calculated according to (29). Equations in (30) can be rewritten as the following forms: Based on the previous algorithm, we can obtain the numerical solution of the fractional differential equations (26) and multiterm fractional differential equations (27), where for and where is the total time of the calculation.

To verify the efficiency of the obtained results in Table 1, the numerical solution for the distributed order fractional Chen system has been computed. In the following calculations, let , with the initial conditions .

Figure 2 shows that the system (22) with parameters and = is asymptotically stable in the equilibrium . Figure 3 demonstrates that the system (22) with parameters and is asymptotically stable in the equilibrium . From Figure 4, we can see that for the distributed order Chen system (21) with parameters and is chaotic. Figure 5 demonstrates that the system (22) with parameters and = has chaotic attractor. Figure 6 shows the chaotic attractor for the distributed order Chen system (22) for parameters and = . From Figure 7, we can see that the system (22) is asymptotically stable in the equilibrium with parameters and = .

#### 6. Conclusion

In this paper, we introduced the nonlinear distributed order fractional differential equations with respect to a nonnegative density function; hence the asymptotical stability for such systems has been investigated. In addition, we presented the distributed order fractional Chen system and then in two special cases the stability for the distributed order fractional Chen is discussed. Numerical solutions were coincident with results of Table 1 described in Section 4. All numerical results are obtained using MATLAB 7.8.

#### References

- A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo,
*Theory and Application of Fractional Differential Equations*, Elsevier, New York, NY, USA, 2006. - K. B. Oldham and J. Spanier,
*The Fractional Calculus*, Acadmic Press, New York, NY, USA, 1974. - E. Reyes-Melo, J. Martinez-Vega, C. Guerrero-Salazar, and U. Ortiz-Mendez, “Application of fractional calculus to the modeling of dielectric relaxation phenomena in polymeric materials,”
*Journal of Applied Polymer Science*, vol. 98, no. 2, pp. 923–935, 2005. View at Publisher · View at Google Scholar · View at Scopus - R. Schumer, D. A. Benson, M. M. Meerschaert, and S. W. Wheatcraft, “Eulerian derivation of the fractional advection-dispersion equation,”
*Journal of Contaminant Hydrology*, vol. 48, no. 1-2, pp. 69–88, 2001. View at Publisher · View at Google Scholar · View at Scopus - B. I. Henry and S. L. Wearne, “Existence of turing instabilities in a two-species fractional reaction-diffusion system,”
*SIAM Journal on Applied Mathematics*, vol. 62, no. 3, pp. 870–887, 2002. View at Publisher · View at Google Scholar · View at Scopus - R. Metzler and J. Klafter, “The random walk's guide to anomalous diffusion: a fractional dynamics approach,”
*Physics Report*, vol. 339, no. 1, pp. 1–77, 2000. View at Google Scholar · View at Scopus - S. Picozzi and B. J. West, “Fractional Langevin model of memory in financial markets,”
*Physical Review E*, vol. 66, no. 4, Article ID 046118, pp. 046118/1–046118/12, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. Ansari, A. Refahi Sheikhani, and S. Kordrostami, “On the generating function ${e}^{xt+y\phi \left(t\right)}$ and its fractional calculus,”
*Central European Journal of Physics*, 2013. View at Publisher · View at Google Scholar - A. Ansari, A. Refahi Sheikhani, and H. Saberi Najafi, “Solution to system of partial fractional differential equations using the fractional exponential operators,”
*Mathematical Methods in the Applied Sciences*, vol. 35, no. 1, pp. 119–123, 2012. View at Publisher · View at Google Scholar · View at Scopus - D. Matignon, “Stability results for fractional differential equations with applications to control processing,” in
*Proceedings of the Computational Engineering in Systems and Application Multiconference (IMACS, IEEE-SMC '96)*, vol. 2, pp. 963–968, Lille, France, 1996. - M. S. Tavazoei and M. Haeri, “A note on the stability of fractional order systems,”
*Mathematics and Computers in Simulation*, vol. 79, no. 5, pp. 1566–1576, 2009. View at Publisher · View at Google Scholar · View at Scopus - V. Daftardar-Gejji and A. Babakhani, “Analysis of a system of fractional differential equations,”
*Journal of Mathematical Analysis and Applications*, vol. 293, no. 2, pp. 511–522, 2004. View at Publisher · View at Google Scholar · View at Scopus - K. Diethelm and N. J. Ford, “Analysis of fractional differential equations,”
*Journal of Mathematical Analysis and Applications*, vol. 265, no. 2, pp. 229–248, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. Sabatier, M. Moze, and C. Farges, “LMI stability conditions for fractional order systems,”
*Computers and Mathematics with Applications*, vol. 59, no. 5, pp. 1594–1609, 2010. View at Publisher · View at Google Scholar · View at Scopus - Y. Li, Y. Chen, and I. Podlubny, “Mittag-Leffler stability of fractional order nonlinear dynamic systems,”
*Automatica*, vol. 45, no. 8, pp. 1965–1969, 2009. View at Publisher · View at Google Scholar · View at Scopus - Y. Li, Y. Chen, and I. Podlubny, “Stability of fractional-order nonlinear dynamic systems: lyapunov direct method and generalized Mittag-Leffler stability,”
*Computers and Mathematics with Applications*, vol. 59, no. 5, pp. 1810–1821, 2010. View at Publisher · View at Google Scholar · View at Scopus - M. Lakshmanan and K. Murali,
*Chaos in Nonlinear Oscillators: Controlling and Synchronization*, World Scientific, Singapore, 1996. - S. K. Han, C. Kurrer, and Y. Kuramoto, “Dephasing and bursting in coupled neural oscillators,”
*Physical Review Letters*, vol. 75, no. 17, pp. 3190–3193, 1995. View at Publisher · View at Google Scholar · View at Scopus - B. Blasius, A. Huppert, and L. Stone, “Complex dynamics and phase synchronization in spatially extended ecological systems,”
*Nature*, vol. 399, no. 6734, pp. 354–359, 1999. View at Publisher · View at Google Scholar · View at Scopus - C. Sparrow,
*The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors*, Springer, New York, NY, USA, 1982. - K. T. Alligood, T. D. Sauer, and J. A. Yorke,
*Chaos: An Introduction to Dynamical Systems*, Springer, New York, NY, USA, 2008. - G. Chen and T. Ueta, “Yet another chaotic attractor,”
*International Journal of Bifurcation and Chaos*, vol. 9, no. 7, pp. 1465–1466, 1999. View at Publisher · View at Google Scholar · View at Scopus - T. Ueta and G. Chen, “Bifurcation analysis of Chen's equation,”
*International Journal of Bifurcation and Chaos in Applied Sciences and Engineering*, vol. 10, no. 8, pp. 1917–1931, 2000. View at Google Scholar · View at Scopus - 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. View at Google Scholar · View at Scopus - C. P. Li and G. J. Peng, “Chaos in Chen's system with a fractional order,”
*Chaos, Solitons and Fractals*, vol. 22, no. 2, pp. 443–450, 2004. View at Publisher · View at Google Scholar · View at Scopus - C. G. Li and G. Chen, “Chaos in the fractional order Chen system and its control,”
*Chaos, Solitons and Fractals*, vol. 22, no. 3, pp. 549–554, 2004. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - C. G. 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. View at Publisher · View at Google Scholar · View at Scopus - M. Shahiri, R. Ghaderi, A. Ranjbar N., S. H. Hosseinnia, and S. Momani, “Chaotic fractional-order Coullet system: synchronization and control approach,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 15, no. 3, pp. 665–674, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. E. Matouk, “Chaos, feedback control and synchronization of a fractional-order modified Autonomous Van der Pol-Duffing circuit,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 2, pp. 975–986, 2011. View at Publisher · View at Google Scholar · View at Scopus - X. Y. Wang and M. J. Wang, “Dynamic analysis of the fractional-order Liu system and its synchronization,”
*Chaos*, vol. 17, no. 3, Article ID 033106, 2007. View at Publisher · View at Google Scholar · View at Scopus - V. Daftardar-Gejji, S. Bhalekar, and P. Gade, “Dynamics of fractional-ordered Chen system with delay,”
*Journal of Physics*, vol. 79, pp. 61–69, 2012. View at Google Scholar - M. Caputo,
*Elasticita e Dissipazione*, Zanichelli, Bologna, Italy, 1969. - M. Caputo, “Mean fractional-order-derivatives differential equations and filters,”
*Annali dell'Università di Ferrara*, vol. 41, no. 1, pp. 73–84, 1995. View at Publisher · View at Google Scholar · View at Scopus - M. Caputo, “Distributed order differential equations modelling dielectric induction and diffusion,”
*Fractional Calculus and Applied Analysis*, vol. 4, pp. 421–442, 2001. View at Google Scholar - L. Bagley and P. J. Torvik, “On the existence of the order domain and the solution of distributed order equations,”
*International Journal of Applied Mathematics*, vol. 1, pp. 865–882, 2000. View at Google Scholar - R. L. Bagley and P. J. Torvik, “On the existence of the order domain and the solution of distributed order equations,”
*International Journal of Applied Mathematics*, vol. 2, pp. 965–987, 2000. View at Google Scholar - K. Diethelm and N. J. Ford, “Numerical analysis for distributed-order differential equations,”
*Journal of Computational and Applied Mathematics*, vol. 225, no. 1, pp. 96–104, 2009. View at Publisher · View at Google Scholar · View at Scopus - H. Saberi Najafi, A. Refahi Sheikhani, and A. Ansari, “Stability analysis of distributed order fractional differential equations,”
*Abstract and Applied Analysis*, vol. 2011, Article ID 175323, 12 pages, 2011. View at Publisher · View at Google Scholar · View at Scopus - A. Refahi Sheikhani, H. Saberi Nadjafi, A. Ansari, and F. Mehrdoust, “Analytic study on linear system of distributed- order fractional differntional equations,”
*Le Matematiche*, vol. 67, p. 313, 2012. View at Google Scholar - Z. Vukic and L. Kuljaca,
*Nonlinear Control Systems*, CRC Press, New York, NY, USA, 2003. - L. Dorciak,
*Numerical Models for Simulation the Fractional-Order Control Systems, UEF-04-94*, The Academy of Sciences, Institute of Experimental Physic, Kosiice, Slovak Republic, 1994. - B. M. Vinagre, Y. Q. Chen, and I. Petráš, “Two direct Tustin discretization methods for fractional-order differentiator/integrator,”
*Journal of the Franklin Institute*, vol. 340, no. 5, pp. 349–362, 2003. View at Publisher · View at Google Scholar · View at Scopus