# A Novel Voltage-Controlled Tri-Valued Memristor and Its Application in Chaotic System

Authors: Xiaoyuan Wang; Xue Zhang; Meng Gao
Journal: Complexity
 (2020)

---

## Abstract

Memristor is a kind of passive nonlinear element, which is widely used in nonlinear systems, especially chaotic systems, because of its nanometer size, nonvolatile property, and good nonlinear characteristics. Compared with general chaotic systems, chaotic systems based on memristors have richer dynamic characteristics. However, the current research mainly focuses on the binary and continuous chaotic systems based on memristors, and studies on the tri-valued and multi-valued memristor chaotic systems are relative scarce. For this reason, a mathematical model of tri-valued memristor is proposed, and the circuit characteristics of the model are studied. Furthermore, based on this model, a new chaotic system is designed and analyzed. This innovation enriches the types of chaotic systems and lays the foundation for the application of tri-valued and multi-valued memristors in nonlinear systems.


---

## Body

## 1. Introduction

In 1971, Chua postulated the concept of memristor according to the mathematical relationship between chargeq and flux φdφ=Mdq and called it the fourth basic circuit element [1]. In 2008, the Hewlett-Packard Lab realized a practical memristor device, which set off an upsurge of research studies on memristor models and their applications [2]. Memristors are widely used in microelectronics, neural network, nonvolatile storage, application and simulation of spontaneous behavior, hard switching, and dynamic storage [3–5]. One of the most typical applications of memristor is to construct a chaotic system. Compared with the general chaotic systems, chaotic oscillators constructed with memristors have more complex and special dynamic characteristics, wider range of parameters, and are extremely sensitive to the initial values. Especially, when appropriate parameters of the memristive chaotic system are given, hidden attractors and coexisting attractors will be behaved [6–8]. In recent years, chaotic oscillators based on memristors are widely used in confidential communication, file encryption, and artificial intelligence [9–11].At present, the research on memristive chaotic systems mainly focuses on the binary and continuous memristors. In 2015, Ma introduced memristor into a four-wing chaotic system and finally obtained line equilibriums by adding a cross product term [12]. Kengne introduced a novel memristor-based oscillator, which was obtained from Shinriki’s circuit by substituting the nonlinear positive conductance with a first-order memristive diode bridge [13]. In 2016, Wu used a memristor to replace the resistor in a parallel RC network and designed a simpler memristor-based Venturi oscillator [14]. By leading memristors into a multiwing chaotic system, Zhou obtained the multiwing hyperchaotic attractor and verified its dynamic characteristics by numerical simulation [15]. In 2017, Wang built a multiscroll chaotic system based on a multisegment memristor [16]. Hu designed two image encryption algorithms based on the chaotic sequences generated by a three-dimensional chaotic circuit based on memristor [17]. In 2018, Fonzin replaced the diode in the original TCMNL (Tamasevicius et al. (1997) oscillator) circuit with a memristor and implemented the hardware of the circuit [18]. In 2019, Min built a hyperchaotic system based on the memristor, which has symmetric coexisting attractors and exhibits complex system characteristics. The correctness of the system was verified by hardware circuit experiments [19]. In 2020, Wang proposed a new hyperchaotic circuit by introducing memristor feedback into a simple Lorenz-like chaotic system. By further analyses, infinite equilibrium points, multiple stability, and symmetric coexisting attractors are found [20].Nevertheless, compared with the continuous and binary memristors, tri-valued and multi-valued memristors have advantages of carrying more information. Constructing a chaotic system based on tri-valued memristors will generate new attractors, expand chaotic types, and broaden ideas for chaotic systems design. Consequently, it is of great significance and application value to propose a tri-valued memristor model and use the model to devise a chaotic system with good performance.This paper is organized as follows. A specific mathematical model of a voltage-controlled tri-valued memristor is proposed and studied in Section2. In Section 3, on the basis of the Lü system, a chaotic system with the proposed voltage-controlled tri-valued memristor is constructed. And, the basic characteristics of the system are analyzed in detail, including dissipative analysis, equilibrium point and stability analyses, and influences of system parameters and initial values on dynamic characteristics, etc. Conclusions are drawn in Section 4.


## 2. Voltage-Controlled Tri-Valued Memristor and its Characteristic Analysis

### 2.1. Mathematical Model of a Voltage-Controlled Tri-Valued Memristor

