Next Article in Journal
Editorial for Special Issue “New Advancements in Pure and Applied Mathematics via Fractals and Fractional Calculus”
Previous Article in Journal
A New Method of Quantifying the Complexity of Fractal Networks
Previous Article in Special Issue
Rational Approximations of Arbitrary Order: A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Oscillators Based on Fractional-Order Memory Elements

Faculty of BERG, Technical University of Kosice, Nemcovej 3, 042 00 Kosice, Slovakia
Submission received: 22 April 2022 / Revised: 16 May 2022 / Accepted: 22 May 2022 / Published: 24 May 2022
(This article belongs to the Special Issue Fractional-Order Circuit Theory and Applications)

Abstract

:
This paper deals with the new oscillator structures that contain new elements, so-called memory elements, known as memristor, meminductor, and memcapacitor. Such circuits can exhibit oscillations as well as chaotic behavior. New mathematical models of fractional-order elements and whole oscillator circuits are proposed as well. An illustrative example to demonstrate the oscillations and the chaotic behavior through the numerical solution of the fractional-order circuit model is provided.

1. Introduction

Although classical electrical elements such as a resistor, capacitor, and inductor have been known for a long time, memristors, memcapacitors, and meminductors are relatively new nonlinear elements with memory [1]. As a novel memory device, the memristor was postulated by Leon Chua in 1971 [2] and manufactured for the first time in 2008 by HP Labs [3]. Until this time, an investigation of the memristor concept was very limited due to the lack of a solid-state implementation of this postulated device.
The memcapacitor and meminductor are also members of a huge family of new circuit elements postulated by Leon Chua in 1978 [4]. The idea of Chua extended the concept of memory elements in an electrical circuit to capacitive and inductive systems, respectively. The memcapacitor and the meminductor were formally defined and described in 2009 [5]. Some elements of the electronic circuit, namely the memcapacitor and the meminductor, require not only the well-known four state variables but also the time integrals of the electric charge and flux [6]. These new state variables, due to integration, lead us to so-called “memory” devices, which are a particular class of higher-order elements (devices) and belong to a broad group of memory systems [4]. They are passive memory devices that can store information without a power supply. Currently, the applications of these memory devices in nonlinear circuits have gained a great deal of attention, and their potential value has attracted many researchers. However, the absence of semiconductor implementation of memcapacitors and meminductors prevents the use of unique functions of these devices in practical implementation. Their properties were so far investigated mainly via mathematical models, equivalent circuits, or emulators [7,8,9,10]. Such investigation is not accurate because it analyzes just their approximation, not fundamental elements.
There are also a considerable number of electrical circuits where non-integer order (or fractional) calculus can be used (see, for example, [11,12,13,14], etc.), with classical electrical circuit theory being limited to variables u (voltage), i (current), q (charge), and ϕ (flux), which are used to describe all four essential components (resistor, capacitor, inductor, and memristor). Models of the non-integer order can also describe the type noted above of memristive systems [15,16,17,18,19,20,21]. In addition, in practice, there is no ideal electrical element, and almost all electrical elements lie between two ideal ones, for example, a fractor (resistor/capacitor) or a fractductor (resistor/inductor) [12,22,23]. Here, we consider it also for the real memristive elements. That means that all real memristive elements should lie in between two ideals, similar to the classical electrical elements noted above. We used a mathematical model of supposed new fractional-order memory elements in the specific oscillator circuit. According to the author’s best knowledge, such an oscillator structure described in this article was used for the first time.
In this article, a new oscillator based on memory devices is studied. The rest of the manuscript is structured as follows: In Section 2 a definition of the fractional calculus and method for the numerical solution of the initial value problem is described. Section 3 presents the fractional-order elements models. In Section 4 the new fractional-order models of the oscillator circuits are proposed. In Section 5 some ideas for further research are discussed. Section 6 concludes this article with some additional comments.

2. Preliminaries

2.1. Definition of Fractional-Order Operator

Fractional calculus has been known since regular calculus, probably with the first evidence dated 30 September 1695, in letter correspondence between Gottfried W. Leibniz and Guillaume de l’Hospital. They mentioned a half order derivative for the first time. The fractional calculus is a generalization of differentiation and integration to common non-integer γ -order, γ R , operator a D t γ , where a and t of interval [ a , t ] are the bounds of the join operation (fractional-order integrals for γ < 0 and derivatives for γ > 0 ).
There are many different definitions for the fractional-order operator a D t γ , but in this paper, we will limit ourselves only to two fundamentals: Caputo’s definition (CD), and Grünwald–Letnikov’s definition (GLD).
The CD can be written as [24]:
a D t γ f ( t ) = 1 Γ ( n γ ) a t f ( n ) ( τ ) ( t τ ) γ n + 1 d τ , n 1 < γ < n ,
where Γ ( . ) denotes Euler gamma function. The CD can be used for electrical circuits where non-integer order derivatives are used in the fractional (non-integer) order model of electrical circuit elements. The main benefit is that the initial conditions for fractional differential equations with Caputo derivatives are the same as for ordinary differential equations, i.e., f ( n ) ( 0 ) = c n , n N .
The GLD is given as follows [24,25]:
a D t γ f ( t ) = lim h s 0 1 h s γ i = 0 t a h s ( 1 ) i γ i f ( t i h s ) ,
where z is the floor function, i.e., the greatest integer smaller than z, and γ i are the binomial coefficients with γ 0 = 1 . This form of the definition is beneficial for obtaining a numerical solution of the fractional differential equation.

