COMPLEXITY Complexity 1099-0526 1076-2787 Hindawi 10.1155/2018/9467435 9467435 Research Article Dynamic Behaviors and the Equivalent Realization of a Novel Fractional-Order Memristor-Based Chaotic Circuit http://orcid.org/0000-0002-7219-484X Yang Ningning12 http://orcid.org/0000-0003-4782-8350 Xu Cheng2 http://orcid.org/0000-0002-3644-4354 Wu Chaojun3 http://orcid.org/0000-0003-3497-6506 Jia Rong12 Liu Chongxin4 Hens Chittaranjan 1 State Key Laboratory Base of Eco-Hydraulic Engineering in Arid Area Xi’an University of Technology Xi’an 710048 Chinaxaut.edu.cn 2 Institute of Water Resources and Hydro-Electric Engineering Xi’an University of Technology Xi’an 710048 Chinaxaut.edu.cn 3 College of Electronics and Information Xi’an Polytechnic University Xi’an 710048 Chinaxpu.edu.cn 4 School of Electrical Engineering Xi’an Jiaotong University Xi’an 710049 Chinaxjtu.edu.cn 2018 2122018 2018 30 05 2018 26 08 2018 17 09 2018 2122018 2018 Copyright © 2018 Ningning Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper proposed a novel fractional-order memristor-based chaotic circuit. A memristive diode bridge cascaded with a fractional-order RL filter constitutes the generalized fractional-order memristor. The mathematical model of the proposed fractional-order chaotic circuit is established by extending the nonlinear capacitor and inductor in the memristive chaotic circuit to the fractional order. Detailed theoretical analysis and numerical simulations are carried out on the dynamic behavior of the proposed circuit by investigating the stability of equilibrium points and the influence of circuit parameters on bifurcations. The results show that the order of the fractional-order circuit has a great influence on the dynamical behavior of the system. The system may exhibit complicated nonlinear dynamic behavior such as bifurcation and chaos with the change of the order. The equivalent circuits of the fractional-order inductor and capacitor are also given in the paper, and the parameters of the equivalent circuits are solved by an undetermined coefficient method. Circuit simulations of the equivalent fractional-order memristive chaotic circuit are carried out in order to validate the correctness of numerical simulations and the practicability of using the integer-order equivalent circuit to substitute the fractional-order element.

Xi’an Science and Technology Bureau201805037 (21)Key Project of Natural Science Basic Research Plan in Shaanxi Province of China2018ZDXM-GY-169Natural Science Foundation of Shaanxi Province2018JM5068Shaanxi Provincial Department of Water Resources2017slkj-15National Natural Science Foundation of China51507134
1. Introduction

Since memristor was first postulated by Leon O. Chua in 1971 , the research on it has been a hot topic in research. Memristor is a nonlinear passive two-terminal component and represents the relationship between a charge and flux . As the fourth fundamental circuit element besides the resistor, inductor, and capacitor, memristor has many unique features that other elements do not have. It can remember its history resistance so be known as a contraction of “memory resistor.” The resistance of memristor depends on the magnitude, direction, and duration of the voltage applied to it . In recent years, memristor has been widely used in many fields such as resistive random access memory (RRAM) [4, 5], neural networks , signal processing , chaotic and control systems , and sample recognition [7, 11].

The fingerprint of the memristor is that the loci in the voltage-current plane exhibit a pinched hysteresis loop which always passes through the origin when driven by any periodic input voltage source or current source . Based on the above statement, a class of diode bridge circuits can be described as generalized memristors because they have been proved to satisfy the fingerprint of a memristor . It is proved that an elementary electronic circuit consisting of a full-wave rectifier and a second-order RLC filter has memory properties . Removing the resistor, another diode bridge circuit cascaded with a second-order filter containing an inductor and a capacitor can also constitute a generalized memristor . Replacing the second-order filter with a first-order parallel RC filter or a series first-order RL filter, a newly generalized memristor is realized, which reduces to a nonlinear fundamental element . Recently, an improved memristive diode bridge circuit was put forward comprising four diodes and an inductor, which has a much simpler circuit realization [18, 19].

Due to the nonlinear characteristics of the memristor, chaotic oscillations can easily emerge. Replacing the Chua’s diode in the canonical Chua’s oscillator with a memristor, a memristor-based chaotic circuit was implemented in 2008 , which is known as the first combination of memristor and chaotic circuit. Since then, the study of memristor chaotic circuits has been widely undertaken . A simple memristor chaotic circuit is proposed in paper , in which the inductor and capacitors are connected in series and one capacitor cascaded with a memristor. In paper , a memristor, an inductor, a capacitor, and a linear negative resistor constitute a chaotic circuit, and paper  puts forward a simplest chaotic circuit which has only three circuit elements: a linear passive inductor, a linear passive capacitor, and a nonlinear active memristor.

