# Effects of Parasitic Fractional Elements to the Dynamics of Memristor

Authors: Rawid Banchuin
Journal: Journal of Electrical and Computer Engineering
 (2019)

---

## Abstract

In this research, we study the effects of the parasitic fractional elements to the dynamic of the memristor where both flux- and charge-controlled memristors have been considered. For doing so, the fractional differential equation-based approach has been used for modeling the memristor and the memristor-based circuits under the effects of the parasitic fractional elements where the resulting equations have been solved both analytically and numerically. From the obtained solutions and simulations, the effects of the parasitic fractional elements to the dynamic of the memristor have been studied. We have found that the parasitic fractional elements cause the charge and flux decay of the memristor similarly to the conventional parasitic elements. Moreover, the impasse points of the phase portraits between flux and charge of the memristor-based circuits can also be broken by the parasitic fractional elements. The effects of the order and the nonlinearity of the parasitic fractional elements have also been studied.


---

## Body

## 1. Introduction

Memristor is the 4th electrical circuit element that has been theoretically found by Leon Chua since 1971 [1]. The memristor can be simply thought of as the resistor with memory. For decades after Chua proposed his original work, the memristor has been practically realized by a research group in Hewlett Packard (HP) labs [2] in 2008. The memristor has been applied in various applications, e.g., parametric oscillator, DRAM, signal processing, neural networks, programmable logic, cross bar switch, passive switch, and control systems. [3–6]. According to the importance of the memristor, the studies in its circuit theoretical aspects have been proposed [7–10]. In [9], the effects of the parasitic element to the memristor have been studied where it has been found that such parasitic element exhibits significant effects to the memristor and memristor-based circuits. However, it was assumed that the orders of the parasitic elements are strictly integers despite that they can be fractional in practice as those parasitic element can also be the fractional elements or fractances [11–13].Therefore, we study the dynamic of the memristor and the memristor-based circuits under the effects of parasitic elements, which are the fractional elements, in this work where comparison with [9] has also been made. For simplicity, we appeal the parasitic elements, which are the fractional elements, as the parasitic fractional elements. Here, both flux- and charge-controlled memristors have been considered. Since the order of the parasitic fractional element can be fractional, the fractional differential equation (FDE), which its order can be fractional, must be used for modeling the memristor and the memristor-based circuits under the effects of such parasitic fractional element instead of the ordinary differential equation (ODE) which its order is strictly integer. The FDE is the extension of the ODE in the domain of fractional calculus which plays a fundamental role in various engineering fields, e.g., bioengineering, control theory, electronics, robotics, and signal processing [14–18]. The analytical solutions of the formulated FDEs, which are nonlinear, have been determined by using the Adomian decomposition method [19]. On the other hand, the numerical solutions have been obtained by using the integer derivative-based approximation of the fractional derivative with the moments of a function [20]. From the obtained solutions and simulations with MATHEMATICA, the effects of the parasitic fractional element to the dynamic of the memristor and the memristor-based circuits have been studied. We have found that the parasitic fractional element also causes the charge and flux decay of the memristor and breaks the impasse points of the phase portraits between flux and charge of the memristor-based circuits. Moreover, the effects of the order and the nonlinearity of the parasitic fractional element have also been studied.In the following section, some background theory on the memristor and the fractional element will be briefly given. The flux and charge decaying of the memristor due to the parasitic fractional element will be illustrated in Section3 followed by the impasse points breaking of the phase portraits between flux and charge of the memristor-based circuits in Section 4 where the effects of the order of the parasitic fractional element will be studied in these sections. The nonlinear effect of the parasitic fractional element will be explored in Section 5, and the conclusion will be finally drawn in Section 6.


## 2. Background Theory

### 2.1. Memristor

Memristor is a nonlinear electrical circuit element. This circuit element relates the flux (φ) and charge (q) through one of the following constitutive relations [9]:(1a)φ=fq,(1b)q=gφ.Noted that bothf(·) and g⋅ must be differentiable and scalar-valued. Moreover, the memristor with its constitutive relation given by (1a) and (1b) can be referred to as the q-controlled memristor and the φ-controlled memristor, respectively. According to [21], the ideal memristor is nonvolatile. Therefore, it can retain the values of φ and q at the moment that it is switched off, i.e., φ and q are always conserved without any change as time goes by. This property of the ideal memristor corresponds to the fundamental circuit law on φ and q, namely, the principle of the conservation of charge and flux [22].From (1a) and (1b), the following mathematical expression of small-signal memristance, M(q), and memductance, W(φ), of the memristor can be obtained [9]:(2)Mq=dfqdq,(3)Wφ=dgφdφ.


