#### Abstract

This paper is presented on the theory and applications of the fractional-order chaotic system described by the Caputo fractional derivative. Considering the new fractional model, it is important to establish the presence or absence of chaotic behaviors. The Lyapunov exponents in the fractional context will be our fundamental tool to arrive at our conclusions. The variations of the model’s parameters will generate chaotic behavior, in general, which will be established using the Lyapunov exponents and bifurcation diagrams. For the system’s phase portrait, we will present and apply an interesting fractional numerical discretization. For confirmation of the results provided in this paper, the circuit schematic is drawn and simulated. As it will be observed, the results obtained after the simulation of the numerical scheme and with the Multisim are in good agreement.

#### 1. Introduction

In the last decade, chaos theory has attracted many researchers. Chaos theory is a field of mathematics that can be applied in many domains: from modeling chaotic financial systems [1], in representing circuit schematics [2, 3], and others. The chaotic systems and their circuit schematics have received many investigations in the literature, see, for example, [2–4]. With the new development of fractional calculus, some researchers were interested in modeling chaotic systems using the new and old fractional operators. The impact of the fractional derivatives is an interesting question which merits investigations. Many problems are opened in fractional calculus. The first question is how to construct the numerical scheme to obtain the phase portraits of chaotic systems. Many numerical techniques are used in fractional differential equations that can be applied naturally to solve chaotic systems. In [5], Danca finds interesting remarks that the characterization of the hyperchaotic system using two positive Lyapunov exponents is not an adequate definition in a fractional context. Therefore, modeling chaos using fractional operators is an open problem that admits some interesting interrogations which are not solved yet. Fractional calculus finds many recent advancements in fractional operators and their applicability in real-world problems. For recent advances of fractional calculus and its applications, see the following investigations [1, 6–10].

The chaos literature in integer version and fractional context is very long, but we recall the investigations we find essential for our study. In [11], Pacheco et al. study a new analysis of a new chaotic system with different families of hidden and self-excited attractors by using the fractional operator. In [12], Rajagopal et al. analyze using fractional-order derivative a memristor no equilibrium chaotic system; they proposed synchronizing their introduced model. In other words, they try to find adequate control to synchronize slave and master chaotic systems. In [2], Avalos-Ruiz et al. propose a control implementation on an FPGA for a class of fractional chaotic systems in the context of fractional variable order derivative with Mittag-Leffler kernel. In [13], Perez et al. propose the fractional-order chaotic system in the context of conformable Liouville–Caputo operator. At the same time, they present a novel numerical scheme for the used fractional operator to depict the phase portraits of their proposed chaotic system. In [14], Gomez et al. investigate the chaos in a calcium oscillation model using the Mittag-Leffler fractional derivative operator. In [15], Atangana and Gomez offer a new hyperchaotic system in the context of the Mittag-Leffler fractional operator. They provided a numerical scheme to depict the considered model’s phase portrait in the fractional operators’ context. They mainly apply their numerical procedure to the fractional Chua model. In [16], Li et al. introduce a new four-wing model with integer-order derivative, which will be subject to investigations in the fractional context. In [3], Akgul gives a chaotic system for modeling the memcapacitor model using a fractional-order derivative. In [17], Rajagopal et al. continue their exciting investigation of the fractional-order chaotic system in the fractional calculus field. They propose a work on the multiscroll chaotic system and give an application to the synchronization of their proposed model. In [18], Akgul et al. address a user interface for the generation of random numbers based on fractional and integer-order chaotic systems. For more investigations in chaos, see in the following investigations [4, 12–14, 19].

The novelties of the present paper are described in the following lines. This paper focuses on the phase space, the bifurcation diagrams, and the chaos’s characterization using the Lyapunov exponents for a fractional version of the Li et al. chaotic system. We analyze the impact of the fractional-order derivative on the dynamics of a fractional version of the system in [16] with Lyapunov exponents and bifurcation diagram aids. Its possible implementation in terms of electrical modeling has been proposed as well. One of the novelties of the present work is the Lyapunov exponents in the fractional context, which is an open problem. Note that the Lyapunov exponents’ properties in the integer context for hyperchaotic and chaotic behaviors are not adequate definitions in the fractional version. This remark was first discussed in [5]. Furthermore, the slight variation of the model’s parameters in the dynamics has been investigated using bifurcation diagrams. We provide in this paper that the initial conditions impact the chaotic behaviors in the fractional context. The circuit schematic of the fractional chaotic system and the results after simulation in Multisim confirm the theoretical findings obtained via the numerical scheme.