2.2. Numerical Solution of Fractional Differential Equation

Based on the fact that both definitions, CD, and GLD, are equivalent for a wide class of the functions, for numerical calculation of the fractional-order derivative, we can use the relation (3) derived from the GLD (2). The relation for the numerical approximation of the γ th derivative at the points k h s , ( k = 1 , 2 , 3 , ) has the following form [24]:
( k L m / h s ) D t k γ f ( t ) h s γ i = 0 k b i ( γ ) f ( t k i ) ,
where L m is the “memory length”, t k = k h s , h s is the time step of calculation (definition (3) is valid only as h s tends towards 0 and that the accuracy of the simulation depends on the value of h s ), and b i ( γ ) ( i = 0 , 1 , 2 , ) are the binomial coefficients. For their calculation we may use the following expression:
b 0 ( γ ) = 1 , b i ( γ ) = 1 1 + γ i b i 1 ( γ ) .
Thus, general numerical solution of the fractional differential equation
0 D t γ u ( t ) = f ( t , u ( t ) ) , 0 t T ,
can be expressed as follows [22]:
u ( t k ) = f t k , u ( t k ) h s γ i = 1 k b i ( γ ) u ( t k i ) .
For the memory term expressed by the sum in (5), a “short memory” principle for various memory lengths of L m can be used. An evaluation of the effect of the memory length and convergence relation of the error between short and long memory was described in [24]. However, this article uses a whole (long) memory to preserve calculation accuracy.

3. Fractional-Order Memristive Elements

Many authors have so far studied the real capacitor and the real inductor, and the experimental evidence of their non-integer order of models has undoubtedly been confirmed [13,26]. The real memristor, memcapacitor, and meminductor have not yet been experimentally studied, only using their simulation models or emulators [10,27,28]. It is because they do not exist as a single component except, for the memristor element constructed in the HP Lab in 2008. However, there are many electrical circuits where memristors, memcapacitors, and meminductors have been used on the theoretical level [29,30,31,32].
Using the relations described in [6] for memcapacitor and meminductor and the well-known connections between the four essential components (resistor, capacitor, inductor, memristor), we can obtain the next floor using the square symmetry shown in Figure 1.
However, as we can observe in Figure 1, there could be some additional electrical element that has not been discovered yet.

3.1. Memristor

The memristor used in this article is a flux-controlled memristor that is characterized by the relation
i ( t ) = W ( ϕ ( t ) ) u ( t ) ,
where W ( ϕ ( t ) ) = d q ( ϕ ) / d ϕ is a memductance of the memristor. Here, we consider the following model:
q ( ϕ ) = m 1 ϕ ( t ) + m 2 ϕ 3 ( t ) ,
where m 1 and m 2 are real constants. The memductance function W that is obtained from the q ( ϕ ) function is [30]:
W ( ϕ ) = d q ( ϕ ) d ϕ = m 1 + 3 m 2 ϕ 2 ( t ) .
Moreover, for such a memristor, a monotone-increasing piecewise-linear characteristic was assumed in [34] and used in [35].
Similarly, as in the case of the real capacitor and real inductor [13,26,36], the memristor is also not an ideal element and we may consider the fractional-order model of this element [37].

3.2. Memcapacitor

An ideal charge-controlled memcapacitor is defined as [5,31]:
u ( t ) = C M 1 ( σ ) q ( t ) , C M 1 ( σ ) = d ϕ ( σ ) / d σ ,
where q ( t ) is the charge on the memcapacitor and u ( t ) is the corresponding voltage across memcapacitor at time t, σ is the time-domain integral of electric charge q passing through the memcapacitor, ϕ ( σ ) is the flux that goes through the memcapacitor, and C M 1 ( σ ) is the inverse memcapacitance, which depends on the state of the device. In this paper, a charge-controlled memcapacitor is formulated in accordance with its definition [31]:
u ( t ) = ( α + β σ 2 ) q ( t ) = α + β t 0 t q ( τ ) d τ 2 q ( t ) ,
where α and β are constants, and their units are F 1 and ( C 2 S 2 F ) 1 , respectively.
Following a linear capacitor model proposed by Westerlund and Ekstam in 1994, for a general input voltage u ( t ) , applied at t = 0 , the current is [26]:
i ( t ) = C 0 D t δ u ( t ) , 0 < δ 1 , for t > 0 ,
where C is the capacitance of the capacitor with unit [F/s 1 δ ]. It is related to the kind of dielectric. Another constant δ (order), δ R , is related to the losses of the capacitor. We may predict a fractional-order model of the real memcapacitor as well [15,17].

