The Scientific World Journal
Volume 2014 (2014), Article ID 943293, 10 pages
http://dx.doi.org/10.1155/2014/943293

Research Article

Multistage Spectral Relaxation Method for Solving the Hyperchaotic Complex Systems

Hassan Saberi Nik¹ and Paulo Rebelo²

¹Department of Mathematics, Islamic Azad University, Mashhad Branch, Mashhad, Iran
²Departamento de Matemática, Universidade da Beira Interior, 6201-001 Covilhã, Portugal

Received 4 August 2014; Accepted 12 September 2014; Published 16 October 2014

Academic Editor: Fazlollah Soleymani

Copyright © 2014 Hassan Saberi Nik and Paulo Rebelo. 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 present a pseudospectral method application for solving the hyperchaotic complex systems. The proposed method, called the multistage spectral relaxation method (MSRM) is based on a technique of extending Gauss-Seidel type relaxation ideas to systems of nonlinear differential equations and using the Chebyshev pseudospectral methods to solve the resulting system on a sequence of multiple intervals. In this new application, the MSRM is used to solve famous hyperchaotic complex systems such as hyperchaotic complex Lorenz system and the complex permanent magnet synchronous motor. We compare this approach to the Runge-Kutta based ode45 solver to show that the MSRM gives accurate results.

1. Introduction

Chaos theory studies the behaviour of dynamical systems that are highly sensitive to initial conditions and have complex and highly unpredictable profiles [1, 2]. Chaotic systems can be observed in a wide variety of applications. In 1982, the complex Lorenz equations were proposed by Fowler et al. [3], which extended nonlinear systems into complex space. After that, some research works in this field have been achieved [4–9]. With in-depth study of complex nonlinear systems, a variety of physical phenomena could be described by the chaotic or hyperchaotic complex systems, for instance, the detuned laser systems and the amplitudes of electromagnetic fields.

The nature of complex chaotic systems precludes the possibility of obtaining closed form analytical solutions of the underlying governing equations. Thus, approximate-analytical methods, which are implemented on a sequence of multiple intervals to increase their radius of convergence, are often used to solve IVPs modelling chaotic systems. Examples of multistage methods that have been developed recently to solve IVPs for chaotic and nonchaotic systems include the multistage homotopy analysis method [10], piecewise homotopy perturbation methods [11, 12], multistage variational iteration method [13], and multistage differential transformation method [14]. Other multistage methods which use numerical integration techniques have also been proposed such as the piecewise spectral homotopy analysis method [15–17] which uses a spectral collocation method to perform the integration process. Accurate solutions of highly chaotic and hyperchaotic systems require resolution over many small intervals. Thus, seeking analytical solutions over the numerous intervals may be impractical or computationally expensive if the solution is sought over very long intervals.

In this paper, we propose a piecewise or multistage spectral relaxation method (MSRM) for solving the hyperchaotic complex systems as an accurate and robust alternative to recent multistage methods. The proposed MSRM was developed using the Gauss-Seidel idea of decoupling systems of equations and using Chebyshev pseudospectral methods to solve the resulting decoupled system on a sequence of multiple intervals. The spectral relaxation method (SRM) was recently proposed in [18, 19].

The rest of the paper is organized as follows. In Section 2, we give a brief description of the proposed MSRM algorithm. In Section 3, we present the numerical implementation of the MSRM on two examples of hyperchaotic complex systems. Finally, the conclusion is given in Section 4.

2. Multistage Spectral Relaxation Method

In this section, we give a brief description of the numerical method of solution used to solve the nonlinear hyperchaotic complex. We employ the multistage spectral relaxation method (MSRM) proposed in [19]. The MSRM algorithm is based on a Gauss-Seidel type of relaxation that decouples and linearises the system and the use of spectral collocation method to solve the linearised equations in a sequential manner. For compactness, we express the system of nonlinear first order differential equations in the form

subject to the initial conditions

where are the unknown variables and are the corresponding initial conditions, are known constant input parameters and is the nonlinear component of the th equation and the dot denotes differentiation with respect to time .

The scheme computes the solution of (1) in a sequence of equal subintervals that makes the entire interval. We define the interval of integration as and divide it into a sequence of nonoverlapping subintervals   , where and . We denote the solution of (1) in the first subinterval as and the solutions in the subsequent subintervals as . For obtaining the solution in the first interval , (2) is used as the initial condition. By using the continuity condition between neighbouring subintervals the obtained solution in the interval is used to obtain the initial condition for the next subinterval . This is applied over the successive subintervals; that is, the obtained solution for each subinterval is used to obtain the initial condition for the next subinterval   . Thus, in each interval we must solve

