Research Article | Open Access

# Numerical Simulation of a Class of Hyperchaotic System Using Barycentric Lagrange Interpolation Collocation Method

**Academic Editor:**Diyi Chen

#### Abstract

Hyperchaotic system, as an important topic, has become an active research subject in nonlinear science. Over the past two decades, hyperchaotic system between nonlinear systems has been extensively studied. Although many kinds of numerical methods of the system have been announced, simple and efficient methods have always been the direction that scholars strive to pursue. Based on this problem, this paper introduces another novel numerical method to solve a class of hyperchaotic system. Barycentric Lagrange interpolation collocation method is given and illustrated with hyperchaotic system () as examples. Numerical simulations are used to verify the effectiveness of the present method.

#### 1. Introduction

Many chaotic systems have been developed such as Lorenz system [1], Rossler system [2], and Chen system [3]. As chaos theory progresses, many new chaotic systems [4–8] have been proposed, specially hyperchaotic systems [9–15]. A hyperchaotic system is usually characterized as a chaotic system with more than one positive Lyapunov exponent, implying that the dynamics expand in more than one direction, giving rise to more complex chaotic dynamics. Barycentric interpolation collocation method [16, 17] is a high precision method. Some authors have used barycentric interpolation collocation method to solve various kinds of problems [16–23]. This paper suggests the barycentric interpolation collocation method to solve a class of hyperchaotic system, and a hyperchaotic system (1) is adopted as an example to elucidate the solution process.

We consider the following 4D butterfly hyperchaotic system with butterfly phenomenon [24]:where are the state variables and are the positive constant parameters of the system which satisfy the following initial conditions:

#### 2. The Numerical Solution of System (1)

First of all, we give initial function and construct the following linear iterative format of system (1):

Next, we use the barycentric Lagrange interpolation collocation method to solve (3).

In the interval takes different nodes, . The barycentric interpolation of can be written as [16, 17]

is, respectively, barycentric Lagrange interpolation primary function and is center of gravity interpolation weight.

Use formula (4), the functions can be expressed as

So, linear iterative format (3) can be written in following partitioned matrix form:The matrix is order matrix. is order unit matrix, diagonal matrix , and diagonal matrix . The vector At last, we use initial conditions (2).

Take formula (4) into initial conditions (2); we can get the following discrete equations of initial conditions:

In this paper, we use displacement method to impose the initial conditions. The detailed procedure is as follows.

The first of (6) are replaced separately by the equation of initial conditions (8) in turn.

So, we can get that are approximate solution of (1) and (2).

#### 3. Numerical Experiment

In this section, six numerical experiments are studied to demonstrate the effectiveness of the present method. All experiments are computed using MatlabR2017a. In Experiments 1–6, we choose Chebyshev nodes, the accuracy of iteration control is , and the initial iteration value ; . Parameters of the numerical Experiments 1–5 are listed in Table 1.

*Experiment 1. *We consider the following hyperchaotic system [25]:where are the state variables and are the positive parameters of the system, which satisfy the following initial conditions:We choose Chebyshev nodes; the number of nodes . Numerical results of Experiment 1 are given in Figures 1 and 2.

**(a)**

**(b)**

Figure 1 is states of the hyperchaotic system for Experiment 1 with , which is obtained by using the current method, and is the states of and and is the states of and . Figure 2 is hyperchaotic attractors of the system for Experiment 1 with , which is obtained by using the current method. Among them, is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph in three-dimensional -space.

*Experiment 2. *We consider the following hyperchaotic system [26]:where are the state variables and are the positive parameters of the system, which satisfy the following initial conditions:

We choose Chebyshev nodes, the number of nodes , and the parameters . Numerical results of Experiment 2 are given in Figures 3–6.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

Figure 3 is phase portraits of a new hyperchaotic system for Experiment 2 with by using the current method. is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the three-dimensional space graph. Figures 4 and 5 are phase portraits of a new hyperchaotic system for Experiment 2 obtained by using the current method with and , respectively. Figure 6 is time series plots of a new hyperchaotic system for Experiment 2 with different parameter value . and represent time series when ; and represent time series when ; and represent time series when .

*Experiment 3. *We consider the following butterfly hyperchaotic system [24]:where are the state variables and are the positive constant parameters of the system, which satisfy the following initial conditions:

We choose Chebyshev nodes; the number of nodes . Numerical results of Experiment 3 are given in Figures 7 and 8.

**(a)**

**(b)**

Figure 7 is states of a novel butterfly hyperchaotic system for Experiment 3 with , which is obtained by using the current method, and is the states of and and is the states of and . Figure 8 is phase portraits of a novel butterfly hyperchaotic system for Experiment 3 with , which is obtained by using the current method. Among them, is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph in three-dimensional -space.

*Experiment 4. *We consider the following hyperchaotic Chen system [27]:where are the state variables and are the positive constant parameters of the system, which satisfy the following initial conditions:

We choose Chebyshev nodes; the number of nodes . Numerical results of Experiment 4 are given in Figures 9 and 10.

**(a)**

**(b)**

