Discrete Dynamics in Nature and Society

Volume 2012 (2012), Article ID 427393, 12 pages

doi:10.1155/2012/427393

Research Article

Application of Multistep Generalized Differential Transform Method for the Solutions of the Fractional-Order Chua's System

Asad Freihat and Shaher Momani

Department of Mathematics, Faculty of Science, The University of Jordan, Amman 1194, Jordan

Received 29 July 2012; Accepted 18 October 2012

Academic Editor: Mingshu Peng

Copyright © 2012 Asad Freihat and Shaher Momani. 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 numerically investigate the dynamical behavior of the fractional-order Chua's system. By utilizing the multistep generalized differential transform method (MSGDTM), we find that the fractional-order Chua's system with “effective dimension” less than three can exhibit chaos as well as other nonlinear behavior. Numerical results are presented graphically and reveal that the multistep generalized differential transform method is an effective and convenient method to solve similar nonlinear problems in fractional calculus.

1. Introduction

The characterization of real dynamical systems using fractional-order dynamical models has proved to be superior to the traditional calculus [1]. Many generalized fundamentals were extracted and were reduced to their known responses when the fractional orders converge to integer values. Fractional derivatives provide an excellent instrument to describe memory and hereditary properties of various materials and processes. Fractional differentiation and integration operators are used to model problems in astrophysics [2–5], chemical physics, signal processing, systems identification, control and robotics [6, 7], and many other areas [8–15].

It is well known that chaos cannot occur in continuous systems of total order less than three. This assertion is based on the usual concepts of order, such as the number of states in a system or the total number of separate differentiations or integrations in the system. The model of system can be rearranged to three single differential equations, where the equations contain the noninteger (fractional) order derivative. The total order of system is changed from 3 to the sum of each particular order. Petrás [16] has studied the dynamics of fractional-order Chua’s system. The key finding of his study is that chaotic behavior exists in the fractional-order Chua’s system with total order less than three. Hartley et al. [17] introduced a particular form of the fractional-order Chua’s system which is known as the Chua-Hartley’s system. This system is different from the usual Chua’s system in that the piecewise-linear nonlinearity is replaced by an appropriate cubic nonlinearity which yields very similar behavior. They found that the lowest value of the total order to have chaos in this system is 2.7.

In this paper, we intend to obtain the approximate solution of the fractional-order Chua's system and Chua-Hartley’s system via the multistep generalized differential transform method (MSGDTM). This method is only a simple modification of the generalized differential transform method (GDTM) [18–21], in which it is treated as an algorithm in a sequence of small intervals (i.e., time step) for finding accurate approximate solutions to the corresponding systems. The approximate solutions obtained by using GDTM are valid only for a short time. The ones obtained by using the MSGDTM [22–24] are more valid and accurate during a long time and are in good agreement with the classical Runge-Kutta method numerical solution when the order of the derivative is one.

This paper is organized as follows. In Section 2, we present some necessary definitions and notations related to fractional calculus. In Section 3, the proposed method is described. In Section 4, the method is applied to the fractional-order Chua’s system and the Chua-Hartley’s system. In Section 5, the numerical simulations are presented graphically. Finally, the conclusions are given in Section 6.

2. Fractional Calculus

In this section, we give some basic definitions and properties of the fractional calculus theory which are used further in this paper [1–5].

Definition 2.1. A function is said to be in the space if it can be written as for some where is continuous in , and it is said to be in the space if .

Definition 2.2. The Riemann-Liouville integral operator of order with is defined as

Properties of the operator can be found in [1–5]. We only need here the following: for , , , and , we have

where is the incomplete beta function which is defined as

The Riemann-Liouville derivative has certain disadvantages when trying to model real-world phenomena with fractional differential equations. Therefore, we will introduce a modified fractional differential operator proposed by Caputo in his work on the theory of viscoelasticity.

Definition 2.3. The Caputo fractional derivative of of order with is defined as

for , , , .

The Caputo fractional derivative was investigated by many authors, for , and , we have

For mathematical properties of fractional derivatives and integrals one can consult the mentioned references.

3. Multistep Generalized Differential Transform Method (MSGDTM)

To describe the multistep generalized differential transform method (MSGDTM) [22–24], we consider the following initial value problem for systems of fractional differential equations

subject to the initial conditions

where is the Caputo fractional derivative of order , where , for . Let be the interval over which we want to find the solution of the initial value problem ((3.1)-(3.2)). In actual applications of the generalized differential transform method (GDTM), the th-order approximate solution of the initial value problem ((3.1)-(3.2)) can be expressed by the finite series

where satisfied the recurrence relation

, and is the differential transform of function for . The basics steps of the GDTM can be found in [18–21].

