Next Article in Journal
An Iterative Algorithm for the Nonlinear MC2 Model with Variational Inequality Method
Previous Article in Journal
The Generalized Viscosity Implicit Midpoint Rule for Nonexpansive Mappings in Banach Space
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Complex Dynamical Behaviors of Lorenz-Stenflo Equations

1
Chongqing Key Laboratory of Social Economy and Applied Statistics, College of Mathematics and Statistics, Chongqing Technology and Business University, Chongqing 400067, China
2
Mathematical Postdoctoral station, College of Mathematics and Statistics, Southwest University, Chongqing 400716, China
3
College of Automation, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(6), 513; https://doi.org/10.3390/math7060513
Submission received: 30 May 2019 / Revised: 30 May 2019 / Accepted: 31 May 2019 / Published: 5 June 2019
(This article belongs to the Section Engineering Mathematics)

Abstract

:
A mathematical chaos model for the dynamical behaviors of atmospheric acoustic-gravity waves is considered in this paper. Boundedness and globally attractive sets of this chaos model are studied by means of the generalized Lyapunov function method. The innovation of this paper is that it not only proves this system is globally bounded but also provides a series of global attraction sets of this system. The rate of trajectories entering from the exterior of the trapping domain to its interior is also obtained. Finally, the detailed numerical simulations are carried out to justify theoretical results. The results in this study can be used to study chaos control and chaos synchronization of this chaos system.

1. Introduction

