# On the Memristances, Parameters, and Analysis of the Fractional Order Memristor.

1. IntroductionRecently, a state-of-the-art electrical circuit element, namely, fractional order memristor, is often cited. This circuit element can be obtained from the generalization of the 4th electrical circuit element, namely, memristor, that has been theoretically found by Leon Chua since 1971 [1], by using the concept of fractional calculus which have been adopted in various disciplines, e.g., biomedical engineering [2,3], control system [4-6], and electronic engineering [7-9]. For decades after Chua proposed his original work, the memristor has been practically realized by a research group in Hewlett Packard (HP) labs [10] in 2008. As a result, the mathematical modelling and analysis attempts of the memristor have been proposed (e.g., [11-15]).

For the fractional order memristor on the other hand, there also exists such modelling and analysis attempts [16-21]. Some of them generalize the memristor by applying concept of the fractional calculus to the voltage-current relationship [16, 17] and termed such generalized memristor as the fracmemristor [17]. On the other hand, others do so by applying the fractional calculus to the memristor's state equation where the often cited HP memristor has been adopted as the basis [18-21]. However, only the analytical expression of the area of voltage-current hysteresis loop has been proposed in [20] and those of the memristances proposed in [18, 19, 21] are in terms of the input voltage despite the fact that fractional order memristor of interested is the generalization of the HP memristor which is actually of a charge/current controlled type. Moreover, these previous works also neglected the dimensional consistency [22, 23] related issues and the boundary effect, which is an important characteristic of the HP memristor [10], has not been analytically modelled.

By this motivation, we generalize the HP memristor in the fractional order domain by also concerning the formerly ignored dimensional consistency and formulate the analytical expression of memristance in term of the input current where boundary effect has also been modelled. We also derive the expressions of those related parameters of the fractional memristor excited by various exciting waveforms including DC with arbitrary delay and sinusoidal and cosinusoidal with arbitrary phase and arbitrary periodic which are the AC waveforms. With these expressions, parameters, and numerical simulations with MATHEMATICA, the behaviors of the fractional order memristor have been thoroughly explored. Unlike [18-21], both transient and asymptotic behaviors of the fractional order memristor excited by AC waveform have been distinguished and analyzed. The effect of phase of AC waveform has also been studied. Moreover, the influence of the fractional order to the areas of voltage-current hysteresis loop and memristance-current lissajous curve has been clearly discussed and the usage of fractional order memristor in the memristor based circuit has also been demonstrated.

In the following section, the overview of memristor will be briefly given followed by the memristor's generalization and derivation of our expressions in Section 3 where the behavioral analysis of the fractional order memristor will also be given. The DC waveform will be firstly treated followed by the AC waveforms where the sinusoidal waveform has been emphasized as it is the most fundamental. This is because the memristances and parameters of the memristor excited by the cosinusoidal and arbitrary periodic waveform can be obtained by using those due to the sinusoidal waveform as the basis as will be shown in Section 3 as well. The usage of fractional order memristor in the memristor based circuit will be presented in Section 4 and the conclusion will be finally drawn in Section 5.

2. The Overview of Memristor

Memristor is a nonlinear electrical circuit element. This circuit element relates the instantaneous flux, [phi](t), and charge, q(t), through the following relationship:

M(t) = d[phi](t)/dq(t) (1)

where M(t) denotes the memristance.

According to [10], M(t) of the HP memristor can be given in terms of the minimum and maximum values of M(t) denoted by [M.sub.on] and [M.sub.off] and the state variable, x(t), as

M(t) = [M.sub.on]x(t) + (1 - x(t))[M.sub.off] (2)

where x(t), which is dimensionless, can be given in terms of the memristor's current, i(t), by

dx/dt = ki(t) (3)

Note that k = [mu][M.sub.on]/[D.sup.2], where [mu] and D, respectively, stand for the ion mobility and semiconductor film of thickness. Therefore, the dimension of k is [(Asec).sup.-1].

As can be seen from (3) and also according to [10], x(t) can be simply given as follows:

x(t) = kq(t) (4)

Therefore, it can be seen that the HP memristor is charge controlled. Since q(t) is a time integration of i(t), it can be stated that the HP memristor is of a current controlled type. Note also that 0 [less than or equal to] x(t) [less than or equal to] 1; thus, [M.sub.on] [less than or equal to] M(t) [less than or equal to] [M.sub.off] as long as the memristor is unsaturated. Otherwise, x(t) will be bounded at either 0 or 1 so M(t) will be equal to either [M.sub.on] or [M.sub.off] according to the boundary effect of the device. Traditionally, such boundary effect can be mathematically modelled by multiplying the RHS of (3) with the window function [24].

3. The Fractional Order Domain Generalization of the Memristor and the Memristances, Related Parameters, and Analysis of the Fractional Order Memristor

By generalizing the memristor in the fractional order domain with the fractional calculus, the fractional order memristor can be obtained. Similarly to [18-21], we perform such generalization by applying the fractional calculus to the memristor's state equation, i.e., (3). In these previous works, [d.sup.[alpha]]x/d[t.sup.[alpha]] = ki(t), where [alpha] stands for the order of the fractional order memristor which can be arbitrary real value and has been obtained from such generalization. However, as x(t) is dimensionless; the dimension of the LHS of this previous generalized equation is given by [sec.sup.-[alpha]] where that of the RHS is [sec.sup.-1] which means that a dimensional inconsistency has always existed.

Therefore, the fractional time component [22], [sigma], which has the dimension of sec, has been introduced for handling this issue. As a result, unlike [18-21], the following generalized state equation has been used instead.

[[sigma].sup.[alpha]-1][[d.sup.[alpha]]x/d[t.sup.[alpha]]] = ki(t) (5)

Similarly to that of the RHS, the dimension of the LHS of (5) is [sec.sup.-1]; thus the dimensional inconsistency issue has been resolved. Note also that (5) is reduced to (3) when [alpha] = 1 despite the presence of a as [[sigma].sup.[alpha]-1] become 1 with such value of [alpha].

