Abstract

This work introduces a three-dimensional, highly nonlinear quadratic oscillator with no linear terms in its equations. Most of the quadratic ordinary differential equations (ODEs) such as Chen, Rossler, and Lorenz have at least one linear term in their equations. Very few quadratic systems have been introduced and all of their terms are nonlinear. Considering this point, a new quadratic system with no linear term is introduced. This oscillator is analyzed by mathematical tools such as bifurcation and Lyapunov exponent diagrams. It is revealed that this system can generate different behaviors such as limit cycle, torus, and chaos for its different parameters’ sets. Besides, the basins of attractions for this system are investigated. As a result, it is revealed that this system’s attractor is self-excited. In addition, the analog circuit of this oscillator is designed and analyzed to assess the feasibility of the system’s chaotic solution. The PSpice simulations confirm the theoretical analysis. The oscillator’s time series complexity is also investigated using sample entropy. It is revealed that this system can generate dynamics with different sample entropies by changing parameters. Finally, impulsive control is applied to the system to represent a possible solution for stabilizing the system.

1. Introduction

Chaos is a complex behavior that has been investigated in nature and mathematics [1]. It refers to the systems’ sensitivity to their initial conditions and parameters [2]. Nonlinear systems such as ordinary differential equations (ODEs) can generate chaotic behavior [3]. Therefore, they can have applications in modeling natural systems with chaotic behaviors such as neurons [4]. Besides, using coupling methods, these systems can investigate collective behaviors of neuronal networks [5]. It is good to mention that sometimes these dynamical systems are used in their nonchaotic mode to model some behaviors of natural systems like central pattern generators (neurons that make rhythms for locomotions) [6]. ODEs chaotic systems have been categorized based on their equilibrium points types and locations [7]. In this way, chaotic systems’ attractors can be divided into two groups: ones that have at least one equilibrium point in their basin of attraction (self-excited attractors) [8] and ones that have no equilibrium point in their basin of attraction (hidden attractor(s)) [9]. Besides, the systems’ equilibrium points are interesting nonlinear dynamics properties for researchers. For example, systems have been introduced and investigated with a line of equilibria [9] or just one stable equilibrium [10]. Besides, two main groups of systems can be defined based on the equations’ time dependency: autonomous systems in which no term as a function of time exists in their equations [11] and nonautonomous systems in which a term dependent to the time can be found in their equations (forced systems) [12]. Besides, forcing systems can lead nonchaotic oscillators to systems that have the capability of generating chaos. For instance, using this technique, two-dimensional systems that cannot generate chaotic time series can demonstrate chaotic attractors [13]. These nonautonomous systems are also capable of generating hidden attractors [14]. Another method to make a system chaotic is considering delays in the equations of the system. For instance, a system with just one equation can make chaos if it has a time delay in its equation [15]. On the other hand, systems with four dimensions or more can generate hyperchaotic behaviors. For instance, a four-dimensional jerk system implicated with memristors has shown the potential of demonstrating hyperchaos [16]. Besides these features, multistability refers to the existence of several attractors (at least two ones) for different initial values for a system without parameters’ changing [17]. In addition to the mentioned properties, some other features are defined for chaotic systems based on the topology and shapes of attractors [18]. Some strange attractors with different symmetries’ types were reported [19]. Besides, the systems with several wings’ attractors (multiscrolls) have grabbed researchers’ interest [20]. For instance, it is investigated how the strange multiscrolls attractor for a system can emerge and how its shape can be preserved [21]. In addition, the systems that their attractors look like known objects were also reported. For instance, chaotic systems have been introduced to look like a Persian carpet [22] or a peanut [23].

Among different ODE systems, quadratic ones are mainly focused on by some researchers interested in finding elegant systems [2]. One reason is that these systems can have simpler equations [24]. Lorenz equations, the first introduced chaotic system, are one of these classes and have just quadratic terms. Some quadratic systems were introduced whose equations’ terms are lower than that of the Lorenz equations [24]. Various dynamics of a quadratic system were studied in [25]. Most of the quadratic systems have at least one linear term in one of their equations [26]. Few systems have been introduced with no linear term in their equations [27]. Xu and Wang introduced such a system built by just nonlinear quadratic terms for the first time [28]. As another example of the pure nonlinear systems category, a multistable system can be mentioned with heterogeneous attractors [29]. Here, an oscillator with absolute nonlinear terms is introduced to generate various types of nonlinear dynamics’ behaviors such as torus and chaos.

