Viability approach to Hamilton-Jacobi-Moskowitz problem involving variable regulation parameters

Anya Désilles; Anya Désilles

doi:10.3934/nhm.2013.8.707

Networks and Heterogeneous Media

2013, Volume 8, Issue 3: 707-726. doi: 10.3934/nhm.2013.8.707

Previous Article Next Article

Viability approach to Hamilton-Jacobi-Moskowitz problem involving variable regulation parameters

Anya Désilles ^{1
,}

1.
ENSTA ParisTech, 828, Boulevard des Maréchaux, 91762 Palaiseau Cedex

Received: 01 April 2012 Revised: 01 July 2013
Primary: 49J40, 34A60; Secondary: 35L65.

A few applications of the viability theory to the solution to the Hamilton-Jacobi-Moskowitz problems are presented. In the considered problem the Hamiltonian (fundamental diagram) depends on time, position and/or some regulation parameters. We study such a problem in its equivalent variational formulation. In this case, the corresponding lagrangian depends on the state of the characteristic dynamical system. As the Lax-Hopf formulae that give the solution in a semi-explicit form for an homogeneous lagrangian do not hold, a capture basin algorithm is proposed to compute the Moskowitz function as a viability solution of the Hamilton-Jacobi-Moskowitz problem with general conditions (including initial, boundary and internal conditions). We present two examples of applications to traffic regulation problems.

Keywords:

Citation: Anya Désilles. Viability approach to Hamilton-Jacobi-Moskowitz problem involving variable regulation parameters[J]. Networks and Heterogeneous Media, 2013, 8(3): 707-726. doi: 10.3934/nhm.2013.8.707

Related Papers:

[1]	Yuwei Cao, Bing Li . Existence and global exponential stability of compact almost automorphic solutions for Clifford-valued high-order Hopfield neutral neural networks with $ D $ operator. AIMS Mathematics, 2022, 7(4): 6182-6203. doi: 10.3934/math.2022344
[2]	Nina Huo, Bing Li, Yongkun Li . Global exponential stability and existence of almost periodic solutions in distribution for Clifford-valued stochastic high-order Hopfield neural networks with time-varying delays. AIMS Mathematics, 2022, 7(3): 3653-3679. doi: 10.3934/math.2022202
[3]	Xiaofang Meng, Yongkun Li . Pseudo almost periodic solutions for quaternion-valued high-order Hopfield neural networks with time-varying delays and leakage delays on time scales. AIMS Mathematics, 2021, 6(9): 10070-10091. doi: 10.3934/math.2021585
[4]	Jin Gao, Lihua Dai . Anti-periodic synchronization of quaternion-valued high-order Hopfield neural networks with delays. AIMS Mathematics, 2022, 7(8): 14051-14075. doi: 10.3934/math.2022775
[5]	Abdulaziz M. Alanazi, R. Sriraman, R. Gurusamy, S. Athithan, P. Vignesh, Zaid Bassfar, Adel R. Alharbi, Amer Aljaedi . System decomposition method-based global stability criteria for T-S fuzzy Clifford-valued delayed neural networks with impulses and leakage term. AIMS Mathematics, 2023, 8(7): 15166-15188. doi: 10.3934/math.2023774
[6]	Qian Cao, Xiaojin Guo . Anti-periodic dynamics on high-order inertial Hopfield neural networks involving time-varying delays. AIMS Mathematics, 2020, 5(6): 5402-5421. doi: 10.3934/math.2020347
[7]	Xiaojin Guo, Chuangxia Huang, Jinde Cao . Nonnegative periodicity on high-order proportional delayed cellular neural networks involving $ D $ operator. AIMS Mathematics, 2021, 6(3): 2228-2243. doi: 10.3934/math.2021135
[8]	Yongkun Li, Xiaoli Huang, Xiaohui Wang . Weyl almost periodic solutions for quaternion-valued shunting inhibitory cellular neural networks with time-varying delays. AIMS Mathematics, 2022, 7(4): 4861-4886. doi: 10.3934/math.2022271
[9]	Qi Shao, Yongkun Li . Almost periodic solutions for Clifford-valued stochastic shunting inhibitory cellular neural networks with mixed delays. AIMS Mathematics, 2024, 9(5): 13439-13461. doi: 10.3934/math.2024655
[10]	Hedi Yang . Weighted pseudo almost periodicity on neutral type CNNs involving multi-proportional delays and D operator. AIMS Mathematics, 2021, 6(2): 1865-1879. doi: 10.3934/math.2021113

Abstract

1. Introduction

As is well known, in the practical application of neural networks (NNs), it is often necessary to use and design NN models with different dynamic characteristics for different application scenarios and purposes. Therefore, the study of the dynamic behavior of NNs has become an important issue that is widely concerned in both theoretical research and practical applications of NNs. Therefore, the dynamics of various types of NNs, especially numerous classical NNs such as recurrent NNs ^[1,2], bidirectional associative memory NNs ^[3], inertial NNs ^[4,5], Hopfield NNs ^[6], Cohen-Grossberg NNs ^[7], etc., have been widely studied. It should be mentioned here that due to the stronger approximation, faster convergence speed, larger storage capacity, and higher fault tolerance of high-order NNs compared to low-order NNs, the dynamics research of high-order NNs has received widespread attention ^{[8,9,10,11,12]}.

Meanwhile, owing to the fact that algebra-valued NNs, such as complex-valued ^[13,14], quaternion-valued ^{[15,16,17,18]}, Clifford-valued ^{[19,20,21,22,23]}, and octonion-valued NNs ^[24,25,26], are extensions of real-valued NNs and have more advantages than real-valued NNs in many application scenarios, research on the dynamics of algebra-valued NNs has gradually become a new hotspot in the field of NN research in recent years. It is worth mentioning here that the Clifford-valued high-order Hopfield fuzzy NN represents a sophisticated integration of Clifford algebra, high-order synaptic connections, and fuzzy logic, enabling it to achieve advanced applications in multidimensional data processing and complex system modeling. For example, it has applications in the fields of multidimensional signal processing, secure communication and image encryption, optimization and control systems, neuroscience, and cognitive modeling ^{[27,28,29,30]}.

On the one hand, from both theoretical and practical perspectives, NN models with time-varying connection weights and time-varying external inputs are more realistic than those with constant connection weights and constant external inputs. Meanwhile, time delay effects are inevitable. As a result, the rate of change in the state of a neuron depends not only on its current state but also on its historical state, and even more so, on the rate of change in its historical state. It is precisely for these reasons that researchers have proposed various neutral-type NN models with $ D $ operators and conducted extensive research on their dynamics ^{[31,32,33,34]}. In addition, fuzzy logic and NNs complement each other: fuzzy systems provide interpretability and handle uncertainty, while NNs offer powerful learning from data. Their integration bridges the gap between data-driven machine learning and human-like reasoning, making systems more adaptable, transparent, and robust in real-world applications. Indeed, fuzzy NNs have been successfully applied in many fields such as signal processing, pattern recognition, associative memory, and image processing ^{[35,36,37,38,39]}.

On the other hand, as is well known, almost periodic oscillation is an important dynamic of NNs with time-varying connection weights and time-varying external inputs. In the past few decades, the almost periodic oscillations of various NNs have been studied by countless scholars ^{[23,31,33,34]}. We know that besides Bohr's concept of almost periodicity, there are also Stepanov almost periodicity, Weyl almost periodicity, Besicovitch almost periodicity, and so on ^[40]. It should be pointed out here that Besicovitch almost periodicity is the most complex almost periodicity among Bohr almost periodicity, Stepanov almost periodicity, and Weyl almost periodicity, and that Stepanov almost periodicity, Weyl almost periodicity, and Besicovitch almost periodicity are referred to as generalized almost periodicity. Meanwhile, it should be noted that the product of two generalized almost periodic functions in the same sense may not necessarily be a generalized almost periodic function in that sense. Because of this reason, the emergence of high-order terms in high-order NNs poses difficulties for studying the generalized almost periodic oscillations of high-order NNs. As a consequence, the results of generalized almost periodic oscillations for high-order NNs are still very rare. Thereupon, it is necessary to further study the generalized almost periodic oscillation problem of high-order NNs.

Inspired by the above observations, this paper considers a class of Clifford-valued high-order Hopfield fuzzy NNs with time-varying delays and $ D $ operators as follows:

$\begin{align} &[x_i(t)-a_i(t)x_i(t-\tau_i(t))]^{\prime} \\ = &-b_i(t)x_i(t)+\sum\limits_{j = 1}^{n}c_{ij}(t)f_j(x_j(t))+\sum\limits_{j = 1}^{n}u_{ij}(t)f_j(x_j(t-\sigma_{ij}(t)))\\ &+\sum\limits_{j = 1}^{n}\gamma_{ij}(t)\mu_j(t)+\sum\limits_{j = 1}^{v}\sum\limits_{k = 1}^{n}\theta_{ijk}(t)g_j(x_j(t-\delta_{ijk}(t)))g_k(x_k(t-\delta_{ijk}(t)))\\ &+\bigwedge\limits_{j = 1}^{n}\alpha_{ij}(t)f_j(x_j(t-\eta_{ij}(t)))+\bigvee\limits_{j = 1}^{n}\beta_{ij}(t)f_j(x_j(t-\eta_{ij}(t)))\\ &+\bigwedge\limits_{j = 1}^{n}\bigwedge\limits_{k = 1}^{n}q_{ijk}(t)g_j(x_j(t-\delta_{ijk}(t)))g_k(x_k(t-\delta_{ijk}(t)))\\ &+\bigvee\limits_{j = 1}^{n}\bigvee\limits_{k = 1}^{n}\nu_{ijk}(t)g_j(x_j(t-\delta_{ijk}(t)))g_k(x_k(t-\delta_{ijk}(t)))\\ &+\bigwedge\limits_{j = 1}^{n}T_{ij}(t)\mu_{j}(t)+\bigvee\limits_{j = 1}^{n}S_{ij}(t)\mu_{j}(t)+I_{i}(t), \end{align}$ $

(1.1)

where $ i \in\mathcal{J}: = \{1, 2, \cdots, n\} $, $ x_i(t) \in\mathcal{A} $ indicates the state of the $ i $th unit at time $ t $; $ \mathcal{A} $ is a real Clifford algebra; $ b_{i}(t)\in\mathcal{A} $ represents the self feedback coefficient at time $ t $; $ \alpha_{ij}(t), \beta_{ij}(t), T_{ij}(t), S_{ij}(t)\in\mathcal{A} $ stand for the elements of the fuzzy feedback MIN template and fuzzy feed forward MAX template, respectively; $ a_{i}(t), c_{ij}(t), u_{ij}(t) $ and $ \theta_{ijk}(t), q_{ijk}(t), \nu_{ijk}(t) \in\mathcal{A} $ represent the first-order and second-order connection weights of the NN; $ \gamma_{ij}(t) $ stands for the element of the feed forward template; $ \bigwedge $ and $ \bigvee $ denote the fuzzy AND and OR operations, respectively; $ \mu_j(t)\in\mathcal{A} $ represents the input of the $ j $th neuron; $ I_{i}(t)\in\mathcal{A} $ corresponds to the external input to the $ i $th unit; $ f_j $ and $ g_j: \mathcal{A}\to\mathcal{A} $ signify the nonlinear activation functions; and $ \tau_i(t), \sigma_{ij}(t), \eta_{ij}(t), \delta_{ijk}(t)\in\mathbb{R^{+}} $ denote the transmission delays.

The initial value condition associated with (1.1) is given as

$\begin{equation} \begin{split} x_{i}(s) = \varphi_{i}(s), s\in[-\varrho, 0], \, \, i\in\mathcal{J}, \end{split} \end{equation}$ $

(1.2)

where $ \varphi_{i}\in BC([-\varrho, 0], \mathcal{A}), \varrho = \max\limits_{i, j, k\in \mathcal{J}}\big\{\sup\limits_{t\in \mathbb{R}}\tau_i(t), \sup\limits_{t\in \mathbb{R}}\sigma_{ij}(t), \sup\limits_{t\in \mathbb{R}}\eta_{ij}(t), \sup\limits_{t\in \mathbb{R}}\delta_{ijk}(t)\big\} $.

The main purpose of this paper is to investigate the existence and stability of Besicovitch almost periodic solutions for system (1.1). The main contributions of this paper are as follows:

1. This paper is the first one to investigate the existence of Besicovitch almost periodic solutions for system (1.1), and the results of this paper still hold true and are new in the following special cases of system (1.1).

(ⅰ) System (1.1) is a real-valued, complex-valued, or quaternion-valued system.

(ⅱ) System (1.1) is a real-valued, complex-valued, or quaternion-valued system without $ D $ operators, i.e. $ a_i (t) = 0 $.

(ⅲ) System (1.1) is a real-valued system without $ D $ operators and fuzzy terms, i.e., $ a_i (t) = \alpha_{ij}(t) = q_{ijk}(t) = \nu_{ijk}(t) = T_{ij}(t) = S_{ij}(t)\equiv 0 $.

2. The research method proposed in the paper can be used to study the generalized almost periodic dynamics for other high-order NNs.

Remark 1.1. The method we propose can be summarized as follows: First, we use the fixed point theorem to prove the existence of solutions for system (1.1) that are bounded and continuous with respect to the Besicovitch seminorm on a closed subset of an appropriate Banach space. Then, we apply the definition and inequality techniques to prove that this solution is Besicovitch almost periodic.

The remaining part of the paper is arranged as follows: In the second section, we review some relevant concepts, introduce some symbols used in this article, cite a useful lemma, and state and prove the completeness of the space we will use. In the third section, we investigate the existence and stability of Besicovitch almost periodic solutions for system (1.1). In the fourth section, we provide an example to demonstrate the correctness of our results. Finally, in the fifth section, we provide a brief conclusion.

2. Preliminaries