Figure 9 is time response of the hyperchaotic Chen system’s variable states for Experiment 4 with , which is obtained by using the current method, and is the states of and and is the states of and . Figure 10 is phase portraits of the hyperchaotic Chen system for Experiment 4 with , which is obtained by using the current method. Among them, is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph in three-dimensional -space; is the graph in three-dimensional -space.

*Experiment 5. *We consider the following hyperchaotic system [28]:where are the state variables and are the positive constant parameters of the system, which satisfy the following initial conditions:

We choose Chebyshev nodes; the number of nodes . Numerical results of Experiment 5 are given in Figures 11 and 12.

**(a)**

**(b)**

Figure 11 is the time series plots of a new hyperchaotic system for Experiment 5 with , which is obtained by using the current method, and is the states of and and is the states of and . Figure 12 is hyperchaotic attractors of a new hyperchaotic system for Experiment 5 with , which is obtained by using the current method. Among them, is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane.

*Experiment 6. *We consider the following 4D hyperchaotic system [29]:where are the state variables and are the positive constant parameters of the system, which satisfy the following initial conditions:

We choose Chebyshev nodes; the number of nodes . Numerical results of Experiment 6 are given in Figures 13–15.

Figure 13 is phase portraits of a 4D hyperchaotic system for Experiment 6 with , , , , , and , which is obtained by using the current method. Among them, is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane; is the graph projected on -plane. Figure 14 is phase portraits of a 4D hyperchaotic system for Experiment 6 with , , , , , and , which is obtained by using the current method. Figure 15 is the time series plots of a 4D hyperchaotic system for Experiment 6. and are obtained by using the current method with , , , , , and . and are obtained by using the current method with , , , , , and .

#### 4. Conclusions and Remarks

In this paper, a class of hyperchaotic system has been solved by using barycentric Lagrange interpolation collocation method. The numerical simulation results are in accord with the theoretical analyses and circuit implementation. Numerical simulations are provided to verify the effectiveness and feasibility of the proposed numerical results, which are in agreement with theoretical analysis. In the further work, we will be devoted to studying fractional-order hyperchaotic system.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

#### Acknowledgments

This paper is supported by the Natural Science Foundation of Inner Mongolia 2017MS0103, Inner Mongolia maker Collaborative Innovation Center of Jining Normal University, and the National Natural Science Foundation of China 11361037.

#### References