3.3. Meminductor

Similar to the definition of a memcapacitor, the ideal flux-controlled meminductor is given as [5,31]:
i ( t ) = L M 1 ( ρ ) ϕ ( t ) , L M 1 ( ρ ) = d q ( ρ ) / d ρ ,
where ϕ ( t ) and i ( t ) denote the flux and current go through a meminductor at time t, ρ is the time-domain integral of electric flux ϕ passing through the meminductor, q ( ρ ) is the charge that goes through the meminductor, which is a function of ρ , and L M 1 ( ρ ) is the inverse meminductance, which depends on their inner variables. By providing a concrete expression of L M 1 ( ρ ) , the flux-controlled mathematic model is shown as follows [31]:
i ( t ) = ( α + β ρ ( t ) ) ϕ ( t ) = α + β t 0 t ϕ ( τ ) d τ ϕ ( t ) ,
where α and β are constants, and their units are H 1 and ( W b S H ) 1 , respectively.
For a linear inductor model suggested in [13,36], the voltage is
u ( t ) = L 0 D t μ i ( t ) , 0 < μ 1 , for t > 0 ,
where L is the inductance of the inductor with unit [H/s 1 μ ]. It is related to the kind of coil core material and depends on the geometry of inductor. Another constant μ (order), μ R , is related to the proximity effect of the inductor. Similarly, we may predict a fractional-order model of the real meminductor as well [15].

4. Models of the Fractional-Order Chaotic Systems

4.1. Memcapacitor–Meminductor Oscillator