On the other hand, fractional calculus, a more than 300-year-old mathematical subject, can describe a real object more accurately than the classical “integer” methods . In addition, due to the introduction of extra adjustable parameters, the fractional-order model will increase the degree of freedom. What is worth mentioning is that the principle of fractional calculus is based on the memory property of the fractional-order integral or derivative , so the connection between fractional calculus and memristor is straightforward, the memristor can be extended to fractional order as well. In paper , a fractional-order memristor model was applied to the Chua’s oscillator for the first time, which constructed a fractional-order memristor-based chaotic circuit. Paper  proposed a fractional-order generalized memristor which consists of four diodes cascaded with a fractional-order parallel RC filter (the order is between 0 and 1) and established a fractional-order chaotic circuit by replacing the traditional Chua’s diode with the proposed memristor. The fractional-order capacitor and inductor are implemented by the corresponding unit circuits. Not only numerical simulations were carried out but also equivalent circuits were realized in the paper.

This paper is organized as follows. Section 2 briefly introduces the fundamental of fractional calculus and constructs a generalized memristor based on a fractional-order series RL filter and then establishes a fractional-order memristive chaotic circuit. In Section 3, detailed theoretical analysis and numerical simulations are carried out on the dynamic behavior of the proposed circuit by investigating the stability of equilibrium points and the influence of system parameters on bifurcations. In Section 4, the equivalent circuits of fractional-order inductor and capacitor are implemented, and circuit simulations of the realized fractional-order memristive chaotic circuit are carried out. Section 5 concludes this article with some additional remarks.

2. Fractional-Order Memristive Chaotic Circuit 2.1. Fundamentals of Fractional Calculus

Fractional calculus, as a generalization of integration and differentiation to integer order, can more accurately characterize the properties of actual objects. The fundamental integral-differential operator Dntλ is defined as [32, 33] (1)Dntλ=dλdtλ,λ>0,1,λ=0,ntdτλ,λ<0,where n and t are the bounds of the operation and λ is the order.

There are three equivalent definitions used for the general fractional differentiation known as the Grünwald-Letnikov (GL) definition, the Riemann-Liouville (RL) definition, and the Caputo definition. The Caputo definition is more convenient for initial condition problems because its physical meaning is clear . In this work, we will use the differential conformal transformation based on Caputo’s derivatives. The Caputo definition can be written as (2)DnCtλft=1Γmλntfmτtτλm+1dτ,m1<λ<m,where Γ is the gamma function and λ is usually greater than zero . The main advantage of Caputo’s approach is that the initial conditions for fractional differential equations with Caputo derivatives are in the same form as for integer-order differential equations, and the Laplace transform of the Caputo derivative allows utilization of initial values of classical integer-order derivatives with known physical interpretations . The Laplace transform of the λ order Caputo differential operator is (3)LD0Ctλft=sλFsk=0m1skdλ1kftdtλ1kt=0,where m is an integer and λm1,m. In this paper, we consider the zero initial conditions, so Laplace transform of Caputo’s fractional derivative reduces to  (4)LD0Ctλft=sλFs.

In this case, the three well-known definitions of fractional derivatives are equivalent .

2.2. Fractional-Order Memristor Model

Paper  has demonstrated that a diode bridge circuit cascaded with an inductor and resistor in series satisfied the definition of a generalized memristor, whose equivalent circuit only contains six fundamental elements, which is depicted in Figure 1. Assuming the voltage across the generalized memristor is um, the current flowing through the generalized memristor is im, and the current passing through the inductor Lm is iL, the mathematical model of the proposed generalized memristor circuit is established as (5)diLdt=umLmRmiLLm1ρLmln2IS+iLeρum2IScoshρum,im=2IS+iLtanhρum,where IS and ρ are diode parameters and ρ=1/2nVT. IS, n, and VT indicate the reverse saturation current, emission coefficient, and the thermal voltage of the diode, respectively. In this paper, diodes 1N4148 are used for circuit simulations, and parameters of the diodes are presented as IS=2.682 nA, n=1.836, and VT=25 mV.

Generalized memristor and its circuit realization: (a) memristive diode bridge with series RL filter and (b) memristor symbol representation.

As nonlinear elements, inductors and capacitors can be extended to fractional order from the point of view of fractional calculus . As a result, a fractional-order inductor-based generalized memristive circuit is realized by replacing the integer-order inductor Lm with a fractional-order inductor Lmλ and whose mathematical model can be expressed as (6)dλiLdtλ=umLmλRmiLLmλ1ρLmλln2IS+iLeρum2IScoshρum,im=2IS+iLtanhρum.

