On the fractional domain generalization of memristive parametric oscillators

Abstract In this research, we generalize a family of electronic parametric oscillators in the fractional domain by using a state of the art circuit element namely fractional memristor. Such family of parametric oscillators is the memristor based Wien family which is an extension of the normal Wien family. Noted that such normal Wien family is one family of the simplest second-order nonparametric oscillators. We derive the equations of the range of oscillating frequency, sustained oscillating frequency, sustained oscillating condition and the output voltage by using our mathematical model of the fractional memristor as the basis. With the obtained results and numerical simulations, the effects of the fractional memristor to the generalized parametric oscillators have been studied where the validation has been performed based on the SPICE HP memristor model. We have found that those oscillators with the fractional memristor of order greater than unity are more preferable.


ABOUT THE AUTHOR
Rawid Banchuin received the B.Eng. degree in electrical engineering from Mahidol University, Bangkok, Thailand in 2000, the degree of M.Eng. in computer engineering and Ph.D. in electrical and computer engineering from King Mongkut's University of Technology Thonburi, Bangkok, Thailand in 2003 and 2008 respectively.
He was with the Department of Electrical Engineering, Rajamangala Institute of Technology, North-Bangkok Campus since 2003.At the present, he is with the Graduated School of Information Technology and Faculty of Engineering, Siam University, Bangkok, Thailand.
Asst.Prof. Dr. Rawid Banchuin is a member of council of engineer (Thailand) and has joined the organizing committee of the international conference on ICT and knowledge engineering, which is jointly organized by IEEE, since 2012.His current research interests include circuit elements with memory, fractance device and fractional order circuits/systems.

PUBLIC INTEREST STATEMENT
The electronic oscillators have been found to be the famous signal generators.Among various oscillators, the Wien oscillator is one of the simplest electronic oscillators.By replacing one of its resistor with the memristor, the oscillator become a parametric one thus it has unique properties which are different from those of the normal oscillator, for example, time-dependent pole locations, oscillating frequency, etc.
On the other hand, a state of the art electronic device namely fractional memristor have been developed.Compared to this recently developed device, the conventional memristor is merely its special case.Motivated by such far more generality of the fractional memristor and the interesting properties of parametric Wien oscillator, we generalize the entire family of this oscillator by applying such newly developed device.