Let us start with the oscillator structure consisting of two memory elements, namely meminductor and memcapacitor, which was designed in [31].
Based on the aforementioned charge-controlled memcapacitor and flux-controlled meminductor, a novel chaotic circuit is depicted in Figure 2. The suggested oscillator consists of a meminductor ( L M ), a memcapacitor ( C M ), a capacitor ( C 1 ), a linear resistor (R), and a negative resistor ( G ). This circuit is also supplied by energy for maintaining the oscillating state.
According to Kirchhoff’s laws for two current nodes and one voltage loop, the following set of differential equations is obtained [31]:
d u 1 ( t ) d t = 1 C 1 i LM ( t ) u 1 ( t ) R , d q CM ( t ) d t = G u CM ( t ) i LM ( t ) , d ϕ LM ( t ) d t = u CM ( t ) u 1 ( t ) .
Substituting Equations (10) and (13) into the above Equation (15), then adding the two well-known relations: d ρ L M ( t ) / d t = ϕ L M ( t ) and d σ C M ( t ) / d t = q C M ( t ) to them ( ρ L M ( t ) is the integral of flux ϕ L M ( t ) passing through the meminductor, σ C M ( t ) is the integral of charge q C M ( t ) passing through the memcapacitor), we obtain the differential equations [31]:
d u 1 ( t ) d t = 1 C 1 ( α ϕ LM ( t ) + β ρ LM ( t ) ϕ LM ( t ) ) u 1 ( t ) R , d q CM ( t ) d t = G α q CM ( t ) + β q CM ( t ) σ CM 2 ( t ) α ϕ LM ( t ) + β ρ LM ( t ) ϕ LM ( t ) , d ϕ LM ( t ) d t = α q CM ( t ) + β q CM ( t ) σ CM 2 ( t ) u 1 ( t ) , d ρ LM ( t ) d t = ϕ LM ( t ) , d σ CM ( t ) d t = q CM ( t ) .
For additional dynamical analysis, by setting τ = t / R C 1 and scale transformations x = u 1 , y = q CM , z = ϕ LM , w = ρ LM , v = σ CM , the dynamical system (16) converts into dimensionless form as follows [31]:
d x ( τ ) d τ = a z ( τ ) + b z ( τ ) w ( τ ) x ( τ ) , d y ( τ ) d τ = d y ( τ ) + e y ( τ ) v 2 ( τ ) f z ( τ ) g z ( τ ) w ( τ ) , d z ( τ ) d τ = h y ( τ ) + j y ( τ ) v 2 ( τ ) k x ( τ ) , d w ( τ ) d τ = n 1 z ( τ ) , d v ( τ ) d τ = n 2 y ( τ ) ,
where its parameters are: a = α R , b = β R , d = G α R C 1 , e = G β R C 1 , f = α R C 1 , g = β R C 1 , h = α R C 1 , j = β R C 1 , k = R C 1 , n 1 = R C 1 , and n 2 = R C 1 .
The nonlinear dynamical system (17), depicted in Figure 2, will exhibit chaotic attractor for the following parameters [31]: a = 1.73 , b = 2.04 , d = 0.46 , e = 0.04 , f = 0.67 , g = 0.19 , h = 0.48 , j = 0.52 , k = n 1 = n 2 = 0.21 , and the initial values setting: x ( 0 ) = 0.2 , y ( 0 ) = 0.5 , z ( 0 ) = 0.45 , w ( 0 ) = 0.1 , v ( 0 ) = 0.5 .
Taking into account the consideration described in the previous section, the fractional-order models of the elements (capacitor, memcapacitor, and meminductor) used in the chaotic oscillator shown in Figure 2 could be applied to Kirchhoff laws. Then, instead of the system (17), we obtain the following set of fractional differential equations:
0 D τ γ 1 x ( τ ) = a z ( τ ) + b z ( τ ) w ( τ ) x ( τ ) , 0 D τ γ 2 y ( τ ) = d y ( τ ) + e y ( τ ) v 2 ( τ ) f z ( τ ) g z ( τ ) w ( τ ) , 0 D τ γ 3 z ( τ ) = h y ( τ ) + j y ( τ ) v 2 ( τ ) k x ( τ ) , 0 D τ γ 4 w ( τ ) = n 1 z ( τ ) , 0 D τ γ 5 v ( τ ) = n 2 y ( τ ) ,
where γ 1 , γ 2 , γ 3 , γ 4 , and γ 5 are the real orders of aforementioned elements used in the chaotic oscillator displayed in Figure 2.
For simulation purposes, a numerical solution of the fractional differential Equation (18) obtained using the relationships (3)–(5), was proposed [38]:
x ( t k ) = [ a z ( t k 1 ) + b z ( t k 1 ) w ( t k 1 ) x ( t k 1 ) ] h s γ 1 i = 1 k b i ( γ 1 ) x ( t k i ) , y ( t k ) = [ d y ( t k 1 ) + e y ( t k 1 ) v ( t k 1 ) 2 f z ( t k 1 ) g z ( t k 1 ) w ( t k 1 ) ] h s γ 2 i = 1 k b i ( γ 2 ) y ( t k i ) , z ( t k ) = [ h y ( t k 1 ) + j y ( t k 1 ) v ( t k 1 ) 2 k x ( t k 1 ) ] h s γ 3 i = 1 k b i ( γ 3 ) z ( t k i ) , w ( t k ) = [ n 1 z ( t k 1 ) ] h s γ 4 i = 1 k b i ( γ 4 ) w ( t k i ) , v ( t k ) = [ n 2 y ( t k 1 ) ] h s γ 5 i = 1 k b i ( γ 5 ) v ( t k i ) ,
where T s i m is the simulation time, h s is the calculation time step, k = 1 , 2 , 3 , N , for N = [ T s i m / h s ] , and initial conditions are ( x ( 0 ) , y ( 0 ) , z ( 0 ) , w ( 0 ) , v ( 0 ) ). The binomial coefficients b i ( γ 1 ) , b i ( γ 2 ) , b i ( γ 3 ) , b i ( γ 4 ) , and b i ( γ 5 ) are calculated according to relation (4), respectively.
Figure 3 depicts the simulation results of the system (18) using numerical solution (19) for the parameters: a = 1.73 , b = 2.04 , d = 0.46 , e = 0.04 , f = 0.67 , g = 0.19 , h = 0.48 , j = 0.52 , k = n 1 = n 2 = 0.21 , orders γ 1 = 0.9 (real capacitor C 1 ), γ 2 = γ 3 = γ 4 = γ 5 = 1 (ideal memcapacitor C M and ideal meminductor L M ), initial conditions x ( 0 ) = 0.2 , y ( 0 ) = 0.5 , z ( 0 ) = 0.45 , w ( 0 ) = 0.1 , v ( 0 ) = 0.5 , step h s = 0.005 , and simulation time 500 s. The Lyapunov exponents of this system with the above parameters, computed according to the algorithm described in [39], have the following values: 0.5807 , 0.1683 , 0.0047 , 0.9647 , 1.0267 , which confirm the system is chaotic because at least one exponent is positive.

4.2. Memristor–Memcapacitor–Meminductor Oscillator

