JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 258948 10.1155/2012/258948 258948 Research Article Application of Piecewise Successive Linearization Method for the Solutions of the Chen Chaotic System Motsa S. S. 1 Khan Y. 2 Shateyi S. 3 Lam Hak-Keung 1 School of Mathematical Sciences University of KwaZulu-Natal Private Bag X01 Pietermaritzburg Scottsville 3209 South Africa ukzn.ac.za 2 Department of Mathematics, Zhejiang University Hangzhou 310027 China zju.edu.cn 3 Department of Mathematics University of Venda Private Bag X5050 Thohoyandou 0950 South Africa univen.ac.za 2012 15 7 2012 2012 18 04 2012 17 06 2012 18 06 2012 2012 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.

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 , which is a three-dimensional system of ordinary differential equations with quadratic nonlinearities, defined as (1.1)dxdt=a(y-x),dydt=(c-a)x-xz+cy,dzdt=xy-bz, where x, y, z are dependent variables and a, b, c 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  have shown that with the parameters a=35 and c=28, system (1.1) exhibits nonchaotic behavior and chaotic behavior when b=12 and b=3, 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 , the multistep homotopy perturbation method  and the multistage variational iteration method .

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., ). 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) . 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 (2.1)x(t)=xs(t)+m=0s-1xm(t),y(t)=ys(t)+m=0s-1ym(t),z(t)=zs(t)+m=0s-1zm(t). Substituting (2.1) in the governing equation (1.1) and neglecting nonlinear terms in xs,ys,zs we obtain (2.2)dxsdt+axs-ays=r1,s-1,(2.3)dysdt-cys+α1,s-1xs+α2,s-1zs=r2,s-1,(2.4)dzsdt+bzs-β1,s-1xs-β2,s-1ys=r3,s-1, subject to the initial conditions (2.5)xs(0)=ys(0)=zs(0)=0, where (2.6)α1,s-1=a-c+m=0s-1zm,α2,s-1=m=0s-1xm,β1,s-1=m=0s-1ym,β2,s-1=m=0s-1xm,r1,s-1=-m=0s-1(dxmdt-aym+axm),r2,s-1=-m=0s-1(dymdt-(c-a)xm-cym)+(m=0s-1xm)(m=0s-1zm),r3,s-1=-m=0s-1(dzmdt+bzm)+(m=0s-1xm)(m=0s-1ym). Starting from a suitable initial approximation x0(t), y0(t), z0(t), the solutions for xs(t), ys(t), zs(t) (for s=1,2,3,) 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 x(0)=x0,y(0)=y0,z(0)=z0 where x0,y0,z0 are given initial values. For this system of first-order equations the initial approximation can easily be solved analytically in any given interval t[t0,tn]. 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 t[t0,tn]. After obtaining solutions for (2.3), the approximate solution x(t), y(t), and z(t) is determined as the series solution (2.7)x(t)=x0(t)+x1(t)+x2(t)+. Similar expressions are obtained for y(t) and z(t). A SLM solution is said to be of order M if the above series is truncated at s=M, that is, if (2.8)x(t)=m=0Mxm(t).

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 t. A simple way of ensuring the validity of the approximations for all values of t 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 Ω=[t0,tn], we divide the interval Ω into sub-intervals Ωi=[ti-1,ti], i=1,2,3,,n where t0t1tn. We solve (2.2)–(2.4) in each subinterval Ωi. When applied to the multiple subintervals the SLM will be referred to as the piecewise SLM (PSLM).

