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][10][11]. Moreover, there is a significant interest in studying hidden attractors because they play an important role in both theo-I. Petráš (B) Faculty of BERG, Technical University of Košice, B. Němcovej 3, Košice 042 00, Slovak Republic e-mail: ivo.petras@tuke.sk retical 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 fractionalorder 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.

Definition of fractional-order operator
Fractional calculus is a generalization of integration and differentiation to non-integer (fractional) order fundamental operator a D α b , where a and b are the bounds of the operation. The standard notation for denoting the left-sided fractional-order integrodifferential operator of a function f (t) defined in the interval [a, b] is a D α t f (t), with α ∈ R. 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 z is the floor function, i.e., the greatest integer smaller than z, and α are the binomial coefficients for α 0 = 1. 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 kh, (k = 1, 2, . . . ) has the following form [13]: where L m is the "memory length", t k = kh, h 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 c (q) j ( j = 0, 1, . . . ) are binomial coefficients. For their calculation we can use the following expression: Then, general numerical solution of the fractional differential equation can be expressed as For the memory term expressed by the sum, a "short memory" principle can be used for various L m .
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 (a > 0). 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 (10 kΩ), 3 capacitors (10 nF), 7 operational amplifiers, 2 diodes, 2 analog multipliers and the power supplies (±15 V DC ) 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 is where h 1 and α are constant, V 0 is the DC voltage applied at t = 0, and 0 < α < 1, (α ∈ R).
For a general input voltage V (t), the current is where 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 is Taking 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 q 1 , q 2 , q 3 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 T sim is the simulation time, k = 1, 2, 3 . . . , N , for N = [T sim / h], and (x(0), y(0), z(0)) is the start

Simulation results and comments
First, let us use proposed numerical solution (10) to simulate chaotic system (6)  Due to the fact that parameter a should be within intervals 1.305 < a < 1.381 or 1.224 < a < 1.303 in order to observe chaotic behavior of the system, let us move from first interval to the second one and consider a = 1.3. Using created MATLAB function, we can call it as follows:  In Figs. 4-6 are depicted chaotic behavior of the system (9) for a = 1.3, real orders q 1 = q 2 = q 3 = 0.98, and the initial conditions (x(0), y(0), z(0)) = (0, 0.1, 0) in state planes (blue color), respectively. We can observe chaotic behavior when the total system order is 2.94.
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].

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][2][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.