subject to

where is the Kronecker delta. As mentioned earlier, the main idea behind the MSRM scheme is decoupling the system of nonlinear IVPs using the Gauss-Seidel idea of decoupling systems of algebraic equations. The proposed MSRM iteration scheme for the solution in the interval is given as

subject to the initial conditions

where is the estimate of the solution after iterations. A suitable initial guess to start the iteration scheme (5) is one that satisfies the initial condition (6). A convenient choice of initial guess that was found to work in the numerical experiments considered in this work is

The Chebyshev spectral method is used to solve (5) on each interval . First, the region is transformed to the interval on which the spectral method is defined by using the linear transformation,

in each interval for . We then discretize the interval using the Chebyshev-Gauss-Lobatto collocation points [20]:

which are the extrema of the th order Chebyshev polynomial:

The Chebyshev spectral collocation method is based on the idea of introducing 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 are the vector functions at the collocation points .

Applying the Chebyshev spectral collocation method in (5) gives

with

where is an identity matrix of order . Thus, starting from the initial approximation (7), the recurrence formula

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

3. Numerical Examples

In this section, we consider two examples which demonstrate the efficiency and accuracy of the proposed method. In particular, we use the MSRM algorithm as an appropriate tool for solving nonlinear IVPs; we apply the method to two complex nonlinear chaotic systems.

Example 1. The hyperchaotic complex Lorenz system can be described as

where , , , , , and are the conjugates of and . When the parameters are chosen as , , , and , the system (16) is hyperchaotic [21].
Replacing the complex variables in system (16) with real and imaginary number variables, one can get an equivalent system as follows:

For (17), the parameters and are defined as

with all other and for .

Through numerical experimentation, it was determined that collocation points and 5 iterations of the MSRM scheme at each interval were sufficient to give accurate results in each interval. Tables 1 and 2 show a comparison of the solutions of the hyperchaotic complex Lorenz system computed by the MSRM and ode45. In Figures 1, 2, and 3, the MSRM graphical results are also compared with ode45 and good agreement is observed. The MRSM phase portraits in Figures 4 and 5 were also found to be exactly the same as those computed using ode45. This shows that the proposed MSRM is a valid tool for solving the hyperchaotic complex Lorenz system.

Table 1: Numerical comparison between MSRM and ode45 for the hyperchaotic complex Lorenz system.

MSRM
ode45
MSRM
ode45
MSRM
ode45

2
−2.91138
−2.91138
21.73155
21.73155
−3.24491
−3.24491

4
−3.63001
−3.63001
6.52144
6.52144
−6.30884
−6.30884

6
2.80571
2.80571
−2.77638
−2.77638
−2.37099
−2.37099

8
0.01134
0.01134
2.09585
2.09585
−0.14880
−0.14880

10
−0.80219
−0.80219
16.48559
16.48560
−0.06690
−0.06690

Table 2: Numerical comparison between MSRM and ode45 for the hyperchaotic complex Lorenz system.

MSRM
ode45
MSRM
ode45
MSRM
ode45

2
23.96851
23.96851
44.32071
44.32071
26.54682
26.54682

4
11.30830
11.30830
14.68007
14.68007
3.25221
3.25221

6
4.65208
4.65208
39.34559
39.34559
12.99055
12.99055

8
−4.99685
−4.99685
33.79560
33.79560
8.02232
8.02232

10
1.98179
1.98179
50.59739
50.59740
24.48234
24.48234

Figure 1: Comparison between the MSRM and ode45 results for the hyperchaotic complex Lorenz system.

Figure 2: Comparison between the MSRM and ode45 results for the hyperchaotic complex Lorenz system.

Figure 3: Comparison between the MSRM and ode45 results for the hyperchaotic complex Lorenz system.

Figure 4: Phase portraits of the hyperchaotic complex Lorenz system.

Figure 5: Phase portraits of the hyperchaotic complex Lorenz system.

Example 2. State equations of a permanent magnet synchronous motor system in a field-oriented rotor can be described as follows [22, 23]:

where , and are the state variables which represent currents and motor angular frequency, respectively; and are the direct-axis stator and quadrature-axis stator voltage components, respectively; is the polar moment of inertia; is the external load torque; is the viscous damping coefficient; is the stator winding resistance; and are the direct-axis stator inductors and quadrature-axis stator inductors, respectively; is the permanent magnet flux; and is the number of pole-pairs; the parameters , , , , , , are all positive.

When the air gap is even, and the motor has no load or power outage, the dimensionless equations of a permanent magnet synchronous motor system can be depicted as

where are both positive parameters. If the current in the system (19) is plural and the variables in the system (20) are complex numbers, by changing cross coupled terms and to conjugate form, Wang and Zhang got a complex permanent magnet synchronous motor system as follows [24]:

where , , , , and are the conjugates of and . Replacing the complex variables in system (21) with real and imaginary number variables, Wang and Zhang got an equivalent system as follows (see [24]):

where are positive parameters determining the chaotic behaviors and bifurcations of system (22). When the parameters satisfy , , there is one positive Lyapunov exponent, two Lyapunov exponents of zero, and two negative Lyapunov exponents for system (22), which means system (22) is chaotic [24]. The values of parameters and initial values are , , and , , , , .

For (21), the parameters and are defined as

with all other and for .

The results obtained were compared to those from the MATLAB inbuilt solver, ode45. The ode45 solver integrates a system of ordinary differential equations using explicit 4th and 5th Runge-Kutta formula. Tables 3 and 4 show a comparison of the solutions of the complex permanent magnet synchronous motor computed by the MSRM and ode45. In Figures 6, 7, and 8, the MSRM graphical results are also compared with ode45 and good agreement is observed. The MRSM phase portraits in Figures 9 and 10 were also found to be exactly the same as those computed using ode45. This shows that the proposed MSRM is a valid tool for solving the complex permanent magnet synchronous motor.

Table 3: Numerical comparison between MSRM and ode45 for the complex permanent magnet synchronous motor.

MSRM
ode45
MSRM
ode45
MSRM
ode45

3
−3.85711
−3.85711
−5.66683
−5.66683
−5.20445
−5.20445

10
−0.33729
−0.33729
−0.49554
−0.49554
−0.49104
−0.49104

17
0.12630
0.12631
0.18555
0.18557
0.15550
0.15551

24
0.05091
0.05105
0.07480
0.07501
0.19500
0.19518

31
−2.55034
−2.54878
−3.74694
−3.74465
−0.79819
−0.79326

38
−3.93154
−3.73551
−5.77619
−5.48818
−5.33693
−5.20595

Table 4: Numerical comparison between MSRM and ode45 for the complex permanent magnet synchronous motor.

MSRM
ode45
MSRM
ode45

3
−7.64635
−7.64635
15.05932
15.05932

10
−0.72144
−0.72143
10.73663
10.73663

17
0.22846
0.22848
14.25582
14.25583

24
0.28649
0.28675
19.33844
19.33921

31
−1.17270
−1.16545
25.34856
25.35739

38
−7.84098
−7.64855
14.98250
14.03140

Figure 6: Comparison between the MSRM and ode45 results for the complex permanent magnet synchronous motor.

Figure 7: Comparison between the MSRM and ode45 results for the complex permanent magnet synchronous motor.

Figure 8: Comparison between the MSRM and ode45 results for the complex permanent magnet synchronous motor.

Figure 9: Phase portraits of the complex permanent magnet synchronous motor.

Figure 10: Phase portraits of the complex permanent magnet synchronous motor.

4. Conclusion

In this paper, we have applied a spectral method called the multistage spectral relaxation method (MSRM) for the solutions of hyperchaotic complex systems. The proposed MSRM was developed using the Gauss-Seidel idea of decoupling systems of equations and using Chebyshev pseudospectral methods to solve the resulting decoupled system on a sequence of multiple intervals. The proposed MSRM was used to solve the hyperchaotic complex Lorenz system and complex permanent magnet synchronous motor. The accuracy and validity of the proposed method was tested against Matlab Runge-Kutta based inbuilt solvers and against previously published results.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

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.

Z. Wei, R. Wang, and A. Liu, “A new finding of the existence of hidden hyperchaotic attractors with no equilibria,” Mathematics and Computers in Simulation, vol. 100, pp. 13–23, 2014.