Let $ \mathcal{A} = \Big\{\sum\limits_{A\in\Omega}x^{A}e_{A}\in\mathbb{R}\Big\} $ indicate a real Clifford algebra over $ \mathbb{R}^m $ ^[41], where $ \Omega = \{\emptyset, 1, 2, \cdots, A, \cdots, 12, \cdots, m\} $, $ e_A = e_{h_1}e_{h_2}\cdots e_{h_v} $, $ 1\le h_1 < h_2 < \cdots < h_v\le m $, and in addition, $ e_{\emptyset} = e_0 = 1 $, and $ e_h $, $ h = 1, 2, \cdots, m $ are said to be Clifford generators and satisfy $ e_p = 1, p = 0, 1, 2, \cdots, s, e_p^{2} = -1, p = s+1, s+2, \cdots, m $, where $ s < m $, and $ e_{p}e_{q}+e_{q}e_{p} = 0, p\neq q, p, q = 1, 2, \cdots, m. $ For every $ x = \sum\limits_{A\in\Omega}x^{A}e_{A}\in\mathcal{A} $ and $ y = (y_{1}, y_{2}, \cdots, y_{n})^{T}\in\mathcal{A}^{n} $, we define $ |x|_{1} = \max\limits_{A\in\Omega}\{|x^{A}|\} $ and $ |y|_{n} = \max\limits_{i\in\mathcal{J}}\{|y_{i}|_{1}\} $, respectively, and then the spaces $ (\mathcal{A}, |\cdot|_1) $ and $ (\mathcal{A}^n, |\cdot|) $ are Banach ones.

Since there is no order relation among Clifford numbers, as in ^[42], for $ x = \sum\limits_{A\in\Omega}x^{A}e_{A}, y = \sum\limits_{A\in\Omega}y^{A}e_{A}, $ we define $ x\bigwedge y = \sum\limits_{A\in\Omega}(\min\{x^{A}, y^{A}\})e_{A} $ and $ x\bigvee y = \sum\limits_{A\in\Omega}(\max\{x^{A}, y^{A}\})e_{A}. $ According to this regulation, for example, regarding the 6th and 7th terms on the right-hand side of Equation (1.1), we have

$\begin{align*} \bigwedge\limits_{j = 1}^{n}\alpha_{ij}(t)f_j(x_j(t-\eta_{ij}(t))) = \sum\limits_{A\in \Omega}\Big(\min\limits_{1\leq j\leq n}\big\{\alpha^A_{ij}(t)f_j^A(x_j(t-\eta_{ij}(t)))\big\}\Big)e_A \end{align*}$ $

and

$\begin{align*} \bigvee\limits_{j = 1}^{n}\beta_{ij}(t)f_j(x_j(t-\eta_{ij}(t))) = \sum\limits_{A\in \Omega}\Big(\max\limits_{1\leq j\leq n}\big\{\alpha^A_{ij}(t)f_j^A(x_j(t-\eta_{ij}(t)))\big\}\Big)e_A. \end{align*}$ $

For $ x = \sum\limits_{A\in\Omega}x^{A}e_{A}\in\mathcal{A} $, we indicate $ x^c = x-x^{\emptyset} $.

For the sake of generality in the subsequent discussion of this section, let $ (\mathbb{X}, \|\cdot\|) $ be a Banach space and $ L_{loc}^{p}(\mathbb{R}, \mathbb{X}) $ with $ 1\le p < +\infty $ be the space consisting of measurable and locally $ p $-integrable functions from $ \mathbb{R} $ into $ \mathbb{X} $. In the next section, we will take $ \mathbb{X} = \mathbb{R} $, $ \mathbb{X} = \mathcal{A} $, or $ \mathbb{X} = \mathcal{A}^n $.

Definition 2.1. ^[40] A bounded continuous function $ \varphi: \mathbb{R} \rightarrow \mathbb{X} $ is said to be almost periodic, if for every $ \varepsilon > 0 $, there exists a number $ \ell(\varepsilon) > 0 $ such that for each $ a\in \mathbb{R} $, there exists a point $ \sigma \in [a, a+\ell] $ satisfying

$ \|\varphi(t+\sigma)-\varphi(t)\| < \varepsilon. $

The family of such functions will be signified by $ AP(\mathbb{R}, \mathbb{X}) $.

For $ \varphi\in L_{loc}^{p}(\mathbb{R}, \mathbb{X}) $, the Besicovitch seminorm is defined as the following:

$\begin{align*} \|\varphi\|_{B^{p}} = \Big\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\|\varphi(t)\|^{p}dt\Big\}^{\frac{1}{p}}. \end{align*}$ $

Definition 2.2. ^[43] A function $ \varphi\in L_{loc}^p(\mathbb{R}, \mathbb{X}) $ is called $ B^{p} $-continuous if $ \lim\limits_{h\to0}\|\varphi(\cdot+h)-\varphi(\cdot)\|_{B^{p}} = 0 $ and is called $ B^{p} $-bounded if $ \|\varphi\|_{B^{p}} < \infty. $

Henceforth, we will denote the set of all functions that are $ B^{p} $-continuous and $ B^{p} $-bounded by $ BC B^p(\mathbb{R}, \mathbb{X}). $

Definition 2.3. ^[40] A function $ \varphi\in L_{loc}^{p}(\mathbb{R}, \mathbb{X}) $ is said to be Besicovitch almost periodic, if for every $ \varepsilon > 0 $, there exists a positive number $ \ell > 0 $ such that for each $ a\in \mathbb{R} $, there exists a point $ \sigma \in [a, a+\ell] $ satisfying

$\|\varphi(\cdot+\sigma)-\varphi(\cdot)\|_{B^{p}} < \varepsilon .$

Denote by $ B^{p}_{AP}(\mathbb{R}, \mathbb{X}) $ the class of such functions and, for simplicity, call them $ B^{p} $-almost periodic functions.

Lemma 2.1. ^[44] If $ \alpha_{ij}, \beta_{ij}\in C(\mathbb{R}, \mathcal{A}), g_{j}\in C(\mathcal{A}, \mathcal{A}), i, j\in\mathcal{J}, $ then one has

$ \Big|\bigwedge\limits_{i = 1}^{n}\alpha_{ij}(t)g_{j}(x)-\bigwedge\limits_{i = 1}^{n}\alpha_{ij}(t)g_{j}(y)\Big|_1 \leq\sum\limits_{i = 1}^{n}|\alpha_{ij}(t)\|\|g_{j}(x)-g_{j}(y)|_1, \\ \Big|\bigvee\limits_{i = 1}^{n}\beta_{ij}(t)g_{j}(x)-\bigvee\limits_{i = 1}^{n}\beta_{ij}(t)g_{j}(y)\Big|_1 \leq\sum\limits_{i = 1}^{n}\|\beta_{ij}(t)\||g_{j}(x)-g_{j}(y)|_1. $

Let $ L^{\infty}(\mathbb{R}, \mathbb{X}) $ be the set of all essentially bounded measurable functions from $ \mathbb{R} $ to $ \mathbb{X} $, then $ (L^{\infty}(\mathbb{R}, \mathbb{X}), \|\cdot\|_\infty) $ is a Banach space, where $ \|\cdot\|_\infty: = ess\sup\limits_{t\in\mathbb{R}}\|\cdot\| $ denotes the essential supremum norm.

Denote

$ \mathbb{Z} = \{x|x\in L^{\infty}(\mathbb{R}, \mathbb{X})\cap BCB^{p}(\mathbb{R}, \mathbb{X})\}. $

Then we have the following lemma which is crucial for the proof of our main result of this paper.

Lemma 2.2. The space $ (\mathbb{Z}, \|\cdot\|_{\infty}) $ is a Banach space.

Proof. Let $ \{{\phi_{n}}\} \subset\mathbb{Z} $ be a Cauchy sequence, and then for every $ \varepsilon > 0 $, there is a positive integer $ N_{1} $ such that for $ n, m > N_{1}, $

$ \|\phi_{n}(\cdot)-\phi_{m}(\cdot)\|_{\infty} < \frac{\varepsilon}{3} . $

Since $ \{\phi_n\}\subset\mathbb{Z}\subset L^{\infty}(\mathbb{R}, \mathbb{X}) $ and $ (L^{\infty}(\mathbb{R}, \mathbb{X}), \|\cdot\|_{\infty}) $ is a Banach space, there exists $ \phi\in L^\infty(\mathbb{R}, \mathcal{A}^n) $ such that $ \phi_n\to \phi $ as $ n\to\infty $ with respect to the norm $ \|\cdot\|_\infty. $ To complete the proof, it suffices to prove that $ \phi\in BCB^p(\mathbb{R}, \mathbb{X}). $ From $ \lim\limits_{n\to\infty} \phi_n = \phi $ in regard to the essential supremum norm, it follows that there exists a positive integer $ N_{2} $ such that for $ n > N_{2} $,

$ \|\phi_n(\cdot)-\phi(\cdot)\|_\infty < \frac{\varepsilon}{3}. $

Now, take $ N_0 = \max\{N_1, N_2\} $, and then, due to the fact that $ \phi_{N_0+1}\in BCB^p(\mathbb{R}, \mathbb{X}) $, there exists a $ \delta = \delta(\varepsilon) > 0 $ such that for any $ h\in \mathbb{R} $ with $ |h| < \varepsilon $, it holds that

$ \|\phi_{N_{0}+1}(\cdot+h)-\phi_{N_{0}+1}(\cdot)\|_{B^{p}} < \frac{\varepsilon}{3}. $

Consequently,

$\begin{align*} &\|\phi(\cdot+h)-\phi(\cdot)\|_{B^{p}}\\ &\leq\|\phi(\cdot+h)-\phi_{N_{0}+1}(\cdot+h)\|_{B^{p}}+\|\phi_{N_{0}+1}(\cdot+h)-\phi_{N_{0}+1}(\cdot)\|_{B^{p}}+\|\phi_{N_{0}+1}(\cdot)-\phi(\cdot)\|_{B^{p}}\\ &\leq\|\phi(\cdot+h)-\phi_{N_{0}+1}(\cdot+h)\|_{\infty}+\|\phi_{N_{0}+1}(\cdot+h)-\phi_{N_{0}+1}(\cdot)\|_{B^{p}}+\|\phi_{N_{0}+1}(\cdot)-\phi(\cdot)\|_{\infty}\\ &\leq\frac{\varepsilon}{3}+\frac{\varepsilon}{3}+\frac{\varepsilon}{3} = \varepsilon, \end{align*}$ $

which implies $ \phi\in BCB^p(\mathbb{R}, \mathbb{X}). $ The proof is completed. □

3. Main results

In this section, for $ x\in L^{\infty}(\mathbb{R}, \mathcal{A}) $, we denote $ |x|_{\infty} = \max\limits_{A\in\Omega}\{ess\sup\limits_{t\in\mathbb{R}}|x^{A}(t)|\} $ and for $ z = (z_{1}, z_{2}, \cdots, z_{n})^{T} = (\sum\limits_{A}z_{1}^{A}e_A, \sum\limits_{A}z_{2}^{A}e_A, \cdots, \sum\limits_{A}z_{n}^{A}e_A)^{T}\in L^{\infty}(\mathbb{R}, \mathcal{A}^{n}) $, we denote $ \|z\|_{\infty} = \max\limits_{i\in\mathcal{J}}\{|z_i|_\infty\} $. Let $ \mathbb{Z} = \{z|z\in L^{\infty}(\mathbb{R}, \mathcal{A}^n)\cap BCB^{p}(\mathbb{R}, \mathcal{A}^n)\}, $ and then, according to Lemma 2.2, $ (\mathbb{Z}, \|\cdot\|_{\infty}) $ is a Banach space. For $ x\in B^{p}_{AP}(\mathbb{R}, \mathcal{A}) $ and $ z\in B^{p}_{AP}(\mathbb{R}, \mathcal{A}^n) $, we will use $ |x|_{B^p} $ and $ \|z\|_{B^p} $ to represent the seminorms of $ x $ and $ z $, respectively.

In what follows, we will employ the following symbols:

$ \bar{g} = \sup\limits_{t\in\mathbb{R}}\|g(t)\|_\mathbb{Y} \; \text{and}\; \underline{g} = \inf\limits_{t\in\mathbb{R}}\|g(t)\|_\mathbb{Y}, $

where $ g:\mathbb{R}\to\mathbb{Y} $ is a bounded function and $ (\mathbb{Y}, \|\cdot\|_\mathbb{Y}) $ is a normed space. Moreover, we will use the following assumptions:

$ (A_1) $ For $ i, j, k\in\mathcal{J} $, functions $ b^\emptyset_i\in AP(\mathbb{R}, \mathbb{R}^+) $ with $ \underline{b}_i^\emptyset > 0 $, $ a_{i}, b^c_i, \mu_{j}, c_{ij}, u_{ij}, \alpha_{ij}, \beta_{ij}, \theta_{ijk}, q_{ijk}, \\ \nu_{ijk}\in AP(\mathbb{R}, \mathcal{A}), \tau_{i}, \sigma_{ij}, \eta_{ij}, \delta_{ijk}\in AP(\mathbb{R}, \mathbb{R})\cap C^{1}(\mathbb{R}, \mathbb{R^{+}}) $ with $ \tau^\prime_{i}(t)\le\bar{\tau}_{i}^{\prime} < 1, \sigma^{\prime}_{ij}(t)\le\bar{\sigma}_{ij}^\prime < 1, \eta^{\prime}_{ij}(t)\le\bar{\eta}_{ij}^\prime < 1, \delta^{\prime}_{ijk}(t)\le\bar{\delta}_{ijk}^\prime < 1 $, where $ \bar{\tau}^\prime, \bar{\sigma}^{\prime}, \bar{\eta}^\prime, \bar{\delta}^{\prime} $ are constants, and $ \gamma_{ij}, T_{ij}, S_{ij}, I_i\in L^{\infty}(\mathbb{R}, \mathcal{A}) \cap B_{AP}^{p}(\mathbb{R}, \mathcal{A}). $