Unlike [18, 19, 21], we derive M(t) of the fractional order memristor as a function i(t) as it has been assumed that the of fractional order memristor is a generalization of the HP memristor which is of a current controlled type as aforementioned. Therefore we directly determine x(t) from (5) by using the Riemann-Liouville fractional order integral [25] as follows:

[mathematical expression not reproducible] (6)

where x(0) and [GAMMA]() denote the initial value of x(t) and the Gamma function [26], respectively.

Since it can be seen from (2) that

M(t) = [M.sub.off] - [M.sub.d]x(t) (7)

where [M.sub.d] = [M.sub.off] - [M.sub.on], the initial memristance value, i.e., M(0), can be immediately given by

M(0) = [M.sub.off] - [M.sub.d]x(0) (8)

Thus by substituting (6) into (7) and keeping (8) in mind, M(t) of the fractional order memristor can be obtained as follows:

[mathematical expression not reproducible] (9)

which shows that M(t) is current-controlled.

If we let [alpha] = 1, (9) will be reduced to

M(t) = M(0) - k[M.sub.d][[integral].sup.t.sub.0]i([tau])d[tau] (10)

Since the integer order integration of i(t) gives q(t), we obtain

M(t) = M(0) - k[M.sub.d]q(t) (11)

By using (8), (11) can be simplified under the assumption that x(0) = 0 and [M.sub.off] [much greater than] [M.sub.on] as follows:

M(t) = [M.sub.off](1 - kq(t)) (12)

which is similar to the original simplified model of the HP memristor [10]. Such correspondence cannot be found in [18, 19, 21] as the integer order integration of the memristor's voltage, v(t), yields [phi](t).

For traditionally including the boundary effect, the window function must be introduced to the state equation as mentioned above. In [18], the linear window function given by f(x(t)) = 1 has been adopted for simplicity as the usage of more accurate yet more complicated window function; e.g., those Joglekar, Biolek, and Prodomakis [24] can be mathematically cumbersome. Unfortunately, using such linear window function is mathematically equivalent to multiplying the RHS of the state equation by 1. As a result, no modification has been made on the state equation; thus the boundary effect modelling has not been performed. Moreover, neither the usage of window function nor alternative boundary effect analytical modelling has been made in both [19] and [21].

In order to model the boundary effect in a simplified manner, we apply two mathematical operators, i.e., max[x, y] and min[x, y], which, respectively, selects the maximum value and minimum value among x and y, to (9). As a result, our expression of M(t) due to arbitrary exciting waveform can be finally given as follows:

[mathematical expression not reproducible], (13)

where K = k[[sigma].sup.1-[alpha]]; thus the dimension of K is [A.sup.-1] [sec.sup.-[alpha]].

Since the dimension of fractional integral of i(t) is [Asec.sup.[alpha]], that of (K[M.sub.d]/[GAMMA]([alpha]))[[integral].sup.t.sub.0] [(t - [tau]).sup.[alpha]-1]i([tau])d[tau] is given by [OMEGA] which is physically measurable, similarly to those of [M.sub.d] and M(0). Therefore, M(0) and (K[M.sub.d]/[GAMMA]([alpha]))[[integral].sup.t.sub.0] [(t - [tau]).sup.[alpha]-1]i([tau])d[tau], which are at the RHS of (13), can be physically combined as they have the same dimensions and the dimension of M(t),which is the LHS of such equation, has also been found to be such physically measurable [OMEGA]; thus it can be seen that our expression of M(t) has dimensional consistency. Moreover, due to the operation of nested max[x, y] and min[x, y], M(t) will be equal to M(0) - (K[M.sub.d]/[GAMMA]([alpha]))[[integral].sup.t.sub.0] [(t - [tau]).sup.[alpha]-1]i([tau])d[tau] if and only if M(0) - (K[M.sub.d]/[GAMMA]([alpha]))[[integral].sup.t.sub.0] [(t - [tau]).sup.[alpha]-1]i([tau])d[tau] lies within [[M.sub.on], [M.sub.off]] which means that the fractional order memristor remains unsaturated. Otherwise, M(t) will be equal to either [M.sub.on] or [M.sub.off] if M(0) - (K[M.sub.d]/[GAMMA]([alpha]))[[integral].sup.t.sub.0][(t - [tau]).sup.[alpha]-1]i([tau])d[tau] is lower than [M.sub.on] or higher than [M.sub.off] which in turn means that the device become saturated at either its on-state or off-state. Therefore it can be seen that the boundary effect has been modelled without any necessity to use the window function and (13) along with its related results is valid to the saturated fractional order memristor. In the following subsections, M(t)'s due to due to various exciting waveforms and the behavioral analysis of fractional order memristor will be presented.

3.1. DC Waveform. Mathematically, the DC waveform with arbitrary delay ([t.sub.d]), which is more generic than the undelay waveform assumed in the previous works [18, 19], can be defined as i(t) = [I.sub.DC]u(t - [t.sub.d]), where [I.sub.DC] and u(t) denote the magnitude of the waveform and the unit step function. Therefore, the resulting M(t) can be straightforwardly obtained by using (13) as

[mathematical expression not reproducible], (14)

By using (14) with [t.sub.d] = 0 sec, K = 100000 [A.sup.-1] [sec.sup.-[alpha]], [M.sub.on] = 1 k[ohm], [M.sub.off] = 100 k[ohm], and M(0) = 80 k[ohm], M(t)'s of the fractional order memristor with various as excited by the DC waveform can be numerically simulated as depicted in Figures 1 and 2 where [I.sub.DC] = 110 [micro]A and [I.sub.DC] = -110 [micro]A have been, respectively, assumed and M(t)'s of the HP memristor simulated by using its SPICE model [27] have also been included.

From these figures, the strong agreements between our M(t)'s obtained by using (14) with [alpha] = 1 and those of the HP memristor can be observed. Since the HP memristor is of order 1 in the context of fractional order domain, such strong agreements verify our expression. These figures also show that M(t) of the fractional order memristor can be either the increasing or decreasing function of t with the final value of [M.sub.on] or [M.sub.off] when the fractional order memristor become saturated if we let [I.sub.DC] > 0 or [I.sub.DC] < 0. Moreover, the rate of change of M(t) has been found to be inversely proportional to [alpha].

