Research article

A new fractional order 6D chaotic model: Study of model dynamics, system structure graph, electronic circuit realization and fractional control

  • M. Higazy a, b, * ,
  • Norah Almalki , a ,
  • Shabbir Muhammad , c, d ,
  • A. Al-Ghamdi , e
Expand
  • a Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
  • b Department of Physics and Engineering Mathematics, Faculty of Electronic Engineering, Menoufia University, Menouf 32952, Egypt
  • c Research Center for Advanced Materials Science (RCAMS), King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
  • d Department of Chemistry, College of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
  • e Department of Physics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*E-mail addresses: (M. Higazy),
(N. Almalki),
(S. Muhammad),
(A. Al-Ghamdi).

Received date: 2022-03-15

  Revised date: 2022-04-04

  Accepted date: 2022-04-07

  Online published: 2022-04-13

Abstract

A new fractional 6D chaotic model is constructed in this paper. The new fractional 6D chaotic model has six positive parameters plus the fractional order with eight nonlinear terms. The complicated chaotic dynamics of the new fractional 6D model is presented and analyzed. The basic properties of this model are studied and its chaotic attractors, dissipative feature, symmetry, equilibrium points, Lyapunov Exponents are investigated. The new dynamics of the 6D fractional model is numerically simulated using Matlab software. In addition, utilizing the graph theory tools certain structural characteristics are calculated. An electrical circuit is built to implement the new 5.4 fractional order 6D model. Finally, an active fractional order controller is proposed to control the new model at different fractional orders. The chaos of the new model is very useful and can be used to produce random keys for data encryption.

Highlights

A new fractional 6D chaotic model is constructed.

New chaotic attractors and its dynamical properties are investigated in context of fractional calculus.

New structural characteristics are calculated for a new 6D fractional chaotic model.

An electronic circuit is implemented to realize the new fractional chaotic model.

A new fractional controller is designed for the new model.

The chaos of the new model is very useful and can be used to produce random keys for data encryption.

Cite this article

M. Higazy , Norah Almalki , Shabbir Muhammad , A. Al-Ghamdi . A new fractional order 6D chaotic model: Study of model dynamics, system structure graph, electronic circuit realization and fractional control[J]. Journal of Ocean Engineering and Science, 2024 , 9(2) : 112 -125 . DOI: 10.1016/j.joes.2022.04.002

1. Introduction

