Comments on "Coexistence of hidden chaotic attractors in a novel no-equilibrium system" (Nonlinear Dyn, doi:10.1007/s11071-016-3170-x).

In this comment, an enhancement of issue published in the paper "Coexistence of hidden chaotic attractors in a novel no-equilibrium system" (Nonlinear Dyn, doi:10.1007/s11071-016-3170-x) is addressed. We have shown that the proposed novel autonomous chaotic system can be extended to its fractional-order version where hidden attractors as well as other dynamical properties of the new no-equilibrium system can be observed. A created MATLAB function for the new fractional-order no-equilibrium system is also presented.

Introduction and motivation

Investigation of chaotic systems has a long tradition. Basically, there are two kinds of attractors observed in chaotic systems: self-excited attractor and hidden attractor. It is well-known that hidden attractors cannot be localized by using the standard computational procedure and therefore it is difficult to predict the existence of them in system. However, hidden attractors appear in numerous dynamical systems [6, 9–11]. Moreover, there is a significant interest in studying hidden attractors because they play an important role in both theoretical problems and practical engineering applications [10]. For more details see also survey paper [5].

On the other hand, the subject of fractional calculus has gained considerable popularity and importance mainly during the past three decades [8, 13].

Fractional (order) calculus, also non-integer order calculus, is known since the classical calculus with the first written reference dated September 1695 in the correspondence between Leibniz and L’Hospital. Nowadays, the fractional calculus has a wide area of applications in various fields [4, 8, 13] as well as in the chaotic systems theory and applications [12].

Based on above consideration, it is possible to combine the fractional calculus into description of the novel chaotic systems where the hidden attractors appeared.

This comment is organized as follows. In Sect. 1, the introduction to the problem and motivation is briefly discussed. Section 2 is focused on fractional calculus fundamentals. Section 3 brings a new fractional-order chaotic system. In Sect. 4, the simulation results and discussion are described. Section 5 concludes this note with some additional remarks for further research investigation.

Preliminaries

Definition of fractional-order operator

Fractional calculus is a generalization of integration and differentiation to non-integer (fractional) order fundamental operator , where a and b are the bounds of the operation. The standard notation for denoting the left-sided fractional-order integro-differential operator of a function f(t) defined in the interval [a, b] is , with .

There exist many definitions for fractional operator (fractional order integrals and derivatives), but in this article we will focus only on the Grünwald–Letnikov definition (GLD), which is given as follows [8, 13]:where is the floor function, i.e., the greatest integer smaller than z, andare the binomial coefficients for This form of definition is very helpful for obtaining numerical solutions of fractional differential equations. However, we must be careful with the initial conditions setting because it could be a complex problem in some systems.

Numerical solution of fractional differential equations

For numerical calculation of fractional-order derivatives, we can use the relation (3) derived from the GLD (1). This approach is based on the fact that for a wide class of functions, the most of definitions are equivalent. The relation for the explicit numerical approximation of q-th derivative at the points has the following form [13]:where is the “memory length”, is the time step of calculation (definition (3) is valid only as h tends toward 0 and that the accuracy of the simulation depends on the value of h) and are binomial coefficients. For their calculation we can use the following expression:Then, general numerical solution of the fractional differential equationcan be expressed asFor the memory term expressed by the sum, a “short memory” principle can be used for various .

An evaluation of the short memory effect and convergence relation of the error between short and long memory were clearly described and also proved in [13].

Model of the new fractional-order system

In paper [11], the simple case of the new chaotic system without equilibrium was described in the following form:in which a is the positive parameter (). The dynamical properties of the new system (6) as for instance equilibrium points (there is no equilibrium in (6)), dissipativity, chaotic behavior for various values of parameter a as for example Poincaré map, bifurcation diagram, Lyapunov exponents, and chaotic attractors have been investigated in paper [11]. Moreover, an electronic circuit implementation of the system (6), which included 15 resistors (), 3 capacitors (10 nF), 7 operational amplifiers, 2 diodes, 2 analog multipliers and the power supplies () was built. In the circuit are three integrators, which are created by the operational amplifiers together with the capacitors in feed-back.

