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

In this work, the analytical expressions of memristances, related parameters, and time domain behavioral analysis of the fractional order memristor have been proposed. Both DC with arbitrary delay and many AC waveforms including arbitrary phase sinusoidal and cosinusoidal waveform along with arbitrary periodic waveform have been taken into account. Unlike the previous works, the formerly ignored dimensional consistency has been taken into account and the analytical modelling of the boundary effect has been performed. Moreover, 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. The influence of the fractional order to the areas of voltage-current hysteresis loop and memristance-current lissajous curve has also been clearly discussed and the usage of fractional order memristor in the memristor based circuit has also been demonstrated.


Introduction
Recently, 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][5][6], and electronic engineering [7][8][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][12][13][14][15]).
For the fractional order memristor on the other hand, there also exists such modelling and analysis attempts [16][17][18][19][20][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][19][20][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][19][20][21], both transient and asymptotic behaviors of the fractional order memristor excited by AC 2 Active and Passive Electronic Components 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 voltagecurrent 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.

The Overview of Memristor
Memristor is a nonlinear electrical circuit element.This circuit element relates the instantaneous flux, (), and charge, (), through the following relationship: where () denotes the memristance.
According to [10], () of the HP memristor can be given in terms of the minimum and maximum values of M(t) denoted by   and   and the state variable, (), as where (), which is dimensionless, can be given in terms of the memristor's current, (), by Note that  =   / 2 , where  and , respectively, stand for the ion mobility and semiconductor film of thickness.Therefore, the dimension of  is (Asec) −1 .
As can be seen from ( 3) and also according to [10], () can be simply given as follows: Therefore, it can be seen that the HP memristor is charge controlled.Since () is a time integration of (), it can be stated that the HP memristor is of a current controlled type.Note also that 0 ≤ () ≤ 1; thus,   ≤ () ≤   as long as the memristor is unsaturated.Otherwise, () will be bounded at either 0 or 1 so () will be equal to either   or   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].

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][19][20][21], we perform such generalization by applying the fractional calculus to the memristor's state equation, i.e., (3).In these previous works,   /  = (), where  stands for the order of the fractional order memristor which can be arbitrary real value and has been obtained from such generalization.However, as () is dimensionless; the dimension of the LHS of this previous generalized equation is given by sec − where that of the RHS is sec −1 which means that a dimensional inconsistency has always existed.Therefore, the fractional time component [22], , which has the dimension of sec, has been introduced for handling this issue.As a result, unlike [18][19][20][21], the following generalized state equation has been used instead.
Similarly to that of the RHS, the dimension of the LHS of ( 5) is sec −1 ; thus the dimensional inconsistency issue has been resolved.Note also that ( 5) is reduced to (3) when  = 1 despite the presence of  as  −1 become 1 with such value of .
Unlike [18,19,21], we derive () of the fractional order memristor as a function () 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 () from (5) by using the Riemann-Liouville fractional order integral [25] as follows: where (0) and D() denote the initial value of () and the Gamma function [26], respectively.Since it can be seen from ( 2) that where   =   −   , the initial memristance value, i.e., (0), can be immediately given by Thus by substituting ( 6) into (7) and keeping (8) in mind, () of the fractional order memristor can be obtained as follows: which shows that () is current-controlled.
If we let  = 1, (9) will be reduced to Since the integer order integration of () gives (), we obtain By using ( 8), (11) can be simplified under the assumption that (0) = 0 and   ≫   as follows: 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(), yields ().
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 (()) = 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[, ] and min[, ], which, respectively, selects the maximum value and minimum value among  and , to (9).As a result, our expression of () due to arbitrary exciting waveform can be finally given as follows: where  =  1− ; thus the dimension of  is A −1 sec − .Since the dimension of fractional integral of () is Asec  , that of (  /Γ()) ∫  0 ( − ) −1 () is given by Ω which is physically measurable, similarly to those of   and (0).Therefore, (0) and (  /Γ()) ∫  0 ( − ) −1 (), which are at the RHS of (13), can be physically combined as they have the same dimensions and the dimension of (), which is the LHS of such equation, has also been found to be such physically measurable Ω; thus it can be seen that our expression of () has dimensional consistency.Moreover, due to the operation of nested max[, ] and min[, ], () will be equal to (0) − (  /Γ()) ∫  0 ( − ) [  ,   ] which means that the fractional order memristor remains unsaturated.Otherwise, () will be equal to either is lower than   or higher than   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, ()'s due to due to various exciting waveforms and the behavioral analysis of fractional order memristor will be presented.