Following the oscillator based on the charge-controlled memcapacitor and the flux-controlled meminductor, which was described in the previous subsection, a novel chaotic circuit containing the flux-controlled memristor (8) instead of the negative resistor was suggested, depicted in Figure 4. The proposed oscillator consists of a meminductor ( L M ), a memristor (M), a memcapacitor ( C M ), a linear resistor (R), and a capacitor ( C 1 ). As in the previous case, this novel circuit is also supplied by energy to maintain the oscillating state.
According to Kirchhoff’s laws for two current nodes and one voltage loop, the following set of differential equations can be obtained:
d u 1 ( t ) d t = 1 C 1 i LM ( t ) u 1 ( t ) R , d q CM ( t ) d t = i M ( t ) i LM ( t ) , d ϕ LM ( t ) d t = u CM ( t ) u 1 ( t ) , d ϕ M ( t ) d t = u CM ( t ) .
Substituting Equations (6), (8), (10) and (13) into the above Equation (20), the same as in the previous case, adding the two well-known relations: d ρ L M ( t ) / d t = ϕ L M ( t ) and d σ C M ( t ) / d t = q C M ( t ) to them ( ρ L M ( t ) is the integral of flux ϕ L M ( t ) passing through the meminductor, σ C M ( t ) is the integral of charge q C M ( t ) passing through the memcapacitor), we get the following set of equations:
d u 1 ( t ) d t = 1 C 1 ( α ϕ LM ( t ) + β ρ LM ( t ) ϕ LM ( t ) ) u 1 ( t ) R , d q CM ( t ) d t = ( m 1 + 3 m 2 ϕ M 2 ( t ) ) α q CM ( t ) + β q CM ( t ) σ CM 2 ( t ) α ϕ LM ( t ) + β ρ LM ( t ) ϕ LM ( t ) , d ϕ LM ( t ) d t = α q CM ( t ) + β q CM ( t ) σ CM 2 ( t ) u 1 ( t ) , d ϕ M ( t ) d t = α q CM ( t ) + β q CM ( t ) σ CM 2 ( t ) , d ρ LM ( t ) d t = ϕ LM ( t ) , d σ CM ( t ) d t = q CM ( t ) .
By setting the same state variables transformation as in the previous example and taking into account the fractional-order models of all three memory elements shown in Figure 4, the dynamical system (21) can be mapped into dimensionless form as follows:
0 D τ γ 1 x ( τ ) = a z ( τ ) + b z ( τ ) w ( τ ) x ( τ ) , 0 D τ γ 2 y ( τ ) = d y ( τ ) + e y ( τ ) v 2 ( τ ) + y ( τ ) u 2 ( τ ) ( p + r v 2 ( τ ) ) f z ( τ ) g z ( τ ) w ( τ ) , 0 D τ γ 3 z ( τ ) = h y ( τ ) + j y ( τ ) v 2 ( τ ) k x ( τ ) , 0 D τ γ 4 u ( τ ) = h y ( τ ) + j y ( τ ) v 2 ( τ ) , 0 D τ γ 5 w ( τ ) = n 1 z ( τ ) , 0 D τ γ 6 v ( τ ) = n 2 y ( τ ) ,
where parameters are: a = α R , b = β R , d = m 1 α R C 1 , e = m 1 β R C 1 , f = α R C 1 , g = β R C 1 , h = α R C 1 , j = β R C 1 , k = R C 1 , p = 3 m 2 α R C 1 , r = 3 m 2 β R C 1 , n 1 = R C 1 , and n 2 = R C 1 , and where γ 1 , γ 2 , γ 3 , γ 4 , γ 5 , and γ 6 are the real orders of the elements used in the oscillator circuit. As in system (18), the equilibrium point of the dynamical system (22) is in the origin.
Numerical solutions of the fractional differential Equations (22) obtained using relationships (3)–(5) are given as:
x ( t k ) = [ a z ( t k 1 ) + b z ( t k 1 ) w ( t k 1 ) x ( t k 1 ) ] h s γ 1 i = 1 k b i ( γ 1 ) x ( t k i ) , y ( t k ) = [ d y ( t k 1 ) + e y ( t k 1 ) v ( t k 1 ) 2 + y ( t k 1 ) u ( t k 1 ) 2 ( p + r v ( t k 1 ) 2 ) f z ( t k 1 ) g z ( t k 1 ) w ( t k 1 ) ] h s γ 2 i = 1 k b i ( γ 2 ) y ( t k i ) , z ( t k ) = [ h y ( t k 1 ) + j y ( t k 1 ) v ( t k 1 ) 2 k x ( t k 1 ) ] h s γ 3 i = 1 k b i ( γ 3 ) z ( t k i ) , u ( t k ) = [ h y ( t k 1 ) + j y ( t k 1 ) v ( t k 1 ) 2 ) ] h s γ 4 i = 1 k b i ( γ 4 ) z ( t k i ) , w ( t k ) = [ n 1 z ( t k 1 ) ] h s γ 5 i = 1 k b i ( γ 5 ) w ( t k i ) , v ( t k ) = [ n 2 y ( t k 1 ) ] h s γ 6 i = 1 k b i ( γ 6 ) v ( t k i ) ,
where T s i m is the simulation time, h s is the calculation time step, k = 1 , 2 , 3 , N , for N = [ T s i m / h s ] , and initial conditions are ( x ( 0 ) , y ( 0 ) , z ( 0 ) , u ( 0 ) , w ( 0 ) , v ( 0 ) ). The binomial coefficients b i ( γ 1 ) , b i ( γ 2 ) , b i ( γ 3 ) , b i ( γ 4 ) , b i ( γ 5 ) and b i ( γ 6 ) are calculated using Equation (4), respectively.
Figure 5 shows simulation results of the system (18) in state space, using numerical solution (19), for the parameters: a = 1.73 , b = 2.04 , d = 0.46 , e = 0.04 , f = 0.67 , g = 0.19 , h = 0.48 , j = 0.52 , k = n 1 = n 2 = 0.21 , m 1 = 1 , m 2 = 0 , orders γ 1 = γ 2 = γ 3 = γ 4 = γ 5 = 0.9 , initial conditions x ( 0 ) = 0.7 , y ( 0 ) = 0.5 , z ( 0 ) = 0.45 , u ( 0 ) = 0.3 , w ( 0 ) = 0.1 , v ( 0 ) = 0.5 , step h s = 0.005 , and simulation time 500 s.
Figure 6 depicts the simulation results of the system (22) in state space, using numerical solution (23), for the parameters: a = 1.73 , b = 2.04 , d = 0.51 , e = 0.04 , f = 0.67 , g = 0.19 , h = 0.48 , j = 0.52 , k = n 1 = n 2 = 0.21 , m 1 = 1 , m 2 = 0.001 , p = 3 d m 2 , r = 3 e m 2 , orders γ 1 = 0.99 (capacitor), γ 2 = 0.98 (memcapacitor), γ 3 = 0.96 (meminductor), γ 4 = 0.97 (memristor), γ 5 = γ 6 = 1 , initial conditions x ( 0 ) = 0.2 , y ( 0 ) = 0.5 , z ( 0 ) = 0.45 , u ( 0 ) = 0.3 , w ( 0 ) = 0.1 , v ( 0 ) = 0.5 , step h s = 0.005 , and simulation time 500 s. The Lyapunov exponents of this system with the above parameters, computed according to the algorithm described in [39], have the following values: 0.0020 , 0.1341 , 0.0002 , 0.7573 , 1.0016 , 0.0019 , which confirm the system is chaotic because at least one exponent is positive.