Generally, real electrical elements as resistor, capacitor, coil and memristor are nonlinear and fractional order; however, their models are very often simplified [12]. This notion is based on fact that there is not ideal electrical element. Now, we already know that all real elements lie in between two. For instance, resistor created by wire has property of resistor and coil too. Moreover, we know element fractor (between resistor and capacitor) as well as fractductor (between resistor and coil).

Westerlund and Ekstam [14] in 1994 proposed a new linear capacitor model. It is based on Curie’s empirical law of 1889 which states that the current through a capacitor iswhere and are constant, is the DC voltage applied at , and .

For a general input voltage V(t), the current iswhere C is the capacitance of the capacitor. It is related to the kind of dielectric. Another constant (order) is related to the losses of the capacitor. Westerlund and Ekstam provided in their work the table of various capacitor dielectrics with appropriate constant which has been obtained experimentally by measurements.

For a current in the capacitor the voltage isTaking into account above consideration on capacitor, the new mathematical model of the real circuit implementation of the chaotic system (6) can be derived. Similarly to method provided in book [12], we can apply Kirchhoff’s circuit laws into the circuit depicted in Fig. 9 in paper [11] and derive a new fractional-order model of the chaotic system (6) in the following dimensionless form:where are the real orders of used capacitors.

For simulation purposes, we will use a numerical solution of equations (9) obtained by using the relationship (3), derived from the GLD (1), which leads to equations in the form:where is the simulation time, , for , and (x(0), y(0), z(0)) is the start point (initial conditions). The binomial coefficients , are calculated according to relation (4). Proposed numerical solution (10) was implemented as the MATLAB function, which syntax is listed in Appendix.

Simulation results and comments

First, let us use proposed numerical solution (10) to simulate chaotic system (6) for the parameter , simulation time s, time step , and initial conditions (). Using created MATLAB function, we can call it as follows:

In Fig. 1 is depicted chaotic behavior of the system (6) for and the initial conditions in state space.

Fig. 1

x–y–z state space of system (6) for

Due to the fact that parameter a should be within intervals or in order to observe chaotic behavior of the system, let us move from first interval to the second one and consider . Using created MATLAB function, we can call it as follows:

In Fig. 2 is depicted chaotic behavior of the system (6) for and the initial conditions in state space.

Fig. 2

x–y–z state space of system (6) for

Now, let us consider the fractional-order system (9), where we assume real order of the identical capacitors , for the parameter , simulation time s, time step , and initial conditions (). Using created MATLAB function, we can call it as follows:

In Figs. 4–6 are depicted chaotic behavior of the system (9) for , real orders , and the initial conditions in state planes (blue color), respectively. We can observe chaotic behavior when the total system order is 2.94.

Fig. 4

x–y state plane of system (9) for

Fig. 6

x–z state plane of system (9) for

In order to have more details, in Fig. 3 is presented the bifurcation diagram of the variable x versus the parameter q within interval [0.93; 1].

Fig. 3

Bifurcation diagram of system (9) when changing the bifurcation parameter q for the selected initial conditions and fixed parameter

The Lyapunov exponents of system (9) for and orders , computed by method described in [15], are: , which confirm chaos.

In addition, let us investigate existence of the coexisting attractors. We will consider the same orders and parameter as in previous example of the fractional-order system (9), but for different initial conditions (). Using created MATLAB function, we can call it as follows:

In Figs. 4, 5 and 6 are depicted chaotic behavior, coexisting attractors, of the system (9) for , real orders , and the initial conditions (blue color) and the initial conditions (red) in state planes (red color), respectively.

Fig. 5

y–z state plane of system (9) for

Conclusions

In this comment paper, we pointed out that the results presented in paper [11] can be enhanced. In spite of the fact that theory of fractional-order nonlinear (chaotic) systems is known and it was described in many papers and books [4, 7, 12], it was just for the system with known equilibrium, except few of them, see, e.g., [1-3].

In this note, a contribution to new phenomenon and almost unexplored area of a fractional-order chaotic no-equilibrium system was described together with the tool (MATLAB function) for its further investigation, even the coexisting attractors. Such nonlinear system without equilibrium is appropriate for practical applications. We have more degree of freedom due to the orders.