$ (A_2) $ For all $ j\in\mathcal{J} $, functions $ f_j, g_j\in C(\mathcal{A}, \mathcal{A}) $ with $ f_j(0) = 0, g_j(0) = 0 $, and there exist positive constants $ L_{j}^{f}, L_{j}^{g}, M_{j}^{g} $, and $ M_{k}^{g} $ such that for any $ u, v\in\mathcal{A} $,

$\begin{align*} |f_j(u)-f_j(v)|_1\leq L_j^f|u-v|_1, |g_j(u)-g_j(v)|_1\leq L_j^g|u-v|_1, |g_{j}(u)|_{1}\leq M_{j}^{g}, |g_{k}(u)|_{1}\leq M_{k}^{g}. \end{align*}$ $

$ (A_3) $ For $ i\in\mathcal{J} $, there exist positive constants $ \vartheta_{i} $ such that

$\begin{align*} \rho: = &\max\limits_{i\in\mathcal{J}}\bigg\{\bar{a}_{i}+\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg[\bar{b}_{i}^{\emptyset}\bar{a}_{i}+\bar{b}_{i}^{c}+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}\bar{c}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{u}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\theta}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\\ &+\sum\limits_{j = 1}^{n}\bar{\alpha}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\beta}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{q}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\nu}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\bigg)\bigg]\bigg\} < 1. \end{align*}$ $

$ (A_4) $ For the constants $ \vartheta_{i}, i\in J $, mentioned in $ (A_3) $, and $ p, q > 1 $ with $ \frac{1}{p}+\frac{1}{q} = 1 $, it holds that

$\begin{align*} P: = \, &2^{p-1}\max _{i \in \mathcal{J}}\bigg\{4^{p-1}(\bar{a}_{i})^{p}\frac{1}{1-\bar{\tau}'_{i}}+70^{p-1}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}} (\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p} \frac{e^{\underline{b}_{i}^{\emptyset}\bar{\tau}_{i}}}{1-\bar{\tau}'_{i}}+70^{p-1} \bigg( \frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{b}_i^c)^{p}\\ &+35^{p-1} \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^{n}(\bar{c}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}+70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+1}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{u}_{ij})^{q}\bigg)^{\frac{p}{q}}\\ &\times \sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\sigma}_{ij}}}{1-\bar{\sigma}'_{ij}} +140^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg[\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\theta}_{i j k}M_k^gL_j^g\vartheta_j)^p\\ &+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\theta}_{ijk}M_j^gL_{k}^{g}\vartheta_{k})^{p}\bigg]\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta}'_{ijk}}+70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n( \bar{\alpha}_{i j})^{q}\bigg)^{\frac{p}{q}}\\ &\times\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p} \frac{e^{\underline{b}_{i}^{\emptyset}\bar{\eta}_{ij}}}{1-\bar{\eta}'_{ij}} +70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n (\bar{\beta}_{i j})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\eta}_{ij}}}{1-\bar{\eta}'_{ij}}\\ &+140^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_k^gL_{j}^{g}\vartheta_{j})^{p}+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_j^gL_k^g\vartheta_k)^p\bigg)\\ &\times \frac{e^{\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta}'_{ijk}}+140^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{k}^{g}L_{j}^{g}\vartheta_{j})^{p}\\ &+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{j}^{g}L_{k}^{g}\vartheta_{k})^{p}\bigg) \frac{e^{\frac{p}{q}\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta}'_{ijk}}\bigg\} < 1. \end{align*}$ $

For $ i\in\mathcal{J} $, let $ y_i(t) = \vartheta_{i}^{-1}x_i(t), Z_i(t) = y_i(t)-a_i(t)y_i(t-\tau_i(t)) $, where $ \vartheta_i > 0 $ are constants, and then system (1.1) turns into

$\begin{align} Z^{\prime}_i(t) = &-b_i^{\emptyset}(t)Z_i(t)-b_i^{\emptyset}(t)a_i(t)y_i(t-\tau_i(t)) -b_i^{c}(t)y_i(t)\\ &+\vartheta_{i}^{-1}\bigg[\sum\limits_{j = 1}^{n}c_{ij}(t)f_j(\vartheta_{j}y_j(t)+\sum\limits_{j = 1}^{n}u_{ij}(t)f_j(\vartheta_{j}y_j(t-\sigma_{ij}(t)))\\ &+\sum\limits_{j = 1}^{n}\gamma_{ij}(t)\mu_{j}(t)+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\theta_{ijk}(t)g_j(\vartheta_{j}y_j(t-\delta_{ijk}(t)))g_k(\vartheta_{k}y_k(t-\gamma_{ijk}(t)))\\ &+\bigwedge\limits_{j = 1}^{n}\alpha_{ij}(t)f_j(\vartheta_{j}y_j(t-\eta_{ij}(t)))+\bigvee\limits_{j = 1}^{n}\beta_{ij}(t)f_j(\vartheta_{j}y_j(t-\eta_{ij}(t)))\\ &+\bigwedge\limits_{j = 1}^{n}\bigwedge\limits_{k = 1}^{n}q_{ijk}(t)g_j(\vartheta_{j}y_j(t-\delta_{ijk}(t)))g_k(\vartheta_{k}y_k(t-\delta_{ijk}(t)))\\ &+\bigvee\limits_{j = 1}^{n}\bigvee\limits_{k = 1}^{n}\nu_{ijk}(t)g_j(\vartheta_{j}y_j(t-\delta_{ijk}(t)))g_k(\vartheta_{k}y_k(t-\delta_{ijk}(t)))\\ &+\bigg.\bigwedge\limits_{j = 1}^{n}T_{ij}(t)\mu_i(t)+\bigvee\limits_{j = 1}^{n}S_{ij}(t)\mu_i(t)+I_i(t)\bigg], \, \, i\in\mathcal{J}. \end{align}$ $

(3.1)

Multiplying both sides of (3.1) with $ e^{\int_{t_0}^{t}b_i^{\emptyset}(u)du} $ and integrating over the interval $ \left[t_0, t\right] $, then it holds that

$\begin{align} y_i(t) = \, & a_i(t)y_i(t-\tau_i(t))+[y_i(t_0)-a_i(t_0)y_i(t_0-\tau_i(t_0))]e^{-\int_{t_0}^{t}b_i^{\emptyset}(u)du}\\ &+\int_{t_0}^{t}e^{-\int_{s}^{t}b_i^{\emptyset}(u)du}(N^{y})_{i}(s)ds, \, \, i\in \mathcal{J}, \end{align}$ $

(3.2)

where

$\begin{align} (N^{y})_{i}(s) = &-b_i^{\emptyset}(s)a_i(s)y_i(s-\tau_i(s)) -b_i^{c}(s)y_i(s)+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}c_{ij}(s)f_j(\vartheta_{j}y_j(s))\bigg.\\ &+\sum\limits_{j = 1}^{n}u_{ij}(s)f_j(\vartheta_{j}y_j(s-\sigma_{ij}(s)))+\sum\limits_{j = 1}^{n}\gamma_{ij}(s)\mu_{j}(s)\bigg.\\ &+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\theta_{ijk}(s)g_j(\vartheta_{j}y_j(s-\delta_{ijk}(s)))g_k(\vartheta_{k}y_k(s-\delta_{ijk}(s)))\\ &+\bigwedge\limits_{j = 1}^{n}\alpha_{ij}(s)f_j(\vartheta_{j}y_j(s-\eta_{ij}(s)))+\bigvee\limits_{j = 1}^{n}\beta_{ij}(s)f_j(\vartheta_{j}y_j(s-\eta_{ij}(s)))\\ &+\bigwedge\limits_{j = 1}^{n}\bigwedge\limits_{k = 1}^{n}q_{ijk}(s)g_j(\vartheta_{j}y_j(s-\delta_{ijk}(s)))g_k(\vartheta_{k}y_k(s-\delta_{ijk}(s)))\\ &+\bigvee\limits_{j = 1}^{n}\bigvee\limits_{k = 1}^{n}\nu_{ijk}(s)g_j(\vartheta_{j}y_j(s-\delta_{ijk}(s)))g_k(\vartheta_{k}y_k(s-\delta_{ijk}(s)))\\ &+\bigg.\bigwedge\limits_{j = 1}^{n}T_{ij}(s)\mu_i(s)+\bigvee\limits_{j = 1}^{n}S_{ij}(s)\mu_i(s)+I_i(s)\bigg). \end{align}$ $

It is easy to verify that if $ y^{\ast}(t) = \left(y_1^{\ast}(t), y_2^{\ast}(t), \cdots, y_n^{\ast}(t)\right) $ solves system (3.1), then $ x^{\ast}(t) = (x_1^{\ast}(t), x_2^{\ast}(t), \cdots, x_n^{\ast}(t)) = (\vartheta_{1}^{-1}y_1^{\ast}(t), \vartheta_{2}^{-1}y_2^{\ast}(t), \cdots, \vartheta_{n}^{-1}y_n^{\ast}(t)) $ solves system (1.1).

Definition 3.1. A function $ x = \left(x_1, x_2, \cdots, x_n\right):\mathbb{R}\rightarrow \mathcal{A}^n $ is called a solution of (1.1) provided that there exist positive numbers $ \vartheta_i $ such that $ y_i(t) = \vartheta_{i}^{-1}x_i(t), i\in \mathcal{J} $ fulfill (3.2).

Set

$ \hat{\varphi} = (\hat{\varphi}_{1}(t), \hat{\varphi}_{2}(t), \cdots, \hat{\varphi}_{n}(t))^{T}, $

where

$\begin{align*} \hat{\varphi}_{i}(t) = \int_{-\infty}^{t}e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\bigg(\sum\limits_{j = 1}^{n}\gamma_{ij}(s)\mu_{j}(s) +\bigwedge\limits_{j = 1}^{n}T_{ij}(s)\mu_{j}(s)+\bigvee\limits_{j = 1}^{n}S_{ij}(s)\mu_{j}(s)+I_{i}(s)\bigg)ds, \, \, i\in\mathcal{J}. \end{align*}$ $

It is easy to see that $ \hat{\varphi} $ is well defined under condition $ (A_1) $. Choose a positive constant $ r $ with $ r > \|\hat{\varphi}\|_{\infty} $.

Then, we are now in a position to present and prove our existence result.

Theorem 3.1. Assume that $ (A_{1})-(A_{4}) $ hold. Then, system (1.1) admits a unique $ B^p $-almost periodic solution in $ \mathbb{Z}^*: = \{\varphi|\varphi\in\mathbb{Z}, \|\varphi-\hat{\varphi}\|_{\infty}\le \frac{\rho r}{1-\rho}\} $.

Proof. Letting $ t_0\rightarrow \infty $, from (3.2), one gets

$\begin{align*} y_i(t) = a_i(t)y_i(t-\tau_i(t)) +\int_{-\infty}^{t}e^{-\int_{s}^{t}b_i^{\emptyset}(u)du}(N^{y})_{i}(s)ds, \, \, i\in \mathcal{J}. \end{align*}$ $

Define a mapping $ T:\mathbb{Z}^*\rightarrow \mathbb{Z}^* $ by setting $ (T\varphi)(t) = ((T\varphi)_{1}(t), (T\varphi)_{2}(t), \cdots, (T\varphi)_{n}(t))^{T} $ for $ \varphi\in\mathbb{Z} $ and $ t\in \mathbb{R} $, where $ (T\varphi)_{i}(t) = a_{i}(t)\varphi_{i}(t-\tau_{i}(t))+\int_{-\infty}^{t}e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}(N^{\varphi})_{i}(s)ds, i\in\mathcal{J}. $

To begin with, we show that $ T(\mathbb{Z}^*)\subset\mathbb{Z}^*. $

Note that, for any $ \varphi\in\mathbb{Z}^* $, it holds that

$ \|\varphi\|_{\infty}\leq\|\varphi-\hat{\varphi}\|_{\infty}+\|\hat{\varphi}\|_{\infty}\le\frac{\rho r}{1-\rho}+r = \frac{r}{1-\rho} $

and $ \|N^\varphi\|_\infty < \infty. $

For every $ \varphi\in\mathbb{Z}^* $, we infer that

$\begin{align*} \|T\varphi-\hat{\varphi}\|_{\infty} = \, & \max\limits_{i\in\mathcal{J}}\bigg\{ess\sup\limits_{t\in\mathbb{R}}\bigg|a_{i}(t)\varphi_{i}(t-\tau_{i}(t))+\int_{-\infty}^{t}e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\bigg[-b_{i}^{\emptyset}(s)a_{i}(s)\varphi_{i}(s-\tau_{i}(s))\\ &-b_{i}^{c}(s)\varphi_{i}(s)+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}c_{ij}(s)f_{j}(\vartheta_{j}\varphi_{j}(s))+\sum\limits_{j = 1}^{n}u_{ij}(s)f_{j}(\vartheta_{j}\varphi_{j}(s-\sigma_{ij}(s)))\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\theta_{ijk}(s)g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\\ &+\bigwedge\limits_{j = 1}^{n}\alpha_{ij}(s)f_{j}(\vartheta_{j}\varphi_{j}(s-\eta_{ij}(s)))+\bigvee\limits_{j = 1}^{n}\beta_{ij}(s)f_{j}(\vartheta_{j}\varphi_{j}(s-\eta_{ij}(s)))\\ &+\bigwedge\limits_{j = 1}^{n}\bigwedge\limits_{k = 1}^{n}q_{ijk}(s)g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\\ &+\bigvee\limits_{j = 1}^{n}\bigwedge\limits_{k = 1}^{n}\nu_{ijk}(s)g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\bigg)\bigg]ds\bigg|_{1}\bigg\}\\ \le&\max\limits_{i\in\mathcal{J}}\bigg\{\bar{a}_{i}\|\varphi\|_{\infty}+\int_{-\infty}^{t}e^{-\underline{b}_{i}^{\emptyset}(t-s)}\bigg[\bar{b}_{i}^{\emptyset}\bar{a}_{i}+\bar{b}_{i}^{c}+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}\bar{c}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{u}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\theta}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\vartheta_{k}+\sum\limits_{j = 1}^{n}\bar{\alpha}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\beta}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\bar{q}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\vartheta_{k}+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\nu}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\vartheta_{k}\bigg)\bigg]\|\varphi\|_{\infty}ds\bigg\}\\ \le&\max\limits_{i\in\mathcal{J}}\bigg\{\bar{a}_{i}+\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg(\bar{b}_{i}^{\emptyset}\bar{a}_{i}+\bar{b}_{i}^{c}+\sum\limits_{j = 1}^{n}\bar{c}_{ij}L_{j}^{f}+\sum\limits_{j = 1}^{n}\bar{u}_{ij}L_{j}^{f}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\theta}_{ijk}L_{j}^{g}M_{k}^{g}\\ &+\sum\limits_{j = 1}^{n}\bar{\alpha}_{ij}L_{j}^{f}+\sum\limits_{j = 1}^{n}\bar{\beta}_{ij}L_{j}^{f}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{q}_{ijk}L_{j}^{g}M_{k}^{g}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\nu}_{ijk}L_{j}^{g}M_{k}^{g}\bigg)\bigg\}\|\varphi\|_{\infty}\\ = \, &\rho\|\varphi\|_{\infty}\\ < \, &\frac{r}{1-\rho}, \, \, i\in\mathcal{J}. \end{align*}$ $

