# The Generation, Analysis, and Circuit Implementation of a New Memristor Based Chaotic System

Authors: Yuxia Li; Xia Huang; Mei Guo
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

We present a new memristor based chaotic circuit, which is generated by replacing the nonlinear resistor in Chua’s circuit with a flux-controlled memristor and a negative conductance. The dynamical behaviors are verified not only by computer simulations but also by Lyapunov exponents, bifurcation analysis, Poincaré mapping, power spectrum analysis, and laboratory experiments.


---

## Body

## 1. Introduction

The memristor was postulated in 1971 by circuit theorist Chua as a missing nonlinear two-terminal electrical component relating electric charge and magnetic flux linkage [1]. In 2008, a team at HP Labs introduced a model for a memristance function based on thin films of titanium dioxide [2]. From then on, memristor has attracted worldwide interest due to its promising applications in nanoelectronic memories, computer logic, and neuromorphic computer architectures [3].Chaos has great potential application in many engineering fields, such as image encryption and secure communications. As a nonlinear circuit element, the memristor is very suitable for the design of chaotic circuits and the generation of chaotic systems. Some memristor based chaotic circuits were presented by replacing the Chua diodes with memristors characterized by monotone-increasing piecewise-linear function [4–6]. Several smooth flux-controlled memristors characterized by cubic or piecewise-quadratic nonlinearities were developed in [7–9]. Moreover, attempts have been made on the generation of memristor based hyperchaotic circuits [10, 11].In this paper, a memristor based chaotic circuit is presented by replacing the nonlinear resistor in Chua’s circuit with a flux-controlled memristor which is characterized by a smooth cubic nonlinear resistor. The rest of this paper is organized as follows. The new chaotic circuit with memristor and the corresponding mathematical model are presented in Section2. The dynamical behaviors of the new system are analyzed in Section 3. Circuit for the implementation of this chaotic circuit is proposed in Section 4. Finally, some conclusions are drawn in Section 5.


## 2. Chaotic Circuit with Memristor

The two nonlinear functionsM(q) and W(ϕ), called the memristance and memductance [6], respectively, are defined by
(1)M(q)=dϕ(q)dq,(2)W(ϕ)=dq(ϕ)dϕ.The memristor used in this paper is a flux-controlled memristor, and we assume that the flux-controlled memristor is characterized by a smooth continuous cubic monotone-increasing nonlinearity as follows:(3)q(ϕ)=a0+a1ϕ+a2ϕ3,
where a1<0,  a2>0.In view of (2), the memductance is given by
(4)W(ϕ)=dq(ϕ)dϕ=a1+3a2ϕ2.
Thus, by replacing Chua’s diode in Chua’s circuit with an active two-terminal memristive circuit consisting of a negative conductance and a flux-controlled memristor characterized by (3), a smooth memristor oscillator can be designed.The memristor based chaotic circuit is shown in Figure1. This circuit consists of four basic circuit elements and has only one negative element in addition to the nonlinearity.Figure 1
Chua’s circuit with a flux-controlled memristor.Applying Kirchhoff’s voltage and current laws and component’s current-voltage relationship to the circuit in Figure1, we can obtain a system of differential equations governing the dynamical behavior of the circuit in Figure 1, described by
(5)C1dv1dt=v2-v1R+Gv1-W(ϕ)v1,C2dv2dt=v1-v2R-i,Ldidt=v2+ri,dϕdt=v1,
where the ϕ-q characteristic curve of the flux-controlled memristor is described by (3) and (4).Equation (5) can be transformed into the dimensionless form as
(6)dxdt=α(y-x+γx-W(w)x),dydt=x-y+z,dzdt=-βy+kz,dwdt=x
by letting x=v1,  y=v2,  z=-i,  w=θ,  α=1/C1,  β=1/L,  γ=G,  k=r/L,  R=1,  C2=1 and defining the nonlinear functions q(w) and W(w) as
(7)q(w)=a0+a1w+a2w3,W(w)=dq(w)dw=a1+3a2w2.Takeα=4,  β=12,  k=0.688,  γ=1.5,  a1=-0.67×10-3,  a2=-0.02×10-3, and k=0.688. For initial conditions (1,1,1,1), the system (6) is chaotic with Lyapunov exponents λ1=0.21383,  λ2=0,  λ3=-0.0002, and λ4=-1.5259,   and the Lyapunov dimension is dL=3.14; the chaotic attractors are shown in Figure 2.Phase portraits of the chaotic system (6).
(a)
x
-
y
(b)
x
-
z
(c)
y
-
z
(d)
x-w
(e)
y-w
(f)
z
-
w
(g)
x
-
y
-
z
(h)
x
-
y
-
w


## 3. Numerical Analysis

Similar to dynamical analysis of general chaotic circuits, by using the conventional dynamical analysis methods, such as bifurcation diagram, Lyapunov exponent spectra, the dynamical behaviors of the chaotic system (6) are studied when the system parameters are varied.


### 3.1. Bifurcation Analysis

