Journal of Applied Mathematics

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

doi:10.1155/2012/258948

Research Article

Application of Piecewise Successive Linearization Method for the Solutions of the Chen Chaotic System

S. S. Motsa,¹ Y. Khan,² and S. Shateyi³

¹School of Mathematical Sciences, University of KwaZulu-Natal, Private Bag X01, Pietermaritzburg, Scottsville 3209, South Africa

²Department of Mathematics, Zhejiang University, Hangzhou 310027, China

³Department of Mathematics, University of Venda, Private Bag X5050, Thohoyandou 0950, South Africa

Received 18 April 2012; Revised 17 June 2012; Accepted 18 June 2012

Academic Editor: Hak-Keung Lam

Copyright © 2012 S. S. Motsa 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.

This paper centres on the application of the new piecewise successive linearization method (PSLM) in solving the chaotic and nonchaotic Chen system. Numerical simulations are presented graphically and comparison is made between the PSLM and Runge-Kutta-based methods. The work shows that the proposed method provides good accuracy and can be easily extended to other dynamical systems including those that are chaotic in nature.

1. Introduction

In this paper, we consider the Chen dynamical system [1], which is a three-dimensional system of ordinary differential equations with quadratic nonlinearities, defined as

where , , are dependent variables and , , are positive constants. The Chen system (1.1) does not admit a closed form solution and it is known to exhibit both chaotic and nonchaotic behaviors for distinct parameter values. Bifurcation studies [2–4] have shown that with the parameters and , system (1.1) exhibits nonchaotic behavior and chaotic behavior when and , respectively. Finding accurate and efficient methods for solving chaotic systems such as the Chen system has been an active research area for many researchers who use such systems as a benchmark for testing the accuracy and performance of methods of solutions under their development. Other variants of the Chen system such as fractional order and four-dimensional Chen systems have been the subject of study by many scholars.

The solution of nonchaotic, chaotic and hyperchaotic Chen system equations has been obtained by different analytical and numerical methods such as Runge-Kutta-based numerical schemes. In recent years, several analytical and seminumerical methods have been utilized in solving the Chen system equations [5, 6]. Because the solutions for the Chen systems rapidly oscillate within small intervals, conventional analytical and numerical algorithms have been implemented in a sequence of small multiple intervals in order to ensure accuracy convergence of the solution techniques. Recent multistage methods that have been reported for solving Chen systems include the multistage differential transform method [7, 8], the multistage homotopy analysis method [9], the multistep homotopy perturbation method [10] and the multistage variational iteration method [11].

The aim of this paper is to obtain the solution of the Chen system using the successive linearization method (SLM) that is implemented in a sequence of intervals. This new algorithm is called the piecewise successive linearization method (PSLM). The standard SLM was developed primarily for solving nonlinear boundary value problems (see, e.g., [12–17]). To the best of our knowledge, the method has not been extended to initial value problem (IVPs) and in particular chaotic IVPs of the type (1.1). We remark that the standard SLM is not suitable for solving dynamical systems of the type (1.1), hence the need to modify the method and implement it on multiple intervals. The paper focuses on the use of the proposed PSLM in solving the Chen system (1.1) using selected parameters that yield both chaotic and nonchaotic results. Numerical results are presented graphically and are compared with Runge-Kutta-based results and excellent agreement is observed.

The remainder of the paper is organized as follows. In Section 2, the basic description of the standard SLM is briefly presented. In Section 3, the PSLM is introduced. The application of the proposed method is numerically investigated in Section 4 to generate results which are discussed therein. Finally, conclusions in Section 5 close the paper.

2. Successive Linearization Method (SLM)

In this section, we give a brief description of the basic idea behind the successive linearization method (SLM) [13–15]. This will be followed by a description of the multistage extension of the SLM algorithm which is applied to the Chen system (1.1).

The SLM approach assumes that the solution of (1.1) can be expressed as

Substituting (2.1) in the governing equation (1.1) and neglecting nonlinear terms in we obtain

subject to the initial conditions

where

Starting from a suitable initial approximation , , , the solutions for , , (for ) can be obtained by iteratively solving the linear system (2.2)–(2.4). A suitable initial guess to start off the SLM algorithm (2.2)–(2.4) is obtained by solving the linear part of (1.1) subject to the problem's initial condition where are given initial values. For this system of first-order equations the initial approximation can easily be solved analytically in any given interval . However, the solution of the SLM algorithm (2.2)–(2.4) cannot be solved exactly. Spectral collocation methods (or any other numerical method) can be used to solve (2.2)–(2.4) in any given interval . After obtaining solutions for (2.3), the approximate solution , , and is determined as the series solution