In addition, by condition $ (A_1) $ and the fact that $ \varphi\in \mathbb{Z} $, we have that for every $ \varepsilon > 0 $, there exists a positive number $ \delta = \delta(\varepsilon) (<\varepsilon) $ such that for any $ h\in \mathbb{R} $ with $ |h| < \delta $, it holds

$\begin{align*} |a_{i}(t+h)-a_{i}(t)|_{1} < \varepsilon, \, \, |\varphi_{i}(\cdot+h)-\varphi_{i}(\cdot)|_{B^p} < \varepsilon, \, \, |\tau_{i}(t+h)-\tau_{i}(t)| < \delta, \, \, i\in \mathcal{J}. \end{align*}$ $

Without loss of generality, in the sequel, we assume that $ h > 0 $, then we deduce that

$\begin{align*} &\|T\varphi(\cdot+h)-T\varphi(\cdot)\|^p_{B^p}\\ \leq\, &2^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t+h)\varphi_{i}(t+h-\tau_{i}(t+h))-a_{i}(t)\varphi_{i}(t-\tau_{i}(t))|^p_{1}dt\bigg\}\\ &+2^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\bigg|\int_{-\infty}^{t+h}e^{-\int_{s}^{t+h}b_{i}^{\emptyset}(u)du}(N^{\varphi})_{i}(s)ds\\ &-\int_{-\infty}^{t}e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}(N^{\varphi})_{i}(s)ds\bigg|^p_{1}dt\bigg\}\\ \leq\, &6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t+h)-a_{i}(t)|^p_{1}|\varphi_{i}(t+h-\tau_{i}(t+h))|^p_{1}dt\bigg\}\\ &+6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t)|^p_{1}|\varphi_{i}(t+h-\tau_{i}(t+h))-\varphi_{i}(t-\tau_{i}(t+h))|^p_{1}dt\bigg\}\\ &+6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t)|^p_{1}|\varphi_{i}(t-\tau_{i}(t+h))-\varphi_{i}(t-\tau_{i}(t))|^p_{1}dt\bigg\}\\ &+4^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\bigg|\int_{-\infty}^{t}\Big|e^{-\int_{s}^{t+h}b_{i}^{\emptyset}(u)du} -e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\Big|(N^{\varphi})_i(s)ds\bigg|^p_{1}dt\bigg\}\\ &+4^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\bigg|\int_{t}^{t+h}e^{-\int_{s}^{t+h}b_{i}^{\emptyset}(u)du}(N^{\varphi})_i(s)ds\bigg|^p_{1}dt\bigg\}\\ \leq\, &6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t+h)-a_{i}(t)|^p_{1}|\varphi_{i}(t+h-\tau_{i}(t+h))|^p_{1}dt\bigg\}\\ &+6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t)|^p_{1}|\varphi_{i}(t+h-\tau_{i}(t+h))-\varphi_{i}(t-\tau_{i}(t+h))|^p_{1}dt\bigg\}\\ &+6^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}|a_{i}(t)|^p_{1}|\varphi_{i}(t-\tau_{i}(t+h))-\varphi_{i}(t-\tau_{i}(t))|^p_{1}dt\bigg\}\\ &+4^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\bigg|\int_{-\infty}^{t}e^{-\underline{b}_{i}^{\emptyset}(t-s)} \bigg|\int_{s}^{t+h}b_{i}^{\emptyset}(u)du-\int_{s}^{t}b_{i}^{\emptyset}(u)du\bigg|(N^{\varphi}(s))_ids\bigg|^p_{1}\bigg\}\\ &+4^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}}\int_{-l}^{l}\bigg|\int_{t}^{t+h}e^{-\int_{s}^{t+h}b_{i}^{\emptyset}(u)du}(N^{\varphi})_i(s)ds\bigg|^p_{1}dt\bigg\}\\ \leq&\max\limits_{i\in\mathcal{J}}\bigg\{6^{p-1}\varepsilon^p\|\varphi\|^p_{\infty} +6^{p-1}\bar{a}_{i}^p\varepsilon^p +6^{p-1}\bar{a}_{i}^p\varepsilon^p+4^{p-1}h^p\bigg[\bigg(\frac{\bar{b}_{i}^{\emptyset}}{\underline{b}_{i}^{\emptyset}}\bigg)^p+\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^p\bigg]\|N^\varphi\|^p_\infty\bigg\}\\ \leq& \max\limits_{i\in\mathcal{J}}\bigg\{6^{p-1}\|\varphi\|^p_{\infty} +6^{p-1}\bar{a}_{i} +6^{p-1}\bar{a}_{i}+4^{p-1}\bigg[\bigg(\frac{\bar{b}_{i}^{\emptyset}}{\underline{b}_{i}^{\emptyset}}\bigg)^p+\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^p\bigg]\|N^\varphi\|^p_\infty\bigg\}\varepsilon^p, \end{align*}$ $

which implies $ T\varphi\in BCB^{p}(\mathbb{R}, \mathcal{A}^{n}) $. Therefore, $ T(\mathbb{Z}^*)\subset \mathbb{Z}^* $.

Next, we will prove that $ T $ is a contraction mapping. In fact, for any $ \varphi, \psi\in\mathbb{Z}^*, i\in\mathcal{J}, $ we have

$\begin{align*} \|T\varphi-T\psi\|_{\infty} \leq\, & ess\sup\limits_{t\in\mathbb{R}}|a_{i}(t)(\varphi_{i}(t-\tau_{i}(t))-\psi_{i}(t-\tau_{i}(t)))|_{1}\\ &+ess\sup\limits_{t\in\mathbb{R}}\bigg|\int_{-\infty}^{t}e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\bigg[-b_{i}^{\emptyset}(s)a_{i}(s)(\varphi_{i}(s-\tau_{i}(s))-\psi_{i}(s-\tau_{i}(s)))\\ &-b_{i}^{c}(s)(\varphi_{i}(s)-\psi_{i}(s))+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}c_{ij}(s)(f_{j}(\vartheta_{j}\varphi_{j}(s))-f_{j}(\vartheta_{j}\psi_{j}(s)))\\ &+\sum\limits_{j = 1}^{n}u_{ij}(s)(f_{j}(\vartheta_{j}\varphi_{j}(s-\sigma_{ij}(s)))-f_{j}(\vartheta_{j}\psi_{j}(s-\sigma_{ij}(s))))\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\theta_{ijk}(s)(g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\\ &-g_{j}(\vartheta_{j}\psi_{j}(s-\delta_{ijk}(s))))g_{k}(\vartheta_{k}\psi_{k}(s-\delta_{ijk}(s))))\\ &+\sum\limits_{j = 1}^{n}\alpha_{ij}(s)(f_{j}(\vartheta_{j}\varphi_{j}(s-\eta_{ij}(s)))-f_{j}(\vartheta_{j}\psi_{j}(s-\eta_{ij}(s))))\\ &+\sum\limits_{j = 1}^{n}\beta_{ij}(s)(f_{j}(\vartheta_{j}\varphi_{j}(s-\eta_{ij}(s)))-f_{j}(\vartheta_{j}\psi_{j}(s-\eta_{ij}(s))))\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}q_{ijk}(s)(g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\\ &-g_{j}(\vartheta_{j}\psi_{j}(s-\delta_{ijk}(s))))g_{k}(\vartheta_{k}\psi_{k}(s-\delta_{ijk}(s))))\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\nu_{ijk}(s)(g_{j}(\vartheta_{j}\varphi_{j}(s-\delta_{ijk}(s)))g_{k}(\vartheta_{k}\varphi_{k}(s-\delta_{ijk}(s)))\\ &-g_{j}(\vartheta_{j}\psi_{j}(s-\delta_{ijk}(s))))g_{k}(\vartheta_{k}\psi_{k}(s-\delta_{ijk}(s)))\bigg]ds\\ \leq\, &\bar{a}_{i}\|\varphi-\psi\|_{\infty}+\int_{-\infty}^{t}e^{\underline{b}_{i}^{\emptyset}(t-s)}\bigg[\bar{b}_{i}^{\emptyset}\bar{a}_{i}+\bar{b}_{i}^{c}+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}\bar{c}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{j = 1}^{n}\bar{u}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\theta}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\alpha}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\beta}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\bar{q}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\nu}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\bigg)\bigg]\|\varphi-\psi\|_{\infty}ds\\ \leq&\max\limits_{i\in\mathcal{J}}\bigg\{\bar{a}_{i}+\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg[\bar{b}_{i}^{\emptyset}\bar{a}_{i}+\bar{b}_{i}^{c}+\vartheta_{i}^{-1}\bigg(\sum\limits_{j = 1}^{n}\bar{c}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{u}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\theta}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\alpha}_{ij}L_{j}^{f}\vartheta_{j}+\sum\limits_{j = 1}^{n}\bar{\beta}_{ij}L_{j}^{f}\vartheta_{j}\\ &+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{q}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}+\sum\limits_{i = 1}^{n}\sum\limits_{k = 1}^{n}\bar{\nu}_{ijk}L_{j}^{g}M_{k}^{g}\vartheta_{j}\bigg)\bigg]\bigg\}\|\varphi-\psi\|_{\infty}\\ = \, &\rho\|\varphi-\psi\|_{\infty}, \, \, i\in\mathcal{J}, \end{align*}$ $

which combined with condition $ (A_{3}) $ means that $ T $ is a contraction mapping. Thereupon, $ T $ has a unique fixed point $ \varphi^{*}\in\mathbb{Z}^*. $

Finally, we will examine that $ \varphi^{*} $ is $ B^p $-almost periodic.

Since $ \varphi^{*}\in\mathbb{Z}^*\subset BCB^p(\mathbb{R}, \mathcal{A}^n) $, for any $ \varepsilon > 0 $, there exists a $ \sigma > 0 (\sigma < \varepsilon) $ such that, for any $ \hbar\in \mathbb{R} $ with $ |\hbar| < \sigma $,

$\begin{align*} |\varphi^*_i(\cdot+\hbar)-\varphi^*_i(\cdot)|_{B^p} < \varepsilon, \, \, i\in \mathcal{J}. \end{align*}$ $

Based on this and $ (A_1) $, there exists $ {\flat} $ such that, for all $ i\in\mathcal{J} $,

$\begin{align} |a_{i}(\cdot+{\flat} )-a_{i}(\cdot)|_{1} < \varepsilon, \, \, \, \, |b_{i}^{\emptyset }(\cdot+{\flat})-b_{i}^{\emptyset }(\cdot)| < \varepsilon, \, \, \, \, |b_{i}^{c }(\cdot+{\flat})-b_{i}^{c}(\cdot)|_{1} < \varepsilon, \end{align}$ $

(3.3)

$\begin{align} |c_{ij}(\cdot+{\flat} )-c_{ij}(\cdot)|_{1} < \varepsilon, \, \, \, \, |u_{ij}(\cdot+{\flat} )-u_{ij}(\cdot)|_{1} < \varepsilon, \, \, \, \, |\gamma_{ij}(\cdot+{\flat} )-\gamma_{ij}(\cdot)|_{B^{p}} < \varepsilon, \end{align}$ $

(3.4)

$\begin{align} |\theta _{ijk}(\cdot+{\flat} )-\theta _{ijk}(\cdot)|_{1} < \varepsilon, \, \, \, \, |\alpha _{ij}(\cdot+{\flat} )-\alpha _{ij}(\cdot)|_{1} < \varepsilon, \, \, \, \, |\beta_{ij}(\cdot+{\flat} )-\beta_{ij}(\cdot)|_{1} < \varepsilon, \end{align}$ $

(3.5)

$\begin{align} |q_{ijk}(\cdot+{\flat} )-q_{ijk}(\cdot)|_{1} < \varepsilon, \, \, \, \, |\nu_{ijk}(\cdot+{\flat} )-\nu_{ijk}(\cdot)|_{1} < \varepsilon, \, \, \, \, |T_{ij}(\cdot+{\flat} )-T_{ij}(\cdot)|_{B^{p}} < \varepsilon, \end{align}$ $

(3.6)

$\begin{align} |\tau_i(t+\flat)-\tau(t)| < \varepsilon, \, \, |\sigma_{ij}(t+\flat)-\sigma_{ij}(t)| < \varepsilon, \, \, |\delta_{ijk}(t+\flat)-\delta_{ijk}(t)| < \varepsilon, \end{align}$ $

(3.7)

$\begin{align} |\eta_{ij}(t+\flat)-\eta_{ij}(t)| < \varepsilon, \, \, |S_{ij}(\cdot+{\flat} )-S_{ij}(\cdot)|_{B^{p}} < \varepsilon, \, \, \, \, |I_{i}(\cdot+{\flat} )-I_{i}(\cdot)|_{B^{p}} < \varepsilon, \end{align}$ $

(3.8)