In [1], the first 3D autonomous chaotic system was found by Lorenz when he investigated atmospheric convection. Lorenz system has been taken as a first model for studying chaos. After that, in [2], Rӧssler built a simpler chaotic 3D system. It is noteworthy that during the past twenty years, the chaos of technological paradigms, like nonlinear electric circuits, has progressively shifted from a pure science to a promising subject that has important applications. It is observed that intentionally generating chaos is a vital case in several technical implementations. In this direction, Chen built a three-dimensional autonomous chaos model [3,4] based on an engineering approach to feedback control, after that Lü has built his system [5] and general model combining all cases as a special case has built and named the unified model [6]. Following Vanecek and Celikovsky's classification of canonical form [7], in the 3D autonomous paradigms that have quadratic nonlinearities, considering the linear term (B=[bij] ), the Lorenz model fulfils b12b21>0 , the Chen model fulfils b12b21<0 , and the Lü model fulfils b12b21=0 although they aren't topologically equivalent. With this direction, they form together a full collection of general Lorenz dynamical models. From another point of view, the Rӧssler system contains a quadratic cross product term and is not part of the general family of Lorenz systems mentioned above. Recently, during the study of a 3D Rӧssler-type autonomous chaotic model with quadratic cross product terms, Liu and Chen [8] constructed a chaotic model with a quadratic product term in every kernel, that could give two attractors at the same time each of which has a double-scroll.
The concept of fractional calculus has been recognized since the evolution of the ordinary one, and its origin probably relates to a 1695 communication of Leibniz and L' Hospital. They discussed the concept of 1.5 order derivative [9]. Despite the long age of the fractional derivation, its applications in physics and engineering are still recent. Much research has been devoted to the creation of chaos using autonomous fractional order nonlinear models. In [10], fractional order chaos and hyper chaos were examined for the Rӧssler model. In [11], system synchronization of fractional order was performed. Control and stability for a collection of fractional-order nonlinear models have been investigated with Caputo derivative in [12]. The fractional order derivative provides many features to the dynamical system. It takes the history of the system dynamics into account. And the order of the derivative increases the number of parameters of the chaotic model. So, we can get new chaotic attractors and increases the dimension of parameters' space which increase its efficiency in secure communication applications. Also, in modelling, the fractional order can be used as a tuner parameter for the response to meet the modelled phenomena. But the classical derivatives and the models based on them do not consider any memories and hereditary effects. Extra fractional and integer order chaotic models and applications can be found in [13], [14], [15], [16], [17], [18], [19], [20] and [30], [31], [32], [33], [39], [40], [41].
Dynamical systems, especially the chaotic ones are widely utilized to model Ocean and atmospheric dynamics. Ocean engineering is interested with large scale wave motions in the Ocean (for example, wind waves, tsunami waves) and air along with the temperature and density of water in the Ocean. Ocean and atmosphere are equally significant in transferring energy between positions in the form of wave propagation. Dynamical systems, especially the chaotic ones help to model various phenomena in Oceans. Several research works have highlighted the values of eddies and chaos inside the Ocean in controlling the response to forcing from climate change. Certain physical processes in the Ocean produce highly non-linear or chaotic fluctuations. In [34], Jianmin Yang and Wenyue Lu numerically studied the generation and evolution of the super-rogue waves. In [35], the reader can see the importance of dynamical systems in modelling and evaluation of dynamics of jack-up platform system under wave, wind, earthquake, and tsunami loads. More numerical and computational studies of dynamical systems can be found in [36], [37], [38].
Encouraged by the previous studies, a new fractional chaotic 6D model is constructed which can be used to model may Oceanic phenomena. We proved the chaos nature of the proposed new fractional 6D model numerically by displaying its dynamics related to the order of the fractional differentiation. The Lyapunov exponents and bifurcation maps of the proposed model is presented. The impact of the order of the fractional differentiation on the Lyapunov exponent are presented. The structural properties of the proposed 6D model are studied via certain graph theory tools. In addition, the proposed 6D 5.4-fractional order chaotic model is realized by constructing an electronic circuit and we present its response. Also, an active fractional controller is produced to control the presented 6D fractional chaotic model.
The remaining of the manuscript is arranged as: the preliminaries and basic definitions of fractional differentiation are given in Section 2. The 6D chaotic model is constructed and its basic characteristics are given in Section 3. The new attractors are displayed in Section 4. Certain structural features of the proposed 6D chaotic model are presented using graph theory tools in Section 5. Section 6 contains the electronic circuit realization of the proposed fractional 6D chaotic model. An active fractional controller is designed and applied for the studied chaotic model in Section 7. The conclusions of our study are put in Section 8.

2. Fractional derivative and its preliminaries and basic definitions