### 2.2. Fractional Element

Fractional element is a circuit element with the impedance function of fractional order,α, where 0 < α < 1. According to [23], the fractional element which the phase of its voltage always leads that of its current is referred to as the fractional inductor. On the other hand, such element which the phase of its voltage always lags that of its current is referred to as fractional capacitor. For the fractional inductor, its s-domain impedance function can be given by(4)ZLα=Lαsα,where Lα denotes the inductivity of the fractional inductor [24]. On the other hand, we have the following impedance function for the fractional capacitor:(5)ZCα=Cαsα−1,where Cα stands for the pseudocapacitance of the fractional capacitor [25].By using (4) and (5), the following flux-charge relationships can be, respectively, obtained for the fractional inductor and fractional capacitor:(6)φLα=LαDαqLα,(7)qCα=CαDαφCα,where Dα[ ] denotes the fractional derivative operator. Moreover, φCα(φLα) and qCα(qLα) stand for the flux and charge of the fractional capacitor (fractional inductor).


## 2.1. Memristor

Memristor is a nonlinear electrical circuit element. This circuit element relates the flux (φ) and charge (q) through one of the following constitutive relations [9]:(1a)φ=fq,(1b)q=gφ.Noted that bothf(·) and g⋅ must be differentiable and scalar-valued. Moreover, the memristor with its constitutive relation given by (1a) and (1b) can be referred to as the q-controlled memristor and the φ-controlled memristor, respectively. According to [21], the ideal memristor is nonvolatile. Therefore, it can retain the values of φ and q at the moment that it is switched off, i.e., φ and q are always conserved without any change as time goes by. This property of the ideal memristor corresponds to the fundamental circuit law on φ and q, namely, the principle of the conservation of charge and flux [22].From (1a) and (1b), the following mathematical expression of small-signal memristance, M(q), and memductance, W(φ), of the memristor can be obtained [9]:(2)Mq=dfqdq,(3)Wφ=dgφdφ.


## 2.2. Fractional Element

Fractional element is a circuit element with the impedance function of fractional order,α, where 0 < α < 1. According to [23], the fractional element which the phase of its voltage always leads that of its current is referred to as the fractional inductor. On the other hand, such element which the phase of its voltage always lags that of its current is referred to as fractional capacitor. For the fractional inductor, its s-domain impedance function can be given by(4)ZLα=Lαsα,where Lα denotes the inductivity of the fractional inductor [24]. On the other hand, we have the following impedance function for the fractional capacitor:(5)ZCα=Cαsα−1,where Cα stands for the pseudocapacitance of the fractional capacitor [25].By using (4) and (5), the following flux-charge relationships can be, respectively, obtained for the fractional inductor and fractional capacitor:(6)φLα=LαDαqLα,(7)qCα=CαDαφCα,where Dα[ ] denotes the fractional derivative operator. Moreover, φCα(φLα) and qCα(qLα) stand for the flux and charge of the fractional capacitor (fractional inductor).


## 3. Flux and Charge Decay by the Parasitic Fractional Elements