Since it can be seen that the fractional order memristor become saturated at a certain time given by [t.sub.sat], by using (14), we have

[mathematical expression not reproducible] (15)

As a result, [t.sub.sat] can be immediately given as follows:

[mathematical expression not reproducible] (16)

which shows that [t.sub.sat] is directly proportional to the size of the difference between M(0) and [M.sub.on]([M.sub.off]). So, [t.sub.sat] reaches its maximum value, i.e., [t.sub.sat(MAX)], if and only if M(0) reaches its possible peak value given by either [M.sub.off] when [I.sub.DC] > 0 or [M.sub.on] when [I.sub.DC] < 0. Thus [t.sub.sat(MAX)] can be found as

[mathematical expression not reproducible] (17)

which can be immediately given in a more compact manner as follows:

[mathematical expression not reproducible] (18)

For confirming the aforesaid observation on the relationship between the rate of change of M(t) and [alpha], [t.sub.sat(MAX)]'s have been simulated by using (18) under the similar assumptions to those of the simulation of M(t)'s shown in Figures 1 and 2 but with varying [I.sub.DC], as depicted in Figure 3 which shows that [t.sub.sat(MAX)] is directly proportional to a. Due to the definition of [t.sub.sat(MAX)], this confirms such observation. Moreover, it can also be seen from Figure 3 that [t.sub.sat(MAX)] is inversely proportional to [absolute value of [I.sub.DC]].

3.2. AC Waveforms. Among various AC waveforms, the sinusoidal waveform has been emphasized as it is the foundation of the others. Unlike those previous works, the sinusoidal waveform with arbitrary phase (0) has been chosen due to its generality. Mathematically, such waveform can be given by i(t) = [I.sub.m]sin([omega]t + [theta]), where [I.sub.m] and [omega], respectively, denote the peak value and angular frequency of i(t). By using (13), M(t) of the fractional order memristor under the arbitrary phase sinusoidal input can be given by

[mathematical expression not reproducible] (19)

where [sub.1][F.sub.2](;, ;) denotes a generalized hypergeometric function with p =1 and q = 2 [28].

By letting K = 100000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 1 k[ohm], [M.sub.off] = 100 k[ohm], and M(0) = 80 k[ohm] similarly to the previous subsection, M(t)'s of the fractional order memristor with various [alpha]'s excited by the sinusoidal input can be simulated as depicted in Figures 4 and 5, where [I.sub.m] = 110 [micro]A and [I.sub.m] = -110 [micro]A have been, respectively, assumed and M(t)'s of the HP memristor simulated by using its SPICE model have also been included. Moreover, we also let [omega] = 1 rad/sec and [theta] = 0 rad. Again, the strong agreements between our expression based M(t)'s with [alpha] =1 and those of the HP memristor which verify our expression can be observed. It can be seen that M(t) of the fractional order memristor with [alpha] [less than or equal to] 1 is periodic with clipped peaks at [M.sub.on] and [M.sub.off] due to the temporary saturation of the fractional order memristor as the unclipped peaks of M(t) lie outside [[M.sub.on], [M.sub.off]]. On the other hand, M(t) of the device with a >1 is not periodic but time independently equal to either [M.sub.on] or [M.sub.off] when t [greater than or equal to] [T.sub.sat] up to the sign of [I.sub.m]. This is because the device becomes permanently saturated after [T.sub.sat] as M(t) contains the time proportional term which starts to lie outside [[M.sub.on], [M.sub.off]] at t = [T.sub.sat] when [alpha] > 1. In order to show that the permanent saturation of the fractional order memristor with [alpha] > 1 excited by sinusoidal input is possible, we simulate M(t) of the device with [alpha] = 1.25 once again by using (20) which is (19) without the boundary effect. The simulation results have been depicted in Figures 6 and 7 which show that M(t)'s eventually lie outside [[M.sub.on], [M.sub.off]] permanently. This yields the saturation when such effect has been included.

[mathematical expression not reproducible] (20)

However, it should be mentioned here that both temporary and permanent saturation do not always occur. Instead, their occurrences are dependent on the conditions on parameters, which makes either the peaks or time proportional term of M(t) be outside [[M.sub.on], [M.sub.off]] such as those of Figures 4 and 5. Moreover the determination of such time proportional term and [T.sub.sat] will be presented later in this subsection.

Similarly to the memristor, the voltage-current lissajous curve is of interest for the fractional order memristor. Therefore, those curves of the fractional order memristor will be simulated by using (19) as the basis for further studying the behavior of this device under the AC input. Now, we let 0 sec < t [less than or equal to] 19 sec, [omega] = 1 rad/sec, and [theta] = 0 rad; the lissajous curves of the fractional order memristor with various a's can be simulated as depicted in Figures 8-13 where [I.sub.m] = 110 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 16 k[ohm], and M(0) = 11 k[ohm] have been assumed in Figures 8-10. On the other hand, [I.sub.m] = -150 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 38 k[ohm], and M(0) = 11.2 k[ohm] have been adopted in Figures 11-13. Based on its SPICE model, the lissajous curves of the HP memristor have also been simulated and compared to those of the fractional order memristor with [alpha] =1 as depicted in Figures 9 and 12 where the strong agreements between the fractional order memristor and HP memristor based curves which are both time invariant shaped symmetric pinched hysteresis loops can be observed. Again, this verifies the accuracy of our expression.

From other figures, it can be seen that the lissajous curve of the fractional order memristor with [alpha] [not equal to] 1 also takes the pinched hysteresis loop shape with pinching point at the origin despite the asymmetricities which means that the fractional order memristor preserves the memristive characteristic [30]. This is unlike those fracmemristor based results previously proposed in [17] which do not display the pinched hysteresis loop at all; therefore such fracmemristor does not employ the memristive characteristic according to [30]. Unlike those of the fractional order memristor with [alpha] = 1 and HP memristor, the shape of the lissajous curve of the fractional order memristor with [alpha] [not equal to] 1 keeps changing. This is because M(t) of such fractional order memristor contains the time proportional term. This can be clearly seen from Figures 14-17 which display M(t)'s of the fractional order memristor with [alpha] [not equal to] 1. Unlike the previous ones depicted in Figures 4 and 5, the peaks of these M(t)'s with [alpha] < 1 are unclipped as can be seen from Figures 14 and 16. This is because fractional order memristor is always unsaturated as such peaks lie within [[M.sub.on], [M.sub.off]]. Moreover, the fractional order memristor with [alpha] > 1 never becomes saturated and can be seen from Figures 15 and 17. This is because the time proportional terms of its M(t)'s are always be within [[M.sub.on], [M.sub.off]]. At this point, it can be seen that both temporary and permanent saturation of the fractional order memristor under sinusoidal excitation do not always occur but depend on the conditions on parameters and input as aforementioned.

