Abstract

This paper investigated the global synchronization of fractional-order memristive neural networks (FMNNs). To deal with the effect of reaction-diffusion and time delay, fractional partial and comparison theorem are introduced. Based on the set value mapping theory and Filippov solution, the activation function is extended to discontinuous case. Adaptive controllers with a compensator are designed owing to the existence of unknown parameters, with the help of Gronwall–Bellman inequality. Numerical simulation examples demonstrate the availability of the theoretical results.

1. Introduction

As a generalization of differential and integral calculus from an integer order to arbitrary order, fractional calculus has more than 300 years history [1]. Since Mandelbort proposed the fact that lots of fractals exist in the field of technology, a whole new array of possibilities for natural science research has been opened [2]. Compared with the conventional integer order calculus, it shows some superiority in the aspect of memory and heredity [3]. That caused its wide range of applications [46].

Neural networks model, as we know, has great potential for controlling highly nonlinear and severe uncertain systems. Due to the ability of studying arbitrary continuous nonlinear functions, it has been fruitfully applied to optimization, secure communication, cryptography, and so on. Taking into account these facts, the combination of fractional calculus and neural networks model is a remarkably great improvement. This kind of model called fractional-order neural networks (FNNs) was proposed by Arena in 1998 [7]. Afterwards, a myriad of studies have focused on the dynamical behaviors of FNNs. One of the most interesting issues is to explore the synchronization problem of FNNs. Various definition types of synchronization have been built [810]. Among them, global synchronization is the most basic and widely applied one, since it guarantees the convergence of the system. From the perspective of engineering applications, it has achieved the required realistic standards [11, 12]. Thus, global synchronization of FNNs has been extensively investigated to date [1315].

Another hot topic of research is memristor-based neural networks, which were first put forward by Chua [16]. In 2008, a nanoscale memory was made by Williams group [17], which was named memristor. Memory characteristics and nanometer dimensions make the memristor show properties the same as the neurons in the human brain. From this standpoint, fractional-order memristive neural networks (FMNNs) have been suggested replacing the FNNs to emulate the human brain. Therefore, the analysis of FMNNs has become increasingly attractive to researchers. At present, there are many excellent results on synchronization of FMNNs [1821]. Velmurugan investigated finite-time synchronization of FMNNs with the method of Laplace transform and generalized Gronwalls inequality [18]. The lag complete synchronization criteria of FMNNs are proposed by period intermittent control [19]. Delay-independent synchronization criteria for FMNNs are established by using the maximum modulus principle, the spectral radii of matrix, and Kakutani fixed point theorem [20]. Chen et al. designed a stabilizing state-feedback controller by employing the FO Razumikhin theorem and linear matrix inequalities [21]. Besides the abovementioned methods, adaptive mechanism has been proposed as a more effective strategy [22]. Compared with the conventional feedback control, the parameters of adaptive controllers will be updated automatically, which makes the energy be saved [23]. It also has significant effect on the analysis of fractional order nonlinear system, especially when the system has unknown parameters [2427]. Thus, it is a natural idea to investigate the synchronization of FMNNs using adaptive control scheme.

It is a remarkable fact that time delays emerge in the communication of neurons [28, 29]. Furthermore, the occurrence of diffusion phenomena has a significant impact in the electrons’ transportation through a nonuniform electromagnetic field. That means, not only the time but also their spatial positions affect the dynamic behaviors of networks. The models covering time delay and reaction-diffusion are good mimicry for the real neural networks in terms of application, but the existence of time delays and reaction-diffusion could bring about some undesirable behaviors [30, 31]. Drawing on the abovementioned concerns, special attention should be paid on time delay and reaction-diffusion. However, the research on FMNNs with both reaction-diffusion and time delay is still in its infancy [32, 33].