Similar expressions are obtained for and . A SLM solution is said to be of order if the above series is truncated at , that is, if

3. Piecewise Successive Linearization Method

It was observed through numerical experimentation that, in its standard form, the SLM will not give accurate solutions for initial value problems of the Chen system type. Thus, it can be remarked that the SLM is ideally suited for boundary value problems whose solutions do not rapidly change in behaviour or oscillate over small regions of the domain of the governing problem. The SLM solution can thus be considered to be local in nature and may not be suitable for initial value problems at very large values of the independent variable . A simple way of ensuring the validity of the approximations for all values of is to determine the solution in a sequence of equal intervals, which are subject to continuity conditions at the end points of each interval. To extend this solution over the interval , we divide the interval into sub-intervals , where . We solve (2.2)–(2.4) in each subinterval . When applied to the multiple subintervals the SLM will be referred to as the piecewise SLM (PSLM).

Let , , be the solution of (2.2)–(2.4) in the first subinterval and , , be the solutions in the subintervals for . The initial conditions used in obtaining the solutions in the subinterval are obtained from the initial conditions of the subinterval . Thus, we solve

subject to the initial conditions

The initial approximations for solving (3.1) are obtained as solutions of the system

subject to the initial condition

An exact solution of (3.3) can easily be found. To solve (3.1) we use the Chebyshev spectral collocation method on each interval . Before applying the spectral method the region is transformed to the interval on which the spectral collocation method is defined. This can be achieved by using the linear transformation

in each interval (for ), After the transformation, the interval is discretized using the Chebyshev-Gauss-Lobatto collocation points [18, 19]

which are the extrema of the th-order Chebyshev polynomial

The basic idea behind the pseudo-spectral collocation method is the introduction of a differentiation matrix which is used to approximate the derivatives of the unknown variables , , at the collocation points as the matrix vector product

where and is the vector function at the collocation points . Similar expressions can be obtained for and . The Chebyshev derivative matrix is of size and its entries are defined [18, 19] as

with

Applying the Chebyshev spectral collocation method in (3.1) gives

where

and is an identity matrix, is an matrix of zeros. The denotes an diagonal matrix corresponding to that is evaluated at the collocation points and is an vector corresponding to the quantities that have been evaluated at the collocation points. Thus, starting from the initial approximation obtained as a solution of (3.3)-(3.4), the recurrence formula

can be used to obtain the solutions , , in the interval . The solution approximating , for example, in the entire interval is given by

Similar expressions can be obtained for and . We observe that when , the proposed piecewise successive linearization method (PSLM) becomes equivalent to the original SLM algorithm.

4. Results and Discussion

In this section, numerical simulations are given to validate the proposed piecewise successive linearization method. The PSLM results were obtained using collocation points and ten iterations (that is ) in each interval. Unless otherwise specified, the width of each interval was taken to be . We also fix the values of the parameters and with for the nonchaotic case and for the chaotic case. The initial conditions are , , and . The results for all the simulations done in this paper are computed for the time span . In order to assess the accuracy and performance of the proposed PSLM approach, the present results are compared with those obtained with the MATLAB built-in solver . The solver integrates a system of ordinary differential equations using explicit 4th and 5th Runge-Kutta (4,5) formula, the Dormand-Prince pair [20].

The results of the PSLM simulation of the Chen system for the nonchaotic case are shown in Figures 1 and 2. In Figure 1 the results of the comparison between the PSLM and Runge-Kutta results for the time series evolution of the state variables are given. It can be seen that there is good agreement between the two results. Figure 2 depicts the phase portraits for the nonchaotic case. We observe that the system converges to a periodic orbit and results in a limit cycle as shown in Figure 2 when the parameters , , and are used.

(a)

(b)

(c)

Figure 1: Time series PSLM (solid line) and (dots) results for nonchaotic case with parameters , , and .

(a)

(b)

(c)

Figure 2: Phase portraits of the Chen system using the PSLM for the nonchaotic case.

The results of the PSLM simulation of the Chen system for the chaotic case are shown in Figures 3 and 4. In Figure 3, the PSLM results are compared with the Runge-Kutta results when parameters , , and are used. For this set of parameters, the system becomes chaotic and the chaotic attractors are shown in Figure 4. We observe that there is good agreement between the PSLM and Runge-Kutta results in the time series illustration depicted in Figure 3. This shows that the proposed method of solution has great potential of being a useful tool for solving complex dynamical systems with chaotic and nonchaotic behaviour.

(a)