This section briefly describes the fractional derivative and its corresponding fractional integral. Several notations were used in the new definition of fractional derivative. Through this article, we use lDtvY(t) and lItvY(t) respectively for the fractional differentiation and its corresponding integration, where l signifies the integration lower limit and v is the order of the fractional derivative or fractional integral. For each vC , and for any function Y:CC , lDtvY(t) and lItvY(t) should meet the following criteria [21]:
If Y(.) is analytic function, then the operator lDtvY(t) and lItvY(t) are analytic w.r.t. vandt .
The fractional differentiation operators satisfy the following linearity rules:
lDtv(c1Y1(t)+c2Y2(t))=c1lDtvY1(t)+c2lDtvY2(t)
lItv(c1Y1(t)+c2Y2(t))=c1lItvY1(t)+c2lItvY2(t)
The fractional operators with zero order must have no effect on Y(t) , i.e.,
lDt0Y(t)=lIt0Y(t)=Y(t).
If v is a positive integer, the fractional operator should give the same result as a normal derivative or integral.
In the coming parts, the Caputo fractional operators of differentiation and integration are introduced in its common sense which defined as following [22] and shall be used to construct the proposed new chaotic model.
lCDtvY(t)=1Γ(nv)ltY(n)(θ)(tθ)nv1dθ,
lCItvY(t)=1Γ(v)ltY(θ)(tθ)v1dθ.
For more details and other definitions as Atangana-Baleanu derivative and Caputo-Fabrizio derivative with its possible applications can be found in [23].
Definition 1. The following formula represents the Laplace transformation of Caputo fractional order operator L{D(ε)f(t)}=sεF(s)i=0k1sεi1f(i)(0);k1<ε<k,kN;F(s)=L{f(t)}.
Definition 2. For the graph G, define the matrix A=A(G)=[aij], such that aij=1 if ui is adjacent to uj in G and aij=0 otherwise. A is called adjacency matrix of G.
Definition 3. The Hermitian matrix of a digraph D is defined as:
Huv(D)={1if(u,v)and(v,u)areedgesinD;iif(u,v)isanedgeinDbut(v,u)isnot;iif(v,u)isanedgeinDbut(u,v)isnot;0otherwise.
The basics of the graph theory and its applications can be found for example in [24]

3. The novel 6D chaotic model and its basic characteristics

In the following, let us write 0CDtυ as Dυ . The 6D nonlinear fractional autonomous novel system is described by
Dυy1=α1(y2y1)+y2y3y4,Dυy2=β1(y2+y1)y1y3y4,Dυy3=γy3+y1y2y4+y5y6y4,Dυy4=δy4+y1y2y3+y5y6y3,Dυy5=α2(y6y5)+y6y3y4,Dυy6=β2(y6+y5)y5y3y4.
where yi (i=1,2,3,,6 ) are the state variables of the proposed model, and α1,α2,β1,β2,γ,δ are the model parameters and all of them are positive real constants.
Constructing this model and determining the values of its parameters such that the model has chaotic dynamic follow some common concepts of chaotification [25,26], to build an autonomous chaotic paradigm or to chaotic a non-chaotic autonomous paradigm, the coming framework rules need to be satisfied:
i) Dissipative of the model, in other words, the energy of the model is diminishing (except if Hamiltonian models are taken into consideration).
ii) The model has an unstable equilibrium. In other words, the equilibrium-evaluated Jacobian value has unstable eigenvalues.
iii) The model consists of one cross product term at least. That is, the dynamic effects between different variables can be considered.
iv) All system trajectories are bounded. In other words, dynamic equilibrium is maintained by increasing and decreasing system energy.
The method of chaotifying the discrete models have some systematic steps [25], but there are no general techniques for non-discrete cases. In the latter case, the above analytical conditions usually need to be combined with trial-and-error computer experiments to accomplish the needed chaotic job. System (6) is also has such a shape. Now, some basic features of (6) are examined.
In the following subsections, some basic characteristics of the proposed system (6) are examined.

3.1. Dissipative feature for the proposed model (6)

The 6D nonlinear fractional autonomous novel system is described by
Dυy1=α1(y2y1)+y2y3y4=v1,Dυy2=β1(y2+y1)y1y3y4=v2,Dυy3=γy3+y1y2y4+y5y6y4=v3,Dυy4=δy4+y1y2y3+y5y6y3=v4,Dυy5=α2(y6y5)+y6y3y4=v5,Dυy6=β2(y6+y5)y5y3y4=v6.
For this model (7), let V=[v1,v2,...,v6]T , we have
·V=i=16yi(Dυyi)=(β1+β2)(α1+α2+δ+γ)=D.
Then, system (7) is dissipative, if
(β1+β2)(α1+α2+δ+γ)<0
Suppose that Θ(t) is a suitable region in R6 which has a smooth boundary and let Θ(t)Λ(t) such that Λ(t) is the flow of the vector field V(t) . Plus, let that the volume of Θ(t) is V(t) . Referring to Liouville's theorem as in [13], we get
V˙(t)=Θ(t)(·V)k=16dvk
Then from (7) and (9), we have
V˙(t)=Θ(t)(D)k=16dvk=DV(t)
Integrating (11), then
V(t)=V(0)eDt
from (12), we conclude that the orbit of this dynamic model (6) forces all contained volumes to be reduced to zero. Therefore, there are attractors that attract all the orbits of the new model. From this, we can show that the newly proposed chaos model (6) is dissipative in the following condition is met:
(β1+β2)(α1+α2+δ+γ)<0

3.2. Model (6) equilibrium points