Assume that the interval is divided into subintervals of equal step size by using the nodes . The main ideas of the MSGDTM are as follows.

First, we apply the GDTM to the initial value problem ((3.1)-(3.2)) over the interval , we will obtain the approximate solution ], using the initial condition , for . For and at each subinterval [,], we will use the initial condition and apply the GDTM to the initial value problem ((3.1)-(3.2)) over the interval [,]. The process is repeated and generates a sequence of approximate solutions , , for . Finally, the MSGDTM assumes the following solution:

The new algorithm, MSGDTM, is simple for computational performance for all values of . As we will see in the next section, the main advantage of the new algorithm is that the obtained solution converges for wide time regions.

4. Solving the Fractional-Order Chua’s System and Chua-Hartley’s System Using the MSGDTM Algorithm

In order to demonstrate the performance and efficiency of the multistep generalized differential transform method for solving linear and nonlinear fractional-order equations, we have applied the method to two examples. In the first example, we consider the fractional-order Chua’s chaotic system, while in the second example, we consider the fractional-order Chua-Hartley’s system.

4.1. The Fractional-Order Chua’s System

The classical Chua’s oscillator is a simple electronic circuit that exhibits nonlinear dynamical phenomena such as bifurcation and chaos [16]. Now, we consider a fractional-order Chua’s system, where integer-order derivatives are replaced by fractional-order ones. This system is expressed as

where are the state variables, are constants, and are parameters describing the order of the fractional time derivatives in the Caputo sense.

Applying the MSGDTM Algorithm to (4.1)–(4.3) gives

where , , and are the differential transformation of , , and , respectively, and equals 1 when and equals 0 otherwise. The differential transform of the initial conditions is given by , , . In view of the differential inverse transform, the differential transform series solution for the system (4.1)–(4.3) can be obtained as

According to the multistep generalized differential transform method, the series solution for the system (4.1)–(4.3) is suggested by

where , , and for satisfy the following recurrence relations:

such that , , and .

Finally, we start with , , and , using the recurrence relation given in (4.8), then we can obtain the multistep solution given in (4.7).

4.2. The Fractional-Order Chua-Hartley’s System

The Chua-Hartley’s system is different from the usual Chua’s system in that the piecewise-linear nonlinearity is replaced by an appropriate cubic nonlinearity which yields a very similar behavior. Derivatives on the left side of the differential equations are also replaced by the fractional derivatives as follows [16, 17]:

where are the state variables, is a positive constant and are parameters describing the order of the fractional time-derivatives in the Caputo sense.

Following the same procedure as the previous system and applying MSGDTM algorithm to (4.9)–(4.11) yields

where , , and are the differential transformation of , , and respectively. The differential transform of the initial conditions is given by , , and . In view of the differential inverse transform, the differential transform series solution for the system (4.9)–(4.11) can be obtained as

Now, according to the MSGDTM algorithm, the series solution for the system (4.8)–(4.11) is suggested by

where ,  , and for satisfy the following recurrence relations:

such that , , and .

Starting with , , and using the recurrence relation given in (4.15), then we can obtain the multistep solution given in (4.14).

5. Numerical Results

The MSGDTM is coded in the computer algebra package Mathemtica. The Mathemtica environment variable digits controlling the number of significant digits is set to 20 in all the calculations done in this paper. The time range studied in this work is and the step size . We take the initial conditions for Chua’s system and Chua-Hartley’s system: and .

We consider the case which corresponds to the classical Chua's system. Figure 1 represents the phase portrait for chaotic solutions. The effective dimension of (4.1)–(4.3) is defined as the sum of orders . We can see that the chaotic attractors of the fractional-order system are similar to that of the integer-order Chua's attractor as shown in Figure 2. Figure 2(b) shows the lowest order we found to yield chaos in this system is 2.62. From the numerical results in Figure 2, it is clear that the approximate solutions depend continuously on the time-fractional derivative .

(a)

(b)

Figure 1: Phase plot of Chua’s chaotic attractor with , : (a) in the - space, (b) in the -- space.

(a)

(b)

Figure 2: Phase plot of Chua’s chaotic attractor in the -- space: (a) , , , (b) , , .

Simulations were performed for the classical integer-order Chua-Hartley’s system in Figure 3. Figure 4(b) shows the lowest order we found to yield chaos in the fractional-order Chua-Hartley’s system is 2.65. From the graphical results in Figures 1–4, it is to conclude that the approximate solutions obtained using the multistep generalized differential transform method are in good agreement with the approximate solutions obtained in [16, 17].

(a)

(b)

Figure 3: Phase plot of Chua-Hartley’s chaotic attractor with and : (a) in the - space, (b) in the -- space.