The paper is structured in the following form. In Section 2, we define the fractional operators’ necessaries for our work. In Section 3, we present the fractional-order chaotic system described by the Caputo derivative. Section 4 offers the numerical scheme used to obtain the considered fractional system’s phase portraits. In Section 5, we illustrate the numerical method by depicting the supposed system’s phase portraits with different fractional-order of the Caputo derivative. We also characterize the nature of the chaos using the Lyapunov exponents. In Section 6, we present the bifurcation diagrams and phase space according to the variation of the considered fractional-order chaotic system’s parameters. In Section 7, we give the possible implementation of the fractional-order chaotic system in modeling electrical circuits; this section will confirm the paper’s theoretical findings. Section 8 provides the stability analysis of the fractional-order chaotic system’s equilibrium points with the Matignon criterion used in fractional calculus. In Section 9, we analyze the sensitivity of the initial conditions in the fractional-order chaotic system. We finish with the final remarks and future perspectives of research.

#### 2. On Fractional Operators in Fractional Calculus

In this paper, we study a fractional differential system; before, it was indispensable to recall the fractional operators used in this work. Many fractional derivatives exist, those with singularities as the Caputo and Riemann–Liouville derivatives, those with nonsingularities as the exponential derivative, and the Mittag-Leffler fractional derivative. Due to space limitations, we provide the Caputo derivative and the Riemann–Liouville derivative and their generalizations.

*Definition 1 (see [20, 21]). *The integral defined by Riemann and Liouville called Riemann–Liouville integral of a given function iswith the function symbolizing the Gamma Euler function and under the order .

There exists a generalization of the Riemann–Liouville operator recently proposed in the literature; we have the following definition.

*Definition 2 (see [22]). *The generalized Riemann–Liouville fractional integral of a given function iswhere the orders and satisfy the following relationship and and the function gamma is , for all .

We can observe that when the order , we get the classical Riemann–Liouville integral operator. Furthermore, when , we have the integer version of the integral.

Caputo proposed another derivative due to the inconvenience of the Riemann–Liouville operator. We recall the Caputo derivative in the following definitions and its generalization.

*Definition 3 (see [20, 21]). *The Caputo fractional derivative of order can be represented as the form of the function :with the Gamma Euler function denoted by .

*Definition 4 (see [22]). *The generalized Caputo derivative of order is described as the form of the function :with the orders and satisfying the relations and and the Gamma function is denoted by , for all .

There exist many fractional operators, and all of them have their advantages. The use of the Caputo derivative in the present work has many motivations. Firstly, all fractional operators take into account the memory effect. Secondly, our main motivation is the initial conditions because the Riemann–Liouville derivative cannot consider them. In the Riemann–Liouville derivative case, the starting condition should be an integral form known in physics to be unrealistic. It cannot be satisfied in many real-world problems. Another basic motivation is in the context of Riemann–Liouville; the derivative of a constant is not null, contrary to the context of the Caputo derivative where the derivative of a constant is null.

#### 3. Fractional-Order Chaotic System

This section provides a fractional version of the chaotic model based on the system presented by Li et al. [16]. Before the investigation, we recall Li et al. chaotic equations represented with the integer-order derivative by the following equation:with initial conditions

The strange attractor is obtained at , , and according to Li et al. work. The Lyapunov exponents of the model presented in system (5) are , , and . In the original paper, the system presented in equation (5) is called a four-wing chaotic equation and admits three real equilibrium points, which are unstable and have two other complex equilibrium points. For more explanations related to the system in equation (5), see [16]. In our present work, we consider the memory effect generated by the fractional operators. Therefore, we investigate the following modified fractional-order chaotic system:with the initial condition defined by the following relation:

Our objective is to give a new fractional-order system; therefore, we set the following values for the parameters of model (7); there are , , , and . In terms of comparison, we can observe that equation (7) is a fractional version of equation (5). Furthermore, we can observe when the order of the Caputo derivative is in equation (7); we recover the constructive equation of original equation (5) presented in [16]. The difference is in the initial conditions. This paper will be to focus on detecting the chaotic behavior when the fractional-order derivative is utilized. The fractional-order operator impact will be analyzed and proved by the comparison of the attractors’ geometries. The variation of the model’s parameters such as and generate strange actuators and will be confirmed using the bifurcation diagrams. We will also establish via the stability analysis of the equilibrium points of our considered chaotic system; there exists a chaotic region according to the variation of the order of the Caputo derivative. Equation (7) is said the system with commensurate fractional-order when ; else, it is said the equation with incommensurate fractional-order.

#### 4. Numerical Scheme for the Fractional-Order Chaotic System

In this section, we apply the numerical scheme proposed by Garrapa [23] for our proposed fractional-order equation (7) necessary to depict the phase portraits. This numerical discretization uses the numerical approach of the Riemann–Liouville integral. In this section, we consider equation (7) with commensurate fractional-order. The first step will be to determine the analytical solution of fractional differential equation (7) using Riemann–Liouville integral:

We set the following function obtained from model (7):

For the graphical representations, we introduce the discretization at the point ; then, equations (9)–(11) can be written in the following form:

Introducing the step-size and , the discretization of the fractional integrals can be represented for the functions , and as the following:where the parameters are given by

For the simplification of the numerical scheme of the Riemann–Liouville integral, we introduce the first-order interpolant polynomial of the functions , , and given by the following relations [23]:

Replacing equations (20)–(22) into equations (16)–(18), we get, after recursive summation, more simple formulas for the numerical approach of the fractional integrals:where the parameters are given explicitly as the following form:and furthermore, when , then the parameter can be expressed as the following forms:

For the numerical schemes which are used to represent the graphics, it is obtained after plugging equations (23)–(27) into equations (13)–(15); we have the following expressions:

For the previous representations, it is necessary to recall the numerical discretization for our functions , , and ; we have the following:

The convergence and stability are not detailed in this section, see [23]; here, we just give information related to these two properties. We assume that , and are the approximate numerical solutions of the fractional system represented in equation (7), and , and are the exact solutions; the residual function with the Caputo derivative is given by the following expression:

We can notice that the convergence of the numerical scheme presented in this section is obtained when the step-size converges to 0 [23]. The stability of our numerical scheme is obtained from the Lipchitz criterion of the functions , , and .

#### 5. Illustration of the Numerical Scheme

This section presents the phase portraits of fractional-order system (7) with commensurate and incommensurate fractional orders. The main objective is to illustrate the numerical scheme presented in previous Section 4. We set the following values for the parameters of the model , , and and the initial conditions given by equation (9). For Figure 1, the considered order is .

**(a)**

**(b)**

**(c)**

**(d)**

For more illustrations, we consider system (7) with commensurate order . The phase portraits are presented in Figure 2.

**(a)**

**(b)**

**(c)**

**(d)**

An interesting procedure of comparison between the geometry of the attractors can be found in [24]. The method consists of comparing the amplitude of the attractors in different orders of the Caputo derivative. The comparison of the attractors’ geometry at the orders and can be analyzed for model (7). We adopt the sketch presented in [24]. We consider the phase portraits in plane in Figures 1 and 2 and construct polygons, calculate their areas, and compare them. We have following Figures 3 and 4.

In Figure 3, amplitudes are and , and in Figure 4, the amplitudes are and ; referring to the dimension and calculating the areas of the polygons and , we notice that the ; thus, there exists a significant difference between the attractors at the orders and . The method adopted in this section proves the influence of the order of the Caputo derivative. To see more the impact of the order of the Caputo derivative, we consider fractional chaotic equation (7) with incommensurate order:

System (31) is obtained when and in equation (7). We have the phase portraits in following Figure 5.

**(a)**

**(b)**

**(c)**

**(d)**

The Caputo derivative’s impact can be observed by comparing the phase portrait in the commensurate scenario and the phase portrait in the incommensurate scenario, which have different geometries. To detect equation (7) chaotic behavior with the orders and , we calculate the Lyapunov exponents. The Lyapunov exponents will permit us, in general, to recognize chaotic behaviors. The Jacobian matrix of fractional differential equation (7) with the same order is given by the following matrix, which is a fundamental tool for finding the Lyapunov exponents:

The Lyapunov exponents in the context of fractional calculus in the chaotic system using Matlab code are proposed in Danca’s work. Following the same proposed algorithm, the Lyapunov exponents for fractional-order equation (7) with commensurate order when are given as follows:

And, its associated Kaplan–Yorke dimension is given as follows:

Following the previous results, we can give the nature of the behaviors. First of all, the fractional system is dissipative because the sum of all Lyapunov exponents gives a negative number. The second property is to detect the chaotic behavior of our presented equation (7) with the same order. We can see with equation (33), we have one positive Lyapunov exponent corresponding to the chaotic behavior of equation (7). We repeat the same analysis with the order ; the Lyapunov exponents are represented as the following form:

And, its associated Kaplan–Yorke dimension is given as follows:

The system is dissipative as well because the addition of all Lyapunov exponents is negative. For the nature of the behaviors, we analyze equation (35); we have the same conclusion as in the previous case, equation (7) has chaotic behavior with the order . It is caused by the existence of one positive Lyapunov exponent. This section’s main conclusion is when the fractional derivative orders vary in the interval , fractional equation (7) with the same order has chaotic behavior. To confirm this assumption, we calculate in following Table 1 the Lyapunov exponents versus the fractional-order derivative’s variation into .

Checking the results on Table 1, we conclude that fractional-order equation (7) is dissipative at all orders in the interval . It is because, at all orders, the addition of the Lyapunov exponents is negative. We observe one positive Lyapunov exponent at all orders , which means the fractional-order system’s solutions have chaotic behavior.

#### 6. Stability of the Equilibrium Points in Fractional Context Calculus

In this part, we discuss the fractional-order system’s stability by using the Matignon criterion. It is a famous criterion used in fractional calculus. Stability analysis is one of the fundamental points in the chaotic and hyperchaotic system. In general, all equilibrium points when the system is chaotic fail to be stable. In this section, we consider equation (7) with commensurate order. It is also interesting to focus on the stability of its equilibrium points. The Jacobian matrix given in the previous section will be used again. We have the matrix described by the following expression:

The following points are the real equilibrium points when you exclude the complex eigenvalues:

The principle of the method consists of evaluating the Jacobian matrix at all points previously mentioned. For the first point , the Jacobian matrix is as follows:

The previous matrix has eigenvalues as the following numbers , , and . We notice that the third eigenvalue satisfies . The first eigenvalue obeys that for all . Thus, the equilibrium point is not stable for all orders of the Caputo derivative.

The matrix in equation (37) at the point is described as the following form:

The previous matrix has eigenvalues as the following numbers , , and . We remark that the third eigenvalue satisfies . The first eigenvalue obeys to for all and the last eigenvalue for all . Thus, the equilibrium point is not stable when exceeds 0.9.

The Jacobian at the point is described by

The previous matrix has eigenvalues as the following numbers , , and . We have that the third eigenvalue satisfies . The first eigenvalue obeys that for all and the last for all . Then, the equilibrium point is not stable when exceeds 0.9.

#### 7. Bifurcation Diagrams and Lyapunov Exponents

This section analyzes the impact caused by the variation of the parameters and on chaotic behavior. Small variations are applied to and , and we depict the responses via the bifurcation diagrams plus phase portraits. In this section, for more clarity in the curve, we consider the step-size . All findings are obtained with equation (7) in commensurate order.

We begin our analysis with the variation of the parameter in the interval . We fix and . In Figure 6, we represent the bifurcation diagram versus the parameter .

We notice that when the parameter varies in the region , fractional system (5)–(8) has chaotic behavior. Note that the variation of the parameter into generates the solutions’ convergence to a stable equilibrium point. For an illustration of the chaotic behavior, we represent the phase portraits of our system at and fixing and at the order , see Figure 7.

**(a)**

**(b)**

**(c)**

**(d)**