In order to verify whether the proposed fractional-order generalized memristor satisfied the three fingerprints of a memristor , numerical simulations are carried out using MATLAB software. The input voltage is given as a sinusoidal voltage um=Amsin2πft, and the parameters of the inductor and the resistor are 280 mH and 580 Ω, respectively. Setting the fractional order λ equals 0.95, the relationship curves of the memristor input voltage and input current versus different frequencies are indicated in Figure 2(a), and the relationship curves with different amplitudes are demonstrated in Figure 2(b). The amplitude of input voltage in Figure 2(a) is set as Am=4 V, and the frequency of input voltage is f=500 Hz in Figure 2(b). Moreover, the hysteresis loops with different orders are shown in Figure 2(c), in which theamplitude and frequency are set as Am=4 V and f=500 Hz, respectively. It can be seen that the hysteresis loops pinched at the origin and the area of the pinched hysteresis loops decrease monotonically as the input voltage frequency increase when the frequency is greater than the critical value. Besides, the pinched hysteresis loop will shrink to a nonlinear single-valued function when the frequency tends to infinity. The areas of pinched hysteresis loops increase as the input voltage amplitude increase as well, as depicted in Figure 2(b). In addition, the pinched hysteresis loops indicated in Figure 2(c) demonstrate that the order of the fractional-order inductor can affect the dynamic characteristics of the memristor. It can be seen that the pinched hysteresis loops with different amplitudes are similar to each other (see Figure 2(b)), but the frequency and the order can bring about an important impact on the memristive characteristics (see Figures 2(a) and 2(c)).

The hysteresis loops of the fractional-order memristor with different sinusoidal input voltages: (a) Am=4V and λ=0.95 with different frequencies, (b) f=500 Hz and λ=0.95 with different amplitudes, and (c) Am=4 V and f=500 Hz with different orders.

2.3. Fractional-Order Memristor-Based Chaotic Circuit

After the realization of the fractional-order generalized memristor, we consider applying the proposed first-order generalized memristor to a classical Chua’s circuit. By replacing the traditional Chua’s oscillator with a fractional-order inductor-based memristor, a fractional-order memristive chaotic circuit is established, which is indicated in Figure 3. Since the practical inductor and capacitor can be extended to fractional order to better express its electricity characteristics, the capacitance and inductance used in this paper are all fractional order. Therefore, the circuit is composed of six elements which are a fractional-order inductor Lλ, two fractional-order capacitors C1λ and C2λ, a fractional-order memristor Mλ, a linear passive resistor, and a negative conductance. Because the proposed fractional-order memristive Chua’s circuit has four dynamic elements, there are four corresponding state variables, which are the voltages across the fractional-order capacitors C1λ and C2λv1 and v2, the current flowing through the fractional-order inductor Lλi3. The last one is the current passing through the fractional-order inductor Lmλ in the fractional-order generalized memristor, MλiL, as shown in (6), which can express the dynamic characteristics of the fractional-order generalized memristor Mλ. Utilizing Kirchhoff voltage law (KVL) and Kirchhoff current law (KCL), the characteristic equation of the fractional-order memristive Chua’s circuit can be derived as (7)dλv1dtλ=v2v1RC1λi3C1λ,dλv2dtλ=v2RG1RC2λ+v1RC2λ2IS+iLtanhρv2C2λ,dλi3dtλ=v1Lλ,dλiLdtλ=v2LmλRmiLLmλ1ρLmλln2IS+iLeρv22IScoshρv2.

Fractional-order memristive Chua’s circuit.

The model of the fractional-order memristive Chua’s circuit is established in MATLAB, and simulations are carried out in numerical. Set capacitance C1λ=100 nF, capacitance C2λ=10 nF, inductance Lλ=12 mH, resistance R=2 kΩ, and conductance G=0.6667 mS, and the parameters of the fractional-order generalized memristor are the same as that mentioned above. When the fractional order λ is taken as 0.98 and the initial values of the state variables are selected as v1=0 V, v2=0.01 V, i3=0 A, and iL=0 A, the phase diagrams of the fractional-order memristive Chua’s circuit are exhibited in Figure 4.

Phase diagrams of the fractional-order memristive Chua’s circuit. Phase diagrams in the (a) v2tv1t plane, (b) v2ti3t plane, (c) v2tiLt plane, and (d) v2timt plane.

3. Chaotic Dynamic Behavior of the Fractional-Order Memristive Chua’s Circuit 3.1. Equilibrium Points and Stability

