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 [25], chemical physics, signal processing, systems identification, control and robotics [6, 7], and many other areas [815].

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) [1821], 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 [2224] 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 [15].

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 [15]. 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) [2224], 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 [1821].

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 .

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 14, 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].

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.