We confirm the chaotic behavior with the existence of one positive and large Lyapunov exponent:

And, its associated Kaplan–Yorke dimension is given as follows:

To validate the convergence of the solutions to a stable equilibrium point, we study the local stability of the nontrivial equilibrium point of equation (7) with commensurate order considering . The equilibrium points of equation (7) with commensurate order with are given by , , and . The Jacobian matrix at the point gives

The eigenvalues of the previous matrix are given by and . We can remark that the first eigenvalue does not obey to Matignon criterion; thus, the point is not local stable. Thus, the convergence of the solutions of equation (7) is not at the point . We now try to a point ; thus, the Jacobian matrix is given by

We obtain as eigenvalues the following values , , and . We can observe that all eigenvalues have negative real parts. Thus, the Matignon criterion is satisfied by all eigenvalues. Thus, the point is a stable equilibrium point. We finish with the point . The Jacobian matrix at the point is described as the following form:

The eigenvalues of the previous Jacobian matrix at this point are given by , , and . We can observe that all eigenvalues have negative real parts; thus, the Matignon criterion is automatically satisfied. Thus, the point coincides with the stable equilibrium point. Referring to the phase portraits, the trajectories converge to the point when .

We finish this section by the bifurcation diagram generated by the variation of the parameter , and we illustrate the impact of the parameter on the phase portraits, see Figure 8, taken at and .

**(a)**

**(b)**

**(c)**

**(d)**

The bifurcation diagram shown in Figure 9 indicates that , in our considered interval , the system has continuously chaotic behavior. Let us ; the chaotic behavior at this point is confirmed by the existence of one positive Lyapunov exponent; we have the following calculations:

And, its associated Kaplan–Yorke dimension is given as follows:For more confirmation of the chaotic region, we depict our fractional system’s phase portraits at and , see Figure 10.

**(a)**

**(b)**

**(c)**

**(d)**

#### 8. Electrical Circuit Schematic of the Fractional System

In this section, we give the circuit representation associated to chaotic system (7) with incommensurate order. For more details in the context of the fractional chaotic circuit, we consider equation (7) with written using electronic tools:where for denotes the value of the capacitors and represents the resistance value. To connect the results in this paper with the experimental results, we draw the circuit schematic and generate the phase portraits using oscilloscopes. We first use our numerical scheme in the previous sections to draw the phase portraits. We consider the order of the Caputo derivative as . The portraits of model (49)–(51) are represented in following Figure 11.

**(a)**

**(b)**

**(c)**

**(d)**

For the circuit schematic, the values of the capacitors are given by , and the values of the resistors are given , , , , , , and . We need to utilize the fractional integrators associated with the order [25]; we have the following formula:

In equation (52), the third-order approximation is used, but more suitable approximations can also be used. This section’s objective is the circuit schematic to validate our theoretical results; therefore, third- or fourth-order approximations can be used for illustration. The method of approximation in equation (52) is not unique. There also exists the 4^{th}-order approximation, which can generate a satisfactory magnitude and phase response over a considerable frequency band [26]. For more related works addressing the approximation of the fractional-order capacitor via fourth-order continued fraction expansion, see Mishra et al.’s investigation [26]. Utilizing the transfer function, the above function can be written as the following formula:

Identifying two above equations (45) and (46), we have the following values for the resistors and the capacitors for the fractional integrator, , , , , , and . The circuit schematic is done with the aids of 14 resistors, five capacitors, four multipliers for the 4 nonlinear functions, and others, see Figure 12.

The results corresponding to the simulation of the electronic circuit in the Multisim platform are shown in following Figures 13–15.

We can observe that the results with Matlab after simulating the numerical scheme in Figure 11 and the findings corresponding to the simulation of the electronic circuit in the Multisim platform shown in Figures 13–15 are in good agreement. This section connects the theoretical results with experimental results.

Before closing this section, we will implement the integer version of model equations (49)–(51). Therefore, the order should be . For the confirmation of our results obtained after the simulation of the electronic circuit in the Multisim platform, we first simulate the following differential equations:

Under the initial conditions given by system (9) and with and , the phase portraits are represented in following Figure 16.

**(a)**

**(b)**

**(c)**