Introduction
Mathematically, the oscillation refers to the continuous recurrence of periodic waveform.The oscillation is parametric if the oscillatory system under consideration is time varying.Many examples of parametric oscillation have been previously presented by Komine et al. (Komine, Galliou, & Makarov, 2003).In electronic engineering, those circuits which generate oscillatory outputs called oscillators have been found to be the famous signal generators.Among various oscillators, Wien oscillator has been found to be one of the simplest second-order oscillators.Such Wien oscillator has four basic structures, which constitute the Wien family.By replacing one of the resistor in the Wien oscillator by the memristor, the memristor-based Wien oscillator can be obtained (Talukdar, Radwan, & Salama, 2011).This memristor-based oscillator has been found to be an electronic parametric oscillatory system because the coefficients of its state space representation are time dependent thus it is a linear time varying system.Moreover, it has been shown that this parametric oscillator has unique properties, which are different from those of a nonparametric Wien oscillator, for example, time-dependent pole locations, oscillating frequency, etc. (Talukdar et al., 2011).
The memristor, which is the key element of this oscillator, is the fourth electrical circuit element that has been theoretically found by Chua since 1971(Chua, 1971).Since the memristor is simply the resistor with memory, it can be thought of as a generalization of resistor.Therefore the memristor-based Wien oscillator has been found to be a generalization of the normal Wien oscillator.For decades after Chua proposed his original work, the memristor has been practically realized by a research group in Hewlett Packard (HP) labs in 2008 (Strukov, Snider, Stewart, & Williams, 2008).Apart from the aforesaid parametric oscillator, the memristor has also been applied in various applications, for example, DRAM, signal processing, neural networks, programmable logic, cross bar switch, passive switch, control systems etc. (Pershin & Di Ventra, 2010), (Mouttet, 2007), (Jeong, Lee, Choi, & Kim, 2009).
Recently, a state of the art electrical circuit element namely fractional memristor have been existed.This circuit element can be obtained from a further generalization of the memristor in the fractional domain by using the concept of fractional calculus, which have been adopted in various disciplines, for example, biomedical engineering (Jesus, Machado & Cunha, 2008), (Tang et al., 2009), control system (Charef, 2006), (Vinagre & Feliu, 2007), (Matušů, 2011) and electronic engineering (Dorčák, Terpák, Petráš, & Dorčáková, 2007), (Krishna, Reddy, & Santa Kumari, 2008), (Radwan & Elwakil, 2012).It can be seen that the memristor is merely a special case of the fractional memristor as it is simply the fractional memristor of first order.Motivated by the far more generality of the fractional memristor and the interesting properties of memristor based Wien oscillators, we generalize this family of parametric oscillators in the fractional domain by applying the fractional memristor.We firstly derive the mathematical model of the fractional memristor for being the basis of our work.Later, we derive the equations of the range of oscillating frequency, sustained oscillating frequency and sustained oscillating condition, which are the important features of the parametric Wien family.We also solve the state space equations of these fractional memristor-based parametric oscillators for determining the output voltage.By using these results and numerical simulations with MATHEMATICA, the effects of the order of fractional memristor to the aforesaid important features and the output voltage have been studied.We also compare our simulation results to those obtained by using the circuit simulations with SPICE model of HP memristor proposed by Biolek et al. (Biolek, Biolek & Biolkova, 2009a, August), (Biolek, Biolek & Biolkova, 2009b), (Biolek, Biolek, Biolkova, & Kolka, 2014).As a result, the strong agreements can be observed.Moreover, we have found that the generalized parametric oscillators with the fractional memristor of order greater than unity are more preferable.
In the following section, some background theory on the memristor and Wien family of oscillators will be briefly given.The mathematical modelling of fractional memristor will be shown in section 3 followed by the derivation of the aforesaid equations in section 4. In section 5, the effects of the order of fractional memristor to the aforementioned important features and output voltage will be studied where some analytical discussions and discussions on other interesting issue will be given in section 6.Finally, the conclusion will be drawn in section 7.

Memristor
Memristor is a nonlinear electrical circuit element.This circuit element relates the instantaneous flux (ϕ(t)) and charge (q(t)) through the following relationship where M(t) denotes the memristance which is time variant.
According to (Strukov et al., 2008), M(t) of the HP memristor can be given in terms of the minimum and maximum values of M(t) denoted by M on and M off and the state variable (x(t)) as where x(t) which is dimensionless, can be given in term of the memristor's current (i(t)) by Noted that k = μM on /D 2 where μ and D respectively stand for the ion mobility and semiconductor film of thickness.Therefore the dimension of k is (A•s) −1 .Moreover, f ðxðtÞÞ stands for the window function, which has been used for modelling the boundary effect of the device, that is, 0

Wien family of oscillators
Wien family of oscillators or Wien family in short are a family of the simplest second-order oscillators.There are four members in this family where each of them has a unique circuital structure.So, they are referred to as the Wien oscillator of type A, B, C and D which can be respectively depicted in Figures 1-4.
Despite their different structures, the oscillating condition and frequency of oscillation (f osc ) of these members of Wien family can be commonly given by ( 4) and ( 5) respectively where C = C 1 = C 2 .Such condition and f osc are time independent as R 1 , R 2 , C 1 and C 2 do.
For the memristor based Wien family, R 1 has been replaced by a memristor (Talukdar et al., 2011).Therefore, the oscillating condition and f osc become As a result, both oscillating condition and f osc are now time variant thus the parametric oscillation can be achieved by the usage of memristor.