5. Discussion

Simulation results presented in this article confirm that new memory elements proposed only a few decades ago could be beneficial in new electronic circuit theory and practice. Moreover, a mathematical model for such elements should be a fractional order due to the memory of these parts. So far, it is possible only for a few well-known classical elements and the memristor [37]. Because the memory elements as memcapacitor and meminductor do not exist as a single electronic part, it is challenging to identify derivative order from experimental data. It will be possible when such elements are on sale. However, we may use an analogy with a classical capacitor and inductor to predict the real order in their mathematical models. Some experiments made by other authors confirm this theory [15,17].
Based on results depicted in Figure 5 and Figure 6, we can see that the new oscillator proposed in this article, depicted in Figure 4, may generate oscillation as well as chaotic behavior. Slight changes in parameters, orders, and initial conditions may produce different results. It could be significant in using such an oscillator with memory in secure communication, coding/decoding, and a new kind of memory circuits for future computers.
Here, we also present an idea for further work. Let us consider the circuit shown in Figure 7.
An important question is: Is it possible to generate chaos or hyperchaos using an electrical circuit, depicted in Figure 7, containing all known elements? It is an open problem.

6. Conclusions

This short article presents a new oscillator based on fractional-order elements with memory. We obtain a new circuit model described by the fractional differential equations.
Such a model was proposed as a numerical solution for simulation purposes as well as further investigation and analysis, for example, calculation of eigenvalues, Lyapunov exponents, Poincaré maps, and bifurcation diagrams. Moreover, the system (22) can work as a chaotic oscillator under appropriate parameters and initial values.
Such a combination of the fractional-order memory elements in the new oscillator circuit is presented for the first time. New is the part where an oscillator with three memory elements was proposed and its mathematical model and the numerical solution were derived. Both practical and theoretical relevance is significant for further circuit analysis.

Funding