$\begin{align} |\varphi_{i}^{*}(\cdot-\tau_{i}(\cdot+{\flat}))-\varphi_{i}^{*}(\cdot-\tau_{i}(\cdot))|_{B^p} < \varepsilon, \, \, \, \, |\varphi_{j}^{*}(\cdot-\sigma_{ij}(\cdot+{\flat}))-\varphi_{j}^{*}(\cdot-\sigma_{ij}(\cdot))|_{B^p} < \varepsilon, \end{align}$ $

(3.9)

$\begin{align} |\varphi_{i}^{*}(\cdot-\eta_{ij}(\cdot+{\flat}))-\varphi_{i}^{*}(\cdot-\eta_{ij}(\cdot))|_{B^p} < \varepsilon, \, \, \, \, |\varphi_{j}^{*}(\cdot-\delta_{ijk}(\cdot+{\flat}))-\varphi_{j}^{*}(\cdot-\delta_{ijk}(\cdot))|_{B^p} < \varepsilon, \end{align}$ $

(3.10)

$\begin{align} |\varphi_{k}^{*}(\cdot-\delta_{ijk}(\cdot+{\flat}))-\varphi_{k}^{*}(\cdot-\delta_{ijk}(\cdot))|_{B^p} < \varepsilon. \end{align}$ $

(3.11)

Then, we deduce that

$\begin{align*} & \|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^p}^p \\ \leq\, & 2^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l | a_i(t+{\flat}) \varphi_i^*(t+{\flat}-\tau_i(t+{\flat}))-a_i(t) \varphi_i^*(t-\tau_i(t))|_1 ^p d t\bigg\} \\ & +2^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg|\int_{-\infty}^{t+{\flat}} e^{-\int_s^{t+{\flat}} b_i^{\emptyset}(u+{\flat}) d u}(N^{\varphi^*}) i(s) d s \\ & -\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u) d u}(N^{\varphi^*})_i(s) d s\bigg|_1 ^p d t\bigg\} \\ \leq\, & 4^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l|a_i(t+{\flat})(\varphi_i^*(t+{\flat}-\tau_i(t+{\flat}))-\varphi_i^*(t-\tau_i(t)))|_1^p d t\bigg\} \\ & +4^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l|(a_i(t+{\flat})-a_i(t)) \varphi_i^*(t-\tau_i(t))|_1^p d t\bigg\} \\ & +2^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} (N^{\varphi^*})_i(s+{\flat}) d s \\ &-\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u) d u}(N^{\varphi^*})_i(s) d s\bigg|_1 ^p d t\bigg\} \\ \leq\, & 4^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}|a_i(t+{\flat})( \varphi_i^*(t+{\flat}-\tau_i(t+{\flat}))-\varphi_i^*(t-\tau_i(t)))|_1 ^p d t\bigg\} \\ & +4^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}|(a_i(t+{\flat})-a_i(t)) \varphi_i^*(t-\tau_i(t))|_1^p d t\bigg\} \\ &+70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}\bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} a_i(s+{\flat}) b_i^{\emptyset}(s+{\flat})\\ &\times(\varphi_i^*(s+{\flat}-\tau_i(s+{\flat}))-\varphi_i^*(s-\tau_i(s))) d s\bigg|_1^p d t\bigg\}\\ & +70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg|\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}(a_i(s+{\flat})-a_i(s)) b_i^{\emptyset}(s+{\flat})\\ &\times\varphi_i^*(s-\tau_i(s)) d s\bigg|_1 ^p d t\bigg\} \\ & +70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} a_i(s)(b_i^{\emptyset}(s+{\flat})\\ &-b_i^{\emptyset}(s))\varphi_i^*(s-\tau_i(s)) d s\bigg|_1 ^p d t\bigg\} \\ &+70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-1}^l\bigg|\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset} (u) du}\Big|\\ &\times a_i(s) b_i^{\emptyset}(s) \varphi_i^*(s-\tau_{i}(s)) d s\bigg|_1 ^p d t\bigg\} \\ & +70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}(b_i^c(s+{\flat}) \varphi_i^*(s+{\flat})-b_i^c(s) \varphi_i^*(s))d s\bigg|_1 ^p d t\bigg\} \\ & +70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big|b_i^c(s) \varphi_i^*(s) d s\bigg|_1 ^p d t \bigg\}\\ \end{align*}$ $

$\begin{align*} &\;\;\;\;+70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \sum\limits_{j = 1}^n c_{i j}(s+{\flat})(f_j(\vartheta_j \varphi_j^*(s+{\flat}))\\ &\;\;\;\;-f_j(\vartheta_j \varphi_i^*(s)))d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \sum\limits_{j = 1}^n(c_{i j}(s+{\flat})\\ &\;\;\;\;-c_{i j}(s))f_j(\vartheta_j \varphi_j^*(s)) d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big| \\ &\;\;\;\;\times\sum\limits_{j = 1}^n c_{i j}(s) f_j(\vartheta_j \varphi_j^*(s)) d s\bigg|_1 ^p d t\bigg\}\\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg|\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \sum\limits_{j = 1}^n u_{i j}(s+{\flat})\\ &\;\;\;\;\times(f_j(\vartheta_j \varphi_j^*(s+{\flat}-\sigma_{i j}(s+{\flat})))-f_j(\vartheta_j \varphi_j^*(s-\sigma_{i j}(s)))) d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg| \int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \sum\limits_{j = 1}^n(u_{i j}(s+{\flat})-u_{i j}(s))\\ &\;\;\;\;\times f_j(\vartheta_j \varphi_j^*(s-\sigma_{i j}(s))) d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big|\\ &\;\;\;\;\times\sum\limits_{j = 1}^n u_{i j}(s) f_j(\vartheta_j \varphi_j^*(s-\sigma_{i j}(s))) d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \bigg(\sum\limits_{j = 1}^n \gamma_{i j}(s+{\flat}) \mu_j(s+{\flat}) \\ &\;\;\;\;-\sum\limits_{j = 1}^n \gamma_{i j}(s) \mu_j(s)\bigg)d s\bigg|_1 ^p d t \bigg\}\\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t \Big|e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} -e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\Big|\sum\limits_{j = 1}^n \gamma_{i j}(s) \mu_j(s)d s\bigg|_1 ^p d t \bigg\}\\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg |\int_{-\infty}^t e^{-\int _s^t b_i^{\emptyset}( u+{\flat}) d u} \sum\limits_{j = 1}^n \sum\limits_{k = 1}^n \theta_{i j k}(s+{\flat})\\ &\;\;\;\;\times(g_j(\vartheta_j \varphi_j^*(s+{\flat}-\delta_{i j k}(s+{\flat}))) g_k ( \vartheta_k \varphi_k^*(s+{\flat}-\delta_{i j k}(s+{\flat}))) \\ &\;\;\;\; -g_j(\vartheta_j \varphi_j^*(s-\delta_{i j k}(s))) g_k(\vartheta_k \varphi_k^*(s-\delta_{i j k}(s )))) d s\bigg|_1 ^p d t\bigg\} \\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg|\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\theta_{i j k}(s+{\flat})\\ &\;\;\;\; -\theta_{i j k}(s)) g_j(\vartheta_j \varphi_j^*(s-\delta_{ijk} k(s))) g_k(\vartheta_k \varphi_k^*(s-\delta_{i j k}(s) )) d s\bigg|_1 ^p d t\bigg\}\\ &\;\;\;\; +70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg|\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big| \\ &\;\;\;\; \times \sum\limits_{j = 1}^n \sum\limits_{k = 1}^n \theta_{i j k}(s) g_j(\vartheta_j \varphi_j^*(s-\delta_{i j k}(s))) g_k(\vartheta_k \varphi_k^*(s-\delta_{i j k}(s))) d s\bigg|_1 ^p d t\bigg\} \\ \end{align*}$ $

Furthermore, based on inequalities (3.3)–(3.11), the Hölder inequality, and the Fubini theorem, we can deduce that

$\begin{align*} K_{1}\leq\, &8^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \bigg[\int_{-l}^{l}|a_i(t+{\flat})|^{p}(| \varphi_i^*(t+{\flat}-\tau_i(t+{\flat}))\\ &-\varphi_{i}^{*}(t-\tau_{i}(t+{\flat}))|_{1}^{p}+|\varphi_{i}^{*}(t-\tau_{i}(t+{\flat}))-\varphi_i^*(t-\tau_i(t))|_{1}^{p}) d t\bigg\}\\ \leq\, &8^{p-1} \max _{i \in \mathcal{J}}\bigg\{(\bar{a}_{i})^{p}\frac{1}{1-\bar{\tau^{\prime}}_{i}}| \varphi_i^*(t+{\flat})-\varphi_{i}^{*}(t)|_{B^{p}}^{p}+(\bar{a}_{i})^{p}\varepsilon^{p}\bigg\}, \\ K_{2}\leq\, &4^{p-1} \max _{i \in \mathcal{J}}\{\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\}, \\ K_{3} \leq\, &140^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}\bigg[\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}ds\bigg]^{\frac{p}{q}}\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}\\ &\times|a_i(s+{\flat})b_i^{\emptyset}(s+{\flat})|_1(|\varphi_i^*(s+{\flat}-\tau_i(s+{\flat}))-\varphi_i^*(s-\tau_i(s+{\flat}))|_1^p\\ &+|\varphi_i^*(s-\tau_i(s+{\flat}))-\varphi_i^*(s-\tau_i(s))|_{1}^{p}) d s d t\bigg\}\\ \leq\, &140^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p}{q}}(\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p}\bigg[\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}\int_{-\infty}^{t-\tau_{i}(t+{\flat})}\frac{e^{-\underline{b}_{i}^{\emptyset}(t-s-\bar{\tau}_{i})}}{1-\bar{\tau}_{i}^{\prime}}\\ &\times|\varphi_{i}^{*}(s+{\flat})-\varphi_{i}^{*}(s)|_{1}^{p}dsdt\bigg]+\frac{\varepsilon^{p}}{\underline{b}_{i}^{\emptyset}}\bigg\}\\ \leq\, &140^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p}{q}}(\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p} \bigg[\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\tau}_{i}}}{1-\bar{\tau^{\prime}}_{i}}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^{l}\int_{-\infty}^{t}e^{-\underline{b}_{i}^{\emptyset}(t-s)}\\ &\times|\varphi_{i}^{*}(s+{\flat})-\varphi_{i}^{*}(s)|_{1}^{p}dsdt\bigg]+\frac{\varepsilon^{p}}{\underline{b}_{i}^{\emptyset}}\bigg\}\\ \leq\, &140^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p}{q}}(\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p} \bigg[\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\tau}_{i}}}{1-\bar{\tau^{\prime}}_{i}}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-\infty}^{l}(2l)^{-1}\int_{s-2l}^{s}e^{-\underline{b}_{i}^{\emptyset}(l-s)}\\ &\times|\varphi_{i}^{*}(t+{\flat})-\varphi_{i}^{*}(t)|_{1}^{p}dtds\bigg]+\frac{\varepsilon^{p}}{\underline{b}_{i}^{\emptyset}}\bigg\}\\ \leq\, &140^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p} \bigg(\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\tau}_{i}}}{1-\bar{\tau^{\prime}}_{i}}\|\varphi^{*}(t+{\flat})-\varphi^{*}(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg)\bigg\}, \\ K_{4}\leq\, &70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg(\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}ds\bigg)^{\frac{p}{q}}\\ &\times\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}|(a_i(s+{\flat})-a_i(s)) b_i^{\emptyset}(s+{\flat})\varphi_i^*(s-\tau_i(s)) |_{1}^{p}d s d t\bigg\} \\ \leq\, & 70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{b}_{i}^{\emptyset})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}.\\ K_{5}\leq\, &70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{a}_{i})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}, \\ K_{6} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-1}^l\bigg[\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset} (u) du}\Big|ds\bigg]^{\frac{p}{q}}\\ &\times \int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset} (u) du}\Big||a_i(s) b_i^{\emptyset}(s) \varphi_i^*(s-\tau_{i}(s)) |_1 ^p d s d t\bigg\} \\ \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} (\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p}\|\varphi^*\|^{p}_{\infty}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{7} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg[\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}ds\bigg]^{\frac{p}{q}}\\ &\times\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}|b_i^c(s+{\flat}) \varphi_i^*(s+{\flat})-b_i^c(s) \varphi_i^*(s)|_{1}^{p}d sd t\bigg\} \\ \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\bigg( \frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p}{q}}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-\infty}^l e^{-\underline{b}_i^{\emptyset}(l-s) }\int_{s-2l}^{s}((\bar{b}_i^c)^{p} \\ \end{align*}$ $

$\begin{align*} &\times|\varphi_i^*(t+{\flat})-\varphi_i^*(t)|_{1}^{p}+|b_i^c(t+{\flat})-b_i^c(t)|_{1}^{p}\|\varphi^*\|_\infty^{p}) d td s\bigg\} \\ \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\bigg( \frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}((\bar{b}_i^c)^{p} \|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\|\varphi^*\|_{\infty}^{p}\varepsilon^{p})\bigg\}, \\ K_{8}\leq\, & 70^{p-1} \max _{i \in\mathcal{J}}\bigg\{(\frac{1}{\underline{b}_{i}^{\emptyset}})^{\frac{2(p+q)}{q}} (\bar{b}_{i}^{\emptyset})^{p}\|\varphi^*\|^{p}_{\infty}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{9}\leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg[ \int_{-\infty}^te^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}ds \bigg]^{\frac{p}{q}}\\ &\times\int_{-\infty}^t e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} \bigg|\sum\limits_{j = 1}^n c_{i j}(s+{\flat})(f_j(\vartheta_j \varphi_j^*(s+{\flat}))-f_j(\vartheta_j \varphi_i^*(s)))\bigg|_1 ^pd sd t\bigg\} \\ \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{ \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p}{q}}\bigg(\sum\limits_{j = 1}^{n}(\bar{c}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-\infty}^l e^{- \underline{b}_i^{\emptyset}(l-s)}\\ &\times\int_{s-2l}^{s} |\varphi_j^*(t+{\flat}))- \varphi_i^*(t)|_1 ^p d t d s\bigg\} \\ \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{ \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^{n}(\bar{c}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*(t+{\flat}))- \varphi^*(t)\|_{B^{p}} ^p\bigg\}, \\ K_{10}\leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}, \\ K_{11}\leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{c}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{12} \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+1}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{u}_{ij})^{q}\bigg)^{\frac{p}{q}} \sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\\ &\times\bigg[\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\sigma}_{ij}}}{1-\bar{\sigma^{\prime}}_{ij}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg]\bigg\}, \\ \end{align*}$ $

