Abstract

Memristor is a nonlinear and memory element that has a future of replacing resistors for nonlinear circuit computation. It exhibits complex properties such as chaos and hyperchaos. A five-dimensional memristor-based circuit in the context of a nonlocal and nonsingular fractional derivative is considered for analysis. The Banach fixed point theorem and contraction principle are utilized to verify the existence and uniqueness of the solution of the five-dimensional system. A numerical method developed by Toufik and Atangana is used to get approximate solutions of the system. Local stability analysis is examined using the Matignon fractional-order stability criteria, and it is shown that the trivial equilibrium point is unstable. The Lyapunov exponents for different fractional orders exposed that the nature of the five-dimensional fractional-order system is hyperchaotic. Bifurcation diagrams are obtained by varying the fractional order and two of the parameters in the model. It is shown using phase-space portraits and time-series orbit figures that the system is sensitive to derivative order change, parameter change, and small initial condition change. Master-slave synchronization of the hyperchaotic system was established, the error analysis was made, and the simulation results of the synchronized systems revealed a strong correlation among themselves.

1. Introduction

In the last decade, fractional differential equations started gaining much attention in modeling several real-world problems in different areas including in mathematical epidemiology, physics, engineering, and many others. Fractional-order operators are either with singular kernels such as the Caputo derivative and the Reimann–Liouville fractional derivatives or with nonsingular kernels such as the Caputo–Fabrizio and Atangana–Baleanu fractional derivatives [14].

One of the differences between integer- and fractional-order derivatives is that the integer-order derivative describes local properties of a certain dynamic system, whereas the fractional-order derivative representation of a dynamic system involves the whole space of the process [5]. That is, applying fractional derivative orders in modeling real-world problems is essential for describing the hereditary specifications and effectiveness of the memory as essential aspects of different mechanisms in the problem [6, 7].

The recent increase in the study of different dynamical systems using fractional-order derivatives is attributed to the fact that most of the dynamic systems associated with complex systems are found to be nonlocal involving long memory in time and intrinsically fractional derivative operators can describe such systems more accurately than the integer derivatives [8]. In other words, important features of many physical systems are best described or exposed by using fractional-order operators.

There are different studies conducted using fractional derivatives for diverse types of systems pertinenet to different problems. For instance, stability analysis of fractional differential equations with unknown parameters was considered using the D-decomposition method by Koksal [9]. Time- and frequency-domain response of the RLC circuit with fractional-order derivative was investigated in [10]. Telegraph equation used for the power transmission line with a nonlocal boundary value problem was deliberated with a different type of finite difference numerical schemes in [11].

Chaos theory has attracted many researchers and has applications in the fields of encryption and secure of communication [12], modeling financial systems, representing circuit diagrams, and many others [13, 14].

Nowadays, there are many pieces of literature dedicated to the analysis of chaotic systems using fractional derivative operators. This includes chaotic systems of Chua’s electrical circuits and memristor-based circuit systems. Some of the literature are reviewed as follows.

Sene applied the Caputo fractional derivative operator in detecting the chaotic behavior of different 3D and 4D chaotic systems. He used Lyapunov exponent and bifurcation diagrams to identify the nature of the chaos and impact of parameter variation for the different chaotic models investigated [1518]. Different chaotic systems including Chua’s electric circuit and several other chaotic systems were analyzed using fractional derivative order mathematical models by Petráš [19]. Bifurcation and chaotic behaviors in fractional-order simplified Lorenz system using Adams–Bashforth–Moulton predictor-corrector scheme were considered in [20]. Atangana–Baleanu fractional derivative operators were used for modeling and analysis of different chaotic and hyperchaotic systems, and solutions were approximated using a two-step Adams–Bashforth numerical scheme in [21].

It was in 1971 that circuit theorist Chua proposed memristor as a missing two-terminal nonlinear electrical component. The three basic components of a circuit are resistor, capacitor, and inductor. Memristor known for its memory effect and nonlinear characteristics is the fourth circuit component. Memristor relates magnetic flux and electric charge linkage in which case it is called a charge-controlled electric model or it models a relationship between charges and flux in which case it is called a flux-controlled model [12, 22].