Different from binary and continuous memristors, this paper presents a mathematical model of a voltage-controlled tri-valued memristor, whoseq−φ relationship is given as follows:(1)q=−1.5+2.5φ+4φ+1−2.5φ−1.According to the mathematical definition of memristors and the derivative of equation (1), the relationship between the memconductance and flux of the tri-valued memristor is given by(2)dqdφ=Gφ=2.5+4sgnφ+1−2.5sgnφ−1=1,φ<−1,9,−1≤φ≤1,4,φ>1,where Gφ denotes memconductance, sgnx represents the symbolic function; when x>0, sgnx=1 and when x≤0, sgnx=−1. The equation (2) shows that the magnitude of the flux φ affects the resistance state of the tri-valued memristor, which means that when the flux φ takes different values, the memristor is in three different memconductances. The flux-charge relationship described by equation (1) and the flux- memconductance relationship described by equation (2) are shown in Figure 1.Figure 1
Characteristic curves of voltage-controlled tri-valued memristor model: (a)φ−q curve; (b) φ−G curve.
(a)(b)Figure1(a) indicates that the q−φ relationship of the proposed memristor is described by a three-segment linear curve that passes through the origin, and each slope of the curve corresponds to the instantaneous conductance value of the memristor. Figure 1(b) intuitively depicts the effect of flux φ on the memconductance of the proposed memristor.


### 2.2. Circuit Characteristics of the Voltage-Controlled Tri-Valued Memristor

To study the circuit characteristics of the voltage-controlled tri-valued memristor, a sinusoidal signalvt=v0sin2πft is applied to this model. According to the mathematical definition of the flux φ, the expression of φt can be shown as(3)φt−φ0=∫0tvτdτ=∫0tv0sin2πfτdτ=−v02πfcos2πft+v02πf.From equation (3), the φt can be written as(4)φt=φ0−v02πfcos2πft+v02πf.It can be seen from equation (4) that the magnitude of the flux φ is not only related to the initial value φ0 but also to the amplitude v0 and the frequency f of the input voltage, and the change interval of φt can be obtained as φ0,φ0+v0/πf. Then, combined with equation (2), we can conclude the specific relationship between φt and the resistance state of the tri-valued memristor in Table 1, which manifests the initial value φ0, the amplitude v0, and the frequency f of the input voltage are the main factors affecting the change of the resistance state of the proposed memristor. In Section 2.2.1, the influence of these three factors on the characteristics of the tri-valued memristor will be studied in detail.Table 1
Relationship between the number of resistance states of tri valued memristors andφt.
CasesValue range ofφtNumber of resistance statesCase1φ0<−1,φ0+v0/πf<−11Case2φ0<−1,−1<φ0+v0/πf<12Case3φ0<−1,φ0+v0/πf>13Case4−1<φ0<1,−1<φ0+v0/πf<11Case5−1<φ0<1,φ0+v0/πf>12Case6φ0>1,φ0+v0/πf>11


#### 2.2.1. Influence of Differentφ0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Table1 indicates that, when the input amplitude v0 and frequency f are fixed, the initial flux φ0 affects the resistance state of the tri-valued memristor. Let v0 = 2 V and f = 0.159 Hz, then φt belongs to the φ0,φ0+4 so when φ0 takes different values, the resistance states and hysteresis curves of the tri-valued memristor are different as shown in Figure 2, and these three curves correspond to Case3, Case5, and Case6 in Table 1. These results indicate that, when φ0 increases gradually, the tri-valued memristor changes from a tri-valued state to a binary state and finally to a single state.Figure 2
v-i curve of tri-valued memristor at different initial values: (a) φ0=−1.5 (Case3), (b) φ0=−0.5 (Case5), and (c) φ0=1.5 (Case6).
(a)(b)(c)


#### 2.2.2. Influence of Differentv0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, f = 0.159 Hz, then we can obtain the range of φt as [−1.5, 2v0 − 1.5]. Then, adjusting the input voltage amplitude v0 to 0.2 V, 1.2 V, and 3 V, which corresponds to Case1, Case2, and Case3 in Table 1 respectively, we can finally obtain the v-i hysteresis curves with different shapes as shown in Figure 3.Figure 3
v-i curve of the tri-valued memristor at different amplitudes: (a) v0 = 0.2 V, (b) v0 = 1.2 V, and (c) v0 = 3 V.
(a)(b)(c)The experimental results in Figure3 illustrate that, when φ0 and f are definite, as the increase of v0, the tri-valued memristor changes from single state to binary state and from binary state to tri-valued, the area of hysteresis curve also increases accordingly.