$\begin{align*} K_{13} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}n^{\frac{p}{q}}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}, \\ K_{14}\leq\, &70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg[\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big|^{\frac{q}{p}}ds\bigg]^{\frac{p}{q}}\\ &\times\int_{-\infty}^t\Big| e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u}-e^{-\int_s^t b_i^{\emptyset}(u) d u} \Big|^{\frac{p}{q}}\bigg|\sum\limits_{j = 1}^n u_{i j}(s) f_j(\vartheta_j \varphi_j^*(s-\sigma_{i j}(s))) \bigg|_1 ^pd s d t\bigg\} \\ \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}}\bigg(\sum\limits_{j = 1}^{n}(\bar{u}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{15} \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{ \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg[\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{p}+\sum\limits_{j = 1}^{n}|\gamma_{ij}|^{p}_\infty\bigg]\varepsilon^{p} \bigg\}, \\ K_{16}\leq\, &70^{p-1} \vartheta_i^{-p} \max _{i \in\mathcal{J}}\bigg\{\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l\bigg[\int_{-\infty}^t \Big|e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} -e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\Big|^{\frac{q}{p}}ds\bigg]^{\frac{p}{q}}\\ &\times\int_{-\infty}^t \Big|e^{-\int_s^t b_i^{\emptyset}(u+{\flat}) d u} -e^{-\int_{s}^{t}b_{i}^{\emptyset}(u)du}\Big|^{\frac{p}{q}}\bigg|\sum\limits_{j = 1}^n \gamma_{i j}(s) \mu_j(s)\bigg|_1 ^pd s d t \bigg\}\\ \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}|\gamma_{ij}|^{q}_\infty\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{17}\leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\overline{\lim\limits_{l\to\infty}}{ \frac{1}{2l}} \int_{-l}^l \bigg(\int_{-\infty}^t e^{-\int _s^t b_i^{\emptyset}( u+{\flat}) d u} ds\bigg)^{\frac{p}{q}}\int_{-\infty}^t e^{-\int _s^t b_i^{\emptyset}( u+{\flat}) d u} \\ &\times\bigg(\sum\limits_{j = 1}^n\sum\limits_{k = 1}^n \bar{\theta}_{i j k}(M^g_k|g_j(\vartheta_j \varphi_j^*(s+{\flat}-\delta_{i j k}(s+{\flat})))-g_j(\vartheta_j \varphi_j^*(s-\delta_{i j k}(s)))|_{1} \\ &+ M_j^g|g_k ( \vartheta_k \varphi_k^*(s+{\flat}-\delta_{i j k}(s+{\flat})))- g_k(\vartheta_k \varphi_k^*(s-\delta_{i j k}(s )))) |_{1})d sd t\bigg)^p\bigg\} \\ \leq\, &280^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg[\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\theta}_{i j k}M_k^gL_j^g\vartheta_j)^p+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\theta}_{ijk}M_j^gL_{k}^{g}\vartheta_{k})^{p}\bigg]\\ &\times\bigg(\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta^{\prime}}_{ijk}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg)\bigg\}, \\ K_{18} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\sum\limits_{k = 1}^{n} \sum\limits_{j = 1}^{n}(M_j^gM_{k}^{g})^{p}\varepsilon^{p}\bigg\}, \\ K_{19}\leq\, &70^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}} \sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\theta}_{ijk}M_j^gM_k^g)^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{20}\leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n( \bar{\alpha}_{i j})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p} \\ &\times\bigg[\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\eta}_{ij}}}{1-\bar{\eta^{\prime}}_{ij}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg]\bigg\}, \\ K_{21} \leq\, &70^{p-1} \max _{i \in \mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}, \\ \end{align*}$ $

$\begin{align*} K_{22} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{\alpha}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{23} \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n (\bar{\beta}_{i j})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\\ &\times\bigg[\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\eta}_{ij}}}{1-\bar{\eta^{\prime}}_{ij}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p} + \varepsilon^{p}\bigg] \bigg\}, \\ K_{24} \leq\, &70^{p-1} \max _{i \in \mathcal{J}}\bigg\{ \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p}\bigg\}, \\ K_{25} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{\beta}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{26} \leq\, &280^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_k^gL_{j}^{g}\vartheta_{j})^{p}+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_j^gL_k^g\vartheta_k)^p\bigg)\\ &\times\bigg(\frac{e^{\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta}_{ijk}^{\prime}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg)\bigg\}, \\ K_{27} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\sum\limits_{k = 1}^{n} \sum\limits_{j = 1}^{n}(M_j^gM_{k}^{g})^{p}\varepsilon^{p}\bigg\}, \\ K_{28} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}(\frac{1}{\underline{b}_{i}^{\emptyset}})^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}}\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{q}_{ijk}M_j^gM_k^g)^p\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{29} \leq\, &280^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{k}^{g}L_{j}^{g}\vartheta_{j})^{p}+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{j}^{g}L_{k}^{g}\vartheta_{k})^{p}\bigg)\\ &\times\bigg(\frac{e^{\frac{p}{q}\underline{b}_{i}^{\emptyset}\bar{\delta}_{ijk}}}{1-\bar{\delta^{\prime}}_{ijk}}\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+\varepsilon^{p}\bigg)\bigg\}, \\ K_{30} \leq\, &70^{p-1}\max\limits_{i\in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}} n^{\frac{2p}{q}}\sum\limits_{k = 1}^{n}\sum\limits_{j = 1}^{n}(M_j^gM_k^g )^p\varepsilon^{p}\bigg\}, \\ K_{31}\leq\, & 70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}}\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\nu}_{ijk}M_j^gM_k^g)^p\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{32} \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg[\bigg(\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{q}\bigg)^{\frac{p}{q}} +\bigg(\sum\limits_{j = 1}^{n}(|T_{ij}|_\infty)^{q}\bigg)^{\frac{p}{q}}\bigg]\varepsilon^{p} \bigg\}, \\ K_{33} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(|T_{ij}|_\infty\bar{\mu}_{j})^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{34} \leq\, &140^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+1}{q}}\bigg[\bigg(\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{q}\bigg)^{\frac{p}{q}} +\bigg(\sum\limits_{j = 1}^{n}(|S_{ij}|_\infty)^{q}\bigg)^{\frac{p}{q}}\bigg]\varepsilon^{p} \bigg\}, \\ K_{35} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{ \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(|S_{ij}|_\infty\bar{\mu}_{j})^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}, \\ K_{36} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\varepsilon^{p} \bigg\}, \\ K_{37} \leq\, &70^{p-1} \max _{i \in\mathcal{J}}\bigg\{\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}}|I_{i}|_\infty^{p}\varepsilon^{\frac{p+q}{q}}\bigg\}. \end{align*}$ $

From the above estimates, it follows that

$\begin{align} &\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p} \leq P\|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p}+Q\varepsilon^{p}, \end{align}$ $

(3.12)

where $ P $ is defined in condition $ (A_4) $ and

$\begin{align*} Q = \, & 2^{p-1}\max _{i \in \mathcal{J}}\bigg\{4^{p-1}(\bar{a}_{i})^{p}\varepsilon^{p-1}+2^{p-1} \bigg( \frac{r}{1-\rho}\bigg)^p\varepsilon^{p-1}+70^{p-1}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p}\varepsilon^{p-1}\\ &+35^{p-1} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{b}_{i}^{\emptyset})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{p-1}+35^{p-1} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}(\bar{a}_{i})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{p-1}\\ &+35^{p-1} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} (\bar{a}_{i}\bar{b}_{i}^{\emptyset})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{\frac{p}{q}}+70^{p-1} \bigg( \frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{p-1}\\ &+35^{p-1} (\frac{1}{\underline{b}_{i}^{\emptyset}})^{\frac{2(p+q)}{q}} (\bar{b}_{i}^{\emptyset})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{\frac{p}{q}}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\\ &\times\bigg(\frac{r}{1-\rho}\bigg)^{p} \varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{c}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{\frac{p}{q}}\\ &+ 70^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+1}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{u}_{ij})^{q}\bigg)^{\frac{p}{q}} \sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\\ &\times \sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}n^{\frac{p}{q}}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{u}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\\ &\times \bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{\frac{p}{q}}+70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg[\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{p}+\sum\limits_{j = 1}^{n}|\gamma_{ij}|^{p}_\infty\bigg]\varepsilon^{p-1} \\ &+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}|\gamma_{ij}|^{q}_\infty\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{p}\varepsilon^{\frac{p}{q}}+140^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\\ &\times\bigg[\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\theta}_{i j k}M_k^gL_j^g\vartheta_j)^p+\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\theta}_{ijk}M_j^gL_{k}^{g}\vartheta_{k})^{p}\bigg]\varepsilon^{p-1}\bigg\} +35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\\ &\times n^{\frac{2p}{q}}\sum\limits_{k = 1}^{n} \sum\limits_{j = 1}^{n}(M_j^gM_{k}^{g})^{p}\varepsilon^{p-1}+35^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}} \sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\theta}_{ijk}M_j^gM_k^g)^{p}\varepsilon^{\frac{p}{q}}\\ &+70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n( \bar{\alpha}_{i j})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\varepsilon^{p-1}\bigg]+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\\ &\times n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} \bigg(\sum\limits_{j = 1}^{n}(\bar{\alpha}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\\ &\times \bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{\frac{p}{q}}+700^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg(\sum\limits_{j = 1}^n (\bar{\beta}_{i j})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\varepsilon^{p-1} \end{align*}$ $

$\begin{align*} &\;\;\;\;\;\;+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\bigg(\frac{r}{1-\rho}\bigg)^{p}\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}}\\ &\;\;\;\;\;\;\times\bigg(\sum\limits_{j = 1}^{n}(\bar{\beta}_{ij})^{q}\bigg)^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(L_{j}^{f}\vartheta_{j})^{p}\|\varphi^*\|_{\infty}^{p}\varepsilon^{\frac{p}{q}}+140^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\\ &\;\;\;\;\;\;\times\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_k^gL_{j}^{g}\vartheta_{j})^{p}+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n(\bar{q}_{ijk}M_j^gL_k^g\vartheta_k)^p\bigg)\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\\ &\;\;\;\;\;\;\times n^{\frac{2p}{q}} \sum\limits_{k = 1}^{n}\sum\limits_{j = 1}^{n}(M_j^gM_{k}^{g})^{p}\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}(\frac{1}{\underline{b}_{i}^{\emptyset}})^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}}\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{q}_{ijk}M_j^gM_k^g)^p\varepsilon^{\frac{p}{q}}\\ &\;\;\;\;\;\;+140^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}n^{\frac{2p}{q}}\bigg(\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{k}^{g}L_{j}^{g}\vartheta_{j})^{p}+\sum\limits_{j = 1}^n \sum\limits_{k = 1}^n( \bar{\nu}_{ijk}M_{j}^{g}L_{k}^{g}\vartheta_{k})^{p}\bigg)\varepsilon^{p-1}\\ &\;\;\;\;\;\;+35^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}} n^{\frac{2p}{q}}\sum\limits_{k = 1}^{n}\sum\limits_{j = 1}^{n}(M_j^gM_k^g )^p\varepsilon^{p-1}+35^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{p}}n^{\frac{2p}{q}}\\ &\;\;\;\;\;\;\times\sum\limits_{j = 1}^{n}\sum\limits_{k = 1}^{n}(\bar{\nu}_{ijk}M_j^gM_k^g)^p\varepsilon^{\frac{p}{q}}+70^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\bigg[\bigg(\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{q}\bigg)^{\frac{p}{q}}\\ &\;\;\;\;\;\;+\bigg(\sum\limits_{j = 1}^{n}(|T_{ij}|_\infty)^{q}\bigg)^{\frac{p}{q}}\bigg]\varepsilon^{p-1}+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(|T_{ij}|_\infty\bar{\mu}_{j})^{p}\varepsilon^{\frac{p}{q}}\\ &\;\;\;\;\;\;+70^{p-1} \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+1}{q}}\bigg[\bigg(\sum\limits_{j = 1}^{n}(\bar{\mu}_{j})^{q}\bigg)^{\frac{p}{q}} +\bigg(\sum\limits_{j = 1}^{n}(|S_{ij}|_\infty)^{q}\bigg)^{\frac{p}{q}}\bigg]\varepsilon^{p-1} \\ &\;\;\;\;\;\;+35^{p-1} \vartheta_i^{-p} \bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}} n^{\frac{p}{q}}\sum\limits_{j = 1}^{n}(|S_{ij}|_\infty\bar{\mu}_{j})^{p}\varepsilon^{\frac{p}{q}}+35^{p-1}\vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{p+q}{q}}\varepsilon^{p-1} \\ &\;\;\;\;\;\;+35^{p-1} \vartheta_i^{-p}\bigg(\frac{1}{\underline{b}_{i}^{\emptyset}}\bigg)^{\frac{2(p+q)}{q}}|I_{i}|_\infty^{p}\varepsilon^{\frac{p}{q}}\bigg\}. \end{align*}$ $

Hence, by (3.12) and $ (A_4) $, it holds that

$\begin{align*} \|\varphi^*(t+{\flat})-\varphi^*(t)\|_{B^{p}}^{p} \leq\frac{Q\varepsilon^{p}}{1-P}, \end{align*}$ $

which implies that $ \varphi^{*} $ is $ B^p $-almost periodic. The proof is finished. □

