
In this paper, we present a 4D hyperchaotic Rabinovich system which obtained by adding a linear controller to 3D Rabinovich system. Based on theoretical analysis and numerical simulations, the rich dynamical phenomena such as boundedness, dissipativity and invariance, equilibria and their stability, chaos and hyperchaos are studied. In addition, the Hopf bifurcation at the zero equilibrium point of the 4D Rabinovich system is investigated. The numerical simulations, including phase diagrams, Lyapunov exponent spectrum, bifurcations, power spectrum and Poincaré maps, are carried out in order to analyze and verify the complex phenomena of the 4D Rabinovich system.
Citation: Junhong Li, Ning Cui. A hyperchaos generated from Rabinovich system[J]. AIMS Mathematics, 2023, 8(1): 1410-1426. doi: 10.3934/math.2023071
[1] | Ning Cui, Junhong Li . A new 4D hyperchaotic system and its control. AIMS Mathematics, 2023, 8(1): 905-923. doi: 10.3934/math.2023044 |
[2] | Junhong Li, Ning Cui . Dynamical behavior and control of a new hyperchaotic Hamiltonian system. AIMS Mathematics, 2022, 7(4): 5117-5132. doi: 10.3934/math.2022285 |
[3] | Figen Kangalgil, Seval Isșık . Effect of immigration in a predator-prey system: Stability, bifurcation and chaos. AIMS Mathematics, 2022, 7(8): 14354-14375. doi: 10.3934/math.2022791 |
[4] | Sukono, Siti Hadiaty Yuningsih, Endang Rusyaman, Sundarapandian Vaidyanathan, Aceng Sambas . Investigation of chaos behavior and integral sliding mode control on financial risk model. AIMS Mathematics, 2022, 7(10): 18377-18392. doi: 10.3934/math.20221012 |
[5] | Minghung Lin, Yiyou Hou, Maryam A. Al-Towailb, Hassan Saberi-Nik . The global attractive sets and synchronization of a fractional-order complex dynamical system. AIMS Mathematics, 2023, 8(2): 3523-3541. doi: 10.3934/math.2023179 |
[6] | Binhao Hong, Chunrui Zhang . Bifurcations and chaotic behavior of a predator-prey model with discrete time. AIMS Mathematics, 2023, 8(6): 13390-13410. doi: 10.3934/math.2023678 |
[7] | Wei Li, Qingkai Xu, Xingjian Wang, Chunrui Zhang . Dynamics analysis of spatiotemporal discrete predator-prey model based on coupled map lattices. AIMS Mathematics, 2025, 10(1): 1248-1299. doi: 10.3934/math.2025059 |
[8] | A. E. Matouk . Chaos and hidden chaos in a 4D dynamical system using the fractal-fractional operators. AIMS Mathematics, 2025, 10(3): 6233-6257. doi: 10.3934/math.2025284 |
[9] | Xiaoyan Hu, Bo Sang, Ning Wang . The chaotic mechanisms in some jerk systems. AIMS Mathematics, 2022, 7(9): 15714-15740. doi: 10.3934/math.2022861 |
[10] | Zhao Li, Shan Zhao . Bifurcation, chaotic behavior and solitary wave solutions for the Akbota equation. AIMS Mathematics, 2024, 9(8): 22590-22601. doi: 10.3934/math.20241100 |
In this paper, we present a 4D hyperchaotic Rabinovich system which obtained by adding a linear controller to 3D Rabinovich system. Based on theoretical analysis and numerical simulations, the rich dynamical phenomena such as boundedness, dissipativity and invariance, equilibria and their stability, chaos and hyperchaos are studied. In addition, the Hopf bifurcation at the zero equilibrium point of the 4D Rabinovich system is investigated. The numerical simulations, including phase diagrams, Lyapunov exponent spectrum, bifurcations, power spectrum and Poincaré maps, are carried out in order to analyze and verify the complex phenomena of the 4D Rabinovich system.
The concept of the hyperchaos was first put forward by Rössler [1] in 1979. Any system with at least two positive Lyapunov exponents is defined as hyperchaotic. Compared to chaotic attractors, hyperchaotic attractors have more complicated dynamical phenomena and stronger randomness and unpredictability. Hyperchaotic systems have aroused wide interest from more and more researchers in the last decades. A number of papers have investigated various aspects of hyperchaos and many valuable results have been obtained. For instance, in applications, in order to improve the security of the cellular neural network system, the chaotic degree of the system can be enhanced by designing 5D memristive hyperchaotic system [2]. For higher computational security, a new 4D hyperchaotic cryptosystem was constructed by adding a new state to the Lorenz system and well used in the AMr-WB G.722.2 codec to fully and partially encrypt the speech codec [3]. In fact, hyperchaos has a wide range applications such as image encryption [4], Hopfield neural network [5] and secure communication [6] and other fields [7,8]. Meanwhile, there are many hyperchaotic systems have been presented so far. Aimin Chen and his cooperators constructed a 4D hyperchaotic system by adding a state feedback controller to Lü system [9]. Based on Chen system, Z. Yan presented a new 4D hyperchaotic system by introducing a state feedback controller [10]. By adding a controlled variable, Gao et al. introduced a new 4D hyperchaotic Lorenz system [11]. Likewise, researchers also formulated 5D Shimizu-Morioka-type hyperchaotic system [12], 5D hyperjerk hypercaotic system [13] and 4D T hyperchaotic system [14] and so on.
In [15], a chaotic Rabinovich system was introduced
{˙x=hy−ax+yz,˙y=hx−by−xz,˙z=−cz+xy, | (1.1) |
where (x,y,z)T∈R3 is the state vector. When (h,a,b,c)=(0.04,1.5,−0.3,1.67), (1.1) has chaotic attractor[15,16]. System (1.1) has similar properties to Lorenz system, the two systems can be considered as special cases of generalized Lorenz system in [17]. Liu and his cooperators formulated a new 4D hyperchaotic Rabinovich system by adding a linear controller to the 3D Rabinovich system [18]. The circuit implementation and the finite-time synchronization for the 4D hyperchaotic Rabinovich system was also studied in [19]. Reference [20] formulated a 4D hyperchaotic Rabinovich system and the dynamical behaviors were studied such as the hidden attractors, multiple limit cycles and boundedness. Based on the 3D chaotic Rabinovich system, Tong et al. presented a new 4D hyperchaotic system by introducing new state variable [21]. The hyperchaos can be generated by adding variables to a chaotic system, which has been verificated by scholars [3,9,10,11,14]. In [18,19,20,21], the hyperchaotic systems were presented by adding a variable to the second equation of system (1.1). In fact, hyperchaos can also be generated by adding a linear controller to the first equation and second equation of system (1.1). Based on it, the following hyperchaotic system is obtained
{˙x=hy−ax+yz+k0u,˙y=hx−by−xz+mu,˙z=−dz+xy,˙u=−kx−ky, | (1.2) |
where k0,h,a,b,d,k,k0,m are positive parameters. Like most hyperchaotic studies (see [14,22,23,24] and so on), the abundant dynamical properties of system (1.2) are investigated by divergence, phase diagrams, equilibrium points, Lyapunov exponents, bifurcation diagram and Poincaré maps. The results show that the new 4D Rabinovich system not only exhibit hyperchaotic and Hopf bifurcation behaviors, but also has the rich dynamical phenomena including periodic, chaotic and static bifurcation. In addition, the 4D projection figures are also given for providing more dynamical information.
The rest of this paper is organized as follows: In the next section, boundedness, dissipativity and invariance, equilibria and their stability of (1.2) are discussed. In the third section, the complex dynamical behaviors such as chaos and hyperchaos are numerically verified by Lyapunov exponents, bifurcation and Poincaré section. In the fourth section, the Hopf bifurcation at the zero equilibrium point of the 4D Rabinovich system is investigated. In addition, an example is given to test and verify the theoretical results. Finally, the conclusions are summarized in the last section.
Theorem 2.1. If k0>m, system (1.2) has an ellipsoidal ultimate bound and positively invariant set
Ω={(x,y,z,u)|mx2+k0y2+(k0−m)[z−h(k0+m)k0−m]2+k0mu2k≤M}, |
where
M={14h2d2(m+k0)2(k0−m)a(d−a),(k0>m,d>a),14h2d2(m+k0)2(k0−m)b(−b+d),(k0>m,d>b),(m+k0)2h2k0−m,(k0>m). | (2.1) |
Proof. V(x,y,z,u)=mx2+k0y2+(k0−m)[z−h(k0+m)k0−m]2+k0mu2k. Differentiating V with respect to time along a trajectory of (1.2), we obtain
˙V(x,y,z,u)2=−amx2−by2k0+dhmz+dhzk0+dmz2−dz2k0. |
When ˙V(x,y,z,u)2=0, we have the following ellipsoidal surface:
Σ={(x,y,z,u)|amx2+bk0y2+d(k0−m)(z−hm+hk02k0−2m)2=dh2(k0+m)24(k0−m)}. |
Outside Σ that is,
amx2+bk0y2+d(k0−m)(z−hm+hk02k0−2m)2<dh2(k0+m)24(k0−m), |
˙V<0, while inside Σ, that is,
amx2+bk0y2+d(k0−m)(z−hm+hk02k0−2m)2>dh2(k0+m)24(k0−m), |
˙V>0. Thus, the ultimate bound for system (1.2) can only be reached on Σ. According to calculation, the maximum value of V on Σ is Vmax=14h2d2(m+k0)2(k0−m)a(d−a),(k0>m,d>a) and
Vmax=14h2d2(m+k0)2(k0−m)b(−b+d),(k0>m,d>b);Vmax=(m+k0)2h2k0−m,(k0>m). |
In addition, Σ⊂Ω, when a trajectory X(t)=(x(t),y(t),z(t),u(t)) of (1.2) is outside Ω, we get ˙V(X(t))<0. Then, limt→+∞ρ(X(t,t0,X0),Ω)=0. When X(t)∈Ω, we also get ˙V(X(t))<0. Thus, any trajectory X(t) (X(t)≠X0) will go into Ω. Therefore, the conclusions of theorem is obtained.
We can see that system (1.2) is invariant for the coordinate transformation
(x,y,z,u)→(−x,−y,z,−u). |
Then, the nonzero equilibria of (1.2) is symmetric with respect to z axis. The divergence of (1.2) is
∇W=∂˙x∂x+∂˙y∂y+∂˙z∂z+∂˙u∂u=−(a+b+d), |
system (1.2) is dissipative if and only if a+b+d>0. It shows that each volume containing the system trajectories shrinks to zero as t→∞ at an exponential rate −(a+b+d).
For m≤k0, system (1.2) only has one equilibrium point O(0,0,0,0) and the Jacobian matrix at O is
J=[−ah0k0h−b0m00−d0−k−k00]. |
Then, the characteristic equation is
(d+λ)(λ3+s2λ2+s1λ+s0), | (2.2) |
where
s2=a+b,s1=ab−h2+km+kk0,s0=akm+bkk0+hkm+hkk0. |
According to Routh-Hurwitz criterion [25], the real parts of eigenvalues are negative if and only if
d>0,a+b>0,(a+b)(ab−h2)+k(ak0+bm−hm−hk0)>0. | (2.3) |
When m>k0, system (1.2) has other two nonzero equilibrium points E1(x∗1,y∗1,z∗1,u∗1) and E2(−x∗1,−y∗1,z∗1,−u∗1), where
x∗1=√d(am+bk0+hm+hk0)m−k0,y∗1=−√d(am+bk0+hm+hk0)m−k0, |
z∗1=−am+bk0+hm+hk0m−k0,u∗1=−a+b+2hm−k0√d(am+bk0+hm+hk0)m−k0. |
The characteristic equation of Jacobi matrix at E1 and E2 is
λ4+δ3λ3+δ2λ2+δ1λ+δ0=0, |
where
δ3=a+b+d,δ2=x∗14d2+ab+ad+bd−h2+km+kk0, |
δ1=1d[x∗12(3x∗12+ad−bd+kk0−km)+k(m+k0)(d+h)]+abd+akm+bkk0−dh2, |
δ0=3kk0x∗12−3kmx∗12+akmd+bkk0d+hkk0d+hkmd. |
Based on Routh-Hurwitz criterion [25], the real parts of eigenvalues are negative if and only if
δ0>0,δ3>0,δ3δ2−δ1>0,δ3δ2δ1−δ21−δ0δ23>0. | (2.4) |
Therefore, we have:
Theorem 2.2. (I) When m≤k0, system (1.2) only has one equilibrium point O(0,0,0,0) and O is asymptotically stable if and only if (2.3) is satisfied.
(II) when m>k0, system (1.2) has two nonzero equilibria E1, E2, (2.4) is the necessary and sufficient condition for the asymptotically stable of E1 and E2.
When (b,d,h,a,k0,k,m)=(1,1,10,10,0.8,0.8,0.8), system (1.2) only has zero-equilibrium point E0(0,0,0,0,0), its corresponding characteristic roots are: −1,0.230,5.232,−16.462, E0 is unstable. The Lyapunov exponents are: LE1=0.199, LE2=0.083, LE3=0.000, LE4=−12.283, system (1.2) is hyperchaotic. Figure 1 shows the hyperchaotic attractors on z−x−y space and y−z−u space. The Poincaré mapping on the x−z plane and power spectrum of time series x(t) are depicted in Figure 2.
In the following, we fix b=1,d=1,h=10,a=10,k0=0.8,m=0.8, Figure 3 indicates the Lyapunov exponent spectrum of system (1.2) with respect to k∈[0.005,1.8] and the corresponding bifurcation diagram is given in Figure 4. These simulation results illustrate the complex dynamical phenomena of system (1.2). When k varies in [0.005,1.8], there are two positive Lyapunov exponents, system (1.2) is hyperchaoic as k varies.
Assume b=1,d=1,h=10,k=0.8,k0=0.8,m=0.8, the different Lyapunov exponents and dynamical properties with different values of parameter a are given in Table 1. It shows that system (1.2) has rich dynamical behaviors including periodic, chaos and hyperchaos with different parameters. The bifurcation diagram of system (1.2) with a∈[0.5,5] is given in Figure 5. Therefore, we can see that periodic orbits, chaotic orbits and hyperchaotic orbits can occur with increasing of parameter a. When a=1.08, Figure 6 indicates the (z,x,y,u) 4D surface of section and the location of the consequents is given in the (z,x,y) subspace and are colored according to their u value. The chaotic attractor and hyperchaotic attractor on y−x−z space and hyperchaotic attractor on y−z−u space are given in Figures 7 and 8, respectively. The Poincaré maps on x−z plane with a=2 and a=4 are depicted in Figure 9.
a | LE1 | LE2 | LE3 | LE4 | Dynamics |
1.08 | 0.000 | −0.075 | −4.661 | −2.538 | Periodic |
2 | 0.090 | 0.000 | −0.395 | −3.696 | Chaos |
4 | 0.325 | 0.047 | −0.000 | −6.373 | Hyperchaos |
Theorem 4.1. Suppose that ab>h2 and k0(a−h)+m(b−h)<0 are satisfied. Then, as m varies and passes through the critical value k=ah2+bh2−a2b−ab2ak0+bm−hm−hk0, system (1.2) undergoes a Hopf bifurcation at O(0,0,0,0).
Proof. Assume that system (1.2) has a pure imaginary root λ=iω,(ω∈R+). From (2.2), we get
s2ω2−s0=0,ω3−s1ω=0, |
then
ω=ω0=√ab−h2+k(k0+m), |
k=k∗=(a+b)(h2−ab)k0(a−h)+m(b−h). |
Substituting k=k∗ into (2.2), we have
λ1=iω,λ2=iω,λ3=−d,λ4=−(a+b). |
Therefore, when ab>h2, k0(a−h)+m(b−h)<0 and k=k∗, the first condition for Hopf bifurcation [26] is satisfied. From (2.2), we have
Re(λ′(k∗))|λ=iω0=h(k0+m)−ak0−bm2(s22+s1)<0, |
Thus, the second condition for a Hopf bifurcation [26] is also met. Hence, Hopf bifurcation exists.
Remark 4.1. When ab−h2+k(k0+m)≤0, system (1.2) has no Hopf bifurcation at the zero equilibrium point.
Theorem 4.2. When ab>h2 and k0(a−h)+m(b−h)<0, the periodic solutions of (1.2) from Hopf bifurcation at O(0,0,0,0) exist for sufficiently small
0<|k−k∗|=|k−(a+b)(h2−ab)k0(a−h)+m(b−h)|. |
And the periodic solutions have the following properties:
(I) if δ1g>0 (resp., δ1g<0), the hopf bifurcation of system (1.2) at (0,0,0,0) is non-degenerate and subcritical (resp. supercritical), and the bifurcating periodic solution exists for m>m∗ (resp., m<m∗) and is unstable (resp., stable), where
δ1g=14dδ(−k2δ01k02+ahδ01s1+kδ01k0s1−2dδ03δ05+2dδ04δ06),δ=ω0[(−ab+h2)(akk0+hkk0−hs1)+k(bm−hk0)(ha+h2+kk0)−k(ak−hm)(a2+ha−kk0+s1)],δ01=kω0(a2bm+a2hk−a2hk0+abhm+ah2k−ah2m−ah2k0+ak2k0−bkmk0−h3m−hkmk0+hkk02+bms1−hk0s1),δ03=ω0(a2bkm−2a2bkk0+a2hk2−a2hkk0+abhkm+2abhkk0+ah2k2−ah2km+ah2kk0−ak3k0+bk2mk0−h3km−2h3kk0+hk2mk0−hk2k02+2a2bs1−2ah2s1−bkms1+hkk0s1),δ04=(−ab+h2)s12+(a3b−a2h2+abh2−2abkm+2abkk0−ahk2+2ahkk0−bhkm−h4+h2km−h2kk0)s1−a2k3k0+abk2mk0−2abk2k02−ahk3k0+ahk2mk0−ahk2k02+bhk2mk0+h2k2mk0+h2k2k02,δ05=12d2+8s1[2ω0(−kk0aω0+hkk0ω0−hω0s1)+d(−hω02a−k2k02+kk0s1],δ06=12d2+8s1[d(−kk0aω0+hkk0ω0−hω0s1)−2ω0(k2k02+ahs1−kk0s1)]. |
(II) The period and characteristic exponent of the bifurcating periodic solution are:
T=2πω0(1+τ2∗ϵ2+O(ϵ4)),β=β2ϵ2+O(ϵ4), |
where ϵ=k−k∗μ2+O[(k−k∗)2] and
μ2=−ReC1(0)α′(0)=−(s22+s1)δ1gh(k0+m)−ak0−bm, |
τ2∗=δ2gω0−δ1g(ams2+bk0s2+hms2+hk0s2+ms1+k0s1)s1(h(k0+m)−ak0−bm), |
β2=δ1g, |
δ2g=14dδ(−k2δ02k02+ahδ02s1+kδ02k0s1−2δ03δ06d−2δ04δ05d). |
(III) The expression of the bifurcating periodic solution is
[xyzu]=[−kk0ϵcos(2πtT)−hω0ϵsin(2πtT)(kk0−s1)ϵcos(2πtT)−aω0ϵsin(2πtT)ϵ2[−kk0(kk0−s1)+hs1a2d+δ05−δ06sin(4πtT)]−k(a+h)ϵcos(2πtT)+kω0ϵsin(2πtT)]+O(ϵ3). |
Proof. Let k=k∗, by straightforward computations, we can obtain
t1=[ihω0−kk0kk0−s1+iaω00−(iω0+a+h)k],t3=[0010],t4=[ak0−hmbm−hk00h2−ab], |
which satisfy
Jt1=iω0t1,Jt3=−dt3,Jt4=−(a+b)t4. |
Now, we use transformation X=QX1, where X=(x,y,z,u)T, X1=(x1,y1,z1,u1)T, and
Q=[−kk0−hω00ak−hmkk0−s1−aω00bm−hk00010−k(a+h)kω00h2−ab], |
then, system (1.2) is transformed into
{˙x1=−ω0y1+F1(x1,y1,z1,u1),˙y1=ω0x1+F2(x1,y1,z1,u1),˙z1=−dz1+F3(x1,y1,z1,u1),˙u1=−(a+b)u1+F4(x1,y1,z1,u1), | (4.1) |
where
δ=ω0 [(−ab+h2)(akk0+hkk0−hs1)+k(bm−hk0)(ha+h2+kk0)−k(ak−hm) (a2+ha−kk0+s1)],
F1(x1,y1,z1,u1) =1δz1(f11x1+f12y1+f13u1),
f11=−kk0ω0(abh−ak2−h3+hkm)+ (kk0−s1)ω0(a2b−ah2−bkm+hkk0),
f12=−kk0ω0(abh−ak2−h3+hkm)+ (kk0−s1)ω0(a2b−ah2−bkm+hkk0),
f13=−kk0ω0(abh−ak2−h3+hkm)+ (kk0−s1)ω0(a2b−ah2−bkm+hkk0),
F2(x1,y1,z1,u1)= −1δz1(f21x1+f22y1+f23u1),
f21=(−ab+h2)(2k2k02−2kk0s1+s12) −k(a+h)(ak2k0−bkmk0−hkmk0+ hkk02+bms1−hk0s1),
f22=−aω0[abk(m−k0)−hk(ak0−bm)] +abs1−h2s1]−hω0[ak(ak+bk0+hk−hm)−h2k(m+k0)],
f23=(h2−ab)[hkk0(m−k0)−kk0(ak−bm)−s1(bm−hk0)]+k(a+h)[ak(ak−2hm) +bm(bm−2hk0)+h2(m2+k02)],
F3(x1,y1,z1,u1)= [−kk0x1−hω0y1+(ak−hm)u1][(kk0−s1)x1−aω0y1+(bm−hk0)u1],
F4(x1,y1,z1,u1)= 1δz1(f41x1+f42y1+f43u1),
f41=kk0 (hω0ka+h2ω0k+k2k0ω0)− (kk0−s1)kω0(a2+ah−kk0+s1),
f42= a3ω02k+a2ω02kh+h2ω02ka−ak2k0ω02+h3 ω02k+k2k0hω02+ akω02s1,
f43=−a2bkmω0−a2hk2ω0+ a2hkk0ω0−abhkmω0−ah2k2 ω0+ah2kmω0+ah2kk0ω0−ak3 k0ω0 +bk2mk0ω0+ h3kmω0+hk2mk0ω0−hk2k02ω0−bkm ω0s1+ hkk0ω0s1.
Furthermore,
g11=14[∂2F1∂x21+∂2F1∂y21+i(∂2F2∂x21+∂2F2∂y21)]=0, |
g02=14[∂2F1∂x21−∂2F1∂y21−2∂2F2∂x1∂y1+i(∂2F2∂x21−∂2F2∂y21+2∂2F1∂x1∂y1)]=0, |
g20=14[∂2F1∂x21−∂2F1∂y21+2∂2F2∂x1∂y1+i(∂2F2∂x21−∂2F2∂y21−2∂2F1∂x1∂y1)]=0, |
G21=18[∂3F1∂x31+∂3F2∂y31+∂3F1∂x1∂y21+∂3F2∂x21∂y1+i(∂3F2∂x31−∂3F2∂y31+∂3F2∂x1∂y21−∂3F1∂x21∂y1)]=0. |
By solving the following equations
[−d00−(a+b)][ω111ω211]=−[h111h211],[−d−2iω000−(a+b)−2iω0][ω120ω220]=−[h120h220], |
where
h111=12[(kk0+ah)s1−k2k02], |
h211=14(∂2F4∂x21+∂2F4∂y21)=0, |
h120=12[−k2k02−hs1a+kk0s1+(hkk0ω0−kk0aω0−hω0s1)i], |
h220=14(∂2F4∂x21−∂2F4∂y21−2i∂2F4∂x1∂y1)=0, |
one obtains
ω111=12d[hs1a−kk0(kk0−s1)],ω211=0,ω220=0, |
ω120=12d2+8s1{2ω0(−kk0aω0+hkk0ω0−hω0s1)+d(−hω02a−k2k02+kk0s1) |
+[d(−kk0aω0+hkk0ω0−hω0s1)−2ω0(k2k02+ahs1−kk0s1)]i}, |
G1110=12[(∂2F1∂x1∂z1+∂2F2∂y1∂z1)+i(∂2F2∂x1∂z1−∂2F1∂y1∂z1)]=12δ(δ01+δ02i), |
where
δ01=kω0(a2bm+a2hk−a2hk0+abhm+ah2k−ah2m−ah2k0+ak2k0−bkmk0−h3m−hkmk0+hkk02+bms1−hk0s1),δ02=(ab−h2)s12+(a3b−a2h2+abh2−2abkk0−ahk2+bhkm−h4+h2km+h2kk0)s1+k2k0(a2k−bma+2abk0+ahk−ahm+ahk0−bhm−h2m−h2k0),G2110=12[(∂2F1∂x1∂u1+∂2F2∂y1∂u1)+i(∂2F2∂x1∂u1−∂2F1∂y1∂u1)]=0,G1101=12[(∂2F1∂x1∂z1−∂2F2∂y1∂z1)+i(∂2F2∂x1∂z1+∂2F1∂y1∂z1)]=−12δ[δ03+δ04i], |
where
δ03=ω0(a2bkm−2a2bkk0+a2hk2−a2hkk0+abhkm+2abhkk0+ah2k2−ah2km+ah2kk0−ak3k0+bk2mk0−h3km−2h3kk0+hk2mk0−hk2k02+2a2bs1−2ah2s1−bkms1+hkk0s1),δ04=(−ab+h2)s12+(a3b−a2h2+abh2−2abkm+2abkk0−ahk2+2ahkk0−bhkm−h4+h2km−h2kk0)s1−a2k3k0+abk2mk0−2abk2k02−ahk3k0+ahk2mk0−ahk2k02+bhk2mk0+h2k2mk0+h2k2k02, |
G2101=12[(∂2F1∂x1∂u1−∂2F2∂y1∂u1)+i(∂2F2∂x1∂u1+∂2F1∂y1∂u1)]=0, |
g21=G21+2∑j=1(2Gj110ωj11+Gj101ωj20)=δ1g+δ2gi, |
where
δ1g=14dδ(−k2δ01k02+ahδ01s1+kδ01k0s1−2dδ03δ05+2dδ04δ06), |
δ2g=14dδ(−k2δ02k02+ahδ02s1+kδ02k0s1−2δ03δ06d−2δ04δ05d), |
δ05=12d2+8s1[2ω0(−kk0aω0+hkk0ω0−hω0s1)+d(−hω02a−k2k02+kk0s1)], |
δ06=12d2+8s1[d(−kk0aω0+hkk0ω0−hω0s1)−2ω0(k2k02+ahs1−kk0s1)]. |
Based on above calculation and analysis, we get
C1(0)=i2ω0(g20g11−2|g11|2−13|g02|2)+12g21=12g21, |
μ2=−ReC1(0)α′(0)=−(s22+s1)δ1gh(k0+m)−ak0−bm, |
τ2∗=δ2gω0−δ1g(ams2+bk0s2+hms2+hk0s2+ms1+k0s1)s1(h(k0+m)−ak0−bm), |
where
ω′(0)=ω0(ams2+bk0s2+hms2+hk0s2+ms1+k0s1)s1s22+s21, |
α′(0)=h(k0+m)−ak0−bm2(s22+s1),β2=2ReC1(0)=δ1g. |
Note α′(0)<0. From ab>h2 and k0(a−h)+m(b−h)<0, we obtain that if δ1g>0 (resp., δ1g<0), then μ2>0 (resp., μ2<0) and β2>0 (resp., β2<0), the hopf bifurcation of system (1.2) at (0,0,0,0) is non-degenerate and subcritical (resp. supercritical), and the bifurcating periodic solution exists for k>k∗ (resp., k<k∗) and is unstable (resp., stable).
Furthermore, the period and characteristic exponent are
T=2πω0(1+τ2∗ϵ2+O(ϵ4)),β=β2ϵ2+O(ϵ4), |
where ϵ=k−k∗μ2+O[(k−k∗)2]. And the expression of the bifurcating periodic solution is (except for an arbitrary phase angle)
X=(x,y,z,u)T=Q(¯y1,¯y2,¯y3,¯y4)T=QY, |
where
¯y1=Reμ,¯y1=Imμ,(¯y3,¯y4,¯y5)T=ω11|μ|2+Re(ω20μ2)+O(|μ|2), |
and
μ=ϵe2itπT+iϵ26ω0[g02e−4itπT−3g20e4itπT+6g11]+O(ϵ3)=ϵe2itπT+O(ϵ3). |
By computations, we can obtain
[xyzu]=[−kk0ϵcos(2πtT)−hω0ϵsin(2πtT)(kk0−s1)ϵcos(2πtT)−aω0ϵsin(2πtT)ϵ2[−kk0(kk0−s1)+hs1a2d+δ05−δ06sin(4πtT)]−k(a+h)ϵcos(2πtT)+kω0ϵsin(2πtT)]+O(ϵ3). |
Based on the above discussion, the conclusions of Theorem 4.2 are proved.
In order to verify the above theoretical analysis, we assume
d=2,h=1,k0=1,m=2,b=1.5,a=0.5. |
According to Theorem 4.1, we get k∗=1. Then from Theorem 4.2, μ2=−4.375 and β2=−0.074, which imply that the Hopf bifurcation of system (1.2) at (0,0,0,0) is nondegenerate and supercritical, a bifurcation periodic solution exists for k<k∗=1 and the bifurcating periodic solution is stable. Figure 10 shows the Hopf periodic solution occurs when k=0.999<k∗=1.
In this paper, we present a new 4D hyperchaotic system by introducing a linear controller to the first equation and second equation of the 3D Rabinovich system, respectively. If k0=0, m=1 and the fourth equation is changed to −ky, system (1.2) will be transformed to 4D hyperchaotic Rabinovich system in [18,19]. Compared with the system in [18], the new 4D system (1.1) has two nonzero equilibrium points which are symmetric about z axis when m>k0 and the dynamical characteristics are more abundant. The complex dynamical behaviors, including boundedness, dissipativity and invariance, equilibria and their stability, chaos and hyperchaos of (1.2) are investigated and analyzed. Furthermore, the existence of Hopf bifurcation, the stability and expression of the Hopf bifurcation at zero-equilibrium point are studied by using the normal form theory and symbolic computations. In order to analyze and verify the complex phenomena of the system, some numerical simulations are carried out including Lyapunov exponents, bifurcations and Poincaré maps, et al. The results show that the new 4D Rabinovich system can exhibit complex dynamical behaviors, such as periodic, chaotic and hyperchaotic. In the real practice, the hyperchaotic Rabinovich system can be applied to generate key stream for the entire encryption process in image encryption scheme [27]. In some cases, hyperchaos and chaos are usually harmful and need to be suppressed such as in biochemical oscillations [8] and flexible shaft rotating-lifting system [28]. Therefore, we will investigate hyperchaos control and chaos control in further research.
Project supported by the Doctoral Scientific Research Foundation of Hanshan Normal University (No. QD202130).
The authors declare that they have no conflicts of interest.
[1] | O. E. Rössler, An equation for hyperchaos, Phys. Lett. A, 71 (1979), 155–157. https://doi.org/10.1016/0375-9601(79)90150-6 |
[2] |
C. Xiu, R. Zhou, S. Zhao, G. Xu, Memristive hyperchaos secure communication based on sliding mode control, Nonlinear Dyn., 104 (2021), 789–805. https://doi.org/10.1007/s11071-021-06302-9 doi: 10.1007/s11071-021-06302-9
![]() |
[3] |
M. Boumaraf, F. Merazka, Secure speech coding communication using hyperchaotic key generators for AMR-WB codec, Multimedia Syst., 27 (2021), 247–269. https://doi.org/10.1007/s00530-020-00738-6 doi: 10.1007/s00530-020-00738-6
![]() |
[4] |
D. Jiang, L. Liu, L. Zhu, X. Wang, X. Rong, H. Chai, Adaptive embedding: A novel meaningful image encryption scheme based on parallel compressive sensing and slant transform, Signal Proc., 188 (2021), 108220. https://doi.org/10.1016/j.sigpro.2021.108220 doi: 10.1016/j.sigpro.2021.108220
![]() |
[5] |
P. C. Rech, Chaos and hyperchaos in a Hopfield neural network, Neurocomputing, 74 (2011), 3361–3364. https://doi.org/10.1016/j.neucom.2011.05.016 doi: 10.1016/j.neucom.2011.05.016
![]() |
[6] |
H. Li, Z. Hua, H. Bao, L. Zhu, M. Chen, B. Bao, Two-dimensional memristive hyperchaotic maps and application in secure communication, IEEE T. Ind. Electron., 68 (2020), 9931–9940. https://doi.org/10.1109/TIE.2020.3022539 doi: 10.1109/TIE.2020.3022539
![]() |
[7] |
Y. Su, X. Wang, Characteristic analysis of new four-dimensional autonomous power system and its application in color image encryption, Multimedia Syst., 28 (2022), 553–571. https://doi.org/10.1007/s00530-021-00861-y doi: 10.1007/s00530-021-00861-y
![]() |
[8] |
Y. Si, H. Liu, Y. Chen, Constructing a 3D exponential hyperchaotic map with application to PRNG, Int. J. Bifurcat. Chaos, 32 (2022), 2250095. https://doi.org/10.1142/S021812742250095X doi: 10.1142/S021812742250095X
![]() |
[9] |
A. Chen, J. Lu, J. Lü, S. Yu, Generating hyperchaotic Lü attractor via state feedback control, Phys. A, 364 (2006), 103–110. https://doi.org/10.1016/j.physa.2005.09.039 doi: 10.1016/j.physa.2005.09.039
![]() |
[10] |
Z. Yan, Controlling hyperchaos in the new hyperchaotic Chen system, Appl. Math. Comput., 168 (2005), 1239–1250. https://doi.org/10.1016/j.amc.2004.10.016 doi: 10.1016/j.amc.2004.10.016
![]() |
[11] |
T. Gao, Z. Chen, Q. Gu, Z. Yuan, A new hyper-chaos generated from generalized Lorenz system via nonlinear feedback, Chaos Soliton. Fract., 35 (2008), 390–397. https://doi.org/10.1016/j.chaos.2006.05.030 doi: 10.1016/j.chaos.2006.05.030
![]() |
[12] |
H. Wang, X. Li X, A novel hyperchaotic system with infinitely many heteroclinic orbits coined, Chaos Soliton. Fract., 106 (2018), 5–15. https://doi.org/10.1016/j.chaos.2017.10.029 doi: 10.1016/j.chaos.2017.10.029
![]() |
[13] |
N. Nguyen, T. Bui, G. Gagnon, P. Giard, G. Kaddoum, Designing a pseudorandom bit generator with a novel five-dimensional-hyperchaotic system, IEEE T. Ind. Electron., 69 (2022), 6101–6110. https://doi.org/10.1109/TIE.2021.3088330 doi: 10.1109/TIE.2021.3088330
![]() |
[14] |
S. Emiroglu, A. Akgül, Y. Adıyaman, T. E. Gümüş, Y. Uyaroglu, M. A. Yalçın, A new hyperchaotic system from T chaotic system: Dynamical analysis, circuit implementation, control and synchronization, Circuit World, 48 (2021), 265–277. https://doi.org/10.1108/CW-09-2020-0223 doi: 10.1108/CW-09-2020-0223
![]() |
[15] | A. S. Pikovski, M. I. Rabinovich, V. Y. Trakhtengerts, Onset of stochasticity in decay confinement of parametric instability, Sov. Phys. JETP, 7 (1978), 715–719. |
[16] |
J. Llibre, M. Messias, P. D. Silva, On the global dynamics of the Rabinovich system, J. Phys. A, 41 (2008), 275210. https://doi.org/10.1088/1751-8113/41/27/275210 doi: 10.1088/1751-8113/41/27/275210
![]() |
[17] | V. A. Boichenko, G. A. Leonov, V. Reitmann, Dimension theory for ordinary differential equations, Vieweg+Teubner Verlag, Wiesbaden, 2005. |
[18] |
Y. Liu, Q. Yang, G. Pang, A hyperchaotic system from the Rabinovich system, J. Comput. Appl. Math., 234 (2010), 101–113. https://doi.org/10.1016/j.cam.2009.12.008 doi: 10.1016/j.cam.2009.12.008
![]() |
[19] |
Y. Liu, Circuit implementation and finite-time synchronization of the 4D Rabinovich hyperchaotic system, Nonlinear Dyn., 67 (2012), 89–96. https://doi.org/10.1007/s11071-011-9960-2 doi: 10.1007/s11071-011-9960-2
![]() |
[20] |
Z. Wei, P. Yu, W. Zhang, M. Yao, Study of hidden attractors, multiple limit cycles from Hopf bifurcation and boundedness of motion in the generalized hyperchaotic Rabinovich system, Nonlinear Dyn., 82 (2015), 131–141. https://doi.org/10.1007/s11071-015-2144-8 doi: 10.1007/s11071-015-2144-8
![]() |
[21] |
X. Tong, Y. Liu, M. Zhang, H. Xu, Z. Wang, An image encryption scheme based on hyperchaotic Rabinovich and exponential chaos maps, Entropy, 17 (2015), 181–196. https://doi.org/10.3390/e17010181 doi: 10.3390/e17010181
![]() |
[22] |
Z. Zhang, L. Huang, A new 5D Hamiltonian conservative hyperchaotic system with four center type equilibrium points, wide range and coexisting hyperchaotic orbits, Nonlinear Dynam., 108 (2022), 637–652. https://doi.org/10.1007/s11071-021-07197-2 doi: 10.1007/s11071-021-07197-2
![]() |
[23] |
S. Yan, X. Sun, Z. Song, Y. Ren, Dynamical analysis and bifurcation mechanism of four-dimensional hyperchaotic system, Eur. Phys. J. Plus, 137 (2022), 734. https://doi.org/10.1140/epjp/s13360-022-02943-w doi: 10.1140/epjp/s13360-022-02943-w
![]() |
[24] |
Z. Li, F. Zhang, X. Zhang, Y. Zhao, A new hyperchaotic complex system and its synchronization realization, Phys. Scripta, 96 (2021), 045208. https://doi.org/10.1088/1402-4896/abdf0c doi: 10.1088/1402-4896/abdf0c
![]() |
[25] |
X. D. Edmund, K. Charles, Routh-Hurwitz criterion in the examination of eigenvalues of a system of nonlinear ordinary differential equations, Phys. Rev. A, 35 (1987), 5288–5290. https://doi.org/10.1103/PhysRevA.35.5288 doi: 10.1103/PhysRevA.35.5288
![]() |
[26] |
J. Guckenheimer, P. Holmes, Nonlinear oscillations, dynamical systems and bifurcations of vector fields, J. Appl. Mech., 51 (1984), 947. https://doi.org/10.1115/1.3167759 doi: 10.1115/1.3167759
![]() |
[27] |
B. Elizabeth, J. Gayathri, S. Subashini, A. Prakash, Hide: Hyperchaotic image encryption using DNA computing, J. Real-Time Image Proc., 19 (2022), 429–443. https://doi.org/10.1007/s11554-021-01194-9 doi: 10.1007/s11554-021-01194-9
![]() |
[28] |
S. Sajjadi, D. Baleanu, A. Jajarmi, H. Pirouz, A new adaptive synchronization and hyperchaos control of a biological snap oscillator, Chaos Soliton. Fract., 138 (2020), 109919. https://doi.org/10.1016/j.chaos.2020.109919 doi: 10.1016/j.chaos.2020.109919
![]() |
1. | Xinna Mao, Hongwei Feng, Maryam A. Al-Towailb, Hassan Saberi-Nik, Dynamical analysis and boundedness for a generalized chaotic Lorenz model, 2023, 8, 2473-6988, 19719, 10.3934/math.20231005 |
a | LE1 | LE2 | LE3 | LE4 | Dynamics |
1.08 | 0.000 | −0.075 | −4.661 | −2.538 | Periodic |
2 | 0.090 | 0.000 | −0.395 | −3.696 | Chaos |
4 | 0.325 | 0.047 | −0.000 | −6.373 | Hyperchaos |