At asymptotic state, it can be seen that such time proportional term of M(t) becomes time independent instead. Therefore, the shape of its lissajous curves of the fractional order memristor with [alpha] [not equal to] 1 is asymptotically unchanged as can be seen from Figures 18 and 19 which depict the lissajous curves of the fractional order memristor with various [alpha]'s simulated by assuming that 100 sec < t [less than or equal to] 119 sec. Apart from being asymmetric when [alpha] = 1, we have found that the lobe area of the lissajous curve of the fractional order memristor is affected by [alpha] similarly to those proposed in [20, 21]. In particular, we have found that the fractional order memristor with lower [alpha] yields the lissajous curve with wider lobe area which refers to more memory effect and less linearity. Besides the voltage-current curve, the memristance-current lissajous curve has been found to be also interesting. Therefore, such curves of the fractional order memristor with various [alpha]'s have also been simulated as depicted in Figures 20 and 21 where the parameters setting similar to those of the voltage-current curve have been assumed. Moreover, we also assume that 100 sec < t [less than or equal to] 119 sec. These figures show that the resulting lissajous curves take the elliptical closed loop shape which are unchanged as M(t) enters the asymptotic state at the assumed time interval. The elliptical shaped lissajous curve means that M(t) is periodic as well as i(t). Similarly to the voltage-current curve, the loop area of the memristance-current lissajous curve is also affected by [alpha]. In particular, the fractional order memristor with lower [alpha] yields the memristance-current lissajous curve with wider loop area.

At this point, we will analytically show that both upper and lower lobes of the voltage-current lissajous curve have equal sizes of areas which means that the fractional order memristor does not store the energy, and such areas are independent of [theta]. Here, we let the area of the upper and lower lobes of the lissajous curve be denoted, respectively, by [A.sub.U] and [A.sub.L]. Since v(t) = [M.sub.a](t)i(t) at asymptotic state where [M.sub.a](t) denotes the asymptotic approximation of M(t), [A.sub.U] and [A.sub.L] can be given by

[mathematical expression not reproducible] (21)

[mathematical expression not reproducible] (22)

With (13) and the asymptotic approximation of the fractional order integration of sinusoidal function [20], [M.sub.a](t) can be found as

[mathematical expression not reproducible] (23)

As i(t) = [I.sub.m] sin([omega]t+[theta]) and the fractional order memristor is unsaturated for the entire simulation period as can be seen from Figures 18 and 19, [A.sub.U] and [A.sub.L] can be finally obtained as follows:

[mathematical expression not reproducible] (24)

[mathematical expression not reproducible] (25)

which show that the upper and lower lobes of the lissajous curve employ equal sizes of areas and such areas are independent of [theta]. Unlike [20] which considered only [A.sub.U], [A.sub.L] has also been formulated in this work. Moreover, our [A.sub.U] which has been derived by using a different approach from that used in [20] totally agrees with the previous result that [theta] = 0 rad has been assumed. This emphasizes the independency from [theta] of the lobe areas. Moreover, (24) and (25) also show that the fractional order memristor does not store the energy as the summation of [A.sub.U] and [A.sub.L] which, respectively, referred to the intake and dissipated energy, is equal to 0.

Besides [A.sub.U] and [A.sub.L], the area within the closed loop of the memristance-current lissajous curve ([A.sub.M]) can be obtained by also using (23) and

[mathematical expression not reproducible] (26)

Since it can be seen from Figures 20 and 21 that the device is unsaturated for the whole simulation period and i(t) = [I.sub.m] sin([omega]t + [theta]), we have

[A.sub.M] = [K[M.sub.d][I.sup.2.sub.m][pi]/[[omega].sup.[alpha]]] sin ([alpha][pi]/2) (27)

which shows that [A.sub.M] is affected by a but independent of [theta] as well as [A.sub.U] and [A.sub.L].

Despite the aforementioned independencies, [theta] does affect the behavior of the fractional order memristor. For illustration, [M.sub.a](t)'s due to various [theta]'s have been simulated by using (20) as depicted in Figures 22 and 23 where [alpha] = 1 and [omega] = 1 rad/sec has been assumed. It should be mentioned here that [I.sub.m] = 110 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 16 k[ohm], and M(0) = 11 k[ohm] have been adopted in Figure 20 where [I.sub.m] = -150 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 38 k[ohm], and M(0) = 11.2 k[ohm] have been assumed in Figure 21. Moreover, [theta] = [pi]/4 rad, [theta] = 3[pi]/4 rad, [theta] = -3[pi]/4 rad, and [theta] = -[pi]/4 rad have been chosen as they are good representatives of those [theta]'s which their coordinates on the Euclidian plane, i.e., (cos([theta]), sin([theta])), are located on the portion of unit circle's arc in quadrants 1,2, 3, and 4 of such plane, respectively. This is because (cos([theta]), sin([theta]))'s of these chosen [theta]'s are exactly located at the middle points of the portion of unit circle's arc. For example, (cos([theta]), sin([theta])) of [theta] = [pi]/4 rad is located at the middle point of the portion of unit circle's arc in quadrant 1, etc.