The equilibrium points of the proposed model (6) can be found by simultaneously solving the next equations:
α1(y2y1)+y2y3y4=0,β1(y2+y1)y1y3y4=0,γy3+y1y2y4+y5y6y4=0,δy4+y1y2y3+y5y6y3=0,α2(y6y5)+y6y3y4=0,β2(y6+y5)y5y3y4=0.
Clearly, P0=(0,0,0,0,0,0) is one of the equilibrium points. The system is very complex and has many other equilibrium points that can be calculated using Matlab. Here, the stability of the zero equilibrium P0 is examined only. The Jaccobian can be obtained by linearizing system (6) at P0 as:
J=(α1α10000β1β1000000γ000000δ000000α2α20000β2β2)
That has the following eigenvalues:
λ1,2=12β112α1±12α12+6α1β1+β12,
λ3=δ , λ4=γ ,
λ5,6=12β212α2±12α22+6α2β2+β22.
Since all model parameters are positive real numbers, it is easy to find that λ2,λ6>0 , proving the instability of the equilibrium P0 . The other nonzero equilibria, it is still possible to numerically evaluate their stabilities. Since two unstable equilibriums has been got at P0 , this long numerical calculation becomes not needed then it is not examined further here.

4. Displaying of new chaotic attractors

In the following, to simplifying the study, the produced model parameters are taken as: α1=α2=a ,β1=β2=b ,γ=c and δ=d . Then proposed model is rewritten as follows.
Dυy1=a(y2y1)+y2y3y4,Dυy2=b(y2+y1)y1y3y4,Dυy3=cy3+y1y2y4+y5y6y4,Dυy4=dy4+y1y2y3+y5y6y3,Dυy5=a(y6y5)+y6y3y4,Dυy6=b(y6+y5)y5y3y4.
Considering condition (12), with initial conditions Y(0)=(1.77,1.66,7.2,1.5,1.77,1.66) and fractional order v=0.8515,0.9010,0.9505,1 , many numerical experiments have been executed, certain discoveries recorded as follows.

4.1. Varying the fractional order with a=30,b=10,c=1,d=10

Fig. 1 displays the spectrum of Lyapunov exponents versus the order v of fractional derivative and parameters are chosen as a=30,b=10,c=1,d=10. From which one can note that LE1 and LE2 are positive for fractional order v[0.8515,1] . In Table 1, some values of Lyapunov exponents are recoded against different values of v . The proposed system (18) has a hyperchaotic behavior for the fractional orderv[0.8515,1] . Fig. 2 shows the bifurcation map of y1 related to the fractional order v where parameters are chosen as a=30,b=10,c=1,d=10. Fig. 3, Fig. 4, show two hyperchaotic attractors at fractional orders 0.9 and 1 respectively. For v outside this period v[0.8515,1] , the system has fixed point behavior as shown in Fig. 5.
Fig.1 Lyapunov exponents spectrum versus the fractional order v where the model parameters are fixed as a=30,b=10,c=1,d=10.
Table 1 Lyapunov exponents of the proposed system (18) for various values of fractional order v.
Fractional order vLE1LE2LE3LE4LE5LE6
0.85155.93890.0031-0.1850-8.7262-40.1658-57.4287
0.90104.98080.0362-0.0060-6.9942-32.1153-46.8727
0.95053.69360.0006-0.0006-5.3108-25.4198-37.3238
10.35140.0007-0.0068-3.8473-19.9668-27.2099
Fig.2 Bifurcation map of y1 related to the fractional order v where the model parameters are fixed as a=30,b=10,c=1,d=10.
Fig.3 Hyperchaotic attractor of the proposed model (18) at fractional order v=1 anda=30,b=10,c=1,d=10.
Fig.4 Hyperchaotic attractor of the proposed model (18) at fractional order v=0.9 anda=30,b=10,c=1,d=10.
Fig.5 Fixed point behavior of the proposed model (18) at fractional order v=0.8 and a=30,b=10,c=1,d=10.
From above experiments, the proposed system is hyperchaotic for fractional order v[0.8515,1] where model parameters are selected as a=30,b=10,c=1,d=10.

4.2. Vary the parameter a and fix b=10,c=1,d=10

As clear from Fig. 6, the spectrum of Lyapunov-exponents show that the proposed model (18) is hyperchaotic for a wide range of the parameter a for fractional order v[0.8515,1] and the system lose its chaotic dynamical behavior outside this period. Fig. 7 displays the bifurcation maps of y1 related to parameter a where other parameters are selected as b=10,c=1,d=10 at different fractional orders.
Fig.6 Lyapunov exponents spectrum versus the parameter a where other parameters are fixed as b=10,c=1,d=10 at different fractional orders v=1,0.9,0.8,0.7 .
Fig.7 Bifurcation maps of y1 related to parameter a where other parameters are selected as b=10,c=1,d=10 at different fractional orders v=1,0.95,0.9,0.88 .

5. Structural characteristic of proposed model (18)

In this section, we display the relations between the different terms of the proposed model via drawing its graph. The proposed model consists of pure state variables and nonlinear product terms. We draw the suggested graph as follows. The graph vertices are the six pure state variables plus the eight nonlinear product terms. There will be a directed edge between two vertices if the corresponding terms affect the other. For example, there will be a directed edge from the vertex corresponds to the termy2y3y4 to the vertex corresponds to the term y1 , but there no edge from the vertex corresponds to y1 towards that corresponds to y2y3y4 , so model (18) can be displayed by the digraph D given in Fig. 8. Many structural characteristics of the model can be calculated from the adjacency and Hermitian matrices.
Fig.8 The digraph D of the studied chaotic model.
The underlying graph G (shown in Fig. 9) of a digraph D is constructed by neglecting the direction of the edges and the loops (edges with same ends). Following Definitions 2, 3, respectively, the adjacency matrix A(G) for the underlying graph of D and the Hermitian matrix H(D) for the model digraph D are constructed as follows in (19) and (20).
Fig.9 The underlying graph G of the model graph D .
A(G)=y1y2y3y4y5y6y2y3y4y1y3y4y1y2y4y4y5y6y1y2y3y3y5y6y3y4y6y3y4y5[0100001110100010000011101000000000111111110000001111111100000100010111000010000101111111000000000011110000000000111100000000000011110000000011110000000000001111000000000011110000000000111100000000]
H(D)=y1y2y3y4y5y6y2y3y4y1y3y4y1y2y4y4y5y6y1y2y3y3y5y6y3y4y6y3y4y5[110000iii0i000110000iii0i000001000iiiiiiii000100iiiiiiii000011000i0iii000011000i0iiiiiii0000000000iiii0000000000iiii000000000000iiii00000000iiii000000000000iiii0000000000iiii0000000000iiii00000000]
The eigenvalues of A(G) are:
{4.755,2.3723,1,1,0,0,0,0,0,0,0,0.6605,3.3723,5.0945} . And the eigenvalues of H(D) are: {3.2361,2,2,1.2361,3.3723,3,2,2.3723,2,2,3.2361,0,0,1.2361} .
The energy of the graph G is defined as E(G)=i=114|λi|=18.2546 , where 14 is the order of G and λ 's are the eigenvalues of its adjacency matrix. This idea was presented by I. Gutman in [24] and has several applications and it can be utilized to estimate the whole energy of graphs. The Hermitian energy is defined as EH(D)=i=114|μi|=27.689, where 14 is the order of D and μ 's are the eigenvalues of its Hermitian matrix.

6. Electric circuit realization of the studied fractional chaotic system

In this section, we show that the proposed chaotic model is applicable by constructing its electric circuit. By constructing an electric circuit for the proposed model, we prove its ability to be applicable in real world applications. Here, we consider the suggested 6D fractional order chaotic paradigm (18). First, we fix the fractional order v=0.9 to design circuit to realize model (18). We follow the method published in [27] to realize 1/s0.9 . With aid of NI Multisim 14.0 package, the electric circuit given in Fig. 10 is designed, where S1, S2, ..., S6 are six summers all of them has 3 inputs and single output; M1, M2, …, M11 are AD633 voltage multiplier with 2 inputs and single output all of them are adjusted with identity gain; six LM741 op. amp's are utilized to simulate the six 0.9-fractional-order voltage integrators with outputs equal the states of the system y1 to y6 ; R1, R2, R3, R4, R5, R15 are 100KΩ resistances and R6, R7, R8, R9 are 10KΩ resistances plus a couple of LM741 op. amps. are used to construct two voltage inverters. The designed values of resistors and capacitors are given in Tables 2 and 3. In Table 4, we record input/output gain's utilized in voltage summers. The electronic circuit implementing the new chaotic 6D model (18) is shown in Fig. 10. Figs. 11 and 12 display the response of 5.4-fractional-order chaotic 6D new model.
Fig.10 Electronic circuit implementing the 5.4-fractional-order 6D new chaotic model (18).
Fig.11 y4 -y1 phase plane from the electronic circuit of the 5.4-fractional-order 6D new chaotic model (18).
Fig.12 y4 -y2 phase plane from the electronic circuit of the 5.4-fractional-order 6D new chaotic model (18).
Table 2 The tree shape 1/s0.9 (2dB) integrator: values of resistors (c.f. [27]).
v Ra Rb Rc
0.91.55MΩ 61.54MΩ 2.5KΩ
Table 3 The tree shape 1/s0.9 (2dB) integrator: values of capacitors (c.f. [27]).
v Ca Cb Cc
0.90.73μF 0.52μF 1.1μF
Table 4 The used voltage summing devices': gains for inputs (A, B, C) and output O.
summer symbolABCO
S1-30-1300.28
S2-101100.28
S3-1-110.015
S4-1-1100.015
S5-1-30300.28
S6-101-100.28

7. Design of fractional order active controller for the proposed chaotic model (18)

In this section, we propose an active control strategy to control the model states and convert it to a fixed point.
For the proposed system (18), the control rule is designed as follows:
C1=G1y1(a(y2y1)+y2y3y4),C2=G2y2b(y2+y1)+y1y3y4,C3=G3y3cy3y1y2y4y4y5y6,C4=G4y4dy4y1y2y3y3y5y6,C5=G5y5a(y6y5)y3y4y6,C6=G6y6b(y6+y5)+y3y4y5.
where Gi 's are control gains. The controlled model becomes:
Dυy1=a(y2y1)+y2y3y4+C1,Dυy2=b(y2+y1)y1y3y4+C2,Dυy3=cy3+y1y2y4+y5y6y4+C3,Dυy4=dy4+y1y2y3+y5y6y3+C4,Dυy5=a(y6y5)+y6y3y4+C5,Dυy6=b(y6+y5)y5y3y4+C6.
Substituting (21) into (22), the dynamics of the controlled model becomes:
Dυy1=G1y1,Dυy2=G2y2,Dυy3=G3y3,Dυy4=G4y4,Dυy5=G5y5,Dυy6=G6y6,
It is obvious that the controlled system (23) has E^0=(0,0,0,0,0,0) as an equilibrium point. The following theorem can be proved for the controlled fractional order model (23).
Theorem 1. If all controller gains Gi 's in the controlled model (23) are chosen to be negative values, then E^0=(0,0,0,0,0,0) of the controlled model (23) is fixed point.
Proof. Taking Laplace transformation for the controlled model (23), then
then
svYi(s)sv1yi(0)=GiYi(s),
where Yi(s) is the Laplace transformation of yi(t), so
Yi(s)=sv1yi(0)svGi
applying inverse taransformation of Laplace for (24), we get
yi(t)=yi(0)Mv(Gitv),
such that Mε(.) is the common function of Mittag-Leffler. Because 0<v1 , so Mv(Gitv)=h=0(Gitv)hΓ(hv+1) approaches zero for all values of Gitε if all Gi 's are negatives, the reader can consult [28], [29] and its cited references for more about function of Mittag-Leffler. Then, the states yi(t) of the model are all stable asymptotically. Figs. 13, 14 and 15 show the time response of the new fractional 6D chaotic model before and after control at fractional order v=1,0.95and0.9 respectivily with all gains Gi=1 . The controller is applied after 30 secs. We can see that the model takes a transient time to reach its final state. Changing the controller gains to be Gi=5 , the transient time is approaches zero as shown in Figs. 16, 17 and 18.
Fig.13 Dynamics of the new fractional 6D chaotic model before and after control at fractional order v=1 , with all gains Gi=1 .
Fig.14 Dynamics of the newfractional 6D chaotic model before and after control at fractional order v=0.95 , with all gains Gi=1 .
Fig.15 Dynamics of the newfractional 6D chaotic model before and after control at fractional order v=0.9 , with all gains Gi=1 .
Fig.16 Dynamics of the newfractional 6D chaotic model before and after control at fractional order v=1 , with all gains Gi=5 .
Fig.17 Dynamics of the newfractional 6D chaotic model before and after control at fractional order v=0.95 , with all gains Gi=5 .
Fig.18 Dynamics of the newfractional 6D chaotic model before and after control at fractional order v=0.9 , with all gains Gi=5 .