(b)

(c)

Figure 3: Time series PSLM (solid line) and (dots) results for nonchaotic case with parameters , , and .

(a)

(b)

(c)

Figure 4: Phase portraits of the Chen system using the PSLM for the chaotic case.

5. Conclusion

In this paper, we carefully develop a novel piecewise successive linearization method (PSLM) which seeks to overcome the difficulties arising when the successive linearization method is extended to initial value problems with complex behaviour. The proposed method of solution is implemented on the Chen system and results are generated using known parameters that give both chaotic and nonchaotic outcomes. The validity of PSLM is verified by comparing its results with the results of built in Matlab-based numerical Runge-Kutta method, . Close agreement of the two sets of results is observed thus demonstrating the accuracy of the PSLM approach for the particular problem considered. In conclusion, the PSLM may be considered as a robust refinement and extension of the existing SLM. The extended method may also be valid for other complex nonlinear initial value problems with chaotic behaviour. This paper can be used as a standard example for other applications in engineering and applied sciences.

References

G. Chen and T. Ueta, “Yet another chaotic attractor,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 9, no. 7, pp. 1465–1466, 1999.

C. P. Li and G. R. Chen, “A note one Hopf bifurcation in Chen’s system,” International Journal of Bifurcation and Chaos, vol. 13, no. 6, pp. 1609–1615, 2003.

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.

T. Ueta and G. Chen, “Bifurcation analysis of Chen's equation,” International Journal of Bifurcation and Chaos, vol. 10, no. 8, pp. 1917–1931, 2000.

A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Homotopy approach for the hyperchaotic Chen system,” Physica Scripta, vol. 81, Article ID 045005, 2010.

M. S. M. Noorani, I. Hashim, R. Ahmad, S. A. Bakar, E. S. Ismail, and A. M. Zakaria, “Comparing numerical methods for the solutions of the Chen system,” Chaos, Solitons & Fractals, vol. 32, p. 1296, 2007.

A. K. Alomari, “A new analytic solution for fractional chaotic dynamical systems using the differential transform method,” Computers & Mathematics with Applications, vol. 61, no. 9, pp. 2528–2534, 2011.

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.

A. K. Alomari, M. S. M. Noorani, and R. Nazar, “Adaptation of homotopy analysis method for the numericanalytic solution of Chen system,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, p. 2336, 2009.

M. S. H. Chowdhury and I. Hashim, “Application of multistage homotopy-perturbation method for the solutions of the Chen system,” Nonlinear Analysis, vol. 10, no. 1, pp. 381–391, 2009.

S. M. Goh, M. S. M. Noorani, and I. Hashim, “On solving the chaotic Chen system: a new time marching design for the variational iteration method using Adomian's polynomial,” Numerical Algorithms, vol. 54, no. 2, pp. 245–260, 2010.

F. G. Awad, P. Sibanda, S. S. Motsa, and O. D. Makinde, “Convection from an inverted cone in a porous medium with cross-diffusion effects,” Computers & Mathematics with Applications, vol. 61, no. 5, pp. 1431–1441, 2011.

Z. G. Makukula, P. Sibanda, and S. S. Motsa, “A note on the solution of the von Kármán equations using series and Chebyshev spectral methods,” Boundary Value Problems, vol. 2010, Article ID 471793, 17 pages, 2010.

S. S. Motsa, “New algorithm for solving non-linear BVPs in heat transfer,” International Journal of Modeling, Simulation & Scientific Computing, vol. 2, no. 3, pp. 355–373, 2011.

S. S. Motsa, P. Sibanda, and S. Shateyi, “On a new quasi-linearization method for systems of nonlinear boundary value problems,” Mathematical Methods in the Applied Sciences, vol. 34, no. 11, pp. 1406–1413, 2011.

S. S. Motsa and P. Sibanda, “A linearisation method for non-linear singular boundary value problems,” Computers & Mathematics with Applications, vol. 63, no. 7, pp. 1197–1203, 2012.

S. Shateyi and S. S. Motsa, “Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with Hall effect,” Boundary Value Problems, vol. 2010, Article ID 257568, 20 pages, 2010.

C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods in Fluid Dynamics, Springer Series in Computational Physics, Springer, Berlin, Germany, 1988.

L. N. Trefethen, Spectral Methods in MATLAB, vol. 10 of Software, Environments, and Tools, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, Pa, USA, 2000.

J. R. Dormand and P. J. Prince, “A family of embedded Runge-Kutta formulae,” Journal of Computational and Applied Mathematics, vol. 6, no. 1, pp. 19–26, 1980.