From both figures, it can be seen that these [M.sub.a](t)'s contain time independent terms which are formerly the time proportional term of M(t) that become time independent at asymptotic state as aforementioned. Since these time independent terms lie within [[M.sub.on], [M.sub.off]], [M.sub.a](t) takes the shape of sinusoidal waveform. However, the minimum peaks of [M.sub.a](t)'s due to the input with [I.sub.m] < 0 and (cos([theta]), sin([theta])) located on the portion of unit circle's arc in quadrants 2 and 3 of the Euclidian plane have been clipped as can be seen from Figure 23 due to the saturation of the fractional order memristor as these peaks are lower than [M.sub.on]. From Figure 22, it can be seen that the input with [I.sub.m] > 0 and [theta] with (cos([theta]), sin([theta])) located on the portion of unit circle's arc in quadrants 2 and 3 yields M(t) with higher time average. If we have assumed that [M.sub.on] [less than or equal to] [M.sub.a](t) [less than or equal to] [M.sub.off] is always satisfied, it can be seen from Figure 23 that M(t) with higher time average can be obtained by using the input with [I.sub.m] < 0 and (cos([theta]), sin([theta])) located on such portion in quadrants 1 and 4.

Now, we will derive the analytical expression of time independent term of [M.sub.a](t). Let such term be denoted by [M.sub.TI]; it can be given by carefully observing (23) as follows:

[mathematical expression not reproducible], (28)

When [alpha] = 1 as assumed in our simulations of [M.sub.a](t), [M.sub.TI] can be reduced to

[M.sub.TI] = min[max [M(0) + [DELTA][M.sub.TI], Mon], [M.sub.off]] (29)

where

[DELTA][M.sub.TI] = -K[[omega].sup.-1][M.sub.d][I.sub.m]cos([theta]) (30)

Therefore, it can be seen that when [I.sub.m] > 0 ([I.sub.m] < 0), if (cos([theta]), sin([theta])) is located on the portion of unit circle's arc in either quadrant 2 or 3 where cos([theta]) < 0, [DELTA][M.sub.TI] > 0 ([DELTA][M.sub.TI] < 0) thus M(0) has been increased (decreased). On the other hand, if (cos([theta]), sin([theta])) is on the portion of unit circle's arc in either quadrant 1 or 4 where cos([theta]) > 0, [DELTA][M.sub.TI] < 0 ([DELTA][M.sub.TI] >0) as [I.sub.m] > 0 ([I.sub.m] < 0) thus M([theta]) has been decreased (increased). As a result, the sinusoidal input with [I.sub.m] > 0 ([I.sub.m] < 0) and (cos([theta]), sin([theta])) located on the portion of unit circle's arc in either quadrant 2 or 3 (1 or 4) of the Euclidian plane yields higher [M.sub.TI] as graphically observed.

At this point, the influence of a to [M.sub.TI] will be explored. By using (28) with [theta] = 0 rad and [omega] = 1 rad/sec, we can simulate [M.sub.TI] as shown in Tables 1 and 2. In Table 1, we assume that [I.sub.m] = 110 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 16 k[ohm], and M(0) = 11 k[ohm], where [I.sub.m] = -150 [micro]A, K = 10000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], [M.sub.off] = 38 k[ohm], and M(0) = 11.2 k[ohm] have been adopted in Table 2. It can be seen from these tables that [M.sub.TI] is proportional to [alpha] when [I.sub.m] > 0 and vice versa when [I.sub.m] < 0.

By subtracting [M.sub.TI] from [M.sub.a](t), we obtain the following purely periodic term of M(t), [M.sub.P](t)

[mathematical expression not reproducible] (31)

which existed in both initial state and asymptotic state.

Therefore, the time proportional term of M(t), [M.sub.TP](t) can be formulated by subtracting [M.sub.P](t) from M(t) as

[mathematical expression not reproducible] (32)

As a result, the aforementioned [T.sub.sat] can be determined by solving the following equation:

[mathematical expression not reproducible] (33)

At this point, we will derive the expressions of M(t)'s due to other AC waveforms by using that due to the sinusoidal waveform, i.e., (19), as the basis. For example, the expression of M(t) due to arbitrary phase cosinusoidal waveform, i.e., i(t) = [I.sub.m] cos([omega]t + [theta]), can be formulated by using (19) and the relationship between the sinusoidal and cosinusoidal functions is given by cos(u) = sin(([pi]/2) - u) where u denotes arbitrary angle. As a result, we have

[mathematical expression not reproducible] (34)

For arbitrary periodic waveform which has never been considered in those previous works, the resulting expression can also be determined based on (19). This is because such waveform can be given as a series of sinusoidal functions, i.e., i(t) = [[summation].sup.[infinity].sub.n=0][[I.sub.mn]sin(n[omega]t + [[theta].sub.n])], where [I.sub.mn] and [[theta].sub.n], respectively, stand for peak value and phase of arbitrary nth term of the series, according to the Fourier theorem. As a result, the expression of M(t) due to arbitrary periodic waveform can be given as follows:

[mathematical expression not reproducible] (35)

4. The Usage of Fractional Order Memristor in the Memristor Based Circuit

In this research, the HP memristor based type A Wien oscillator [29] has been chosen as the candidate memristor based circuit. The unique characteristic of such circuit, which is either [R.sub.1] or [R.sub.2] replaced by the memristive device as depicted in Figure 22, is the fluctuated frequency of oscillation. The smaller range of fluctuation refers to the better chance that the system has sustained oscillation which can be obtained if and only if all poles of the system are fixed in the s-plane [29]. For studying the usage of fractional order memristor, we replace [R.sub.1] of the circuit by such fractional order device instead of the HP memristor as traditionally did [29] and analyze the effect of [alpha] to the fluctuation in frequency of oscillation which determines the chance that the system has sustained oscillation, as mentioned above. Let the range of such fluctuation be denoted by [DELTA][f.sub.osc]; it can be mathematically defined as given by (36) where [f.sub.up] and [f.sub.low] stand for the upper bound and lower boundary.

[DELTA][f.sub.osc] = [f.sub.up] - [f.sub.low] (36)