In this section, equilibrium points and eigenvalues of the corresponding Jacobian matrix are calculated to qualitatively analyze the dynamic behavior of the proposed fractional-order memristive chaotic circuit. The equilibrium points of such a circuit can be obtained by solving the following equations: (8)v1=0,v2RG1=R2IS+iLtanhρv2,v2=RmiL+1ρln2IS+iLeρv22IScoshρv2,v2=i3R.

Obviously, O1=0,0,0,0 is one equilibrium point of the fractional-order memristive Chua’s circuit, but the analytic solutions of v2 and iL cannot be computed in MATLAB, so we consider using the graphic analytic method to obtain the other equilibrium points . The relationship between v2 and iL can be written as (9)v2=±1ρacosh2IS+iLeρRmiL2IS,iL=RG1v2Rtanhρv22IS.

Utilizing (9), we can draw two functions between v2 and iL; then the intersection coordinates of two functions can easily be obtained, as shown in Figure 5. Thus, the other two equilibrium points are obtained as O2=0,1.134,0.000567,0.000189 and O3=0,1.134,0.000567,0.000189.

Two functions and their intersection points.

The Jacobian matrix of the characteristic (7) can be described as (10)1RC1λ1RC1λ1C1λ01RC2λΛ10tanhρv2C2λ1Lλ0000Λ20RmLmλ1ρLmλ2IS+iL,where Λ1 and Λ2 can be expressed as (11)Λ1=RG1RC2λ+ρC2λ2IS+iLtanh2ρv21,(12)Λ2=1Lmλ2ISeρv2coshρv2ρLmλ2IS+iLρ2IS+iLeρv22IScoshρv2ρ2IS+iLeρv2sinhρv22IScosh2ρv2.

Then the eigenvalues at equilibrium points O1,2,3 can be calculated as (13)O1:ζ1=61,124,000,ζ2=20,000,ζ3,4=4000±26,000i,O2,3:ζ1,2=6106±26,054i,ζ3,4=10,038±15,524i.

It can be seen that equilibrium point O1 has a positive real root and equilibrium point O2,3 have two complex conjugate roots with the positive real part. As a result, the system is unstable and has chaotic attractors.

3.2. Bifurcation with the Changing of Order <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M92"><mml:mi>λ</mml:mi></mml:math></inline-formula>

The fractional-order models can increase the flexibility and degrees of freedom by means of the fractional parameters. In this section, the change of the fractional order λ is applied in the fractional-order memristive Chua’s circuit, and the influence of the order on the performance of the fractional-order memristive chaotic circuit is discussed. Taking the same circuit parameters as mentioned above and the initial conditions are set as (0 V, 0.01 V, 0 A, 0 A), the bifurcation diagram of the fractional-order chaotic circuit described in (7) can be obtained, as depicted in Figure 6. It shows that the dynamical behaviors of the system can be divided into three states: periodic state, bifurcation, and chaos. The system has two blocks of chaos around λ=0.93 and λ=0.99 when the order λ varies from 0.9 to 1. With the increase of the order, the system enters into the chaotic state from the periodic state abruptly when λ>0.92. However, the chaotic behavior does not exist for a long time; the system runs back to the periodic state when order λ is just over 0.94. When the order is higher than 0.97, the system enters into chaos once more. Furthermore, the system displays complicated nonlinear dynamic behavior such as bifurcation and chaos when order λ is greater than 0.96. In order to further illustrate the influence of the order on the dynamic behaviors of the system, a comparative analysis is carried out. The phase diagrams of the fractional-order memristive Chua’s circuit when order equals to 0.955 and 0.965 are indicated in Figure 7. It can be seen that the circuit is a single cycle when λ=0.955, but when λ=0.965, the circuit has two or more cycles, bifurcation occurs, which corresponds well with the results in the bifurcation diagram as shown in Figure 6. In addition, it is worth noting that the orbits of state variable v1t with order λ do not change smoothly with order λ, which is different from most of other bifurcation diagrams. From the above analysis, we know that different orders will cause different dynamical behaviors of the memristive chaotic circuit, such as bifurcation and chaos, which indicates that the order plays a significant role in the performance of the fractional-order memristive chaotic circuit.

Bifurcation diagram of the fractional-order memristive chaotic circuit with the different order.

The phase diagrams of the fractional-order memristive Chua’s circuit with the different order.

3.3. Bifurcation with the Changing of Inductor <inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M105"><mml:msup><mml:mrow><mml:mi>L</mml:mi></mml:mrow><mml:mrow><mml:mi>λ</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>

The bifurcation diagram of v2t with inductor Lλ is depicted in Figure 8.

Bifurcation diagrams of the fractional-order memristive chaotic circuit and integer-order memristive chaotic circuit with different inductance Lλ.