At present, there are several studies conducted on memristor-based chaotic circuits using both integer and fractional-order derivatives. In [23], a conformal fractional-order simplest memristor-based chaotic circuit was investigated based on conformable Adomian decomposition method, Lyapunov exponent, bifurcation diagram, and Poincare sections. Buscarino et al. introduced a chaotic circuit based on a realistic model of the HP memristor, and numerical results showed a generation of chaotic attractors [22]. A novel 5D chaotic system with flux-controlled memristor and integer-order derivative, extracted from Wang’s 4D hyperchaotic system, was proposed by Wang et al. [24]. A memristor-based chaotic circuit modified by replacing the nonlinear resistor in Chua’s circuit with a flux-controlled memristor was analyzed in [12].

This research aims to study the memory effect properties and detection of chaos in a five-dimensional memristor-based system. Accordingly, an integer model memristor-based circuit is represented by Atangana–Baleanu fractional derivatives in the Caputo sense (ABC), and the existence and uniqueness of solution of the ABC fractional model are analyzed based on Banach fixed point theorem for contraction principle. From the different concepts of fractional-order operators introduced above, the purpose of choosing the ABC fractional derivative is due to the fact that it possesses nonlocal kernels and of course it allows the inclusion of traditional initial conditions in the formulation of a mathematical model. Numerical approximation of the ABC fractional model is made using the newly developed numerical approximation for fractional derivative developed by Toufik and Atangana in [25]. Local stability of the fractional model representation is accomplished using the Matignon stability criterion. The existence and nature of chaos in the fractional model are investigated using Lyapunov exponents. A bifurcation diagram for different fractional derivatives and parameter variation is given. Several phase portraits are depicted as a verification for the impact of different parameter values and different fractional derivative orders. The impact of initial conditions on the solution trajectory of the system is also investigated using simulation of the trajectories of the system for different initial conditions. Lastly, master-slave synchronization of the system is performed accompanied by the corresponding error analysis and simulation of several cases of the synchronized system. All the phase portraits and solution trajectories in this work were obtained from the numerical scheme of Toufik and Atangana adapted for the memristor chaotic model considered in this work. A computing software application called Matlab 2019a is used for the simulation of different results.

A Matlab code for Lyapunov exponents and bifurcation diagram of fractional-order systems called Danca algorithm [26] is used to quantify the chaos by calculating Lyapunov exponents and obtaining bifurcation diagrams for different fractional orders and different parameter values of the model. Some of the evidence for the originality of this work includes the application of Atangana–Baleanu fractional operator to the memristor-based system considered in this study, application of the newly developed numerical approximation by Toufik and Atangana for the fractional-order systems, obtaining the phase portraits of the system from the numerical scheme, and performing synchronization of the five-dimensional system using the numerical approximation.

The remaining part of this paper is arranged as follows. In Section 2, the memristor-based circuit model considered in this study is described in the context of integer derivative. Section 3 is devoted to the fractional derivative representation of the memristor-based systems following some recapping of preliminary concepts and definitions of Atangana–Baleanu fractional derivatives. The existence and uniqueness of the solution for the fractional derivative representation of the model are portrayed in Section 4 of the paper. The numerical scheme applied to get all the phase portraits and approximate time-series solution of the memristor-based system are developed in Section 5. Section 6 is concerned with the local stability analysis of the fractional model followed by Lyapunov exponents, bifurcation diagrams with different fractional orders, and parameter variation in Section 7. Investigation of the impact of small change in the initial conditions on the dynamics of the system is considered in Section 8 followed by synchronization of the hyperchaotic system in Section 9. Finally, in Section 10, conclusions are given followed by list of references.

2. Mathematical Model Description of Memristor-Based Circuit