According to [29], these boundaries can be obtained by solving (37) and (38) which have been formulated by assuming that [C.sub.1] = [C.sub.2] = C, the memristor is unsaturated, and the frequency of oscillation has been found to lie within the range that the sustained oscillation can be assured. Note also that [I.sub.mem] = [V.sub.mem]/M(0), where [V.sub.mem] denotes the peak value of the voltage dropped across the memristor which can be determined from the oscillator output voltage, [V.sub.out](t), and also depends on the initial voltages of [C.sub.1] and [C.sub.2] [29]. Moreover, [V.sub.out](t) can be obtained from the state space representation and output equation of the fractional memristor based Wien oscillator which are, respectively, given by (39) and (40) where [V.sub.C1](t) and [V.sub.C2](t) denote the voltage dropped across [C.sub.1] and [C.sub.2] and M(t) stands for the memristance of the conventional HP memristor.

[mathematical expression not reproducible] (37)

[mathematical expression not reproducible] (38)

[mathematical expression not reproducible] (39)

[mathematical expression not reproducible] (40)

However, this is not the case when the fractional order memristor has been used as M(t) of such device must be adopted. Therefore [V.sub.out](t) must be determined based on our derived M(t) instead where [f.sub.up] and [f.sub.low] must be evaluated from (41) and (42). As a result, the corresponding [DELTA][f.sub.osc] will be different from that of the original conventional memristor based circuit and the different chance of obtaining sustained oscillation can be expected.

[mathematical expression not reproducible] (41)

[mathematical expression not reproducible] (42)

Moreover, the condition for ensuring the occurrence of sustained oscillation can be given by () where [f.sub.avr] which stands for the average oscillating frequency is given by ().

[mathematical expression not reproducible] (43)

[mathematical expression not reproducible] (44)

By assuming 0.1 V and -0.95 V as the approximate initial voltages of [C.sub.1] and [C.sub.2] and also assuming that [R.sub.2] = 5 k[ohm], [C.sub.1] = [C.sub.2] = 3.2[micro]F, K = 100000 [A.sup.-1][sec.sup.-[alpha]], [M.sub.on] = 100 [ohm], and [M.sub.off] = 16 k[ohm], [DELTA][f.sub.osc] of the fractional order memristor based Wien oscillator with [alpha] < 1 and [alpha] > 1 can be obtained for various M(0)'s by numerically solving (41) and (42). as tabulated in Table 3 where [DELTA][f.sub.osc] of the original HP memristor based circuit which is equivalent to the fractional order memristor with [alpha] =1 in the context of this work, determined by solving (37) and (38), has also been included.

It has been found that [DELTA][f.sub.osc] is inversely proportional to M(0) which is in agreement with [29]. Since 4.1 k[ohm] [less than or equal to] M(0) [less than or equal to] 5.9 k[ohm] [29] and the probability of obtaining sustained oscillation is inversely proportional to [DELTA][f.sub.osc], M(0) = 5.9 k[ohm] is recommended as it minimizes [DELTA][f.sub.osc] and thus maximizes such probability. It can also be seen that [DELTA][f.sub.osc] is inversely proportional to a which means that the probability of obtaining sustained oscillation is directly proportional to [alpha]. Therefore, the fractional order memristor with larger [alpha] is recommended and the fractional order memristor with [alpha] > 1 should be adopted as it can increase such probability from that of the original circuit which its memristive device has [alpha] = 1. On the other hand, the fractional order memristor with [alpha] < 1 should be avoided as it decreases such probability.

Finally, by further assuming that M(0) = 5 k[ohm], [R.sub.3] = 20.2 k[ohm], and [R.sub.4] = 10 k[ohm], we can simulate [V.sub.out](t) at asymptotic state and the lissajous patterns of [V.sub.out](t) and M(t) for various [alpha]'s as depicted in Figures 23-28 where the units of 100 [V.sub.out](t), t, and M(t) are V, sec, and [ohm], respectively. The comparison of the results with [alpha] = 1 to the SPICE HP memristor model based counterparts has been made for verification where a strong agreement can be observed. We have found that there exists neither temporary nor permanent saturation of the fractional order memristor. This is because M(t) is always within [[M.sub.on], [M.sub.off]] in this scenario due to the assumed conditions on parameters and input, as can be seen from the simulated lissajous patterns. Such conditions have been adopted for ensuring the unsaturation which yields the proper operation of the oscillator [29]. From these lissajous patterns, we have also found that [V.sub.out](t) is of the same frequency as M(t) and the phase difference between [V.sub.out](t) and M(t), which is less than 90[degrees], is inversely proportional to

5. Conclusion

In this work, the HP memristor has been generalized in a fractional order domain by applying the fractional calculus to its state equation. Unlike [18-21], the dimensional consistency has been taken into account. Moreover, the boundary effect has also been modelled. Therefore the analytical expression of M(t) which has been derived as a function of i(t) and its related results are dimensional consistent and valid to such generalized device in its saturation states. By using such expression, M(t)'s due to various waveforms and the related parameters have been formulated. With the simulations by using these M(t)'s and parameters, the behaviors of the fractional order memristor have been thoroughly explored. Therefore this research gives a precise understanding on the characteristics of such up to date nonlinear electrical circuit element which has been recently applied as the basis of the net grid type fracmemristor [31].

https://doi.org/10.1155/2018/3408480

Data Availability

The simulated data used to support the findings of this study are included within the article.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The author would like to acknowledge Mahidol University, Thailand, for the online database service which is our primary information resource.

References

[1] L. O. Chua, "Memristor--the missing circuit element," IEEE Transactions on Circuit Theory, vol. 18, no. 5, pp. 507-519, 1971.

[2] I. S. Jesus, J. A. Tenreiro MacHado, and J. Boaventure Cunha, "Fractional electrical impedances in botanical elements," Journal of Vibration and Control, vol. 14, no. 9-10, pp. 1389-1402, 2008.

[3] C. Tang, F. You, G. Cheng, D. Gao, F. Fu, and X. Dong, "Modeling the frequency dependence of the electrical properties of the live human skull," Physiological Measurement, vol. 30, no. 12, pp. 1293-1301, 2009.

[4] A. Charef, "Analogue realisation of fractional-order integrator, differentiator and fractional PIAD[micro] controller," IEE Proceedings--Control Theory and Applications, vol. 153, no. 6, pp. 714-720, 2006.