In 1963, E.N. Lorenz [1] obtained the famous Lorenz chaotic system to describe weather changes. The Lorenz chaotic system stimulated the interest of researchers to study chaotic systems and chaotic phenomena. Since then, some new chaotic systems, which are not equivalent to the Lorenz system, were found, i.e., Rössler system [2,3], Chua system [4], Chen system [5], Lu system [6], Shimizu-Morioka system [7], hyperchaos Lorenz system [7], etc. Chaotic systems have potential applications in different scientific areas, such as secure communications [8], image encryption [9], laser technology [10], physical systems [11], circuits system [12] and chaos synchronization [13,14,15,16,17].
In 1996, the physicist Lennart Stenflo [18] obtained a new four-dimensional continuous-time dynamical chaotic system by adding a new variable w to the Lorenz system to describe the complex dynamical behaviors of the atmospheric acoustic-gravity waves. The Lorenz-Stenflo system can be described by the following ordinary differential equations [18]:
{ d x d t = a ( y x ) + r w , d y d t = c x y x z , d z d t = x y b z , d w d t = x a w ,
where a , b , c , r are parameters of system (1) and a > 0 , b > 0 , c > 0 , r > 0 . a is the Prandtl number of system (1), c is the Rayleigh number of system (1), b is the geometric parameter of system (1), and r is the rotation parameter of system (1). The Lorenz-Stenflo system is a four-dimensional continuous-time dynamical system that can describe the dynamical behaviors of the atmospheric acoustic-gravity waves in a rotating atmosphere [18]. It is important for us to study acoustic gravity waves because they are associated with minor weather changes and large-scale phenomena. Dynamical behaviors of the Lorenz-Stenflo system, such as periodicity [19,20], bifurcation phenomenon [21,22], synchronization behaviors [23], and chaotic behaviors, have been studied by many researchers. System (1) has chaotic attractors when a = 1 , b = 0.7 , c = 26 , r = 1.5 . Figure 1 shows the chaotic attractors of the Lorenz-Stenflo system (1) in three-dimensional (3D) spaces.
Boundedness of chaotic systems is an important concept in dynamical systems [24,25], which plays an important role in chaos control and chaos synchronization. Boundedness of the Lorenz system has been investigated by Leonov et al. in a series of articles [26,27,28]. One can show that there is a bounded ellipsoid in R 3 , which all orbits of the Lorenz system with an exponential rate will eventually enter.
Despite the fact that many qualitative results on the Lorenz-Stenflo system have been obtained [19,20,21,22,23], there is a fundamental question that has not been completely answered so far: Is there a global trapping region for the Lorenz-Stenflo system? How to get the boundedness of a chaotic system is particularly significant both for theoretical research and engineering applications [13,24,29,30]. Motivated by the above discussion, we will investigate the ultimate bound set and the global attractive sets of the Lorenz-Stenflo system. The novelty of this paper is that this paper not only proves the Lorenz-Stenflo system is globally bounded for all parameters of the system by means of the generalized Lyapunov function method, but it also gives a family of mathematical expressions of globally exponential attractive sets for the Lorenz-Stenflo system, with respect to all parameters of this system. The innovation of this paper is that this paper not only proves the Lorenz-Stenflo system is globally bounded but also provides a series of global attraction sets of this Lorenz-Stenflo system. The rate of trajectories entering from the exterior of the trapping domain to its interior is also obtained.

2. Boundedness and Global Attraction

The ultimate bound set and global domain of attraction of the Lorenz-Stenflo system will be studied in the following part. We can get the following Theorem 1 and Theorem 2.
Theorem 1.
For any a > 0 , b > 0 , c > 0 , r > 0 , there exists a positive number L > 0 , such that
Ψ λ , m = { ( x , y , z , w ) | λ ( x m 1 ) 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r ( w m 3 ) 2 L }
is the ultimate bound set and positively invariant set of the Lorenz-Stenflo system (1), where X ( t ) = ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) , λ > 0 , m > 0 , and m 1 , m 3 are arbitrary real numbers.
Proof. 
Define the generalized Lyapunov function
V ( X ) = λ ( x m 1 ) 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r ( w m 3 ) 2 ,
where λ > 0 , m > 0 , X ( t ) = ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) , and m 1 , m 3 are arbitrary real numbers.
And we can get
d V ( X ( t ) ) d t | ( 1 ) = 2 λ ( x m 1 ) d x d t + 2 m y d y d t + 2 m ( z λ a + m c m ) d z d t + 2 λ r ( w m 3 ) d w d t , = 2 λ ( x m 1 ) ( a y a x + r w ) + 2 m y ( c x y x z ) + 2 m ( z λ a + m c m ) ( x y b z ) + 2 λ r ( w m 3 ) ( x a w ) , = 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 λ ( r m 3 + a m 1 ) x 2 a λ m 1 y + 2 b ( λ a + m c ) z + 2 λ r ( a m 3 m 1 ) w .
Obviously, the set Γ defined by
{ X | 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 λ ( r m 3 + a m 1 ) x 2 a λ m 1 y + 2 b ( λ a + m c ) z + 2 λ r ( a m 3 m 1 ) w    = 0 . }
is an ellipsoid in R 4 for a > 0 , b > 0 , c > 0 , r > 0 , λ > 0 , m > 0 . Outside Γ , d V ( X ( t ) ) d t < 0 , while inside Γ , d V ( X ( t ) ) d t > 0 . Thus, boundedness of system (1) can only be reached on Γ . Obviously, the continuous function (3) can reach its maximum value max X Γ V ( X ) = L on the bounded closed set Γ . Clearly, { X | V ( X ) max X Γ V ( X ) = L , X Γ } contains solutions to system (1). It is easy to prove that Ψ λ , m is the ultimate bound set and positively invariant set for system (1).
The proof is completed. □
Although Theorem 1 claims that there exist the ultimate bound set and positively invariant set for system (1), it is very difficult to calculate the mathematical analytic expression of the positive number L > 0 , since there are too many parameters in (2).We can get the expression of (2) easily for the special case m 1 = 0 , m 3 = 0 . We can get the following theorem for the special case m 1 = 0 , m 3 = 0 .
Theorem 2.
For any a > 0 , b > 0 , c > 0 , r > 0 , such that
Ω λ , m = { ( x , y , z , w ) | λ x 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r w 2 R 2 }
is the ultimate bound set and positively invariant set of Lorenz-Stenflo system (1), where λ > 0 , m > 0 ,
R 2 = { b 2 ( λ α + m c ) 2 4 m a ( b a ) , a 1 , b 2 a , b 2 ( λ a + m c ) 2 4 m ( b 1 ) , a > 1 , b 2 , ( λ α + m c ) 2 m , b < 2 a , b < 2 .
Proof. 
Define the generalized Lyapunov function
V 1 ( X ) = λ x 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r w 2 ,
where λ > 0 , m > 0 , X ( t ) = ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) .
And we can get
d V 1 ( X ( t ) ) d t | ( 1 ) = 2 λ x d x d t + 2 m y d y d t + 2 m ( z λ a + m c m ) d z d t + 2 λ r w d w d t , = 2 λ x ( a y a x + r w ) + 2 m y ( c x y x z ) + 2 m ( z λ a + m c m ) ( x y b z ) + 2 λ r w ( x a w ) , = 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 b ( λ a + m c ) z .
Obviously,
Γ 1 = { ( x , y , z , w ) | 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 b ( λ a + m c ) z = 0 }
is an ellipsoid in R 4 for a > 0 , b > 0 , c > 0 , r > 0 , λ > 0 , m > 0 . Outside Γ 1 , d V 1 ( X ( t ) ) d t < 0 , while inside Γ 1 , d V 1 ( X ( t ) ) d t > 0 . Thus, the ultimate boundedness for system (1) can only be reached on the bounded closed set Γ 1 . Obviously, the continuous function V 1 ( X ) can reach its maximum value max X Γ V 1 ( X ) = R 2 on the bounded closed set Γ 1 . In order to calculate mathematical analytic expression of max X Γ V 1 ( X ) = R 2 , we have to solve the following optimization problem:
{ max   V 1 ( X ) = max { λ x 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r w 2 } , s . t . 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 b ( λ a + m c ) z = 0 ,
The above optimization problem is equivalent to
{ max   V 1 ( X ) = max { λ x 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r w 2 } , s . t . λ x 2 b ( λ a + m c ) 2 4 a m + m y 2 b ( λ a + m c ) 2 4 m + m ( z λ a + m c 2 m ) 2 ( λ a + m c ) 2 4 m + λ r w 2 b ( λ a + m c ) 2 4 a m = 1 .
Let us take
λ x = x 1 , m y = y 1 , m z = z 1 , λ r w = w 1 ,
Then, the above optimization problem (9) is changed into
{ max   V 1 ( X ) = max { x 1 2 + y 1 2 + ( z 1 λ a + m c m ) 2 + w 1 2 } , s . t . x 1 2 b ( λ a + m c ) 2 4 a m + y 1 2 b ( λ a + m c ) 2 4 m + ( z 1 λ a + m c 2 m ) 2 ( λ a + m c ) 2 4 m + w 1 2 b ( λ a + m c ) 2 4 a m = 1 ,
We can get the solution of optimization problem (10) as
R 2 = { b 2 ( λ α + m c ) 2 4 m a ( b a ) , a 1 , b 2 a , b 2 ( λ a + m c ) 2 4 m ( b 1 ) , a > 1 , b 2 , ( λ α + m c ) 2 m , b < 2 a , b < 2 .
The proof is completed. □
Remark 1.
(i) Let us take λ > 0 , m > 0 in Theorem 2, then we can get a series of ultimate bound and positively invariant sets of the Lorenz-Stenflo system (1), according to Theorem 2;
(ii) Theorem 2 gives a family of mathematical expressions of the ultimate bound sets of the Lorenz-Stenflo system with respect to the parameters of the system. Particularly, we can get the following conclusion of Theorem 2for the special case λ = 1 , m = 1 :
For any a > 0 , b > 0 , c > 0 , r > 0 , such that
Ω 2 = { ( x , y , z , w ) | x 2 + y 2 + ( z a c ) 2 + r w 2 l 2 }
is the ultimate bound set and positively invariant set of the Lorenz-Stenflo system (1), where
l 2 = { b 2 ( α + c ) 2 4 a ( b a ) , a 1 , b 2 a , b 2 ( a + c ) 2 4 ( b 1 ) , a > 1 , b 2 , ( α + c ) 2 , b < 2 a , b < 2 .
Figure 2 shows the chaotic attractor of the Lorenz-Stenflo system (1) in xOyz space, defined by Ω 2 in (11) with a = 1 , b = 0.7 , c = 26 , r = 1.5 .
Remark 2.
From Figure 2, we can see that numerical simulation is consistent with theoretical analysis results (11).
Though Theorem 1 and Theorem 2 give the ultimate bound and positively invariant set of the Lorenz-Stenflo system (1), it does not give the rate of the trajectories entering from the exterior of the trapping domain to its interior. The rate of trajectories entering from the exterior of the trapping domain to its interior of the Lorenz-Stenflo system (1) is described by the following Theorem 3:
Theorem 3.
For any a > 0 , b > 0 , c > 0 , r > 0 , and let
M 2 = 1 θ [ a λ ( m 1 ) 2 + λ r 2 ( m 3 ) 2 a + a 2 λ 2 ( m 1 ) 2 m + b ( λ a + m c ) 2 m + a λ r ( m 3 ) 2 + λ r ( m 1 ) 2 a ] ,
V ( X ) = λ ( x m 1 ) 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r ( w m 3 ) 2 ,
λ > 0 , m > 0 , m 1 R , m 3 R ,
θ = min ( a , b , 1 ) > 0 ,   X ( t ) = ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) .
Then we can obtain the following exponential estimate for the trajectory of system (1)
[ V ( X ( t ) ) M 2 ] [ V ( X ( t 0 ) ) M 2 ] e θ ( t t 0 ) .
Thus, Δ λ , m = { X | V ( X ) M 2 } is the globally exponential attractive set of system (1), i.e.,
lim ¯ t + V ( X ( t ) ) M 2 .
Proof. 
Define
f ( x ) = a λ x 2 + 2 λ r m 3 x , h ( y ) = m y 2 2 a λ m 1 y , g ( w ) = a λ r w 2 2 λ r m 1 w ,
then we can get
max x R f ( x ) = λ r 2 ( m 3 ) 2 a , max y R h ( y ) = a 2 λ 2 m 1 2 m , max w R g ( w ) = λ r ( m 1 ) 2 a .
Define
V ( X ) = λ ( x m 1 ) 2 + m y 2 + m ( z λ a + m c m ) 2 + λ r ( w m 3 ) 2 ,
where λ > 0 , m > 0 , X ( t ) = ( x ( t ) , y ( t ) , z ( t ) , w ( t ) ) , and m 1 , m 3 are arbitrary real numbers.
And we have
d V ( X ( t ) ) d t | ( 1 ) = 2 λ ( x m 1 ) d x d t + 2 m y d y d t + 2 m ( z λ a + m c m ) d z d t + 2 λ r ( w m 3 ) d w d t , = 2 λ ( x m 1 ) ( a y a x + r w ) + 2 m y ( c x y x z ) + 2 m ( z λ a + m c m ) ( x y b z )     + 2 λ r ( w m 3 ) ( x a w ) , = 2 a λ x 2 2 m y 2 2 m b z 2 2 a λ r w 2 + 2 λ ( r m 3 + a m 1 ) x 2 a λ m 1 y     + 2 b ( λ a + m c ) z + 2 λ r ( a m 3 m 1 ) w , = 2 a λ x 2 + 2 λ ( r m 3 + a m 1 ) x 2 m y 2 2 a λ m 1 y 2 m b z 2   + 2 b ( λ a + m c ) z     2 a λ r w 2 + 2 λ r ( a m 3 m 1 ) w , = a λ x 2 + 2 λ a m 1 x a λ x 2 + 2 λ r m 3 x 2 m y 2 2 a λ m 1 y 2 m b z 2   + 2 b ( λ a + m c ) z a λ r w 2 + 2 λ r a m 3 w a λ r w 2 2 λ r m 1 w , = a λ x 2 + 2 λ a m 1 x + f ( x ) 2 m y 2 2 a λ m 1 y 2 m b z 2   + 2 b ( λ a + m c ) z a λ r w 2 + 2 λ r a m 3 w + g ( w ) , a λ x 2 + 2 λ a m 1 x + max x R f ( x ) 2 m y 2 2 a λ m 1 y 2 m b z 2   + 2 b ( λ a + m c ) z a λ r w 2 + 2 λ r a m 3 w + max w R g ( w ) , = a λ x 2 + 2 λ a m 1 x + λ r 2 ( m 3 ) 2 a 2 m y 2 2 a λ m 1 y 2 m b z 2   + 2 b ( λ a + m c ) z a λ r w 2 + 2 λ r a m 3 w + λ r ( m 1 ) 2 a , = a λ ( x m 1 ) 2 + a λ m 1 2 + λ r 2 ( m 3 ) 2 a m y 2 + h ( y ) 2 m b z 2   + 2 b ( λ a + m c ) z a λ r w 2 + 2 λ r a m 3 w + λ r ( m 1 ) 2 a , a λ ( x m 1 ) 2 + a λ m 1 2 + λ r 2 ( m 3 ) 2 a m y 2 + max y R h ( y ) m b z 2   + 2 b ( λ a + m c ) z a λ r ( w m 3 ) 2 + a λ r m 3 2 + λ r ( m 1 ) 2 a , = a λ ( x m 1 ) 2 + a λ m 1 2 + λ r 2 ( m 3 ) 2 a m y 2 + a 2 λ 2 m 1 2 m m b z 2   + 2 b ( λ a + m c ) z a λ r ( w m 3 ) 2 + a λ r m 3 2 + λ r ( m 1 ) 2 a , = a λ ( x m 1 ) 2 + a λ m 1 2 + λ r 2 ( m 3 ) 2 a m y 2 + a 2 λ 2 m 1 2 m m b ( z λ a + m c m ) 2 + b ( λ a + m c ) 2 m a λ r ( w m 3 ) 2 + a λ r m 3 2 + λ r ( m 1 ) 2 a , = a λ ( x m 1 ) 2 m y 2 m b ( z λ a + m c m ) 2 a λ r ( w m 3 ) 2 + a λ m 1 2 + λ r 2 ( m 3 ) 2 a + a 2 λ 2 ( m 1 ) 2 m + b ( λ a + m c ) 2 m + a λ r ( m 3 ) 2 + λ r ( m 1 ) 2 a , = θ V ( X ) + a λ m 1 2 + λ r 2 ( m 3 ) 2 a + a 2 λ 2 ( m 1 ) 2 m + b ( λ a + m c ) 2 m + a λ r ( m 3 ) 2 + λ r ( m 1 ) 2 a , = θ [ V ( X ( t ) ) M 2 ] .
Thus, we have the exponential estimate for the trajectory of system (1):
[ V ( X ( t ) ) M 2 ] [ V ( X ( t 0 ) ) M 2 ] e θ ( t t 0 ) .
and
lim ¯ t + V ( X ( t ) ) M 2 .
Hence,
Δ λ , m = { X | V ( X ) M 2 }
is the globally exponential attractive set of system (1).
The proof is completed. □
Remark 3.
(i) Let us take λ > 0 , m > 0 in Theorem 3, then we can get a series of global exponential attractive sets of the Lorenz-Stenflo system (1), according to Theorem 3;
(ii) Theorem 3 not only gives a family of mathematical expressions of global exponential attractive sets for the Lorenz-Stenflo system with respect to the parameters of this system, but also provides an explicit exponential rate of trajectories entering from the exterior of the trapping domain to its interior.

3. Conclusions

The main goal of this paper is to study the Lorenz-Stenflo chaos system, describing the dynamical behavior of atmospheric acoustic-gravity waves. A family of hyperelliptic estimates of the ultimate bound and positively invariant set for the Lorenz-Stenflo system is obtained by means of the generalized Lyapunov function method. Furthermore, the rate of trajectories entering from the exterior of the trapping domain to its interior is also obtained. This generalized Lyapunov function method can be extended to study other chaotic systems. The results in this study provide a new and in-depth understanding of the Lorenz-Stenflo chaotic system. The results in this paper can be used to study the Lyapunov dimension of attractors, the Hausdorff dimension of attractors, chaos control, and chaos synchronization of this Lorenz-Stenflo chaos system in the future.

Author Contributions

The authors have made the same contribution. All authors read and approved the final manuscript.

Funding

This work is supported by the Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant Nos. KJQN201800818, KJ1500605), the National Natural Science Foundation of China (Grant No. 11871122), China Postdoctoral Science Foundation (Grant No. 2016M590850), and Chongqing Postdoctoral Science Foundation Special Funded Project (Grant No. Xm2017174) and the Program for University Innovation Team of Chongqing (Grant No: CXTDX201601026).

Acknowledgments

The authors express their sincere gratitude to the editors for the careful reading of the original manuscript and useful comments, which helped to improve the presentation of the results and accentuate important details. The authors will thank Gennady A. Leonov in Russian Academy of Sciences, Jinhu Lü in Chinese Academy of Sciences, Qigui Yang in South China University of Technology, Xiaofeng Liao in Chongqing University, Gaoxiang Yang in Ankang University, and Min Xiao in Nanjing University of Posts and Telecommunications very much for their help.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Lorenz, E.N. Deterministic non-periods flows. J. Atmos. Sci. 1963, 20, 130–141. [Google Scholar] [CrossRef]
  2. Rössler, O.E. An equation for hyperchaos. Phys. Lett. A 1979, 71, 155–157. [Google Scholar] [CrossRef]
  3. Kuznetsov, N.; Mokaev, T.; Vasilyev, P. Numerical justification of Leonov conjecture on Lyapunov dimension of Rössler attractor. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 1027–1034. [Google Scholar] [CrossRef]
  4. Chua, L.O.; Komura, M.; Matsumoto, T. The double scroll family. IEEE Trans. Circuits Syst. 1986, 33, 1072–1097. [Google Scholar] [CrossRef]
  5. Chen, G.; Ueta, T. Yet another chaotic attractor. Int. J. Bifurc. Chaos 1999, 9, 1465–1466. [Google Scholar] [CrossRef]
  6. Lu, J.; Chen, G. A new chaotic attractor coined. Int. J. Bifurc. Chaos 2002, 12, 659–661. [Google Scholar] [CrossRef]
  7. Leonov, G. General existence conditions of homoclinic trajectories in dissipative systems. Shimizu-Morioka, Lu and Chen systems. Phys. Lett. A 2012, 376, 3045–3050. [Google Scholar] [CrossRef]
  8. Wang, X.Y.; Wang, M.J. A hyperchaos generated from Lorenz system. Phys. A 2008, 387, 3751–3758. [Google Scholar] [CrossRef]
  9. Wu, X.; Zhu, C.; Kan, H. An improved secure communication scheme based passive synchronization of hyperchaotic complex nonlinear system. Appl. Math. Comput. 2015, 252, 201–214. [Google Scholar] [CrossRef]
  10. Zhang, Y.; Wang, X. A new image encryption algorithm based on non-adjacent coupled map lattices. Appl. Soft Comput. 2015, 26, 10–20. [Google Scholar] [CrossRef]
  11. Pisarchik, A.N.; Arecchi, F.T.; Meucci, R.; Garbo, A.D. Synchronization of Shilnikov chaos in CO2 laser with feedback. Laser Phys. 2014, 11, 1235–1239. [Google Scholar]
  12. Zhang, F.C.; Liao, X.F.; Zhang, G.Y. Qualitative behaviors of the continuous-time chaotic dynamical systems describing the interaction of waves in plasma. Nonlinear Dyn. 2017, 88, 1623–1629. [Google Scholar] [CrossRef]
  13. Kuznetsov, N.V.; Leonov, G.A.; Yuldashev, M.V.; Yuldashev, R.V. Hidden attractors in dynamical models of phase-locked loop circuits: Limitations of simulation in MATLAB and SPICE. Commun. Nonlinear Sci. Numer. Simul. 2017, 51, 39–49. [Google Scholar] [Green Version]
  14. Zhou, P.; Zhu, P. A practical synchronization approach for fractional-order chaotic systems. Nonlinear Dyn. 2017, 89, 1719–1726. [Google Scholar] [CrossRef]
  15. Zhou, P.; Cai, H.; Yang, C.D. Stabilization of the unstable equilibrium points of the fractional-order BLDCM chaotic system in the sense of Lyapunov by a single-state variable. Nonlinear Dyn. 2016, 84, 2357–2361. [Google Scholar] [CrossRef]
  16. Ren, F.; Cao, J. Anti-synchronization of stochastic perturbed delayed chaotic neural networks. Neural Comput. Appl. 2009, 18, 515–521. [Google Scholar] [CrossRef]
  17. Ahmad, I.; Saaban, A.; Ibrahim, A.; Shahzad, M. Robust Finite-Time Anti-Synchronization of Chaotic Systems with Different Dimensions. Mathematics 2015, 3, 1222–1240. [Google Scholar] [CrossRef] [Green Version]
  18. Zhu, X.H.; Du, W.S. A New Family of Chaotic Systems with Different Closed Curve Equilibrium. Mathematics 2019, 7, 94. [Google Scholar] [CrossRef]
  19. Stenflo, L. Generalized Lorenz equations for acoustic-gravity waves in the atmosphere. Phys. Scr. 1996, 53, 83–84. [Google Scholar] [CrossRef]
  20. Yu, M.Y.; Yang, B. Periodic and chaotic solutions of the generalized Lorenz equations. Phys. Scr. 1996, 54, 140–142. [Google Scholar] [CrossRef]
  21. Park, J.; Lee, H.; Jeon, Y.L.; Baik, J.J. Periodicity of the Lorenz–Stenflo equations. Phys. Scr. 2015, 90, 065201. [Google Scholar] [CrossRef]
  22. Yu, M.Y.; Zhou, C.T.; Lai, C.H. The bifurcation characteristics of the generalized Lorenz equations. Phys. Scr. 1996, 54, 321–324. [Google Scholar] [CrossRef]
  23. Zhou, C.; Lai, C.H.; Yu, M.Y. Bifurcation behavior of the generalized Lorenz equations at large rotation numbers. J. Math. Phys. 1997, 38, 5225–5239. [Google Scholar] [CrossRef]
  24. Chen, Y.; Shi, Z.; Lin, C. Some criteria for the global finite-time synchronization of two Lorenz–Stenflo systems coupled by a new controller. Appl. Math. Model. 2014, 38, 4076–4085. [Google Scholar] [CrossRef]
  25. Zhang, F.C.; Mu, C.L.; Zhou, S.M.; Zheng, P. New results of the ultimate bound on the trajectories of the family of the Lorenz systems. Discrete Contin. Dyn. Syst. Ser. B 2015, 20, 1261–1276. [Google Scholar]
  26. Zhang, F.C.; Liao, X.F.; Mu, C.L.; Zhang, G.Y.; Chen, Y.A. On global boundedness of the Chen system. Discret. Contin. Dyn. Syst. Ser. B 2017, 22, 1673–1681. [Google Scholar] [CrossRef] [Green Version]
  27. Leonov, G. Bounds for attractors and the existence of homoclinic orbits in the Lorenz system. J. Appl. Math. Mech. 2001, 65, 19–32. [Google Scholar] [CrossRef]
  28. Leonov, G.; Bunin, A.; Koksch, N. Attractor localization of the Lorenz system. Z. Angew. Math. Mech. 1987, 67, 649–656. [Google Scholar] [CrossRef]
  29. Zhang, F.C.; Zhang, G.Y. Further results on ultimate bound on the trajectories of the Lorenz system. Qual. Theory Dyn. Syst. 2016, 15, 221–235. [Google Scholar] [CrossRef]
  30. Leonov, G.; Kuznetsov, N. Hidden attractors in dynamical systems. From hidden oscillations in Hilbert–Kolmogorov, Aizerman, and Kalman problems to hidden chaotic attractor in Chua circuits. Int. J. Bifurc. Chaos 2013, 23, 1330002. [Google Scholar] [CrossRef]
Figure 1. Chaotic attractors of the Lorenz-Stenflo system (1) in three-dimensional (3D) spaces with a = 1 , b = 0.7 , c = 26 , r = 1.5 .
Figure 1. Chaotic attractors of the Lorenz-Stenflo system (1) in three-dimensional (3D) spaces with a = 1 , b = 0.7 , c = 26 , r = 1.5 .
Mathematics 07 00513 g001
Figure 2. Chaotic attractor of the Lorenz-Stenflos ystem (1) in xOyz space, defined by Ω 2 in (11) with a = 1 , b = 0.7 , c = 26 , r = 1.5 .
Figure 2. Chaotic attractor of the Lorenz-Stenflos ystem (1) in xOyz space, defined by Ω 2 in (11) with a = 1 , b = 0.7 , c = 26 , r = 1.5 .
Mathematics 07 00513 g002

Share and Cite

MDPI and ACS Style

Zhang, F.; Xiao, M. Complex Dynamical Behaviors of Lorenz-Stenflo Equations. Mathematics 2019, 7, 513. https://doi.org/10.3390/math7060513

AMA Style

Zhang F, Xiao M. Complex Dynamical Behaviors of Lorenz-Stenflo Equations. Mathematics. 2019; 7(6):513. https://doi.org/10.3390/math7060513

Chicago/Turabian Style

Zhang, Fuchen, and Min Xiao. 2019. "Complex Dynamical Behaviors of Lorenz-Stenflo Equations" Mathematics 7, no. 6: 513. https://doi.org/10.3390/math7060513

APA Style

Zhang, F., & Xiao, M. (2019). Complex Dynamical Behaviors of Lorenz-Stenflo Equations. Mathematics, 7(6), 513. https://doi.org/10.3390/math7060513

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop