Finite Time Synchronization for Stochastic Fractional-Order Memristive BAM Neural Networks with Multiple Delays

1. Introduction

Artificial intelligence has been an active field of research, and neural networks have emerged as a prominent branch due to their intelligence characteristics and potential for real-world applications. Neural networks have revolutionized the field of artificial intelligence by enabling computers to process and analyze large volumes of complex data with remarkable accuracy. These models are based on the structure and processes of the brain, where neurons are interconnected and communicated by using electrical signals. Neural networks are characterized by their remarkable ability to learn from data and improve their performance over time, without being explicitly programmed. Neural networks have evolved over time, with various models developed to address different types of problems. For example Cohen-Grossberg neural networks [1], Hopfield neural networks [2] and cellular neural networks [3].

Kosko’s bidirectional associative memory neural networks (BAMNNs) are noteworthy extension of the traditional single-layer neural networks [4]. The BAMNNs consist of two layers of neurons that are not interconnected within their own layer. In contrast, the neurons in different layers are fully connected, allowing for bidirectional information flow between the two layers. This unique structure enables the BAMNNs to function as both input and output layers, providing powerful information storage and associative memory capabilities. In signal processing, the BAMNNs can be used to filter signals or extract features, while in pattern recognition, it can classify images or recognize speech. In optimization, the BAMNNs can be used to identify the optimal solution, while in automatic control, it can be used to regulate or stabilize a system. The progress of artificial intelligence and the evolution of neural networks have created novel opportunities to tackle intricate issues in diverse domains. In summary, these advancements have paved the way for innovative problem-solving approaches that were previously unattainable. The BAMNNs’ unique architecture and capabilities make its a powerful tool for engineering applications, and its is anticipated that this technology will maintain its importance in future research and development, and continue to make substantial contributions to various fields.

Due to the restriction of network size and synaptic elements, the functions of artificial neural networks are greatly limited. If common connection weights and the self-feedback connection weights of BAMNNs are established by memristor [5], then its model can be built in a circuit. Memristor [6,7] is a circuit element in electronic circuit theory that behaves nonlinearly and has two terminals. Its unique feature has led to its widespread use and potential application in a variety of fields, including artificial intelligence, data storage, and neuromorphic computing. Adding memristor to neural networks makes it possible for artificial neural networks to simulate human brain on the circuit, which makes the research on memristor neural networks more meaningful. Therefore, the resistance in the traditional neural networks is replaced by memristor, and the BAMNNs based on memristor is formed (MBAMNNs). Compared with traditional neural networks, Memristor neural networks have stronger learning and associative memory abilities, allowing for more efficient processing and storage of information, thereby improving the efficiency and accuracy of artificial intelligence. Additionally, due to the nonlinear characteristics of memristors and their applications in circuits, Memristor neural networks also have lower energy consumption and higher speed. Therefore, the development of Memristor neural networks have broad application prospects in the field of artificial intelligence.

However, due to the limitation of amplifier conversion speed, the phenomenon of time delay in neural network systems is inevitable. Research indicates that the presence of time delay is a significant factor contributing to complex dynamic behaviors such as system instability and chaos [8]. To enhance the versatility and efficiency of the BAMNNs, Ding and Huang [9] developed a novel BAMNNs model in 2006. The focus of their investigation was on analyzing the global exponential stability of the equilibrium point and studying its characteristics in this model. The time delay BAMNNs has been advanced by its positive impact on its development [10,11,12,13].

Fractional calculus extends the traditional differentiation and integration operations to non-integer orders [14], and has been introduced into neural networks to capture the characteristics of memory and inheritance [15,16,17]. The emergence of fractional-order calculus has spurred the development of neural networks [6,7,18,19], which have found applications in diverse areas, including signal detection, fault diagnosis, optimization analysis, associative memory, and risk assessment. The fractional-order memristive neural networks (FMNN) are the specific type of fractional-order neural networks, which have been widely studied for their stability properties. For instance, scholars have investigated the asymptotic stability of FMNNs with delay using Caputo fractional differentiation and Filippov solution properties [20], and also have investigated the asymptotic stability of FMNNs with delay by leveraging the properties of Filippov solutions and Leibniz theorem [21].

As one of the significant research directions in the nonlinear systems field, synchronization includes quasi-consistent synchronization [22], projective synchronization [23], full synchronization [24], Mittag-Leffler synchronization [25], global synchronization [26] and many other types. And it is widely used in cryptography [27], image encryption [28], secure communication [29]. In engineering applications, people want to realize synchronization as soon as possible, so the concept of finite time synchronization is proposed. Due to its ability to achieve faster convergence speed in network systems, finite-time synchronization has become a crucial aspect in developing effective control strategies for realizing system stability or synchronization [30,31].

This paper addresses the challenge of achieving finite-time synchronization. The definition of finite-time synchronization used in this article is that the synchronization error be kept within a certain range for a limited time interval. However, dealing with the time delay term in this context is challenging. Previous studies have utilized the Hölder inequality [32] and the generalized Gronwall inequality [33,34] to address the finite-time synchronization problem of fractional-order time delay neural networks, providing valuable insights into the problem. However, this paper proposes a new criterion based on the quadratic fractional-order Gronwall inequality with time delay and comparison principle, offering a fresh perspective on the problem.

This paper presents significant contributions towards the study of finite-time synchronization in fractional-order stochastic MBAMNNs with time delay. The key contributions are as follows:

(1) We improved Lemma 2 in [34] by deriving a quadratic fractional-order Gronwall inequality with time delay, which is a crucial tool for analyzing the finite-time synchronization problem in stochastic neural networks.

(2) A novel criterion for achieving finite-time synchronization is proposed, which allows for the computation of the required synchronization time T. This criterion provides a new approach to analyze finite-time synchronization and has the potential to be widely applicable in the field of neural networks research. The paper is structured as follows: Section 2 introduces relevant concepts and presents the neural networks model used in this study. Section 3 proposes a novel quadratic fractional-order Gronwall inequality that takes into account time delay. This inequality is useful for studying the finite-time synchronization problem in fractional-order stochastic MBAMNNs with time delay, and by utilizing differential inclusion and set-valued mapping theory, a new criterion for determining the required time T for finite-time synchronization is derived. Section 4 provides a numerical example that demonstrates the effectiveness of the proposed results. Finally, suggestions for future research are presented.

2. Preliminaries and model

This section provides an overview of the necessary preliminaries related to fractional-order derivatives and the model of fractional-order stochastic MBAMNNs. We begin by introducing the fundamental concepts related to fractional-order derivatives and then move on to describe the fractional-order stochastic MBAMNNs model. And the definition of finite-time synchronization is provided.

2.1. Preliminaries

Notations: The norm and absolute value of vectors and matrices are defined as follows. Let N and $\mathbb{R}$ denote the sets of positive integers and real numbers, respectively.

The norm of a vector $e_x\left(\mathfrak{t}\right)=(e_{x,1}\left(\mathfrak{t}\right),e_{x,2}\left(\mathfrak{t}\right),\cdots ,e_{x,\mathfrak{n}}\left(\mathfrak{t}\right))\in C([0,+\infty ),{\mathbb{R}}^{\mathfrak{n}})$ is given by $\Vert e_x\left(\mathfrak{t}\right)\Vert ={\sum }_{\kappa =1}^{\mathfrak{n}}\mid e_{x,\kappa }\left(\mathfrak{t}\right)\mid$. Similarly, the norm of a vector $e_y\left(\mathfrak{t}\right)=(e_{y,1}\left(\mathfrak{t}\right),e_{y,2}\left(\mathfrak{t}\right),\cdots ,$$e_{y,\mathfrak{m}}\left(\mathfrak{t}\right))\in C([0,+\infty ),{\mathbb{R}}^{\mathfrak{m}})$ is defined as $\Vert e_y\left(\mathfrak{t}\right)\Vert ={\sum }_{\iota =1}^{\mathfrak{m}}\mid e_{y,\iota }\left(\mathfrak{t}\right)\mid ,\mathfrak{n},\mathfrak{m}\in N$.The induced norm of a matrix $\mathfrak{A}$ is denoted by $\Vert \mathfrak{A}\Vert =\underset{1\le \iota \le \mathfrak{n}}{max}{\sum }_{\kappa =1}^{\mathfrak{n}}\mid {\mathfrak{a}}_{\kappa \iota }\mid$.The absolute value of a vector $x\left(\mathfrak{t}\right)\in {\mathbb{R}}^{\mathfrak{n}}$ is defined as $\mid x\left(\mathfrak{t}\right)\mid ={(\mid x_1\left(\mathfrak{t}\right)\mid ,\mid x_2\left(\mathfrak{t}\right)\mid ,\cdots ,\mid x_{\mathfrak{n}}\left(\mathfrak{t}\right)\mid )}^T$.

Following that, we provide a review and introduction of several definitions and lemmas related to fractional calculus.

Definition 1.

[35] A fractional-order integral of a function $\kappa \left(\mathfrak{t}\right)$ with order α can be defined as

$$I_{{\mathfrak{t}}_0}^{\alpha }\kappa \left(\mathfrak{t}\right)=\frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{\left(\mathfrak{t}-\theta \right)}^{\alpha -1}\kappa \left(\theta \right)d\theta ,$$

where $\alpha >0,\mathfrak{t}\ge {\mathfrak{t}}_0$ and $\Gamma \left(\alpha \right)=\in {\mathfrak{t}}_0^{\infty }{\mathfrak{t}}^{\alpha -1}{e}^{-\mathfrak{t}}d\mathfrak{t}$.

Especially, $I_{{\mathfrak{t}}_0}^{\alpha }\kappa \left(\mathfrak{t}\right)=\kappa \left(\mathfrak{t}\right)-\kappa \left({\mathfrak{t}}_0\right)$ for $0<\alpha <1.$

Definition 2.

[35] Suppose $\kappa \left(\mathfrak{t}\right)\in C^{\iota }([0,+\infty ),\mathbb{R})$, where ι is a positive integer. The Caputo derivative of order α of the function $\kappa \left(\mathfrak{t}\right)$ can be expressed as:

$${}^c{D}_0^{\alpha }\kappa \left(\mathfrak{t}\right)=\in {\mathfrak{t}}_0^{\mathfrak{t}}\frac{{(\mathfrak{t}-\theta )}^{-\alpha }}{\Gamma (\iota -\alpha )}{\kappa }^{\left(\iota \right)}\left(\theta \right)d\theta ,$$

where $\kappa -1<\alpha <\kappa ,\kappa \in N.$

For the convenience, we use $D^{\alpha }$ to represent ${}^c{D}_0^{\alpha }.$

Lemma 1.

[36] Let $\mathfrak{t}\in [{\mathfrak{t}}_0,T]$, $u\left(\mathfrak{t}\right),v\left(\mathfrak{t}\right)\in \mathbb{R},\kappa \left(\mathfrak{t}\right)\in C^1([{\mathfrak{t}}_0,T],\mathbb{R})$ and satisfies

$$\left\{\begin{array}{cc}{\kappa }^{{}^{\prime }}\left(\mathfrak{t}\right)\le u\left(\mathfrak{t}\right)\kappa \left(\mathfrak{t}\right)+v\left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0,T],\hfill \\ \kappa \left({\mathfrak{t}}_0\right)\le {\kappa }_0,\hfill & {\kappa }_0\in \mathbb{R}.\hfill \end{array}\right.$$

Then, we have

$$\kappa \left(\mathfrak{t}\right)\le {\kappa }_{0}exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}u\left(\theta \right)d\theta \right)+{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}exp\left({\int }_{\theta }^{\mathfrak{t}}u\left(s\right)ds\right)v\left(\theta \right)d\theta .$$

2.2. Model

We investigate a kind of fractional-order differential equations that capture the dynamics of fractional-order stochastic MBAMNNs with time delays. These equations are viewed as the driving system (1), which models the interactions between neurons in MBAMNNs and accounts for the influence of discontinuous jumps and time delays. Through examining the stability and analytical solutions of these equations, this study aims to enhance the comprehension of the behavior of MBAMNNs, ultimately leading to more comprehensive analysis for practical applications of this model.

$$\left\{\begin{array}{cc}{D}^{\alpha }{x}_{\kappa }\left(\mathfrak{t}\right)=\hfill & -u_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)x_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{a}}_{\iota \kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{b}}_{\iota \kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)\hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)+I_{\kappa },\kappa =1,2,\dots ,\mathfrak{n},\hfill \\ D^{\alpha }{y}_{\iota }\left(\mathfrak{t}\right)=\hfill & -v_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)y_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{c}}_{\kappa \iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{d}}_{\kappa \iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota },\iota =1,2,\dots ,\mathfrak{m},\hfill \end{array}\right.$$

(1)

where $0<\alpha <1,t\in [0,+\infty ).$

In this system, the positive parameters ${\mu }_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)$ and $v_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)$ represent the rates of neuron self-inhibition, while $x_{\kappa }\left(\mathfrak{t}\right)$ and $y_{\iota }\left(\mathfrak{t}\right)$ denote the state variables of the $\kappa$-th and $\iota$-th neuron, respectively. The activation functions without time delay are denoted by ${\mathfrak{f}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)$ and ${\mathfrak{g}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)$, and those with time delay are denoted by ${\mathfrak{f}}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)$ and ${\mathfrak{g}}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)$. The neural connection memristive weights matrices are represented by ${\mathfrak{a}}_{\iota \kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)$, ${\mathfrak{b}}_{\iota \kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)$, ${\mathfrak{c}}_{\kappa \iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)$, and ${\mathfrak{d}}_{\kappa \iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)$. Stochastic terms representing Brownian motion are denoted by $r_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)$ and $s_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)$. The constant input vectors are represented by $I_{\kappa }$ and $J_{\iota }$. The time delay parameters ${\tau }_1$ and ${\tau }_2$ satisfy $0\le {\tau }_1\le \tau$ and $0\le {\tau }_2\le \tau$, where $\tau$ is a constant.

The initial conditions of fractional-order stochastic MBAMNNs (1) are given by $x\left(\mathfrak{t}\right)={(x_1\left(\mathfrak{t}\right),x_2\left(\mathfrak{t}\right),\cdots ,x_{\mathfrak{n}}\left(\mathfrak{t}\right))}^T,$$y\left(\mathfrak{t}\right)={(y_1\left(\mathfrak{t}\right),y_2\left(\mathfrak{t}\right),\cdots ,y_{\mathfrak{m}}\left(\mathfrak{t}\right))}^T$, where $x_{\kappa }\left(r\right)={\varphi }_{\kappa }\left(r\right)\in C([-{\tau }_1,0],$${\mathbb{R}}^{\mathfrak{n}})$ and $y_{\iota }\left(r\right)={\psi }_{\iota }\left(r\right)\in C([-{\tau }_2,0],{\mathbb{R}}^{\mathfrak{m}})$. Here, ${\varphi }_{\kappa }\left(r\right)$ and ${\psi }_{\iota }\left(r\right)$ are continuous functions on $[-{\tau }_1,0]$ and $[-{\tau }_2,0]$.

Then, the corresponding system of drive system (1) is gived by

$$\left\{\begin{array}{cc}{D}^{\alpha }{\breve{x}}_{\kappa }\left(\mathfrak{t}\right)=\hfill & -u_{\kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right){\breve{x}}_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{a}}_{\iota \kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right){\mathfrak{f}}_{\iota }\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{b}}_{\iota \kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right){\mathfrak{f}}_{\iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right)\hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)+I_{\kappa }+{\gamma }_{\kappa },\kappa =1,2,\dots ,\mathfrak{n},\hfill \\ D^{\alpha }{\breve{y}}_{\iota }\left(\mathfrak{t}\right)=\hfill & -v_j\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right){\breve{y}}_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{c}}_{\kappa \iota }\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right){\mathfrak{g}}_{\kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{d}}_{\kappa \iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right){\mathfrak{g}}_{\kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota }+{\eta }_{\iota },\iota =1,2,\dots ,\mathfrak{m},\hfill \end{array}\right.$$

(2)

The initial conditions of the corresponding system (2) are ${\breve{x}}_{\kappa }\left(r\right)={\breve{\varphi }}_{\kappa }\left(r\right)\in C([-{\tau }_1,0],{\mathbb{R}}^{\mathfrak{n}}),{\breve{y}}_{\iota }\left(r\right)={\breve{\psi }}_{\iota }\left(r\right)\in C([-{\tau }_2,0],{\mathbb{R}}^{\mathfrak{m}})$; ${\gamma }_{\kappa }$ and ${\eta }_{\iota }$ are the following controllers:

$$\left\{\begin{array}{c}{\gamma }_{\kappa }=-{\zeta }_{\kappa }({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)-x_{\kappa }\left(\mathfrak{t}\right)),\hfill \\ {\eta }_{\iota }=-{\zeta }_{\iota }({\breve{y}}_{\iota }\left(\mathfrak{t}\right)-y_{\iota }\left(\mathfrak{t}\right)),\hfill \end{array}\right.$$

where ${\zeta }_{\kappa }$ and ${\zeta }_{\iota }$ are both positive numbers called the control gain.

The synchronization error, as defined by systems (1) and (2), can be expressed as:

$$\left\{\begin{array}{cc}{e}_{x,\kappa }\left(\mathfrak{t}\right)={\breve{x}}_{\kappa }\left(\mathfrak{t}\right)-x_{\kappa }\left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0,T],\hfill \\ e_{y,\iota }\left(\mathfrak{t}\right)={\breve{y}}_{\iota }\left(\mathfrak{t}\right)-y_{\iota }\left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0,T],\hfill \\ {\phi }_{x,\kappa }\left(r\right)={\breve{\varphi }}_{\kappa }\left(r\right)-{\varphi }_{\kappa }\left(r\right),\hfill & r\in [{\mathfrak{t}}_0-{\tau }_1,{\mathfrak{t}}_0],\hfill \\ {\breve{\phi }}_{y,\iota }\left(r\right)={\breve{\psi }}_{\iota }\left(r\right)-{\psi }_{\iota }\left(r\right),\hfill & r\in [{\mathfrak{t}}_0-{\tau }_2,{\mathfrak{t}}_0].\hfill \end{array}\right.$$

Then we obtain synchronization error is $\Vert e_x\left(\mathfrak{t}\right)\Vert +\Vert e_y\left(\mathfrak{t}\right)\Vert ;$ And denote $\phi =\underset{s\in [-{\tau }_1,0]}{sup}\Vert {\phi }_{x,\kappa }\left(r\right)\Vert ,\breve{\phi }=\underset{s\in [-{\tau }_2,0]}{sup}\Vert {\breve{\phi }}_{y,\iota }\left(r\right)\Vert .$

Definition 3.

If there is a real number $T>0$ such that for any $\mathfrak{t}>T$ and $ϵ>0$, the synchronization error satisfies

$$\Vert e_x\left(\mathfrak{t}\right)\Vert +\Vert e_y\left(\mathfrak{t}\right)\Vert \le ϵ,T<\mathfrak{t}<+\infty$$

then it can be inferred that the drive system (1) and the response system (2) achieve finite-time synchronization at T.

Remark 1.

To obtain a sufficient condition for achieving finite-time synchronization between systems (1) and (2), it is necessary to identify an appropriate evaluation function $\chi \left(\mathfrak{t}\right)$. This function should satisfy $\Vert e_x\left(\mathfrak{t}\right)\Vert +\Vert e_y\left(\mathfrak{t}\right)\Vert \le \chi \left(\mathfrak{t}\right)\le ϵ.$