- E. N. Lorenz, “Deterministic non-periodic flow,”
*Journal of the Atmospheric Sciences*, vol. 20, pp. 130–141, 1963. View at: Publisher Site | Google Scholar - O. E. Rössler, “An equation for continuous chaos,”
*Physics Letters A*, vol. 57, no. 5, pp. 397-398, 1976. View at: Publisher Site | Google Scholar - G. Chen and T. Ueta, “Yet another chaotic attractor,”
*International Journal of Bifurcation and Chaos*, vol. 9, no. 7, pp. 1465-1466, 1999. View at: Publisher Site | Google Scholar | MathSciNet - D. Chen, R. Zhang, J. C. Sprott, H. Chen, and X. Ma, “Synchronization between integer-order chaotic systems and a class of fractional-order chaotic systems via sliding mode control,”
*Chaos: An Interdisciplinary Journal of Nonlinear Science*, vol. 22, no. 2, Article ID 023130, 2012. View at: Publisher Site | Google Scholar - A. E. Matouk, A. A. Elsadany, E. Ahmed, and H. N. Agiza, “Dynamical behavior of fractional-order Hastings-Powell food chain model and its discretization,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 27, no. 1–3, pp. 153–167, 2015. View at: Publisher Site | Google Scholar | MathSciNet - D. Y. Chen, R. F. Zhang, J. C. Sprott, and X. Y. Ma, “Synchronization between integer-order chaotic systems and a class of fractional-order chaotic system based on fuzzy sliding mode control,”
*Nonlinear Dynamics*, vol. 70, no. 2, pp. 1549–1561, 2012. View at: Publisher Site | Google Scholar - D. Y. Chen, W. L. Zhao, J. C. Sprott, and X. Ma, “Application of Takagi-Sugeno fuzzy model to a class of chaotic synchronization and anti-synchronization,”
*Nonlinear Dynamics*, vol. 73, no. 3, pp. 1495–1505, 2013. View at: Publisher Site | Google Scholar | MathSciNet - D. Y. Chen, C. Wu, H. H. C. Iu, and X. Ma, “Circuit simulation for synchronization of a fractional-order and integer-order chaotic system,”
*Nonlinear Dynamics*, vol. 73, no. 3, pp. 1671–1686, 2013. View at: Publisher Site | Google Scholar | MathSciNet - Z. Chen, Y. Yang, G. Qi, and Z. Yuan, “A novel hyperchaos system only with one equilibrium,”
*Physics Letters A*, vol. 360, no. 6, pp. 696–701, 2007. View at: Publisher Site | Google Scholar | MathSciNet - X. Wang and M. Wang, “A hyperchaos generated from Lorenz system,”
*Physica A: Statistical Mechanics and Its Applications*, vol. 387, no. 14, pp. 3751–3758, 2008. View at: Publisher Site | Google Scholar | MathSciNet - A. E. Matouk, “Dynamics and control in a novel hyperchaotic system,”
*International Journal of Dynamics and Control*, 2018. View at: Publisher Site | Google Scholar - N. Yujun, W. Xingyuan, W. Mingjun, and Z. Huaguang, “A new hyperchaotic system and its circuit implementation,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 15, no. 11, pp. 3518–3524, 2010. View at: Publisher Site | Google Scholar - A. M. A. El-Sayed, H. M. Nour, A. Elsaid, A. E. Matouk, and A. Elsonbaty, “Circuit realization, bifurcations, chaos and hyperchaos in a new 4D system,”
*Applied Mathematics and Computation*, vol. 239, pp. 333–345, 2014. View at: Publisher Site | Google Scholar | MathSciNet - S. Pang and Y. Liu, “A new hyperchaotic system from the Lü system and its control,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 8, pp. 2775–2789, 2011. View at: Publisher Site | Google Scholar | MathSciNet - X. Liu, X. Shen, and H. Zhang, “Multi-scroll chaotic and hyperchaotic attractors generated from Chen system,”
*International Journal of Bifurcation and Chaos*, vol. 22, Article ID 1250033, 2012. View at: Publisher Site | Google Scholar - S. P. Li and Z. Q. Wang,
*Barycentric Interpolation Collocation Method for Nonlinear Problems*, National Defense Industry Press, Beijing, China, 2015. - S. P. Li and Z. Q. Wang,
*High-Precision Non-Grid Center of Gravity Interpolation Collocation Method: Algorithm, Program And Engineering Application*, Science Press, Beijing, China, 2012. - H. Liu, J. Huang, Y. Pan, and J. Zhang, “Barycentric interpolation collocation methods for solving linear and nonlinear high-dimensional Fredholm integral equations,”
*Journal of Computational and Applied Mathematics*, vol. 327, pp. 141–154, 2018 (Basque). View at: Publisher Site | Google Scholar | MathSciNet - W.-H. Luo, T.-Z. Huang, X.-M. Gu, and Y. Liu, “Barycentric rational collocation methods for a class of nonlinear parabolic partial differential equations,”
*Applied Mathematics Letters*, vol. 68, pp. 13–19, 2017. View at: Publisher Site | Google Scholar | MathSciNet - F. F. Liu, Y. L. Wang, and S. G. Li, “Barycentric interpolation collocation method for solving the coupled viscous Burgers' equations,”
*International Journal of Computer Mathematics*, vol. 95, pp. 2162–2173, 2018. View at: Publisher Site | Google Scholar | MathSciNet - Y. L. Wang, D. Tian, and Z. Y. Li, “Numerical method for singularly perturbed delay parabolic partial differential equations,”
*Thermal Science*, vol. 21, no. 4, pp. 1595–1599, 2017. View at: Publisher Site | Google Scholar - S. P. Li and Z. Q. Wang, “Barycentric interpolation collocation method for solving elastic problems,”
*Zhongnan Daxue Xuebao (Ziran Kexue Ban)/Journal of Central South University (Science and Technology)*, vol. 44, pp. 2031–2040, 2013. View at: Google Scholar - S. P. Li and Z. Q. Wang, “Barycentric interpolation collocation method for solving nonlinear vibration problems,”
*Noise Vibration Control*, vol. 28, pp. 49–52, 2018. View at: Google Scholar - L. Zhang, “A novel 4-D butterfly hyperchaotic system,”
*Optik - International Journal for Light and Electron Optics*, vol. 131, pp. 215–220, 2017. View at: Publisher Site | Google Scholar - J. P. Singh and B. K. Roy, “A novel hyperchaotic system with stable and unstable line of equilibria and sigma shaped Poincare map,”
*IFAC-PapersOnLine*, vol. 49, no. 1, pp. 526–531, 2016. View at: Publisher Site | Google Scholar - A. E. Matouk, “On the periodic orbits bifurcating from a fold Hopf bifurcation in two hyperchaotic systems,”
*Optik - International Journal for Light and Electron Optics*, vol. 126, no. 24, pp. 4890–4895, 2015. View at: Publisher Site | Google Scholar - D. Sadaoui, A. Boukabou, and S. Hadef, “Predictive feedback control and synchronization of hyperchaotic systems,”
*Applied Mathematics and Computation*, vol. 247, pp. 235–243, 2014. View at: Publisher Site | Google Scholar | MathSciNet - C.-L. Li, J.-B. Xiong, and W. Li, “A new hyperchaotic system and its generalized synchronization,”
*Optik-International Journal for Light and Electron Optics*, vol. 125, no. 1, pp. 575–579, 2014. View at: Publisher Site | Google Scholar - J. P. Singh, B. K. Roy, and S. Jafari, “New family of 4-D hyperchaotic and chaotic systems with quadric surfaces of equilibria,”
*Chaos, Solitons & Fractals*, vol. 106, pp. 243–257, 2018. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2019 Xiaofei Zhou 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.