# Study of a Fractional-Order Chaotic System Represented by the Caputo Operator

Authors: Ndolane Sene
Journal: Complexity
 (2021)

---

## 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.


---

## Body

## 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. InSection 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 functionz:0,+∞⟶ℝ is(1)Iαzt=1Γα∫0tt−sα−1zsds,with the function Γ… symbolizing the Gamma Euler function and under the order α>0.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 functionz:0,+∞⟶ℝ is(2)Iα,ρzt=1Γα∫0ttρ−sρρα−1zsdss1−ρ,where the orders α and ρ satisfy the following relationship α and ρ>0 and the function gamma is Γ…, for all t>0.We can observe that when the orderρ=1, we get the classical Riemann–Liouville integral operator. Furthermore, when α=ρ=1, 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α∈0,1 can be represented as the form of the function z:0,+∞⟶ℝ:(3)Dcαzt=1Γ1−α∫0tt−s−αz′sds,with the Gamma Euler function denoted by Γ….Definition 4 (see [22]).
The generalized Caputo derivative of orderα∈0,1 is described as the form of the function z:0,+∞⟶ℝ:(4)Dα,ρzt=1Γ1−α∫0ttρ−sρρ−αz′sds,with the orders α and ρ satisfying the relations α∈0.1 and ρ>0 and the Gamma function is denoted by Γ…, for all t>0.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:(5)x′=ax+y+yz,y′=−xz+yz,z′=−bz−cxy+m,with initial conditions(6)x0=1,y0=−1,z0=1.The strange attractor is obtained ata=1, b=1, and m=1 according to Li et al. work. The Lyapunov exponents of the model presented in system (5) are LYE1=0.409, LYE2=0, and LYE3=−1.773. 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:(7)Dcκx=ax+y+yz,Dcβy=−xz+yz,Dcαz=−bz−cxy+m,with the initial condition defined by the following relation:(8)x0=x0=0.2,y0=y0=0.2,z0=z0=0.2.Our objective is to give a new fractional-order system; therefore, we set the following values for the parameters of model (7); there are a=1, b=1, c=1, and m=0. 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 κ=β=α=1 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 a and c 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:(9)xt=x0+Iαut,x1,(10)yt=y0+Iαvt,x1,(11)zt=z0+Iαwt,x1.We set the following function obtained from model (7):(12)ut,x1=ax+y+yz,vt,x1=−xz+yz,wt,x1=−bz−cxy.For the graphical representations, we introduce the discretization at the pointtn; then, equations (9)–(11) can be written in the following form:(13)xtn=x0+Iαutn,x1,(14)ytn=y0+Iαvtn,x1,(15)ztn=z0+Iαwtn,x1.Introducing the step-sizeh and tn=nh, the discretization of the fractional integrals can be represented for the functions u,v, and w as the following:(16)Iαutn,x1=hα∑j=1ndn−jutj,x1j,(17)Iαvtn,x1=hα∑j=1ndn−jvtj,x1j,(18)Iαwtn,x1=hα∑j=1ndn−jwtj,x1j,where the parameters are given by(19)dn−j=n−j+1α−n−jα1/Γ1+α,dn=n+1α−nα1/Γ1+α.For the simplification of the numerical scheme of the Riemann–Liouville integral, we introduce the first-order interpolant polynomial of the functionsuτ, vτ, and wτ given by the following relations [23]:(20)uτ=utj+1,x1j+1+τ−tj+1hutj+1,x1j+1−utj,x1j,(21)vτ=vtj+1,x1j+1+τ−tj+1hvtj+1,x1j+1−vtj,x1j,(22)wτ=wtj+1,x1j+1+τ−tj+1hwtj+1,x1j+1−wtj,x1j.Replacing equations (20)–(22) into equations (16)–(18), we get, after recursive summation, more simple formulas for the numerical approach of the fractional integrals:(23)Iαutn,x1=hαd¯nαu0+∑j=1ndn−jαutj,x1j,(24)Iαvtn,x1=hαd¯nαv0+∑j=1ndn−jαvtj,x1j,(25)Iαwtn,x1=hαd¯nαw0+∑j=1ndn−jαwtj,x1j,where the parameters are given explicitly as the following form:(26)d¯nα=n−1α−nαn−α−1Γ2+α,and furthermore, when n=1,2,…, then the parameter d can be expressed as the following forms:(27)d0α=1Γ2+α,dnα=n−1α+1−2nα+1+n+1α+1Γ2+α.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:(28)xtn=x0+hαd¯nαu0+∑j=1nκn−jαutj,x1j,ytn=y0+hαd¯nαv0+∑j=1nκn−jαvtj,x1j,ztn=z0+hαd¯nαw0+∑j=1nκn−jαwtj,x1j.For the previous representations, it is necessary to recall the numerical discretization for our functionsu, v, and w; we have the following:(29)utj,x1j=axj+yj+yjzj,vtj,x1j=−xjzj+yjzj,wtj,x1j=−bzj−cxjyj.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 xtn,ytn, and ztn are the approximate numerical solutions of the fractional system represented in equation (7), and xn,yn, and zn are the exact solutions; the residual function with the Caputo derivative is given by the following expression:(30)xtn−xn=Ohminα+1,2,ytn−yn=Ohminα+1,2,ztn−zn=Ohminα+1,2.We can notice that the convergence of the numerical scheme presented in this section is obtained when the step-sizeh converges to 0 [23]. The stability of our numerical scheme is obtained from the Lipchitz criterion of the functions u, v, and w.


## 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 a=1, b=1, and c=1 and the initial conditions given by equation (9). For Figure 1, the considered order is α=0.95.Figure 1
Phase portraits with the orderα=0.95.
(a)(b)(c)(d)For more illustrations, we consider system (7) with commensurate order α=0.90. The phase portraits are presented in Figure 2.Figure 2
Phase portraits with the orderα=0.9.
(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 α=0.9 and α=0.95 can be analyzed for model (7). We adopt the sketch presented in [24]. We consider the phase portraits in plane x,y in Figures 1 and 2 and construct polygons, calculate their areas, and compare them. We have following Figures 3 and 4.Figure 3
Phase spaces with polygon with the orderα=0.9.Figure 4
Phase spaces with polygon with the orderα=0.95.In Figure3, amplitudes are EF=12.2cm and AB=CD=4cm, and in Figure 4, the amplitudes are KL=10.88cm and GH=IJ=3cm; referring to the dimension and calculating the areas of the polygons AFBDEC and GLHJKI, we notice that the areaAFBDEC>areaGLHJKI; thus, there exists a significant difference between the attractors at the orders α=0.9 and α=0.95. 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:(31)x′=ax+y+yz,y′=−xz+yz,Dcαz=−bz−cxy.System (31) is obtained when κ=1 and β=1 in equation (7). We have the phase portraits in following Figure 5.Figure 5
Phase portraits with the orderα=0.9 in equations (32)–(34).
(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 α=0.9 and α=0.95, 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:(32)J=a1+zy−zz−x+y−cy−cx−b.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 α=0.95 are given as follows:(33)LYE1=0.3449,LYE3=0,LYE3=−2.5326.And, its associated Kaplan–Yorke dimension is given as follows:(34)dimLYE=2+LYE1+LYE2LYE3=2.1362.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 α=0.9; the Lyapunov exponents are represented as the following form:(35)LYE1=0.3775,LYE2=0,LYE3=−3.0042.And, its associated Kaplan–Yorke dimension is given as follows:(36)dimLYE=2+LYE1+LYE2LYE3=2.1257.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 α=0.9. 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 0.9,1, 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 0.9,1.Table 1
Lyapunov exponents.
αLE1LE2LE30.900.37750−3.00420.910.35810−2.78580.920.30830−2.83790.930.27010−2.73370.940.29790−2.55460.950.34490−2.53260.960.21770−2.30490.970.21070−2.16710.980.21710−2.09310.990.24020−2.07010.99500.22150−1.9927Checking the results on Table1, we conclude that fractional-order equation (7) is dissipative at all orders in the interval 0.9,1. 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:(37)J=a1+zy−zz−x+y−cy−cx−1.The following points are the real equilibrium points when you exclude the complex eigenvalues:(38)P1=0,0,0,P2=−1.41,−1.41,−2,P3=1.41,1.41,−2.The principle of the method consists of evaluating the Jacobian matrix at all points previously mentioned. For the first pointP1, the Jacobian matrix is as follows:(39)JP1=11000000−1.The previous matrix has eigenvalues as the following numbersλ1=1, λ2=0, and λ3=−1. We notice that the third eigenvalue satisfies argλ3=π>απ/2. The first eigenvalue obeys that argλ1≈0<απ/2 for all α. Thus, the equilibrium point is not stable for all orders α of the Caputo derivative.The matrix in equation (37) at the point P2 is described as the following form:(40)JP2=1−1−1.412−201.411.41−1.The previous matrix has eigenvalues as the following numbersλ1=0.1230+1.8776i, λ2=0.1230−1.8776i, and λ3=−2.246. We remark that the third eigenvalue satisfies argλ3=π>απ/2. The first eigenvalue obeys to argλ1≈43π/90>απ/2 for all α<0.9 and the last eigenvalue argλ2≈43π/90>απ/2 for all α. Thus, the equilibrium point is not stable when α exceeds 0.9.The Jacobian at the pointP3 is described by(41)JP2=1−11.412−20−1.41−1.41−1.The previous matrix has eigenvalues as the following numbersλ1=0.1230+1.8776i, λ2=0.1230−1.8776i, and λ3=−2.246. We have that the third eigenvalue satisfies argλ3=π>απ/2. The first eigenvalue obeys that argλ1≈43π/90>απ/2 for all α<0.9 and the last argλ2≈43π/90>απ/2 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 parametersa and c on chaotic behavior. Small variations are applied to a and c, 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 h=0.005. All findings are obtained with equation (7) in commensurate order.We begin our analysis with the variation of the parametera in the interval 0,5. We fix b=1 and c=1. In Figure 6, we represent the bifurcation diagram versus the parameter a.Figure 6
Bifurcation diagram versus the parametera at order α=0.95.We notice that when the parametera varies in the region 0.6,5, fractional system (5)–(8) has chaotic behavior. Note that the variation of the parameter a into 0,0.5 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 a=1.5 and fixing b=1 and c=1 at the order α=0.95, see Figure 7.Figure 7
Phase portraits with the orderα=0.95 and a=1.5.
(a)(b)(c)(d)We confirm the chaotic behavior with the existence of one positive and large Lyapunov exponent:(42)LYE1=0.1918,LYE2=0,LYE3=−2.4806.And, its associated Kaplan–Yorke dimension is given as follows:(43)dimLYE=2+LYE1+LYE2LYE3=2.0773.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 a=0.4. The equilibrium points of equation (7) with commensurate order with a=0.4 are given by O0,0,0, A−1.18,−1.18,−1.4, and B1.18,1.18,−1.4. The Jacobian matrix at the point O0,0,0 gives(44)JO=0.41000000−1.The eigenvalues of the previous matrix are given byλ1=0.4,λ2=0, and λ3=−1. We can remark that the first eigenvalue does not obey to Matignon criterion; thus, the point O0,0,0 is not local stable. Thus, the convergence of the solutions of equation (7) is not at the point O0,0,0. We now try to a point A−1.18,−1.18,−1.4; thus, the Jacobian matrix is given by(45)JA=0.4−0.4−1.181.4−1.401.181.18−1.We obtain as eigenvalues the following valuesλ1=−0.0763+1.4507i, λ2=−0.0763−1.4507i, and λ3=−1.8474. We can observe that all eigenvalues have negative real parts. Thus, the Matignon criterion is satisfied by all eigenvalues. Thus, the point A−1.18,−1.18,−1.4 is a stable equilibrium point. We finish with the point B1.18,1.18,−1.4. The Jacobian matrix at the point B is described as the following form:(46)JC=0.4−0.41.181.4−1.40−1.18−1.18−1.The eigenvalues of the previous Jacobian matrix at this point are given byλ1=−0.0763+1.4507i, λ2=−0.0763−1.4507i, and λ3=−1.8474. We can observe that all eigenvalues have negative real parts; thus, the Matignon criterion is automatically satisfied. Thus, the point C coincides with the stable equilibrium point. Referring to the phase portraits, the trajectories converge to the point B when a=0.4.We finish this section by the bifurcation diagram generated by the variation of the parameterc∈0,25, and we illustrate the impact of the parameter c on the phase portraits, see Figure 8, taken at a=1 and c=4.Figure 8
Phase portraits with the orderα=0.95.
(a)(b)(c)(d)The bifurcation diagram shown in Figure9 indicates that , in our considered interval 0,25, the system has continuously chaotic behavior. Let us c=4; the chaotic behavior at this point is confirmed by the existence of one positive Lyapunov exponent; we have the following calculations:(47)LYE1=0.2485,LYE2=0,LYE3=−2.4311.Figure 9
Bifurcation diagram versus the parameterc at order κ=0.95.And, its associated Kaplan–Yorke dimension is given as follows:(48)dimLYE=2+LYE1+LYE2LYE3=2.1022.For more confirmation of the chaotic region, we depict our fractional system’s phase portraits at a=1 and c=0.5, see Figure 10.Figure 10
Phase portraits with the orderα=0.95.
(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 κ=β=1 written using electronic tools:(49)x′=1R1C1x+1R2C1y+1R3C1yz,(50)y′=−1R4C1xz+1R5C2yz,(51)Dcαz=−1R6C3z−1R7C3xy,where Ci for i=1,…,n denotes the value of the capacitors and Ri 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 α=0.95. The portraits of model (49)–(51) are represented in following Figure 11.Figure 11
Phase portraits with incommensurate fractional system (42)–(44).
(a)(b)(c)(d)For the circuit schematic, the values of the capacitors are given byC1=C2=C3=1nF, and the values of the resistors are given R1=250kΩ, R2=250kΩ, R3=6.25kΩ, R4=6.25kΩ, R5=6.25kΩ, R6=250kΩ, and R7=6.25kΩ. We need to utilize the fractional integrators associated with the order α=0.95 [25]; we have the following formula:(52)1s0.95≈1.2862s2+18.6004s+2.0833s3+18.4738s2+2.6547s+0.003.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 4th-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:(53)Ts=1/C3s+1/R9C3+1/C4s+1/R10C4+1/C5s+1/R15C5.Identifying two above equations (45) and (46), we have the following values for the resistors and the capacitors for the fractional integrator, C3=1.27nF, C4=4.7nF, C5=3.63nF, R9=700MΩ, R10=1.5MΩ, and R15=15kΩ. 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.Figure 12
Circuit schematic for fractional-order system (42)–(44).The results corresponding to the simulation of the electronic circuit in the Multisim platform are shown in following Figures13–15.Figure 13
x,y space portrait of the Multisim simulation of equations (49)–(51).Figure 14
x,z space portrait of the Multisim simulation of equations (49)–(51).Figure 15
y,z space portrait of the Multisim simulation of equations (49)–(51).We can observe that the results with Matlab after simulating the numerical scheme in Figure11 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 α=1. 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:(54)x′=ax+y+yz,(55)y′=−xz+yz,(56)z′=−z−cxy.Under the initial conditions given by system (9) and with a=1 and c=1, the phase portraits are represented in following Figure 16.Figure 16
Phase portraits with system (54)–(56).
(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:(57)x′=1R1C1x+1R2C1y+1R3C1yz,(58)y′=−1R4C1xz+1R5C2yz,(59)z′=−1R6C3z−1R7C3xy,where the values of capacitors are given C1=C2=C3=1nF and resistors are given by R1=250kΩ, R2=250kΩ, R3=6.25kΩ, R4=6.25kΩ, R5=6.25kΩ, R6=250kΩ, and R7=6.25kΩ. The schematic of the above system is represented in following Figure 17.Figure 17
Circuit schematic for system (57)–(59).The phase portraits obtained after the simulation of the circuit in Figure17 are in following Figures 18–20.Figure 18
x,y space portrait of the Multisim simulation of equations (54)–(56).Figure 19
x,z space portrait of the Multisim simulation of equations (54)–(56).Figure 20
y,z space portrait of the Multisim simulation of equations (54)–(56).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 α=0.95. The comparisons are made with the conditions given by(60)x0=x0=0.2,y0=y0=0.2,z0=z0=0.2.We first consider the variation of the parametera, and we influence the last initial condition z0. The bifurcation map21 associated with the initial conditions 0.2,0.2,0.2 (blue) and 0.2,0.2,0.5 (red) versus the variation of the parameter a is given in following Figure 21.Figure 21
Bifurcation diagram versus the parametera at 0.2,0.2,0.2 and 0.2,0.2,0.5.The present bifurcation diagram, see Figure21, informs us that the chaotic regions do not change significantly when the parameter a varies. Note that, for the initial condition 0.2,0.2,0.2, the chaotic region is 0.6,5, and for the initial condition 0.2,0.2,0.5, the chaotic region is 0.7,5. The initial conditions’ impact can be observed at the beginning of the bifurcation diagram, in the interval 0,1, where system (7) trajectories converge to a stable equilibrium point. We notice that, with the first initial condition 0.2,0.2,0.2, 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 a, we illustrate by representing the phase portraits at a=0.6 and a=0.7. When a=0.6, we can observe that, in Figure 22, the phase space in x,y plane (blue) corresponding to 0.2,0.2,0.2 is chaotic, contrary to the phase space22 in plane x,y (red) corresponding to 0.2,0.2,0.5, where clearly the solutions of the system converge to an stable equilibrium point. Let us now a=0.7; in Figure 22, the phase space in the plane x,y (blue) corresponding to the condition 0.2,0.2,0.2 is chaotic, and the phase space22 in the plane x,y (red) corresponding to 0.2,0.2,0.5 is also chaotic.Figure 22
(a) Portrait ata=0.6 in the plane x,y. (b) Portrait at a=0.6 in the plane x,y.
(a)(b)In the second, we supposey0=0.5, x0=0.2, and z0=0.2 are fixed. The variation of the parameter a 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 0.2,0.2,0.2 (blue) and 0.2,0.2,0.5 (red) versus the variation of the parameter a is given in following Figure 23.Figure 23
Bifurcation diagram versus the parametera at 0.2,0.2,0.2 and 0.2,0.5,0.2.The bifurcation diagram, see Figure23, informs us the chaotic regions 0.6,5 for the starting condition of equation (7) are given by 0.2,0.2,0.2 and 0.7,5 for the new condition 0.2,0.5,0.2. The influence of the initial conditions can be noticed in the interval 0.6,0.7, after this region equation (7) has strange attractors with the conditions 0.2,0.2,0.2 and 0.2,0.5,0.2. For illustration, see Figure 24, we consider a=0.65, and we represent the phase space in the plane x,z for more clarity.Figure 24
(a) Portrait ata=0.65 in the plane x,z. (b) Portrait at a=0.65 in the plane x,y.
(a)(b)We finish this section with the perturbation of the first initial condition. The bifurcation diagram, see Figure25, associated to the initial conditions 0.2,0.2,0.2 (blue) and 0.5,0.2,0.2 (red) versus the variation of the parameter a is given in the following figure.Figure 25
Bifurcation diagram versus the parametera at 0.2,0.2,0.2 and 0.5,0.2,0.2.The bifurcation diagram presented in Figure25 can be interpreted as with the condition 0.2,0.2,0.2, the chaotic region does not change, and with 0.5,0.2,0.2, the new chaotic region is also the interval 0.7,5. In bifurcation diagrams, see Figures 21, 23, and 25, the general observation is that the chaotic regions do not significantly change when the parameter a varies. We also observe that the chaotic region 0.7,5 is conserved for all new conditions 0.2,0.2,0.5, 0.2,0.5,0.2, and 0.5,0.2,0.2. 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 into0.9,1. 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?



---

*Source: 5534872-2021-06-14.xml*