Despite the conservation ofφ and q in the ideal phenomena, the memristor with parasitic elements exhibits the decaying of both φ and q after being switched off as proposed in [9]. In this work, we will show that such decaying of both φ and q can also be occurred by the parasitic fractional elements.Firstly, the parasitic fractional capacitor will be considered. The memristor with parasitic fractional capacitor can be modeled as an ideal memristor connected in parallel with a fractional capacitor as depicted in Figure1; thus, we have φCα=φ. If we assume that the memristor is φ-controlled, the following relationship can be obtained by using the conventional circuit analysis after the memristor is switched off:(8)qCα=−gφ.Figure 1
Equivalent circuit model of the memristor affected by the parasitic fractional capacitor.By using (1b), (7), and (8), the following FDE can be obtained:(9)Dαφ=−qCα.Therefore, we have found that(10)q∝−Dαφ.On the other hand, the relationships of the memristor affected by parasitic conventional capacitor,C, have been proposed in [9] as follows:(11)dϕdt=−qC,(12)q∝−dϕdt.Since the time domain fractional derivative is also referred to as the rate of change similarly to ordinary derivative, it can be found from (10) that φ of the memristor affected by the parasitic fractional capacitor is decreased with respect to time when q > 0 and vice versa similarly to that of the parasitic conventional capacitor-affected device as proposed in [9] and can be seen in this work from (11). Here, we assume that q=gφ lies in the first and third quadrants of φ-q plane as it is a strictly increasing function and g0=0 as also assumed in [9]. After applying the dynamic route theory [26] to (10), we have found that(13)limt⟶∞φ=0.Thus, it can be stated that the parasitic fractional capacitor causes the decaying ofφ after the memristor is switched off similarly to the parasitic conventional capacitor as proposed by [9]. Since q=gφ is a strictly increasing function, we have(14)limt⟶∞q=0,which states the decaying of q is occurred by the parasitic fractional capacitor as well as by the parasitic conventional capacitor [9]. It should be mentioned here that both (13) and (14) can also be obtained by applying the dynamic route theory to (11) and keeping in mind that q=gφ.Now, consider the parasitic fractional inductor. The memristor with parasitic fractional inductor can be modeled as an ideal memristor connected in series with a fractional inductor as depicted in Figure2; therefore, qLα=q. By assuming that the memristor is q-controlled and applying (1a) and (6), the following equation can be obtained after the memristor is switched off:(15)Dαq=−φLα.Figure 2
Equivalent circuit model of the memristor affected by the parasitic fractional inductor.From (15), it can be seen that(16)φ∝−Dαq,which shows that q is decreased with respect to time when φ > 0 and vice versa. If we assume that φ = f(q) lies in the first and third quadrants of q-φ plane and f(0) = 0, we have found that(17)limt⟶∞q=0.by using the dynamic route theory.Asφ = f(q) is a strictly increasing function, we also have(18)limt⟶∞φ=0.Therefore, it can be seen that the parasitic fractional inductor induces the decaying ofq which in turn causes the decaying of φ when the memristor is switched off similarly to the parasitic conventional inductor as proposed in [9] by using the intuition without the mathematical proof.As a practical illustration, let the constitutive relation of the memristor be specifically given byq=φ3. Therefore, the effect of the parasitic fractional capacitor can be analytically described by the following nonlinear FDE:(19)Dαφ=φ3Cα,which can be derived by using (9) and the above constitutive relation.If the parasitic conventional capacitor has been considered on the other hand similarly to the previous study [9], its effect can be described by the following ODE:(20)dφdt=φ3C.By using the Adomian decomposition, the following analytical solution of (17) can be obtained:(21)φ=φ0+φ1+φ2+φ3⋯,where φ0 denotes the initial value of φ (at t = 0). On the other hand, the higher order terms can be given in terms of φ0 as follows:(22)φ1=−φ03tαCαΓα+1,φ2=3φ05t2αCα2Γ2α+1,φ3=−φ07t3α3Γα+12+Γ2α+1Cα3αΓα+12Γ3α,⋮where Г(·) stands for the gamma function [26].For determining the numerical solution of (20), the integer derivative-based approximation of Dαφ by using the moments of a function has been used. As a result, the following system of ODE can be obtained:(23)dφdt=−φ3/Cα+∑p=2NΓp−1+αt1−p−αVp/Γ2−αΓα−1p−1!−1+∑p=2NΓp−1+α/Γαp−1!/Γ1−αt−αφ1+∑p=1NΓp−1+α/Γα−1p!/Γ2−αt1−α,dV2dt=−φ,dV3dt=−2tφ,⋮dVNdt=1−NtN−2φ,where the initial values of V2, V3, …, VN have been automatically assumed to be 0 due to the adopted approximation methodology [21].Here, we letN = 7 in (23) as this value of N gives a satisfied approximation of the fractional derivative [20] while being not too large; thus, the numerical solution of (17), i.e., φ, can be obtained by also maintaining a reasonable amount of computational effort. The obtained φ can be simulated against t for various values of α as depicted in Figure 3 where φ0 = 1 and Cα = 0.01 have been assumed. Similarly to [7–10], all related quantities have been given in dimensionless manner without loss of generality as this work also focuses on the circuit theoretical aspect of memristor. Since q=φ3, q can be simulated as depicted in Figure 4 where the similar assumption has been adopted. From these figures which φ and q of the memristor affected by the conventional parasitic element [9] which can be simulated here by using the numerical solution of (20) have also been included for comparison, the decaying of both φ and q due to the parasitic fractional element can be clearly observed. This confirms the aforementioned dynamic route theory based on generic results. Moreover, we have also found from Figures 3 and 4 that the decaying rates of φ and q, which are much faster than those of the φ and q decaying caused by the conventional parasitic element, are inversely proportional to α. This is because the discharging rate of the fractional capacitor and the flux decaying rate of the fractional inductor which serve as the parasitic elements in this scenario are inversely proportional to α. For a demonstration, consider a resistor connected in parallel with the fractional capacitor as shown in Figure 5. After the resistor being switched off, the following linear FDE can be obtained by applying the conventional circuit analysis where R stands for the resistance of the resistor:(24)DαqCα+qCαRCα=0.Figure 3
φ vs t of the memristor affected by the parasitic fractional capacitor and by the conventional parasitic capacitor [9].Figure 4
q vs t of the memristor affected by the parasitic fractional capacitor and by the conventional parasitic capacitor [9].Figure 5
The source-free resistor-fractional capacitor circuit.After the Laplace transformation and inverse Laplace transformation, we have(25)qCα=q0CαEα−tαRCα,where q0Cα denotes the initial value of qCα (at t = 0) and Eα stands for the Mittag–Leffler function [26] which can be defined in terms of arbitrary variable, x, as(26)Eαx=∑k=0∞xkΓαk+1.On the other hand, the charge of the conventional capacitor-resistor circuit, i.e.,qC, obeys(27)dqCdt+qCRC=0.Thus, it can be analytically given by(28)qC=q0Cexp−tRC,where q0C denotes the initial value of qC.With (25), qCα can be simulated for various values of α by assuming that q0Cα = 1, Cα = 0.01, and R = 1 as depicted in Figure 6 where the charge of the conventional capacitor-resistor circuit, i.e., qC, simulated by using (28) has also been included. It can be seen from Figure 6 that the decreasing rate of qCα, which is much faster than that of qC and referred to the discharging rate of the fractional capacitor, is inversely proportional to α. Noted that the aforesaid relationship between the flux decaying rate of the fractional inductor and α can also be demonstrated in a similar manner.Figure 6
q
C
α
/
q
C vs t of a source-free fractional capacitor-resistor circuit/a source-free conventional capacitor-resistor circuit.


