Characteristic Analysis and Circuit Implementation of a Novel Fractional-Order Memristor-Based Clamping Voltage Drift

Abstract

The ideal magnetic flux-controlled memristor was introduced into a four-dimensional chaotic system and combined with fractional calculus theory, and a novel four-dimensional commensurate fractional-order system was proposed and solved using the Adomian decomposition method. The system orders, parameters, and initial values were studied as independent variables in the bifurcation diagram and Lyapunov exponents spectrum, and it was discovered that changing these variables can cause the system to exhibit more complex and rich dynamical behaviors. The system had an offset boosting, which was discovered by adding a constant term after the decoupled linear term. Finally, the results of the numerical simulation were verified through the use of analog circuits and FPGA designs, and a control scheme for the system circuit was also suggested.

1. Introduction

The first chaotic system was accidentally proposed by Lorenz [1], and scholars began to investigate various chaotic systems with different types of attractors [2,3,4,5,6,7] and equilibria [8,9,10,11]. The exceptional qualities of chaotic systems, such as their initial sensitivity and long-term unpredictability, allow for their cross-combination with other scientific disciplines, including image encryption [12], finance [13], and so forth. It has therefore become a priority to figure out how to establish new chaotic systems with improved physical properties, and the establishment of memristive chaotic system is an effective method [14]. Leon Chua once proposed the fourth fundamental element from a symmetry arguments memristor, which serves as the link between magnetic flux and charge [15]. In recent years, memristors have received much attention due to their unique non-volatile performance and variety of nonlinear characteristics, such as their Lissajous curve’s variation with frequency and the appearance of the voltage–current’s hysteresis curve shrinking at the origin [16]. By introducing memristors with different properties into the existing chaotic circuit systems, scholars have obtained a new class of memristive chaotic circuit systems, studied their abundantly rich dynamical behaviors [17,18,19,20,21], and used them to establish nonlinear circuits [22], neural networks [23,24], secure communications [25], and other engineering fields. However, most of these studies are based on the integer-order calculus system. Fractional-order calculus can accurately describe the physical model as a generalized description. Therefore, it has drawn the attention of scholars and has become the focus of nonlinear theory and engineering application studies.

The theoretical development of fractional-order calculus has been very slow due to a lack of practical engineering experience and the substantial amount of calculations using traditional pens. However, with the advancement of computer technology in recent years and the inclusion of variable differences in the range of $(0,1]$, numerous fields increasingly focus on fractional-calculus operators. Meanwhile, the combination of fractional-order calculus and memristors expands the influence of control parameters on the system performance and exhibits richer dynamical behaviors, including high-dimensional [26] and variable-order memristive chaotic systems [27]. With an increased study depth, scholars have begun to shift their original focus from the dynamical behaviors of chaotic memristive circuits with different parameters [28,29,30] to the study of memristive chaotic systems with changing initial values [31,32]. As a result, we can investigate a novel fractional-order memristor chaotic system and, by altering the initial value, can make the system exhibit extreme multistability. This can improve the complexity of chaotic pseudo-random sequences and better satisfy the requirements of code division multiple access and the spread spectrum in secure communication [33,34,35].

Meanwhile, a number of non-ideal issues with memristor operation have emerged in recent years along with the development of memristors, including random variability, voltage-drift and the aging-induced degradation of the memristor threshold voltage [36,37]. When we designed the memristive chaotic circuit, the integrator was made by the operational amplifier, and the positive input was grounded. The integration function was obtained by using the virtual short concept of the two inputs of the ideal operational amplifier, where the negative terminal is also zero. However, the input of the operational amplifier exists as offset; that is, when the positive input terminal is grounded, the negative input terminal is not zero, the input offset is also integrated, and the output terminal has voltage. Thus, the DC voltage integration drift occurs. Consequently, the current study will concentrate on how to design memristive chaotic circuits to clamp voltage integration drift. The paper is organized as follows: In Section 2, the mathematical expression of a fractional memristor is derived through fractional calculus theory. The performance of the memristor is studied through numerical research and circuit implementation. In Section 3, when orders, system parameters, and system initial values are selected as variables, the bifurcation diagram and Lyapunov exponents (LEs) spectrum are used to illustrate the characteristics of the system and then the offset boosting phenomenon of the system is researched. In Section 4, the analog circuit and FPGA design of the multi-type chaotic attractor are implemented, and the control circuit of the system is designed through the stability analysis of the equilibrium to satisfy the control requirements. The last section summarizes the work of this paper and gives prospects for future work.

2. Fractional-Order Chaotic Memristor Circuit and Model

2.1. Fractional-Order Calculus

The operator ${}_0^C{D}_t^q$ is generally used to represent fractional-order calculus. D is a differential operator, d stands for the differential definition, $\alpha$ and t are the boundaries of the operator, and q is a real number. When q is positive, it is a differentiation, and when q is negative, it is an integration. Different from integer-order calculus, there are many definitions of fractional-order calculus, among which, the Caputo definition [38] is more common, which is as follows:

$${}_0^C{D}_{t_0}^{q}f\left(t\right)=\left\{\begin{array}{c}{\left[\mathsf{\Gamma }\left(k-q\right)\right]}^{-1}{\int }_{t_0}^t\frac{f^{\left(k\right)}\left(\tau \right)}{{(t-\tau )}^{q-k+1}}\mathrm{d}\tau ,k-1<q<k,\hfill \\ \frac{{\mathrm{d}}^{k}f\left(t\right)}{\mathrm{d}{t}^k},q=k.\hfill \end{array}\right.$$

(1)

where $k\in N$, $q\in R^{+}$, $\Gamma \left(\mathcal{L}\right)={\int }_0^{+\infty }{\mathrm{e}}^{-t}{t}^{\mathcal{L}-1}dt$ is a gamma function, and C is the Caputo differential definition. When $t\in \left(t_0,t_1\right)$, $k-1<q<k$, the properties of Caputo fractional-order calculus are

$${}_0^{C}J{}_0^{C}Df\left(t\right)=f\left(t\right)-\sum _{k=0}^{k-1}{f}^{\left(k\right)}\left(t_0\right)\frac{{\left(t-t_0\right)}^k}{k!}.$$

(2)

$J_{t_0}^q$ is the q-order integral operator, $D_{t_0}^q$ is the q-order differential operator, and $D_{t_0}^{-q}=J_{t_0}^q$. When $q=0$, $D_{t_0}^{q}f\left(t\right)=J_{t_0}^{q}f\left(t\right)=f\left(t\right)$. In practical application, the definition of the fractional-order derivative of time is often defined by Caputo. From [39], we can obtain ${}_{-\infty }^{\mathrm{C}}{D}_t^q{\mathrm{e}}^{\lambda t}={\lambda }^q{\mathrm{e}}^{\lambda t}$; then, according to Euler’s equation ${\mathrm{e}}^{\mathrm{i}t}=cost+\mathrm{i}·sint$, we can expand it as follows:

$${}_{-\infty }^{\mathrm{C}}{D}_t^q{\mathrm{e}}^{\mathrm{i}\omega t}={\left(\mathrm{i}\omega \right)}^q{\mathrm{e}}^{\mathrm{i}\omega t}={\omega }^q{\mathrm{e}}^{\mathrm{i}\omega t+\mathrm{i}\frac{q\pi }{2}}={\omega }^{q}cos\left(\omega t+\frac{q\pi }{2}\right)+\mathrm{i}{\omega }^{q}sin\left(\omega t+\frac{q\pi }{2}\right).$$

(3)

Thus, Equation (3) can be

$$\left\{\begin{array}{c}{}_{-\infty }^{\mathrm{C}}{D}_t^{q}sin\left(\omega t\right)={\omega }^{q}sin\left(\omega t+\frac{q}{2}\pi \right),\\ {}_{-\infty }^{\mathrm{C}}{D}_t^{q}cos\left(\omega t\right)={\omega }^{q}cos\left(\omega t+\frac{q}{2}\pi \right).\end{array}\right.$$

(4)

and, for Equation (4) above, the fractional-order derivative formulae of sine and cosine functions are obtained by expanding the corresponding real and imaginary parts of Euler’s Equation (3). Then, from the short-term memory principle [39], when $t\gg 1$, Equation (4) can be approximated as

$$\left\{\begin{array}{c}{}_0^{\mathrm{C}}{D}_t^{q}sin\left(\omega t\right){\approx }_{-\infty }^{\mathrm{C}}{D}_t^{q}sin\left(\omega t\right)={\omega }^{q}sin\left(\omega t+\frac{q}{2}\pi \right),\hfill \\ {}_0^{\mathrm{C}}{D}_t^{q}cos\left(\omega t\right){\approx }_{-\infty }^{\mathrm{C}}{D}_t^{q}cos\left(\omega t\right)={\omega }^{q}cos\left(\omega t+\frac{q}{2}\pi \right).\hfill \end{array}\right.$$

(5)

the corresponding fractional-order integral is as follows:

$${}_0^C{J}_t^{q}f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(q\right)}{\int }_0^t\frac{f\left(\tau \right)}{{(t-\tau )}^{1-q}}\mathrm{d}\tau .$$

(6)

Therefore, when $q\in (0,1)$, we can obtain

$${}_0^C{D}_t^{1-q}\left({}_0^C{D}_t^{q}f\left(t\right)\right)=\dot{f}\left(t\right).$$

(7)

the subsequent derivation of the memristor model will make extensive use of many of the equations mentioned above that are related to this section.

2.2. Fractional-Order Memristor Model and Circuit

2.2.1. Magnetic Flux-Controlled Memristor Model

The model of the nonlinear fractional-order quadratic memristor is

$$\left\{\begin{array}{c}i\left(t\right)=W\left(\varphi \right)v\left(t\right),\\ {}_0^c{D}_t^{q}f\left(t\right)=v\left(t\right).\end{array}\right.$$

(8)

where $q>0$, $i\left(t\right)$ is the input current, $v\left(t\right)$ is the output voltage, $W\left(\varphi \right)$ is the memristor value, $\zeta$ and $\xi$ are parameters, $f\left(t\right)$ is the memristor internal variable, and $W\left(\varphi \right)=\zeta +\xi f{\left(t\right)}^2$. The driving voltage on both sides of the memristor is $v\left(t\right)=V_{m}sin\left(\omega t\right)$, and through Equations (7) and (8) simultaneously, we have

$$\dot{f}\left(t\right)={V_m}_0^C{D}_t^{1-q}sin\left(\omega t\right).$$

(9)

where $V_m$ is the amplitude and $\omega$ is the angular frequency. Then, from the short-term memory principle [39], Equation (9) can be approximated as

$$\dot{f}\left(t\right)\approx V_m{\omega }^{1-q}sin\left(\omega t+\frac{1-q}{2}\pi \right).$$

(10)

and by integrating the two sides of Equation (10), we obtain

$$f\left(t\right)\approx f\left(0\right)+\frac{V_m}{{\omega }^q}\left[cos\left(\frac{1-q}{2}\pi \right)-cos\left(\omega t+\frac{1-q}{2}\pi \right)\right].$$

(11)

When $q=1$, the above expression is simplified to an integer-order memristor equation and, when $0<q<1$, Equations (11) and (8) simultaneously become

$$\begin{array}{c}i\left(t\right)=\left(\zeta +\xi f{\left(t\right)}^2\right)v\left(t\right)\hfill \\ =\left(\zeta +\xi {\left(f\left(0\right)+\frac{V_m}{{\omega }^q}\left[cos\left(\frac{1-q}{2}\pi \right)-cos\left(\omega t+\frac{1-q}{2}\pi \right)\right]\right)}^2\right)V_{m}sin\left(\omega t\right).\hfill \end{array}$$

(12)

According to Equation (12), the system parameters of the memristor, the initial magnetic flux of the memristor, and the amplitude and frequency of the sinusoidal voltage excitation are all related to the current flowing through the quadratic nonlinear magnetic flux-controlled memristor. When the input voltage is $v\left(t\right)=V_{m}sin\left(\omega t\right)$, for different amplitudes and frequencies of sinusoidal voltage excitation and different internal initial conditions of the memristor, the quadratic nonlinear flux memristor hysteresis loop on the $v-i$ plane is shown in Figure 1a–c, and Figure 1d–f displays the memristor’s time-domain waveform.

From above, Figure 1a shows the alter of the area of the hysteresis loop when the order q of Equation (12) is changed; Figure 1b illustrates how the area of the hysteresis loop changes when the frequency is altered, and, when the increment is large enough, it becomes a single-valued function; Figure 1c shows that there is a positive correlation between the variation in the excitation voltage amplitude $V_m$ at the memristor terminal and the area of the hysteresis loop. Then, in the subsequent group, Figure 1d is the effect of order on the memristor, and the memristance $W\left(\varphi \right)$ will be positive or negative. When the fractional-order memristor degenerates to integer-order, it is found that the integer-order memristance is always negative; Figure 1e demonstrates that the dynamic range of the current memristor variation increases with a decreasing order; Figure 1f is the time-domain waveform of the memristor output power, and when the order changes, the memristor power changes between positive and negative. When the fractional-order memristor degenerates into an integer-order, the integer-order memristor power is always less than zero.

2.2.2. Implementation of Chaotic Circuit

According to the nonlinear magnetic flux-controlled memristor model described in Equation (12), parameters were set and the equivalent circuit of the fractional-order memristor was established, as shown in Figure 2, where, different from the traditional active magnetic flux-controlled memristor, a shunt resistor $R_2$ was added to the integrating capacitor to clamp the DC voltage integral drift. Furthermore, the fractional-order memristor model is

$$\left\{\begin{array}{c}i=\left(-\frac{R_5}{R_4}\frac{R_7}{R_6}\frac{R_9}{R_8}+\frac{R_9}{R_3}{\varphi }^2\right),\hfill \\ {}_0{D}_t^q\phi =-\frac{v}{R_1{C}_0}-\frac{\varphi }{R_2{C}_0}.\hfill \end{array}\right.$$

(13)

where v and i represent the input voltage and current, respectively, and $\varphi$ is the voltage at both ends of the integrating capacitor $C_0$. The circuit unit values involved in the equation are shown in

Table 1. Value of circuit elements of fractional-order memristor.

Circuit Element	Physical Meaning of Parameters	Parameter Value
$R_1,R_2,R_4,R_5,R_8,R_9$	resistance	10k
$R_3$	resistance	62 k
$R_6$	resistance	100 k
$R_7$	resistance	125 k
$C_0$	capacitance	1 nF

.

Following that, different excitation amplitudes and frequencies of a sinusoidal voltage source can be chosen in accordance with the parameters in Table 1 to simulate various hysteresis loops shrinking at the origin in Multism, as shown in Figure 3.

Figure 3a shows that the voltage excited by sinusoidal voltages has not changed, but the frequencies have changed: this is a “hard” switching condition (large voltage excursions or long-term bias), and any symmetrical AC voltage bias will cause double-loop $v-i$ hysteresis and shrink to a line at high frequency (blue line) as long as it remains in the memristor mechanism. Figure 3b shows that the frequencies have not changed and that the voltages have changed, and as the voltage decreases, the hysteresis loop of the memristor tends to be linear.

2.2.3. Dynamical Chaotic Circuit Model

According to Figure 2, the dynamical elements, capacitance $C_1$ and $C_2$, inductance L, and magnetic flux-controlled memristor $W\left(\varphi \right)$ correspond to state variables $v_1$, $v_2$, $i_L$, and $\varphi$, respectively. Therefore, the equations of state can be described by Kirchhoff voltage law as

$$\left\{\begin{array}{c}\frac{{\mathrm{d}}^q{v}_1}{\mathrm{d}t}=\frac{v_2-v_1}{R{C}_1}-\frac{W\left(\varphi \right)v_1}{C_1},\hfill \\ \frac{{\mathrm{d}}^q{v}_2}{\mathrm{d}t}=\frac{v_1-v_2}{R{C}_2}-\frac{i_L}{C_2},\hfill \\ \frac{{\mathrm{d}}^q{i}_L}{\mathrm{d}t}=\frac{v_2}{L},\hfill \\ \frac{{\mathrm{d}}^q\varphi }{\mathrm{d}t}=-\frac{v_1}{R_1{C}_0}-\frac{\varphi }{R_2{C}_0}.\hfill \end{array}\right.$$

(14)

where $q>0$. Let $x_1=v_1$, $x_2=v_2$, $x_3=R{i}_L$, and $x_4=\varphi$. Define nonlinear function $W\left(\varphi \right)=-\alpha +\beta {\varphi }^2$, and then let $a=\frac{\alpha -1}{R{C}_1}$, $b=\frac{1}{R{C}_1}$, $c=\frac{\beta }{R{C}_1}$, $R=1$, $R{C}_2=1$, $\frac{1}{L}=d$, $-\frac{1}{R_1{C}_0}=e$, $\frac{1}{R_2{C}_0}=f$, $\zeta =1.25$, and $\xi =0.16$. The normalized representation of Equation (14) can be rewritten as

$$\left\{\begin{array}{c}{D}_t^q{x}_1=a{x}_1+b{x}_2-c{x}_1{x_4}^2,\hfill \\ D_t^q{x}_2=x_1-x_2-x_3,\hfill \\ D_t^q{x}_3=d{x}_2,\hfill \\ D_t^q{x}_4=e{x}_1-f{x}_4.\hfill \end{array}\right.$$

(15)

where q represents the fractional order of the commensurate system, and a, b, c, d, e, and f are the system parameters.

2.3. Adomian Decomposition Algorithm and Solution

2.3.1. Adomian Decomposition Algorithm

The Adomian decomposition algorithm is a numerical analysis algorithm with a high accuracy and fast convergence. Its general method is to decompose the differential equation into nonlinear terms, linear terms, and constant terms, convert the nonlinear terms into a special polynomial, and then use the inverse operator method to derive step by step. Finally, the approximate solution of the differential equation is obtained. It is the best candidate method for solving fractional-order chaotic systems $D_{t_0}^{q}x\left(t\right)=f\left(x\left(t\right)\right)$ at present. $x\left(t\right)$ is the function variable, where $x\left(t\right)={\left(x_1\left(t\right),x_2\left(t\right),x_3\left(t\right)\cdots x_n\left(t\right)\right)}^{\mathrm{T}}$, $D_{t_0}^q$ represents the Caputo differential operator of q-order, whereupon its expression, based on the Adomian decomposition algorithm, is as follows:

$$\left\{\begin{array}{c}{D}_{t_0}^{q}x\left(t\right)=Lx\left(t\right)+Nx\left(t\right)+o\left(x\left(t\right)\right),\hfill \\ x^{\left(k\right)}\left(t_0^{+}\right)=b_k.\hfill \end{array}\right.$$

(16)

where, $m\in N$, $k=0,1,2,\cdots ,m-1$ and $m-1<q⩽m$, the linear operator of the system can be represented by L and the nonlinear operator can be represented by N, the $o\left(t\right)={\left[o_1\left(t\right),o_2\left(t\right),\cdots ,o_n\left(t\right)\right]}^{\mathrm{T}}$ representative constant and the representative initial value. From Equation (2), fractional integration operator $J_{t_0}^q$ on both sides of Equation (16), and we get

$$x=J_{t_0}^{q}Lx+J_{t_0}^{q}Nx+J_{t_0}^{q}o+\frac{{\left(t-t_0\right)}^k}{k!}\sum _{k=0}^{m-1}{b}_k.$$

(17)

where $t\in \left(t_0,t_1\right)$. From the properties of fractional integral operators [38], the nonlinear term in Equation (17) can be decomposed according to Equation (18):

$$A_j^i=\frac{1}{i!}{\left[\frac{d^i}{d{\lambda }^i}N\left(\sum _{k=0}^i{\left(\lambda \right)}^k{x}_j^k\right)\right]}_{\lambda =0}.$$

(18)

where $i=0,1,\cdots \infty$, $j=0,1,\cdots ,n$. Then, the nonlinear term is

$$Nx=\sum _{i=0}^{\infty }{A}^i\left(x^0,x^1,x^2,\cdots ,x^i\right).$$

(19)

Therefore, the numerical solution of Equation (16) is as follows:

$$\left\{\begin{array}{c}{x}^0=J_{t_0}^{q}o+\frac{{\left(t-t_0\right)}^k}{k!}\sum _{k=0}^{m-1}{b}_k,\hfill \\ x^1=J_{t_0}^{q}L{x}^0+J_{t_0}^q{A}^0\left(x^0\right),\hfill \\ x^2=J_{t_0}^{q}L{x}^1+J_{t_0}^q{A}^1\left(x^0,x^1\right),\hfill \\ \cdots \hfill \\ x^i=J_{t_0}^{q}L{x}^{i-1}+J_{t_0}^q{A}^{i-1}\left(x^0,x^1,\cdots ,x^{i-1}\right),\hfill \\ \cdots \hfill \end{array}\right.$$

(20)

The derivation times of Equation (20) approach infinity; that is, when the derivation terms approach infinity, its value is closer to the analytic solution.

2.3.2. Decomposition Form of Solution

According to the decomposition algorithm in Section 2.3.1 above, the fractional-order chaotic system (15) is decomposed, and the linear terms, nonlinear terms, and constant terms are obtained as $L{\left(x_1,x_2,x_3,x_4\right)}^{\mathrm{T}}={\left(a{x}_1+b{x}_2,x_1-x_2-x_3,d{x}_2,e{x}_1-f{x}_4\right)}^{\mathrm{T}}$, $N{\left(x_1,x_2,x_3,x_4\right)}^{\mathrm{T}}={\left(c{x}_1{x}_4{x}_4,0,0,0\right)}^{\mathrm{T}}$, ${\left(0,0,0,0\right)}^{\mathrm{T}}=\mathbf{O}$. Thus, from Equation (17), the solution of this system is expressed by

$$\left(\begin{array}{c}{x}_1\left(t\right)\\ x_2\left(t\right)\\ x_3\left(t\right)\\ x_4\left(t\right)\end{array}\right)=J_{t_0}^q\left(\begin{array}{c}\begin{array}{c}a{x}_1+b{x}_2\\ x_1-x_2-x_3\\ d{x}_2\end{array}\\ e{x}_1-f{x}_4\end{array}\right)+J_{t_0}^q\left(\begin{array}{c}c{x}_1{x}_4{x}_4\\ 0\\ 0\\ 0\end{array}\right)+\mathbf{O}+\left(\begin{array}{c}{x}_1\left(t_0\right)\\ x_2\left(t_0\right)\\ x_3\left(t_0\right)\\ x_4\left(t_0\right)\end{array}\right).$$

(21)

The nonlinear term of Equation (21) is $x_1{x}_4{x}_4$, and we can obtain the nonlinear term decomposition from Equation (18), as shown in the Appendix A section. As a consequence, the form of the solution of the system after taking the finite-term iteration is as follows:

$$\begin{array}{c}{\tilde{x}}_i\left(t\right)=x_i^0+x_i^1+x_i^2+x_i^3+x_i^4+x_i^5\hfill \\ =c_i^0+c_i^1\frac{{\left(t-t_0\right)}^q}{\Gamma (q+1)}+c_i^2\frac{{\left(t-t_0\right)}^{2q}}{\Gamma (2q+1)}+c_i^3\frac{{\left(t-t_0\right)}^{3q}}{\Gamma (3q+1)}+c_i^4\frac{{\left(t-t_0\right)}^{4q}}{\Gamma (4q+1)}+c_i^5\frac{{\left(t-t_0\right)}^{5q}}{\Gamma (5q+1)}.\hfill \end{array}$$

(22)

where $i=1,2,3,4$. Based on this approximate analytical solution, the parameter f is selected as the system variable. The others parameter values are $a=3$, $b=12$, $c=2$, $d=32$, $e=-37$, $q=0.9$, and the initial condition $\left(0.1,0.1,0.1,0.1\right)$. By changing the parameter f, the system phase diagram of Figure 4 can be obtained.

Figure 4d depicts a single-scroll symmetric chaotic attractor with $f=9$, and when $f=23$, $f=19$, and $f=13$, we can obtain the periodic-1 attractor, periodic-2 attractor, and single-scroll chaotic attractor as shown in Figure 4a–c, respectively.

3. Dynamical Analysis

3.1. Dynamical Analysis with Fractional Orders q

To study the dynamical behaviors of system (15), we set $q\in [0.7,1]$, and the other system parameters were $a=3$, $b=12$, $c=2$, $d=32$, $e=-37$, $f=12$, and $(0.1,0.1,0.1,0.1)$. Figure 5 shows the bifurcation diagram, LEs diagram, and that order q is the variable.

From Figure 5a, there are three distinct periodic windows, i.e., $q\in \left[0.70,0.76\right]$, $q\in \left[0.77,0.78\right]$, and $q\in \left[0.91,1\right]$ can be found, and, except for the periodic, the others regions are in a chaotic state. The same result can be obtained from LEs in Figure 5b. We list the state of the system (15) under different orders q in

Table 2. Dynamical behaviors of system (15) as the order q changes.

q	Dynamical Properties	Phase Diagrams
(0.700,0.748)	periodic-1	Figure 6a
(0.748,0.754)	periodic-2	Figure 6b
(0.754,0.780)	periodic-4	Figure 6c
(0.780,0.850)	chaotic	Figure 6d
(0.850,0.910)	chaotic	Figure 6e
(0.910,1.000)	periodic-2	Figure 6f

, and the corresponding dynamical behaviors are shown in Figure 6.

Figure 6 shows the dynamical behaviors with different orders q of the system (15), which are consistent with the corresponding dynamical properties in Table 2. When fractional-order q is taken as a variable, the system can produce richer and more complex attractors.

3.2. Dynamical Analysis with Parameters

3.2.1. Parameter a as the Variable

We fixed the order $q=0.9$, the parameter a is used as a variable and selected the system parameters as $b=12,c=2,d=32,e=-37,f=12$, the initial value as $(0.1,0.1,0.1,0.1)$, similarly. Then we obtained the bifurcation diagram and LEs diagram as show in Figure 7.

In Figure 7, it can be found that the system changes back and forth between a chaotic and periodic state, and that the dynamical behaviors are very rich.

3.2.2. Parameter d as the Variable

Likewise, we fixed the order q, used the parameter d as a variable, and selected the system parameters as $a=3,b=12,c=2,e=-37,f=12$ and the initial value as $(0.1,0.1,0.1,0.1)$, similarly. Then, we obtained the bifurcation diagram and LEs diagram, as shown in Figure 8.

Regarding studying the dynamical behaviors of $d\in \left(28,32\right)$, Figure 8 shows the bifurcation diagram and LEs when selecting parameter d as the variable, and the numerical research of the system (15) indicates that, in addition to three periodic windows, i.e., $d\in \left(28.61,28.82\right)$, $d\in \left(29.23,29.43\right)$ and $d\in \left(30.68,31.01\right)$, the system is in a chaotic state in most of the other parameter ranges.

3.2.3. Parameter f as the Variable

Under the same parameters and conditions as the above analysis, $f\in \left(8,12\right)$ is selected as the variable, and the bifurcation diagram and LEs diagram are shown in Figure 9.

From the above Figure 9, it can be found thatm when parameter f is a variable, the dynamical behaviors of the system (15) exhibit alternating chaotic and periodic states.

3.3. Dynamical Analysis with Initial Value

Like before, we fixed the order $q=0.9$, the parameter was selected as $a=3,$$b=12,c=2,d=32,e=-37$, and $f=12$, and then we changed the initial value $x_1\in [0.2,0.2]$ and $x_2\in [0.2,0.2]$, whereas $x_3$ and $x_4$ were not changed. Figure 10 and Figure 11 show the bifurcation diagram and LEs diagram of system (15) with $x_1$ and $x_2$ as variables.

Figure 10 and Figure 11 show that the system (15) has symmetrical changes, i.e., when other conditions remain constant and the initial value changes, the system can emerge with coexisting attractors, as shown in Figure 12.

When some initial values of system (15) change, it can be found that the position of the attractor changes. This illustrates the multiple stability of the system, where it is just a dynamical property exhibited by the initial value space.

3.4. Offset Boosting

For a nonlinear system, adding a constant term after the uncoupled linear term of the variable that occurs only once in the system can offset the system (15). The constant is used as the control quantity, and the offset boosting can control the system offset. In Equation (15), variable $x_2$ satisfies the conditions of phase offset boosting. The constant term $\overline{p}$ is added to the second dimension so we can obtain

$$\left\{\begin{array}{c}{D}_t^q{x}_1=a{x}_1+b{x}_2-c{x}_1{x_4}^2,\hfill \\ D_t^q{x}_2=x_1-x_2-\left(x_3+\overline{p}\right),\hfill \\ D_t^q{x}_3=d{x}_2,\hfill \\ D_t^q{x}_4=e{x}_1-f{x}_4.\hfill \end{array}\right.$$

(23)

For the study of the offset boosting phenomenon, the parameters are consistent with the above, the initial values are $\left(0.1,0.1,0.1,0.1\right)$, and the offset $\overline{p}$ is set to 0, 0.2, and 0.4, respectively. When we change the fractional-orders q, we can find the offset boosting shown in Figure 13.

The system exhibits rich dynamical behaviors. To further study the offset boosting of system (15), under the same conditions, we only changed the value of f, in which, the offset $\overline{p}$ was set to 0, 0.2, 0.4, and 0.6, and Figure 14 shows the position offset of the attractor under different parameters f.

As shown in Figure 14, when $f=9$, the system is shown as being in a double-scroll chaotic state. When $f=11$, the system remains in a double-scroll periodic state. When $f=14$, the system is shown as being in a single-scroll chaotic state. When $f=17$, the system remains in a single-scroll periodic state. Similarly, if we change the other parameters, we can find more offset boosting, indicating the rich and complex dynamical behaviors of the system.

4. System Circuit Implementation

4.1. Laplace Transform of Fractional-Order Memristor Circuit

Since it is impossib;le to perform a fractional differential calculation of a function or system directly in the time-domain, we can approximate the fractional-order operator using an integer-order operator as follows:

$$L\left\{{}_0{D}_t^{q}x\left(t\right)\right\}=\underset{0}{\overset{\infty }{\int }}{{\mathrm{e}}^{-z}}_0{D}_t^{q}x\left(t\right)dt=s^{q}X\left(s\right)-\sum _{k=0}^{n-1}{s}^k\left({\left.{}_0{D}_t^{q-k-1}x\left(t\right)\right|}_{t=0}\right)$$

(24)

where the order $n-1\le q<n$, the initial value of the fractional-order system (15) is zero, and, in the time-domain transformation, the transfer function $H\left(s\right)=\frac{H\left(0\right)}{s^q}$ can be used to replace the fractional operator q. The Laplace transform can be used for system identification by using $q=0.9$, where the maximum approximation error is 2 dB, and by considering the $\frac{1}{s^q}$ Bode diagram of the system (15). Thus, we can obtain

$$\left\{\begin{array}{c}{s}^{0.9}{X}_1\left(s\right)=a{X}_1\left(s\right)+b{X}_2\left(s\right)-c{X}_1\left(s\right){X_4}^2\left(s\right)\hfill \\ s^{0.9}{X}_2\left(s\right)=X_1\left(s\right)-X_2\left(s\right)-X_3\left(s\right)\hfill \\ \begin{array}{c}{s}^{0.9}{X}_3\left(s\right)=d{X}_2\left(s\right)\hfill \\ s^{0.9}{X}_4\left(s\right)=e{X}_1\left(s\right)-f{X}_4\left(s\right)\hfill \end{array}\hfill \end{array}\right.$$

(25)

According to the $\frac{1}{s^q}$ Bode diagram (the maximum approximation error is 2 dB), we have

$$\frac{1}{s^{0.9}}\approx \frac{2.2675\left(s+1.292\right)\left(s+215.4\right)}{\left(s+0.01292\right)\left(s+2.154\right)\left(s+359.4\right)}$$

(26)

Here, let $g=2.2675$, $k=216.692$, $l=278.2968$, $m=361.567$, $n=778.819$, $p=10$, i.e., $\frac{g\left(s^2+ks+l\right)}{s^3+m{s}^2+ns+p}=\frac{2.2675\left(s+1.292\right)\left(s+215.4\right)}{\left(s+0.01292\right)\left(s+2.154\right)\left(s+359.4\right)}$, and, by Equation (25), we can obtain

$$\left\{\begin{array}{c}{X}_1\left(s\right)=\frac{g\left(s^2+ks+l\right)}{s^3+m{s}^2+ns+p}\left(a{X}_1\left(s\right)+b{X}_2\left(s\right)-c{X}_1\left(s\right){X_4}^2\left(s\right)\right)\hfill \\ X_2\left(s\right)=\frac{g\left(s^2+ks+l\right)}{s^3+m{s}^2+ns+p}(X_1\left(s\right)-X_2\left(s\right)-X_3\left(s\right))\hfill \\ \begin{array}{c}{X}_3\left(s\right)=\frac{g\left(s^2+ks+l\right)}{s^3+m{s}^2+ns+p}d{X}_2\left(s\right)\hfill \\ X_4\left(s\right)=\frac{g\left(s^2+ks+l\right)}{s^3+m{s}^2+ns+p}(e{X}_1\left(s\right)-f{X}_4\left(s\right))\hfill \end{array}\hfill \end{array}\right.$$

(27)

Let ${\chi }_1=x$, ${\chi }_2=\dot{x}$, ${\chi }_3=\ddot{x}$, ${\chi }_4=\stackrel{⃛}{x}$, ${\chi }_5=y$, ${\chi }_6=\dot{y}$, ${\chi }_7=\ddot{y}$, ${\chi }_8=\stackrel{⃛}{y}$, ${\chi }_9=z$, ${\chi }_{10}=\dot{z}$, ${\chi }_{11}=\ddot{z}$, ${\chi }_{12}=\stackrel{⃛}{z}$, ${\chi }_{13}=w$, ${\chi }_{14}=\dot{w}$, ${\chi }_{15}=\ddot{w}$, ${\chi }_{16}=\stackrel{⃛}{w}$; then, Equation (27) can be further integrated into a 16-order integer-order differential equation, as shown in

Table 3. Sixteen-order integer-order differential equation of system (27).

$\frac{\mathrm{d}{\chi }_1}{\mathrm{d}t}={\chi }_2$	$\frac{\mathrm{d}{\chi }_2}{\mathrm{d}t}={\chi }_3$	$\frac{\mathrm{d}{\chi }_3}{\mathrm{d}t}={\chi }_4$	$\begin{array}{c}\frac{\mathrm{d}{\chi }_4}{\mathrm{d}t}=\left(ag-m\right){\chi }_3+\left(agk-n\right){\chi }_2\hfill \\ +\left(agl-p\right){\chi }_1+bg{\chi }_7+bgk{\chi }_6+bgl{\chi }_5\hfill \\ -cg{\chi }_3{\chi }_{15}^2-cgkcg{\chi }_2{\chi }_{14}^2-cg{\chi }_1{\chi }_{13}\hfill \end{array}$
$\frac{\mathrm{d}{\chi }_5}{\mathrm{d}t}={\chi }_6$	$\frac{\mathrm{d}{\chi }_6}{\mathrm{d}t}={\chi }_7$	$\frac{\mathrm{d}{\chi }_7}{\mathrm{d}t}={\chi }_8$	$\begin{array}{c}\frac{\mathrm{d}{\chi }_8}{\mathrm{d}t}=\mathrm{g}{\chi }_3-\left(g-m\right){\chi }_7-g{\chi }_{11}+gk{\chi }_2\hfill \\ -\left(gk-n\right){\chi }_6-gk{\chi }_{10}+g{\chi }_1-\left(g-p\right){\chi }_5-g{\chi }_9\hfill \end{array}$
$\frac{\mathrm{d}{\chi }_9}{\mathrm{d}t}={\chi }_{10}$	$\frac{\mathrm{d}{\chi }_{10}}{\mathrm{d}t}={\chi }_{11}$	$\frac{\mathrm{d}{\chi }_{11}}{\mathrm{d}t}={\chi }_{12}$	$\frac{\mathrm{d}{\chi }_{12}}{\mathrm{d}t}=dg{\chi }_7+dk{\chi }_6+dl{\chi }_5-m{\chi }_{11}-n{\chi }_{10}-p{\chi }_9$
$\frac{\mathrm{d}{\chi }_{13}}{\mathrm{d}t}={\chi }_{14}$	$\frac{\mathrm{d}{\chi }_{14}}{\mathrm{d}t}={\chi }_{15}$	$\frac{\mathrm{d}{\chi }_{15}}{\mathrm{d}t}={\chi }_{16}$	$\begin{array}{c}\frac{\mathrm{d}{\chi }_{16}}{\mathrm{d}t}=ge{\chi }_3+gke{\chi }_2+gle{\chi }_1\hfill \\ -\left(gf-m\right){\chi }_{15}-\left(gkf-n\right){\chi }_{14}-\left(glf-p\right){\chi }_3\hfill \end{array}$

.

The initial values of $x\left(t\right)$,$\frac{\mathrm{d}x\left(t\right)}{\mathrm{d}t}$; $y\left(t\right)$, $\frac{\mathrm{d}y\left(t\right)}{\mathrm{d}t}$; $z\left(t\right)$, $\frac{\mathrm{d}z\left(t\right)}{\mathrm{d}t}$; $w\left(t\right)$, $\frac{\mathrm{d}w\left(t\right)}{\mathrm{d}t}$ are considered as zero throughout time-domain transformations, i.e., ${\chi }_1\left(0\right)=0$, ${\chi }_2\left(0\right)=0$, ${\chi }_3\left(0\right)=0$, ${\chi }_4\left(0\right)=0.1$, ${\chi }_5\left(0\right)=0$, ${\chi }_6\left(0\right)=0$, ${\chi }_7\left(0\right)=0$, ${\chi }_8\left(0\right)=0.1$, ${\chi }_9\left(0\right)=0$, ${\chi }_{10}\left(0\right)=0$, ${\chi }_{11}\left(0\right)=0$, ${\chi }_{12}\left(0\right)=0.1$, ${\chi }_{13}\left(0\right)=0$, ${\chi }_{14}\left(0\right)=0$, ${\chi }_{15}\left(0\right)=0$, ${\chi }_{16}\left(0\right)=0.1$ in Table 3.

Therefore, a fractional-order system (15) can be converted into an integer-order system by time-frequency domain conversion and setting the step size to 0.1. It is a useful method for calculating fractional-order circuit systems using the Laplace transform. The fractional-order memristor analog circuit will then be implemented using the aforementioned method to validate the previous numerical simulation.

4.2. Circuit Simulation and FPGA Implementation

4.2.1. Analog Circuit Simulation

The analog circuit was designed as shown in Figure 15 in order to verify that the system (15) is functioning correctly. Multiplication was implemented by the AD633 multiplier, and the integral and addition operational amplifiers were implemented by TL082. All of the chips needed in the circuit were powered by $\pm 15$ V DC.

The principle equation of fractional-order memristor circuit is given according to the equivalent approximate fractional-order operator. According to the system (15), the equivalent circuit equation can be obtained as follows:

$$\left\{\begin{array}{c}\frac{d^q{x}_1}{dt}=\frac{RS}{R_1}\frac{1}{RB{C}_0}\frac{RA}{R{A}_1}{x}_1+\frac{RS}{R_2}\frac{1}{RB{C}_0}\frac{RA}{R{A}_1}{x}_2+\frac{RS}{R_3}\frac{1}{RB{C}_0}\frac{RA}{R{A}_1}{x_4}^2(-x_1),\hfill \\ \frac{d^q{x}_2}{dt}=\frac{R{S}_1}{R_4}\frac{1}{R{B}_1{C}_0}\frac{R{A}_2}{R{A}_3}{x}_1+\frac{R{S}_1}{R_5}\frac{1}{R{B}_1{C}_0}\frac{R{A}_2}{R{A}_3}(-x_2)+\frac{R{S}_1}{R_6}\frac{1}{R{B}_1{C}_0}\frac{R{A}_2}{R{A}_3}(-x_3),\hfill \\ \frac{d^q{x}_3}{dt}=\frac{R{S}_2}{R_7}\frac{1}{R{B}_2{C}_0}\frac{R{A}_4}{R{A}_5}{x}_2,\hfill \\ \frac{d^q{x}_4}{dt}=\frac{R{S}_3}{R_8}\frac{1}{R{B}_3{C}_0}\frac{R{A}_6}{R{A}_7}(-x_1)+\frac{R{S}_3}{R_9}\frac{1}{R{B}_3{C}_0}\frac{R{A}_4}{R{A}_5}(-x_4).\hfill \end{array}\right.$$

(28)

where $R_i$, $RA$, $R{A}_i$, $RB$, $R{B}_i$, $RS$, $R{S}_i$, and so on are variable resistors corresponding to variable parameters in system (15). In Equation (28), $C_0$ represents the fractional-order modules in parallel on operational amplifiers U2, U5, U8, and U11. The circuit simulation results under different resistances $R_9$ are shown in Figure 16, and the resistance values, capacitance values, and corresponding parameters are shown in Figure 15. The simulation results of the circuit verify the symmetric and multistability of the system (15) well.

The numerical simulation results under different parameters f are consistent with the simulation results of the analog circuit under different resistances $R_9$, which proves that the digital circuit can verify the chaotic behaviors of this system, including symmetric multi-stability.

4.2.2. FPGA System Implementation

The selected system’s parameters are consistent with those of the aforementioned dynamics study: the step is $h=0.01$, and the initial value is $(x_1\left(0\right),x_2\left(0\right),x_3\left(0\right),$$x_4\left(0\right))=(0.1,0.1,0.1,0.1)$. The FPGA digital circuit of the system (15) was realized using the FPGA chip EPLC6Q240 of the Altera Cyclone4 series. To capture the system’s phase diagram on the digital oscilloscope, the only change was to the system parameter f, as shown in Figure 17.

In practice, the FPGA-based system design has received more and more attention. During the design of a fractional-order memristor system, modules are divided according to their functions at the top-level, and the relationship and cooperation mode between modules are planned. Verilog hardware is used at the bottom to realize the functions of each module. The clock frequency is 51.13 MHz, and the proportion of registers used is $11\%$. The attractor phase diagram displayed by the oscilloscope in Figure 17 further verifies the existence of the system and the possibility of engineering applications.

4.3. Fractional-Order Memristor Control Circuit

Consider the master system of fractional-order memristor circuit (28). The slave system can be written as

$$\left\{\begin{array}{cc}\hfill \frac{d^q{y}_1}{dt}& =\frac{RS}{R_1}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_1+\frac{RS}{R_2}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_2+\frac{RS}{R_3}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_4^2\left(-y_1\right)+U_1,\hfill \\ \hfill \frac{d^q{y}_2}{dt}& =\frac{R{S}_1}{R_4}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}{y}_1+\frac{R{S}_1}{R_5}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(-y_2\right)+\frac{R{S}_1}{R_6}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(-y_3\right)+U_2,\hfill \\ \hfill \frac{d^q{y}_3}{dt}& =\frac{R{S}_2}{R_7}\frac{1}{R{B}_2·C_0}\frac{R{A}_4}{R{A}_5}{y}_2+U_3,\hfill \\ \hfill \frac{d^q{y}_4}{dt}& =\frac{R{S}_3}{R_8}\frac{1}{R{B}_3·C_0}\frac{R{A}_6}{R{A}_7}\left(-y_1\right)+\frac{R{S}_3}{R_9}\frac{1}{R{B}_3·C_0}\frac{R{A}_4}{R{A}_5}\left(-y_4\right)+U_4.\hfill \end{array}\right.$$

(29)

Then, the error system can be written as $e_i=y_i-x_i(i=1,2,3,4)$. Thus, we can obtain

$$\left\{\begin{array}{c}\frac{d^q{e}_1}{dt}=\frac{RS}{R_1}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}\left(y_1-x_1\right)+\frac{RS}{R_2}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}\left(y_2-x_1\right)-\frac{RS}{R_3}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}\left(y_1{y_4}^2-x_1{x_4}^2\right)+U_1,\hfill \\ \frac{d^q{e}_2}{dt}=\frac{R{S}_1}{R_4}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(y_1-x_1\right)-\frac{R{S}_1}{R_5}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(y_2-x_2\right)-\frac{R{S}_1}{R_6}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(y_3-x_3\right)+U_2,\hfill \\ \frac{d^q{e}_3}{dt}=\frac{R{S}_2}{R_7}\frac{1}{R{B}_2·C_0}\frac{R{A}_4}{R{A}_5}\left(y_2-x_2\right)+U_3,\hfill \\ \frac{d^q{e}_4}{dt}=-\frac{R{S}_3}{R_8}\frac{1}{R{B}_3·C_0}\frac{R{A}_6}{R{A}_7}\left(y_1-x_1\right)-\frac{R{S}_3}{R_9}\frac{1}{R{B}_3·C_0}\frac{R{A}_4}{R{A}_5}\left(y_4-x_4\right)+U_4.\hfill \end{array}\right.$$

(30)

In this case, we give $U_1=\frac{RS·RA·y_1}{C_0·RB·R{A}_1·R_r}$ while setting the other controls $U_2,U_3,U_4$ to 0. Consequently, the system equation of the circuit controlled by a fractional-order memristor is as follows:

$$\left\{\begin{array}{c}\frac{d^q{y}_1}{dt}=\frac{RS}{R_1}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_1+\frac{RS}{R_2}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_2+\frac{RS}{R_3}\frac{1}{RB·C_0}\frac{RA}{R{A}_1}{y}_4^2\left(-y_1\right)+\frac{RS·RA}{C_0·RB·R{A}_1·R_r}{y}_1,\hfill \\ \frac{d^q{y}_2}{dt}=\frac{R{S}_1}{R_4}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}{y}_1+\frac{R{S}_1}{R_5}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(-y_2\right)+\frac{R{S}_1}{R_6}\frac{1}{R{B}_1·C_0}\frac{R{A}_2}{R{A}_3}\left(-y_3\right),\hfill \\ \frac{d^q{y}_3}{dt}=\frac{R{S}_2}{R_7}\frac{1}{R{B}_2·C_0}\frac{R{A}_4}{R{A}_5}{y}_2,\hfill \\ \frac{d^q{y}_4}{dt}=\frac{R{S}_3}{R_8}\frac{1}{R{B}_3·C_0}\frac{R{A}_6}{R{A}_7}\left(-y_1\right)+\frac{R{S}_3}{R_9}\frac{1}{R{B}_3·C_0}\frac{R{A}_4}{R{A}_5}\left(-y_4\right).\hfill \end{array}\right.$$

(31)

The values of each device in Figure 15 can be determined using the system parameters $a=3$, $b=12$, $c=2$, $d=32$, $e=-37$, and $q=0.9$. Naturally, $O=(0,0,0,0)$ is the system’s equilibrium, and the Jacobian matrix following the circuit system’s linearization is

$$J=\left(\begin{array}{cccc}\frac{RS}{R_1}\frac{1}{RB{C}_0}\frac{RA}{R{A}_1}-\frac{RS·RA}{C_0·RB·R{A}_1·R_r}& \frac{RS}{R_2}\frac{1}{RB{C}_0}\frac{RA}{R{A}_1}& 0& 0\\ 1& -1& -1& 0\\ 0& \frac{R{S}_2}{R_7}\frac{1}{R{B}_2{C}_0}\frac{R{A}_4}{R{A}_5}& 0& 0\\ \frac{R{S}_3}{R_8}\frac{1}{R{B}_3{C}_0}\frac{R{A}_6}{R{A}_7}& 0& 0& \frac{R{S}_3}{R_9}\frac{1}{R{B}_3{C}_0}\frac{R{A}_4}{R{A}_5}\end{array}\right).$$

(32)

where $U_1=\frac{RS·RA}{C_0·RB·R{A}_1·R_r}=3.3$. Then, we can determine that the eigenvalues are ${\lambda }_1=13.7868$, ${\lambda }_2=-1.2396+3.4949i$, ${\lambda }_3=-1.2396-3.4949i$, and ${\lambda }_4=-0.6076$, respectively.

Lemma 1.

The stability lemma of fractional-order systems [40].

The commensurate system is stable if the following condition is satisfied:

$$\left|arg\left({\lambda }_i\right)\right|=\left|arg\left(\mathrm{Re}\left({\lambda }_i\right)+\mathrm{Im}\left({\lambda }_i\right)\right)\right|>\frac{q\pi }{2}.$$

(33)

where ${\lambda }_i\left(i=1,2,\cdots ,n\right)$ are eigenvalues, and the order q must satisfy

$$\begin{array}{c}q<\frac{2}{\pi }arctan\left(\frac{\left|\mathrm{Im}\left({\lambda }_i\right)\right|}{\left|\mathrm{Re}\left({\lambda }_i\right)\right|}\right).\end{array}$$

(34)

Accordingly, from Lemma 1, $\left|arg\left({\lambda }_i\right)\right|=2.8$, i.e., $\left|arg\left({\lambda }_i\right)\right|>\frac{q\pi }{2}=1.41$, which satisfies Equation (33), occurs when the system is at zero equilibrium O. As a result, the system with $U_i(i=1,2,3,4)$ can gradually stabilize to the equilibrium, and the fractional-order memristor feedback control circuit is designed as shown in Figure 18.

From the first channel of the feedback control circuit of Figure 18, when the resistance $R_r=30$ k$\mathsf{\Omega }$ and the switch $S1$ on the control circuit is closed at time $t=30$, the time-domain waveform of the state variable of the circuit experiment simulation of the memristor chaotic system feedback control circuit system is designed as shown in Figure 19.

According to the circuit experiment results of the time-domain waveform diagram of state variables (time responses of error system) in Figure 19, it can be found that its trajectory asymptotically converges to the origin by adding a feedback controller, which shows the control effect of the designed controller.

5. Conclusions

In this work, a novel fractional-order memristor system was established by integrating the fractional-order memristor into Chua’s circuit. Following that, the bifurcation diagram and LEs with system order, parameters, and initial values as independent variables were studied. The study revealed that altering these variables can cause the system to produce dynamical characteristics that are richer and more complex. In particular, the system exhibits periodic limit cycles and a single-scroll chaotic attractor; when the parameters and order of the system are changed, the system can generate a symmetric single-scroll chaotic attractor when the initial value is changed. Then, by adding a constant term after the decoupled linear term of the system, it was found that there is an offset boosting phenomenon in the system. Finally, the designed FPGA and analog circuit were established in accordance with the circuit equation to verify the numerical calculation results. The control circuit of the system was designed through the stability analysis of the equilibrium to satisfy the control requirements.

Based on the results obtained, the dynamics of a novel fractional-order memristor chaotic system studied in this paper are very rich due to the fractional-order operator, memristor, and offset boosting control, which provide a theoretical foundation for spread spectrum communication, pseudo-random sequence generators, and other practical engineering applications. Meanwhile, the possibility of the next work was further verified through the implementation of FPGA and analog circuits. Therefore, in the following work, we will summarize the shortcomings of this work, design and implement a fractional-order chaotic signal generator (floating point multiplier, floating point adder, fractional-order chaotic operation controller), and explore a novel fractional-order memristor chaotic system control scheme.