The next section will apply Lemma 1, which provides a quadratic fractional Gronwall inequality with time delay. This inequality is used to analyze the behavior of the MBAMNNs system.

3. Main results

This section presents a novel approach to obtain the evaluation function ${\chi }_1\left(\mathfrak{t}\right)$ by improving quadratic fractional Gronwall inequality with time delay. We then utilize Theorem 1 to convert this inequality into a format that is consistent with Lemma 2, enabling us to derive a synchronization criterion for the drive system (1) and corresponding system (2) in finite time. Specifically, the application of Lemma 2 leads to the novel criterion for finite-time synchronization.

Quadratic fractional Gronwall inequality with time delay is given below.

Lemma 2.

Let $T\in (0,+\infty )$ and $\omega \left(\mathfrak{t}\right)$, $\mathfrak{a}\left(\mathfrak{t}\right)$, $\mathfrak{b}\left(\mathfrak{t}\right)$, $\sigma \left(\mathfrak{t}\right)$ and $\upsilon \left(\mathfrak{t}\right)$ be continuous functions that are nonnegative and defined on $[{\mathfrak{t}}_0,T]$. Let $\varphi \left(\mathfrak{t}\right)$ be a nonnegative continuous function defined on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and suppose

$$\left\{\begin{array}{cc}{\omega }^2\left(\mathfrak{t}\right)\le {\upsilon }^2\left(\mathfrak{t}\right){\left\{{\int }_{{\mathfrak{t}}_0}^T{[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]}^{p}d\theta \right\}}^{\frac{1}{p}}+2{\sigma }^2\left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0,T],\hfill \\ \omega \left(\mathfrak{t}\right)\le \varphi \left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0].\hfill \end{array}\right.$$

(3)

Assume $\sigma \left(\mathfrak{t}\right)$ and $\upsilon \left(\mathfrak{t}\right)$ are nondecreasing on $[{\mathfrak{t}}_0,T]$, $\varphi \left(\mathfrak{t}\right)$ is nondecreasing on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and $\varphi \left({\mathfrak{t}}_0\right)=\sigma \left({\mathfrak{t}}_0\right).$

(1) If $\sigma \left(\mathfrak{t}\right)>\upsilon \left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\sigma }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}}.\hfill \end{array}$$

(2) If $\sigma \left(\mathfrak{t}\right)\le \upsilon \left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\upsilon }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}},\hfill \end{array}$$

where $\mathfrak{t}\in [{\mathfrak{t}}_0,T],\Omega >0,p\ge 1$ and they are all constants.

Proof 

(Proof). Define a function $d\left(\mathfrak{t}\right)$ by

$$d\left(\mathfrak{t}\right)={\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]}^{p}d\theta ,\mathfrak{t}\in [{\mathfrak{t}}_0,T].$$

As shown in Equation (3), we get

$${\omega }^2\left(\mathfrak{t}\right)\le 2{\sigma }^2\left(\mathfrak{t}\right)+{\upsilon }^2\left(\mathfrak{t}\right)d^{\frac{1}{p}}\left(\mathfrak{t}\right),\mathfrak{t}\in [{\mathfrak{t}}_0,T].$$

(4)

When $\mathfrak{t}\in [{\mathfrak{t}}_0,{\mathfrak{t}}_0+\Omega ]$ and the function $d\left({\mathfrak{t}}_0\right)=0$, then

$$\begin{array}{cc}\hfill d^{{}^{\prime }}\left(\mathfrak{t}\right)& ={[\mathfrak{a}\left(\mathfrak{t}\right){\omega }^2\left(\mathfrak{t}\right)+\mathfrak{b}\left(\mathfrak{t}\right){\omega }^2(\mathfrak{t}-\Omega )]}^p\hfill \\ \hfill & \le {\{\mathfrak{a}\left(\mathfrak{t}\right)[2{\sigma }^2\left(\mathfrak{t}\right)+{\upsilon }^2\left(\mathfrak{t}\right)d^{\frac{1}{p}}\left(\mathfrak{t}\right)]+\mathfrak{b}\left(\mathfrak{t}\right){\varphi }^2(\mathfrak{t}-\Omega )\}}^p\hfill \\ \hfill & \le 2^{p-1}\left\{{[2\mathfrak{a}\left(\mathfrak{t}\right){\sigma }^2\left(\mathfrak{t}\right)+\mathfrak{b}\left(\mathfrak{t}\right){\varphi }^2(\mathfrak{t}-\Omega )]}^p+{\mathfrak{a}}^p\left(\mathfrak{t}\right){\upsilon }^{2p}\left(\mathfrak{t}\right)d\left(\mathfrak{t}\right)\right\}\hfill \\ \hfill & =2^{p-1}{[2\mathfrak{a}\left(\mathfrak{t}\right){\sigma }^2\left(\mathfrak{t}\right)+\mathfrak{b}\left(\mathfrak{t}\right){\varphi }^2(\mathfrak{t}-\Omega )]}^p+2^{p-1}{\mathfrak{a}}^p\left(\mathfrak{t}\right){\upsilon }^{2p}\left(\mathfrak{t}\right)d\left(\mathfrak{t}\right).\hfill \end{array}$$

From Lemma 1, we obtain

$$\begin{array}{cc}\hfill d\left(\mathfrak{t}\right)\le & {\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^T{2}^{p-1}{\mathfrak{a}}^p\left(s\right){\upsilon }^{2p}\left(s\right)ds\right)d\theta .\hfill \end{array}$$

(5)

By using inequalities (4) and (5), we have

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^t{2}^{p-1}{\mathfrak{a}}^p\left(s\right){\upsilon }^{2p}\left(s\right)ds\right)d\theta \right\}}^{\frac{1}{p}}\hfill \\ \hfill & +2{\sigma }^2\left(\mathfrak{t}\right).\hfill \end{array}$$

Assume $\sigma \left(\mathfrak{t}\right),\upsilon \left(\mathfrak{t}\right)$ are nondecreasing on $[{\mathfrak{t}}_0,{\mathfrak{t}}_0+\Omega ]$, $\varphi \left(\mathfrak{t}\right)$ is nondecreasing on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and $\sigma \left({\mathfrak{t}}_0\right)=\varphi \left({\mathfrak{t}}_0\right)$.

(1) If $\sigma \left(\mathfrak{t}\right)>\upsilon \left(\mathfrak{t}\right)$, then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^{\mathfrak{t}}{2}^{p-1}{\mathfrak{a}}^p\left(s\right){\upsilon }^{2p}\left(s\right)ds\right)d\theta \right\}}^{\frac{1}{p}}\hfill \\ \hfill & +2{\sigma }^2\left(\mathfrak{t}\right)\hfill \\ \hfill \le & {\sigma }^2\left(\mathfrak{t}\right){\left\{{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}exp\left({\int }_{\theta }^{\mathfrak{t}}{2}^{p-1}{[2\mathfrak{a}\left(s\right)+\mathfrak{b}\left(s\right)]}^p{\sigma }^{2p}\left(\mathfrak{t}\right)ds\right)d\theta \right\}}^{\frac{1}{p}}\hfill \\ \hfill & +2{\sigma }^2\left(\mathfrak{t}\right)\hfill \\ \hfill \le & {\sigma }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+2\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}}.\hfill \end{array}$$

(2) If $\sigma \left(\mathfrak{t}\right)\le \upsilon \left(\mathfrak{t}\right)$, similar to case (1), we obtain

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\upsilon }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+2\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}}.\hfill \end{array}$$

For $\mathfrak{t}\in [{\mathfrak{t}}_0+\Omega ,T]$, we get

$$\begin{array}{cc}\hfill d^{{}^{\prime }}\left(\mathfrak{t}\right)& ={[\mathfrak{a}\left(\mathfrak{t}\right){\omega }^2\left(\mathfrak{t}\right)+\mathfrak{b}\left(\mathfrak{t}\right){\omega }^2(\mathfrak{t}-\Omega )]}^p\hfill \\ \hfill & \le {\{\mathfrak{a}\left(\mathfrak{t}\right)[2{\sigma }^2\left(\mathfrak{t}\right)+d^{\frac{1}{p}}\left(\mathfrak{t}\right){\upsilon }^2\left(\mathfrak{t}\right)]+\mathfrak{b}\left(\mathfrak{t}\right)[2{\sigma }^2(\mathfrak{t}-\Omega )+d^{\frac{1}{p}}(\mathfrak{t}-\Omega ){\sigma }^2(t-\Omega )]\}}^p\hfill \\ \hfill & \le {\{[{\upsilon }^2\left(\mathfrak{t}\right)\mathfrak{a}\left(\mathfrak{t}\right)+{\sigma }^2(\mathfrak{t}-\Omega )\mathfrak{b}\left(\mathfrak{t}\right)]d^{\frac{1}{p}}\left(\mathfrak{t}\right)+2{\sigma }^2\left(\mathfrak{t}\right)\mathfrak{a}\left(\mathfrak{t}\right)+2{\sigma }^2(\mathfrak{t}-\Omega )\mathfrak{b}\left(\mathfrak{t}\right)\}}^p\hfill \\ \hfill & \le 2^{p-1}{[{\upsilon }^2\left(\mathfrak{t}\right)\mathfrak{a}\left(\mathfrak{t}\right)+{\sigma }^2(\mathfrak{t}-\Omega )\mathfrak{b}\left(\mathfrak{t}\right)]}^{p}d\left(\mathfrak{t}\right)+2^p{[{\sigma }^2\left(\mathfrak{t}\right)\mathfrak{a}\left(\mathfrak{t}\right)+{\sigma }^2(\mathfrak{t}-\Omega )\mathfrak{b}\left(\mathfrak{t}\right)]}^p.\hfill \end{array}$$