The orbits painted in red indicate the integer-order system with order λ=1, and the orbits painted in blue express the fractional-order system whose order λ=0.92. For the integer-order system, as the increase of inductor value, the orbits of the memristive Chua’s circuit start from the periodic state and then enter into chaos abruptly. After increasing for a period of time, the system returns back to periodic behavior. The phenomenon is not the same as that in the fractional-order system. The inductance value required for the fractional-order system to enter into chaotic behavior is larger than that of the integer-order system. Moreover, the fractional-order system is still in chaos when the inductance is greater than 33 mH, but the orbits are in a periodic state for the integer-order systems. In order to verify whether the theoretical analysis of the bifurcation diagram is correct, the phase diagrams of the integer-order system and fractional-order system with order λ=0.92 are given in Figure 9 when inductance Lλ=20 mH and Lλ=40 mH, respectively. It can be seen that the results of the phase diagrams well demonstrate the correctness of the bifurcation diagram. The integer-order system is chaotic but the fractional-order system is periodic when inductance value Lλ=20 mH. Meanwhile, the integer-order system has a periodic behavior but the fractional-order system has a chaotic behavior when inductance value Lλ=40 mH.

Phase diagrams of the integer-order system and fractional-order system: (a) Lλ=20 mH and (b) Lλ=40 mH.

4. Circuit Realization of the Fractional-Order Memristive Chaotic Circuit

In this section, equivalent circuits are used to simulate the fractional-order memristive chaotic circuit. Since the inductor and capacitor can be extended to fractional order, the Oustaloup filter algorithm is used to obtain the approximate transfer function of fractional Laplace transform operator sλ , then simplifying the transfer functions into zero-pole forms and using chain circuits to realize the equivalent circuit of fractional-order inductor and capacitor .

4.1. Equivalent Circuit of Fractional-Order Generalized Memristor

The fractional Laplace transform operator of inductor in fractional-order generalized memristor can be physically realized using the RL chain circuits in Figure 10. The transfer function of the fractional-order inductor can be approximately expressed as (14)1Lsλ=i=1n1/Lis+Ri/Li+1R0,where λ is the order of the fractional-order inductor. Selecting L=280 mH, λ=0.95,0.98 and n=7, resistance and inductance can be calculated by the method of undetermined coefficients. The parameters of the equivalent chain circuit of fractional-order inductor are listed in Tables 1 and 2. The magnitude and phase of the ideal 1/Lsλ and 7th order Oustaloup approximation and fractance circuit approximation with L=280 mH and order λ=0.95,0.98 are presented in Figure 11. From the figure, we observe that the approximation is relatively good from 10−1 Hz to 105 Hz.

Equivalent chain circuit of fractional-order inductor.

The resistance parameters of the equivalent chain circuit of fractional-order inductor.

λ R0 (Ω) R1 (Ω) R2 (Ω) R3 (Ω) R4 (Ω) R5 (Ω) R6 (Ω) R7 (Ω)
0.95 1.25e6 7.74e4 964.42 12.14 1.53e1 1.92e3 2.52e6 7.76e1
0.98 2.03e6 2.55e5 2.77e3 30.39 3.33e1 3.65e3 1.78e6 4.18e1

The inductance parameters of the equivalent chain circuit of fractional-order inductor.

λ L1 (H) L2 (H) L3 (H) L4 (H) L5 (H) L6 (H) L7 (H)
0.95 0.6896 0.8595 1.0820 1.3621 1.7147 0.4440 1.3688e5
0.98 2.4385 2.6473 2.9025 3.1825 3.4893 0.3367 7.9025e4

Bode diagram of Oustaloup approximation and fractance circuit approximation.

4.2. Equivalent Circuit of Fractional-Order Memristive Chua’s Circuit

After realizing the equivalent circuit of fractional-order inductor, we consider using the same method to achieve the equivalent circuit of the fractional-order capacitor. The approximate equivalence of the fractional-order capacitor can be actualized utilizing RC ladder topologies, as shown in Figure 12. The transfer function of the RC ladder network is expressed as (15)1Csλ=i=1n1/Cis+1/RiCi+R0.

Equivalent ladder circuit of the fractional-order capacitor.

In consideration of capacitors C1λ=100 nF and C2λ=10 nF in the fractional-order memristive Chua’s circuit, choosing λ=0.98 and n=7, resistance and capacitance can be calculated as well, as listed in Tables 3 and 4. The bode diagram of the ideal 1/Csλ and 7th order Oustaloup approximation and fractance circuit approximation with Cλ=10 nF and 100 nF and order λ=0.98 is indicated in Figure 13. It can be seen that the approximation is pretty good from 10−2 Hz to 105 Hz. The chain equivalent parameters of fractional-order inductor Lλ in the fractional-order memristive Chua’s circuit can be calculated as well, as listed in Tables 5 and 6, and the bode diagram shows a good approximation from 10−2 Hz to 105 Hz, as shown in Figure 14. Thus, the equivalent circuit of the fractional-order memristive Chua’s circuit has been completed as shown in Figure 15, which has four equivalent circuit units, and then the circuit simulation can be carried out on the fractional-order memristive equivalent circuit.

