Although the globally attractive sets of a hyperchaotic system have important applications in the fields of engineering, science, and technology, it is often a difficult task for the researchers to obtain the globally attractive set of the hyperchaotic systems due to the complexity of the hyperchaotic systems. Therefore, we will study the globally attractive set of a generalized hyperchaotic Lorenz–Stenflo system describing the evolution of finite amplitude acoustic gravity waves in a rotating atmosphere in this paper. Based on Lyapunov-like functional approach combining some simple inequalities, we derive the globally attractive set of this system with its parameters. The effectiveness of the proposed methods is illustrated via numerical examples.
National Natural Science Foundation of China11501064Chongqing Municipal Education CommissionKJ1500605Chongqing Technology and Business University2014-56-11China Postdoctoral Science Foundation2016M590850Chongqing Postdoctoral Science Foundation Special Funded ProjectXm2017174Chongqing Technology and Business University17520731. Introduction
In 1963, Lorenz found the well-known three-dimensional Lorenz model when he studied the dynamics of the atmosphere [1]. Since then, various complex dynamical behaviors of the Lorenz system have been studied by mathematicians, physicists, and engineers from various fields due to various applications in the fields of engineering, science, and technology [2–14]. In order to improve the stability or predictability of the Lorenz system, Stenflo and Leonov derived the following four-dimensional Lorenz–Stenflo system with four parameters to describe the dynamics of the atmosphere [15, 16]:(1)dxdt=ay-x+dw,dydt=cx-y-xz,dzdt=xy-bz,dwdt=-x-aw.
In order to give a better description of the atmosphere, Chen and Liang propose a generalized Lorenz–Stenflo system with six parameters according to the Lorenz–Stenflo system [17]:(2)dxdt=ay-x+sw,dydt=cx-dy-xz,dzdt=xy-bz,dwdt=-x-rw,where x, y, z, and w are state variables and a, b, c, d, r, and s are positive parameters of system (2). System (2) can describe the dynamic behavior of finite amplitude acoustic gravity waves in a rotating atmosphere.
The Lyapunov exponents of the dynamical system (2) are calculated numerically for the parameter values a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64 with the initial state x0,y0,z0,w0=2.2,2.0,10.5,20. System (2) has Lyapunov exponents as λLE1=0.0696, λLE2=0.0359, λLE3=0.0002, and λLE4=-24.5176 for the parameters listed above (see [18, 19] for a detailed discussion of Lyapunov exponents of strange attractors in dynamical systems). Thus, system (2) has two positive Lyapunov exponents and the strange attractor, which means system (2) can exhibit a variety of interesting and complex chaotic behaviors. System (2) has a hyperchaotic attractor with a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64, as shown in Figures 1–4.
Projection of hyperchaotic attractor of system (2) onto the xOyz space with a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64.
Projection of hyperchaotic attractor of system (2) onto the xOyw space with a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64.
Projection of hyperchaotic attractor of system (2) onto the xOzw space with a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64.
Projection of hyperchaotic attractor of system (2) onto the yOzw space with a=19.42, b=1.91, c=29.45, d=2.86, r=0.23, and s=9.64.
In this paper, all the simulations are carried out by using fourth-order Runge-Kutta Method with time-step h=0.005.
The rest of this paper is organized as follows. In Section 2, the globally attractive set for the chaotic attractors in (2) is studied using Lyapunov stability theory. To validate the ultimate bound estimation, numerical simulations are also provided. Finally, the conclusions are drawn in Section 3.
2. Bounds for the Chaotic Attractors in System (2)Theorem 1.
For any a>0, b>0, c>0, d>0, r>0, s>0, there exists a positive number M>0, such that(3)Ωλ,m=X∣λx-m12+my2+mz-aλ+cmm2+λsw-m32≤M,∀λ>0,∀m>0is the ultimate bound set of system (2), where Xt=xt,yt,zt,wt.
Proof.
Define the following Lyapunov-like function:(4)Vλ,mX=Vλ,mx,y,z,w=λx-m12+my2+mz-aλ+cmm2+λsw-m32,where ∀λ>0,∀m>0,Xt=xt,yt,zt,vt,ut,ωt, and m1∈R, m3∈R are arbitrary constants.
And we can get(5)dVλ,mXtdt2=2λx-m1dxdt+2mydydt+2mz-aλ+cmmdzdt+2λsw-m3dwdt=2λx-m1ay-ax+sw+2mycx-dy-xz+2mz-aλ+cmmxy-bz+2λsw-m3-x-rw=-2aλx2+2aλm1+λsm3x-2dmy2-2aλm1y-2bmz2+2bmm2z-2λsrw2-2λsm1-λsrm3w.Let dVXt/dt=0. Then, we can get the surface Γ:(6)X∣-aλx2+aλm1+λsm3x-dmy2-aλm1y-bmz2+bmm2z-λsrw2-λsm1-λsrm3w=0is an ellipsoid in R4 for ∀λ>0,∀m>0,a>0,b>0,c>0,d>0,r>0,s>0. Outside Γ, dVλ,mXt/dt<0, while inside Γ, dVλ,mXt/dt>0. Thus, the ultimate boundedness for system (2) can only be reached on Γ. Since the Lyapunov-like function Vλ,mX is a continuous function and Γ is a bounded closed set, then the function (4) can reach its maximum value maxX∈ΓVλ,mX=M on the surface Γ that is defined in (6). Obviously, X∣Vλ,mX≤maxX∈ΓVλ,mX=M,X∈Γ contains solutions of system (2). It is obvious that the set Ωλ,m is the ultimate bound set for system (2).
This completes the proof.
Theorem 2.
Suppose that ∀a>0,b>0,d>0,r>0,c>0,s>0,λ>0,m>0.
Let xt,yt,zt,wt be an arbitrary solution of system (2) and (7)Lλ,m2=1θa2λ2md+λsr+aλm12+bmm22+λs2a+λsrm32,θ=mina,b,d,r>0,Vλ,mX=Vλ,mx,y,z,w=λx-m12+my2+mz-aλ+cmm2+λsw-m32,∀λ>0,∀m>0,∀m1∈R,∀m3∈R.Then the estimation(8)Vλ,mXt-Lλ,m2≤Vλ,mXt0-Lλ,m2e-θt-t0holds for system (2), and thus Ωλ,m=X∣Vλ,mX≤Lλ,m2 is the globally exponential attractive set of system (2); that is, lim-t→+∞Vλ,mXt≤Lλ,m2.
Proof.
Define the following functions:(9)fx=-aλx2+2λsm3x,hy=-dmy2-2aλm1y,gw=-λsrw2-2λsm1w.then we can get(10)maxx∈Rfx=λs2m32a,maxy∈Rhy=a2λ2m12dm,maxw∈Rgw=λsm12r.Construct the Lyapunov-like function(11)Vλ,mX=Vλ,mx,y,z,w=λx-m12+my2+mz-aλ+cmm2+λsw-m32,∀λ>0,∀m>0,∀m1∈R,∀m3∈R.Differentiating the above Lyapunov-like function Vλ,mX in (11) with respect to time t along the trajectory of system (2) yields(12)dVλ,mXtdt2=2λx-m1dxdt+2mydydt+2mz-aλ+cmmdzdt+2λsw-m3dwdt=2λx-m1ay-ax+sw+2mycx-dy-xz+2mz-aλ+cmmxy-bz+2λsw-m3-x-rw=-2aλx2+2aλm1+λsm3x-2dmy2-2aλm1y-2bmz2+2bmm2z-2λsrw2-2λsm1-λsrm3w=-aλx2+2aλm1x-aλx2+2λsm3x-dmy2-dmy2-2aλm1y-bmz2-bmz2+2bmm2z-λsrw2-2λsm1w-λsrw2+2λsrm3w=-aλx2-2m1x-aλx2+2λsm3x-dmy2-dmy2-2aλm1y-bmz2-2m2z-bmz2-λsrw2-2m3w-λsrw2-2λsm1w=-aλx-m12+aλm12+fx-dmy2+hy-bmz-m22+bmm22-bmz2-λsrw-m32+λsrm32+gw=-aλx-m12-dmy2-bmz-m22-λsrw-m32+fx+hy+gw-bmz2+aλm12+bmm22+λsrm32≤-θVλ,mX+maxx∈Rfx+maxy∈Rhy+maxw∈Rgw+aλm12+bmm22+λsrm32=-θVλ,mX+λs2m32a+a2λ2m12md+λsm12r+aλm12+bmm22+λsrm32=-θVλ,mX-Lλ,m2.Thus, we have(13)Vλ,mXt-Lλ,m2≤Vλ,mXt0-Lλ,m2e-θt-t0,limt→+∞¯Vλ,mXt≤Lλ,m2,which clearly shows that Ωλ,m=X∣Vλ,mX≤Lλ,m2 is the globally exponential attractive set of system (2).
The proof is complete.
Remark 3.
(i) In particular, let us take m1=0, m3=0 in Theorem 2, we can get the conclusions below.
Suppose that ∀a>0,b>0,d>0,r>0,c>0,s>0,λ>0,m>0.
Let xt,yt,zt,wt be an arbitrary solution of system (2) and (14)Mλ,m2=baλ+cm2θm,θ=mina,b,d,r>0,Vλ,mX=Vλ,mx,y,z,w=λx2+my2+mz-aλ+cmm2+λsw2,∀λ>0,∀m>0.Then the estimation(15)Vλ,mXt-Mλ,m2≤Vλ,mXt0-Mλ,m2e-θt-t0holds for system (2), and thus (16)Σλ,m=x,y,z,w∣λx2+my2+mz-aλ+cmm2+λsw2≤Mλ,m2,∀λ>0,∀m>0is the globally exponential attractive set and positive invariant set of system (2); that is, (17)limt→+∞¯Vλ,mXt≤Mλ,m2.
(ii) Let us take m1=0,m3=0,λ=1,m=1; then we can get(18)Σ1,1=x,y,z,w∣x2+y2+z-a-c2+sw2≤ba+c2mina,b,d,ras the globally exponential attractive set and positive invariant set of system (2) according to Theorem 2.
When a=19.42,b=1.91,c=29.45,d=2.86,r=0.23,s=9.64, we can get that (19)Ω1,1=x,y,z,w∣x2+y2+z-48.872+9.64w2≤140.82as the globally exponential attractive set and positive invariant set of system (2).
Figure 5 shows hyperchaotic attractor of system (2) in the xOyz space defined by Ω1,1. Figure 6 shows hyperchaotic attractor of system (2) in the xOyw space defined by Ω1,1. Figure 7 shows hyperchaotic attractor of system (2) in the xOzw space defined by Ω1,1. Figure 8 shows hyperchaotic attractor of system (2) in the yOzw space defined by Ω1,1.
Localization of hyperchaotic attractor of system (2) in the xOyz space defined by Ω1,1.
Localization of hyperchaotic attractor of system (2) in the xOyw space defined by Ω1,1.
Localization of hyperchaotic attractor of system (2) in the xOzw space defined by Ω1,1.
Localization of hyperchaotic attractor of system (2) in the yOzw space defined by Ω1,1.
3. Conclusions
In this paper, we have investigated some global dynamics of a generalized Lorenz–Stenflo system describing the evolution of finite amplitude acoustic gravity waves in a rotating atmosphere. Based on the Lyapunov method, the globally attractive sets were formulated combining simple inequalities. Finally, numerical examples were presented to show the effectiveness of the proposed method.
Disclosure
All authors have read and approved the final manuscript.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported by National Natural Science Foundation of China (Grant no. 11501064), the Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant no. KJ1500605), the Research Fund of Chongqing Technology and Business University (Grant no. 2014-56-11), China Postdoctoral Science Foundation (Grant no. 2016M590850), Chongqing Postdoctoral Science Foundation Special Funded Project (Grant no. Xm2017174), and the Research Fund of Chongqing Technology and Business University (Grant no. 1752073). The authors thank Professors Jinhu Lu in Institute of Systems Science, Chinese Academy of Sciences, Gaoxiang Yang in Ankang University, Ping Zhou in Chongqing University of Posts and Telecommunications, and Min Xiao in Nanjing University of Posts and Telecommunications for their help.
LorenzE. N.Deterministic nonperiodic flow196320213014110.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2LeonovG. A.KuznetsovN. V.On differences and similarities in the analysis of Lorenz, Chen, and LU systems2015256334343MR331607210.1016/j.amc.2014.12.132Zbl1338.370462-s2.0-84922423521WangX.WangM.A hyperchaos generated from Lorenz system20083871437513758MR258681110.1016/j.physa.2008.02.0202-s2.0-41649108176LeonovG. A.General existence conditions of homoclinic trajectories in dissipative systems. Lorenz, Shimizu-Morioka, LU and Chen systems20123764530453050MR299718010.1016/j.physleta.2012.07.003Zbl1266.340792-s2.0-84868088152ZhangF.ZhangG.Further results on ultimate bound on the trajectories of the Lorenz system2016151221235MR348401110.1007/s12346-015-0137-0Zbl1338.652752-s2.0-84962791905LeonovG. A.Bounds for attractors and the existence of homoclinic orbits in the lorenz system200165119322-s2.0-003564098410.1016/S0021-8928(01)00004-1Zbl1025.34048LeonovG. A.BuninA. I.KokschN.Attraktorlokalisierung des Lorenz-Systems1987671264965610.1002/zamm.19870671215MR928581ZhangF.MuC.ZhouS.ZhengP.New results of the ultimate bound on the trajectories of the family of the Lorenz systems2015204126112762-s2.0-8498304430910.3934/dcdsb.2015.20.1261Zbl1348.37024LeonovG. A.KuznetsovN. V.Hidden attractors in dynamical systems: from hidden oscillations in Hilbert-Kolmogorov, Aizerman, and KALman problems to hidden chaotic attractor in Chua circuits20132316910.1142/S02181274133000241330002MR3038624Zbl1270.340032-s2.0-84874642094LeonovG. A.KuznetsovN. V.KorzhemanovaN. A.KusakinD. V.Lyapunov dimension formula for the global attractor of the Lorenz system20164184103MR350810010.1016/j.cnsns.2016.04.0322-s2.0-84969764618LeonovG. A.KuznetsovN. V.MokaevT. N.Homoclinic orbits, and self-excited and hidden attractors in a Lorenz-like system describing convective fluid motion201522481421145810.1140/epjst/e2015-02470-32-s2.0-84937868918KuznetsovN. V.MokaevT. N.VasilyevP. A.Numerical justification of Leonov conjecture on LYApunov dimension of Rossler attractor201419410271034MR311927810.1016/j.cnsns.2013.07.0262-s2.0-84886303693KuznetsovN. V.LeonovG. A.YuldashevM. V.YuldashevR. V.Hidden attractors in dynamical models of phase-locked loop circuits: limitations of simulation in MATLAB and SPICE201751394910.1016/j.cnsns.2017.03.0102-s2.0-85016520472ElsayedE. M.AhmedA. M.Dynamics of a three-dimensional systems of rational difference equations201639510261038MR347413910.1002/mma.3540Zbl1339.390142-s2.0-84961145342StenfloL.Generalized Lorenz equations for acoustic-gravity waves in the atmosphere199653183842-s2.0-000280438210.1088/0031-8949/53/1/015LeonovG. A.Generalized Lorenz equations for acoustic-gravity waves in the atmosphere. Attractors dimension, convergence and homoclinic trajectories201716622532267MR369388110.3934/cpaa.2017111Zbl06790273ChenY.-M.LiangH.-H.Zero-zero-Hopf bifurcation and ultimate bound estimation of a generalized Lorenz-Stenflo hyperchaotic system2017401034243432MR366458610.1002/mma.4236Zbl1371.340542-s2.0-84996567131FredericksonP.KaplanJ. L.YorkeE. D.YorkeJ. A.The Liapunov dimension of strange attractors1983492185207MR70864210.1016/0022-0396(83)90011-62-s2.0-48749146954WolfA.SwiftJ. B.SwinneyH. L.VastanoJ. A.Determining Lyapunov exponents from a time series1985163285317MR80570610.1016/0167-2789(85)90011-92-s2.0-0008494528