Wang et al. [27] developed a new four-dimensional hyperchaotic circuit system using an improved modularized design and proposed implementation of the circuit. The hyperchaotic circuit system mentioned in [24] is described under an integer-order system of differential equations. Later on, Wang et al. [24] proposed a five-dimensional flux-controlled memristor-based circuit system and demonstrated that the system exhibits hyperchaotic character under an integer-order system of differential equations. The five-dimensional nondimensionalized flux-controlled memristor-based circuit developed by Wang et al. is described using integer-order differential as given in the following equation:where the memductance used in this study is

In this study, motivated by the contribution of Wang et al. [24] in developing the five-dimensional memristor-based circuit, we analyzed system (1) under ABC fractional derivatives.

3. Fractional Derivative Representation of the Model

In this section, we recall the definitions and basic properties of nonsingular and nonlocal Atangana–Baleanu fractional derivative of Caputo type (ABC) used for our analysis of the circuit.

Definition 1. Let and let . The AB fractional derivative of order is defined as [7, 28, 29]where and the Mittag-Leffler function is

Definition 2. The Atangana–Baleanu (AB) fractional integral of the function is given by [7, 28, 29]

Lemma 1. The AB fractional derivative and AB fractional integral of fulfill [30]

Lemma 2. For , the AB fractional derivative in the Caputo sense satisfies the Lipschitz condition [29]:

Now we continue with the reformulation of (1) in terms of ABC fractional derivatives.

Lemma 3. Let . Then, , admits one solution given by

We can describe the result in Lemma 3 in the form of the Banach fixed theorem point as follows. We begin by defining a Banach spacewith a norm defined as .

Now define the operators by

We are now ready to describe the dynamic equation for the memristor-based chaotic circuit given in (1) using ABC fractional derivative:

The initial condition is given as and .

4. Existence Theory on the Model

Here, the existence and uniqueness of the solution for the ABC model given in (15) are shown using Banach fixed point theorem for contraction mapping. The following two theorems are worth recalling before proceeding further.

Theorem 1. Let be any nonempty closed subset of a Banach space . Then, any contraction has a unique fixed point [31, 32].

Theorem 2. Assume that are continuous functions satisfying the following conditions:(C1) There exist constants such that .(C2) , , and ; then, the ABC fractional derivative system given by (15) has a unique solution in the region .

Proof. Let us show that the operator defined in (10) is well defined in the sense that and is continuous on , wherefor .
Now for any and from (10), we haveConsequently, we have .
To show continuous differentiability on we proceed from which is continuous on , and then we conclude that is continuous on ; as a result, .
To show that the operator has a fixed point based on Theorem 1, it is enough to show that is a contraction mapping. Indeed, let , and it is not difficult to show thatwhere . Since by hypothesis of Theorem 2, we conclude that is a contraction mapping.
To show that the remaining operators , and are contraction mappings, we proceed in the same manner.
Let us show that the operator defined in (14) is well defined in the sense that and is continuous on , wherefor .
Now for any and from (11), we haveConsequently, we have . To show continuous differentiability on , we proceed from which is continuous on , and we conclude that is continuous on ; as a result, .
To show that the operator has a fixed point, we apply Theorem 1. Following the theorem, it is sufficient to show that is a contraction mapping. Indeed, let . Then,where Since , we conclude that is a contraction mapping.
Let us show that the operator defined in (12) is well defined in the sense that and is continuous on , wherefor .
Now for any and from (12), we haveConsequently, we have . To show continuous differentiability on , we proceed from which is continuous on , and we conclude that is continuous on ; as a result, .
To show that the operator has a fixed point, we apply Theorem 1. Following the theorem, it is sufficient to show that is a contraction mapping. Indeed, let . Then,where Since by hypothesis , we conclude that is a contraction mapping.
We have then proved that the operators , and are well defined and are contraction mappings. The case for the operators and follows immediately. Hence, by the Banach fixed point theorem, system (17) has a unique solution in .

5. Numerical Solutions