The mathematical model of fractional memristor
By generalizing the memristor in the fractional domain, the fractional memristor can be obtained.Such generalization can be performed by applying the fractional calculus concept to the memristor's state equation.According to the previous fractional memristor modelling attempts (Fouda & Radwan, 2013), (Fouda & Radwan, 2015), such fractional domain state equation can be given by where α stands for the order of the fractional derivative which is not strictly integer but arbitrary real value.
However, as x(t) is dimensionless, the dimension of the LHS of this previous generalized is given by sec −α where that of the RHS is sec −1 , which means that a dimensional inconsistency has been existed.Moreover, the boundary effect has also been ignored as f(x(t)) has been excluded.Therefore, the fractional time component, σ which has the dimension of s (Gómez-Aguilar, Rosales-García, Bernal-Alvarado, Córdova-Fraga, & Guzmán-Cabrera, 2012), (Podlubny, 2002), (Banchuin, 2017), must be introduced and (3) must be used as the basis for handling these issues.As a result, the following fractional domain state equation which has been adopted in this study, can be obtained Similarly to that of the RHS, the dimension of the LHS of ( 7) is sec −1 thus the dimensional inconsistency issue has been resolved.Moreover, the boundary effect has been now taken into account by the inclusion of f ðxðtÞÞ.Noted also that (9) can be reduced to (3) when α = 1 despite the presence of σ as σ α −1 become 1 with such value of α.
Since the linear dopant drift has been assumed by Talukdar et al. according to their adopted mathematical model of the basis memristor (Radwan, Zidan & Salama, 2010b), (Radwan, Zidan & Salama, 2010a, August), we use the following f ðxðtÞÞ in our work for ceteris paribus.By using the Riemann-Liouville fractional order integral (Sabatier, Agrawal, & Machado, 2007), x(t) can be found as where x(0) denotes the initial value of x(t) and Г() stands for the Gamma function (Beyer, 2018), Moreover, we define K = kσ 1-α for simplicity thus the dimensions of Since it can be seen from ( 2) that where M d = M off -M on , the initial memristance value i.e.M(0), can be immediately defined as Therefore M(t) of the fractional memristor can be obtained by using ( 2), ( 10) and ( 12) as remained Ω similarly to those of the traditional resistance and memristance of the conventional memristor.
By using (15), M(t) of the fractional memristor due to any exciting signal can be determined.If we assume that such signal is a sinusoidal one with I m , f and θ as its arbitrary peak value, frequency and phase angle, M(t) can be obtained by applying the approximation of the fractional order integration of sinusoidal function (Kilbas, Srivastava, & Trujillo, 2006) as follows

The fractional memristor based Wien family
This generalized parametric Wien family can be obtained by using the fractional memristor for R 1 in the Wien oscillator of all types.As a result, the state space representations of all members of the family can be commonly given by where Noted also that V C1 (t) and V C2 (t) denote the voltage dropped across C 1 and C 2 and where M(t) is given by ( 16) as the fractional memristor is applied by the sinusoidal signal because the Wien oscillator is a sinusoidal oscillator.
Since M(t) is in terms of trigonometric function, it exhibits an oscillatory behavior and so does f osc as f osc is depended on M(t).Similarly to its conventional memristor based prototype, f osc of the fractional memristor based Wien family is also related to M(t) by ( 7) thus it is inversely proportional to M(t).As a result, f osc reaches its maximum when M(t) reaches its minimum and vice versa.By using (7) with M(t) as given by ( 16), it has been found that such maximum and minimum value of f osc , that is, f max and f min , can be respectively determined by solving (17) and ( 18) which are applicable if and only if I m > 0.
For I m < 0 on the other hand, ( 19) and ( 20) must be applied instead.
Before we proceed further, it should be mentioned here that ( 17)-( 20) and the ongoing ones have been derived by assuming that M min < MðtÞ < M max has been satisfied in order to ensure the proper operations of the oscillators as previously did by Talukdar et al.At this point, it can be stated that the range of f osc (Δf osc ) can be simply obtained in terms of f max and f min as follows Similarly to its conventional memristor-based progenitor, the fractional memristor-based Wien family finally reach the sustained oscillation state despite the fluctuation in f osc .At such state, f osc become constant at its average value (f avr ) where we have found that Moreover, the following condition must be satisfied for ensuring the occurrence of sustained oscillation.
In order to determine the output voltage (V out (t)), ( 17) must be firstly solved.The solution of (17) can be found as where Noted also that I is a 2 dimensional identity matrix.
Despite the similar V C1 (t) and V C2 (t) as the state space representations of all members of the family can be commonly given by ( 17), each member of the fractional memristor-based Wien family employs different V out (t) due to its unique circuit topology.For type A and type B fractional memristor-based Wien oscillator, V out (t) can be given by either ( 26) or ( 27) where that of type C and type D can be given by either ( 28) or (29).