#### 2.2.3. Influence of Differentf Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, v0  = 2 V, then the range of φt is [−1.5, 0.637/f − 1.5]. When the input voltage frequency f is adjusted to 0.1 Hz, 1 Hz, and 1.5 Hz (corresponding to Case3, Case2, and Case1 in Table 1, respectively), hysteresis curves at different input frequencies can be obtained, as shown in Figure 4, which shows that, when the frequency of the input signal increases, the tri-valued memristor changes from tri-valued memristor state to binary memristor state and finally to single memristor state, which is consistent with the theoretical derivation in Table 1.Figure 4
v-i curve of the tri-valued memristor at different frequencies: (a) f = 0.1 Hz, (b) f = 1 Hz, and (c) f = 1.5 Hz.
(a)(b)(c)


#### 2.2.4. Influence of Different Input Signals on the Characteristics of Tri-Valued Memristor

The above three influencing factors are all discussed under the condition of sinusoidal signal as the input. In this section, the circuit characteristics under different kinds of input signals are studied. By applying sine, square, and triangular waves with amplitude of 4 V and frequency of 0.159 Hz to the tri-valued memristor model, respectively, and letφ0=−1.5, the hysteresis curves of the three input signals can be obtained, respectively, as shown in Figure 5.Figure 5
v-i curve of tri-valued memristor at different input signals: (a) sine wave, (b) square wave, and (c) triangular wave.
(a)(b)(c)From the above simulation results, it can be concluded that when the periodic signal of any zero DC component acts on the tri-valued memristor, its input and output responses can be represented as a hysteresis curve across the origin in thev-i plane, and all of these curves have three value characteristics.


## 2.1. Mathematical Model of a Voltage-Controlled Tri-Valued Memristor

Different from binary and continuous memristors, this paper presents a mathematical model of a voltage-controlled tri-valued memristor, whoseq−φ relationship is given as follows:(1)q=−1.5+2.5φ+4φ+1−2.5φ−1.According to the mathematical definition of memristors and the derivative of equation (1), the relationship between the memconductance and flux of the tri-valued memristor is given by(2)dqdφ=Gφ=2.5+4sgnφ+1−2.5sgnφ−1=1,φ<−1,9,−1≤φ≤1,4,φ>1,where Gφ denotes memconductance, sgnx represents the symbolic function; when x>0, sgnx=1 and when x≤0, sgnx=−1. The equation (2) shows that the magnitude of the flux φ affects the resistance state of the tri-valued memristor, which means that when the flux φ takes different values, the memristor is in three different memconductances. The flux-charge relationship described by equation (1) and the flux- memconductance relationship described by equation (2) are shown in Figure 1.Figure 1
Characteristic curves of voltage-controlled tri-valued memristor model: (a)φ−q curve; (b) φ−G curve.
(a)(b)Figure1(a) indicates that the q−φ relationship of the proposed memristor is described by a three-segment linear curve that passes through the origin, and each slope of the curve corresponds to the instantaneous conductance value of the memristor. Figure 1(b) intuitively depicts the effect of flux φ on the memconductance of the proposed memristor.


## 2.2. Circuit Characteristics of the Voltage-Controlled Tri-Valued Memristor

To study the circuit characteristics of the voltage-controlled tri-valued memristor, a sinusoidal signalvt=v0sin2πft is applied to this model. According to the mathematical definition of the flux φ, the expression of φt can be shown as(3)φt−φ0=∫0tvτdτ=∫0tv0sin2πfτdτ=−v02πfcos2πft+v02πf.From equation (3), the φt can be written as(4)φt=φ0−v02πfcos2πft+v02πf.It can be seen from equation (4) that the magnitude of the flux φ is not only related to the initial value φ0 but also to the amplitude v0 and the frequency f of the input voltage, and the change interval of φt can be obtained as φ0,φ0+v0/πf. Then, combined with equation (2), we can conclude the specific relationship between φt and the resistance state of the tri-valued memristor in Table 1, which manifests the initial value φ0, the amplitude v0, and the frequency f of the input voltage are the main factors affecting the change of the resistance state of the proposed memristor. In Section 2.2.1, the influence of these three factors on the characteristics of the tri-valued memristor will be studied in detail.Table 1
Relationship between the number of resistance states of tri valued memristors andφt.
CasesValue range ofφtNumber of resistance statesCase1φ0<−1,φ0+v0/πf<−11Case2φ0<−1,−1<φ0+v0/πf<12Case3φ0<−1,φ0+v0/πf>13Case4−1<φ0<1,−1<φ0+v0/πf<11Case5−1<φ0<1,φ0+v0/πf>12Case6φ0>1,φ0+v0/πf>11