This work was supported in part by the Slovak Grant Agency for Science under grant VEGA 1/0365/19 and by the Slovak Research and Development Agency under Contracts No. APVV-14-0892 and No. APVV-18-0526.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Yin, Z.; Tian, H.; Chen, G.; Chua, L.O. What are Memristor, Memcapacitor, and Meminductor? IEEE Trans. Circuits Syst. II Express Briefs 2015, 62, 402–406. [Google Scholar] [CrossRef] [Green Version]
  2. Chua, L.O. Memristor: The missing circuit element. IEEE Trans. Circuit Theory 1971, 18, 507–519. [Google Scholar] [CrossRef]
  3. Strukov, D.B.; Snider, G.S.; Stewart, D.R.; Williams, R.S. The missing memristor found. Nature 2008, 453, 80–83. [Google Scholar] [CrossRef] [PubMed]
  4. Chua, L.O.; Sung, M.K. Memristive devices and systems. Proc. IEEE 1976, 64, 209–223. [Google Scholar] [CrossRef]
  5. Ventra, M.D.; Pershin, Y.V.; Chua, L.O. Circuit elements with memory: Memristors, memcapacitors and meminductors. Proc. IEEE 2009, 97, 1717–1724. [Google Scholar] [CrossRef] [Green Version]
  6. Biolek, D.; Biolek, Z.; Biolkova, V. SPICE modeling of memristive, memcapacitative and meminductive systems. In Proceedings of the 2009 European Conference on Circuit Theory and Design, Antalya, Turkey, 23–27 August 2009; pp. 249–252. [Google Scholar] [CrossRef]
  7. Romero, F.J.; Ohata, A.; Toral-Lopez, A.; Godoy, A.; Morales, D.P.; Rodriguez, N. Memcapacitor and Meminductor Circuit Emulators: A Review. Electronics 2021, 10, 1225. [Google Scholar] [CrossRef]
  8. Khalil, N.A.; Fouda, M.E.; Said, L.A.; Radwan, A.G.; Soliman, A.M. A general emulator for fractional-order memristive elements with multiple pinched points and application. AEU—Int. J. Electron. Commun. 2020, 124, 153338. [Google Scholar] [CrossRef]
  9. Khalil, N.A.; Hezayyin, H.G.; Said, L.A.; Madian, A.H.; Radwan, A.G. Active emulation circuits of fractional-order memristive elements and its applications. AEU—Int. J. Electron. Commun. 2021, 138, 153855. [Google Scholar] [CrossRef]
  10. Cam Taskiran, Z.; Sagbas, M.; Ayten, U.; Sedef, H. A New Universal Mutator Circuit for Memcapacitor and Meminductor Elements. AEU—Int. J. Electron. Commun. 2020, 119, 153180. [Google Scholar] [CrossRef]
  11. Arena, P.; Caponetto, R.; Fortuna, L.; Porto, D. Nonlinear Noninteger Order Circuits and Systems: An Introduction; World Scientific: Singapore, 2000. [Google Scholar] [CrossRef]
  12. Elwakil, A.S. Fractional-order circuits and systems: An emerging interdisciplinary research area. IEEE Circuits Syst. Mag. 2010, 10, 40–50. [Google Scholar] [CrossRef]
  13. Schafer, I.; Kruger, K. Modelling of lossy coils using fractional derivatives. J. Phys. D Appl. Phys. 2008, 41, 045001. [Google Scholar] [CrossRef]
  14. Khalil, N.A.; Said, L.A.; Radwan, A.G.; Soliman, A.M. General fractional order mem-elements mutators. Microelectron. J. 2019, 90, 211–221. [Google Scholar] [CrossRef]
  15. Abdelouahab, M.S.; Lozi, R.; Chua, L.O. Memfractance: A Mathematical Paradigmfor Circuit Elements with Memory. Int. J. Bifurc. Chaos Appl. Sci. Eng. 2014, 24, 1430023–1430029. [Google Scholar] [CrossRef] [Green Version]
  16. Coopmans, C.; Petráš, I.; Chen, Y. Analogue Fractional-Order Generalized Memristive Devices. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, San Diego, CA, USA, 30 August–2 September 2009; Volume 4, pp. 1127–1136. [Google Scholar] [CrossRef]
  17. Guo, Z.; Si, G.; Diao, L.; Jia, L.; Zhang, Y. Generalized modeling of the fractional-order memcapacitor and its character analysis. Commun. Nonlinear Sci. Numer. Simul. 2018, 59, 177–189. [Google Scholar] [CrossRef]
  18. Yang, N.N.; Xu, C.; Wu, C.; Jia, R.; Liu, C. Fractional-order cubic nonlinear flux-controlled memristor: Theoretical analysis, numerical calculation and circuit simulation. Nonlinear Dyn. 2019, 97, 33–44. [Google Scholar] [CrossRef]
  19. Zhou, P.; Tang, J. Clarify the physical process for fractional dynamical systems. Nonlinear Dyn. 2020, 100, 2353–2364. [Google Scholar] [CrossRef]
  20. Tenreiro Machado, J.A.; Lopes, A.M. Multidimensional scaling locus of memristor and fractional order elements. J. Adv. Res. 2020, 25, 147–157. [Google Scholar] [CrossRef]
  21. Tenreiro Machado, J. Fractional generalization of memristor and higher order elements. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 264–275. [Google Scholar] [CrossRef] [Green Version]
  22. Petráš, I. Fractional-Order Nonlinear Systems; Springer: New York, NY, USA, 2011. [Google Scholar] [CrossRef] [Green Version]
  23. Caponetto, R.; Dongola, G.; Fortuna, L.; Petráš, I. Fractional Order Systems: Modeling and Control Applications; World Scientific: Singapore, 2010. [Google Scholar] [CrossRef]
  24. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  25. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA, 1974. [Google Scholar]
  26. Westerlund, S.; Ekstam, L. Capacitor theory. IEEE Trans. Dielectr. Electr. Insul. 1994, 1, 826–839. [Google Scholar] [CrossRef]
  27. Lopez, C.S.; Gomez, V.H.C.; Aguilar, M.A.C.; Lopez, F.E.M. PID controller design based on memductor. Int. J. Electron. Commun. 2019, 101, 9–14. [Google Scholar] [CrossRef]
  28. Pershin, Y.V.; Ventra, M.D. Memristive circuits simulate memcapacitors and meminductors. Electron. Lett. 2010, 46, 517–518. [Google Scholar] [CrossRef] [Green Version]
  29. Ma, X.; Mou, J.; Liu, J.; Ma, C.; Yang, F.; Zhao, X. A novel simple chaotic circuit based on memristor–memcapacitor. Nonlinear Dyn. 2020, 100, 2859–2876. [Google Scholar] [CrossRef]
  30. Innocenti, G.; Liu, X.; Bi, X.; Yan, H.; Mou, J. A Chaotic Oscillator Based on Meminductor, Memcapacitor, and Memristor. Complexity 2021, 2021, 7223557. [Google Scholar] [CrossRef]
  31. Wang, X.; Yu, J.; Jin, C.; Iu, H.H.C.; Yu, S. Chaotic oscillator based on memcapacitor and meminductor. Nonlinear Dyn. 2019, 96, 161–173. [Google Scholar] [CrossRef]
  32. Radwan, A.G.; Fouda, M.E. On the Mathematical Modeling of Memristor, Memcapacitor, and Meminductor; Springer: Berlin/Heidelberg, Germany, 2015. [Google Scholar]
  33. Petráš, I.; Chen, Y. Fractional-order circuit elements with memory. In Proceedings of the 13th International Carpathian Control Conference (ICCC), High Tatras, Slovakia, 28–31 May 2012; pp. 552–558. [Google Scholar] [CrossRef]
  34. Itoh, M.; Chua, L.O. Memristor oscillation. Int. J. Bifurcat. Chaos Appl. Sci. Eng. 2008, 18, 3183–3206. [Google Scholar] [CrossRef]
  35. Petráš, I. Fractional-order memristor-based Chua’s circuit. IEEE Trans. Circuits Syst. II Express Briefs 2010, 57, 975–979. [Google Scholar] [CrossRef]
  36. Westerlund, S. Dead Matter Has Memory! Causal Consulting: Kalmar, Sweden, 2002. [Google Scholar]
  37. Wang, S.F.; Ye, A. Dynamical Properties of Fractional-Order Memristor. Symmetry 2020, 12, 437. [Google Scholar] [CrossRef] [Green Version]
  38. Petráš, I. “Comments on “Chaotic oscillator based on memcapacitor and meminductor” (Nonlinear Dyn, DOI: 10.1007/s11071-019-04781-5). Nonlinear Dyn. 2020, 102, 2945–2950. [Google Scholar] [CrossRef]
  39. Wolf, A.; Swift, J.B.; Swinney, H.L.; Vastano, J.A. Determining Lyapunov exponents from a time series. Phys. D Nonlinear Phenom. 1985, 16, 285–317. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Connections between all known electrical elements (inspired by Refs. [6,33]).