## 4. Breaking the Impasse Points of the Phase Portraits between Flux and Charge

According to [9], the phase portraits between φ and q of the φ-controlled memristor-inductor circuit exhibit discontinuous jumps at the impasse points which can be broken by the effect of parasitic capacitance. By using intuition, it can be stated that such impasse points breaking of q-controlled memristor-capacitor circuit can be occurred by the parasitic inductor albeit this is not mentioned in [9]. In this work, we study the effect of the parasitic fractional element to the phase portraits between φ and q of the memristor-inductor/capacitor circuits. Firstly, the memristor-inductor circuit affected by the parasitic fractional capacitor will be considered. Since we assume that the memristor and inductor are connected in parallel, the parasitic fractional capacitor included circuit can be modeled as depicted in Figure 7. Therefore, we have(29)dQLdt=−φL,(30)Dαφ=QL−gφCα,where QL = −qL which stands for the charge flowing through the inductor. Noted also that QL depends on q due to Kirchhoff’s charge law (KCHL) [8] and we assume that the memristor is φ-controlled.Figure 7
The memristor-inductor circuit with the parasitic fractional capacitor.If we let(31)gφ=φ33−φ.Equations (29) and (30) become(32)dQLdt=−φL,(33)Dαφ=QL−φ3/3+φCα.After applying the Adomian decomposition by keeping in mind thatdQL/dt≈D0.999…QL and the fractional integration of order 0.999… can be approximated by the integer one, we obtained(34)φ=φ0+φ1+φ2⋯,(35)QL=QL0+QL1+QL2⋯,where QL0 denotes the initial value of QL (at t = 0). Noted also that(36)φ1=−tαQL0+φ0−φ03/3CαΓα+1,φ2=tα3Cα2φ02−1φ03−3φ0−3QL0tαΓ2α+1−3Cαφ0tΓα+2,⋮(37)
Q
L
1
=
−
φ
0
t
L
,
Q
L
2
=
t
α
+
1
φ
0
3
−
3
φ
0
−
3
Q
L
0
3
L
C
α
Γ
α
+
2
.
⋮For determining the numerical solution of (32) and (33), the integer derivative-based approximation of Dαφ by using the moments of a function has also been used. As a result, we have(38)dQLdt=−φL,dφdt=QL/Cα−φ3/3Cα+φ/Cα+∑p=2NΓp−1+αt1−p−αVp/Γ2−αΓα−1p−1!−1+∑p=2NΓp−1+α/Γαp−1!/Γ1−αt−α1+∑p=1NΓp−1+α/Γα−1p!/Γ2−αt1−α,dV2dt=−φ,dV3dt=−2tφ,⋮dVNdt=1−NtN−2φ.In [9] where the conventional parasitic capacitor has been considered on the other hand, the following conventional ODEs will be obtained instead [9]:(39)dQLdt=−φL,(40)dφdt=QL−φ3/3+φCα.By using (38) with N = 7, the numerical solution of a system of FDE composed of (32) and (33), i.e., φ and QL, can be obtained and simulated against t for various values of α as depicted in Figures 8 and 9 where φ0  = 0.1, QL0  = 0.1, L = 1 [9], and Cα = 0.01 have also been assumed. Moreover, φ and QL of the memristor-inductor circuit affected by the conventional parasitic capacitor have also been included in these figures. Noted that the initial values of V2, V3, …, VN have been automatically assumed to be 0 due to the adopted approximation of the fractional derivative. From these figures, it can be seen that both φ and QL of both the parasitic fractional capacitor and conventional parasitic capacitor-affected circuits exhibit the oscillations, i.e., the memristor-inductor circuit with parasitic capacitor/fractional capacitor become an oscillator/a fractional oscillator. Noted that the simulation of the parasitic conventional capacitor-affected circuit studied in the previous work has been performed by using (39) and (40). As proposed in [9], the memristor-inductor circuit under the effect of the conventional parasitic capacitor becomes the conventional Van der Pol oscillator. In this work, it has been found that the parasitic fractional capacitor-affected memristor-inductor circuit becomes a fractional Van der Pol oscillator [27–29]. This is not surprising as (32) and (33) resemble the state space representation of such oscillator. We have also found that the amplitude of QL and the oscillating frequencies of both φ and QL are affected by α. This is also reasonable as the significant effect of the order on the output’s magnitude and frequency is a typical characteristic of any fractional oscillator [30–32].Figure 8
φ vs t of the memristor-inductor circuit with the parasitic fractional capacitor/conventional parasitic capacitor.Figure 9
Q
L vs φ of the memristor-inductor circuit with the parasitic fractional capacitor/conventional parasitic capacitor [9].Now, consider the phase portraits betweenφ and QL which can be simulated for various values of α by also assuming that φ0  = 1 and Cα = 0.01 as depicted in Figure 10 where the phase portrait of the conventional parasitic capacitor-affected memristor-inductor circuit proposed in [9] has also been included. Equation (38) has been used for the simulation of the circuit with the parasitic fractional capacitor where (39) and (40) have been adopted in this work in order to simulate the conventional parasitic capacitor-affected circuit for comparison. Similarly to that of the circuit affected by the conventional parasitic capacitor, there is none of any impasse point on each phase portrait of the memristor-inductor circuit with the parasitic fractional capacitor as there is no jump behavior occurred. Therefore, it can be stated that such impasse points on the phase portraits can also be broken by the effect of the parasitic fractional capacitor as well as the conventional parasitic capacitor [9] despite its fractional order. Such impasse point breakings yield the continuous oscillations of φ and QL as shown in Figures 8 and 9. Moreover, it has been found that α significantly affects both size and shape of the phase portraits. This is because the memristor-inductor circuit with the parasitic fractional capacitor becomes a fractional Van der Pol oscillator which the significant effects of α on the phase portraits between outputs are typical [27, 28], as stated above. As expected, the simulated phase portrait with α approaching 1 becomes more similar to that of the memristor-inductor circuit with the conventional parasitic capacitor. This is not surprising as the charge-flux relationship of the fractional capacitor with α = 1 is similar to that of the conventional capacitor as can be seen from (7). From the depicted phase portraits, it can also be seen that the memory and the linearity of the memristor-inductor-parasitic fractional capacitor circuit are, respectively, inversely and directly proportional to α.Figure 10
Q
L vs φ of the memristor-inductor circuit with the parasitic fractional capacitor/conventional parasitic capacitor [9].At this point, we will analyze the effects of the parasitic fractional inductor to the memristor-capacitor circuit. Since we assume that the memristor and capacitor are connected in parallel, the parasitic fractional inductor included circuit can be modeled as depicted in Figure11. As a result, the system of FDE given by (41) and (42) can be obtained if we let the memristor be of a q-controlled type with fq=q3/3−q.(41)dΦCdt=−qC,(42)Dαq=ΦC−q3/3+qLα.Figure 11
The memristor-capacitor circuit with parasitic fractional inductor.Noted thatΦC=φC which denotes the flux of the capacitor, and ΦC depends on φ due to Kirchhoff’s flux law (KFLL) [8]. By solving (41) and (42) analytically with the Adomian decomposition, we have(43)q=q0+q1+q2⋯,(44)ΦC=ΦC0+ΦC1+ΦC2⋯,where ΦC0 denotes the initial value of ΦC (at t = 0) and(45)q1=−tαΦC0+q0−q03/3LαΓα+1,q2=tα3Lα2q02−1q03−3q0−3ΦC0tαΓ2α+1−3Lαq0tΓα+2,⋮(46)ΦC1=−q0tC,ΦC2=tα+1q03−3q0−3ΦC03CLαΓα+2.⋮For determiningq and ΦC numerically, the integer derivative-based approximation of Dαq by using the moments of a function with N = 7 has been used. Besides other initial values which have been assumed to be 0 due to the adopted approximation of Dαq, we let q0 = 0.1, ΦC0  = 0.1, C = 1, and Lα = 0.01. The resulting q and ΦC and their phase portraits have been simulated for various values of α as depicted in Figures 12–14 where those of the conventional parasitic inductor-affected circuit have also been included. It should be mentioned here that the conventional parasitic inductor-related results can be simulated by using the following ODEs: (47)dΦCdt=−qC,(48)dqdt=ΦC−q3/3+qL.Figure 12
:q vs t of the memristor-capacitor circuit with the parasitic fractional inductor/conventional parasitic inductor.Figure 13
Φ
C vs t the memristor-capacitor circuit with the parasitic fractional inductor/conventional parasitic inductor.Figure 14
ΦC vs q the memristor-capacitor circuit with the parasitic fractional inductor/conventional parasitic inductor.From Figures12 and 13, α-controlled continuous oscillations of both q and ΦC can be observed for the circuit with the parasitic fractional inductor. This is because the entire circuit acts as a fractional oscillator. From the phase portraits between q and ΦC depicted in Figure 14, it can be observed that the impasse points on the phase portraits can also be broken by the parasitic fractional inductor similarly to the conventional parasitic inductor. It can also be seen that α employs significant effects to both size and shape of the phase portrait of the parasitic fractional inductor-affected circuit as such circuit becomes a fractional Van der Pol oscillator. Moreover, the phase portrait of the parasitic fractional inductor-affected circuit with α approaching 1 becomes more similar to that of the circuit with conventional parasitic inductor which is actually the memristor-capacitor-inductor-based Van der Pol oscillator [8]. This is because the fractional inductor with α = 1 employs a similar charge-flux relationship to that of the conventional inductor as can be seen from (6).