In this part, the numerical scheme applied to get the phase portrait of the dynamic system (15) is introduced. In the context of chaotic or hyperchaotic fractional differential equations, the use of analytical methods such as the Sudumu transform method, the Laplace transform method, the homotopy analysis method, and the homotopy perturbation method cannot easily be applied because of the nonlinearities of the system [13]. This leads to the need for using numerical methods to approximate the solutions of systems of fractional differential equations. Some of the numerical methods that can be applied for this case include Adams–Bashforth and Toufik–Atangana numerical schemes [8]. Both of these methods are based on Lagrange interpolation polynomials. In this study, the newly developed numerical approximation for fractional derivatives by Toufik and Atangana is employed. The numerical scheme is particularly developed for approximation of Atangana–Baleanu fractional derivative considered in this study, and it is proved to be convergent, stable, and consistent [8].

For convenience, let us write (15) in the following form:where

Now from Lemma 3 and the first equation of (25), we have

The solution for (27) is given as follows:

Applying Lagrange’s interpolation polynomial on to equality leads towhere .

Substituting (28) into (29), we obtainwhere

Substituting in (31) and (32) results in

The expression (30) can be expressed in terms of (31) and (34) as shown in the following equation:

Similarly, we get the following equations for the rest of the state variables:

6. Local Stability Analysis

In this section of the study, the local stability analysis of the fractional model represented in (25) is performed. It is known that the equilibrium points of chaotic systems are not generally stable. Some of the standard methods of stability analysis in fractional calculus are the Matignon criterion and the Laplace transform methods. In this work, the Matignon method is used for its simplicity and is most commonly used in the literature for the same purpose [1, 7, 13, 17].

The Matignon criterion is given bywhere represents the Jacobian matrix, is the set of the eigenvalues of , and is the fractional-order derivative. In the context of fractional derivative, an equilibrium point of (25) is said to be locally stable provided that the Matignon criterion (40) is satisfied for each of the eigenvalues of the Jacobian matrix.

To determine if the equilibrium points of (25) are stable or not, we proceed as follows:(i)The equilibrium points: the equilibrium points of (25) are given as , and line equilibrium is given by .(ii)The Jacobian matrix: in this study, analysis is made regarding the trivial equilibrium point for simplicity. The parameter values used, unless otherwise mentioned, are set as , . Thus, the Jacobian matrix of (25) evaluated at the trivial equilibrium point is given by

The eigenvalues of the Jacobian matrix are given by , , , , and . To check the local stability of the trivial equilibrium point, using the Matignon criteria, we need to verify if all the eigenvalues satisfy condition (40). Indeed, which is not possible since . is true for all . It can then be concluded that the equilibrium point is locally unstable for the parameter values considered.

Furthermore, since one of the eigenvalues of the Jacobian matrix is with positive real part, it can be inferred that system (25) satisfies the necessary condition for showing the double scroll attractor [19].

7. Lyapunov Exponents, Bifurcation, and Chaos with Different Fractional Orders of and Different Parameter Values

In this section, the level of chaos in system (25) is quantified using the Lyapunov exponent method. Bifurcation diagrams of the system (25) related to the fractional derivative order and three of the parameters in the model named, , and , are depicted.

A Matlab code for Lyapunov exponents of fractional-order systems named the Danca algorithm [26] is used to quantify the chaos by calculating Lyapunov exponents for different fractional orders of model (25). The initial conditions used in this part of the work are given by . The parameter values used are similar to the ones used above for calculating the Jacobian matrix. The corresponding Lyapunov exponents (LEs) for different fractional orders are shown in Table 1.

It is then possible to conclude that the 5D system considered in this study is dissipative since the sum of the LEs in each column of Table 1 is negative and the system exhibits a hyperchaotic behavior since there are at least two positive LEs in each column of the table and the largest LE is positive. Moreover, Kaplan–Yorke dimension coresponding to the fractional derivatives considered in Table 1can be calculated. For instance, the dimension of two of the fractional orders is given as follows.

For q = 0.99,

For q = 0.98,

7.1. Bifurcations due to Variation of the Fractional Order