### 2.2.1. Influence of Differentφ0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Table1 indicates that, when the input amplitude v0 and frequency f are fixed, the initial flux φ0 affects the resistance state of the tri-valued memristor. Let v0 = 2 V and f = 0.159 Hz, then φt belongs to the φ0,φ0+4 so when φ0 takes different values, the resistance states and hysteresis curves of the tri-valued memristor are different as shown in Figure 2, and these three curves correspond to Case3, Case5, and Case6 in Table 1. These results indicate that, when φ0 increases gradually, the tri-valued memristor changes from a tri-valued state to a binary state and finally to a single state.Figure 2
v-i curve of tri-valued memristor at different initial values: (a) φ0=−1.5 (Case3), (b) φ0=−0.5 (Case5), and (c) φ0=1.5 (Case6).
(a)(b)(c)


### 2.2.2. Influence of Differentv0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, f = 0.159 Hz, then we can obtain the range of φt as [−1.5, 2v0 − 1.5]. Then, adjusting the input voltage amplitude v0 to 0.2 V, 1.2 V, and 3 V, which corresponds to Case1, Case2, and Case3 in Table 1 respectively, we can finally obtain the v-i hysteresis curves with different shapes as shown in Figure 3.Figure 3
v-i curve of the tri-valued memristor at different amplitudes: (a) v0 = 0.2 V, (b) v0 = 1.2 V, and (c) v0 = 3 V.
(a)(b)(c)The experimental results in Figure3 illustrate that, when φ0 and f are definite, as the increase of v0, the tri-valued memristor changes from single state to binary state and from binary state to tri-valued, the area of hysteresis curve also increases accordingly.


### 2.2.3. Influence of Differentf Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, v0  = 2 V, then the range of φt is [−1.5, 0.637/f − 1.5]. When the input voltage frequency f is adjusted to 0.1 Hz, 1 Hz, and 1.5 Hz (corresponding to Case3, Case2, and Case1 in Table 1, respectively), hysteresis curves at different input frequencies can be obtained, as shown in Figure 4, which shows that, when the frequency of the input signal increases, the tri-valued memristor changes from tri-valued memristor state to binary memristor state and finally to single memristor state, which is consistent with the theoretical derivation in Table 1.Figure 4
v-i curve of the tri-valued memristor at different frequencies: (a) f = 0.1 Hz, (b) f = 1 Hz, and (c) f = 1.5 Hz.
(a)(b)(c)


### 2.2.4. Influence of Different Input Signals on the Characteristics of Tri-Valued Memristor

The above three influencing factors are all discussed under the condition of sinusoidal signal as the input. In this section, the circuit characteristics under different kinds of input signals are studied. By applying sine, square, and triangular waves with amplitude of 4 V and frequency of 0.159 Hz to the tri-valued memristor model, respectively, and letφ0=−1.5, the hysteresis curves of the three input signals can be obtained, respectively, as shown in Figure 5.Figure 5
v-i curve of tri-valued memristor at different input signals: (a) sine wave, (b) square wave, and (c) triangular wave.
(a)(b)(c)From the above simulation results, it can be concluded that when the periodic signal of any zero DC component acts on the tri-valued memristor, its input and output responses can be represented as a hysteresis curve across the origin in thev-i plane, and all of these curves have three value characteristics.


## 2.2.1. Influence of Differentφ0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Table1 indicates that, when the input amplitude v0 and frequency f are fixed, the initial flux φ0 affects the resistance state of the tri-valued memristor. Let v0 = 2 V and f = 0.159 Hz, then φt belongs to the φ0,φ0+4 so when φ0 takes different values, the resistance states and hysteresis curves of the tri-valued memristor are different as shown in Figure 2, and these three curves correspond to Case3, Case5, and Case6 in Table 1. These results indicate that, when φ0 increases gradually, the tri-valued memristor changes from a tri-valued state to a binary state and finally to a single state.Figure 2
v-i curve of tri-valued memristor at different initial values: (a) φ0=−1.5 (Case3), (b) φ0=−0.5 (Case5), and (c) φ0=1.5 (Case6).
(a)(b)(c)