For realistic neural networks systems, one must include uncertainties. Uncertainty could occur in parameters in the mathematical model, for they are hardly to be exactly known in advance [3436]. Therefore, the study on synchronization of the drive-response system considering unknown parameters is one of the most challenging tasks. Expanding on the existing literature, the global synchronization problem for reaction-diffusion FMNNs with time delay and unknown parameters is proposed firstly. The main contents of this paper are the following. (1) Based on fractional comparison theorem, the problem of synchronization for the integral delayed neural networks is generalized to the fractional order case. (2) By utilizing Caputo partial fractional derivative, reaction-diffusion is emphasized for the FMNNs, which extends the problem from the time domain to time-spatial dimension. (3) With the help of functional differential inclusions, the solution of the discontinuous systems is guaranteed in the sense of Filippov. (4) Employing Gronwall–Bellman inequality, the adaptive controllers with compensator are designed to achieve the synchronization, despite the presence of parametric uncertainties.

This paper is structured as follows. In Section 2, we briefly introduce preliminaries and model description. Two types of adaptive controllers for global synchronization of the drive-response system are designed in Section 3. In Section 4, two simulations are shown to verify the correctness of the proposed results. A short conclusion and some discussions are given in Section 5.

Notations. let be the space of -dimensional Euclidean space. The norm of vector is defined as . A bounded open-set containing the origin equips with the smooth boundary and .

2. Preliminaries and Model Description

Preliminaries about fractional calculus and the model of reaction-diffusion FMNNs with time delay are introduced in this section. Concept of Filippov solution, fractional comparison theorem, and Gronwall–Bellman inequality need to be emphasized.

2.1. Preliminaries

Definition and the properties about Caputo Riemann–Liuville fractional calculus are given in the following.

Definition 1. (see [37]). For any , the time Caputo fractional derivative of order for a function is defined bywhere . Particularly, when , thenIt should be noted that the initial conditions of Caputo fractional derivative are same as integer-order derivative, which makes Caputo fractional derivative equip specific physical sense in the real world. In this paper, simple notation is used to replace .

Definition 2. (see [37]). Riemann–Liuville fractional integral of order for a function is defined bywhere .

Proposition 1. (see [37]). For arbitrary constants and , there is

Proposition 2. (see [37]). If and , then .
Several important conclusions about the fractional-order system are listed below.

Lemma 1. (see [33]). Let , be a continuous and differentiable function on ; then, for any , ,where . Particularly, when , there is

Lemma 2. (see [38]). Consider the system as follows:where and . The solution of the system is asymptotically stable, when .

Lemma 3 (fractional comparison theorem; see [38]). Consider the following fractional order inequality:and fractional order system:where , , , and are continuous and nonnegative functions in , on . If and , then the inequality holds:for all .

Lemma 4. (see [37, 39]). Suppose that is integrable on and is derivable respect to . Denote , then

Lemma 5. (see [40]). If with , then there is .

Lemma 6. (see [41]). Let be a cube , and let be a real-valued function belonging to which vanishes on the boundary of the , i.e., , then

Lemma 7. (Gronwall–Bellman inequality; see [42]). If satisfies , where are the known real functions, then

If is a constant,

2.2. System Description

Reaction-diffusion FMNNs with time delay are described aswhere , is the quantity of units; refers to the order of the equation; represents the state of the th unit at time and in space ; represents the time variable; is a Laplace operator, which means ; the smooth constants work as the transmission diffusion operator along the th unit at th position; , , where and stand for the strength of the th unit on the th unit at time in space and at time in space , respectively; is the self-feedback connection weight; and denote the discontinuous activation functions; signifies the external input; and stands for time delay. In addition, the memristor-based connection weights are given bywhere , , switching jumps and are all constants.

Consider the following Dirichlet-type boundary conditions of (15)and initial conditionswhere are real-value continuous functions defined on .

Take (15) as the drive system, then the corresponding response system is defined bywhere and is the controller which will be designed later.

Similarly, there areand , which have the same property to .

Since the right-hand sides of equations are discontinuous, the solutions of the equations in the conventional sense have been shown to be invalid. In this case, Filippov solutions [43] are introduced to deal with the discontinuous equations in functional differential inclusions framework. The solution’s global existence for fractional differential inclusions with time delay is given by [44].