Then, by utilizing Lemma 1 and inequality (5), we arrive at

$$\begin{array}{cc}\hfill d\left(\mathfrak{t}\right)\le & {\int }_{{\mathfrak{t}}_0}^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{\mathfrak{a}}^p\left(s\right){\upsilon }^{2p}\left(s\right)ds\right)d\theta \hfill \\ \hfill & \times exp\left({\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}\times 2^{p-1}{[\mathfrak{a}\left(\theta \right){\upsilon }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^{p}d\theta \right)\hfill \\ & +{\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^p{[\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^{\mathfrak{t}}{2}^{p-1}{[\mathfrak{a}\left(s\right){\upsilon }^2\left(s\right)+\mathfrak{b}\left(s\right){\upsilon }^2(s-\Omega )]}^{p}ds\right)d\theta .\hfill \end{array}$$

By using (4), we have

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right)\left\{{\int }_{{\mathfrak{t}}_0}^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp({\int }_{\theta }^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{\mathfrak{a}}^p\left(s\right)\times {\upsilon }^{2p}\left(s\right)ds)d\theta \right.\hfill \\ \hfill & \times exp\left({\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^{p-1}{[\mathfrak{a}\left(\theta \right){\upsilon }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^{p}d\theta \right)\hfill \\ \hfill & +{\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^p{[\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^p\hfill \\ \hfill & {\left.\times exp\left({\int }_{\theta }^{\mathfrak{t}}{2}^{p-1}{[\mathfrak{a}\left(s\right){\upsilon }^2\left(s\right)+\mathfrak{b}\left(s\right){\upsilon }^2(s-\Omega )]}^{p}ds\right)d\theta \right\}}^{\frac{1}{p}}+2{\sigma }^2\left(\mathfrak{t}\right).\hfill \end{array}$$

Similarly, assume $\sigma \left(\mathfrak{t}\right),\upsilon \left(\mathfrak{t}\right)$ are nondecreasing on $[{\mathfrak{t}}_0+\Omega ,T]$, $\varphi \left(\mathfrak{t}\right)$ is nondecreasing on $t\in [{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and $\sigma \left({\mathfrak{t}}_0\right)=\varphi \left({\mathfrak{t}}_0\right)$.

(1) If $\sigma \left(\mathfrak{t}\right)>\upsilon \left(\mathfrak{t}\right)$, then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right)\left\{{\int }_{{\mathfrak{t}}_0}^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{[2\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\varphi }^2(\theta -\Omega )]}^{p}exp\left({\int }_{\theta }^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{\mathfrak{a}}^p\left(s\right){\upsilon }^{2p}\left(s\right)ds\right)d\theta \right.\hfill \\ \hfill & \times exp\left({\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^{p-1}{[\mathfrak{a}\left(\theta \right){\upsilon }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^{p}d\theta \right)\hfill \\ \hfill & +{\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^p{[\mathfrak{a}\left(\theta \right){\sigma }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\sigma }^2(\theta -\Omega )]}^p\hfill \\ \hfill & {\left.\times exp\left({\int }_{\theta }^{\mathfrak{t}}{2}^{p-1}{[\mathfrak{a}\left(s\right){\upsilon }^2\left(s\right)+\mathfrak{b}\left(s\right){\upsilon }^2(s-\Omega )]}^{p}ds\right)d\theta \right\}}^{\frac{1}{p}}+2{\sigma }^2\left(\mathfrak{t}\right)\hfill \\ \hfill \le & {\sigma }^2\left(\mathfrak{t}\right)\left\{exp[\left({\int }_{{\mathfrak{t}}_0}^{{\mathfrak{t}}_0+\Omega }{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{b}\left(\theta \right)+\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1]\right.\hfill \\ \hfill & \times exp\left({\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)\hfill \\ \hfill & {\left.+exp\left({\int }_{{\mathfrak{t}}_0+\Omega }^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1\right\}}^{\frac{1}{p}}+2{\sigma }^2\left(\mathfrak{t}\right)\hfill \\ \hfill \le & {\sigma }^2\left(\mathfrak{t}\right){\left\{exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+2\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1\right\}}^{\frac{1}{p}}+2{\sigma }^2\left(\mathfrak{t}\right).\hfill \end{array}\end{array}$$

(2) If $\sigma \left(\mathfrak{t}\right)\le \upsilon \left(\mathfrak{t}\right)$, then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\upsilon }^{2p}\left(\mathfrak{t}\right){[\mathfrak{b}\left(\theta \right)+2\mathfrak{a}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}}.\hfill \end{array}$$

Based on the above analysis, from Lemma 2 and when $\sigma \left(\mathfrak{t}\right)$ and $\upsilon \left(\mathfrak{t}\right)$ are non-decreasing functions for $\mathfrak{t}\in [{\mathfrak{t}}_0,T]$, $\varphi \left(\mathfrak{t}\right)$ is a non-decreasing function for $t\in [{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ with $\varphi \left({\mathfrak{t}}_0\right)=\sigma \left({\mathfrak{t}}_0\right),$ we get the following results.

(1) If $\sigma \left(\mathfrak{t}\right)>\upsilon \left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\sigma }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}},\mathfrak{t}\in [{\mathfrak{t}}_0,T].\hfill \end{array}$$

(2) If $\sigma \left(\mathfrak{t}\right)\le \upsilon \left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\upsilon }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\upsilon }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}},\mathfrak{t}\in [{\mathfrak{t}}_0,T].\hfill \end{array}$$

□

Theorem 1.

Assume $T\in (0,+\infty )$ and non-negative continuous functions $\omega \left(\mathfrak{t}\right)$, $\sigma \left(\mathfrak{t}\right)$, $\upsilon \left(\mathfrak{t}\right)$, $\mathfrak{a}\left(\mathfrak{t}\right)$, and $\mathfrak{b}\left(\mathfrak{t}\right)$ defined on $[{\mathfrak{t}}_0,T]$. Let $\varphi \left(\mathfrak{t}\right)$ be a non-negative continuous function defined on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and suppose

$$\left\{\begin{array}{cc}{\omega }^2\left(\mathfrak{t}\right)\le 2{\sigma }^2\left(\mathfrak{t}\right)+{\displaystyle \frac{\upsilon \left(\mathfrak{t}\right)}{{\Gamma }^2\left(\alpha \right)}}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2(\alpha -1)}[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]d\theta ,\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0,T],\hfill \\ \omega \left(\mathfrak{t}\right)\le \varphi \left(\mathfrak{t}\right),\hfill & \mathfrak{t}\in [{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0].\hfill \end{array}\right.$$

Assume $\sigma \left(\mathfrak{t}\right)$ and $\upsilon \left(\mathfrak{t}\right)$ are nondecreasing functions on $[{\mathfrak{t}}_0,T]$, and $\varphi \left(\mathfrak{t}\right)$ is a nondecreasing function on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ with $\varphi \left({\mathfrak{t}}_0\right)=\sigma \left({\mathfrak{t}}_0\right)$.

(1) If $\sigma \left(\mathfrak{t}\right)>H\left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & {\sigma }^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\sigma }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}}.\hfill \end{array}$$

(2) If $\sigma \left(\mathfrak{t}\right)\le H\left(\mathfrak{t}\right),$ then

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & H^2\left(\mathfrak{t}\right){\left\{2+[exp\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{2}^{p-1}{\upsilon }^{2p}\left(\mathfrak{t}\right){[2\mathfrak{a}\left(\theta \right)+\mathfrak{b}\left(\theta \right)]}^{p}d\theta \right)-1]\right\}}^{\frac{1}{p}},\hfill \end{array}$$

where $\mathfrak{t}\in [{\mathfrak{t}}_0,T]$ and $p,q>0$ such that $\frac{1}{q}+\frac{1}{p}=1,\alpha >\frac{1+p}{2p}$ and $H^2\left(\mathfrak{t}\right)=\frac{\upsilon \left(\mathfrak{t}\right){(\mathfrak{t}-{\mathfrak{t}}_0)}^{2(\alpha -1)+\frac{1}{q}}}{{\Gamma }^2\left(\alpha \right){(2q(\alpha -1)+1)}^{\frac{1}{q}}}.$

Proof 

(Proof). By using the Hölder inequality, it follows that

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le & 2{\sigma }^2\left(\mathfrak{t}\right)+{\displaystyle \frac{\upsilon \left(\mathfrak{t}\right)}{{\Gamma }^2\left(\alpha \right)}}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2(\alpha -1)}[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]d\theta \hfill \\ \hfill \le & 2{\sigma }^2\left(\mathfrak{t}\right)+{\displaystyle \frac{\upsilon \left(\mathfrak{t}\right)}{{\Gamma }^2\left(\alpha \right)}}{\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2q(\alpha -1)}d\theta \right)}^{\frac{1}{q}}{\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]}^{p}d\theta \right)}^{\frac{1}{p}}\hfill \\ \hfill \le & 2{\sigma }^2\left(\mathfrak{t}\right)+{\displaystyle \frac{\upsilon \left(\mathfrak{t}\right)}{{\Gamma }^2\left(\alpha \right)}}{\displaystyle \frac{{(\mathfrak{t}-{\mathfrak{t}}_0)}^{2(\alpha -1)+\frac{1}{q}}}{{[2q(\alpha -1)+1]}^{\frac{1}{q}}}}{\left({\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]}^{p}d\theta \right)}^{\frac{1}{p}}.\hfill \end{array}$$

Let $\alpha >\frac{1+p}{2p},H^2\left(\mathfrak{t}\right)=\frac{\upsilon \left(\mathfrak{t}\right){(\mathfrak{t}-{\mathfrak{t}}_0)}^{2(\alpha -1)+\frac{1}{q}}}{{\Gamma }^2\left(\alpha \right){(2q(\alpha -1)+1)}^{\frac{1}{q}}},$ then from Lemma 2 we obtain