## 5. Nonlinear Effect of the Parasitic Fractional Element

The conventional parasitic element can be nonlinear [9] and so does the parasitic fractional element. Since the ODEs of the memristor-inductor circuit affected by the nonlinear conventional parasitic capacitor can be given by [9](49)dQLdt=−φL,(50)dφdt=QL−gφCφ,where Cφ which is a function of φ is the capacitance of the nonlinear conventional parasitic capacitor, the system of FDE of the memristor-inductor circuit affected by nonlinear parasitic fractional capacitor can be given as follows:(51)dQLdt=−φL,(52)Dαφ=QL−gφCαφ,where Cαφ which is a function of φ stands for the pseudocapacitance of the nonlinear parasitic fractional capacitor. Noted also that φCα=φ due to the circuit topology.By assuming thatCαφ=Cα1+φ2 where Cα means Cα0 in this scenario and using the Adomian decomposition, the analytical solution of (51) and (52) can also be given by (34) and (35) but with the following terms instead of those given by (36) and (37).(53)φ1=−tαQL0+φ0−φ03/3Cα1+φ02Γα+1,φ2=tα9Cα21+φ023φ04+6φ02+6φ0QL0−3QL0tαφ03−3φ0−3QL0−1Γ2α+1−9Cαφ0t1+φ022LΓα+2,⋮(54)
Q
L
1
=
−
φ
0
t
L
,
Q
L
2
=
t
α
+
1
φ
0
3
−
3
φ
0
−
3
Q
L
0
3
L
C
α
1
+
φ
0
2
Γ
α
+
2
.
⋮For determining the numerical solution of (51) and (52), the integer derivative-based approximation of Dαφ by using the moments of a function with N = 7 has been adopted and φ0 = 0.1, QL0 = 0.1, L = 1, and Cα = 0.01 have also been assumed. The other initial values have been assumed to be 0 by default according to the adopted approximation of Dαφ. The obtained solutions, i.e., φ and QL, and their phase portraits have been simulated for various values of α as depicted in Figures 15–17 where the nonlinear conventional parasitic capacitor-related results including the phase portraits [9] simulated here by using (49) and (50) with Cφ=C01+φ2 have also been included. From Figures 15 and 16, it can be seen that the nonlinear parasitic fractional capacitor also yields the continuous oscillations of φ and QL similarly to the nonlinear conventional parasitic capacitor as reported by [9]. Moreover, the properties of the oscillation caused by the nonlinear parasitic fractional capacitor are α dependent similarly to those caused by its linear counterpart.Figure 15
φ vs t of the memristor-inductor circuit with the nonlinear parasitic fractional capacitor/nonlinear conventional parasitic capacitor.Figure 16
Q
L vs t of the memristor-inductor circuit with the nonlinear parasitic fractional capacitor/nonlinear conventional parasitic capacitor.Figure 17
Q
L vs φ of the memristor-inductor circuit with the nonlinear parasitic fractional/nonlinear conventional parasitic capacitor [9].From Figure17, it can be observed that the nonlinear parasitic fractional capacitor also causes the impasse points breaking of the phase portraits between φ and QL similarly to the nonlinear conventional parasitic capacitor as also reported by [9]. Noted that both sizes and shapes of the phase portraits of the circuit with the nonlinear parasitic fractional capacitor are also significantly affected by α despite its nonlinearity. This is not surprising because the dependency on α of the phase portraits between outputs can be typically found in any nonlinear fractional system [33, 34] including our memristor-inductor-nonlinear parasitic fractional capacitor circuit. Moreover, it can be seen that the phase portraits between φ and QL depicted in Figure 17 have wider loop areas which refers to less linearity and more memory, than those of the previous case depicted in Figure 6. Such wider loop areas can be clearly observed when α approaches 1. On the other hand, the phase portraits depicted in Figure 12 also have more roundish shapes compared to those depicted in Figure 10 when α approaches 1. These differences in the sizes and shapes of the phase portraits are caused by the nonlinearities which have been taken into account in this scenario. Therefore, it can be seen that the parasitic fractional capacitor also has the nonlinear effect similarly to the conventional parasitic capacitor but with dependency on α. When α approaches 1, such effect yields merely a little loop area enlargement plus more roundish phase portraits similarly to the nonlinear effect of the conventional parasitic capacitor [9]. When α approaches 0, the nonlinear effect of the parasitic fractional capacitor yields the significant enlargement in the loop area. Therefore, it can be stated that the nonlinear effect of the parasitic element also exists in the fractional order domain where the fractional element can be assumed.Now we let the parasitic fractional inductor be nonlinear. With the nonlinear conventional parasitic inductor, the dynamic of memristor-capacitor circuit can be analytically described by the following ODEs:(55)dΦCdt=−qC,(56)dqdt=ΦC−fqLq,where Lq is the inductance of the nonlinear parasitic inductor.Therefore, the system of FDE of the memristor-capacitor circuit affected by such nonlinear parasitic fractional inductor can be given as follows:(57)dΦCdt=−qC,(58)Dαq=ΦC−fqLαq,where Lαq, which is a function of q, denotes the inductivity of the nonlinear parasitic fractional inductor. Noted also that qLα=q due to the circuit topology.If we letLαq=Lα1+q2 where Lα = Lα0, both q and ΦC can also be analytically given by (43) and (44) but with(59)q1=−tαΦC0+q0−q03/3Lα1+q02Γα+1,q2=tα9Lα21+q023q04+6q02+6q0ΦC0−3ΦC0tαq03−3q0−3ΦC0−1Γ2α+1−9Lαq0t1+q022CΓα+2,⋮(60)ΦC1=−q0tC,ΦC2=tα+1q03−3q0−3ΦC03CLα1+q02Γα+2.⋮after applying the Adomian decomposition.By numerically solving (57) and (58) via the moments of a function-based integer derivative-based approximation of Dαq with N = 7 where q0 = 0.1, ΦC0 = 0.1, C = 1, and Lα = 0.01 have been assumed, q, ΦC, and their phase portraits have been simulated for various values of α as depicted in Figures 18–20 where the nonlinear conventional parasitic inductor-related results have also been included for comparison. Noted also that other initial values are automatically assumed to be 0 due to the adopted approximation of Dαq. Moreover, (55) and (56) with Lq=L01+q2 have been used for simulating such nonlinear conventional parasitic inductor-related results.Figure 18
q vs t of the memristor-capacitor circuit with the nonlinear parasitic fractional inductor/nonlinear conventional parasitic inductor.Figure 19
Φ
C vs t the memristor-capacitor circuit with the nonlinear parasitic fractional inductor/nonlinear conventional parasitic inductor.Figure 20
Φ
C vs q the memristor-capacitor circuit with the nonlinear parasitic fractional inductor/nonlinear conventional parasitic inductor.From these figures, the continuous oscillations ofq and ΦC along with the impasse points breaking of the phase portraits between q and ΦC can be observed from both nonlinear conventional parasitic inductor and nonlinear fractional parasitic inductor-related results albeit the properties of the latter are α dependent because our memristor-capacitor-nonlinear parasitic fractional inductor circuit acts as a nonlinear fractional oscillator which is a nonlinear fractional system. Therefore, it can be seen that the parasitic fractional inductor also employs the nonlinear effect similarly to the conventional parasitic inductor but with dependency on α. By this and the aforementioned nonlinear effect of the parasitic fractional capacitor, the existence of the nonlinear effect of parasitic element in the fractional order domain has been validated. The phase portraits between q and ΦC with wider loop areas compared to those depicted in Figure 14 can be seen where the significant loop area enlargement can be observed when α approaches 0 and more roundish shapes of the phase portraits can also be found when α approaches 1. This is because the formerly neglected nonlinearities have been now considered.