Figure 1. Connections between all known electrical elements (inspired by Refs. [6,33]).
Fractalfract 06 00283 g001
Figure 2. Chaotic oscillator circuit based on two memory elements [31].
Figure 2. Chaotic oscillator circuit based on two memory elements [31].
Fractalfract 06 00283 g002
Figure 3. One-scroll attractors of the system (18) in state space (left) and state plane (right).
Figure 3. One-scroll attractors of the system (18) in state space (left) and state plane (right).
Fractalfract 06 00283 g003
Figure 4. Chaotic oscillator circuit based on three memory elements.
Figure 4. Chaotic oscillator circuit based on three memory elements.
Fractalfract 06 00283 g004
Figure 5. Orbits of the fractional-order oscillator (22) in state spaces x y z , u w v , respectively.
Figure 5. Orbits of the fractional-order oscillator (22) in state spaces x y z , u w v , respectively.
Fractalfract 06 00283 g005
Figure 6. One-scroll attractors of the system (22) in state space x y z (left) and u w v (right).
Figure 6. One-scroll attractors of the system (22) in state space x y z (left) and u w v (right).
Fractalfract 06 00283 g006
Figure 7. New oscillator circuit based on all known electrical elements.
Figure 7. New oscillator circuit based on all known electrical elements.
Fractalfract 06 00283 g007
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Petráš, I. Oscillators Based on Fractional-Order Memory Elements. Fractal Fract. 2022, 6, 283. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6060283

AMA Style

Petráš I. Oscillators Based on Fractional-Order Memory Elements. Fractal and Fractional. 2022; 6(6):283. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6060283

Chicago/Turabian Style

Petráš, Ivo. 2022. "Oscillators Based on Fractional-Order Memory Elements" Fractal and Fractional 6, no. 6: 283. https://0-doi-org.brum.beds.ac.uk/10.3390/fractalfract6060283

Article Metrics

Back to TopTop