$$\begin{array}{cc}\hfill {\omega }^2\left(\mathfrak{t}\right)\le 2{\sigma }^2\left(\mathfrak{t}\right)+& H^2\left(\mathfrak{t}\right){\left\{{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{[\mathfrak{a}\left(\theta \right){\omega }^2\left(\theta \right)+\mathfrak{b}\left(\theta \right){\omega }^2(\theta -\Omega )]}^{p}d\theta \right\}}^{\frac{1}{p}}.\hfill \end{array}$$

The proof is completed. □

To analyze the solutions of the discontinuous systems represented by equations (1) and (2), Filippov regularization is used. This involves transforming the equations into differential inclusions and set-valued maps.

The drive system represented by equation (1) can be expressed in terms of a differential inclusion, which is a powerful tool in the theory of differential inclusions. By using this approach, we can study the behavior of the system even when it experiences discontinuities or impulses.

Overall, Filippov regularization allows us to analyze the solutions of discontinuous systems like Equations (1) and (2) in a rigorous and systematic way, providing insights into their behavior and enabling us to make informed decisions about their design and operation.

$$\left\{\begin{array}{cc}{D}^{\alpha }{x}_{\kappa }\left(\mathfrak{t}\right)\in \hfill & -co\left[u_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)\right]x_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}co\left[{\mathfrak{a}}_{\iota \kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)\right]{\mathfrak{f}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}co\left[{\mathfrak{b}}_{\iota \kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\right]\hfill \\ \hfill & \times {\mathfrak{f}}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)+\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)+I_{\kappa },\hfill \\ D^{\alpha }{y}_{\iota }\left(\mathfrak{t}\right)\in \hfill & -co\left[v_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)\right]y_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}co\left[{\mathfrak{c}}_{\kappa \iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)\right]{\mathfrak{g}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}co\left[{\mathfrak{d}}_{\kappa \iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)\right]\hfill \\ \hfill & \times {\mathfrak{g}}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota }.\hfill \end{array}\right.$$

According to the definition of set-valued maps, we obtain

$$\begin{array}{ccccccccccc}& & & co\left[u_{\kappa }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{u}}_{\kappa },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{u}}_{\kappa },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{u}}_{\kappa },{\ddot{u}}_{\kappa }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.& & & co\left[{\mathfrak{a}}_{\iota \kappa }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{\mathfrak{a}}}_{\iota \kappa },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{\mathfrak{a}}}_{\iota \kappa },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{\mathfrak{a}}}_{\iota \kappa },{\ddot{\mathfrak{a}}}_{\iota \kappa }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.& & \\ \\ & & & co\left[{\mathfrak{b}}_{\iota \kappa }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{\mathfrak{b}}}_{\iota \kappa },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{\mathfrak{b}}}_{\iota \kappa },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{\mathfrak{b}}}_{\iota \kappa },{\ddot{\mathfrak{b}}}_{\iota \kappa }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.& & & co\left[v_{\iota }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{v}}_{\iota },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{v}}_{\iota },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{v}}_{\iota },{\ddot{v}}_{\iota }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.& & \\ \\ & & & co\left[{\mathfrak{c}}_{\kappa \iota }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{\mathfrak{c}}}_{\kappa \iota },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{\mathfrak{c}}}_{\kappa \iota },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{\mathfrak{c}}}_{\kappa \iota },{\ddot{\mathfrak{c}}}_{\kappa \iota }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.& & & co\left[{\mathfrak{d}}_{\kappa \iota }\left(\varsigma \right)\right]=\hfill & \hfill \left\{\begin{array}{cccc}& {\dot{\mathfrak{d}}}_{\kappa \iota },\hfill & & \mid \varsigma \mid <\mathring{T},\hfill \\ & {\ddot{\mathfrak{d}}}_{\kappa \iota },\hfill & & \mid \varsigma \mid >\mathring{T},\hfill \\ & co\left\{{\dot{\mathfrak{d}}}_{\kappa \iota },{\ddot{\mathfrak{d}}}_{\kappa \iota }\right\},\hfill & & \mid \varsigma \mid =\mathring{T},\hfill \end{array}\right.\end{array}$$

where $\varsigma \in \mathbb{R},$ the switching jumps ${\mathring{T}}_{\kappa },T_{\iota }^{**}$ are positive contants, ${\dot{u}}_i$, ${\ddot{u}}_i$, ${\dot{\mathfrak{a}}}_{\iota \kappa }$, ${\ddot{\mathfrak{a}}}_{\iota \kappa }$, ${\dot{\mathfrak{b}}}_{\iota \kappa }$, ${\ddot{\mathfrak{b}}}_{\iota \kappa }$, ${\dot{v}}_j$, ${\ddot{v}}_j$, ${\dot{\mathfrak{c}}}_{\kappa \iota }$, ${\ddot{\mathfrak{c}}}_{\kappa \iota }$, ${\dot{\mathfrak{d}}}_{\kappa \iota }$ and ${\ddot{\mathfrak{d}}}_{\kappa \iota }$ are all contant numbers. $co\left[u_{\kappa }\left(\varsigma \right)\right]$, $co\left[{\mathfrak{a}}_{\iota \kappa }\left(\varsigma \right)\right]$, $co\left[{\mathfrak{b}}_{\iota \kappa }\left(\varsigma \right)\right]$, $co\left[v_j\left(\varsigma \right)\right]$, $co\left[{\mathfrak{c}}_{\kappa \iota }\left(\varsigma \right)\right]$, $co\left[{\mathfrak{d}}_{\kappa \iota }\left(\varsigma \right)\right]$ are all compact, closed and convex.

Let

$$\begin{array}{ccccccc}& {\acute{u}}_{\kappa }\left(\varsigma \right)\in co\left[u_{\kappa }\left(\varsigma \right)\right],\hfill & & {\acute{\mathfrak{a}}}_{\iota \kappa }\left(\varsigma \right)\in co\left[{\mathfrak{a}}_{\iota \kappa }\left(\varsigma \right)\right],\hfill & & & \hfill {\acute{\mathfrak{b}}}_{\iota \kappa }\left(\varsigma \right)\in co\left[{\mathfrak{b}}_{\iota \kappa }\left(\varsigma \right)\right],\\ & {\acute{v}}_{\iota }\left(\varsigma \right)\in co\left[v_{\iota }\left(\varsigma \right)\right],\hfill & & {\acute{\mathfrak{c}}}_{\kappa \iota }\left(\varsigma \right)\in co\left[{\mathfrak{c}}_{\kappa \iota }\left(\varsigma \right)\right],\hfill & & & \hfill {\acute{\mathfrak{d}}}_{\kappa \iota }\left(\varsigma \right)\in co\left[{\mathfrak{d}}_{\kappa \iota }\left(\varsigma \right)\right].\end{array}$$

By modifying the drive system (1), we can achieve

$$\left\{\begin{array}{cc}{D}^{\alpha }{x}_{\kappa }\left(\mathfrak{t}\right)=\hfill & -{\acute{u}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)x_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}{\acute{\mathfrak{a}}}_{\iota \kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}{\acute{\mathfrak{b}}}_{\iota \kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)\hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)+I_{\kappa },\hfill \\ D^{\alpha }{y}_{\iota }\left(\mathfrak{t}\right)=\hfill & -{\acute{v}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)y_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\acute{\mathfrak{c}}}_{\kappa \iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\acute{\mathfrak{d}}}_{\kappa \iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota }.\hfill \end{array}\right.$$

(6)

Similarly, let

$$\begin{array}{ccccccc}& {\stackrel{`}{u}}_{\kappa }\left(\varsigma \right)\in co\left[u_{\kappa }\left(\varsigma \right)\right],\hfill & & {\stackrel{`}{\mathfrak{a}}}_{\iota \kappa }\left(\varsigma \right)\in co\left[{\mathfrak{a}}_{\iota \kappa }\left(\varsigma \right)\right],\hfill & & & \hfill {\stackrel{`}{\mathfrak{b}}}_{\iota \kappa }\left(\varsigma \right)\in co\left[{\mathfrak{b}}_{\iota \kappa }\left(\varsigma \right)\right],\\ & {\stackrel{`}{v}}_{\iota }\left(\varsigma \right)\in co\left[v_{\iota }\left(\varsigma \right)\right],\hfill & & {\stackrel{`}{\mathfrak{c}}}_{\kappa \iota }\left(\varsigma \right)\in co\left[{\mathfrak{c}}_{\kappa \iota }\left(\varsigma \right)\right],\hfill & & & \hfill {\stackrel{`}{\mathfrak{d}}}_{\kappa \iota }\left(\varsigma \right)\in co\left[{\mathfrak{d}}_{\kappa \iota }\left(\varsigma \right)\right],\end{array}$$

by employing the similar method, we can modify the corresponding system (2) as follows:

$$\left\{\begin{array}{cc}{D}^{\alpha }{\breve{x}}_{\kappa }\left(\mathfrak{t}\right)=\hfill & -{\stackrel{`}{u}}_{\kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right){\breve{x}}_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}{\stackrel{`}{\mathfrak{a}}}_{\iota \kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right){\mathfrak{f}}_{\iota }\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}{\stackrel{`}{\mathfrak{b}}}_{\iota \kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right){\mathfrak{f}}_{\iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right)\hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right)dB\left(\mathfrak{t}\right)+I_{\kappa }+{\gamma }_{\kappa },\hfill \\ D^{\alpha }{\breve{y}}_{\iota }\left(\mathfrak{t}\right)=\hfill & -{\stackrel{`}{v}}_{\iota }\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right){\breve{y}}_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\stackrel{`}{\mathfrak{c}}}_{\kappa \iota }\left({\breve{y}}_{\iota }\left(\mathfrak{t}\right)\right){\mathfrak{g}}_{\kappa }\left({\breve{x}}_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\stackrel{`}{\mathfrak{d}}}_{\kappa \iota }\left({\breve{y}}_{\iota }(\mathfrak{t}-{\tau }_2)\right){\mathfrak{g}}_{\kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left({\breve{x}}_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota }+{\eta }_{\iota }.\hfill \end{array}\right.$$

(7)

Assumption A1.

Let function ${\mathfrak{f}}_{\iota }$ satisfy the Lipschitz condition, there exist positive constants ${\mathfrak{f}}_{\iota }^{*}$ such that

$$\mid {\mathfrak{f}}_{\iota }\left({\mathfrak{t}}_1\right)-{\mathfrak{f}}_{\iota }\left({\mathfrak{t}}_2\right)\mid \le {\mathfrak{f}}_{\iota }^{*}\mid {\mathfrak{t}}_1-{\mathfrak{t}}_2\mid$$

where ${\mathfrak{t}}_1,{\mathfrak{t}}_2\in \mathbb{R},{\mathfrak{f}}_{\iota }^{*}$ are positive Lipschitz constants and assume functions ${\mathfrak{g}}_{\kappa },s_{\kappa }$ and $r_{\iota }$ satisfy this condition equally.

Assumption A2.

${\mathfrak{f}}_{\iota }(\pm \mathring{T}={\mathfrak{g}}_{\kappa }(\pm \mathring{T})=0.$

Lemma 3.

[37] Under Assumption A1 and Assumption A2, we know for any ${\stackrel{`}{\mathfrak{c}}}_{\kappa \iota },{\acute{\mathfrak{c}}}_{\kappa \iota }\in co\left[{\mathfrak{c}}_{\kappa \iota }\left(\varrho \right)\right]$

$$\mid {\stackrel{`}{\mathfrak{c}}}_{\kappa \iota }\left({\mathfrak{t}}_2\right){\mathfrak{g}}_{\kappa }\left({\mathfrak{t}}_2\right)-{\acute{\mathfrak{c}}}_{\kappa \iota }\left({\mathfrak{t}}_1\right){\mathfrak{g}}_{\kappa }\left({\mathfrak{t}}_1\right)\mid \le {\mathfrak{c}}_{\kappa \iota }^{*}{\mathfrak{g}}_{\kappa }^{*}\mid {\mathfrak{t}}_2-{\mathfrak{t}}_1\mid ,$$

where ${\mathfrak{c}}_{\kappa \iota }^{*}=max\left\{\mid {\ddot{\mathfrak{c}}}_{\kappa \iota }\mid ,\mid {\dot{\mathfrak{c}}}_{\kappa \iota }\mid \right\}.$

The synchronization error system, by Assumption A1 and Lemma 3, can be expressed as:

$$\left\{\begin{array}{cc}{D}^{\alpha }{e}_{x,\kappa }\left(\mathfrak{t}\right)\le \hfill & -(u_{\kappa }^{*}+{\zeta }_{\kappa })\mid e_{x,\kappa }\left(\mathfrak{t}\right)\mid +\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{a}}_{\iota \kappa }^{*}{f}_{\iota }^{*}\mid e_{y,\iota }\left(\mathfrak{t}\right)\mid +\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{b}}_{\iota \kappa }^{*}{f}_{\iota }^{*}\mid {\tilde{e}}_{y,\iota }\left(\mathfrak{t}\right)\mid \hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_{\iota }^{*}\mid {\tilde{e}}_{y,\iota }\left(\mathfrak{t}\right)\mid dB\left(\mathfrak{t}\right),\hfill \\ D^{\alpha }{e}_{y,\iota }\left(\mathfrak{t}\right)\le \hfill & -(v_{\iota }^{*}+{\zeta }_{\iota })\mid e_{y,\iota }\left(\mathfrak{t}\right)\mid +\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{c}}_{\kappa \iota }^{*}{g}_{\kappa }^{*}\mid e_{x,\kappa }\left(\mathfrak{t}\right)\mid +\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{d}}_{\kappa \iota }^{*}{g}_{\kappa }^{*}\mid {\tilde{e}}_{x,\kappa }\left(\mathfrak{t}\right)\mid \hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }^{*}\mid {\tilde{e}}_{x,\kappa }\left(\mathfrak{t}\right)\mid dB\left(\mathfrak{t}\right),\hfill \end{array}\right.$$

(8)

where denote ${\tilde{e}}_{y,\iota }\left(\mathfrak{t}\right)=e_{y,\iota }(\mathfrak{t}-{\tau }_2),$${\tilde{e}}_{x,\kappa }\left(\mathfrak{t}\right)=e_{x,\kappa }(\mathfrak{t}-{\tau }_1),$$u_{\kappa }^{*}=max\left\{\mid {\ddot{u}}_{\kappa }\mid ,\mid {\dot{u}}_{\kappa }\mid \right\},$$v_{\iota }^{*}=max\left\{\mid {\ddot{v}}_{\iota }\mid ,\mid {\dot{v}}_{\iota }\mid \right\},$ ${\mathfrak{a}}_{\iota \kappa }^{*}=max\left\{\mid {\ddot{\mathfrak{a}}}_{\iota \kappa }\mid ,\mid {\dot{\mathfrak{a}}}_{\iota \kappa }\mid \right\},$${\mathfrak{b}}_{\iota \kappa }^{*}=max\left\{\mid {\ddot{\mathfrak{b}}}_{\iota \kappa }\mid ,\mid {\dot{\mathfrak{b}}}_{\iota \kappa }\mid \right\},$${\mathfrak{c}}_{\kappa \iota }^{*}=max\left\{\mid {\ddot{\mathfrak{c}}}_{\kappa \iota }\mid ,\mid {\dot{\mathfrak{c}}}_{\kappa \iota }\mid \right\}$ and ${\mathfrak{d}}_{\kappa \iota }^{*}=max\left\{\mid {\ddot{\mathfrak{d}}}_{\kappa \iota }\mid ,\mid {\dot{\mathfrak{d}}}_{\kappa \iota }\mid \right\}.$

For the sake of convenience, we can express inequality (8) as

