Abstract
Although some numerical methods of the Lorenz 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 a novel numerical method to solve the Lorenz-type chaotic system which is based on barycentric Lagrange interpolation collocation method (BLICM). The system is adopted as an example to elucidate the solution process. Numerical simulations are used to verify the effectiveness of the present method.
1. Introduction
In 1963, Edward Lorenz developed a simplified mathematical model for atmospheric convection. The model is a system of three ordinary differential equations now known as the Lorenz equations (see [1–3]):with the initial conditions
where is proportional to the rate of convection, to the horizontal temperature variation, and to the vertical temperature variation. The constants are system parameters proportional to the Prandtl number, Rayleigh number, and certain physical dimensions of the layer itself.
As chaos theory progresses, many new Lorenz-type systems [4–6] have been proposed, specially Lorenz hyperchaotic systems [7–10]. The Lorenz system is widely used in electric circuits, chemical reactions, and forward osmosis. Although some numerical methods of the Lorenz system have been announced, simple and efficient methods have always been the direction that scholars strive to pursue.
With the development of numerical analysis, there are some high-precision methods, such as variational iteration method [11–13], BLICM [14–22], and so on [23]. J.P. Berrut [24–26] introduced barycentric Lagrange interpolation, [27, 28] studied numerical stability of barycentric Lagrange interpolation, and [15, 16] give algorithm of BLICM. Some authors [14, 17–22] have used BLICM to solve all sorts of problems and show the BLICM is a high precision numerical method. This paper suggests the BLICM to solve the Lorenz system. The system (1) is adopted as an example to elucidate the solution process.
2. The Numerical Solution of the System (1)
First of all, we give initial function and construct following linear iterative format of system (1)Next, we use BLICM to solve (3).
Using the barycentric Lagrange interpolation functions [14–16, 24–26], we can get following:
where is respectively barycentric interpolation primary function, is center of gravity Lagrange interpolation weight, and .
Substitute formulae (4) and (5) into iterative format (3) and let . So, linear iterative format (3) can be written in following partitioned matrix form:where is order unit matrix, is order matrix, and the vectorThe vector
The first line of (7) is replaced separately by the equation of initial conditions (6) in turn.
So, we can get that are approximate solution of (1) and (2).
3. Numerical Experiment
In this section, some numerical examples are studied to find some new chaotic behaviors and verify the existing chaotic dynamic behaviors. In Experiments 1–5, the accuracy of iteration control is , the initial iteration value , and for parameters , and see Table 1.
Experiment 1. We consider the model (1) with and the initial conditions [4].
We choose Chebyshev nodes, and the number of nodes . Figure 1 is obtained by using the current method with . Among them, (a) is the time series plot; (b) is the phase diagram of ; (c) is the three-dimensional space graph; (d) is the graph projected on -plane; (e) is the graph projected on -plane. Figures 2 and 3 are obtained by using the current method at and , respectively. We can see that the fluctuation amplitude of and increases, while the fluctuation amplitude of decreases with the increase of . The corresponding graphs , and also have obvious changes.

(a)

(b)

(c)

(d)

(e)

(a)

(b)

(c)

(d)

(e)

(a)

(b)

(c)

(d)

(e)
Experiment 2. We consider the Lorenz-type system [6]
We choose Chebyshev nodes, the number of nodes , and the parameters and the initial conditions .
Figure 4 is strange attractors of chaotic system for Experiment 2 at by using the current method. (a) is the graph projected on -plane; (b) is the graph projected on -plane; (c) is the graph projected on -plane. Figure 7 is time series plots of chaotic system for Experiment 2 at different parameter value . We can see that the frequency of fluctuations of , and accelerates obviously with the increase of . When , the fluctuations of , and change obviously, and their fluctuations become smaller and smaller and finally stop at a certain value. Figures 5 and 6 are strange attractors of chaotic system for Experiment 2 by using the current method at and respectively.

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)
Experiment 3. We consider the 3D autonomous chaotic Lorenz-type system [7]
We choose Chebyshev nodes and the number of nodes and the parameters and the initial conditions .
Figure 8 is phase portraits of the 3D chaotic Lorenz type system for Experiment 3 at by using the current method. (a) is the graph projected on -plane; (b) is the graph projected on -plane; (c) is the graph projected on -plane. Figure 11 is time series plots of the 3D chaotic Lorenz type system for Experiment 3 at different parameter value . We can see that the fluctuation range of changes obviously with the increase of . When , the fluctuation range of is to , and when , the fluctuation range of is to . Figures 9 and 10 are phase portraits of the 3D chaotic Lorenz type system for Experiment 3 by using the current method at and , respectively.

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)

(a)

(b)

(c)
Experiment 4. We consider the Lorenz system [5]where , and are real parameters, which satisfy the following initial conditions:
We choose Chebyshev nodes and the number of nodes . Figure 12 is obtained by using the current method with the parameters , and . In Figure 12, (a) is the time series plot of ; (b) is the time series plot of ; (c) is the three-dimensional space graph; (d) is the graph projected on -plane; (e) is the graph projected on -plane; (f) is the graph projected on -plane.

(a)

(b)

(c)

(d)

(e)

(f)
Experiment 5. We consider the new chaotic system [8]where , and are state variables and , and are real parameters, which satisfy the following initial conditions:
We choose Chebyshev nodes and the number of nodes . Figure 13 is obtained by using the current method with the parameters , , and . In Figure 13, (a) is the time series plot; (b) is the three-dimensional space graph; (c) is the graph projected on -plane; (d) is the graph projected on -plane; (e) is the graph projected on -plane.

(a)

(b)

(c)

(d)

(e)
4. Conclusions and Remarks
In this paper, the Lorenz System has solved by using BLICM. These numerical experiments illustrate that the numerical results of the present method are the same as the experimental results.
All computations are performed by the MatlabR2017b software packages.
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 , Inner Mongolia maker collaborative innovation center of Jining Normal University, and the National Natural Science Foundation of China .