The fractional domain behavioral analysis
In order to do so, the effects of α which is the fractional domain parameter, to f avr , sustained oscillating condition, Δf osc and V out (t) which are the key features of this parametric Wien family, must be analyzed by using those equations derived in the previous section and numerical simulations.Noted also that all members of the family have been studied and the analytical discussion on the obtained results have also been made as will be seen in the subsequent section.This is unlike the previous analysis attempt which only type A oscillator has been analyzed where merely some features has been considered and θ = 0 rad has been assumed even in the mathematical formulations for simplicity (Banchuin, 2018).Throughout this work, we have assumed that C 1 = C 2 = C = 3.2 μF and R 2 = 5 kΩ for all members of the family (Talukdar et al., 2011).However, we use R 4 = 10 kΩ unlike Talukdar et al.We have also assumed that D = 10 nm, σ = 1 s, M off = 16 kΩ, M on = 100 Ω and μ = 10 −14 m 2 (s•V) −1 .For analyzing the effects of α, we let it take three possible values, that is, 0.75 which is lower than 1, 1 which makes the fractional memristor be equivalent to an conventional one and 1.25 which is greater than 1.As results, f avr , R 3 /R 4 for ensuring sustained oscillation and Δf osc have been tabulated against α and M(0) for all types of fractional memristor-based Wien oscillator as shown in Tables 1-9.Apart from being inversely proportional to M(0) similarly to that of the memristor-based Wien family, f avr is directly proportional to α for the fractional memristor based oscillators of type A, C and D as can be seen from Table 1 and III.On the other hand, it can be seen from Table 2 that the opposite relationship can be observed from type B circuit.From Table 4 and VI, it can be observed that R 3 /R 4 for sustained oscillation guaranteed is directly proportional to M(0) but inversely proportional to α for the oscillators of type A, C and D. On the other hand, such R 3 /R 4 become directly proportional to α for type B which yields lower values of R 3 /R 4 than those of the other types when α ≤ 1 as can be seen from Table 5.Therefore, type B oscillator requires smaller R 3 than those of the others for similarly obtaining sustained oscillation.It can also be seen that the effects of α to f avr and R 3 /R 4 for ensuring sustained oscillation are small due to their insignificant dependencies of both quantities on α.This is not surprising because we have found in this scenario that Therefore, both f avr and R 3 /R 4 are largely depended on M(0), which is independent of α, as can be seen from ( 20) and ( 21).
However, we have found that α significantly affects Δf osc as Tables 7-9 shows that Δf osc is inversely proportional to α at significant degree beside being inversely proportional to M(0).Therefore, we promote the usage of fractional memristor with α > 1 for increasing the chance of obtaining sustained oscillation according to its lowest Δf osc .These tables also show that fractional memristor-based oscillator of type A employs minimum Δf osc thus it has been found to be more preferable than the others due to its best chance of obtaining the sustained oscillation.Moreover, we simulate V out (t) for various α's as depicted in Figures 5-10 where we assume that M(0) = 5 kΩ.In Figures 5-7, 0 s < t < 1 s which refers to the early state of V out (t) have been assumed.On the other hand, 99 s < t < 100 s which refers to the asymptotic state, have been assumed for the rests.
For the validation of our results, the comparison of our results with α = 1 to their corresponding SPICE HP memristor model-based benchmarks have been made as also shown in the figures.The schematic diagram of SPICE HP memristor model can be shown as depicted in Figure 11 (Biolek, Biolek & Biolkova, 2009b), (Biolek et al., 2014), (Kozma, Pino, & Pazienza, 2012) where TE and BE represent the top and bottom electrodes of the memristor (Kozma et al., 2012), i(t) denotes the current flowing through the memristor and x(t), which is physically the normalized width of the doped layer is now emulated in term of the voltage dropped across the parallel combination of conventional capacitor and resistor thus x(0) can be emulated as the capacitor's initial voltage.Those results with α = 1 have been chosen as the fractional memristor with such order is equivalent to the conventional one.Moreover, the circuit simulation by using the SPICE HP memristor model has been found to be convenient and also more preferable than using other HP memristor emulator circuits, for example, those proposed by Kim et al. (Kim, Sah, Yang, Cho, & Chua, 2012), Yang et al. (Yang, Cho, Sah, Kim, & Jung, 2012), Sah et al. (Sah, Yang, Kim, & Chua, 2012), Olumodeji (Olumodeji & Gottardi, 2017), etc.This is because the HP memristor which has been adopted in this work as our basis of the fractional memristor, is a specifically fabricated electronic device at nanometer level that employs a unique characteristic namely the boundary effect which makes the memristance be fixed either at on-state value or off-state value when become saturated (Anusudha & Prabaharan, 2018).This effect has been taken into account by the SPICE HP memristor model despite being ignored by those other HP memristor emulator circuits.Therefore, the SPICE HP memristor model has been found to be one of the best so far HP memristor emulator for computer simulation to the knowledge of the author.From Figures 5-10   the strong agreements between our results and the SPICE HP memristor model-based counterparts which verify the validity of our results, can be observed.Noted also that our asymptotic V out (t) of type A oscillator depicted in Figure 9 bears a strong resemblance to the previous result (Banchuin, 2018) as we let θ = 0 rad for our numerical simulation and different outcomes can be expected if we assume nonzero θ. 5-7 to Figures 8-10, we have found that the magnitude of V out (t) become significantly greater at asymptotic state which sustained oscillation has occurred, and the increasing in such magnitude is inversely proportional to α.This promotes the usage of the fractional memristor with α > 1 for enhancing the constancy in magnitude of the oscillator's output.It can also be observed that the magnitudes of V out (t) of the oscillators of type C and D are significantly larger than to those of the others.This is because R 3 > M(t) for all t according to the assumed parameters stated above.Therefore the magnitudes of V out (t) of type C and D which are directly proportional to R 3 as can be seen from ( 29), are significantly larger than to those of types A and B, which in turn are proportional to M(t) as shown by ( 27).Before proceed further, it should be mentioned here that there exist phase differences between V out (t)'s for different α's at the asymptotic state.These phase differences have been caused by the accumulation of the differences in f osc of V out (t)'s for different α's occurred during the early state as the phase angle is mathematically the integration of frequency.Such differences in f osc occurred because the sustained oscillation is not yet achieved during the early state thus f osc is still fluctuated.Since α significantly affects Δf osc , it affects the fluctuation of f osc therefore different α's causes different fluctuations and the differences in f osc thus the phase differences have been yielded.Now, we simulate V out (t) at asymptotic state against M(t) for various α's as depicted in Figures 12-20, where M(0) = 5 kΩ has also been assumed and the vertical axes have been multiplied by 100 for visibility.Again, the comparison of those results with α = 1 to their         corresponding SPICE HP memristor model-based benchmarks have also been made where strong agreements which validate our results, can also be observed.Noted also that the simulation results of type A oscillator depicted in Figures 12-14 are very similar to the previous results (Banchuin, 2018) as we let θ = 0 rad.From these lissajous curves, we have found that V out (t) is of the same frequency as M(t).Such frequency is given by f avr as the sustained oscillation has been already achieved.It has also been found that the phase difference between V out (t) and M(t) which is less than π/2 rad, is inversely proportional to α for the type A, C and D oscillators.For type B oscillator on the other hand, we have found that such phase difference which lies between π/2 rad to π rad, is directly proportional to α.It can also be seen that the memory effects and the linearity in the relationship between M(t) and V out (t) are inversely and directly proportional to α respectively where the oscillators of type A and B have been found to be more linear than the rests assuming the same α.Therefore, we promote the usage of either type A or type B oscillator and the fractional memristor with α > 1 in those applications which linearity issues are crucial.