## 2.2.2. Influence of Differentv0 Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, f = 0.159 Hz, then we can obtain the range of φt as [−1.5, 2v0 − 1.5]. Then, adjusting the input voltage amplitude v0 to 0.2 V, 1.2 V, and 3 V, which corresponds to Case1, Case2, and Case3 in Table 1 respectively, we can finally obtain the v-i hysteresis curves with different shapes as shown in Figure 3.Figure 3
v-i curve of the tri-valued memristor at different amplitudes: (a) v0 = 0.2 V, (b) v0 = 1.2 V, and (c) v0 = 3 V.
(a)(b)(c)The experimental results in Figure3 illustrate that, when φ0 and f are definite, as the increase of v0, the tri-valued memristor changes from single state to binary state and from binary state to tri-valued, the area of hysteresis curve also increases accordingly.


## 2.2.3. Influence of Differentf Values on the Characteristics of Voltage-Controlled Tri-Valued Memristor

Letφ0=−1.5, v0  = 2 V, then the range of φt is [−1.5, 0.637/f − 1.5]. When the input voltage frequency f is adjusted to 0.1 Hz, 1 Hz, and 1.5 Hz (corresponding to Case3, Case2, and Case1 in Table 1, respectively), hysteresis curves at different input frequencies can be obtained, as shown in Figure 4, which shows that, when the frequency of the input signal increases, the tri-valued memristor changes from tri-valued memristor state to binary memristor state and finally to single memristor state, which is consistent with the theoretical derivation in Table 1.Figure 4
v-i curve of the tri-valued memristor at different frequencies: (a) f = 0.1 Hz, (b) f = 1 Hz, and (c) f = 1.5 Hz.
(a)(b)(c)


## 2.2.4. Influence of Different Input Signals on the Characteristics of Tri-Valued Memristor

The above three influencing factors are all discussed under the condition of sinusoidal signal as the input. In this section, the circuit characteristics under different kinds of input signals are studied. By applying sine, square, and triangular waves with amplitude of 4 V and frequency of 0.159 Hz to the tri-valued memristor model, respectively, and letφ0=−1.5, the hysteresis curves of the three input signals can be obtained, respectively, as shown in Figure 5.Figure 5
v-i curve of tri-valued memristor at different input signals: (a) sine wave, (b) square wave, and (c) triangular wave.
(a)(b)(c)From the above simulation results, it can be concluded that when the periodic signal of any zero DC component acts on the tri-valued memristor, its input and output responses can be represented as a hysteresis curve across the origin in thev-i plane, and all of these curves have three value characteristics.


## 3. Chaotic System Based on a Voltage-Controlled Tri-Valued Memristor

In 2001, Professor Lü proposed the famous Lü system [21], which can realize the conversion between Lorenz system [22] and Chen system [23]. In this paper, by adding a voltage-controlled tri-valued memristor model to the Lü system, a chaotic system based on a tri-valued memristor is generated as follows:(5)x˙=ay−x,y˙=cy−xz,z˙=xy−bz−dGw,w˙=z,where a, b, c, and d are system parameters, all of them are real constants, and Gw represents the memconductance of the voltage-controlled tri-valued memristor shown in equation (2); here, the variable w is selected as a dimensionless mathematical representation of the flux φ. When the parameters are set as a = 40, b = 5, c = 24.4, and d = 50 and initial value x0,y0,z0,w0=0.01,0.01,0.01,0.01, the phase diagrams of the system are as shown in Figure 6. The corresponding Lyapunov exponent values are calculated as LE1 = 4.2486, LE2 = 0.0025, LE3 = −0.004, and LE4 = −24.8471 by the Jacobi method. Figure 7 shows the Poincare mappings obtained when the cross planes are selected as x = −10 and z = 50, which are dense points with hierarchical structures. All the above results indicate that the system can be working as a chaotic oscillator under appropriate parameters and initial values.Figure 6
Phase diagrams of the system based on the tri-valued memristor: (a)x-y, (b) x-z, (c) y-z, and (d) x-y-z.
(a)(b)(c)(d)Figure 7
Poincare sections of the system: (a)x = −10 and (b) z = 50.
(a)(b)