$$\left\{\begin{array}{cc}{D}^{\alpha }{e}_x\left(\mathfrak{t}\right)\le \hfill & -\mathfrak{U}\mid e_x\left(\mathfrak{t}\right)\mid +\mathfrak{A}F\mid e_y\left(\mathfrak{t}\right)\mid +\mathfrak{B}F\mid {\tilde{e}}_y\left(\mathfrak{t}\right)\mid +R\mid {\tilde{e}}_y\left(\mathfrak{t}\right)\mid dB\left(\mathfrak{t}\right),\hfill \end{array}{D}^{\alpha }{e}_y\left(\mathfrak{t}\right)\le -\mathfrak{V}\mid e_y\left(\mathfrak{t}\right)\mid +\mathfrak{C}G\mid e_x\left(\mathfrak{t}\right)\mid +\mathfrak{D}G\mid {\tilde{e}}_x\left(\mathfrak{t}\right)\mid +S\mid {\tilde{e}}_x\left(\mathfrak{t}\right)\mid dB\left(\mathfrak{t}\right),\right.$$

(9)

where $\mathfrak{U}=diag\{u_1^{*}+{\zeta }_1,u_2^{*}+{\zeta }_2,\cdots ,u_{\mathfrak{n}}^{*}+{\zeta }_{\mathfrak{n}}\},$$\mathfrak{V}=diag\{v_1^{*}+{\zeta }_1,v_2^{*}+{\zeta }_2,\cdots ,v_{\mathfrak{m}}^{*}+{\zeta }_{\mathfrak{m}}\},$$\mathfrak{A}={\left({\mathfrak{a}}_{\iota \kappa }^{*}\right)}_{\mathfrak{m}\times \mathfrak{n}},$$\mathfrak{B}={\left({\mathfrak{b}}_{\iota \kappa }^{*}\right)}_{\mathfrak{m}\times \mathfrak{n}},\mathfrak{C}={\left({\mathfrak{c}}_{\kappa \iota }^{*}\right)}_{\mathfrak{n}\times \mathfrak{m}},\mathfrak{D}={\left({\mathfrak{d}}_{\kappa \iota }^{*}\right)}_{\mathfrak{n}\times \mathfrak{m}},F=\underset{1\le l\le \mathfrak{m}}{max}\left\{{\mathfrak{f}}_{\iota }^{*}\right\},G=\underset{1\le \kappa \le \mathfrak{n}}{max}\left\{{\mathfrak{g}}_{\kappa }^{*}\right\},R=\underset{1\le \iota \le \mathfrak{m}}{max}\left\{r_{\iota }^{*}\right\},$ and $S=\underset{1\le \kappa \le \mathfrak{n}}{max}\left\{s_{\kappa }^{*}\right\}.$

Remark 2.

To ensure that $\Vert e_x\left(\mathfrak{t}\right)\Vert +\Vert e_y\left(\mathfrak{t}\right)\Vert \le \chi \left(\mathfrak{t}\right)\le ϵ$, we can simply find an evaluation function ${\chi }_1\left(\mathfrak{t}\right)$ that satisfies $max\left\{\Vert e_x\left(\mathfrak{t}\right)\Vert ,\Vert e_y\left(\mathfrak{t}\right)\Vert \right\}\le {\chi }_1\left(\mathfrak{t}\right)$ with ${\chi }_1\left(\mathfrak{t}\right)\le \frac{ϵ}{2}.$ This will guarantee that $\Vert e_x\left(\mathfrak{t}\right)\Vert +\Vert e_y\left(\mathfrak{t}\right)\Vert$ remains below ϵ, while it also will ensure that $\chi \left(\mathfrak{t}\right)$ is never less than $max\left\{\Vert e_x\left(\mathfrak{t}\right)\Vert ,\Vert e_y\left(\mathfrak{t}\right)\Vert \right\}.$

Lemma 4.

[38] [Burkholder-Davis-Gundy inequality] For any $\Phi \in L_{F_{\mathfrak{t}}}^p([-\tau ,0];H),$ then

$$\begin{array}{c}\hfill \mathbb{E}\left[\underset{{\mathfrak{t}}_0\le \mathfrak{t}\le T}{sup}\mid {\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}\Phi \left(u\right)d\omega \left(u\right){\mid }^p\right]\le c_p{\left[{\int }_{{\mathfrak{t}}_0}^T\mathbb{E}\mid \Phi \left(u\right){\mid }^{2}du\right]}^{\frac{p}{2}}\end{array}$$

where $c_p={\left(\frac{p^{p+1}}{2{(p-1)}^{p-1}}\right)}^{\frac{p}{2}}(p\ge 2)$.

Theorem 2.

Assume Assumption A1, Assumption A2 and Remark 2 hold and the following conditions are satisfied.

(1) If $\varphi >H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill & {\varphi }^2\left\{2+{\left(exp(2^{p-1}{\varphi }^{2p}[4{\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2+{2\Vert \mathfrak{B}\Vert }^2{F}^2\Vert +R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\}\le {\displaystyle \frac{{ϵ}^2}{4}}.\hfill \end{array}\end{array}$$

(2) If $\varphi <H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill & H^2\left(\mathfrak{t}\right)\left\{2+{\left(exp(2^{p-1}{H}^{2p}\left(\mathfrak{t}\right)[4{\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2+{2\Vert \mathfrak{B}\Vert }^2{F}^2\Vert +R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\}\le {\displaystyle \frac{{ϵ}^2}{4}},\hfill \end{array}\end{array}$$

$$H^2\left(\mathfrak{t}\right)={\displaystyle \frac{4{\mathfrak{t}}^{2\alpha -1-\frac{1}{p}}}{{\Gamma }^2\left(\alpha \right){(2q(\alpha -1)+1)}^{\frac{1}{q}}}},$$

where $\mathfrak{t}\in \left[0,T\right],0<ϵ,\frac{1+p}{2p}<\alpha <1,$ and $\frac{1}{p}+\frac{1}{q}=1$ ($p,q\in N$).

Then the drive system (1) and the corresponding system (2) are finite-time synchronized.

Proof 

(Proof).

By Definition 1, for $0<\alpha <1$, we can obtain the following integral inequalities:

$$\left\{\begin{array}{cc}{e}_x\left(\mathfrak{t}\right)\le \hfill & e_x\left(0\right)+\frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}R\mid {\tilde{e}}_y\left(\theta \right)\mid dB\left(\theta \right)+\frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{U}\mid e_x\left(\theta \right)\mid \right.\hfill \\ \hfill & \left.+\mathfrak{A}F\mid e_y\left(\theta \right)\mid +\mathfrak{B}F\mid {\tilde{e}}_y\left(\theta \right)\mid \right)d\theta ,\hfill \\ e_y\left(\mathfrak{t}\right)\le \hfill & e_y\left(0\right)+\frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}S\mid {\tilde{e}}_x\left(\theta \right)\mid dB\left(\theta \right)+\frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{V}\mid e_y\left(\theta \right)\mid \right.\hfill \\ \hfill & \left.+\mathfrak{C}G\mid e_x\left(\theta \right)\mid +\mathfrak{D}G\mid {\tilde{e}}_x\left(\theta \right)\mid \right)d\theta .\hfill \end{array}\right.$$

(10)

By taking the norm on both sides of inequality (10) simultaneously, we obtain

$$\begin{array}{cc}\hfill \Vert e_x\left(\mathfrak{t}\right)\Vert \le & \Vert e_x\left(0\right)\Vert +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}R\mid {\tilde{e}}_y\left(\theta \right)\mid dB\left(\theta \right)\Vert \hfill \\ \hfill & +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{U}\mid e_x\left(\theta \right)\mid +\mathfrak{A}F\mid e_y\left(\theta \right)\mid +\mathfrak{B}F\mid {\tilde{e}}_y\left(\theta \right)\mid \right)d\theta \Vert ,\hfill \end{array}$$

$$\begin{array}{cc}\hfill \Vert e_y\left(\mathfrak{t}\right)\Vert \le & \Vert e_y\left(0\right)\Vert +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}S\mid {\tilde{e}}_x\left(\theta \right)\mid dB\left(\theta \right)\Vert \hfill \\ \hfill & +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{V}\mid e_y\left(\theta \right)\mid +\mathfrak{C}G\mid e_x\left(\theta \right)\mid +\mathfrak{D}G\mid {\tilde{e}}_x\left(\theta \right)\mid \right)d\theta \Vert .\hfill \end{array}$$

Assuming $\Vert e_x\left(\mathfrak{t}\right)\Vert >\Vert e_y\left(\mathfrak{t}\right)\Vert ,$ we can square both sides to obtain $\Vert e_x{\left(\mathfrak{t}\right)\Vert }^2$. By using Lemma 4, we can then rewrite the expression as follows

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Vert e_x{\left(\mathfrak{t}\right)\Vert }^2\le & (\Vert e_x\left(0\right)\Vert +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}R\times \mid {\tilde{e}}_y\left(\theta \right)\mid dB\left(\theta \right)\Vert \hfill \\ \hfill & +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{U}\mid e_x\left(\theta \right)\mid +\mathfrak{A}F\mid e_y\left(\theta \right)\mid +\mathfrak{B}F\mid {\tilde{e}}_y\left(\theta \right)\mid \right)d\theta {\Vert )}^2\hfill \\ \hfill \le & 2\Vert e_x{\left(0\right)\Vert }^2+2(\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}R\mid {\tilde{e}}_y\left(\theta \right)\mid dB\left(\theta \right)\Vert \hfill \\ \hfill & +\Vert \frac{1}{\Gamma \left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{\alpha -1}\left(-\mathfrak{U}\mid e_x\left(\theta \right)\mid +\mathfrak{A}F\mid e_y\left(\theta \right)\mid +\mathfrak{B}F\mid {\tilde{e}}_y\left(\theta \right)\mid \right)d\theta {\Vert )}^2\hfill \\ \hfill \le & 2\Vert e_x{\left(0\right)\Vert }^2+\frac{16}{{\Gamma }^2\left(\alpha \right)}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2(\alpha -1)}{R}^2{\Vert {\tilde{e}}_x\left(\theta \right)\Vert }^{2}d\theta \hfill \\ \hfill & +{\displaystyle \frac{8}{{\Gamma }^2\left(\alpha \right)}}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2(\alpha -1)}[{\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2\Vert e_x{\left(\theta \right)\Vert }^2+{\Vert \mathfrak{B}\Vert }^2{F}^2\Vert {\tilde{e}}_x\left(\theta \right){\Vert }^2]d\theta \hfill \end{array}\end{array}$$

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \le & 2\Vert e_x{\left(0\right)\Vert }^2+{\displaystyle \frac{8}{{\Gamma }^2\left(\alpha \right)}}{\int }_{{\mathfrak{t}}_0}^{\mathfrak{t}}{(\mathfrak{t}-\theta )}^{2(\alpha -1)}[{\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2{\Vert e_x\left(\theta \right)\Vert }^2\hfill \\ \hfill & +{(\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2)\Vert {\tilde{e}}_x\left(\theta \right){\Vert }^2]d\theta .\hfill \end{array}\end{array}$$

From Lemma 2, we define the following functions: $\omega \left(\mathfrak{t}\right)=\Vert e_x\left(\mathfrak{t}\right)\Vert$, $\sigma \left(\mathfrak{t}\right)=\Vert e_x\left(0\right)\Vert =\phi$, $\upsilon \left(\mathfrak{t}\right)=8$, $\mathfrak{a}\left(\mathfrak{t}\right)={\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2$ and $\mathfrak{b}\left(\mathfrak{t}\right)={\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2$. The initial value ${\mathfrak{t}}_0=0$ and $\Omega =\tau$. It is easy to see that all of these functions are non-negative and continuous. Additionally, $\sigma \left(\mathfrak{t}\right)$ and $\upsilon \left(\mathfrak{t}\right)$ are non-decreasing on $[{\mathfrak{t}}_0,T]$, $\varphi \left(\mathfrak{t}\right)$ is non-decreasing on $[{\mathfrak{t}}_0-\Omega ,{\mathfrak{t}}_0]$ and $\sigma \left({\mathfrak{t}}_0\right)=\varphi \left({\mathfrak{t}}_0\right)$.

By using Lemma 2, we can obtain the following results.

(1) If $\phi >H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Vert e_x{\left(\mathfrak{t}\right)\Vert }^2\le & {\phi }^2\left\{2+{\left(exp(2^{p-1}{\phi }^{2p}\left[2\right(\Vert \mathfrak{U}\Vert +{\Vert \mathfrak{A}\Vert F)}^2+{\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\}.\hfill \end{array}\end{array}$$

(2) If $\phi \le H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \Vert e_x{\left(\mathfrak{t}\right)\Vert }^2\le & H^2\left(\mathfrak{t}\right)\left\{2+{\left(exp(2^{p-1}{H}^{2p}\left(\mathfrak{t}\right)\left[2\right(\Vert \mathfrak{U}\Vert +{\Vert \mathfrak{A}\Vert F)}^2+{\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\},\hfill \end{array}\end{array}$$

where $H^2\left(\mathfrak{t}\right)={\displaystyle \frac{4t^{2\alpha -1-\frac{1}{p}}}{{\Gamma }^2\left(\alpha \right){(2q(\alpha -1)+1)}^{\frac{1}{q}}}}.$

Therefore, based on the hypothesis conditions, it can be concluded that systems (1) and (2) can achieve synchronization. □

So when $\Vert e_x\left(\mathfrak{t}\right)\Vert =max\left\{\Vert e_x\left(\mathfrak{t}\right)\Vert ,\Vert e_y\left(\mathfrak{t}\right)\Vert \right\}$, Remark 2 indicates that the evaluation function ${\chi }_1\left(\mathfrak{t}\right)$ can be determined as follows.

(1) If $\phi >H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\chi }_1^2\left(\mathfrak{t}\right)=& {\phi }^2\left\{2+{\left(exp(2^{p-1}{\phi }^{2p}[2{\left(\Vert \mathfrak{U}\Vert +\Vert \mathfrak{A}\Vert F\right)}^2+{\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\}.\hfill \end{array}\end{array}$$

(2) If $\phi \le H\left(\mathfrak{t}\right),$ then

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill {\chi }_1^2\left(\mathfrak{t}\right)=& H^2\left(\mathfrak{t}\right)\left\{2+{\left(exp(2^{p-1}{H}^{2p}\left(\mathfrak{t}\right)\left[2\right(\Vert \mathfrak{U}\Vert +{\Vert \mathfrak{A}\Vert F)}^2+{\Vert \mathfrak{B}\Vert }^2{F}^2+2R^2{]}^p\mathfrak{t})-1\right)}^{\frac{1}{p}}\right\}.\hfill \end{array}\end{array}$$

In Section 4, we will present a numerical example to provide a more visual demonstration of the finite-time synchronization achieved between systems (1) and (2).

4. Numerical examples

Compared to conventional neural networks, neural networks incorporating stochasticity possess greater adaptability and robustness in achieving finite-time synchronization. This is because stochasticity can increase the complexity of the system, endowing it with enhanced fault-tolerance and adaptability, thus facilitating more efficient adaptation to diverse environments and application scenarios. Moreover, neural networks with stochasticity exhibit advantageous characteristics in handling nonlinear and complex problems. In practical applications, the parameters and states of the neural network systems are often uncertain, owing to the presence of uncertainty. Stochasticity can more effectively model such uncertainty and bolster the reliability of the neural network systems, thereby elevating its performance in practical applications. Therefore, investigating neural networks with stochasticity may contribute to enhancing the application capabilities and performance of neural networks.

We illustrate the practical application of Theorem 2 in achieving finite-time synchronization between the systems (1) and (2) through a numerical example. By showing this example, we can validate the effectiveness of the proposed synchronization method. Specifically, the example involves simulating the behavior of the systems with varying initial conditions, and analyzing the resulting trajectories. The insights gained from this example are used to serve as evidence to show the practical relevance of the finite-time synchronization approach presented in this paper.

Example

Consider the fractional-order stochastic MBAMNNs with time delay

$$\left\{\begin{array}{cc}{D}^{\alpha }{x}_{\kappa }\left(\mathfrak{t}\right)=\hfill & -u_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)x_{\kappa }\left(\mathfrak{t}\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{a}}_{\iota \kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)+\sum _{\iota =1}^{\mathfrak{m}}{\mathfrak{b}}_{\iota \kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right){\mathfrak{f}}_{\iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right)\hfill \\ \hfill & +\sum _{\iota =1}^{\mathfrak{m}}{r}_j(y_{\iota }(\mathfrak{t}-{\tau }_2)dB\left(\mathfrak{t}\right)+I_{\kappa },\hfill \\ D^{\alpha }{y}_{\iota }\left(\mathfrak{t}\right)=\hfill & -v_{\iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right)y_{\iota }\left(\mathfrak{t}\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{c}}_{\kappa \iota }\left(y_{\iota }\left(\mathfrak{t}\right)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }\left(\mathfrak{t}\right)\right)+\sum _{\kappa =1}^{\mathfrak{n}}{\mathfrak{d}}_{\kappa \iota }\left(y_{\iota }(\mathfrak{t}-{\tau }_2)\right){\mathfrak{g}}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)\hfill \\ \hfill & +\sum _{\kappa =1}^{\mathfrak{n}}{s}_{\kappa }\left(x_{\kappa }(\mathfrak{t}-{\tau }_1)\right)dB\left(\mathfrak{t}\right)+J_{\iota }.\hfill \end{array}\right.$$

where

$$\begin{array}{cccccc}\hfill u_1\left(x_1\right)=& \left\{\begin{array}{cccc}& 0.6,\hfill & & \mid x_1\mid \le \frac{7}{10},\hfill \\ & -0.6,\hfill & & \mid x_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill u_2\left(x_2\right)=& \left\{\begin{array}{cccc}& 0.4,\hfill & & \mid x_2\mid \le \frac{7}{10},\hfill \\ & -0.4,\hfill & & \mid x_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill v_1\left(y_1\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid y_1\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid y_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill v_2\left(y_2\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid y_2\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid y_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{a}}_{11}\left(x_1\right)=& \left\{\begin{array}{cccc}& 0.1,\hfill & & \mid x_1\mid \le \frac{7}{10},\hfill \\ & -0.1,\hfill & & \mid x_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{a}}_{12}\left(x_2\right)=& \left\{\begin{array}{cccc}& 0.1,\hfill & & \mid x_2\mid \le \frac{7}{10},\hfill \\ & -0.1,\hfill & & \mid x_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{a}}_{21}\left(x_1\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid x_1\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid x_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{a}}_{22}\left(x_2\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid x_2\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid x_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{b}}_{11}\left({\tilde{x}}_1\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid x_1\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid x_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{b}}_{12}\left({\tilde{x}}_2\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid x_2\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid x_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{b}}_{21}\left({\tilde{x}}_1\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid x_1\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid x_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{b}}_{22}\left({\tilde{x}}_2\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid x_2\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid x_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{c}}_{11}\left(y_1\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid y_1\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid y_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{c}}_{12}\left(y_2\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid y_2\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid y_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{c}}_{21}\left(y_1\right)=& \left\{\begin{array}{cccc}& 0.1,\hfill & & \mid y_1\mid \le \frac{7}{10},\hfill \\ & -0.1,\hfill & & \mid y_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{c}}_{22}\left(y_2\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid y_2\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid y_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{d}}_{11}\left({\tilde{y}}_1\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid y_1\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid y_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{d}}_{12}\left({\tilde{y}}_2\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid y_2\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid y_2\mid >\frac{7}{10},\hfill \end{array}\right.\hfill \\ \hfill {\mathfrak{d}}_{21}\left({\tilde{y}}_1\right)=& \left\{\begin{array}{cccc}& 0.2,\hfill & & \mid y_1\mid \le \frac{7}{10},\hfill \\ & -0.2,\hfill & & \mid y_1\mid >\frac{7}{10},\hfill \end{array}\right.\hfill & & & \hfill {\mathfrak{d}}_{22}\left({\tilde{y}}_2\right)=& \left\{\begin{array}{cccc}& 0.3,\hfill & & \mid y_2\mid \le \frac{7}{10},\hfill \\ & -0.3,\hfill & & \mid y_2\mid >\frac{7}{10}.\hfill \end{array}\right.\hfill \end{array}$$

Let $\mathfrak{n}=\mathfrak{m}=2,{\tau }_1={\tau }_2=0.2,{\mathfrak{g}}_{\kappa }\left(x_{\kappa }\right)=tanh(\mid x_{\kappa }\mid -\frac{7}{10}),{\mathfrak{f}}_{\iota }\left(y_{\iota }\right)=tanh(\mid y_{\iota }\mid -\frac{7}{10}),r_{\iota }\left(y_{\iota }\right)=tanh(\mid r_{\iota }\mid -\frac{7}{10}),s_{\kappa }\left(x_{\kappa }\right)=tanh(\mid s_{\kappa }\mid -\frac{7}{10}),I_1=I_2=J_1=J_2=2,x\left(0\right)={(0.1,0.2)}^T,y\left(0\right)={(0.1,0.2)}^T.$

Then, assume ${\zeta }_{\kappa }={(-0.1,-0.2)}^T,{\zeta }_{\iota }={(-0.1,-0.2)}^T,\breve{x}\left(0\right)$$={(0.2,0.3)}^T$ and $\breve{y}\left(0\right)={(0.2,0.3)}^T$ in corresponding system (2).

Let $p=2,q=2,\alpha =0.8>{\displaystyle \frac{1+p}{2p}},ϵ=2$ such that ${\displaystyle \frac{ϵ}{2}}=1.$ It is easy to calculate that $\phi =0.2,\breve{\phi }=0.2,\Vert \mathfrak{U}\Vert =0.5,\Vert \mathfrak{V}\Vert =0.1,F=G=S=R=1,\Vert \mathfrak{A}\Vert =0.3,\Vert \mathfrak{B}\Vert =0.5,\Vert \mathfrak{C}\Vert =0.4$ and $\Vert \mathfrak{D}\Vert =0.6.$ By analyzing Figure 1 and Figure 2, we can obtain the error components $e_{x,1}\left(\mathfrak{t}\right),e_{y,1}\left(\mathfrak{t}\right)$ and $e_{x,2}\left(\mathfrak{t}\right),e_{y,2}\left(\mathfrak{t}\right)$ of systems (1) and (2) when $\kappa ,\iota =1$ and $\kappa ,\iota =2$. By using Figure 3 and Figure 4, we can calculate the error and square of error of systems (1) and (2). Finally, based on Figure 5 and Theorem 2, we can determine the finite-time synchronization time $T=0.3177$.

5. Conclusions

We have enhanced the fractional-order Gronwall inequality for studying finite-time synchronization in fractional-order stochastic MBAMNNs systems, based on the work originally proposed in [34]. Then we presented an illustrative example to show the effectiveness of our proposed approach. However, it should be noted that we have only analyzed the finite-time synchronization of continuous neural network systems, and have not provided a detailed description of discontinuous neural networks with impulses. Hence, we will investigate the dynamic behaviors of fractional-order neural networks with impulses in the future.