Remark 3.1. Although we can prove that $ W = (L^\infty(\mathbb{R}, \mathcal{A})\cap B_{AP}^p(\mathbb{R}, \mathcal{A}), \|\cdot\|_\infty) $ is a Banach space, we still cannot directly use the fixed point theorem to determine the existence of almost periodic solutions for (1.1). Because there are higher-order terms in system (1.1), and $ W $ is not an algebra, we cannot prove that operator $ T $ is a self mapping.

It is easy to prove the following stability results using the same method as the proof of Theorem 4.1 in ^[33] or the proof of Theorem 15 in ^[31].

Theorem 3.2. Assume that $ (A_{1})-(A_{4}) $ hold. Then system (1.1) possesses a unique Besicovitch almost periodic solution, which is globally exponentially stable, i.e., if $ \bar{x} $ is the Besecovitch almost periodic solution with initial value $ \bar{\varphi} $ and $ x(t) $ is an arbitrary solution of system (1.1) with initial value $ \varphi $, then there exist positive numbers $ \zeta > 0 $ and $ N > 0 $ satisfying

$ |x(t)-\bar{x}(t)|_{1}\leq N\|\varphi-\bar{\varphi}\|_{\varrho}e^{-\zeta t}, \, \, \, t > 0, $

in which $ \|\varphi-\bar{\varphi}\|_{\varrho} = \max\limits_{i\in\mathcal{J}}\bigg\{\sup\limits_{t\in[- \varrho, 0]}|\varphi_{i}(t)-a_{i}(t)\varphi_{i}(t)-(\bar{\varphi}_{i}(t)-a_{i}(t)\bar{\varphi}_{i}(t))|_{1}\bigg\} $.

4. Numerical examples

In this section, we provide an example to demonstrate the validity of the results obtained in this paper.

Example 4.1. In system (1.1), let $ m = 3, n = 2 $, and for $ i, j, k = 1, 2 $, take the coefficients are as follows:

$\begin{align*} x_i(t) = \, &e_0x_i^0(t)+e_1x_i^1(t)+e_2x_i^2(t)+e_3x_i^3(t)+e_{12}x_i^{12}(t)+e_{13}x_i^{13}(t)+e_{23}x_i^{23}(t)+e_{123}x_i^{123}(t), \\ f_j(x) = \, &\frac{1}{100}e_{0}\sin{x_j^{12}}+\frac{3}{250}\sin(x_j^{12}+x_j^{123})e_{1}+\frac{1}{168}\sin(x_j^{1}+x_j^{13})e_{2}+\frac{1}{53}\sin(x_j^{2}+x_j^{123})e_{3}\\ &+\frac{1}{125}e_{12}\arctan{x_j^{3}}+\frac{1}{156}\sin(x_j^{3}+x_j^{123})e_{13}+\frac{1}{150}e_{23}\tanh{x_j^{12}}+\frac{1}{40}\sin(x_j^{1}+x_j^{12}\\ &+x_j^{123})e_{123}, \\ g_j(x) = \, &\frac{1}{48}e_{0}\sin{x_j^{13}}+\frac{1}{153}\sin(x_j^{12}+x_j^{23})e_{1}+\frac{1}{150}e_{2}\arctan{x_j^{13}}+\frac{1}{120}\sin(x_j^{12}+x_j^{123})e_{3}\\ &+\frac{1}{150}\sin(x_j^{12}-x_j^{23})e_{12}+\frac{1}{250}\sin(x_j^{2}+x_j^{13})e_{13}+\frac{1}{57}e_{23}\sin{x_j^{13}}+\frac{1}{60}\sin(x_j^{0}+x_j^{3}\\ &+x_j^{23})e_{123}, \\ a_1(t) = \, &(0.01+0.004\sin{t})e_0+(0.01+0.001\sin{\sqrt{6}t})e_1+(0.01+0.001\cos{\sqrt{3}t})e_2\\ &+(0.01+0.002\sin{\sqrt{5}t})e_3+(0.01+0.003\sin{t})e_{12}+(0.01+0.002\cos{t})e_{13}\\ &+(0.01+0.002\cos{\sqrt{2}t})e_{23}+(0.01+0.002\sin{\sqrt{2}t})e_{123}, \\ a_2(t) = \, &(0.01+0.002\sin{t})e_0+(0.01+0.002\cos{3t})e_1+(0.01+0.001\sin{\sqrt{2}t})e_2\\ &+(0.01+0.001\sin{\sqrt{7}t})e_3+(0.01+0.001\sin{\sqrt{5}t})e_{12}+(0.01+0.002\cos{t})e_{13}\\ &+(0.01+0.002\cos{\sqrt{5}t})e_{23}+(0.01+0.002\sin{\sqrt{3}t})e_{123}, \\ b_{1}(t) = \, & (10+0.05\sin t)e_0+(0.2+0.01\cos\sqrt{2}t)e_1+(0.2+0.02\sin t)e_2\\ &+(0.2+0.01\sin\sqrt{3}t)e_3+(0.2+0.06\sin3t)e_{12}+(0.2+0.05\sin2t)e_{13}\\ &+(0.2+0.01\sin t)e_{23}+(0.2+0.01\cos\sqrt{3}t)e_{123}, \\ b_2(t) = \, & (0.2+0.01\cos\sqrt{3}t)e_0+(0.2+0.02\sin t)e_1+(0.2+0.07\cos t)e_2\\ &+(0.2+0.05\cos\sqrt{5}t)e_3+(10+0.05\sin\sqrt{5}t)e_{12}+(0.2+0.02\cos\sqrt{5}t)e_{13}\\ &+(0.2+0.01\sin t)e_{23}+(0.2+0.01\sin7t)e_{123}, \\ c_{11}(t) = \, &0.01e_{0}\sin{2t}+0.02e_{3}\sin{2t}+0.02e_{23}\cos{\sqrt{2}t}+0.03e_{123}\cos{11t}, \\ c_{12}(t) = \, &0.01e_{0}\sin{\sqrt{5}t}+0.02e_{2}\cos^{2}{3t}+0.01e_{3}\sin{5t}+0.03e_{12}\sin{\sqrt{3}t}, \\ c_{21}(t) = \, &0.01e_{0}\sin{6t}+0.02e_{3}\cos{\sqrt{2}t}+0.03e_{23}\cos{\sqrt{3}t}+0.03e_{123}\cos^{2}{2t}, \\ c_{22}(t) = \, &0.01e_{0}\sin^{2}{7t}+0.04e_{2}\cos{6t}+0.04e_{3}\sin{\sqrt{5}t}+0.03e_{12}\cos{7t}, \\ u_{11}(t) = \, &0.02e_{0}\sin{4t}+0.01e_{3}\cos{\sqrt{2}t}+0.04e_{23}\cos{\sqrt{3}t}+0.03e_{123}\sin^{2}{2t}, \\ u_{12}(t) = \, &0.02e_{0}\cos{9t}+0.03e_{2}\cos^{2}{3t}+0.04e_{3}\sin{5t}+0.01e_{12}\sin{3t}, \\ u_{21}(t) = \, &0.02e_{0}\sin{3t}+0.03e_{3}\cos{\sqrt{3}t}+0.01e_{23}\cos{\sqrt{3}t}+0.02e_{123}\sin^{2}{7t}, \\ u_{22}(t) = \, &0.02e_{0}\cos{t}+0.03e_{2}\sin^{2}{5t}+0.03e_{3}\sin{3t}+0.04e_{12}\cos{2t}, \\ \alpha_{11}(t) = \, &0.01e_{3}\sin{3t}+0.04e_{12}\cos{5t}+0.02e_{13}\cos{\sqrt{5}t}+0.03e_{123}\sin^{2}{3t}, \\ \alpha_{12}(t) = \, &0.01e_{2}\sin{\sqrt{3}t}+0.02e_{3}\sin{4t}+0.03e_{12}\sin{7t}+0.01e_{123}\cos^{2}{5t}, \\ \alpha_{21}(t) = \, &0.03e_{3}\cos{4t}+0.01e_{12}\sin{\sqrt{5}t}+0.03e_{23}\cos{4t}+0.01e_{123}\cos^{2}{5t}, \\ \alpha_{22}(t) = \, &0.01e_{0}\cos{\sqrt{5}t}+0.04e_{3}\sin{3t}+0.03e_{12}\sin^{2}{3t}+0.02e_{123}\cos{4t}, \\ \beta_{11}(t) = \, &0.01e_{0}\cos{5t}+0.02e_{1}\sin{3t}+0.01e_{2}\sin{\sqrt{7}t}+0.03e_{123}\cos{3t}, \\ \beta_{12}(t) = \, &0.03e_{0}\sin{7t}+0.02e_{1}\sin{t}+0.03e_{2}\sin{\sqrt{3}t}+0.02e_{23}\sin{3t}, \\ \beta_{21}(t) = \, &0.03e_{0}\cos{7t}+0.02e_{1}\sin{\sqrt{5}t}+0.04e_{2}\sin{\sqrt{6}t}+0.01e_{123}\sin{5t}, \\ \beta_{22}(t) = \, &0.03e_{0}\cos{2t}+0.02e_{1}\cos{3t}+0.04e_{2}\cos{\sqrt{2}t}+0.02e_{23}\cos{3t}, \\ \theta_{111}(t) = \, &0.03e_{0}\sin{\sqrt{2}t}+0.02e_{1}\cos{\sqrt{5}t}+0.01e_{2}\sin{\sqrt{7}t}+0.02e_{12}\cos{2t}, \\ \theta_{112}(t) = \, &0.04e_{0}\cos{\sqrt{5}t}+0.02e_{1}\sin{3t}+0.03e_{2}\sin{\sqrt{3}t}+0.03e_{13}\cos{\sqrt{3}t}, \\ \theta_{121}(t) = \, &0.02e_{0}\sin{\sqrt{3}t}+0.02e_{2}\cos{3t}+0.04e_{12}\sin{\sqrt{6}t}+0.01e_{23}\cos{3t}, \\ \theta_{122}(t) = \, &0.03e_{0}\sin{\sqrt{5}t}+0.04e_{2}\cos{5t}+0.03e_{12}\sin{\sqrt{5}t}+0.02e_{23}\cos{4t}, \\ \theta_{211}(t) = \, &0.03e_{0}\sin{3t}+0.03e_{1}\cos{t}+0.01e_{2}\cos{4t}+0.01e_{12}\sin{5t}, \\ \theta_{212}(t) = \, &0.04e_{0}\cos{3t}+0.02e_{1}\sin{t}+0.01e_{2}\cos{3t}+0.01e_{13}\sin{3t}, \\ \theta_{221}(t) = \, &0.03e_{0}\sin{4t}+0.03e_{2}\cos{\sqrt{2}t}+0.02e_{12}\sin{5t}+0.03e_{23}\cos{3t}, \\ \theta_{222}(t) = \, &0.01e_{0}\sin{4t}+0.01e_{2}\cos{3t}+0.01e_{12}\sin{4t}+0.04e_{123}\cos{2t}, \\ q_{111}(t) = \, &0.06e_{0}\sin{5t}+0.04e_{1}\cos{6t}+0.03e_{12}\sin{\sqrt{3}t}+0.03e_{23}\sin^{2}{2t}, \\ q_{112}(t) = \, &0.05e_{0}\sin{2t}+0.02e_{1}\cos{\sqrt{2}t}+0.04e_{2}\sin{3t}+0.02e_{23}\cos{3t}, \\ q_{121}(t) = \, &0.02e_{0}\cos{4t}+0.02e_{1}\cos^{2}{3t}+0.05e_{2}\sin{4t}+0.03e_{23}\cos{4t}, \\ q_{122}(t) = \, &0.05e_{0}\sin{4t}+0.03e_{1}\cos{\sqrt{3}t}+0.03e_{12}\cos{2t}+0.04e_{23}\sin{3t}, \\ q_{211}(t) = \, &0.01e_{0}\cos{\sqrt{5}t}+0.02e_{1}\sin{3t}+0.04e_{12}\cos{\sqrt{3}t}+0.02e_{23}\sin{t}, \\ q_{212}(t) = \, &0.01e_{0}\cos{4t}+0.06e_{1}\cos{\sqrt{2}t}+0.03e_{12}\sin{5t}+0.04e_{23}\cos{2t}, \\ q_{221}(t) = \, &0.01e_{0}\cos{3t}+0.02e_{1}\sin{t}+0.02e_{2}\cos{\sqrt{3}t}+0.01e_{23}\sin{2t}, \\ q_{222}(t) = \, &0.01e_{0}\cos{9t}+0.05e_{1}\cos^{2}{3t}+0.03e_{12}\sin{t}+0.05e_{23}\cos{t}, \\ \nu_{111}(t) = \, &0.02e_{0}\sin{3t}+0.03e_{1}\cos{\sqrt{3}t}+0.01e_{2}\cos{5t}+0.01e_{12}\sin{4t}, \\ \nu_{112}(t) = \, &0.02e_{0}\cos{t}+0.03e_{1}\sin^{2}{5t}+0.02e_{12}\cos{3t}+0.03e_{23}\sin{2t}, \\ \nu_{121}(t) = \, &0.05e_{0}\sin{3t}+0.04e_{2}\sin{3t}+0.03e_{12}\cos{2t}+0.02e_{23}\sin{3t}, \\ \nu_{122}(t) = \, &0.03e_{0}\cos{2t}+0.04e_{2}\cos{\sqrt{3}t}+0.05e_{12}\sin{\sqrt{2}t}+0.04e_{23}\cos{\sqrt{3}t}, \\ \nu_{211}(t) = \, &0.03e_{0}\sin{7t}+0.04e_{1}\sin{3t}+0.02e_{2}\cos{\sqrt{3}t}+0.03e_{12}\sin{2t}, \\ \nu_{212}(t) = \, &0.05e_{0}\cos{5t}+0.02e_{2}\cos{4t}+0.03e_{12}\sin{7t}+0.04e_{23}\cos{t}, \\ \nu_{221}(t) = \, &0.02e_{0}\sin{3t}+0.03e_{2}\cos{2t}+0.04e_{12}\sin{t}+0.01e_{23}\cos{2t}, \\ \nu_{222}(t) = \, &0.03e_{0}\sin{6t}+0.03e_{2}\cos{t}+0.01e_{12}\sin{2t}+0.04e_{23}\cos{3t}, \\ I_1(t) = \, &0.32e_{0}(\sin{\sqrt{5}t}+\frac{1}{1+t^2})+0.5e_{1}\cos{\sqrt{2}t}+0.36e_{2}\sin{\sqrt{3}t}+0.23e_{3}\cos{\sqrt{3}t}\\ &+0.49e_{12}\sin{\sqrt{2}t}+0.25e_{13}\cos{\sqrt{5}t}+0.42e_{23}\cos{\sqrt{5}t}+0.15e_{123}(\sin{\sqrt{3}t}+e^{-|t|}), \\ I_2(t) = \, &0.25e_{0}\cos{\sqrt{3}t}+0.42e_{1}(\sin{\sqrt{2}t}+e^{-|t|})+0.28e_{2}\sin{\sqrt{3}t}+0.45e_{3}\cos{\sqrt{3}t}\\ &+0.32e_{12}(\cos{\sqrt{2}t}+\frac{1}{1+t^{2}})+0.46e_{13}\sin{\sqrt{3}t}+0.15e_{23}\sin{\sqrt{5}t}+0.26e_{123}\cos{\sqrt{3}t}, \\ \gamma_{ij}(t) = \, &0.07e_{0}\cos{2t}+0.03e_{1}\sin{2t}+0.04e_{2}(\sin{2t}+\frac{1}{1+t^{2}})+0.06e_{3}\cos{t}\\ &+0.08e_{12}(\sin{t}+e^{-|t|})+0.02e_{13}\cos{t}+0.05e_{23}\sin{2t}+0.04e_{123}\cos{t}, \\ \mu_{j}(t) = \, &0.2e_{0}\sin{\sqrt{2}t}+0.3e_{1}\sin{\sqrt{3}t}+0.4e_{2}\sin{\sqrt{2}t}+0.6e_{3}\cos{\sqrt{2}t}\\ &+0.8e_{12}\cos{\sqrt{3}t}+0.3e_{13}\cos{\sqrt{5}t}+0.4e_{23}\sin{\sqrt{3}t}+0.2e_{123}\cos{\sqrt{6}t}, \\ T_{ij}(t) = \, &0.006e_{0}\sin{t}+0.004e_{1}\sin{t}+0.003e_{2}(\sin{t}+\frac{1}{1+t^{2}})+0.001e_{3}\cos{t}\\ &+0.002e_{12}\sin{t}+0.003e_{13}(\cos{t}+e^{-|t|})+0.002e_{23}\cos{t}+0.005e_{123}\sin{t}, \\ S_{ij}(t) = \, &0.002e_{0}\cos{t}+0.003e_{1}\sin{t}+0.001e_{2}\cos{t}+0.001e_{3}(\sin{t}+\frac{1}{2+t^{2}})\\ &+0.003e_{12}(\sin{t}+e^{-|t+1|})+0.002e_{13}\cos{t}+0.001e_{23}\sin{t}+0.002e_{123}\cos{t}, \\ \sigma_{ij}(t) = \, &1-0.3\sin{t}, \quad \eta_{ij}(t) = 1-0.8\cos{3t}, \quad \delta_{ijk}(t) = 1-0.6\sin{2t}, \\ \tau_1(t) = \, &1-0.1\sin{t}, \quad \tau_2(t) = 1-0.3\cos{t} . \end{align*}$ $