DC Waveform.
Mathematically, the DC waveform with arbitrary delay (  ), which is more generic than the undelay waveform assumed in the previous works [18,19], can be defined as () =   ( −   ), where   and () denote the magnitude of the waveform and the unit step function.Therefore, the resulting () can be straightforwardly obtained by using (13) as By using (14) with   = 0 sec,  = 100000 A −1 sec − ,   = 1 kΩ,   = 100 kΩ, and (0) = 80 kΩ, ()'s of the fractional order memristor with various 's excited by the DC waveform can be numerically simulated as depicted in Figures 1 and 2 where   = 110 A and   = -110 A have been, respectively, assumed and ()'s of the HP memristor simulated by using its SPICE model [27] have also been included.
From these figures, the strong agreements between our ()'s obtained by using (14) with  = 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 () of the fractional order memristor can be either the increasing or decreasing function of  with the final value of    or   when the fractional order memristor become saturated if we let   > 0 or   < 0.Moreover, the rate of change of () has been found to be inversely proportional to .
Since it can be seen that the fractional order memristor become saturated at a certain time given by   , by using ( 14), we have As a result,   can be immediately given as follows: which shows that   is directly proportional to the size of the difference between (0) and   (  ).So,   reaches its maximum value, i.e.,  () , if and only if (0) reaches its possible peak value given by either   when   > 0 or   when   < 0. Thus  () can be found as which can be immediately given in a more compact manner as follows: For confirming the aforesaid observation on the relationship between the rate of change of () and ,  () 's have been simulated by using ( 18) under the similar assumptions to those of the simulation of ()'s shown in Figures 1 and 2 but with varying   , as depicted in Figure 3 which shows that  () is directly proportional to .Due to the definition of t () , this confirms such observation.Moreover, it can also be seen from Figure 3 that  () is inversely proportional to |  |.

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 () has been chosen due to its generality.Mathematically, such waveform can be given by () =   sin( + ), where   and , respectively, denote the peak value and angular frequency of ().By using (13), () of the fractional order memristor under the arbitrary phase sinusoidal input can be given by where 1 F 2 ( ; , ; ) denotes a generalized hypergeometric function with  = 1 and  = 2 [28].
By letting  = 100000 A −1 sec − ,   = 1 kΩ,   = 100 kΩ, and (0) = 80 kΩ similarly to the previous subsection, ()'s of the fractional order memristor with various 's excited by the sinusoidal input can be simulated as depicted in Figures 4 and 5, where   = 110A and   = -110 A have been, respectively, assumed and ()'s of the HP memristor simulated by using its SPICE model have also been included.Moreover, we also let  = 1 rad/sec and  = 0 rad.Again, the strong agreements between our expression based ()'s with  = 1 and those of the HP memristor which verify our expression can be observed.It can be seen that () of the fractional order memristor with  ≤ 1 is periodic with clipped peaks at   and   due to the temporary saturation of the fractional order memristor as the unclipped peaks of () lie outside [  ,   ].On the other hand, () of the device with  > 1 is not periodic but time independently equal to either   or   when  ≥   up to the sign of   .This is because the device becomes permanently saturated after   as () contains the time proportional term which starts to lie outside [  ,   ] at  =   when  > 1.In order to show that the permanent saturation of the fractional order memristor with  > 1 excited by sinusoidal input is possible, we simulate () of the device with  = 1.25 once again by using (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 () be outside [  ,   ] such as those of Figures 4  and 5. Moreover the determination of such time proportional term and   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 <  ≤ 19 sec,  = 1 rad/sec, and  = 0 rad; the lissajous curves of the fractional order memristor with various 's can be simulated as depicted in Figures 8-13 where   = 110 A, K = 10000 A −1 sec − ,   = 100 Ω, M  = 16 kΩ, and M(0) = 11 kΩ have been assumed in Figures 8-10.On the other hand,   = -150 A,  = 10000 A −1 sec − ,   = 100 Ω,   = 38 kΩ, and M(0) = 11.2 kΩ 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  = 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  ̸ = 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  = 1 and HP memristor, the shape of the lissajous curve of the fractional order memristor with  ̸ = 1 keeps changing.This is because () of such fractional order memristor contains the time proportional term.This can be clearly seen from Figures 14-17 which display ()'s of the fractional order memristor with  ̸ = 1.Unlike the previous ones depicted in Figures 4 and 5, the peaks of these ()'s with  < 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 [  ,   ].Moreover, the fractional order memristor with  > 1 never becomes saturated and can be seen from Figures 15 and 17.This is because the time proportional terms of its ()'s are always be within [  ,   ].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 () becomes time independent instead.Therefore, the shape of its lissajous curves of the fractional order memristor with  ̸ = 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 's simulated by assuming that 100 sec <  ≤ 119 sec.Apart from being asymmetric when  ̸ = 1, we have found that the lobe area of the lissajous curve of the fractional order memristor is affected by  similarly to those proposed in [20,21].In particular, we have found that the fractional order memristor with lower  yields the lissajous curve with wider lobe area which refers to more memory effect and less linearity.Besides the voltage-current curve, the memristancecurrent lissajous curve has been found to be also interesting.Therefore, such curves of the fractional order memristor with various '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 <  ≤ 119 sec.These figures show that the resulting lissajous curves take the elliptical closed loop shape which are unchanged as () enters the asymptotic state at the assumed time interval.The elliptical shaped lissajous curve means that () is periodic as well as

𝑖(𝑡).
Similarly to the voltage-current curve, the loop area of the memristance-current lissajous curve is also affected by .In particular, the fractional order memristor with lower  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 .Here, we let the area of the upper and lower lobes of the lissajous curve be denoted, respectively, by   and   .Since V() =   ()() at asymptotic state where   () denotes the asymptotic approximation of (),   and   can be given by With ( 13) and the asymptotic approximation of the fractional order integration of sinusoidal function [20],   () can be found as As () =   sin(+) and the fractional order memristor is unsaturated for the entire simulation period as can be seen from Figures 18 and 19,   and   can be finally obtained as follows: which show that the upper and lower lobes of the lissajous curve employ equal sizes of areas and such areas are independent of .Unlike [20] which considered only   ,   has also been formulated in this work.Moreover, our   which has been derived by using a different approach from that used in [20] totally agrees with the previous result that  = 0 rad has been assumed.This emphasizes the independency from  of the lobe areas.Moreover, ( 24) and ( 25) also show that the fractional order memristor does not store the energy as the summation of   and   which, respectively, referred to the intake and dissipated energy, is equal to 0.
Besides   and   , the area within the closed loop of the memristance-current lissajous curve (  ) can be obtained by also using (23) and Since it can be seen from Figures 20 and 21 that the device is unsaturated for the whole simulation period and () =   sin( + ), we have which shows that   is affected by  but independent of  as well as   and   .
Despite the aforementioned independencies,  does affect the behavior of the fractional order memristor.For illustration, M  (t)'s due to various 's have been simulated by using (20) as depicted in Figures 22 and 23 where  = 1 and  = 1 rad/sec has been assumed.It should be mentioned here that kΩ, and M(0) = 11 kΩ have been adopted in Figure 20 where I  = -150 A, K = 10000 A −1 sec − , M  = 100 Ω, M  = 38 kΩ, and M(0) = 11.2 kΩ have been assumed in Figure 21.Moreover,  = /4 rad,  = 3/4 rad,  = -3/4 rad, and  = -/4 rad have been chosen as they are good representatives of those 's which their coordinates on the Euclidian plane, i.e., (cos(), sin()), 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(), sin())'s of these chosen 's are exactly located at the middle points of the portion of unit circle's arc.For example, (cos(), sin()) of  = /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   ()'s contain time independent terms which are formerly the time proportional term of () that become time independent at asymptotic state as aforementioned.Since these time independent terms lie within [  ,   ],   () takes the shape of sinusoidal waveform.However, the minimum peaks of   ()'s due to the input with   < 0 and (cos(), sin()) 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   .From Figure 22, it can be seen that the input with   > 0 and  with (cos(), sin()) located on the portion of unit circle's arc in quadrants 2 and 3 yields () with higher time average.If we have assumed that   ≤   () ≤   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   < 0 and (cos(), sin()) located on such portion in quadrants 1 and 4. Now, we will derive the analytical expression of time independent term of   ().Let such term be denoted by   ; it can be given by carefully observing (23) as follows: When  = 1 as assumed in our simulations of   (),   can be reduced to where Therefore, it can be seen that when   > 0 (  < 0), if (cos(), sin()) is located on the portion of unit circle's arc in either quadrant 2 or 3 where cos() < 0, Δ  > 0 (Δ  < 0) thus (0) has been increased (decreased).On the other hand, if (cos(), sin()) is on the portion of unit circle's arc in either quadrant 1 or 4 where cos() > 0, Δ  < 0 (Δ  > 0) as   > 0 (  < 0) thus (0) has been decreased (increased).As a result, the sinusoidal input with   > 0 (  < 0) and (cos(), sin()) located on the portion of unit circle's arc in either quadrant 2 or 3 (1 or 4) of the Euclidian plane yields higher   as graphically observed.
By subtracting   from   (), we obtain the following purely periodic term of (),   () which existed in both initial state and asymptotic state.Therefore, the time proportional term of (),   () can be formulated by subtracting   () from () as As a result, the aforementioned   can be determined by solving the following equation:  At this point, we will derive the expressions of ()'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 () due to arbitrary phase cosinusoidal waveform, i.e., () =   cos( + ), can be formulated by using (19) and the relationship between the sinusoidal and cosinusoidal functions is given by cos() = sin((/2) − ) where  denotes arbitrary angle.As a result, we have 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., () = ∑ (35)

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  1 or  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  1 of the circuit by such fractional order device instead of the HP memristor as traditionally did [29] and analyze the effect of  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 Δ  ; it can be mathematically defined as given by (36) where   and   stand for the upper bound and lower boundary.
According to [29], these boundaries can be obtained by solving (37) and (38) which have been formulated by assuming that  1 =  2 = , 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   =   /(0), where   denotes the peak value of the voltage dropped across the memristor which can be determined from the oscillator output voltage,   (), and also depends on the initial voltages of  1 and  2 [29].Moreover,   () 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  1 () and  2 () denote the voltage dropped across  1 and  2 and () stands for the memristance of the conventional HP memristor.

𝑓
However, this is not the case when the fractional order memristor has been used as () of such device must be adopted.Therefore   () must be determined based on our derived () instead where   and   must be evaluated from ( 41) and (42).As a result, the corresponding Δ  will be different from that of the original conventional memristor based circuit and the different chance of obtaining sustained oscillation can be expected.
Moreover, the condition for ensuring the occurrence of sustained oscillation can be given by () where  V which stands for the average oscillating frequency is given by ().
By assuming 0.1 V and -0.95 V as the approximate initial voltages of  1 and  2 and also assuming that  2 = 5 kΩ,  1 =  2 = 3.2F,  = 100000 A −1 sec − ,   = 100 Ω, and   = 16 kΩ, Δ  of the fractional order memristor based Wien oscillator with  < 1 and  > 1 can be obtained for various (0)'s by numerically solving (41) and (42).as tabulated in Table 3 where Δ  of the original HP memristor based circuit which is equivalent to the fractional order memristor with  = 1 in the context of this work, determined by solving (37) and (38), has also been included.It has been found that Δ  is inversely proportional to (0) which is in agreement with [29].Since 4.1 kΩ ≤ (0) ≤ 5.9 kΩ [29] and the probability of obtaining sustained oscillation is inversely proportional to Δ  , (0) = 5.9 kΩ is recommended as it minimizes Δ  and thus maximizes such probability.It can also be seen that Δ  is inversely proportional to  which means that the probability of obtaining sustained oscillation is directly proportional to .Therefore, the fractional order memristor with larger  is recommended and the fractional order memristor with  > 1 should be adopted as it can increase such probability from that of the original circuit which its memristive device has  = 1.On the other hand, the fractional order memristor with  < 1 should be avoided as it decreases such probability.
Finally, by further assuming that (0) = 5 kΩ,  3 = 20.2kΩ, and  4 = 10 kΩ, we can simulate   () at asymptotic state and the lissajous patterns of   () and () for various 's as depicted in Figures 23-28 where the units of 100   (), , and (t) are V, sec, and Ω, respectively.The comparison of the results with  = 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 () is always within [  ,   ] 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   () is of the same frequency as () and the phase difference between   () and (), which is less than 90 ∘ , is inversely proportional to .

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][19][20][21], the dimensional consistency has been taken into account.Moreover, the boundary effect has also been modelled.Therefore the analytical expression of () which has been derived as a function of () and its related results are dimensional consistent and valid to such generalized device in its saturation states.By using such expression, ()'s due to various waveforms and the related parameters have been formulated.With the simulations by using these ()'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].

Figure 7 :
Figure 7: () of the boundary effect free fractional order memristor under the sinusoidal input with   = -110 A and  = 1.25.
(19)h is(19)without the boundary effect.The simulation results have been depicted in Figures6 and 7which show that ()'s eventually lie outside [  ,   ] permanently.This yields the saturation when such effect has been included.() of the boundary effect free fractional order memristor under the sinusoidal input with   = 110 A and  = 1.25. −