By comparing Figures
Before proceed further, it is worthy to mention here that it is not surprising that the simulation results of type C and D oscillator are very much similar so that they are coincided as can be seen from the above tables and figures in spite of their different circuit structures.This is because such similarities can also be observed from the conventional memristor based type C and D oscillator as proposed by Talukdar et al.Therefore, it can be stated that the fractional domain generalization inherits the similarities in the results of types C and D memristor-based Wien oscillators.

Discussions
At this point, some further discussions on the obtained simulated results should be given.In order to do so, the approximated closed form expression of V out (t) should be derived.By ignoring the high order terms of (24) and using the linear approximations of 1/M(t) and its integral, V out (t) of for A and B fractional memristor-based Wien oscillator can be approximately determined by baring (26) in mind as given by (31).On the other hand, that of type C and D counterparts can be obtained with (28) be kept in mind as given by (32).Noted also that X(t) and Y(t) can be respectively given by ( 33) and (34) and f osc become f avr at the asymptotic state.
From ( 31) to (34), it can be seen that V out (t) contains sinusoidal term with time increasing magnitude.This explains the asymptotically larger magnitude of V out (t) as graphically displayed above.By carefully observing these equations, such sinusoidal term which is symbolized by Z(t), can be given as follows Therefore the aforesaid time increasing magnitude of Z(t) (A Z (t)) can be immediately given by As a result, we have the following rate of change of A Z (t) From (37), it has been found that Since the cosine function is non monotonic, we have found that dA Z (t)/dt is inversely proportional to α due to (38).This explains the above graphically portrayed relationship between the increasing in magnitude of V out (t) and α.
Moreover, as f avr is robust to α as already demonstrated, so does the time at which sustained oscillation is occurred (T sus ).This is because T sus can be determined given f avr and M(t) at t = T sus , that is, M(T sus ), by using (40) therefore any α related behavior of f avr is immediately inherited by T sus including the robustness to α. Noted also that (40) has been formulated by using ( 14) with M min < MðtÞ < M max satisfied and keeping in mind that f osc become steady at f avr at t = T sus as the sustained oscillation is occurred Now some other interesting issues will also be discussed.Firstly, if more degree of freedom is required then we can replace C 1 and C 2 by the fractional capacitors, that is, C 1β and C 2γ .Noted that the usage of fractional capacitors in the nonparametric Wien oscillators has been previously presented by Radwan et al. (Radwan, Soliman, & Elwakil, 2008).Since we let these fractional capacitors be of incommensurate orders for obtaining maximum degree of freedom, ( 22) and ( 26)-( 29) become ( 41)-( 45) where V C1β ðtÞ and V C 2γ ðtÞ denote the voltage dropped across C 1β and C 2γ .
From these equations, it can be seen that more degree of freedom has been obtained as V out (t) can also be now controlled via β and γ apart from α.As an illustration, we will demonstrate the controlling of V out (t) of type A oscillator by simply changing β and γ.Such type A oscillator has been chosen due to its higher linearity mentioned above.By using the above fractional capacitor based equations with M(0) = 5 kΩ, C 1β = 3.2 μF•s β−1 and C 2γ = 3.2 μF•s γ−1 , the resulting V out (t) can be simulated for various combinations of β and γ as depicted in Figures 21-23 where α = 1.25 has been assumed as we promote the usage of the fractional memristor with α > 1 due to its many pros.The significant different transient behaviors and steady state magnitudes of V out (t) can be observed in spite of their same α.This demonstrates the controllability of V out (t) by using β and γ.
For validation, the comparison with benchmarks simulated by using the SPICE HP memristor model with order of 1.25 have been made as also shown in these figures where the strong agreements, which verify such controllability and also validates our results with other α apart from 1 that has been verified in the previous section, can be observed.Noted that a similar controllability can also be achieved by the oscillator of other types.It should be mentioned here that the order of SPICE HP memristor model which is conventionally 1, can be made nonunity for emulating the fractional memristor by simply let the value of current controlled current source in its schematic be k α i(t)f(x(t)) and replacing the conventional capacitor with the fractional one where such fractional capacitor along with C 1β and C 2γ can be simply emulated by using the resistive-conventional capacitor network for circuit simulation purpose.This is because the fractional impedance function can be approximately given in term of the conventional rational polynomial impedance function (Matsuda & Fujii, 1993), (Oustaloup, Levron, Mathieu, & Nanot, 2000), (Valsa & Vlach, 2013).However, order of the resulting resistive-capacitor network can be very high for obtaining acceptably low phase ripple within the acceptably wide frequency range (Valsa & Vlach, 2013), which is necessary for accurately emulating the fractional capacitor as its phase is ideally constant with respected to the frequency.
Secondly, if nonlinearity near boundaries of the memristive device must also be taken into account, we must assume the nonlinear dopant drift thus another f(x(t)) beside ( 8) must be applied.The examples of such f(x(t)) are those proposed by Joklecar and Wolf (Joglekar & Wolf, 2009), Biolek et al. (Biolek, Biolek, & Biolkova, 2009b), Prodomakis et al. (Prodromakis, Peh, Papavassiliou, & Toumazou, 2011), etc.Without regarding to any specific f(x(t)), the state space representation of the oscillator of any type can be now given by