Then, it is easy to see that conditions $ (A_1) $ and $ (A_2) $ are satisfied.

Moreover, take $ \vartheta_{1} = \vartheta_{2} = 1, p = 3, q = \frac{3}{2} $, then through simple calculations, we obtain

$\begin{align*} &\bar{a}_{1} = 0.014, \bar{a}_{2} = 0.012, \bar{b}_{1}^{\emptyset} = 10.05, \bar{b}_{1}^{c} = 0.26, \underline{b}_{1}^{\emptyset} = 9.95, \bar{b}_{2}^{\emptyset} = 10.05, \bar{b}_{2}^{c} = 0.27, \underline{b}_{2}^{\emptyset} = 9.95, \\ &\bar{c}_{11} = 0.03, \bar{c}_{12} = 0.03, \bar{c}_{21} = 0.03, \bar{c}_{22} = 0.04, \bar{u}_{11} = 0.04, \bar{u}_{12} = 0.04, \bar{u}_{21} = 0.03, \bar{u}_{22} = 0.04, \\ &\bar{\alpha}_{11} = 0.04, \bar{\alpha}_{12} = 0.03, \bar{\alpha}_{21} = 0.03, \bar{\alpha}_{22} = 0.04, \bar{\beta}_{11} = 0.03, \bar{\beta}_{12} = 0.03, \bar{\beta}_{21} = 0.04, \bar{\beta}_{22} = 0.04, \\ &\bar{\theta}_{111} = 0.03, \bar{\theta}_{112} = 0.04, \bar{\theta}_{121} = 0.04, \bar{\theta}_{122} = 0.04, \bar{\theta}_{211} = 0.03, \bar{\theta}_{212} = 0.04, \bar{\theta}_{221} = 0.03, \\ &\bar{\theta}_{222} = 0.04, \bar{q}_{111} = 0.06, \bar{q}_{112} = 0.05, \bar{q}_{121} = 0.05, \bar{q}_{122} = 0.05, \bar{q}_{211} = 0.04, \bar{q}_{212} = 0.06, \\ &\bar{q}_{221} = 0.02, \bar{q}_{222} = 0.05, \bar{\nu}_{111} = 0.03, \bar{\nu}_{112} = 0.03, \bar{\nu}_{121} = 0.05, \bar{\nu}_{122} = 0.05, \bar{\nu}_{211} = 0.04, \\ &\bar{\nu}_{212} = 0.05, \bar{\nu}_{221} = 0.04, \bar{\nu}_{222} = 0.04, \bar{\tau}_{i} = \bar{\tau}'_{i} = 0.1, \bar{\sigma}_{ij} = \bar{\sigma}'_{ij} = 0.3, \bar{\eta}_{ij} = \bar{\eta}'_{ij} = 0.8, \\ &\bar{\delta}_{ijk} = \bar{\delta}'_{ijk} = 0.6, L_{1}^{f} = L_{2}^{f} = \frac{1}{40}, L_{1}^{g} = L_{2}^{g} = \frac{1}{48}, M_{1}^{g} = M_{2}^{g} = \frac{1}{48} , \rho \approx0.054972 < 1, \\ &P\approx0.518973 < 1 . \end{align*}$ $

Hence, $ (A_{3}) $ and $ (A_{4}) $ are also satisfied. Consequently, in view of Theorem 3.2, we know that system (1.1) has a unique Besicovitch almost periodic solution that is globally exponentially stable (see Figures 1–4).

Figure 1. Curves of $ x_{1}^{0}(t), x_{2}^{0}(t), x_{1}^{1}(t) $, and $ x_{2}^{1}(t) $ of system (1.1) with two different initial values.

DownLoad: Full-Size Img PowerPoint

Figure 2. Curves of $ x_{1}^{2}(t), x_{2}^{2}(t), x_{1}^{3}(t) $, and $ x_{2}^{3}(t) $ of system (1.1) with two different initial values.

DownLoad: Full-Size Img PowerPoint

Figure 3. Curves of $ x_{1}^{12}(t), x_{2}^{12}(t), x_{1}^{13}(t) $, and $ x_{2}^{13}(t) $ of system (1.1) with two different initial values.

DownLoad: Full-Size Img PowerPoint

Figure 4. Curves of $ x_{1}^{23}(t), x_{2}^{23}(t), x_{1}^{123}(t) $, and $ x_{2}^{123}(t) $ of system (1.1) with two different initial values.

DownLoad: Full-Size Img PowerPoint

Remark 4.1. Even when the system considered in Example 4.1 degenerates into a real-valued system, there are no existing results to derive the results of Example 4.1.

5. Conclusions

This article introduces a new method to establish the existence and global exponential stability of Besicovitch almost periodic solutions for Clifford-valued high-order Hopfield fuzzy NNs with $ D $ operators. The methods and results of this article can be applied to study the generalized almost periodic and almost automorphic dynamics of high-order NNs.

Author contributions

Bing Li: Methodology, Conceptualization, Writing - review and editing; Yuan Ning: Writing - original draft, Visualization; Yongkun Li: Methodology, Conceptualization, Funding acquisition, Writing - review and editing. All authors have read and approved the final version of the manuscript for publication.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was supported by the National Natural Science Foundation of China, grant number 12261098.

Conflict of interest

The authors declare no conflict of interest.

References

[1]	J.-P. Aubin, "Viability Theory," Systems & Control: Foundations & Applications, Birkhäuser, Boston, Inc., Boston, MA, 1991.
[2]	J.-P. Aubin, A. M. Bayen and P. Saint-Pierre, Dirichlet problems for some hamilton-jacobi equations with inequality constraints, SIAM Journal on Control and Optimization, 47 (2008), 2348-2380. doi: 10.1137/060659569
[3]	J.-P. Aubin, A. M. Bayen and P. Saint-Pierre, "Viability Theory: New Directions," New directions. Second edition. Springer, Heidelberg, 2011. doi: 10.1007/978-3-642-16684-6
[4]	J.-P. Aubin and A. Cellina, "Differential Inclusions," Set-valued maps and viability theory. Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences], 264, Springer-Verlag, Berlin, 1984. doi: 10.1007/978-3-642-69512-4
[5]	E. N. Barron and R. Jensen, Semicontinuous viscosity solutions for Hamilton-Jacobi equations with convex Hamiltonians, Comm. Partial Differential Equations, 15 (1990), 1713-1742. doi: 10.1080/03605309908820745
[6]	A. M. Bayen and C. G. Claudel, Solutions to switched Hamilton-Jacobi equations and conservation laws using hybrid components, Hybrid systems: Computation and control, 101–115, Lecture Notes in Computer Science, 4981, Springer, Berlin, 2008. doi: 10.1007/978-3-540-78929-1_8
[7]	C. Canudas de Wit, Best-effort highway traffic congestion control via variable speed limits, In "Decision and Control and European Control Conference (CDC-ECC), 2011 50th IEEE Conference Proceedings," (2001), 5959-5964.
[8]	R. C. Carlson, I. Papamichail, M. Papageorgiou and A. Messmer, Optimal motorway traffic flow control involving variable speed limits and ramp metering, Transportation Science, 44 (2010), 238-253. doi: 10.1287/trsc.1090.0314
[9]	C. G. Claudel and A. M. Bayen, Convex formulations of data assimilation problems for a class of Hamilton-Jacobi equations, SIAM Journal on Control and Optimization, 49 (2011), 383-402. doi: 10.1137/090778754
[10]	C. G. Claudel and A. M. Bayen, Lax-Hopf based incorporation of internal boundary conditions into Hamilton-Jacobi equation. Part I: Theory, IEEE Transactions on Automatic Control, 55 (2010), 1142-1157. doi: 10.1109/TAC.2010.2041976
[11]	C. G. Claudel and A. M. Bayen, Lax-Hopf based incorporation of internal boundary conditions into Hamilton-Jacobi equation. Part II: Computational methods, IEEE Transactions on Automatic Control, 55 (2010), 1158-1174. doi: 10.1109/TAC.2010.2045439
[12]	C. Daganzo, A variational formulation of kinematic waves: Basic theory and complex boundary conditions, Transporation Research B, 39 (2005), 187-196. doi: 10.1016/j.trb.2004.04.003
[13]	C. Daganzo, On the variational theory of traffic flow: Well-posedness, duality and applications, Networks and Heterogeneous Media, 1 (2006), 601-619. doi: 10.3934/nhm.2006.1.601
[14]	L. C. Edie, Car following and steady state theory for non-congested traffic, Operations Research, 9 (1961), 66-76. doi: 10.1287/opre.9.1.66
[15]	H. Frankowska, Lower semicontinuous solutions of Hamilton-Jacobi-Bellman equations, SIAM J. Control Optim., 31 (1993), 257-272. doi: 10.1137/0331016
[16]	H. Greenberg, An analysis of traffic flow, Operations Research, 7 (1959), 79-85. doi: 10.1287/opre.7.1.79
[17]	B. D. Greenshields, A study of traffic capacity, HRB Proc., 14 (1934), 448-481.
[18]	M. J. Lighthill and G. B. Whitham, On kinematic waves: II. A theory of traffic flow on long crowded roads, Proc. Royal Society, Ser. A, 229 (1955), 317-345. doi: 10.1098/rspa.1955.0089
[19]	G. Newell, Nonlinear effects in the dynamics of car following, Operations Research, 9 (1961), 209-229. doi: 10.1287/opre.9.2.209
[20]	M. Papageorgiou, E. Kosmatopoulos and I. Papamichail, Effects of variable speed limits on motorway traffic flow, Transportation Research Record: Journal of the Transportation Research Board, 2047 (2008), 37-48. doi: 10.3141/2047-05
[21]	P. I. Richards, Shock waves on the highway, Operations Research, 4 (1956), 42-51. doi: 10.1287/opre.4.1.42
[22]	P. Saint-Pierre, Approximation of the viability kernel, Applied Mathematics and Optimisation, 29 (1994), 187-209. doi: 10.1007/BF01204182
[23]	R. T. Underwood, Speed, volume and density relationships, quality and theory of traffic flow, in "Yale Bureau of Highway Traffic," 1961, 141-88.

Reader Comments

Your name:*

Email:*
© 2013 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

Networks and Heterogeneous Media

1.3 1.9

Metrics

Article views(3364) PDF downloads(97) Cited by(3)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

Networks and Heterogeneous Media

Viability approach to Hamilton-Jacobi-Moskowitz problem involving variable regulation parameters

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Numerical examples

5. Conclusions

Author contributions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Catalog

Networks and Heterogeneous Media

Viability approach to Hamilton-Jacobi-Moskowitz problem involving variable regulation parameters

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Numerical examples

5. Conclusions

Author contributions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog