Dynamical Properties of Fractional-Order Memristor

: The properties of a fractional-order memristor is studied, and the inﬂuences of parameters are analyzed and compared. The results reﬂect that the resistance value of a fractional-order memristor can be a ﬀ ected by fraction-order, frequency, the switch resistor ratio, average mobility and so on. In addition, the circuit of a fractional-order memristor that is serially connected and connected in parallel with inductance and capacitance are studied. Then, the current–voltage characteristics of a simple series one-port circuits that are composed of a fractional-order memristor and a capacitor, or composed of a fractional-order memristor and a inductor are studied separately. The results demonstrate that at the periodic excitation, the memristor in the series circuits will have capacitive properties or inductive properties as the fractional order changes, the dynamical properties can be used in a memristive circuit.


Introduction
Fractional calculus, an important branch of mathematics, was born in 1695 and appeared almost simultaneously to classic calculus. Fractional calculus, in a narrow sense, mainly includes fractional differentials and fractional integrals, and it broadly includes fractional differences and fractional sum quotients. Since the theory of fractional calculus has been successfully applied to various fields in recent years, people have gradually discovered that fractional calculus can describe some non-classical phenomena in the fields of natural science and engineering applications. The current popular areas of fractional calculus include fractional numerical algorithms and fractional synchronization.
The successful development of a memristor has provided a new avenue for electronic technology and information technology, and it is expected to realize new functions. Their non-volatility makes memristors play a key role in memory, neural networks, and pattern recognition. Chua firstly defined a memristor (MR) in 1971 [1], and HP Labs reported the successful fabrication of nanoscale memristive devices [2]. A memristor is the fourth two-terminal fundamental circuit element with information storage ability and, as such, has attracted immense worldwide interest from both industry and academia [3]. Applications of MRs are used in many fields such as filter design, programmable logic, biological systems, and neural systems such as neural synaptic weighting with a pulse-based memristor circuit [4], Boolean logic operations and computing circuits based on memristors [5], and the voltage-current relationship of active memristors and frequency [6]. A generalized boundary condition memristor model was proposed in [7]. Research on a coupling behavior-based series-parallel flux-controlled memristor has been also conducted [8][9][10]. Two types of nanoscale nonlinear memristor models and their series-parallel circuit have been investigated [11], as have the characteristics of a memristor and its application in the circuit design [12][13][14][15][16][17]. First order mem-circuits have been where m = [α] + 1 and Γ(m) is the Euler's gamma function; when α ∈ (0, 1): Then, we obtain: The real part and imaginary part of the sine and cosine functions can be obtained by separating Equation (3).
Some basic properties of Caputo fractional calculus are as follows.
(1) a D α t [u f (t) + vg(t)] = u a D α t f (t) + v a D α t g(t) (2) a D α t a D β t f (t) = a D β t a D α t f (t) = a D α+β t f (t), In which the constant K 0, and the Caputo FD is C a D α t K = 0, α > 0.

The Model of Fractional-Order Memristor
The structure of an MR and its symbol are shown in Figure 1.
Symmetry 2019, 11, x FOR PEER REVIEW 2 of 12 synchronization [20][21][22], and coupling as a third relation in memristive systems have been proposed [23][24][25][26][27][28][29][30]. However, there is less research on fractional-order memristor, and there has been no literature physical background. In this paper, the analytical solution of parameter expression is derived by using a fractional-order memristor model in Section 2. The properties of fractional-order memristor is studied in Section 3, The results reflect that the resistance value of a fractional-order memristor can be affected by fraction-order, frequency, the switch resistor ratio, average mobility, and so on; additionally, the circuit of α M serially connected and connected in parallel with α M , L and C are also studied separately. Finally, the conclusion is presented in Section 4.

The Fractional Derivative
The α order Caputo derivative is defined as: (2) Then, we obtain: The real part and imaginary part of the sine and cosine functions can be obtained by separating Equation (3).
Some basic properties of Caputo fractional calculus are as follows. The mathematical model are defined as follows  The mathematical model are defined as follows

The Model of Fractional-Order Memristor
In which D −α t is the α integral of x(t), M α (x) is a fractional-order memristor when supposing the input current i(t) = I m sin(ωt), switch resistor R ON = 100 Ω and R OFF = 10 kΩ, the current I m = 0.2 mA D = 10 nm, the average mobility u v = 10 −14 m.s −1 V −1 , the length x 0 = 0.01, and the frequency ω = 5 rad/s. The transient current curve and the voltage curve of the memristor are shown in Figure 2a,b respectively. (4c) In which is a fractional-order memristor when supposing the input current  The simulation result for when α has different values is shown in Figure 3a. It was found the an MR possesses memristive properties, with pinched hysteresis loops forming inclined "8", and this property can be used to realize signal storage or computing. From Figure 3b, it can be seen that the smaller the value of fractional order, the greater the dynamic range amplitude of the resistance was. The simulation result for when α has different values is shown in Figure 3a. It was found the an MR possesses memristive properties, with pinched hysteresis loops forming inclined "8", and this property can be used to realize signal storage or computing. From Figure 3b, it can be seen that the smaller the value of fractional order, the greater the dynamic range amplitude of the resistance was. (4c) In which   The simulation result for when α has different values is shown in Figure 3a. It was found the an MR possesses memristive properties, with pinched hysteresis loops forming inclined "8", and this property can be used to realize signal storage or computing. From Figure 3b, it can be seen that the smaller the value of fractional order, the greater the dynamic range amplitude of the resistance was.  The characteristic curves of a fractional-order memristor with different parameters are shown in Figure 4. From Figure 4a, it can be seen that the difference of resistance decreased as the frequency ω increased. From Figure 4b, it can be seen that the switch resistor increased and the curves of v − i were  The characteristic curves of a fractional-order memristor with different parameters are shown in Figure 4. From Figure 4a, it can be seen that the difference of resistance decreased as the frequency ω increased. From Figure 4b, it can be seen that the switch resistor increased and the curves of i v − were inclined to right. From Figure 4c, it can be seen that the difference of resistance increased as the value of v μ increased.

The Properties of Fractional-Order Memristor
The fractional-order memristor in six cases of connection, serially and in parallel, are discussed in this section.

The Two Fractional-Order Memristors in Serial
The serial circuit of two fractional-order memristors is shown in Figure 5.

The Properties of Fractional-Order Memristor
The fractional-order memristor in six cases of connection, serially and in parallel, are discussed in this section.

The Two Fractional-Order Memristors in Serial
The serial circuit of two fractional-order memristors is shown in Figure 5.
The characteristic curves of a fractional-order memristor with different parameters are shown in Figure 4. From Figure 4a, it can be seen that the difference of resistance decreased as the frequency ω increased. From Figure 4b, it can be seen that the switch resistor increased and the curves of i v − were inclined to right. From Figure 4c, it can be seen that the difference of resistance increased as the value of v μ increased.

The Properties of Fractional-Order Memristor
The fractional-order memristor in six cases of connection, serially and in parallel, are discussed in this section.

The Two Fractional-Order Memristors in Serial
The serial circuit of two fractional-order memristors is shown in Figure 5.   By choosing the current i(t) = I m sin(ωt), based on the Caputo differential, Euler's formulas, and the separating variables method, we obtain: When t >> 1, Equation (11) can be simplified as: .
Then, the two sides of Equation (12) can be integrated to get: Symmetry 2020, 12, 437

of 11
By choosing the parameters M1(u , and α = 0.98, one can obtain the curves of v − i that are are shown in Figure 6. Then, the two sides of Equation (12) can be integrated to get: By choosing the parameters M1(  Figure 6. It was found the two serial MRs also possessed the memristive properties with pinched hysteresis loops behaving as inclined "8", as shown in Fig. 6. It can be seen that the values of M1, M2 and M12 are increased with the value of v(t), and the value of M12 is bigger than those of M1 and M2.

The Circuit of Fractional-Order MR in parallel
The circuit of two fractional-order MRs in parallel is shown in Figure 7. According to Equations (4b) and (4c), it be written as: When 1 >> t ,Equation (11) can be simplified as: and Figure 6. Simulation results of two fractional-order, serially connected memristors v(t) − i(t) curves.
It was found the two serial MRs also possessed the memristive properties with pinched hysteresis loops behaving as inclined "8", as shown in Figure 6. It can be seen that the values of M1, M2 and M12 are increased with the value of v(t), and the value of M12 is bigger than those of M1 and M2.

The Circuit of Fractional-Order MR in Parallel
The circuit of two fractional-order MRs in parallel is shown in Figure 7. According to Equations (4b) and (4c), it be written as: Symmetry 2019, 11, x FOR PEER REVIEW 6 of 12 Then, one can set and we can obtain: The value of ) (x M α can be calculated. and ω = 5. The simulation result when using these parameters is shown in Figure 8. By choosing the voltage v(t) = V m sin(ωt), one can obtain: When t >> 1,Equation (11) can be simplified as: Then, one can set and we can obtain: The value of M α (x) can be calculated.
To further study the dynamic behaviors of this MR circuit in parallel, the parameters were configured as the following: the input voltage v(t) = V m sin ωt, and the parameters R ON1 = 100 Ω, R OFF1 = 10 kΩ, V m = 3 V, D = 10 nm, u v = 10 −14 m.s −1 V −1 , R ON2 = 120 Ω, R OFF2 = 18 kΩ, and ω = 5. The simulation result when using these parameters is shown in Figure 8.
The value of ) (x M α can be calculated. and ω = 5. The simulation result when using these parameters is shown in Figure 8. It was found the two MRs that were connected in parallel also possessed the memristive properties with pinched hysteresis loops which are hown in Figure 8.

The Circuit of Fractal-Order Memristor and Capacitor That Are Serially Connected
The fractional-order memristor M α (t) and serially connected capacitor C are shown in Figure 9. By assume the current i(t) = I m sin ωt, M expresses the memristor, and C is the capacitor.
It was found the two MRs that were connected in parallel also possessed the memristive properties with pinched hysteresis loops which are hown in Figure 8.

The Circuit of Fractal-Order Memristor and Capacitor That Are Serially Connected
The fractional-order memristor ) (t M α and serially connected capacitor C are shown in  we can obtain:  We can obtain: where The simulation result is shown in Figure 10. From Figure 10a, it can be seen as the parameter α decreased, the area of hysteresis loops increased and the difference of resistance increased with same current. From Figure 10b, it can be seen that if the parameters α and C were not varied, as ω increased, the area of hysteresis loops decreased. From Figure 10c, it can be seen that if the parameters α and ω were not varied, the area of hysteresis loops was not almost changed and the capacitor has little effect on the circuit.

The Circuit of Fractal-Order Memristor and Capacitor That Were Connected in Parallel
A memristor and a capacitor which were connected in parallel are shown in Figure 11, we can assume that By applying Kirchhoff's current law (KCL), the current of circuit can be written as: (17)

The Circuit of Fractal-Order Memristor and Capacitor That Were Connected in Parallel
A memristor and a capacitor which were connected in parallel are shown in Figure 11, we can assume that v(t) = V m sin(ωt), ω(0) = 0,V m is the voltage magnitude.

The Circuit of Fractal-Order Memristor and Capacitor That Were Connected in Parallel
A memristor and a capacitor which were connected in parallel are shown in Figure 11, we can assume that By applying Kirchhoff's current law (KCL), the current of circuit can be written as: (17) Figure 11. The circuit of a memristor and a capacitor that were connected in parallel. By applying Kirchhoff's current law (KCL), the current of circuit can be written as: Additionally, when the parameters α and ω were chosen as different values, it can be seen the circuit of a memristor and a capacitor that were connected in parallel also possessed the memristive properties with pinched hysteresis loops which are shown in simulations in Figure 12a,b.

The Circuit of Fractal-Order Memristor and Inductor That Are Serially Connected
The circuit of the memristor and serially connected inductor are shown in Figure 13, where we assumed ) sin( The voltage was: and The simulation result is shown in Figure 14. From Figure 14a, it can be seen that as the parameters ω and L were not varied, because as α decreased, the area of hysteresis loops increased and the difference of resistance increased with the same current. From Figure 14b, it can be seen that if the parametersα and L were not varied, as ω increased, the area of the hysteresis loops decreased. From Figure 14c, it can be seen that if the parameters α and ω were not varied, the area of the hysteresis loops was not almost changed, and it is shown that the inductor barely affected the circuit.

The Circuit of Fractal-Order Memristor and Inductor That Are Serially Connected
The circuit of the memristor and serially connected inductor are shown in Figure 13, where we assumed ) sin( The voltage was: and The simulation result is shown in Figure 14. From Figure 14a, it can be seen that as the parameters ω and L were not varied, because as α decreased, the area of hysteresis loops increased and the difference of resistance increased with the same current. From Figure 14b, it can be seen that if the parametersα and L were not varied, as ω increased, the area of the hysteresis loops decreased. From Figure 14c, it can be seen that if the parameters α and ω were not varied, the area of the hysteresis loops was not almost changed, and it is shown that the inductor barely affected the circuit. Figure 13. The circuit of the memristor and serially connected inductor. Figure 13. The circuit of the memristor and serially connected inductor.
The voltage was: and The simulation result is shown in Figure 14. From Figure 14a, it can be seen that as the parameters ω and L were not varied, because as α decreased, the area of hysteresis loops increased and the difference of resistance increased with the same current. From Figure 14b, it can be seen that if the parameters α and L were not varied, as ω increased, the area of the hysteresis loops decreased. From Figure 14c, it can be seen that if the parameters α and ω were not varied, the area of the hysteresis loops was not almost changed, and it is shown that the inductor barely affected the circuit.

The Circuit of Fractal-Order Memristor and Inductor Connected in Parallel.
A memristor and a capacitor that were connected in parallel are shown in Figure 15. By V can be found as the voltage amplitude.
By applying Kirchhoff's current law (KCL), the current of circuit can be written as: The simulation is shown in Figure 16a,b. It can be seen the circuit of the memristor and inductor that were connected in parallel also possessed memristive properties with pinched hysteresis loops. From Figure 16a, it can be seen the area of pinched hysteresis loops increased withα increased.
From Figure 16b, it can be seen the area of the pinched hysteresis loops decreased with the ω increased. Figure 15. The circuit of a memristor and a inductor that were connected in parallel.

The Circuit of Fractal-Order Memristor and Inductor Connected in Parallel
A memristor and a capacitor that were connected in parallel are shown in Figure 15. By assuming the voltage v(t) = V m sin(ωt), ω(0) = 0, V m can be found as the voltage amplitude. By applying Kirchhoff's current law (KCL), the current of circuit can be written as:

The Circuit of Fractal-Order Memristor and Inductor Connected in Parallel.
A memristor and a capacitor that were connected in parallel are shown in Figure 15. By V can be found as the voltage amplitude.
By applying Kirchhoff's current law (KCL), the current of circuit can be written as: The simulation is shown in Figure 16a,b. It can be seen the circuit of the memristor and inductor that were connected in parallel also possessed memristive properties with pinched hysteresis loops. From Figure 16a, it can be seen the area of pinched hysteresis loops increased withα increased.
From Figure 16b, it can be seen the area of the pinched hysteresis loops decreased with the ω increased. Figure 15. The circuit of a memristor and a inductor that were connected in parallel. Figure 15. The circuit of a memristor and a inductor that were connected in parallel.
The simulation is shown in Figure 16a,b. It can be seen the circuit of the memristor and inductor that were connected in parallel also possessed memristive properties with pinched hysteresis loops. From Figure 16a, it can be seen the area of pinched hysteresis loops increased with α increased. From Figure 16b, it can be seen the area of the pinched hysteresis loops decreased with the ω increased. Symmetry 2019, 11, x FOR PEER REVIEW 11 of 12

Conclusion
In summary, this paper has presented a fractional-order memristor model and verifies the three essential characteristics of a fractional-order memristor. In addition, the properties of fractional-order memristor have been described. In a simple memristive series circuit, with a change of the fractional derivative order, the series circuit of a fractional-order memristor and a capacitor or inductor shows a conversion from a pure capacitor circuit to a memristive circuit. The series circuit shows conversions of purely inductive and memristive circuits. Here, analytical solutions were derived by a fractional-order memristor, the properties of fractional-order memristor model parameters were obtained, and simulation results were given. The results showed that material properties determine the order of the fractional derivative, so the best memory capacity of a physical memristor can be achieved by finding materials that are compatible with the excitation frequency.

Conclusions
In summary, this paper has presented a fractional-order memristor model and verifies the three essential characteristics of a fractional-order memristor. In addition, the properties of fractional-order memristor have been described. In a simple memristive series circuit, with a change of the fractional derivative order, the series circuit of a fractional-order memristor and a capacitor or inductor shows a conversion from a pure capacitor circuit to a memristive circuit. The series circuit shows conversions of purely inductive and memristive circuits. Here, analytical solutions were derived by a fractional-order memristor, the properties of fractional-order memristor model parameters were obtained, and simulation results were given. The results showed that material properties determine the order of the fractional derivative, so the best memory capacity of a physical memristor can be achieved by finding materials that are compatible with the excitation frequency.