Assumption 1. For each , is piecewise continuous function. Moreover, has at most limited the number of discontinuous points on any compact interval of . The abovementioned properties are the same to . Based on the definition of Filippov solution, (15) can be rewritten aswhere , are absolutely continuous on any compact subinterval of :Or equivalently, there are measurable functions , , , and satisfyfor , , .
For writing convenience, we denote as and as .
The IVP of the driven system (15) can be written asThe IVP of the response system (19) under the controller could be given bywhere . , , , and , .
Now, let be the synchronization error between (15) and (19). The IVP about the error system is of the formwhere , the symbol is written as and which stands for .
Denote as the norm of vector .

Definition 3. The drive system (15) and response system (19) is said to be globally synchronized, if .

Remark 1. Owing to the zero Dirichlet boundary conditions of (24) and (25), a zero equilibrium state of error system (26) exists. The global synchronization of systems (24) and (25) corresponds to the global asymptotic stability of the zero equilibrium of system (26).

3. Main Results

The abovementioned synchronization problem in Section 2 are considered in two cases, namely, the drive system (15) has known or unknown parameters.

3.1. Global Synchronization for Systems with Known Parameters

First of all, we assume that the parameters in (15) are known. Then, the following adaptive controllers are proposed to implement global synchronization.

Assumption 2. For , there are positive constants and such thatwhere for any ,For the convenience of description, the following notations are introduced:

Theorem 1. Under Assumption 1 and Assumption 2, the global synchronization of (24) and (25) can be achieved based on the following adaptive controllers:where , , and and are constants which satisfy the following inequalities:

Proof. Employ the following auxiliary function:By the means of Lemma 1, the inequality can be established:According to Lemma 4, taking the derivative of (32) yieldsUsing Green’s formula and Lemma 6, it can be obtained thatRecalling Assumption 1 and Assumption 2, the inequalities hold:Combining controller (30) designed in Theorem 1, it follows thatFrom Definition 2 and Property 2, one obtainsSubstituting (38) into (37) and noting that and meet inequality (31), it can be deduced thatwhere and .
According to Lemma 2 and Lemma 3, it can be assert that the global synchronization between (24) and (25) can be realized.
The proof of Theorem 1 completes.

Remark 2. When , the error system (26) be changed into the integer form. The conclusion is the same as the result in Theorem 3.1 in article [45]. If the reaction-diffusion in space disappears, then the states of the system’s variable are only relevant to time but independent of space. Thus, the error system (26) can be converted toFor convenience, the inequalities are given by

Corollary 1. Under Assumption 1 and Assumption 2, the global synchronization between (24) and (25) can be achieved based on the following adaptive controllers:where , , and and are constants which satisfy inequality (41).

Remark 3. In fact, synchronization of FMNNs with time delay has been investigated in [20]. The results about synchronization of the systems with single delay in Theorem 3 in [20] are the same as Corollary 1 obtained in this paper. Suppose that there is no time delay in the course of communication between the different neurons; then, error system (26) can be reduced toFor convenience, the inequalities are given by

Corollary 2. Under Assumption 1 and Assumption 2, the global synchronization between (24) and (25) can be realized based on the following adaptive controllers:where , , and and satisfy inequality (44).

Remark 4. Comparing the coefficients of the adaptive controllers (31) and (44), it can be concluded that the controllers for global synchronization of the system with time delay are more consumptive.

3.2. Global Synchronization for Systems with Unknown Parameters

In real neural networks, the parameters of the system may be unknown, so it is difficult to achieve synchronization by the conventional adaptive control schemes. In the following text, adaptive controllers with the dynamic compensators are introduced.

Here, suppose that self connection weight and input in the drive system (15) are unknown. Then, the IVP of the drive system (15) can be written as

The IVP of the response system (19) can be shown aswhere and and are treated as the compensation of and .

Then, for , the IVP of the synchronization error can be described as

Assumption 3. There are positive constants , , and such thatwhere , , and the notation represents .

Theorem 2. Under Assumption 1 and Assumption 3, the global asymptotic synchronization between (46) and (47) can realized based on the following adaptive controllers:where , , , , and are positive constants, and and are constants which satisfy the following inequalities:with and .