8. Conclusions

The new fractional 6D chaotic model presented above has complicated behaviors and all characteristics of the chaotic system. The proposed fractional 6D chaotic system with different topology and structure than the existing 6D systems is constructed in this work. The basic features and dynamical response of the new model have been investigated to prove it's chaotic. Certain graph theory tools have been utilized to display some hidden properties of the proposed model. In addition, an electronic circuit realizing the proposed model has been designed. Also, an active fractional controller is proposed to control the proposed model. In the future work, we suggest utilizing this fractional 6D model to produce random keys for data encryption and in secure communication applications.
Data Availability Statement
Not applicable.

Author contributions

Conceptualization, M.H.; methodology, M.H., N.A.; software, M.H., S.M., A.A., N.A.; validation, M.H., S.M., A.A., N.A; formal analysis, M.H., S.M., A.A., N.A.; investigation, M.H., S.M., A.A., N.A; resources, M.H., S.M., A.A., N.A.; data curation: M.H.; Funding acquisition, M.H., S.M., A.A., N.A.; writing—original draft preparation, M.H.; writing—review and editing, M.H., S.M., A.A., N.A; visualization, M.H.; supervision, M.H.; project administration, M.H., S.M., All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge the support and funding of Research Center for Advanced Material Science (RCAMS) at King Khalid University through Grant No. RCAMS/KKU/009-21.

Declaration of Competing Interest

The authors declare no conflict of interests.

Acknowledgements

The authors acknowledge the support and funding of Research Center for Advanced Material Science (RCAMS) at King Khalid University through Grant No. RCAMS/KKU/009-21.
[1]
C. Sparrow. The Lorenz Equations: Bifurcations, Chaos, and Strange Attractors. Springer-Verlag, NY (1982)

[2]
OE. Rosslor. Phys. Lett. A, 57 (1976), pp. 397-398

[3]
G Chen, T. Ueta. Int. J. Bifurcat Chaos, 9 (1999), pp. 1465-1466

[4]
T Ueta, G. Chen. Int. J. Bifurcat Chaos, 10 (2000), pp. 1917-1931

[5]
JH Lu G. Chen. Int. J. Bifurcat Chaos, 12 (2002), pp. 659-661

[6]
JH Lu, G Chen. Int. J. Bifurcat Chaos, 12 (2002), pp. 2917-2926

[7]
A Vanecek. Celikovsky's. Control Systems: from Linear Analysis to Synthesis of Chaos. Prentice-Hall, London (1996)

[8]
WB Liu G. Chen. Int J Bifurcat Chaos, 13 (2003), pp. 261-267

[9]
I. Podlubny. Fractional Differential Equations. Academic Press, New York (1999)

[10]

[11]
H. Zhu, S. Zhou, J. Zhang. Chaos Solit. Fract., 39 (2009), pp. 1595-1603

[12]
L. Chen, Y. Chai, R. Wu, J. Yang. IEEE Trans. Circuits Syst. Express Briefs, 59 (9) (Sept.2012), pp. 602-606, DOI: 10.1109/TCSII.2012.2206936.

[13]
E.E. Mahmoud, M. Higazy, T.M. Al-Harthi. Mathematics, 7 (2019), p. 877 DOI: 10.3390/math7100877

[14]

[15]
W.M. Ahmad, J.C. Sprott. Chaos Solitons & Fractals, 16 (2003), pp. 339-351

[16]
R. Gao. Chaos Solitons Fractals, 122 (2019), pp. 236-244

[17]
EE Mahmoud, M Higazy, OA. Althagafi. Symmetry 12 (11) (2020), p. 1765, DOI: 10.3390/sym12111765.

[18]
M. Higazy, Emad E. Mahmoud, E.M. Khalil, S. Abdel-Khalek, S.M. Abo-Dahab, Hammad Alotaibi. complex., 2021 (2021), p. 17, DOI: 6692369

[19]