The chaotic feasibility of nonlinear ODEs systems always has been a question. Designing analog circuits for chaotic systems has been a hot topic recently. Electrical circuits simulated with PSpice or implemented physically are tools to assess ODE systems’ chaotic behaviors. For example, an electrical circuit was introduced to regenerate the chaotic signals with a multiscroll dynamic [30]. In another instance, analog electrical circuits of a system with multistability were impacted [23]. Using memristors to model chaotic dynamics is one of the hot topics; for instance, a five-dimensional system with three linear dimensions was implicated using two memristors [1]. Besides, chaotic systems’ implication using digital circuits like field-programmable gate array (FPGA) has been carried out to assess the possibility of implicating chaotic systems. For instance, a jerk system feasibility with strange coexisting attractors was assessed with FPGA [31]. In another example, the chaotic time series of a system with coexisting attractors and strange fixed points’ curves was regenerated using FPGA [23]. One of the applications of these circuits is random number generation [32]. Other applications can be secure communications [33] and image encryption [34]. In this work, the system’s analog circuit is designed with PSpice, and the results of simulations are reported.

The complexity of chaotic systems’ signals has recently become an exciting subject for researchers [35]. For instance, the complexities of a system with hidden attractors (for time series of its different parameters’ values) were calculated and discussed [36]. Sample entropy is a feature for comparing the complexity of time series repetitively [37]. In this method, the philosophy of calculating complexity is based on the possibility of predicting the future of the signals based on their previous samples [37]. This method has some advantages in comparison with other methods of measuring complexity. For instance, it is less dependent on the length of time series than approximated entropy [37]. Here, sample entropy is used for calculating the complexities of the oscillator’s signals for different ranges of the introduced system’s parameters.

Controlling chaotic oscillators has been an interesting topic [38]. Various methods have been proposed to control the chaotic dynamics [39, 40]. Impulsive control is a method of stabilizing nonlinear systems such as the ones with infinite [41, 42] or finite delays [43], delayed neural networks [44] (that also includes exponentially stabilization [45] and fixed time control [46]), stochastic delayed systems [47], or singularly perturbed models [48]. For instance, it was used for stabilizing systems whose states are not measurable [49]. In another example, an event-based version of this method was used for controlling Chua-coupled systems [50]. This method also has been used for synchronization among nonlinear systems [51], switched complex networks [52, 53], high-dimensional Kuramoto systems [54], and fuzzy neural networks [55]. Some advanced methods of impulsive control have been introduced, for instance, versions with adaptable frequencies [56]. In this paper, an impulsive-based method for controlling the introduced pure nonlinear system is implicated as a possible solution for stabilizing its equilibrium points.

In the next section, the system’s equations whose terms are all nonlinear quadratic are presented (Section 2). Also, the oscillator’s bifurcation and Lyapunov diagrams for different parameters’ values are analyzed. Besides, the basin of attractions of the pure nonlinear oscillator is plotted and discussed. Section 3 explains the design of the introduced pure nonlinear oscillator’s analog circuit and its simulations with PSpice. The next part assesses the complexity of the oscillator’s signals for various parameter values (Section 4). Applying the impulsive control method (Section 5) to the proposed system helps to enhance its applications. The simulations’ results are concluded in the final part (Section 6).

2. The Highly Nonlinear System: Analytical and Numerical Analysis

The construction of chaotic dynamics is an unknown subject that attracted lots of attention [3, 57]. After revealing some counterexamples for the hypothesis of a relation between saddle equilibrium points and chaotic attractors [58, 59], many works have been focused on studying chaotic flows with unique properties [60, 61]. They have tried to understand the construction of chaotic attractors. Some examples are chaotic flows with different equilibrium points [62, 63] and special attractors [64]. So, a pure nonlinear chaotic flow is proposed here, and its various dynamics are investigated. The oscillator can be described by three-dimension equations that are coupled as follows:where and are the system’s variables when and are considered parameters. The system is symmetric with the change of variables . So any attractor of system (1) has a twin in reversed time and is symmetric to the origin of the main attractor. The system’s equilibrium points are as follows:

Considering these eight fixed points, the system’s Jacobian and eigenvalues are as follows:

The types of equilibria when and are shown in Table 1 (considering eigenvalues for each equilibrium).

The system’s attractors for different parameters’ values have been presented in Figure 1. Figures 1(a)1(d) demonstrate periods 1, 2, 4 and chaotic behaviors of the oscillator.

The Lyapunov exponent and bifurcation diagrams for different parameters’ set are calculated to investigate more about possible behaviors that the introduced system can present. Firstly the parameter is fixed (), and Lyapunov and bifurcation diagrams for a range of are plotted (Figure 2). Figure 2(a) demonstrates two Lyapunov exponents that have higher values than the rest. The third Lyapunov exponent’s values are always negative and have a higher absolute value than the two others. For the two Lyapunov with higher values, the system’s behavior is periodic when one is zero and the other is negative. For the situation that one of them is zero and another is positive, the system’s behavior is chaotic. When both are zero, the system’s behavior is the torus. Figure 2(a) demonstrates all of the mentioned situations; therefore, the system has the capability of having limit cycles, torus, and chaotic solutions. Figure 2(b) is the bifurcation diagram for the same range of . Period windows can be seen in Figure 2(b). In the bifurcation diagram, a period-doubling route to chaos can be observed by decreasing parameter .

In the next step, parameter is fixed, and the oscillator’s behaviors for various are investigated. Figure 3 reveals the Lyapunov exponent and bifurcation diagrams when and the ’s value changes. For better visualization, the Lyapunov exponent with the largest absolute value (its value is always negative) is not plotted in Figure 3(a). The system’s different behaviors from different limit cycles’ periods to torus and chaos can be seen based on the previously explained situations of the two larger Lyapunov exponents (Figure 3(a)). An inverse route of the period-doubling route to chaos can be observed in the bifurcation diagram by increasing (Figure 3(b)).

The basin of attractions when the oscillator’s parameters are set and are plotted for a range of initial values (Figure 4). Two surfaces each containing four equilibrium points are plotted. Studying the basin of attraction of the oscillator shows that the oscillator has only one attractor. Figures 4(a) and 4(b) show the parts of plates that and , respectively. The equilibrium points and are located at the edge of the unstable region and the basin of attraction. The type of both of them is unstable (spiral). It can be seen that some equilibrium points exist in the system’s attractor’s basin of attraction. Therefore, the system’s attractor is self-excited.

In the next section, an analog circuit of the system is implicated for the system when it is in its chaotic mode.

3. Analog System’s Circuit, Design, and PSpice Implication

The pure nonlinear system’s analog circuit in the chaotic mode is designed. Simple elements such as resistors and Op-Amps are used in its designed circuit. Its circuit’s schematic is demonstrated in Figure 5. AD633/AD as an analog device is used for multiplying variables together. The values of capacitors and resistors are tuned to compensate for the mentioned coefficient. To avoid the analog devices’ saturation, , and are considered. Therefore, the system’s equations can be rewritten as follows:

The new version of the system’s equation (Eq. (4)) is designed by analog elements (Figure 5). The values of the analog elements are as follows: and . Finally, the implicated system’s equation to simulate in PSpice can be written as follows:

The circuit simulation in PSpice when and is demonstrated in Figure 5. All elements that are used are analog. The outputs of the designed analog circuit compared to Matlab simulations are demonstrated in Figure 6.

4. The Pure Nonlinear System’s Complexity Analysis