(a)

(b)

Figure 4: Phase plot of Chua-Hartley's chaotic attractor in the -- space: (a) , , , (b) , , , .

6. Conclusions

In this present work, the multistep generalized differential transform method was introduced to obtain the solutions of the fractional-order Chua's system and Chua-Hartley's system by time discretization. This method has the advantage of giving an analytical form of the solution within each time interval which is not possible using purely numerical techniques like the fourth-order Runge-Kutta method (RK4). We conclude that MSGDTM is a very reliable method in solving a broad array of dynamical problems in fractional calculus due to its consistency used in a longer time frame.

References

J. Sabatier, O. P. Agrawal, and J. A. Tenreiro Machado, Advances in Fractional Calculus; Theoretical Developments and Applications in Physics and Engineering, Springer, Dordrecht, The Netherlands, 2007.

K. B. Oldham and J. Spanier, The Fractional Calculus, Academic Press, New York, NY, USA, 1974.

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

I. Podlubny, Fractional Differential Equations, vol. 198, Academic Press, San Diego, Calif, USA, 1999.

C. Li and W. Deng, “Remarks on fractional derivatives,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 777–784, 2007.

G. M. Zaslavsky, “Chaos, fractional kinetics, and anomalous transport,” Physics Reports, vol. 371, no. 6, pp. 461–580, 2002.

D. Baleanu, K. Diethelm, E. Scalas, and J. J. Trujillo, Fractional Calculus Models and Numerical Methods, World Scientific Publishing, Singapore, 2009.

R. L. Bagley and P. J. Torvik, “A theoretical basis for the application of fractional calculus,” Journal of Rheology, vol. 27, pp. 201–210, 1983.

R. L. Bagley and R. A. Calico, “Fractional order state equations for the control of viscoelastically damped structures,” Journal of Guidance, Control, and Dynamics, vol. 14, no. 2, pp. 304–311, 1991.

E. J. S. Pires, J. A. T. Machado, and P. B. de Moura, “Fractional order dynamics in a GA planner,” Signal Processing, vol. 83, pp. 2377–2386, 2003.

K. S. Hedrih and V. N. Stanojević, “A model of gear transmission: fractional order system dynamics,” Mathematical Problems in Engineering, vol. 2010, Article ID 972873, 2010.

J. Cao, C. Ma, H. Xie, and Z. Jiang, “Nonlinear dynamics of duffing system with fractional order damping,” Journal of Computational and Nonlinear Dynamics, vol. 5, no. 4, pp. 041012–041018, 2010.

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.

K. Moaddy, S. Momani, and I. Hashim, “The non-standard finite difference scheme for linear fractional PDEs in fluid mechanics,” Computers & Mathematics with Applications, vol. 61, no. 4, pp. 1209–1216, 2011.

A. G. Radwan, K. Moaddy, and S. Momani, “Stability and non-standard finite difference method of the generalized Chua's circuit,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 961–970, 2011.

I. Petrás, “A note on the fractional-order Chua's system,” Chaos Solution & Fractals, vol. 38, pp. 140–147, 2008.

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

Z. Odibat, S. Momani, and V. S. Erturk, “Generalized differential transform method: application to differential equations of fractional order,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 467–477, 2008.

S. Momani and Z. Odibat, “A novel method for nonlinear fractional partial differential equations: combination of DTM and generalized Taylor's formula,” Journal of Computational and Applied Mathematics, vol. 220, no. 1-2, pp. 85–95, 2008.

Z. Odibat and S. Momani, “A generalized differential transform method for linear partial differential equations of fractional order,” Applied Mathematics Letters, vol. 21, no. 2, pp. 194–199, 2008.

V. S. Erturk, S. Momani, and Z. Odibat, “Application of generalized differential transform method to multi-order fractional differential equations,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 8, pp. 1642–1654, 2008.

Z. M. Odibat, C. Bertelle, M. A. Aziz-Alaoui, and G. H. E. Duchamp, “A multi-step differential transform method and application to non-chaotic or chaotic systems,” Computers & Mathematics with Applications, vol. 59, no. 4, pp. 1462–1472, 2010.

V. S. Ertürk, Z. M. Odibat, and S. Momani, “An approximate solution of a fractional order differential equation model of human T-cell lymphotropic virus I (HTLV-I) infection of T-cells,” Computers & Mathematics with Applications, vol. 62, no. 3, pp. 996–1002, 2011.

V. Erturk, G. Zaman, and S. Momani, “A numericanalytic method for approximating a giving up smoking model containing fractional derivatives,” Computers & Mathematics With Applications, vol. 64, no. 10, pp. 3065–3074, 2012.