For obtaining bifurcation diagrams due to the variation of the fractional order , the values of all the parameters are kept fixed and the order of the fractional derivative is varied in the interval with an increment of . The other parameter values used in this simulation are , and the initial condition is . The bifurcation diagram is shown in Figure 1.

As shown in Figure 1, when , the system is hyperchaotic. As a verification of the observation, the phase portraits corresponding to q = 0.98 and q = 0.99 using the numerical approximations for fractional-order systems indicated by equations (35)–(39) are shown in Figures 2 and 3, respectively. In both Figures 2 and 3, the dependence of the hyperchaotic system on the fractional derivative is observable as the two figures are seen to be different from each other for a fractional order difference of 0.01. If we consider the trajectories in Figure 3, the orbit due to q = 0.98 converged in less than 13.5 seconds which is not the case for q = 0.99. That means different fractional derivative orders generate reasonably different dynamics of the system (25).

7.2. Bifurcations due to Variation of the Parameter

The values of the parameters used are , the fractional derivative order is , and the initial condition is . The value of the parameter is made to vary in the interval (29.5, 30.5) with an increment of 0.001. The bifurcation diagram is shown in Figure 4. It can be observed from the diagram that system (25) exhibits significant hyperchaos throughout the interval under consideration. In support of this conclusion, some of the phase portraits of hyperchaotic system (25) are depicted in Figure 5 projected on different planes.

It can be inferred from Figures 5(a)–5(d) that the orbits seem to overlap initially, but as time increases, their divergence from each other increases. That is, as time goes, the orbit due to showed more contraction than the trajectory due to ; this is because of the sum of all the Lyapunov exponents (LEs) of the respective parameter values considered in this case. The sum of LEs for is −22.9526 (the LEs are shown in Table 1), the sum of LEs for is −23.4693 (LEs are 15.3796, 0.7298, 0.0051, −8.4066, and −31.1772), and the sum of LEs for is −23.2106 (LEs are 15.3817, 0.7298, 0.0051, −8.4066, and −30.9209). Since the system under investigation is dissipative, the contraction rate of volumes due to is stronger than the contraction rate of volumes due to which is inferred from the magnitude of the sum of all LEs for the corresponding parameter values.

7.3. Bifurcations due to Variation of the Parameter

The values of the parameters used are , the fractional derivative order is , and the initial condition is . The value of the parameter is made to vary in the interval [7, 8.5] with an increment of 0.001. The bifurcation diagram is shown in Figure 6. It can be observed from the diagram that system (25) exhibits significant hyperchaos throughout the interval under consideration. In support of this conclusion, some of the phase portraits of the hyperchaotic system (25) are depicted in Figure 7 projected on different planes and some of the time-series trajectories are depicted in Figure 8.

It can be inferred from Figures 7(a)–7(c) that the orbits seem to overlap initially, but as time increases, their divergence from each other increases. That is, as time goes, the orbit due to showed more contraction than the trajectory due to , as can be observed from the sum of all LEs of the respective parameter values considered in this case. The sum of LEs for is −22.9526 (the LEs are shown in Table 1) and the sum of LEs for is −22.4277 (LEs are 15.3837, 0.7298, −7.8817, 0.0050, and −30.6645). Since the system under investigation is dissipative, the contraction rate of volumes due to is stronger than the contraction rate of volumes due to which is inferred by the magnitude of the sum of all LEs corresponding to the parameter values.

The trajectories in Figure 8 seem to overlap for the first few seconds and then begin diverging from each other; the divergence increases with time.

8. Impact of Initial Condition

In this section, the impact of different initial conditions on the dynamics of the system (25) is addressed. It is well known that one of the properties of chaotic systems is sensitivity to initial conditions, and thus it seems relevant to verify the impact of applying different initial conditions on the phase-space and time-series solutions of system (25). Accordingly, the parameter values and derivative order used are , and . In Figure 9, some of the phase portraits of hyperchaotic system (25) are depicted corresponding to initial conditions (0.11,0.11,0.11,0.11.0.11) by varying initial coordinates of is made to run for 5 seconds with a time step of h = 0.001. The simulation result is shown in Figure 9.