[5] B. M. Vinagre and V. Feliu, "Optimal fractional controllers for rational order systems: a special case of the Wiener-Hopf spectral factorization method," Institute of Electrical and Electronics Engineers Transactions on Automatic Control, vol. 52, no. 12, pp. 2385-2389, 2007.

[6] R. Matusu, "Application of fractional order calculus to control theory," International Journal of Mathematical Models and Methods in Applied Sciences, vol. 5, pp. 1162-1169, 2011.

[7] L. Dork, J. Terpk, I. Petr, and F. Dorkov, "Electronic realization of the fractional-order systems," Acta Montanistica Slovaca, vol. 12, pp. 231-237, 2007.

[8] B. T. Krishna, K. V. V. S. Reddy, and S. Santha Kumari, "Time domain response calculations of fractance device of order 1/2," Journal of Active & Passive Electronic Devices, vol. 3, pp. 355-367, 2008.

[9] A. G. Radwan and A. S. Elwakil, "An expression for the voltage response of a current-excited fractance device based on fractional-order trigonometric identities," International Journal of Circuit Theory and Applications, vol. 40, no. 5, pp. 533-538, 2012.

[10] D. B. Strukov, G. S. Snider, D. R. Stewart, and R. S. Williams, "The missing memristor found," Nature, vol. 453, pp. 80-83, 2008.

[11] A. G. Radwan, M. A. Zidan, and K. N. Salama, "HP Memristor mathematical model for periodic signals and DC," in Proceedings of the 2010 53rd IEEE International Midwest Symposium on Circuits and Systems (MWSCAS), pp. 861-864, Seattle, Wash, USA, August 2010.

[12] A. G. Radwan, M. A. Zidan, and K. N. Salama, "On the mathematical modeling of Memristors," in Proceedings of the 2010 International Conference on Microelectronics, ICM'10, pp. 284-287, Cairo, Egypt, December 2010.

[13] Y. N. Joglekar and S. J. Wolf, "The elusive memristor: properties of basic electrical circuits," European Journal of Physics, vol. 30, no. 4, pp. 661-675, 2009.

[14] S. Shin, K. Kim, and S.-M. Kang, "Compact models for memristors based on charge-flux constitutive relationships," IEEE Transactions on Computer-Aided Design of Integrated Circuits and Systems, vol. 29, no. 4, pp. 590-598, 2010.

[15] C. Yakopcic, T. M. Taha, G. Subramanyam, R. E. Pino, and S. Rogers, "A memristor device model," IEEE Electron Device Letters, vol. 32, no. 10, pp. 1436-1438, 2011.

[16] J. Tenreiro Machado, "Fractional generalization of memristor and higher order elements," Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 2, pp. 264-275, 2013.

[17] Y.-F. Pu and X. Yuan, "Fracmemristor: fractional-order memristor," IEEE Access, vol. 4, pp. 1872-1888, 2016.

[18] M. E. Fouda and A. G. Radwan, "On the fractional-order memristor model," Fractional Calculus and Applied Analysis, vol. 4, pp. 1-7, 2013.

[19] M. E. Fouda and A. G. Radwan, "Fractional-order Memristor Response Under DC and Periodic Signals," Circuits, Systems and Signal Processing, vol. 34, no. 3, pp. 961-970, 2015.

[20] Y. Yu, B. Bao, H. Kang, and M. Shi, "Calculating area of fractional-order memristor pinched hysteresis loop," The Journal of Engineering, vol. 2015, no. 11, pp. 325-327, 2015.

[21] M. Shi and S. Hu, "Pinched hysteresis loop characteristics of a fractional-order HP Ti[O.sub.2] memristor," in Intelligent Computing, Networked Control, and Their Engineering Applications, vol. 762 of Communications in Computer and Information Science, pp. 705-713, Springer, Singapore, Singapore, 2017.

[22] J. F. Gomez-Aguilar, J. J. Rosales-Garcia, J. J Bernal-Alvarado, T. Cordova-Fraga, and R. Guzman-Cabrera, "Fractional mechanical oscillators," Revista Mexicana de Fisica, vol. 58, pp. 348-352, 2012.

[23] R. Banchuin, "Novel expressions for time domain responses of fractance device," Cogent Engineering, vol. 4, no. 1, 2017.

[24] A. G. Radwan and E. Mohammed, "Memristor: models, types, and applications," in On the Mathematical Modeling of Memristor, Memcapacitor, and Meminductor, pp. 13-49, Springer, Cham, Switzerland, 2015.

[25] J. Sabatier, O. P. Agrawal, and J. A. Machado, Advance in Fractional Calculus: Theoretical Developments and Applications in Physics and Engineering, Springer, Dordrecht, The Netherlands, 2007.

[26] W. H. Beyer, CRC Handbook of Mathematical Sciences, CRC Press, Boca Raton, Fla, USA, 1987.

[27] Z. Biolek, D. Biolek, and V. Biolkova, "SPICE model of memristor with nonlinear dopant drift," Radioengineering, vol. 18, no. 2, pp. 210-214, 2009.

[28] B. Dwork, Generalized Hypergeometric Functions, Clarendon Press, Oxford, UK, 1990.

[29] A. Talukdar, A. G. Radwan, and K. N. Salama, "Generalized model for Memristor-based Wien family oscillators," Microelectronics Journal, vol. 42, no. 9, pp. 1032-1038, 2011.

[30] L. Chua, "If it's pinched it's a memristor," Semiconductor Science and Technology, vol. 29, no. 10, Article ID 104001, 2014.

[31] L. Xu, G. Huang, and Y. Pu, "Numerical Simulation Research of Fracmemristor Circuit Based on HP Memristor," Journal of Circuits, Systems and Computers, vol. 27, no. 14, Article ID 1850227, 2018.

Rawid Banchuin (ID)

Faculty of Engineering and Graduated School of Information Technology, Siam University, Bangkok, Thailand

Correspondence should be addressed to Rawid Banchuin; rawid_b@yahoo.com

Received 26 July 2018; Revised 5 September 2018; Accepted 24 September 2018; Published 1 November 2018

Academic Editor: Stephan Gift

Caption: Figure 1: M(t) of the fractional order memristor excited by a 110 [micro]A DC waveform: [alpha] = 0.75 (red), [alpha] = 1 (green), [alpha] = 1.25 (blue), and HP memristor (black dots).