The resistance parameters of the equivalent ladder circuit of the fractional-order capacitor.

Cλ (nF) R0 (Ω) R1 (Ω) R2 (Ω) R3 (Ω) R4 (Ω) R5 (Ω) R6 (Ω) R7 (Ω)
10 13.80 109.66 1.01e4 9.21e5 8.40e7 7.66e9 1.57e13 6.70e7
100 1.38 10.97 1.01e3 9.21e4 8.40e6 7.66e8 1.57e12 6.70e6

The capacitance parameters of the equivalent ladder circuit of the fractional-order capacitor.

Cλ (nF) C1 (μF) C2 (μF) C3 (μF) C4 (μF) C5 (μF) C6 (μF) C7 (μF)
10 8.71e2 9.45e2 1.04e1 1.14e1 1.25e1 1.20e2 2.82e3
100 8.71e1 9.45e1 1.04 1.14 1.25 1.20e1 2.82e4

Bode diagram of Oustaloup approximation and fractance circuit approximation.

The resistance parameters of the equivalent chain circuit of fractional-order inductor.

R0 (Ω) R1 (Ω) R2 (Ω) R3 (Ω) R4 (Ω) R5 (Ω) R6 (Ω) R7 (Ω)
8.69e4 1.09e4 118.80 1.30 1.43e2 1.57e4 7.63e8 1.79e2

The inductance parameters of the equivalent chain circuit of fractional-order inductor.

L1 (H) L2 (H) L3 (H) L4 (H) L5 (H) L6 (H) L7 (H)
0.1045 0.1135 0.1244 0.1364 0.1495 0.0144 3.3868e3

Bode diagram of Oustaloup approximation and fractance circuit approximation.

The equivalent realization of the fractional-order memristor-based chaotic circuit.

4.3. Circuit Simulations of the Fractional-Order Memristive Equivalent Circuit

In this section, we consider utilizing PSpice to simulate the fractional-order memristive equivalent circuit. Considering order λ=0.98, the parameters of resistance Rm and diodes in fractional-order generalized memristor are the same as mentioned above, and fractional-order inductor Lm=280 mH is realized by a chain unit, whose equivalent parameters are listed in Tables 1 and 2. In addition, resistance R and negative conductance G in the fractional-order memristive chaotic circuit are constant, while fractional-order capacitance C1λ=100 nF and C2λ=10 nF and fractional-order inductance Lλ=12 mH are implemented by corresponding equivalent circuits, whose parameters are given in Tables 36. Thus, an equivalent circuit of the fractional-order memristor-based fractional-order Chua’s circuit is achieved. Figure 16 shows the phase diagrams of the realized fractional-order memristive Chua’s circuit by PSpice, and the results well verify the correctness of numerical simulations and the practicability of the fractional-order equivalent circuit.

Phase diagrams of the realized fractional-order memristive Chua’s circuit. Phase diagrams in the (a) v2tv1t plane, (b) v2ti3t plane, (c) v2tiLt plane, and (d) v2timt plane.

5. Conclusions

In this paper, a fractional-order memristor-based chaotic circuit is presented. The inductor in generalized memristor and the capacitors and inductors in the memristive chaotic circuit are all fractional order. Firstly, the mathematical model of fractional-order inductor-based generalized memristor is established, and then the characteristic equations of the fractional-order memristive chaotic circuit are derived. After that, theoretical analysis and numerical simulations are carried out such as the analysis of equilibrium point and stability and the influence of the order on the stability of the system. The results show that the fractional-order circuit exhibits different dynamic behaviors such as bifurcation and chaos with the change of the order, which indicates the importance of order effects on the dynamic behaviors of the system. In order to verify the aforementioned analysis, the equivalent circuit of the fractional-order memristive chaotic circuit is presented. The nonlinear elements inductor and capacitor are approximately equivalent through a unit circuit. Using the method of undetermined coefficients to solve circuit parameters, a fractional-order memristor-based equivalent circuit is constructed in PSpice. The results of circuit simulations well confirm the investigation on the fractional-order memristive chaotic circuit both in theoretical and numerical.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant No. 51507134), Scientific Research Program Funded by Shaanxi Provincial Department of Water Resources (Grant No. 2017slkj-15), Natural Science Foundation of Shaanxi Province (Grant No. 2018JM5068), Key Project of Natural Science Basic Research Plan in Shaanxi Province of China (Grant No. 2018ZDXM-GY-169), and Xi’an Science and Technology Bureau innovation project (Grant No. 201805037 (21)).