Holding the values ofα=4,  β=12,  γ=1.5,  a1=-0.67×10-3,  a2=-0.02×10-3 fixed and changing the value of k, the spectra of Lyapunov exponents on the (k-λ) plane are obtained as shown in Figure 3, and the corresponding bifurcation diagram is given in Figure 4.Figure 3
The Lyapunov exponents in the (k-λ) plane.Figure 4
Bifurcation diagram for increasingk.It can be seen from Figures4 and 5 that the bifurcation diagram is consistent well with the Lyapunov exponent spectrum. From Figure 5, we observe that system (6) undergoes the following routes as k increases.(1)
When0.2<k≤0.321,  λ1=0,  λ2,λ3,λ4<0, system (6) is periodic (Figure 5(a)).(2)
When0.322<k≤0.508,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(b)).(3)
When0.509<k≤0.631, λ1>0,  λ2=0,  λ3,λ4<0, system (6) is chaotic (Figure 5(c)).(4)
When0.632<k≤0.638,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(d)).(5)
When0.639<k≤0.66,λ1>0,  λ2=0,  λ3,λ4<0, system (6) is chaotic (Figure 5(e)).(6)
When0.661<k≤0.682,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(f)).(7)
When0.683<k≤0.7,λ1>0,  λ2=0,  λ3,λ4<0,  system (6) is chaotic (Figure 5(g)).Figure 5
Phase portraits of system (6) in (i) y-w plane and (ii) z-w plane with different k. (a) Period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.569; (d) quasi-period orbit with k=0.635; (e) chaotic attractor with k=0.659; (f) quasi-period orbit with k=0.675; (g) chaotic attractor with k=0.688.


### 3.2. Poincaré Maps and Power Spectrum Analysis

Poincaré mapping is shown in Figure6, and system parameters are selected as α=4,  β=12, a1=-0.67×10-3,  a2=0.02×10-3,   and 0.2<k≤0.7.Poincaré mapping of system (6) with different k: (a) period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.688.
(a)
(b)
(c)The power spectrum of system (6) with different k is displayed in Figure 7.Power spectrum of system (6) with different k: (a) period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.688.
(a)
(b)
(c)


## 3.1. Bifurcation Analysis

Holding the values ofα=4,  β=12,  γ=1.5,  a1=-0.67×10-3,  a2=-0.02×10-3 fixed and changing the value of k, the spectra of Lyapunov exponents on the (k-λ) plane are obtained as shown in Figure 3, and the corresponding bifurcation diagram is given in Figure 4.Figure 3
The Lyapunov exponents in the (k-λ) plane.Figure 4
Bifurcation diagram for increasingk.It can be seen from Figures4 and 5 that the bifurcation diagram is consistent well with the Lyapunov exponent spectrum. From Figure 5, we observe that system (6) undergoes the following routes as k increases.(1)
When0.2<k≤0.321,  λ1=0,  λ2,λ3,λ4<0, system (6) is periodic (Figure 5(a)).(2)
When0.322<k≤0.508,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(b)).(3)
When0.509<k≤0.631, λ1>0,  λ2=0,  λ3,λ4<0, system (6) is chaotic (Figure 5(c)).(4)
When0.632<k≤0.638,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(d)).(5)
When0.639<k≤0.66,λ1>0,  λ2=0,  λ3,λ4<0, system (6) is chaotic (Figure 5(e)).(6)
When0.661<k≤0.682,  λ1=λ2=0,  λ3,λ4<0, system (6) is quasi-periodic (Figure 5(f)).(7)
When0.683<k≤0.7,λ1>0,  λ2=0,  λ3,λ4<0,  system (6) is chaotic (Figure 5(g)).Figure 5
Phase portraits of system (6) in (i) y-w plane and (ii) z-w plane with different k. (a) Period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.569; (d) quasi-period orbit with k=0.635; (e) chaotic attractor with k=0.659; (f) quasi-period orbit with k=0.675; (g) chaotic attractor with k=0.688.


## 3.2. Poincaré Maps and Power Spectrum Analysis

Poincaré mapping is shown in Figure6, and system parameters are selected as α=4,  β=12, a1=-0.67×10-3,  a2=0.02×10-3,   and 0.2<k≤0.7.Poincaré mapping of system (6) with different k: (a) period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.688.
(a)
(b)
(c)The power spectrum of system (6) with different k is displayed in Figure 7.Power spectrum of system (6) with different k: (a) period orbit with k=0.321; (b) quasi-period orbit with k=0.486; (c) chaotic attractor with k=0.688.
(a)
(b)
(c)


## 4. Circuit Design for Chaotic System

In this section, we present experimental confirmation of the above numerical results through Multisim modeling. The circuit is designed as in Figure8.Figure 8
The chaotic circuit.The two-terminal memristor consists of five Op-Amps TL084, two analog multipliers AD633, thirteen resistors, and three capacitances. The parameter values of the circuit areR1,R2,  R7,  R8,  R11,  R12=1kΩ; R3=12.5kΩ; R4=2.5kΩ; R5=36kΩ; R6=1.5kΩ; R9=3.1kΩ; R10=0.5kΩ; C1=10nF; C2=100nF,  C3=47nF. The power supplies are Vp=+15 V and Vn=-15 V. For different parameter R13, the attractors are observed as in Figure 9.Figure 9
Phase portraits of system (6) with different R13, measured from the experimental setup: (a) period orbit,  R13=0.305kΩ; (b) quasi-period orbit, R13=0.329kΩ; (c) chaotic attractor, R13=0.527kΩ; (d) quasi-period orbit, R13=0.612kΩ; (e) chaotic attractor, R13=0.639 kΩ; (f) quasi-period orbit, R13=0.667 kΩ; (g) chaotic attractor, R13=0.689kΩ.


## 5.  Conclusions

In this paper, we present a memristor chaotic circuit, which is derived from Chua’s oscillator by replacing Chua’s diode with a charge-controlled memristor. This chaotic circuit uses only four basic circuit elements and has only one negative element in addition to the nonlinearity. The resulting chaotic system is demonstrated by computer simulations and verified by Lyapunov exponents, bifurcation, poincaré mapping, and power spectrum analysis.



---

*Source: 398306-2013-11-13.xml*