[20]
Emad E. Mahmoud, M. Higazy, Hammad Alotaibi, 4 S.M. Abo-Dahab, S. Abdel-Khalek, E.M. Khalil. Chaos Solitons Fractals, 144 (2021), Article 110715, DOI: https://www.sciencedirect.com/science/article/pii/S0960077921000680/pdfft?md5=836c3a0eac067a0cd30723798922875d&pid=1-s2.0-S0960077921000680-main.pdf

[21]
S. Das. Functional Fractional Calculus for System Identifcation and Controls. Springer Verlag, Berlin (2007)

[22]
S.G. Samko, A.A. Kilbas, O.I. Marichev. Fractional Integrals and Derivatives:Theory and Applications. CRC Press (1993)

[23]
A. Atangana, D. Baleanu. Therm. Sci., 20 (2) (2016), pp. 763-769

[24]
JX. Liu, XL. Li. Linear Algebra Appl, 466 (2015), pp. 182-207

[25]
G. Chen. GChen, XYu (Eds.), Chaos Control: Theory and Applications, Springer, Berlin (2003), pp. 159-177

[26]
XF. Wang. GChen, XYu (Eds.), Chaos Control: Theory and Applications, Springer, Berlin (2003), pp. 179-204

[27]
W.M. Ahmad, J.C. Sprott. Chaos Solitons & Fractals, 16 (2003), pp. 339-351

[28]
L.A. Quezada-Téllez, S. Carrillo-Moreno, O. Rosas-Jaimes, J.J. Flores-Godoy, G. Fernndez-Anaya. Int. J. Nonlinear Sci. Numerical Simul., 18 (2017), pp. 371-384

[29]
Saxena R. K. Mathai A. M. Haubold, H. J., (2009), arXiv: 0909.0230v2.

[30]
A.M.S. Mahdy, E.S.M. Youssef. Int. J. Modern Phys. C, 32 (01) (2021), Article 2150002, DOI: 10.1142/S0129183121500029.

[31]
Y.A. Amer, A.M.S. Mahdy, R.T. Shwayaa, E.S.M. Youssef. J. Eng. Appl. Sci., 13 (17) (2018), pp. 7388-7394, DOI: 10.36478/jeasci.2018.7388.7394.

[32]
Missing reference

[33]
Y.A. Amer, A.M.S. Mahdy, H.A.R. Namoos. J. Eng. Appl. Sci., 13 (20) (2018), pp. 8489-8493, DOI: 10.36478/jeasci.2018.8489.8493.

[34]

[35]
Zaid Mohammed Ghazi, Imad Shakir Abbood, Farzad Hejazi. J. Ocean Eng. Sci., 7 (1) (2022), pp. 41-57, DOI: https://www.sciencedirect.com/science/article/pii/S2468013321000401/pdfft?md5=f6ec8a115d329fcf634a111a834b50cc&pid=1-s2.0-S2468013321000401-main.pdf

[36]
Mostafa MA Khater, Shabbir Muhammad, A. Al-Ghamdi, M. Higazy. J. Ocean Eng. Sci. (2022), DOI: in press

[37]
Lanre Akinyemi, Udoh Akpan, Pundikala Veeresha, Hadi Rezazadeh. J. Ocean Eng. Sci. (2022), DOI: in press

[38]
M. Higazy, Shabbir Muhammad, A. Al-Ghamdi, Mostafa M.A. Khater. J. Ocean Eng. Sci.(2022), DOI: in press

[39]
Fuzhang Wang, Shabbir Muhammad, A. Al-Ghamdi, M. Higazy, Mostafa M.A. Khater. In (1 + 1)-dimension; inelastic interaction of long-surface gravity waves of small-amplitude unidirectional propagation. J. Ocean Eng. Sci. (2022), In Press

[40]
M. Higazy, A. El-Mesady, A.M.S. Mahdy, Sami Ullah, A. Al-Ghamdi. Numerical, Approximate Solutions, and Optimal Control on the Deathly Lassa Hemorrhagic Fever Disease in Pregnant Women. J. Funct. Spaces, 2021 (2021), Article 2444920,

[41]
Amr M.S. Mahdy, Mahmoud Higazy, Mohamed S. Mohamed. Optimal and Memristor-Based Control of A Nonlinear Fractional Tumor-Immune Model. Comput. Mater. Contin, 67 (3) (2021), pp. 3463-3486,

Outlines

/