Caption: Figure 2: M(t) of the fractional order memristor excited by a -110 [micro]A DC waveform: [alpha] = 0.75 (red), [alpha] = 1 (green), [alpha] = 1.25 (blue), and HP memristor (black dots).

Caption: Figure 3: [t.sub.sat(MAX)] of the current-controlled fractional order memristor against [absolute value of [I.sub.DC]]: [alpha] = 0.75 (red), [alpha] = 1 (green), and [alpha] = 1.25 (blue).

Caption: Figure 4: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] = 110 [micro]A: [alpha] = 0.75 (red), [alpha] = 1 (green), [alpha] = 1.25 (blue), and HP memristor (black dots).

Caption: Figure 5: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] = -110 [micro]A: [alpha] = 0.75 (red), [alpha] = 1 (green), [alpha] = 1.25 (blue), and HP memristor (black dots).

Caption: Figure 6: M(t) of the boundary effect free fractional order memristor under the sinusoidal input with [I.sub.m] = 110 [micro]A and [alpha] = 1.25.

Caption: Figure 7: M(t) of the boundary effect free fractional order memristor under the sinusoidal input with [I.sub.m] = -110 [micro]A and [alpha] = 1.25.

Caption: Figure 8: v(t)-i(t) of the fractional order memristor: [alpha] = 0.75 and [I.sub.m] > 0.

Caption: Figure 9: v(t)-i(t) of the fractional order memristor (green) and HP memristor (black dots): [alpha] = 1 and [I.sub.m] > 0.

Caption: Figure 10: v(t)-i(t) of the fractional order memristor: [alpha] = 1.25 and [I.sub.m] > [degrees].

Caption: Figure 11: v(t)-i(t) of the fractional order memristor: [alpha] = 0.75 and [I.sub.m] < [degrees].

Caption: Figure 12: v(t)-i(t) of the fractional order memristor (green) and HP memristor (black dots): [alpha] = 1 and [I.sub.m] < 0.

Caption: Figure 13: v(t)-i(t) of the fractional order memristor: [alpha] = 1.25 and [I.sub.m] < [degrees].

Caption: Figure 14: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] > 0: [alpha] = 0.75.

Caption: Figure 15: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] > 0: [alpha] = 1.25.

Caption: Figure 16: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] < 0: [alpha] = 0.75.

Caption: Figure 17: M(t) of the fractional order memristor under the sinusoidal input with [I.sub.m] < 0: [alpha] = 1.25.

Caption: Figure 18: v(t)-i(t) of the fractional order memristor at asymptotic state ([I.sub.m] > 0): [alpha] = 0.75 (red), [alpha] = 1 (green), and [alpha] = 1.25 (blue).

Caption: Figure 19: v(t)-i(t) of the fractional order memristor at asymptotic state ([I.sub.m] < 0): [alpha] = 0.75 (red), [alpha] = 1 (green), and [alpha] = 1.25 (blue).

Caption: Figure 20: M(t)-i(t) of the fractional order memristor ([I.sub.m] > 0): [alpha] = 0.75 (red), [alpha] = 1 (green), and [alpha] = 1.25 (blue).

Caption: Figure 21: M(t)-i(t) of the fractional order memristor ([I.sub.m] < 0): [alpha] = 0.75 (red), [alpha] = 1 (green), and [alpha] = 1.25 (blue).

Caption: Figure 22: [M.sub.a](t) of the fractional order memristor under the sinusoidal input with [I.sub.m] > 0: [theta] = 45[degrees] (red), [theta] = 135[degrees] (green), [theta] = 225[degrees] (blue), and [theta] = 315[degrees] (magenta).

Caption: Figure 23: [M.sub.a](t) of the fractional order memristor under the sinusoidal input with [I.sub.m] < 0: [theta] = 45[degrees] (red), [theta] = 135[degrees] (green), [theta] = 225[degrees] (blue), and [theta] = 315[degrees] (magenta).

Caption: Figure 24: The memristor based Wien oscillator [29].

Caption: Figure 25: [V.sub.out](t) at asymptotic state of Type A Wien oscillator: fractional memristor with [alpha] = 0.75 (blue), fractional memristor with [alpha] = 1 (green), fractional memristor with [alpha] = 1.25 (red), and HP memristor (black dots).

Caption: Figure 26: [V.sub.out](t) versus M(t) of Type A Wien oscillator: fractional memristor with [alpha] = 0.75.

Caption: Figure 27: [V.sub.out] (t) versus M(t) of Type A Wien oscillator: fractional memristor with [alpha] = 1 (green line) and HP memristor (orange dashed).

Caption: Figure 28: [V.sub.out](t) versus M(t) of Type A Wien oscillator: fractional memristor with [alpha] = 1.25.

Table 1: [M.sub.TI] due to [I.sub.M] > 0. [alpha] [M.sub.TI] (k[ohm]) 0.75 6.9296 1.00 8.2175 1.25 9.7365 Table 2: [M.sub.TI] due to [I.sub.M] < 0. [alpha] [M.sub.TI] (k[ohm]) 0.75 24.4306 1.00 20.2443 1.25 16.4772 Table 3: [DELTA][f.sub.osc] (Hz) for various M(0)'s. M(0) (k[ohm]) [alpha] = 0.75 HP ([alpha] = 1) 4.1 9.16 3.32 4.4 7.64 2.88 4.7 6.44 2.52 5.0 5.49 2.23 5.3 4.73 1.98 5.6 4.11 1.78 5.9 3.59 1.60 M(0) (k[ohm]) [alpha] = 1.25 4.1 1.49 4.4 1.32 4.7 1.19 5.0 1.07 5.3 0.97 5.6 0.89 5.9 0.81

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Banchuin, Rawid |

Publication: | Active and Passive Electronic Components |

Article Type: | Report |

Date: | Jan 1, 2018 |

Words: | 8583 |

Previous Article: | A New Fractal Multiband Antenna for Wireless Power Transmission Applications. |

Next Article: | DRV Evaluation of 6T SRAM Cell Using Efficient Optimization Techniques. |

Topics: |