### 3.1. Dissipative Analysis

From the perspective of the dissipation of the chaotic system, to generate chaotic attractors, it is necessary for the system to be dissipative. So, we calculate the dissipativity of the system as follows:(6)∇V=∂x˙∂x+∂y˙∂y+∂z˙∂z+∂w˙∂w=−a+c−b.Settinga = 40, b = 5, c = 24.4, and d = 50, we can get ∇V=−20.6<0, implying that all trajectories are ultimately confined to a specific subset of zero volume and the system satisfies the conditions for chaos.


### 3.2. Equilibrium Point and Stability Analysis

Letx˙=y˙=z˙=w˙=0 in equation (5), because all the system parameters a, b, c, and d are not nonzero, so we can obtain x = 0, y = 0, z = 0, and Gw=0. But, as we know from the mathematical definition of the tri-valued memristor in equations (1) and (2), Gw is not equal to zero, so we can draw the conclusion that the chaotic system based on the novel tri-valued memristor has no equilibrium point, and the attractors generated by the system are hidden attractors.


### 3.3. Influence of System Parameters and Initial Values on Dynamic Characteristics

#### 3.3.1. Influence of Parameterc on System Dynamics

Givenx0,y0,z0,w0=0.01,0.01,0.01,0.01, a = 40, b = 5, and d = 50, the Lyapunov exponent spectrum and the bifurcation diagram of the system with respect to the parameter c can be obtained as shown in Figures 8 and 9. Table 2 presents the states of the system when the parameter c is in different intervals. And, Figure 10 shows the attractor phase diagrams of the system on the x-z plane when c takes different values.Figure 8
Lyapunov exponent spectrum corresponding toc.Figure 9
Bifurcation diagram corresponding toc.Table 2
Under different system parametersc, the system is in different states.
Range of system parametercSystem state[28.9, 29.2] [31, 33]Period-1 state[28.1, 28.5]Period-2 state[24, 27.5]Chaotic stateFigure 10
Influence of parameterc on the attractors of the x-z plane: (a) c = 29, period-1, (b) c = 28.1, period-2, and (c) c = 24.4, chaotic state.
(a)(b)(c)


#### 3.3.2. Influence of the Initial Value on the System

Chaos is extremely sensitive to initial values, and different initial values will eventually produce different trajectories. It is necessary to estimate the impact of initial conditions on system (5) for fixed sets of parameter values. The sensitivity of the sequences can be analyzed by measuring the correlation of the two sequences, which is defined by(7)Co=EXt−μXYt−μYσXσY,where Xt and Yt are two sequences generated by the system (5) with slight changes in its initial value, μ and σ are the mean value and standard deviation, and E[·] is the expectation function [24]. The closer the correlation value is to 0, the higher the sensitivity of the system is, and the greater the influence of the initial value on the system is.In this paper, we slightly change each variable in the initial valuex0,y0,z0,w0 with a 10−8 difference. Concretely taking the x0 case as an example, let x0′=x0+10−8, then we can obtain a different sequence pair X1,X2 under initial values x0,y0,z0,w0=0.01,0.01,0.01,0.01 and x0′,y0,z0,w0=0.01+10−8,0.01,0.01,0.01. In the same way, we can get the other sequence pairs Y1,Y2, Z1,Z2, and W1,W2 generated by applying a tiny change to y0, z0, and w0 in the initial value. Table 3 shows the correlation results in each case. As can be seen, the correlation values of each sequence pair are extremely close to 0. To further illustrate the system’s sensitivity to initial states, the timing diagrams of the sequence pairs X1,X2, Y1,Y2, and Z1,Z2 under x0−x0′ case are shown in Figure 11 as an example, which intuitively demonstrate the differences between the two output sequences generated by the system under tiny different initial conditions. So, we can conclude the system is extremely sensitive to the initial values.Table 3
Correlation values of different sequences under tiny different initial values.
Initial valuesCorrelation ofX1, X2Correlation ofY1, Y2Correlation ofZ1, Z2x0−x0′ case−0.0743−0.0220−0.0348y0−y0′ case0.0084−0.0002−0.0799z0−z0′ case0.09090.04750.0731w0−w0′ case0.12780.1042−0.1441Figure 11
The timing diagrams of different sequence pairs of underx0,y0,z0,w0 and x0′,y0,z0,w0: (a) X1,X2, (b) Y1,Y2, and (c) Z1,Z2.
(a)(b)(c)