Conclusion
The fractional domain generalization of the memristor-based parametric Wien oscillators by using fractional memristor have been performed in this work.We have found that the effects of α to f avr and sustained oscillating condition are insignificant.However, it has been found that α significantly affects Δf osc , the phase difference between M(t) and V out (t), the asymptotic magnitude of V out (t) and the difference between such asymptotic value and its early state counterpart.From the analysis results which have been verified by using the SPICE HP memristor model-based simulation, we promote the usage of the fractional memristor with α > 1 due to its many pros and the further generalization by using fractional capacitors if more degree of freedom is desired.We also give a discussion on the effect of nonlinear dopant drift which has been previously ignored by Talukdar et al. (Talukdar et al., 2011).Therefore, this research gives a simple, yet effective design methodology of the parametric Wien family and clearly demonstrates an interesting application of the state of the art fractional memristor.You are free to: Sharecopy and redistribute the material in any medium or format.Adaptremix, transform, and build upon the material for any purpose, even commercially.The licensor cannot revoke these freedoms as long as you follow the license terms.
Under the following terms: Attribution -You must give appropriate credit, provide a link to the license, and indicate if changes were made.You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.

No additional restrictions
You may not apply legal terms or technological measures that legally restrict others from doing anything the license permits.

Figure 11 .
Figure 11.The schematic of SPICE HP memristor model.
the Joklecar's window function which can be given by(Joglekar & Wolf, 2009) f ðxðtÞÞ ¼ 1 À ð2xðtÞ À 1Þ 2p (47)where p ≥ 1 and p 2 fIg, the resulting V out (t) of type A oscillator can be simulated as depicted in Figure24where p = 1, α = 1.25, β = γ = 0.95, C 1β = 3.2 μF•s β−1 and C 2γ = 3.2 μF•s γ−1 have been assumed.Since it can be seen that V out (t) become eventually distorted, it can be stated such nonlinearity near boundaries causes the oscillator to suffer a nonlinear distortion.For validation, comparison with the order 1.25 SPICE HP memristor model-based benchmark has been made as also shown in these figures where the strong agreements, which verify the above statement and further validates of our results with other α apart from 1, can be observed.It should be mentioned here that the SPICE HP memristor model with Joklecar's window function has been adopted in this scenario unlike the previous simulations which such model with linear window function has been used as the nonlinear dopant drift has been assumed.

Cogent
Engineering (ISSN: 2331-1916) is published by Cogent OA, part of Taylor & Francis Group.Publishing with Cogent OA ensures: • Immediate, universal access to your article on publication • High visibility and discoverability via the Cogent OA website as well as Taylor & Francis Online • Download and citation statistics for your article • Rapid online publication • Input from, and dialog with, expert editors and editorial boards • Retention of full copyright of your article • Guaranteed legacy preservation of your article • Discounts and waivers for authors in developing regions Submit your manuscript to a Cogent OA journal at www.CogentOA.com

Table 1 .
F AVR (HZ) of fractional memristor-based Type A Wien oscillator

Table 2 .
F AVR (HZ) of fractional memristor-based Type B Wien oscillator

Table 3 .
F AVR (HZ) of fractional memristor-based Types C and D Wien oscillator

Table 6 .
R3/R4 of fractional memristor-based Types C and D Wien oscillator

Table 9 .
ΔFOSC (HZ) of fractional memristor-based Types C and D Wien oscillator