**(d)**

For simulation of system (54)–(56), we rewrite equations (54)–(56) by using resistors and capacitors; we have the following representation:where the values of capacitors are given and resistors are given by , , , , , , and . The schematic of the above system is represented in following Figure 17.

The phase portraits obtained after the simulation of the circuit in Figure 17 are in following Figures 18–20.

We can observe that the theoretical results 16 and the experimental results 18, 19, and 20 are in good agreement.

#### 9. Sensitivity to Initial Conditions

In this section, we study the sensitivity of our model when the initial conditions change permanently. By the definition of chaos, the chaotic systems are very sensitive to the starting conditions. In this section, we try to measure the impact of the initial conditions via the bifurcation diagrams. For the rest of the paper, we consider equation (7) with commensurate order . The comparisons are made with the conditions given by

We first consider the variation of the parameter , and we influence the last initial condition . The bifurcation map21 associated with the initial conditions (blue) and (red) versus the variation of the parameter is given in following Figure 21.

The present bifurcation diagram, see Figure 21, informs us that the chaotic regions do not change significantly when the parameter varies. Note that, for the initial condition , the chaotic region is , and for the initial condition , the chaotic region is . The initial conditions’ impact can be observed at the beginning of the bifurcation diagram, in the interval , where system (7) trajectories converge to a stable equilibrium point. We notice that, with the first initial condition , the system enters first in the chaotic region. The small perturbation on the starting condition can generate changes but not significantly when the parameter varies. To see the initial condition’s influence according to the variation of , we illustrate by representing the phase portraits at and . When , we can observe that, in Figure 22, the phase space in plane (blue) corresponding to is chaotic, contrary to the phase space22 in plane (red) corresponding to , where clearly the solutions of the system converge to an stable equilibrium point. Let us now ; in Figure 22, the phase space in the plane (blue) corresponding to the condition is chaotic, and the phase space22 in the plane (red) corresponding to is also chaotic.

**(a)**

**(b)**

In the second, we suppose , , and are fixed. The variation of the parameter is also considered in the bifurcation diagram. The sensitivity to the starting condition can also be observed in the bifurcation diagram. The bifurcation diagram, see Figure 23, associated with the conditions (blue) and (red) versus the variation of the parameter is given in following Figure 23.

The bifurcation diagram, see Figure 23, informs us the chaotic regions for the starting condition of equation (7) are given by and for the new condition . The influence of the initial conditions can be noticed in the interval , after this region equation (7) has strange attractors with the conditions and . For illustration, see Figure 24, we consider , and we represent the phase space in the plane for more clarity.

**(a)**

**(b)**

We finish this section with the perturbation of the first initial condition. The bifurcation diagram, see Figure 25, associated to the initial conditions (blue) and (red) versus the variation of the parameter is given in the following figure.

The bifurcation diagram presented in Figure 25 can be interpreted as with the condition , the chaotic region does not change, and with , the new chaotic region is also the interval . In bifurcation diagrams, see Figures 21, 23, and 25, the general observation is that the chaotic regions do not significantly change when the parameter varies. We also observe that the chaotic region is conserved for all new conditions , , and . Thus, the small variation in the starting conditions does not affect the chaotic region of our present model.

#### 10. Conclusion

In this paper, we have proposed a new system of differential equations that exhibit chaotic behavior and then studied with the Caputo fractional derivative. We have discussed the nature of the chaos for different orders of the fractional derivative using Lyapunov exponents. It is noticed that the fractional-order system considered in this paper has chaotic behavior when the order of the fractional derivative is into . The small variation of the model’s parameters has been studied using bifurcation diagrams, which play an essential role in the chaos characterization. We also find that the initial condition’s choice is crucial because our fractional-order system is very sensitive to the initial conditions. The investigations related to the Lyapunov exponents in the context of the fractional derivative are an open problem. The incommensurate fractional-order chaotic system can be investigated for future work. For example, how the phase portraits can be obtained in the context of incommensurate order. How can the Lyapunov exponents be calculated for a fractional-order chaotic system with incommensurate order? What about the variations of the parameters? What will be the impact of the different fractional-order derivatives?

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The author declares that there are no conflicts of interest.