## 3.1. Dissipative Analysis

From the perspective of the dissipation of the chaotic system, to generate chaotic attractors, it is necessary for the system to be dissipative. So, we calculate the dissipativity of the system as follows:(6)∇V=∂x˙∂x+∂y˙∂y+∂z˙∂z+∂w˙∂w=−a+c−b.Settinga = 40, b = 5, c = 24.4, and d = 50, we can get ∇V=−20.6<0, implying that all trajectories are ultimately confined to a specific subset of zero volume and the system satisfies the conditions for chaos.


## 3.2. Equilibrium Point and Stability Analysis

Letx˙=y˙=z˙=w˙=0 in equation (5), because all the system parameters a, b, c, and d are not nonzero, so we can obtain x = 0, y = 0, z = 0, and Gw=0. But, as we know from the mathematical definition of the tri-valued memristor in equations (1) and (2), Gw is not equal to zero, so we can draw the conclusion that the chaotic system based on the novel tri-valued memristor has no equilibrium point, and the attractors generated by the system are hidden attractors.


## 3.3. Influence of System Parameters and Initial Values on Dynamic Characteristics

### 3.3.1. Influence of Parameterc on System Dynamics

Givenx0,y0,z0,w0=0.01,0.01,0.01,0.01, a = 40, b = 5, and d = 50, the Lyapunov exponent spectrum and the bifurcation diagram of the system with respect to the parameter c can be obtained as shown in Figures 8 and 9. Table 2 presents the states of the system when the parameter c is in different intervals. And, Figure 10 shows the attractor phase diagrams of the system on the x-z plane when c takes different values.Figure 8
Lyapunov exponent spectrum corresponding toc.Figure 9
Bifurcation diagram corresponding toc.Table 2
Under different system parametersc, the system is in different states.
Range of system parametercSystem state[28.9, 29.2] [31, 33]Period-1 state[28.1, 28.5]Period-2 state[24, 27.5]Chaotic stateFigure 10
Influence of parameterc on the attractors of the x-z plane: (a) c = 29, period-1, (b) c = 28.1, period-2, and (c) c = 24.4, chaotic state.
(a)(b)(c)


### 3.3.2. Influence of the Initial Value on the System

Chaos is extremely sensitive to initial values, and different initial values will eventually produce different trajectories. It is necessary to estimate the impact of initial conditions on system (5) for fixed sets of parameter values. The sensitivity of the sequences can be analyzed by measuring the correlation of the two sequences, which is defined by(7)Co=EXt−μXYt−μYσXσY,where Xt and Yt are two sequences generated by the system (5) with slight changes in its initial value, μ and σ are the mean value and standard deviation, and E[·] is the expectation function [24]. The closer the correlation value is to 0, the higher the sensitivity of the system is, and the greater the influence of the initial value on the system is.In this paper, we slightly change each variable in the initial valuex0,y0,z0,w0 with a 10−8 difference. Concretely taking the x0 case as an example, let x0′=x0+10−8, then we can obtain a different sequence pair X1,X2 under initial values x0,y0,z0,w0=0.01,0.01,0.01,0.01 and x0′,y0,z0,w0=0.01+10−8,0.01,0.01,0.01. In the same way, we can get the other sequence pairs Y1,Y2, Z1,Z2, and W1,W2 generated by applying a tiny change to y0, z0, and w0 in the initial value. Table 3 shows the correlation results in each case. As can be seen, the correlation values of each sequence pair are extremely close to 0. To further illustrate the system’s sensitivity to initial states, the timing diagrams of the sequence pairs X1,X2, Y1,Y2, and Z1,Z2 under x0−x0′ case are shown in Figure 11 as an example, which intuitively demonstrate the differences between the two output sequences generated by the system under tiny different initial conditions. So, we can conclude the system is extremely sensitive to the initial values.Table 3
Correlation values of different sequences under tiny different initial values.
Initial valuesCorrelation ofX1, X2Correlation ofY1, Y2Correlation ofZ1, Z2x0−x0′ case−0.0743−0.0220−0.0348y0−y0′ case0.0084−0.0002−0.0799z0−z0′ case0.09090.04750.0731w0−w0′ case0.12780.1042−0.1441Figure 11
The timing diagrams of different sequence pairs of underx0,y0,z0,w0 and x0′,y0,z0,w0: (a) X1,X2, (b) Y1,Y2, and (c) Z1,Z2.
(a)(b)(c)