Let x1(t), y1(t), z1(t) be the solution of (2.2)–(2.4) in the first subinterval [t0,t1] and xi(t), yi(t), zi(t) be the solutions in the subintervals Ωi for 2in. The initial conditions used in obtaining the solutions in the subinterval Ωi(2in) are obtained from the initial conditions of the subinterval Ωi-1. Thus, we solve (3.1)dxsidt+axsi-aysi=r1,i-1i,dysidt-cysi+α1,s-1ixsi+α2,s-1izsi=r2,s-1i,dzsidt+bzsi-β1,s-1ixsi-β2,s-1iysi=r3,s-1i, subject to the initial conditions (3.2)xsi(ti-1)=ysi(ti-1)=zsi(ti-1)=0. The initial approximations for solving (3.1) are obtained as solutions of the system (3.3)dx0idt=a(y0i-x0i),dy0idt=(c-a)x0i+cy0i,dz0idt=-bz0i, subject to the initial condition (3.4)x0i(ti-1)=xi-1(ti-1),y0i(ti-1)=yi-1(ti-1),z0i(ti-1)=zi-1(ti-1). An exact solution of (3.3) can easily be found. To solve (3.1) we use the Chebyshev spectral collocation method on each interval [ti-1,ti]. Before applying the spectral method the region [ti-1,ti] is transformed to the interval [-1,1] on which the spectral collocation method is defined. This can be achieved by using the linear transformation (3.5)t=(ti-ti-1)τ2+(ti+ti-1)2,τ[-1,1] in each interval [ti-1,ti] (for i=1,,n), After the transformation, the interval [ti-1,ti] is discretized using the Chebyshev-Gauss-Lobatto collocation points [18, 19] (3.6)τji=cos(πjN),j=0,1,,N which are the extrema of the Nth-order Chebyshev polynomial (3.7)TN(τ)=cos(Ncos-1τ). The basic idea behind the pseudo-spectral collocation method is the introduction of a differentiation matrix D which is used to approximate the derivatives of the unknown variables xsi(t), ysi(t), zsi(t) at the collocation points as the matrix vector product (3.8)dxsidt=k=0NDjkxsi(τki)=DXsi,j=0,1,,N, where D=2D/(ti-ti-1) and Xsi=[xsi(τ0i),xsi(τ2i),,xsi(τNi)]T is the vector function at the collocation points τji. Similar expressions can be obtained for ysi(t) and zsi(t). The Chebyshev derivative matrix D is of size (N+1)×(N+1) and its entries are defined [18, 19] as (3.9)Djk=cjck(-1)j+kzj-zk,jk;j,k=0,1,,N,Dkk=-zk2(1-zk2),k=1,2,,N-1,D00=2N2+16=-DNN, with (3.10)ck={2k=0,N1-1kN-1. Applying the Chebyshev spectral collocation method in (3.1) gives (3.11)As-1iFsi=Rs-1i,Fsi(τNi-1)=0, where (3.12)As-1i=[D+aI-aIOα~1,s-1iD-cIα~2,s-1i-β~1,s-1i-β~2,s-1iD+bI],Rs-1i=[r1,s-1ir2,s-1ir3,s-1i],Fsi=[XsiYsiZsi] and I is an (N+1)×(N+1) identity matrix, O is an (N+1)×(N+1) matrix of zeros. The ()~ denotes an (N+1)×(N+1) diagonal matrix corresponding to () that is evaluated at the collocation points and rk,s-1i(k=1,2,3) is an (N+1)×1 vector corresponding to the quantities rk,s-1i(t) 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 (3.13)Fsi=(As-1i)-1Rs-1i can be used to obtain the solutions xi(t), yi(t), zi(t) in the interval [ti,ti-1]. The solution approximating x(t), for example, in the entire interval [t0,tn] is given by (3.14)x(t)={x1(t),t[t0,t1]x2(t),t[t1,t2]xn(t),t[tn-1,tn]. Similar expressions can be obtained for y(t) and z(t). We observe that when n=1, 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 N=10 collocation points and ten iterations (that is M=10) in each [ti-1,ti] interval. Unless otherwise specified, the width of each interval Δt=ti-ti-1 was taken to be Δt=0.1. We also fix the values of the parameters a=35 and c=28 with b=12 for the nonchaotic case and b=3 for the chaotic case. The initial conditions are x(0)=-10, y(0)=0, and z(0)=37. The results for all the simulations done in this paper are computed for the time span t[0,10]. 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 ode45. The ode45 solver integrates a system of ordinary differential equations using explicit 4th and 5th Runge-Kutta (4,5) formula, the Dormand-Prince pair .

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 a=35, b=12, and c=28 are used.

Time series PSLM (solid line) and ode45 (dots) results for nonchaotic case with parameters a=35, b=12, and c=28.

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 a=35, c=28, and b=3 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.

Time series PSLM (solid line) and ode45 (dots) results for nonchaotic case with parameters a=35, b=3, and c=28.

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, ode45. 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.

Chen G. Ueta T. Yet another chaotic attractor International Journal of Bifurcation and Chaos in Applied Sciences and Engineering 1999 9 7 1465 1466 10.1142/S0218127499001024 1729683 ZBL0962.37013 Li C. P. Chen G. R. A note one Hopf bifurcation in Chen’s system International Journal of Bifurcation and Chaos 2003 13 6 1609 1615 10.1142/S0218127403007394 J. Zhou T. Chen G. Zhang S. Local bifurcations of the Chen system International Journal of Bifurcation and Chaos 2002 12 10 2257 2270 10.1142/S0218127402005819 1941281 ZBL1047.34044 Ueta T. Chen G. Bifurcation analysis of Chen's equation International Journal of Bifurcation and Chaos 2000 10 8 1917 1931 10.1142/S0218127400001183 1787214 ZBL1090.37531 Alomari A. K. Noorani M. S. M. Nazar R. Homotopy approach for the hyperchaotic Chen system Physica Scripta 2010 81 045005 ZBL1190.37029 Noorani M. S. M. Hashim I. Ahmad R. Bakar S. A. Ismail E. S. Zakaria A. M. Comparing numerical methods for the solutions of the Chen system Chaos, Solitons & Fractals 2007 32 1296 ZBL1131.65101 Alomari A. K. A new analytic solution for fractional chaotic dynamical systems using the differential transform method Computers & Mathematics with Applications 2011 61 9 2528 2534 10.1016/j.camwa.2011.02.043 2795001 ZBL1221.65191 Odibat Z. M. Bertelle C. Aziz-Alaoui M. A. Duchamp G. H. E. A multi-step differential transform method and application to non-chaotic or chaotic systems Computers & Mathematics with Applications 2010 59 4 1462 1472 10.1016/j.camwa.2009.11.005 2591936 ZBL1189.65170 Alomari A. K. Noorani M. S. M. Nazar R. Adaptation of homotopy analysis method for the numericanalytic solution of Chen system Communications in Nonlinear Science and Numerical Simulation 2009 14 2336 Chowdhury M. S. H. Hashim I. Application of multistage homotopy-perturbation method for the solutions of the Chen system Nonlinear Analysis 2009 10 1 381 391 10.1016/j.nonrwa.2007.09.014 2451718 ZBL1154.65350 Goh S. M. Noorani M. S. M. Hashim I. On solving the chaotic Chen system: a new time marching design for the variational iteration method using Adomian's polynomial Numerical Algorithms 2010 54 2 245 260 10.1007/s11075-009-9333-9 2646079 ZBL1190.65189 Awad F. G. Sibanda P. Motsa S. S. Makinde O. D. Convection from an inverted cone in a porous medium with cross-diffusion effects Computers & Mathematics with Applications 2011 61 5 1431 1441 10.1016/j.camwa.2011.01.015 2773416 ZBL1217.76071 Makukula Z. G. Sibanda P. Motsa S. S. A note on the solution of the von Kármán equations using series and Chebyshev spectral methods Boundary Value Problems 2010 2010 17 471793 10.1155/2010/471793 2728258 ZBL1207.35248 Motsa S. S. New algorithm for solving non-linear BVPs in heat transfer International Journal of Modeling, Simulation & Scientific Computing 2011 2 3 355 373 10.1142/S1793962311000499 Motsa S. S. Sibanda P. Shateyi S. On a new quasi-linearization method for systems of nonlinear boundary value problems Mathematical Methods in the Applied Sciences 2011 34 11 1406 1413 10.1002/mma.1449 2839382 Motsa S. S. Sibanda P. A linearisation method for non-linear singular boundary value problems Computers & Mathematics with Applications 2012 63 7 1197 1203 10.1016/j.camwa.2011.12.035 Shateyi S. Motsa S. S. Variable viscosity on magnetohydrodynamic fluid flow and heat transfer over an unsteady stretching surface with Hall effect Boundary Value Problems 2010 2010 20 10.1155/2010/257568 257568 2719288 ZBL1195.76438 Canuto C. Hussaini M. Y. Quarteroni A. Zang T. A. Spectral Methods in Fluid Dynamics 1988 Berlin, Germany Springer xiv+557 Springer Series in Computational Physics 917480 Trefethen L. N. Spectral Methods in MATLAB 2000 10 Philadelphia, Pa, USA Society for Industrial and Applied Mathematics (SIAM) xviii+165 Software, Environments, and Tools 10.1137/1.9780898719598 1776072 Dormand J. R. Prince P. J. A family of embedded Runge-Kutta formulae Journal of Computational and Applied Mathematics 1980 6 1 19 26 10.1016/0771-050X(80)90013-3 568599 ZBL0448.65045