Defining the complexity of the time series based on their predictability results in the definition of sample entropy (SaEn). Accepting this definition, SaEn is applied to estimate the complexity of the system’s time series, as reviewed in the Introduction section. SaEn tries to measure the predictability of samples of time series when the previous samples are observed. The algorithm of calculating SaEn can be read in [37]. To calculate the algorithm of SaEn, and are considered. The algorithm is applied to the oscillator’s attractors (the variable time series) for ranges of the parameters ( and ). The initial conditions are considered and the transient time parts of the time series are emitted before calculating SaEn. The results of SaEn values can be observed in Figure 7. The attractor is a fixed point in parameters that SaEn values are zero. A trend can be seen that increasing , at first, causes an increase in SaEn and then decreases it. In comparison with Figure 2, generally chaotic states of the system have more sample entropy values than periodic ones. Besides, a trend also can be observed that decreasing parameter values increases SaEn values. Comparing this trend with the bifurcation diagram reveals that the chaotic regions generally have more complexity.

5. Impulsive Control

Here, the pure nonlinear oscillator is controlled using impulsive control. In the first step, the system under impulsive control can be described as follows [6567]:

When is continuous, is continuous; is the vector of state variables; and , in general, represents linear terms of systems when contains nonlinear terms.

Definition 1. Assuming then is said to belong to the class , if(1) is continuous in and for each , exists(2) is locally Lipschitzian in

Definition 2. For , it is considered:

Definition 3. (comparison system). Let and assume thatwhere is continuous and is nondecreasing. Then the following system is the comparison system of Eq. (6):

Theorem 1. These three conditions are assumed:(1), when is continuous in for each exists. exists, (2)(3) and on , when (continuous strictly increasing function class so that ) are satisfied. Next, the global asymptotic stability for the trivial solution of the comparison system implies global asymptotic stability of impulsive system (6) trivial solution

Theorem 2. Let for all k ≥ 1. consequently, system (6) origin is global asymptotically stable if Theorem 1 conditions and the following conditions are held:(1) is nondecreasing, exists, for all (2) (3)There exists a such that is held for all , or there exists an so that for all (4) and there exists in class such that

Theorem 3. The origin of the introduced pure nonlinear chaotic system is asymptotically stable if there exists a and a differentiable at , and nonincreasing function which satisfies the following: is the largest eigenvalue assuming is a positive definite symmetric matrix and and are the smallest and the largest eigenvalues of , respectively. denote the spectral radius of for the pure nonlinear chaotic system considered so that . It is as in Theorem 1, are varying but satisfy the following:Furthermore, for a given constant ,This theorem’s proof can be seen in [66].

Remark 1. Theorem 3 also gives an estimate for the upper bound. and of impulsive intervals are given byThe introduced pure nonlinear system when and are set is considered. According to the second section, this system has eight equilibrium points. Assuming as an equilibrium point of the system, the system equations considering can be rewritten as follows:Equilibrium is considered to be stabilized. Without losing generosity for stabilizing equilibrium points of the system, the same method can be applied to the other equilibria of the pure nonlinear system. In this way, considering (6) and (12), for the equilibrium point , equations can be rewritten as follows:Considering and , then The stabilized system numerical simulations are plotted in Figure 8. Figures 8(a)8(c) demonstrate the time series of and , respectively, for the oscillator described in Eq. 12. The time series of and for the stabilized system using the impulsive controller (based on Eq. 13) are demonstrated in Figures 8(d)8(f), respectively.

6. Conclusion

Here, a pure nonlinear 3D system was presented. It was observed that the system could generate periodic, torus, and chaotic time series. Analytical analysis revealed that the oscillator has eight unstable equilibrium points for a set of parameters. The basin of attraction diagrams showed for this set of parameters that the system attractor is self-excited. The pure nonlinear system’s feasibility was investigated with an analog circuit built by simple elements like capacitors and Op-Amps. Changing parameters’ values revealed that the system could generate time series with a wide range of complexities. A possible solution to system stabilization was described by using the impulsive controller on the system. For this system, when both constants (parameters and ) were equal to zero, the system had an unbounded solution. According to the authors’ best knowledge, no pure nonlinear quadratic system has been introduced before with no constant values in its equations. Therefore, searching for such a system can be interesting for future research.

Data Availability

All the numerical simulation parameters are mentioned in the respective text part, and there are no additional data requirements for the simulation results.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was funded by the Center for Nonlinear Systems, Chennai Institute of Technology, India, vide funding number CIT/CNS/2021/RP-015.