
This paper employed the well-known Riccati transformation method to deduce a Kneser-type oscillation criterion for second-order dynamic equations. These results are considered an extension and improvement of the known Kneser results for second-order differential equations and are new for other time scales. We have included examples to highlight the significance of the results we achieved.
Citation: Taher S. Hassan, Amir Abdel Menaem, Yousef Jawarneh, Naveed Iqbal, Akbar Ali. Oscillation criterion of Kneser type for half-linear second-order dynamic equations with deviating arguments[J]. AIMS Mathematics, 2024, 9(7): 19446-19458. doi: 10.3934/math.2024947
[1] | Erxi Zhu . Time-delayed feedback control for chaotic systems with coexisting attractors. AIMS Mathematics, 2024, 9(1): 1088-1102. doi: 10.3934/math.2024053 |
[2] | Ning Cui, Junhong Li . A new 4D hyperchaotic system and its control. AIMS Mathematics, 2023, 8(1): 905-923. doi: 10.3934/math.2023044 |
[3] | 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 |
[4] | Kexin Zhang, Caihui Yu, Hongbin Wang, Xianghong Li . Multi-scale dynamics of predator-prey systems with Holling-IV functional response. AIMS Mathematics, 2024, 9(2): 3559-3575. doi: 10.3934/math.2024174 |
[5] | 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 |
[6] | 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 |
[7] | Sulasri Suddin, Fajar Adi-Kusumo, Mardiah Suci Hardianti, Gunardi . Bifurcation analysis of a diffuse large b-cell lymphoma growth model in germinal center. AIMS Mathematics, 2025, 10(5): 12631-12660. doi: 10.3934/math.2025570 |
[8] | 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 |
[9] | Hany A. Hosham, Thoraya N. Alharthi . Bifurcation and chaos in simple discontinuous systems separated by a hypersurface. AIMS Mathematics, 2024, 9(7): 17025-17038. doi: 10.3934/math.2024826 |
[10] | Na Min, Hongyang Zhang, Xiaobin Gao, Pengyu Zeng . Impacts of hunting cooperation and prey harvesting in a Leslie-Gower prey-predator system with strong Allee effect. AIMS Mathematics, 2024, 9(12): 34618-34646. doi: 10.3934/math.20241649 |
This paper employed the well-known Riccati transformation method to deduce a Kneser-type oscillation criterion for second-order dynamic equations. These results are considered an extension and improvement of the known Kneser results for second-order differential equations and are new for other time scales. We have included examples to highlight the significance of the results we achieved.
Chaos theory is a field of study that involves the mathematical analysis of complex deterministic systems exhibiting unstable aperiodic behavior. In chaotic systems, even a tiny difference in the initial conditions can lead to significantly different outcomes over time. As a way of understanding and describing the long-term behavior of a chaotic system, a chaotic attractor has at least one positive Lyapunov exponent [1]. Bifurcation diagrams can provide insights into the mechanisms behind the chaotic behavior. After transients decay, bifurcation diagrams can be constructed by recording either the Poincaré section intersections or the local peaks of the system's trajectory. By understanding the underlying mechanisms, it is possible to control and predict the behavior of the system. Amplitude control refers to the ability to regulate the magnitude of a signal, typically in the context of chaotic systems. It was shown in [2] that symmetric chaotic systems can be both multistable and also have an independent amplitude control parameter. Chaos has found applications across various fields, including economics [3], cryptography [4], and robotics [5], among others.
In their work on Chua circuits, Leonov, Kuznetsov et al. [6,7,8] introduced the classification of attractors as hidden or self-excited, and they proposed an effective procedure for localizing hidden attractors based on homotopy and numerical continuation. The classification is a crucial concept, which represents a contemporary extension of Andronov's theory of oscillations and provides a fundamental understanding of the existence and characteristics of hidden oscillatory behavior in chaotic models [9,10].
A self-excited attractor refers to an attractor whose basin of attraction intersects an open neighborhood of an equilibrium. It can be numerically found by plotting the trajectory of the system with random initial conditions near the equilibria. On the other hand, a hidden attractor is an attractor whose basin of attraction does not intersect with small neighborhoods of any equilibrium points. This implies that it cannot be numerically found using random initial conditions and demands special analytical-numerical procedures to locate its basin of attraction. In practice, there is a need to control the hidden attractors, because some structures (such as bridges or airplane wings) can display catastrophic responses under disturbances [11,12]. According to [13], a saddle focus of index 2 is defined by an equilibrium with three eigenvalues: one real γ and a pair of complex conjugates σ±jω, where γ<0 and σ>0. Shilnikov homoclinic and heteroclinic theorems provide sufficient conditions for the existence of self-excited chaotic attractors in 3D systems with unstable hyperbolic equilibria. Saddle-foci with an index of 2 are frequently encountered in such systems. As noted by Zhou and Chen [14], the chaotic attractors of Lorenz, Rössler, Chua, Chen, and many other systems are in the sense of Shilnikov, but not all chaotic attractors are of Shilnikov sense. Beyond the scope of Shilnikov theorems, many chaotic flows with various types of non-hyperbolic equilibria were introduced [15,16,17]. Li and Hai [18] introduced a three-dimensional robust chaotic system with two non-hyperbolic equilibrium points, exhibiting features of amplitude and position modulation. It is shown that the presence of linear terms is not essential in producing complex behaviors [19].
In a review paper, Jafari et al. [20] noted that dissipative hidden chaotic flows can be categorized into three types: those with a stable equilibrium [21,22], those without equilibrium [23], and those with an infinite number of equilibria [24]. Molaie et al. [25] proposed twenty-three chaotic flows with a single stable equilibrium point. Jafari and Sprott [26] proposed nine chaotic flows with hidden attractors that have a line equilibrium. Kumarasamy et al. [27] investigated how saddle-node bifurcations contribute to the emergence of hidden attractors, both in systems with stable equilibrium points and those without equilibrium points. Antimonotonicity refers to the intertwined creation and destruction of periodic orbits in nonlinear dynamical systems, leading to inevitable reversals of bifurcation cascades [28]. For a family of jerk systems with a hyperbolic tangent function, Li et al. [29] illustrated that antimonotonicity is important in elucidating the generation mechanism of hidden attractors with a stable equilibrium. Recently, Dong et al. [30] investigated a chaotic system that has no equilibrium through numerical simulations and circuit implementations, and discussed engineering applications such as synchronization, offset boosting control, and image encryption.
A chameleon system is a chaotic system that, as parameter values vary, can exhibit various types of chaotic attractors, including self-excited attractors and the three types of hidden attractors. This feature provides flexibility in applications needing different types of chaotic signals. In a review paper, Sprott [31] categorized chaotic flows from a quadratic family based on equilibrium points, in which both self-excited and hidden attractors are involved. Jafari et al. [32] introduced a chameleon system with quadratic nonlinearities, contributing to an enhanced understanding of chaotic attractors, especially those that are hidden. Recall that a hyperchaotic system is defined as a dynamical system that exhibits chaotic behavior, characterized by the presence of multiple positive Lyapunov exponents. Motivated by [32], Rajagopal et al. [33] presented a hyperchaotic chameleon system, contributing to a deeper understanding of hidden chaotic flows of higher dimensions. Natiq et al. [34] presented a chameleon system with a trigonometric nonlinear controller. Cang et al. [35] investigated the hidden and self-excited coexisting attractors in the Yang system [36]. Signing et al. [37] introduced a chameleon cryptosystem that combines chameleon chaotic systems with dynamic DNA coding to enhance image encryption security through the use of chaotic sequences. Fan et al. [38] investigated the two-parameter space bifurcation of the Wei chameleon system. Recently, Tiwari et al. [39] introduced a class of quadratic chameleon systems that allow for the coexistence of hidden chaotic attractors and conservative tori.
A Hopf bifurcation is a local bifurcation where an equilibrium of a dynamical system changes stability as a parameter is varied. This phenomenon occurs when the equilibrium possesses a pair of purely imaginary eigenvalues, accompanied by specific transversality conditions and an absence of zero eigenvalues. The bifurcation can be classified as either supercritical or subcritical, resulting in either a stable or an unstable limit cycle within an invariant two-dimensional manifold, respectively. The latter situation presents a potential risk, as stable, large-amplitude limit cycles can coexist with the stable equilibrium point [40]. Stankevich et al. [41] and Zhao et al.[42] studied the dynamics of the Chua systems, establishing a connection between subcritical Hopf bifurcations of equilibria and the emergence of hidden attractors. Using a family of quadratic jerk systems, Liu et al. [43] also showed that subcritical Hopf bifurcations are closely related to the emergence of hidden attractors. Yang et al. [44] proposed a 7D hyperchaotic system with five positive exponents and investigated Hopf bifurcation to understand how complex dynamics emerge in this system. Li and Chen [45] demonstrated the complex neuromorphic behaviors of a novel current-controlled Chua corsage memristor-based neuron circuit through Hopf bifurcation and edge-of-chaos analysis. Liu et al. [46] explored a 4D hyperchaotic system through Hopf bifurcation analysis and the linear control method.
The design of a chaotic system with independent non-bifurcation parameters faces a significant challenge due to the broadband random-like nature of chaos [47,48]. Signal polarity refers to whether a signal oscillates around zero (bipolar) or around a non-zero reference level (unipolar). Offset boosting involves adjusting the average value of a state variable by introducing a constant term into the system model, resulting in signal offset and enabling control of signal polarity. Li et al. [49] proposed an offset boosting-based strategy for the identification of multistability in dynamical systems, employing nonbifurcation operations for diagnosing multistability. In chaotic and hyperchaotic systems [50,51], offset boosting controls have been implemented, enabling direct offset regulation through the use of a single constant.
The paper is organized as follows. In Section 2, we introduce a novel chameleon system and conduct a detailed study on its symmetry, dissipation, stability analysis, Hopf bifurcation, and chaotic dynamics with μ varying. In Section 3, we perform a bifurcation analysis for scenarios 1 to 4 with μ=0. We investigate the system's behavior under various parameter values, focusing on the emergence of different dynamical phenomena. In Section 4, we investigate hidden chaotic bistability with offset boosting through the variation of μ. In the final section, we synthesize our key findings and discuss their implications.
Consider the quadratic system
{dxdt=y,dydt=−x+y(z+μ),dzdt=f(x,y,z), | (2.1) |
where
f(x,y,z)=a−bz+cy2−0.7x2−3xy. |
It is symmetric with respect to 180-degree rotation about the z-axis. This means that the system is invariant with respect to the transformation (x,y,z)→(−x,−y,z), corresponding to a 180-degree rotation about the z-axis.
Let h=(h1,h2,h3) be the vector field of system (2.1). Thus the divergence function is
∇⋅V=h1x+h2y+h3z=−b+(z+μ). |
It is clear that system (2.1) is not always dissipative as the divergence depends on z. However, in specific regions of the phase space where the condition ∇⋅V<0 holds on average, we have a dissipative flow. Due to the symmetry, system (2.1) can have a symmetrical attractor or a symmetric pair of coexisting attractors.
If a=b=0, system (2.1) has a line equilibrium along the z-axis. This means that every point (0,0,z0) on the axis is an equilibrium point.
If a≠0 and b=0, we have f(0,0,z)=a≠0. Hence, there is no equilibrium point. When x=y=0, we have dxdt=dydt=0 and dzdt=a≠0. This indicates that the z-axis is invariant under the system's flow. The direction of the flow along the z-axis is determined by the sign of a: if a>0, the flow is in the forward direction; if a<0, the flow is in the backward direction.
If b≠0, there exists a unique equilibrium P:(0,0,a/b). The Jacobian matrix of system (2.1) at P is
A=[010−1μ+ab000−b]. |
The corresponding characteristic polynomial is
φA(λ)=λ3+p1λ2+p2λ+p3=(λ+b)(λ2−(μ+ab)λ+1). | (2.2) |
Proposition 2.1. Consider system (2.1) with b≠0. It is locally asymptotically stable around P if
b>0anda<−bμ. | (2.3) |
Proof. From (2.2) and Routh-Hurwitz criterion, the equilibrium P of system (2.1) is asymptotically stable if b>0 and μ+ab<0. Such constraints are equivalent to the inequalities stated in (2.3). Thus the conclusion holds.
Corollary 2.1. Consider system (2.1) with b>0anda<−bμ, and let
Δ=(μ+ab)2−4. |
We have the following conclusions.
(i) If Δ<0, then the equilibrium P is a stable node-focus.
(ii) If Δ≥0, then the equilibrium P is a stable node.
Proof. (ⅰ) Assume that Δ<0. From (2.2), we know that the characteristic polynomial has a negative root and a pair of complex conjugate roots with a negative real part. Thus the equilibrium P is a stable node-focus.
(ⅱ) Assume that Δ≥0. From (2.2), we know that the characteristic polynomial has three negative roots, counted with multiplicity. Therefore, the equilibrium P is a stable node.
Assume a,b≠0 and μ=0. Then system (2.1) can assume all types of hyperbolic equilibrium points at P:(0,0,a/b).
Corollary 2.2. Consider system (2.1) with b≠0 and μ=0. The following statements hold.
(i) It is locally asymptotically stable around P if a<0 and b>0. In this case, if a>−2b, then the equilibrium P is a stable node-focus; if a≤−2b, then the equilibrium is a stable node.
(ii) If a>0 and b>0, the equilibrium P is unstable. In this case, if a<2b, the equilibrium P is a saddle-focus of index 2; if a≥2b, the equilibrium P is a saddle of index 2.
(iii) If b<0, the equilibrium P is unstable. In this case, if 0<a<−2b, the equilibrium P is a saddle-focus of index 1; if a≥−2b, the equilibrium P is a saddle of index 1; if 2b<a<0, the equilibrium P is an unstable node-focus; and if a≤2b, the equilibrium P is an unstable node.
Proof. (ⅰ) By setting μ=0, the statements follow from Proposition 2.1 and Corollary 2.1.
(ⅱ) Setting μ=0, then the characteristic polynomial (2.2) becomes
φA(λ)=(λ+b)(λ2−abλ+1). | (2.4) |
Since a>0 and b>0, the quadratic factor of (2.4) can have a root with a positive real part. Thus the equilibrium P is unstable.
Assume that a,b>0. If a<2b, the cubic polynomial (2.4) has a negative root and a pair of complex conjugate roots with a positive real part. Therefore, the equilibrium P is a saddle-focus of index 2. If a≥2b, the cubic polynomial (2.4) has a negative root and two positive roots. Therefore, the equilibrium P is a saddle of index 2.
(ⅲ) If b<0, the characteristic polynomial (2.4) has a positive root λ=−b. Thus the equilibrium P is unstable. The rest of the proof for classifying the equilibrium is similar to that of case (ⅱ).
When a=0, system (2.1) becomes
{dxdt=y,dydt=−x+y(z+μ),dzdt=−bz+cy2−0.7x2−3xy. | (2.5) |
It has a unique equilibrium point at the origin. Here we assume that b>0.
Proposition 2.2. When μ=0, a Hopf bifurcation occurs at the origin for system (2.5).
Proof. The characteristic polynomial of this system is
φA(λ)=λ3+p1(μ)λ2+p2(μ)λ+p3(μ), |
where
p1=−μ+b,p2=−bμ+1,p3=b. |
Note that
(p1p2−p3)(μ)=bμ2−(1+b2)μ. |
So
(p1p2−p3)(0)=0, | (2.6) |
(p1p2−p3)′(0)=−(1+b2)<0. | (2.7) |
Thus, a Hopf bifurcation occurs from the origin at μ=0.
When μ=0, system (2.5) becomes
{dxdt=y,dydt=−x+F2,dzdt=−bz+F3, | (2.8) |
where
F2=yz,F3=cy2−0.7x2−3xy. |
According to [52], the first Lyapunov number a1 at the Hopf bifurcation can be obtained as
a1=−216λ3F2yz(F3xx+F3yy)+F2yz16(4ω2+λ23)(λ3(F3xx−F3yy)+4ωF3xy), | (2.9) |
where λ3=−b and ω=−1.
By a simple computation, we have
a1=(30c−7)b2+60b+(80c−56)80b(b2+4). | (2.10) |
Let us consider the following system:
{dxdt=y,dydt=−x+y(z+μ),dzdt=−bz+y2−0.7x2−3xy, | (2.11) |
where b>0.
Corollary 2.3. When μ=0, a subcritical Hopf bifurcation occurs at the origin for system (2.11), which generates one unstable limit cycle when μ is negative but not too small. For μ=0, the origin is unstable.
Proof. Note that b>0 and system (2.11) is a special case of (2.5) with c=1.
By setting c=1 in (2.10), we have
a1=23b2+60b+2480b(b2+4)>0, | (2.12) |
where a1 is the first Lyapunov number. According to Proposition 2.2, conditions (2.6) and (2.7), and the sign of (2.12), we have the conclusions.
The following results give us insights into the stability of the system (2.5).
Proposition 2.3. Consider system (2.5) with μ≠0 and b>0. The stability properties of the origin are characterized as follows:
(i) If −2<μ<0, the origin is a stable node-focus.
(ii) If 0<μ<2, the origin is a saddle-focus of index 2.
(iii) If μ≤−2, the origin is a stable node.
(iv) If μ≥2, the origin is a saddle of index 2.
Proof. The Jacobian matrix of the system at the origin is
A=[010−1μ000−b]. |
The corresponding characteristic polynomial is
φA(λ)=(λ+b)(λ2−μλ+1). | (2.13) |
By analyzing of the roots of (2.13), we can draw the four conclusions.
Consider the following three systems:
{dxdt=y,dydt=−x+y(z+μ),dzdt=y2−0.7x2−3xy, | (2.14) |
{dxdt=y,dydt=−x+y(z+μ),dzdt=0.001+y2−0.7x2−3xy, | (2.15) |
{dxdt=y,dydt=−x+y(z+μ),dzdt=−0.01z+y2−0.7x2−3xy, | (2.16) |
where μ∈[−0.08,0.02]. These systems can be obtained from system (2.1) by setting (a,b,c)=(0,0,1), (a,b,c)=(0.001,0,1), and (a,b,c)=(0,0.01,1), respectively. Figure 1 shows the bifurcation diagrams for these systems with respect to the parameter μ.
Let (L(k)1,L(k)2,L(k)3),k=1,2,3, be the Lyapunov exponent spectra of systems (2.14), (2.15), and (2.16), respectively. For each k, the exponents satisfy L(k)1>0, L(k)2=0, and L(k)3<0. The variations of L(k)1 and L(k)3 with respect to the parameter μ for 1≤k≤3 are shown in Figure 2.
Based on the two figures above and the nature of equilibrium points, we classify the dynamical behaviors of these three systems into four scenarios:
Scenario 1: Hidden chaotic dynamics with a line equilibrium Note that system (2.14) has a line equilibrium x=y=0. Thus, from panel (a) of Figure 1, we see that system (2.14) exhibits a symmetric pair of hidden chaotic attractors for μ∈[−0.08,0.02], coexisting with a line equilibrium.
Scenario 2: Hidden chaotic dynamics without equilibrium It is clear that there is no equilibrium in system (2.15). Thus, from panel (b) of Figure 1, we see that system (2.15) exhibits a symmetric pair of hidden chaotic attractors without an equilibrium point for μ∈[−0.08,0.02].
Scenario 3: Hidden chaotic dynamics with a stable equilibrium By Proposition 2.3, we know that system (2.16) has a stable equilibrium at the origin when μ∈[−0.08,0). Thus, from panel (c) of Figure 1, we see that system (2.16) exhibits a symmetrical hidden chaotic attractor for μ∈[−0.08,0).
Scenario 4: Self-excited chaotic dynamics with an unstable equilibrium By Corollary 2.3 and Proposition 2.3, we know that system (2.16) has an unstable equilibrium at the origin when μ∈[0,0.02]. Thus, from panel (c) of Figure 1, we see that system (2.16) exhibits a symmetrical self-excited chaotic attractor for μ∈[0,0.02].
Using two fixed initial conditions and varying parameters c and a, we demonstrate the system's period-doubling routes to chaos when μ=0. The analysis reveals that system (2.1) can exhibit either self-excited or hidden chaotic attractors under different parameter settings.
Let
f1(x,y)=y2−0.7x2−3xy, | (3.1) |
f2(x,y)=1.008y2−0.7x2−3xy. | (3.2) |
When μ=0 and f(x,y,z)=f1(x,y), system (2.1) becomes
{dxdt=y,dydt=−x+yz,dzdt=y2−0.7x2−3xy. | (3.3) |
It has a symmetric pair of hidden chaotic attractors with a line equilibrium, which is shown in Figure 3 by initiating with (x,y,z)=(±0.8,±0.8,0). Figure 4 shows various projections of the coexisting hidden chaotic attractors. The corresponding Poincaré maps for the two attractors are shown in Figure 5. The set of dense points covering a certain region on the plane y=0 shows that each attractor is chaotic. When the initial conditions are set to (0.8,0.8,0), system (3.3) has Lyapunov exponents of (L1,L2,L3)=(0.0548,0,−0.6914) and a Kaplan-Yorke dimension of DKY=2−L1/L3≈2.079. The sum of the three exponents L1+L2+L3=−0.6366 signifies that the system exhibits dissipative behavior, characterized by a chaotic attractor.
When μ=0 and f(x,y,z)=f2(x,y), system (2.1) becomes
{dxdt=y,dydt=−x+yz,dzdt=1.008y2−0.7x2−3xy. | (3.4) |
It has a symmetric hidden chaotic attractor with respect to the z-axis, which is visualized in Figure 6 by starting from the initial conditions (x,y,z)=(0.8,0.8,0). Moreover, the line equilibrium z-axis is denoted by a green line. The Poincaré maps of the attractor are shown in Figure 7, in which forward and backward directions are considered, respectively. When the initial conditions are set to (0.8,0.8,0), system (3.4) has Lyapunov exponents of (L1,L2,L3)=(0.1179,0,−0.7641) and a Kaplan-Yorke dimension of DKY=2−L1/L3≈2.154.
Consider the following system:
{dxdt=y,dydt=−x+yz,dzdt=cy2−0.7x2−3xy, | (3.5) |
where c∈[0.8,1.008] is the coefficient of the nonlinear term y2. It contains systems (3.3) and (3.4) as special cases.
For the two initial conditions (±0.8,±0.8,0), the bifurcation diagram is presented in Figure 8(a) for c∈[0.94,1.008]. It shows the development of the hidden chaotic motions (with a line equilibrium) via period-doubling cascades. The periodic orbit undergoes a symmetry-breaking bifurcation around c≈0.967, transitioning into a saddle orbit, and giving rise to a pair of asymmetric stable periodic orbits. For each initial data, the first period-doubling bifurcation occurs at c≈0.991. The period and amplitude of the two asymmetric attractors grow with increasing c until they touch the saddle orbit and merge into a broader symmetrical hidden chaotic attractor at c≈1.0027 in an attractor merging crisis. Recall that an attractor merging crisis refers to two or more chaotic attractors merging to form a single attractor as the critical parameter value is passed [53]. For the initial data (x(0),y(0),z(0))=(0.8,0.8,0), the variations of Lyapunov exponents and their sum with respect to the parameter c are shown in Figure 8(b). The Lyapunov exponent spectrum depicted in Figure 8(b) aligns with the bifurcation diagram illustrated in Figure 8(a). In Figure 8(b), we also plotted the variation of Lsum=L1+L2+L3 with respect to c. For each value of c, the sum Lsum is negative, indicating that the system is dissipative. Hence, the system is always dissipative for c within the interval [0.94,1.008].
Letting b=μ=0 in system (2.1), we have the following system:
{dxdt=y,dydt=−x+yz,dzdt=a+(y2−0.7x2−3xy). | (3.6) |
When a≠0, it has no equilibrium but does exhibit an invariant straight line along the z-axis. Recall that when a=0, the system can have a symmetric pair of hidden chaotic attractors with a line equilibrium (z-axis).
For the two initial conditions (±0.8,±0.8,0), the bifurcation diagram is presented in Figure 9(a) for a∈[−0.02,0.004]. It shows the development of the hidden chaotic motions (with no equilibrium) via period-doubling cascades. The left panel of Figure 9 with the change of a is qualitative similar to Figure 8(a) with the change of c. For the initial data (x(0),y(0),z(0))=(0.8,0.8,0), the variation of Lyapunov exponents with respect to the parameter a is shown in Figure 9(b). Note that at a≈0.0013, a symmetric pair of hidden chaotic attractors combine into a symmetrical hidden chaotic attractor. At a=0.004, there exists a symmetrical hidden chaotic attractor with a Lyapunov exponent spectrum of (L1,L2,L3)=(0.1162,0,−0.7621) and a Kaplan-Yorke dimension of 2.152.
Letting b=0.01,μ=0 in system (2.1), we have the following system:
{dxdt=y,dydt=−x+yz,dzdt=(a−0.01z)+(y2−0.7x2−3xy), | (3.7) |
which has a single equilibrium at P:(0,0,100a). According to Corollary 2, the equilibrium P is locally asymptotically stable for a<0.
For the two initial conditions (±0.8,±0.8,0), the bifurcation diagram is presented in Figure 10(a) for a∈[−0.03,0]. It shows the development of chaotic motions via period-doubling cascades. For the initial data (x(0),y(0),z(0))=(0.8,0.8,0), the variation of Lyapunov exponents with respect to the parameter a is show in Figure 10(b). From Figure 10, there is a narrow range near a=0 with a<0 that corresponds to hidden chaotic dynamics with a stable node-focus.
For two initial conditions (±0.8,±0.8,0), the bifurcation diagram is presented in Figure 11 for a∈[0,2×10(−4)]. When a=0, from (2.12), we know that the first Lyapunov number a1 of system (3.7) is positive, thus the origin is unstable. Recalling from Corollary 2.2 that the equilibrium P is a saddle-focus of index 2 for a∈(0,2×10(−4)], hence, we have that the chaotic range [0,2×10(−4)] corresponds to self-excited dynamics.
With two initial conditions (∓0.8,∓0.8,−μ), the variations of the average values of the state variables in systems (2.14) and (2.15) with respect to μ are shown in Figures 12 and 13, respectively. Since the difference between the two system is caused by a small constant of 0.001, the figures are almost the same. For both systems, the average value of z decreases smoothly as the parameter μ increases.
Let (L(k)1,L(k)2,L(k)3),k=1,2, be the Lyapunov exponent spectra of systems (2.14) and (2.15), respectively. For each k, we have L(k)1>0, L(k)2=0, and L(k)3<0. For 1≤k≤2, the variations of L(k)1 and L(k)3 with respect to the parameter μ are shown in Figure 14.
Figure 14 indicates that as the initial conditions vary with the parameter μ, both of these systems exhibit chaotic dynamics. For system (2.14), there exists a symmetric pair of hidden chaotic attractors accompanied by an infinite number of equilibrium points. As for system (2.15), there also exists a symmetric pair of hidden chaotic attractors, but there are no equilibrium points. As shown in panel (a) of Figure 14, the largest Lyapunov exponent L(2)1 (red curve) of system (2.15) consistently exceeds the largest Lyapunov exponent L(1)1 (blue curve) of system (2.14), indicating that system (2.15) demonstrates stronger chaotic dynamics than system (2.14).
For systems (2.14) and (2.15) with μ=−5,5, the coexistence of hidden chaotic attractors is illustrated in Figures 15 and 16, respectively.
In this paper, we have conducted a comprehensive investigation of a novel chameleon system, examining its symmetry, dissipation, stability analysis, Hopf bifurcation, and various chaotic dynamics. The parameter μ serves as both the Hopf bifurcation parameter and the offset boosting parameter. The other parameters are also critical control parameters, with variation of which, the system undergoes period-doubling bifurcations leading to various types of chaos, including hidden chaos and self-excited chaos, with hidden chaos being predominant. By offset boosting the variable z, we also investigated the coexistence of hidden chaotic attractors. The results of this study enrich our understanding of the dynamical behavior of chameleon systems, providing important insights for research in related fields such as information security, cryptography, and nonlinear dynamics. In the future, it is hoped that there will be studies about chameleon systems subjected to small random perturbations.
Jie Liu: Writing-first draft. Bo Sang: Conceptualization, Writing-second and third drafts. Lihua Fan: Software. Chun Wang: Software. Xueqing Liu: Software. Ning Wang: Conceptualization, Software. Irfan Ahmad: Conceptualization. All authors have read and approved the final version of the manuscript for publication.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
This research is supported by the Shandong Provincial Natural Science Foundation (No. ZR2018MA025).
The authors declare that they have no competing interests in this paper.
[1] | S. Hilger, Analysis on measure chains–A unified approach to continuous and discrete calculus, Results Math., 18 (1990), 18–56. https://doi.org/10.1007/BF03323153 |
[2] | V. Kac, P. Chueng, Quantum calculus, New York: Springer, 2002. https://doi.org/10.1007/978-1-4613-0071-7 |
[3] | M. Bohner, A. Peterson, Dynamic equations on time scales: an introduction with applications, Boston: Birkhäuser, 2001. https://doi.org/10.1007/978-1-4612-0201-1 |
[4] | M. Bohner, A. Peterson, Advances in dynamic equations on time scales, Boston: Birkhäuser, 2003. https://doi.org/10.1007/978-0-8176-8230-9 |
[5] | R. P. Agarwal, S. L. Shieh, C. C. Yeh, Oscillation criteria for second-order retarded differential equations, Math. Comput. Modell., 26 (1997), 1–11. https://doi.org/10.1016/S0895-7177(97)00141-6 |
[6] | L. Erbe, T. S. Hassan, A. Peterson, S. H. Saker, Oscillation criteria for sublinear half-linear delay dynamic equations on time scales, Int. J. Differ. Equ., 3 (2008), 227–245. |
[7] |
B. Baculikova, Oscillation of second-order nonlinear noncanonical differential equations with deviating argument, Appl. Math. Lett., 91 (2019), 68–75. https://doi.org/10.1016/j.aml.2018.11.021 doi: 10.1016/j.aml.2018.11.021
![]() |
[8] |
O. Bazighifan, E. M. El-Nabulsi, Different techniques for studying oscillatory behavior of solution of differential equations, Rocky Mountain J. Math., 51 (2021), 77–86. https://doi.org/10.1216/rmj.2021.51.77 doi: 10.1216/rmj.2021.51.77
![]() |
[9] |
T. S. Hassan, Oscillation criteria for half-linear dynamic equations on time scales, J. Math. Anal. Appl., 345 (2008), 176–185. https://doi.org/10.1016/j.jmaa.2008.04.019 doi: 10.1016/j.jmaa.2008.04.019
![]() |
[10] |
S. H. Saker, Oscillation criteria of second-order half-linear dynamic equations on time scales, J. Comput. Appl. Math., 177 (2005), 375–387. https://doi.org/10.1016/j.cam.2004.09.028 doi: 10.1016/j.cam.2004.09.028
![]() |
[11] |
B. Almarri, A. H. Ali, A. M. Lopes, O. Bazighifan, Nonlinear differential equations with distributed delay: some new oscillatory solutions, Mathematics, 10 (2022), 995. https://doi.org/10.3390/math10060995 doi: 10.3390/math10060995
![]() |
[12] | I. Jadlovská, Iterative oscillation results for second-order differential equations with advanced argument, Electron. J. Differ. Eq., 2017 (2017), 1–11. |
[13] |
M. Bohner, K. S. Vidhyaa, E. Thandapani, Oscillation of noncanonical second-order advanced differential equations via canonical transform, Constructive Math. Anal., 5 (2022), 7–13. https://doi.org/10.33205/cma.1055356 doi: 10.33205/cma.1055356
![]() |
[14] |
G. E. Chatzarakis, J. Džurina, I. Jadlovská, New oscillation criteria for second-order half-linear advanced differential equations, Appl. Math. Comput., 347 (2019), 404–416. https://doi.org/10.1016/j.amc.2018.10.091 doi: 10.1016/j.amc.2018.10.091
![]() |
[15] |
T. S. Hassan, C. Cesarano, R. A. El-Nabulsi, W. Anukool, Improved Hille-type oscillation criteria for second-order quasilinear dynamic equations, Mathematics, 10 (2022), 3675. https://doi.org/10.3390/math10193675 doi: 10.3390/math10193675
![]() |
[16] |
G. E. Chatzarakis, O. Moaaz, T. Li, B. Qaraad, Some oscillation theorems for nonlinear second-order differential equations with an advanced argument, Adv. Differ. Equ., 2020 (2020), 160. https://doi.org/10.1186/s13662-020-02626-9 doi: 10.1186/s13662-020-02626-9
![]() |
[17] |
R. P. Agarwal, C. Zhang, T. Li, New Kamenev-type oscillation criteria for second-order nonlinear advanced dynamic equations, Appl. Math. Comput., 225 (2013), 822–828. https://doi.org/10.1016/j.amc.2013.09.072 doi: 10.1016/j.amc.2013.09.072
![]() |
[18] |
S. Frassu, G. Viglialoro, Boundedness in a chemotaxis system with consumed chemoattractant and produced chemorepellent, Nonlinear Anal., 213 (2021), 112505. https://doi.org/10.1016/j.na.2021.112505 doi: 10.1016/j.na.2021.112505
![]() |
[19] |
T. Li, G. Viglialoro, Boundedness for a nonlocal reaction chemotaxis model even in the attraction-dominated regime, Differ. Integral Equ., 34 (2021), 315–336. https://doi.org/10.57262/die034-0506-315 doi: 10.57262/die034-0506-315
![]() |
[20] |
R. P. Agarwal, M. Bohner, T. Li, Oscillatory behavior of second-order half-linear damped dynamic equations, Appl. Math. Comput., 254 (2015), 408–418. https://doi.org/10.1016/j.amc.2014.12.091 doi: 10.1016/j.amc.2014.12.091
![]() |
[21] | M. Bohner, T. S. Hassan, T. Li, Fite-Hille-Wintner-type oscillation criteria for second-order half-linear dynamic equations with deviating arguments, Indag. Math., 29 (2018), 548–560. https://doi.org/10.1016/j.indag.2017.10.006 |
[22] |
M. Bohner, T. Li, Oscillation of second-order p-Laplace dynamic equations with a nonpositive neutral coefficient, Appl. Math. Lett., 37 (2014), 72–76. https://doi.org/10.1016/j.aml.2014.05.012 doi: 10.1016/j.aml.2014.05.012
![]() |
[23] |
T. Li, N. Pintus, G. Viglialoro, Properties of solutions to porous medium problems with different sources and boundary conditions, Z. Angew. Math. Phys., 70 (2019), 86. https://doi.org/10.1007/s00033-019-1130-2 doi: 10.1007/s00033-019-1130-2
![]() |
[24] |
T. S. Hassan, M. Bohner, I. L. Florentina, A. A. Menaem, M. B. Mesmouli, New criteria of oscillation for linear Sturm-Liouville delay noncanonical dynamic equations, Mathematics, 11 (2023), 4850. https://doi.org/10.3390/math11234850 doi: 10.3390/math11234850
![]() |
[25] |
G. V. Demidenko, I. I. Matveeva, Asymptotic stability of solutions to a class of second-order delay differential equations, Mathematics, 9 (2021), 1847. https://doi.org/10.3390/math9161847 doi: 10.3390/math9161847
![]() |
[26] | A. Kneser, Untersuchungen über die reellen nullstellen der integrale linearer differentialgleichungen, Math. Ann., 42 (1893), 409–435. |
[27] |
I. Jadlovská, J. Džurina, Kneser-type oscillation criteria for second-order half-linear delay differential equations, Appl. Math. Comput., 380 (2020), 125289. https://doi.org/10.1016/j.amc.2020.125289 doi: 10.1016/j.amc.2020.125289
![]() |
[28] | O. Došlý, P. Řehák, Half-linear differential equations, Elsevier, 2005. |
[29] | R. P. Agarwal, S. R. Grace, D. O'Regan, Oscillation theory for second-order linear, half-linear, superlinear and sublinear dynamic equations, Springer Dordrecht, 2002. https://doi.org/10.1007/978-94-017-2515-6 |
[30] | T. S. Hassan, Y. Sun, A. A. Menaem, Improved oscillation results for functional nonlinear dynamic equations of second-order, Mathematics, 8 (2020), 1897. https://doi.org/10.3390/math8111897 |