## 3.3.1. Influence of Parameterc on System Dynamics

Givenx0,y0,z0,w0=0.01,0.01,0.01,0.01, a = 40, b = 5, and d = 50, the Lyapunov exponent spectrum and the bifurcation diagram of the system with respect to the parameter c can be obtained as shown in Figures 8 and 9. Table 2 presents the states of the system when the parameter c is in different intervals. And, Figure 10 shows the attractor phase diagrams of the system on the x-z plane when c takes different values.Figure 8
Lyapunov exponent spectrum corresponding toc.Figure 9
Bifurcation diagram corresponding toc.Table 2
Under different system parametersc, the system is in different states.
Range of system parametercSystem state[28.9, 29.2] [31, 33]Period-1 state[28.1, 28.5]Period-2 state[24, 27.5]Chaotic stateFigure 10
Influence of parameterc on the attractors of the x-z plane: (a) c = 29, period-1, (b) c = 28.1, period-2, and (c) c = 24.4, chaotic state.
(a)(b)(c)


## 3.3.2. Influence of the Initial Value on the System

Chaos is extremely sensitive to initial values, and different initial values will eventually produce different trajectories. It is necessary to estimate the impact of initial conditions on system (5) for fixed sets of parameter values. The sensitivity of the sequences can be analyzed by measuring the correlation of the two sequences, which is defined by(7)Co=EXt−μXYt−μYσXσY,where Xt and Yt are two sequences generated by the system (5) with slight changes in its initial value, μ and σ are the mean value and standard deviation, and E[·] is the expectation function [24]. The closer the correlation value is to 0, the higher the sensitivity of the system is, and the greater the influence of the initial value on the system is.In this paper, we slightly change each variable in the initial valuex0,y0,z0,w0 with a 10−8 difference. Concretely taking the x0 case as an example, let x0′=x0+10−8, then we can obtain a different sequence pair X1,X2 under initial values x0,y0,z0,w0=0.01,0.01,0.01,0.01 and x0′,y0,z0,w0=0.01+10−8,0.01,0.01,0.01. In the same way, we can get the other sequence pairs Y1,Y2, Z1,Z2, and W1,W2 generated by applying a tiny change to y0, z0, and w0 in the initial value. Table 3 shows the correlation results in each case. As can be seen, the correlation values of each sequence pair are extremely close to 0. To further illustrate the system’s sensitivity to initial states, the timing diagrams of the sequence pairs X1,X2, Y1,Y2, and Z1,Z2 under x0−x0′ case are shown in Figure 11 as an example, which intuitively demonstrate the differences between the two output sequences generated by the system under tiny different initial conditions. So, we can conclude the system is extremely sensitive to the initial values.Table 3
Correlation values of different sequences under tiny different initial values.
Initial valuesCorrelation ofX1, X2Correlation ofY1, Y2Correlation ofZ1, Z2x0−x0′ case−0.0743−0.0220−0.0348y0−y0′ case0.0084−0.0002−0.0799z0−z0′ case0.09090.04750.0731w0−w0′ case0.12780.1042−0.1441Figure 11
The timing diagrams of different sequence pairs of underx0,y0,z0,w0 and x0′,y0,z0,w0: (a) X1,X2, (b) Y1,Y2, and (c) Z1,Z2.
(a)(b)(c)


## 4. Conclusion

In this paper, a voltage-controlled tri-valued memristor is proposed for the first time. To demonstrate its unique characteristics, the circuit characteristics, parameters properties, and influence factors on its properties have been studied in detail. And, then a novel chaotic system is successfully built based on the Lü system; by introducing the tri-valued memristor to the system, the maximum Lyapunov exponent has been improved. In addition, some conventional analyses like the influences of parameters and initial values on the system are considered. This study shows that the tri-valued memristor is suitable for building chaotic systems, which will enrich the types of nonlinear system, widen the application of tri-valued memristor, and lay a foundation for the subsequent application of tri-valued and multi-valued memristors.



---

*Source: 6949703-2020-07-25.xml*