A. C. Fowler, J. D. Gibbon, and M. J. McGuinness, “The complex Lorenz equations,” Physica D: Nonlinear Phenomena, vol. 4, no. 2, pp. 139–163, 1981/82.

A. G. Vladimirov, V. Y. Toronov, and V. L. Derbov, “Properties of the phase space and bifurcations in the complex Lorenz model,” Technical Physics, vol. 43–48, pp. 877–884, 1998.

A. Rauh, L. Hannibal, and N. B. Abraham, “Global stability properties of the complex Lorenz model,” Physica D: Nonlinear Phenomena, vol. 99, no. 1, pp. 45–58, 1996.

R. L. Lang, “A stochastic complex model with random imaginary noise,” Nonlinear Dynamics, vol. 62, no. 3, pp. 561–565, 2010.

G. M. Mahmoud, M. E. Ahmed, and E. E. Mahmoud, “Analysis of hyperchaotic complex Lorenz systems,” International Journal of Modern Physics C, vol. 19, no. 10, pp. 1477–1494, 2008.

G. M. Mahmoud, E. E. Mahmoud, and M. E. Ahmed, “On the hyperchaotic complex Lü system,” Nonlinear Dynamics, vol. 58, no. 4, pp. 725–738, 2009.

G. M. Mahmoud, M. A. Al-Kashif, and S. A. Aly, “Basic properties and chaotic synchronization of complex Lorenz system,” International Journal of Modern Physics C, vol. 18, no. 2, pp. 253–265, 2007.

A. K. Alomari, M. S. Noorani, and R. Nazar, “Explicit series solutions of some linear and nonlinear Schrodinger equations via the homotopy analysis method,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1196–1207, 2009.

S. Effati, H. Saberi Nik, and M. Shirazian, “An improvement to the homotopy perturbation method for solving the Hamilton-Jacobi-Bellman equation,” IMA Journal of Mathematical Control and Information, vol. 30, no. 4, pp. 487–506, 2013.

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

S. M. Goh, M. S. 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.

A. Gökdoğan, M. Merdan, and A. Yildirim, “The modified algorithm for the differential transform method to solution of Genesio systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 1, pp. 45–51, 2012.

H. Saberi Nik, S. Effati, S. S. Motsa, and M. Shirazian, “Spectral homotopy analysis method and its convergence for solving a class of nonlinear optimal control problems,” Numerical Algorithms, vol. 65, no. 1, pp. 171–194, 2014.

H. Saberi Nik, S. Effati, S. S. Motsa, and S. Shateyi, “A new piecewise-spectral homotopy analysis method for solving chaotic systems of initial value problems,” Mathematical Problems in Engineering, vol. 2013, Article ID 583193, 13 pages, 2013.

S. Effati, H. Saberi Nik, and A. Jajarmi, “Hyperchaos control of the hyperchaotic Chen system by optimal control design,” Nonlinear Dynamics, vol. 73, no. 1-2, pp. 499–508, 2013.

S. S. Motsa, P. G. Dlamini, and M. Khumalo, “Solving hyperchaotic systems using the spectral relaxation method,” Abstract and Applied Analysis, vol. 2012, Article ID 203461, 18 pages, 2012.

S. S. Motsa, P. Dlamini, and M. Khumalo, “A new multistage spectral relaxation method for solving chaotic initial value systems,” Nonlinear Dynamics, vol. 72, no. 1-2, pp. 265–283, 2013.

L. N. Trefethen, Spectral Methods in MATLAB, SIAM, 2000.

C. Luo and X. Wang, “Hybrid modified function projective synchronization of two different dimensional complex nonlinear systems with parameters identification,” Journal of the Franklin Institute, vol. 350, no. 9, pp. 2646–2663, 2013.

J. Yu, B. Chen, H. Yu, and J. Gao, “Adaptive fuzzy tracking control for the chaotic permanent magnet synchronous motor drive system via backstepping,” Nonlinear Analysis: Real World Applications, vol. 12, no. 1, pp. 671–681, 2011.

F. Zhang, C. Mu, X. Wang, I. Ahmed, and Y. Shu, “Solution bounds of a new complex PMSM system,” Nonlinear Dynamics, vol. 74, no. 4, pp. 1041–1051, 2013.

X. Wang and H. Zhang, “Backstepping-based lag synchronization of a complex permanent magnet synchronous motor system,” Chinese Physics B, vol. 22, no. 4, Article ID 048902, 2013.