Proof. Employ the following auxiliary function:whereSimilarly to the proof of Theorem 1, one obtainsSince and fulfill condition (51), it followswhere .
By Lemma 5, it can be concluded thatAccording to the definition of , thusBy Lemma 7, it reveals thatTherefore,From (59), it is easy to see that the global synchronization between (46) and (47) can be realized.
The proof of Theorem 2 completes.

Remark 5. The result obtained in Theorem 2 is more practical for considering the situation where the parameters of the drive system (24) are unknown. Through the proof of Theorem 2, it can be asserted that the dynamical compensator is necessary for the unknown system. If the reaction-diffusion in space disappears, the error system (48) can be converted to

Corollary 3. Under Assumption 1 and Assumption 3, the global synchronization between (46) and (47) can be achieved based on the following adaptive controllers:where , , , , and .

Suppose that there is no time delay in the course of communication; then, the error system (48) can be changed into

Corollary 4. Under Assumption 1 and Assumption 3, (46) and (47) can achieve global synchronization based on the following adaptive controllers:where , , , , and .

4. Numerical Simulations

Several numerical simulations are presented to illustrate the main results. The first example is used to validate Theorem 1. The second one is used to text Theorem 2. Now, let us consider the 2D () reaction-diffusion FMNNs with time delay and zero Direclet boundary values. Here, , , and .

Example 1. Consider error system (26) with such parameters:where the memristive coefficients are defined bySuppose that time delay and discontinuous activation functions are . Let initial values of (24) and (25) bewhere and . Given the fractional order and . Under these parameters, the dynamical behaviors of the synchronization errors between (24) and (25) are shown in Figure 1. From the figure, it can be observed that the synchronization errors gradually increase as time goes on. Notice that and have only one discontinuous point 0, Assumption 1 is valid. According to Assumption 2, we can choose and . By calculation, the parameters of controllers are selected as , , , , , and . It is easy to verify that condition (9) in Theorem 1 is satisfied. Figure 2 suggests that the synchronization errors asymptotically converge to zero.
Theorem 1 is not only dependent on order but also related with reaction-diffusion coefficients. Figure 3 shows the trajectories of the error system at spatial position with different orders when the reaction-diffusion coefficients . While if the order is fixed , the trajectories of the error system at spatial position with different reaction-diffusion coefficients are shown in Figure 4. From the numerical results, we could note that the error system would achieve the asymptotic stability faster with the higher order. It also holds true with increasing reaction-diffusion coefficient.

Example 2. Consider the error system (48) with the following parameters:whereAssume , , and (47) with unknown parameters and , and initial values of compensators in (47) are given byActivation functions are and and time delay . Initial values of (46) and (47) are chosen aswhere and .
As seen in Figure 5, the error state is instable. Choosing , , and , we obtain , , and and , , , and , . Condition (51) in Theorem 2 is content clear. Figure 6 shows that the synchronization errors achieve stability rapidly under the adaptive controller (50).
For the error system (48), under different orders with the same reaction-diffusion coefficients and different reaction-diffusion coefficients with the same order simulations are carried out. Corresponding error state trajectories are depicted in Figures 7 and 8. It reveals that the asymptotic stability of the system varies with orders and reaction-diffusion coefficients. Compared with the former case, the figures show little differences with different parameters.

5. Conclusion and Discussion

The global synchronization problem for reaction-diffusion FMNNs with time delay was emphasized firstly in this paper. Both the known and unknown parameters cases were investigated, separately. Based on fractional comparison theorem and Gronwall–Bellman inequality, two kinds of adaptive control schemes were proposed. Numerical simulations were implemented to examine the effectiveness of the obtained results.

Conventional adaptive controllers have mostly been used for controlling neural networks system, while they cannot guarantee the synchronization of unknown systems. An adaptive controller with a dynamic compensator such as the one described in this paper can eliminate this problem. However, the true values of the unknown parameters could not be obtained. Further research is still in progress to design an adaptive controller with dynamic compensator to achieve the identity of parameters.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (no. 61772063) and Beijing Natural Science Foundation (no. Z180005).