## 6. Conclusion

The dynamic of the memristor under the effects of the parasitic fractional elements has been studied in this work based on the FDE approach where bothφ-controlled and q-controlled memristors have been considered. The comparison to the previous work which studied such device affected by the parasitic conventional elements [9] has also been made. The analytical solutions of the formulated FDEs have been determined by using the Adomian decomposition method where the numerical solutions have been obtained by using the integer derivative-based approximation of the fractional derivative with the moments of a function. From the obtained solutions and simulations, we have found that the decaying of φ and q and the impasse point breaking of the phase portraits are also occurred by the effect of the parasitic fractional elements as well as the parasitic conventional elements as reported by [9]. It has also been found that the decaying rates of φ and q are inversely proportional to α which also affects both size and shape of phase portraits and the properties of the oscillatory outputs of the memristor-based circuits which acts as the fractional Van der Pol oscillators. Moreover, we have found that the memory and the linearity are, respectively, inversely and directly proportional to α and the nonlinear effect of parasitic element, which occurs in the conventional integer domain, remains valid in the fractional order domain where the influence of such effect is α dependent. From these results, it can be seen that this research gives a further understanding on the characteristics of memristor which is a state-of-the-art electrical circuit element. Since this work limitedly focuses on the theoretical memristor similarly to [9], a similar study on the practical memristor, e.g., NTC/PTC thermistor [35] and HP memristor, which also suffers the adverse effects of parasitic elements as reported in previous studies, e.g., [35–37], should be considered as a tentative future work.



---

*Source: 4190641-2019-06-27.xml*