As can be observed from the figures, hyperchaotic system (25) is sensitive to changes in initial conditions. As the variation in the initial condition increases, the dynamics of the system become different from each other at least for this example.

As can be observed from Figures 9 and 10, a small difference of the initial condition generates a significantly observable change in the dynamics of system (25).

9. Synchronization of the Hyperchaotic Model

Synchronization between two chaotic systems is one of the most interesting phenomena in the study of dynamic systems. Synchronization is an occurrence in which two or more chaotic/hyperchaotic systems express a strong correlation among themselves. It is very interesting to see that two hyperchaotic systems, being highly sensitive to initial conditions and exponential divergence of nearby orbits, get synchronized starting from two different initial conditions. Pecora and Caroll [33, 34] in 1990 showed for the first time the possibility of synchronizing two chaotic systems starting from two different initial conditions. Pecora and Caroll used a technique called replacement synchronization. Chaotic synchronization has a range of applications including securing communications based on phase synchronization and controlling insulin production by beta cells. In this study, we applied the master-slave synchronization technique. That is, we considered two identical copies of system (25) and related them with a coupling function. We considered different initial conditions for the coupled systems, showed that the error dynamics are asymptotically stable, and portrayed some of the simulation results. The coupling function is added to the slave system to make it respond to the master system.

Now the master hyperchaotic fractional model is given by

Let the slave model be given by

Define the error terms by

The error dynamics are then given by

The vector is chosen aswhere

The choice of the matrix H is arbitrary except that all of its eigenvalues must satisfy the Matignon criteria of stability (40).

Then, the error dynamics become

The eigenvalues of matrix H are −40.8876, −13.0000, −12.1143, −4.9900, and −3.0081, and all the eigenvalues satisfy Matignon criteria, and thus the error dynamics are asymptotically stable.

Now the slave system is given by

The numerical simulation that verified a strong correlation between the master and the slave system in the form of phase space, time-series orbits, and the error dynamics are depicted in Figures 1113. The numerical approximation for fractional order system (25) given in equations (35)–(39) is adapted to simulate each of the error dynamics (50), the master system (44), and slave system (51). The simulation results are shown in Figures 1113. The parameter values and derivative order used are , and . The initial condition used for master systems is (0.11, 0.11, 0.11, 0.11, 0.11), and it is (−10, 10, −5, −0.5, −0.5) for the slave system. The error graphs indicate fast synchronization of the slave and master systems.

10. Conclusion

In this study, a 5-dimensional memristor-based hyperchaotic circuit in the context of the fractional operator was considered and different qualitative and quantitative analyses are made: numerical approximation of the solution, bifurcation diagrams, Lyapunov exponents for different fractional orders, and parameter values are used to investigate the nature of the solution and chaotic system. The system is found to be hyperchaotic. Sensitivity to initial conditions and sensitivity to parameter value changes were found to cause a significant effect on the dynamics of the system. Furthermore, in the course of making several simulation results on the dynamics of the system, the hyperchaotic system is sensitive to almost any change made in the system including time scale increment change used in the study. The sensitivity of the system to fractional derivative order is very strong in the sense that the simulation results for q = 0.98 and 0.99 are significantly different (see Figures 2, 3, and 14). This is perhaps related to the contraction rate of volumes due to the different fractional-order derivatives. The magnitude of the sum of all the LEs corresponding to different fractional-order derivatives decreases as the order of the derivative increases for the case of the hyperchaotic system considered in this study (see Table 1). This requires further investigation. Synchronization of the system is also established and exhibited a strong agreement. The dependence of synchronization on the fractional derivative orders and variation of parameters needs further investigation.

Data Availability

No datasets were generated or analyzed during the current study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

The authors declare that the study was realized in collaboration with equal responsibility. All authors read and approved the final manuscript.

Acknowledgments

The first and third authors would like to thank Azarbaijan Shahid Madani University.