Chua L. O. Memristor—the missing circuit element IEEE Transactions on Circuit Theory 1971 18 5 507 519 10.1109/TCT.1971.10833372-s2.0-0015127532 Di Ventra M. Pershin Y. V. Chua L. O. Circuit elements with memory: memristors, memcapacitors, and meminductors Proceedings of the IEEE 2009 97 10 1717 1724 10.1109/JPROC.2009.20210772-s2.0-70349684961 Mohanty S. P. Memristor: from basics to deployment IEEE Potentials 2013 32 3 34 39 10.1109/MPOT.2012.22162982-s2.0-84877740461 Dongale T. D. Patil K. P. Mullani S. B. More K. V. Delekar S. D. Patil P. S. Gaikwad P. K. Kamat R. K. Investigation of process parameter variation in the memristor based resistive random access memory (RRAM): effect of device size variations Materials Science in Semiconductor Processing 2015 35 174 180 10.1016/j.mssp.2015.03.0152-s2.0-84925622473 Thangamani V. Memristor-based resistive random access memory: hybrid architecture for low power compact memory design Control Theory & Informatics 2014 4 Thomas A. Memristor-based neural networks Journal of Physics D: Applied Physics 2013 46 9, article 093001 10.1088/0022-3727/46/9/0930012-s2.0-84874057890 Duan S. Dong Z. Hu X. Wang L. Li H. Small-world Hopfield neural networks with weight salience priority and memristor synapses for digit recognition Neural Computing and Applications 2016 27 4 837 844 10.1007/s00521-015-1899-72-s2.0-84926155964 Rakkiyappan R. Chandrasekar A. Laksmanan S. Park J. H. State estimation of memristor-based recurrent neural networks with time-varying delays based on passivity theory Complexity 2014 19 4 43 10.1002/cplx.214822-s2.0-84897614854 Mouttet B. Proposal for memristors in signal processing International ICST Conference 2008 Berlin, Heidelberg Springer 11 13 El-Sayed A. M. A. Elsaid A. Nour H. M. Elsonbaty A. Dynamical behavior, chaos control and synchronization of a memristor-based ADVP circuit Communications in Nonlinear Science and Numerical Simulation 2013 18 1 148 170 10.1016/j.cnsns.2012.06.0112-s2.0-84865637481 Sheridan P. Ma W. Lu W. Pattern recognition with memristor networks 2014 IEEE International Symposium on Circuits and Systems June 2014 Melbourne, VIC, Australia 1078 1081 Chua L. Everything you wish to know about memristors but are afraid to ask Radioengineering 2015 24 2 319 368 10.13164/re.2015.03192-s2.0-84932605254 Corinto F. Ascoli A. Memristive diode bridge with LCR filter Electronics Letters 2012 48 14 824 825 10.1049/el.2012.14802-s2.0-84863656814 Wu H. Bao B. Liu Z. Xu Q. Jiang P. Chaotic and periodic bursting phenomena in a memristive Wien-bridge oscillator Nonlinear Dynamics 2016 83 1-2 893 903 10.1007/s11071-015-2375-82-s2.0-84940921367 Bao B. Yu J. Hu F. Liu Z. Generalized memristor consisting of diode bridge with first order parallel RC filter International Journal of Bifurcation and Chaos 2014 24 11, article 1450143 10.1142/S02181274145014302-s2.0-84920900228 Chen M. Li M. Yu Q. Bao B. Xu Q. Wang J. Dynamics of self-excited attractors and hidden attractors in generalized memristor-based Chua’s circuit Nonlinear Dynamics 2015 81 1-2 215 226 10.1007/s11071-015-1983-72-s2.0-84930758110 Wu H. G. Bao B. C. Xu Q. First order generalized memristor emulator based on diode bridge and series RL filter Sexualities 2015 16 16 651 664 Bao B. C. Wu P. Bao H. Chen M. Xu Q. Chaotic bursting in memristive diode bridge-coupled Sallen-Key lowpass filter Electronics Letters 2017 53 16 1104 1105 10.1049/el.2017.16472-s2.0-85027450623 Xu Q. Zhang Q. Wang N. Wu H. Bao B. An improved memristive diode bridge-based band pass filter chaotic circuit Mathematical Problems in Engineering 2017 2017 11 2461964 10.1155/2017/24619642-s2.0-85042169332 Itoh M. Chua L. O. Memristor oscillators International Journal of Bifurcation and Chaos 2008 18 11 3183 3206 10.1142/S02181274080223542-s2.0-58449129665 Muthuswamy B. Kokate P. P. Memristor-based chaotic circuits IETE Technical Review 2009 26 6 417 429 10.4103/0256-4602.578272-s2.0-73249115972 Muthuswamy B. Implementing memristor based chaotic circuits International Journal of Bifurcation and Chaos 2010 20 5 1335 1350 10.1142/S02181274100265142-s2.0-77954578392 Bao B. Ma Z. Xu J. Liu Z. Xu Q. A simple memristor chaotic circuit with complex dynamics International Journal of Bifurcation and Chaos 2011 21 9 2629 2645 10.1142/S02181274110299992-s2.0-81055149518 Wu R. Wang C. A new simple chaotic circuit based on memristor International Journal of Bifurcation and Chaos 2016 26 9, article 1650145 10.1142/S02181274165014552-s2.0-84984994773 Teng L. Iu H. H. C. Wang X. Wang X. Chaotic behavior in fractional-order memristor-based simplest chaotic circuit using fourth degree polynomial Nonlinear Dynamics 2014 77 1-2 231 241 10.1007/s11071-014-1286-42-s2.0-84902242203 Rajagopal K. Guessas L. Karthikeyan A. Srinivasan A. Adam G. Fractional order memristor no equilibrium chaotic system with its adaptive sliding mode synchronization and genetically optimized fractional order PID synchronization Complexity 2017 2017 19 1892618 10.1155/2017/18926182-s2.0-85018636379 Zou L. Peng Y. Feng Y. Tu Z. Stabilization and synchronization of memristive chaotic circuits by impulsive control Complexity 2017 2017 10 5186714 10.1155/2017/51867142-s2.0-85042232288 Petras I. Chen Y. Q. Coopmans C. Fractional-order memristive systems 2009 IEEE Conference on Emerging Technologies & Factory Automation September 2009 Mallorca, Spain 1 8 10.1109/etfa.2009.53471422-s2.0-77949897523 Petras I. Chen Y. Q. Fractional-order circuit elements with memory Proceedings of the 13th International Carpathian Control Conference (ICCC) May 2012 High Tatras, Slovakia 552 558 10.1109/carpathiancc.2012.62287062-s2.0-84864231388 Petras I. Fractional-order memristor-based Chua's circuit IEEE Transactions on Circuits and Systems II: Express Briefs 2010 57 12 975 979 10.1109/TCSII.2010.20831502-s2.0-78650415784 Yang N. Xu C. Wu C. Jia R. Liu C. Modeling and analysis of a fractional-order generalized memristor-based chaotic system and circuit implementation International Journal of Bifurcation and Chaos 2017 27 13, article 1750199 10.1142/s02181274175019912-s2.0-85040076896 Oldham K. B. Spanier J. The fractional calculus Mathematical Gazette 1974 56 247 396 400 Podlubny I. Fractional differential equations vol. 198 of Mathematics in Science and Engineering 1999 San Diego, CA, USA Academic Press Hartley T. T. Lorenzo C. F. The error incurred in using the caputo-derivative Laplace-transform Volume 4: 7th International Conference on Multibody Systems, Nonlinear Dynamics, and Control, Parts A, B and C August-September 2009 San Diego, California, USA 271 278 10.1115/DETC2009-876482-s2.0-77953704972 Wang F. Z. Shi L. Wu H. Helian N. Chua L. O. Fractional memristor Applied Physics Letters 2017 111 24, article 243502 10.1063/1.50009192-s2.0-85038579069 Podlubny I. Petráš I. Vinagre B. M. O'Leary P. Dorčák Ľ. Analogue realizations of fractional-order controllers Nonlinear Dynamics 2002 29 1/4 281 296 10.1023/A:10165566043202-s2.0-0036650890 Heymans N. Podlubny I. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives Rheologica Acta 2006 45 5 765 771 10.1007/s00397-005-0043-52-s2.0-33745869026 Abdelouahab M.-S. Lozi R. Chua L. Memfractance: a mathematical paradigm for circuit elements with memory International Journal of Bifurcation and Chaos 2014 24 9, article 1430023 10.1142/S02181274143002372-s2.0-84928485676 Adhikari S. P. Sah M. P. Kim H. Chua L. O. Three fingerprints of memristor IEEE Transactions on Circuits and Systems I: Regular Papers 2013 60 11 3008 3021 10.1109/TCSI.2013.22561712-s2.0-84887136051 Oustaloup A. Levron F. Mathieu B. Nanot F. M. Frequency-band complex noninteger differentiator: characterization and synthesis IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications 2000 47 1 25 39 10.1109/81.8173852-s2.0-0033888365 Freeborn T. Design and Implementation of Fractional Step Filters 2010 Diss. University Of Calgary