Existence and Stability of a Nonlinear Distributed Delayed Periodic AG-Ecosystem with Competition on Time Scales

Abstract

The Ayala-Gilpin (AG) kinetics system is one of the famous mathematical models of ecosystem. This model has been widely concerned and studied since it was proposed. This paper stresses on a nonlinear distributed delayed periodic AG-ecosystem with competition on time scales. In the sense of time scale, our model unifies and generalizes the discrete and continuous cases. Firstly, with the aid of the auxiliary function having only two zeros in the real number field, we apply inequality technique and coincidence degree theory to obtain some sufficient criteria which ensure that this model has periodic solutions on time scales. Meanwhile, the global asymptotic stability of the periodic solution is founded by employing stability theory in the sense of Lyapunov. Eventually, we provide an illustrative example and conduct numerical simulation by means of MATLAB tools.

1. Introduction

This paper mainly deals with the following nonlinear competitive periodic Ayala-Gilpin ecosystem with distributed lags on time scale

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{\mathcal{U}}^{\Delta }\left(\varsigma \right)=r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right){\left[e^{\mathcal{U}\left(\varsigma \right)}\right]}^{{\theta }_1}-c_{12}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}\Delta s-{\varphi }_1\left(\varsigma \right)e^{-\mathcal{U}\left(\varsigma \right)},\varsigma \in \mathbb{T},\hfill \\ {\mathcal{V}}^{\Delta }\left(\varsigma \right)=r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right){\left[e^{\mathcal{V}\left(\varsigma \right)}\right]}^{{\theta }_2}-c_{21}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}\Delta s-{\varphi }_2\left(\varsigma \right)e^{-\mathcal{V}\left(\varsigma \right)},\varsigma \in \mathbb{T},\hfill \\ \mathcal{U}\left(\varsigma \right)={\phi }_1\left(\varsigma \right),\mathcal{V}\left(\varsigma \right)={\phi }_2\left(\varsigma \right),\varsigma \in (-\xi ,0]\cap \mathbb{T},\hfill \end{array}\right.\end{array}$$

(1)

where $\mathbb{T}$ is a time scale, $\Delta$ is the delta (or Hilger) derivative on $\mathbb{T}$, $\mathcal{U}\left(\varsigma \right)$ and $\mathcal{V}\left(\varsigma \right)$ are the quantity densities of two vying species at moment $\varsigma$, $r_i\left(\varsigma \right)(i=1,2)$ stands for the intrinsic birth rate, $c_{11}\left(\varsigma \right)$ and $c_{22}\left(\varsigma \right)$ are the intraspecific competition rates, $c_{12}\left(\varsigma \right)$ and $c_{21}\left(\varsigma \right)$ are the interspecific competition rates, the kernel function of distributed-lags is given by $k_i\left(\varsigma \right)(i=1,2)$, ${\xi }_i\left(\varsigma \right)>0(i=1,2)$ is the distributed-lag function, ${\varphi }_i\left(\varsigma \right)(i=1,2)$ is the manual control term, the constants ${\theta }_i>0(i=1,2)$ measures the nonlinear interferences within species, ${\phi }_i\left(\varsigma \right)(i=1,2)$ is the initial function, $\xi =max\{{sup}_{\varsigma \in \mathbb{T}}{\xi }_1\left(\varsigma \right),{sup}_{\varsigma \in \mathbb{T}}{\xi }_2\left(\varsigma \right)\}$.

The proposal of the model is related to the famous experimental study of Drosophila competition. Combined with the experimental study of Drosophila competition, Ayala, Gilpin and Eherenfeld [1] put forward the following nonlinear dynamic model in 1973.

$$\begin{array}{c}\hfill \left\{\begin{array}{c}\frac{d\mathcal{U}\left(\varsigma \right)}{d\varsigma }=r_1\mathcal{U}\left(\varsigma \right)\left[1-{\left(\frac{\mathcal{U}\left(\varsigma \right)}{K_1}\right)}^{{\theta }_1}-c_{12}\frac{\mathcal{V}\left(\varsigma \right)}{K_2}\right],\hfill \\ \frac{d\mathcal{V}\left(\varsigma \right)}{d\varsigma }=r_2\mathcal{V}\left(\varsigma \right)\left[1-{\left(\frac{\mathcal{V}\left(\varsigma \right)}{K_2}\right)}^{{\theta }_2}-c_{21}\frac{\mathcal{U}\left(\varsigma \right)}{K_1}\right],\hfill \end{array}\right.\end{array}$$

(2)

where $r_i>0(i=1,2)$ represents inherent growth ratios. $K_i>0(i=1,2)$ stands for the maximum capacity of environment to species. The constant ${\theta }_i(i=1,2)$ measures the nonlinear interferences within species. $c_{12}>0$ and $c_{21}>0$ mean the vying rate among species.

The nonlinear interferences within species are measured by constants ${\theta }_1>0$ and ${\theta }_2>0$. $c_{12}>0$ and $c_{21}>0$ are the measures of competition between species.

When ${\theta }_1={\theta }_2=1$, (2) becomes the below Lotka-Volterra competitive model

$$\begin{array}{c}\hfill \left\{\begin{array}{c}\frac{d\mathcal{U}\left(\varsigma \right)}{d\varsigma }=r_1\mathcal{U}\left(\varsigma \right)\left[1-\frac{\mathcal{U}\left(\varsigma \right)}{K_1}-c_{12}\frac{\mathcal{V}\left(\varsigma \right)}{K_2}\right],\hfill \\ \frac{d\mathcal{V}\left(\varsigma \right)}{d\varsigma }=r_2\mathcal{V}\left(\varsigma \right)\left[1-\frac{\mathcal{V}\left(\varsigma \right)}{K_2}-c_{21}\frac{\mathcal{U}\left(\varsigma \right)}{K_1}\right].\hfill \end{array}\right.\end{array}$$

(3)

Therefore, the Ayala-Gilpin model is a generalization of Lotka-Volterra model. The parameters ${\theta }_1$ and ${\theta }_2$ can take any positive real number. Therefore, Ayala-Gilpin model has been widely favored since it was proposed. The kinetics properties of Ayala-Gilpin model have been extensively explored. In [2,3,4,5,6,7], the authors dealt with the persistence, extinction and attraction of AG-system. Amdouni et al. [8] investigated the existence and global exponential stability of pseudo almost periodic solutions for a generalized competitive AG-model. Korobenko et al. [9] handled the evolutionary stability of a diffusion AG-model. Zhao [10,11] studied the multiplicity of almost periodic solution and local exponential stability for two AG-systems with harvest term. If the time-lag, impulse and random effects are considered in the AG-system, many excellent results have been achieved (see [12,13,14,15,16,17,18,19,20]). with the exception of the classical AG-system, some extended AG-systems have also been studied (see [4,6,8,9,14,16,20,21,22]).

Moreover, the model (1) is more suitable for the actual situation of ecosystem. For example, the predation process is not instantaneous, but takes a period of time to complete. At the same time, the number of predators can not be increased immediately after predators prey. Therefore, time lag is common in the whole process of predation and transformation. As we all know, the environment of the ecosystem often presents certain periodic changes over time. For the harmonious coexistence of people and ecosystem, sometimes it is necessary to protect and intervene the ecosystem manually.

From the perspective of mathematical theory, it is also of great value to study model (1). According to the definition of time scale, the model (1) contains difference case and differential case. Indeed, when $\mathbb{T}={\mathbb{N}}^{+}$, the model (1) becomes the following difference equation

$$\begin{array}{c}\hfill \left\{\begin{array}{c}\mathcal{U}(\varsigma +1)-\mathcal{U}\left(\varsigma \right)=r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right){\left[e^{\mathcal{U}\left(\varsigma \right)}\right]}^{{\theta }_1}-c_{12}\left(\varsigma \right)\hfill \\ \times {\displaystyle \sum _{s=-{\xi }_1\left(\varsigma \right)}^0}{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}-{\varphi }_1\left(\varsigma \right)e^{\mathcal{U}\left(\varsigma \right)},\varsigma \in {\mathbb{N}}^{+},\hfill \\ \mathcal{V}(\varsigma +1)-\mathcal{V}\left(\varsigma \right)=r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right){\left[e^{\mathcal{V}\left(\varsigma \right)}\right]}^{{\theta }_2}-c_{21}\left(\varsigma \right)\hfill \\ \times {\displaystyle \sum _{s=-{\xi }_2\left(\varsigma \right)}^0}{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}-{\varphi }_2\left(\varsigma \right)e^{\mathcal{V}\left(\varsigma \right)}],\varsigma \in {\mathbb{N}}^{+},\hfill \\ \mathcal{U}\left(\varsigma \right)={\phi }_1\left(\varsigma \right),\mathcal{V}={\phi }_2\left(\varsigma \right),\varsigma \in {\mathbb{N}}^{+}.\hfill \end{array}\right.\end{array}$$

(4)

When $\mathbb{T}=\mathbb{R}$, let $x\left(\varsigma \right)=e^{\mathcal{U}\left(\varsigma \right)},y\left(\varsigma \right)=e^{\mathcal{V}\left(\varsigma \right)}$, then the model (1) becomes the following differential equation

$$\begin{array}{c}\hfill \left\{\begin{array}{c}\frac{dx\left(\varsigma \right)}{d\varsigma }=x\left(\varsigma \right)[r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right){\left[x\left(\varsigma \right)\right]}^{{\theta }_1}-c_{12}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)y(\varsigma +s)ds]-{\varphi }_1\left(\varsigma \right),\varsigma \in \mathbb{R},\hfill \\ \frac{dy\left(\varsigma \right)}{d\varsigma }=y\left(\varsigma \right)[r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right){\left[y\left(\varsigma \right)\right]}^{{\theta }_2}-c_{21}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)x(\varsigma +s)ds]-{\varphi }_2\left(\varsigma \right),\varsigma \in \mathbb{R},\hfill \\ x\left(\varsigma \right)=e^{{\phi }_1\left(\varsigma \right)},y\left(\varsigma \right)=e^{{\phi }_2\left(\varsigma \right)},\varsigma \in \mathbb{R}.\hfill \end{array}\right.\end{array}$$

(5)

In addition, a complex number set like $\mathbb{T}={\bigcup }_{n=-\infty }^{+\infty }[n+\frac{1}{4},n+\frac{1}{2}]\cup \left\{n\right\}$ is also a time scale. In 1988, Hilger first raised the time scale theory in his Ph.D. thesis [23], aiming at unifying difference and differential. For further study of time scale theory, please refer to the monographs [24,25] To the best my knowledge, no one has studied the periodic solution and stability of AG-system on time scales. Consequently, it is worthwhile to study the periodic dynamic behavior of (1).

The highlights of our work and the differences from previous published works are mainly manifested in two aspects. (i) We study the dynamic properties of AG-system in the sense of time scales, which can unify the differential and difference forms of AG-systems within the same framework. However, most previous studies on AG-systems have focused on the case of continuous differential equations (see some of them [6,7,8,9,10,11,12,13,14,15,16,22]). Previous papers dealing with pure difference AG-systems are rare. some of their papers involved discrete AG-systems containing impulsive terms (see [2,17,21]). (ii) Our model (1) is the closest to classical model (2). The model (1) only adds control harvest terms and distribution lags based on classical model (2). Compared to model (1), some previous papers on generalized AG-systems have made many modifications to classical model (2). For example, Wang et al. [4] added random terms and multiple species to their model. In [6,14,16], the authors not only added multiple species to their models, but also added nonlinear measures ${\theta }_{ij}(i\ne j)$ to the inter species competition term. In [20,21], the authors added Lévy and Markovian jumps to their models.

The surplus framework of the manuscript is arranged as follows. We mainly introduce the basic concepts and results of time scales, important assumptions and necessary propositions in Section 2. Section 3 focuses on finding the sufficient conditions for the existence of periodic solutions to system (1). In Section 4, the Lyapunov functional is applied to discuss the global asymptotic stability of system (1). Section 5 provides an example and numerical simulation to verify our main results. In the last Section 6, we made a brief conclusion and outlook.

2. Preliminaries

This section first concisely reviews the elementary knowledge of calculus on time scale. The following statements are taken from Refs. [24,25].

$\varphi \ne \mathbb{T}\subset \mathbb{R}$ is called a time scale if $\mathbb{T}$ is closed. Define some operators as follows:

$$\sigma \left(\varsigma \right)=inf\{\iota \in \mathbb{T}:\iota >\varsigma \},\rho \left(\varsigma \right)=sup\{\iota \in \mathbb{T}:\iota <\varsigma \},\mu \left(\varsigma \right)=\sigma \left(\varsigma \right)-\varsigma ,\forall \varsigma \in \mathbb{T},$$

then, we call that $\sigma$ is jumping forward, $\rho$ is jumping backward, and $\mu$ is the graininess.

For $\varsigma \in \mathbb{T}$, if $\varsigma >inf\mathbb{T}$ and $\rho \left(\varsigma \right)=\varsigma$ ($\varsigma <sup\mathbb{T}$ and $\sigma \left(\varsigma \right)=\varsigma$), then we call that $\varsigma$ is left-dense (right-dense). If $\rho \left(\varsigma \right)<\varsigma$ ($\sigma \left(\varsigma \right)>\varsigma$), then we call that $\varsigma$ is left-scattered (right-scattered). In addition, we have

$$\begin{array}{c}\hfill {\mathbb{T}}^k=\left\{\begin{array}{cc}\mathbb{T}\backslash \left\{M\right\},\hfill & \mathbb{T}\mathrm{achieves}\mathrm{a}\mathrm{left}\text{-}\mathrm{scattered}\mathrm{maximum}M,\hfill \\ \mathbb{T},\hfill & \mathrm{otherwise},\hfill \end{array}\right.\end{array}$$

and

$$\begin{array}{c}\hfill {\mathbb{T}}^k=\left\{\begin{array}{cc}\mathbb{T}\backslash \left\{m\right\},\hfill & \mathbb{T}\mathrm{reachs}\mathrm{a}\mathrm{right}\text{-}\mathrm{scattered}\mathrm{minimum}m,\hfill \\ \mathbb{T},\hfill & \mathrm{otherwise}.\hfill \end{array}\right.\end{array}$$

For $ɷ>0$, we call that a time scale $\mathbb{T}$ is $ɷ$-periodic, provided that $\forall \varsigma \in \mathbb{T}\Rightarrow \varsigma +ɷ\in \mathbb{T}$. Obviously, if $\mathbb{T}$ is $ɷ$-periodic, then $\mathbb{T}$ is unbounded above.

Definition 1.

We call that $u:\mathbb{T}\to \mathbb{R}$ is regulated iff ${lim}_{t\to {\varsigma }^{+}}u\left(t\right)$ and ${lim}_{t\to {\varsigma }^{-}}u\left(t\right)$ all exist (finite), $\forall \varsigma \in \mathbb{T}$.

Definition 2.

We call that $u:\mathbb{T}\to \mathbb{R}$ is $rd$-continuous iff, for all right-dense point. $\varsigma \in \mathbb{T}$, we have ${lim}_{t\to {\varsigma }^{+}}u\left(t\right)=u\left(\varsigma \right)$, and for all left-dense point $\iota \in \mathbb{T}$, we have that ${lim}_{t\to {\iota }^{-}}u\left(t\right)=u\left({\iota }^{-}\right)$ exists (finite). We denote the collection of all $rd$-continuous functions $u:\mathbb{T}\to \mathbb{R}$ as $C_{rd}(\mathbb{T},\mathbb{R})$.

Definition 3.

Let $u:\mathbb{T}\to \mathbb{R}$ and $\varsigma \in {\mathbb{T}}^k$. We call that a number $u^{\Delta }\left(\varsigma \right)$ (if exists) is the Δ-derivative of u at ς iff, for a given $\forall ϵ>0$, there has $\delta >0$ such that

$$|[u\left(\sigma \left(\varsigma \right)\right)-u\left(\iota \right)]-u^{\Delta }\left(\varsigma \right)[\sigma \left(\varsigma \right)-\iota ]|<ϵ|\sigma \left(\varsigma \right)-\iota |,\forall \iota \in (\varsigma -\delta ,\varsigma +\delta )\cap \mathbb{T}.$$

The collection of first-order Δ-differentiable functions is denoted by

$$C_{rd}^1(\mathbb{T},\mathbb{R})=\{u:\mathbb{T}\to \mathbb{R}:u^{\Delta }\left(\varsigma \right)\mathrm{exists},\mathrm{and}{u}^{\Delta }\left(\varsigma \right)\in C_{rd}(\mathbb{T},\mathbb{R}),\forall \varsigma \in \mathbb{R}\}.$$

According to the above definitions, one easily knows that u is $\Delta$-differentiable $\Rightarrow u$ is continuous $\Rightarrow u$ is $rd$-continuous $\Rightarrow u$ is regulated.

Lemma 1.

If u is regulated, then there is a Δ-differentiable function U with differentiable region D that contents

$$U^{\Delta }\left(\varsigma \right)=u\left(\varsigma \right),\forall \varsigma \in D.$$

Definition 4.

Assume that $u:\mathbb{T}\to \mathbb{R}$ is regulated, then we have the following concepts.

$\left(1\right)$

We call that any function U as in Lemma 1 is a Δ-antiderivative of u.

$\left(2\right)$

We define the Δ-indefinite integral of u as

$$\int u(\varsigma )\Delta \varsigma =U(\varsigma )+c,$$

where an arbitrary constant c is a Δ-integral constant.

$\left(3\right)$

We define the Δ-definite integral as

$${\int }_{\alpha }^{\beta }u\left(\varsigma \right)\Delta \varsigma =U\left(\beta \right)-U\left(\alpha \right),\forall \alpha ,\beta \in \mathbb{T}.$$

Lemma 2.

Let $\alpha ,\beta \in \mathbb{T}$, ${\lambda }_1,{\lambda }_2\in \mathbb{R}$, $u,v\in C_{rd}(\mathbb{T},\mathbb{R})$, then we have the followings

$\left(1\right)$

${\int }_{\alpha }^{\beta }[{\lambda }_{1}u\left(\varsigma \right)+{\lambda }_{2}v\left(\varsigma \right)]\Delta \varsigma ={\lambda }_1{\int }_{\alpha }^{\beta }u\left(\varsigma \right)\Delta \varsigma +{\lambda }_2{\int }_{\alpha }^{\beta }v\left(\varsigma \right)\Delta \varsigma$.

$\left(2\right)$

$\forall \alpha \le \varsigma <\beta$, $u\left(\varsigma \right)\ge 0\Rightarrow {\int }_{\alpha }^{\beta }u\left(\varsigma \right)\Delta \varsigma \ge 0$.

$\left(3\right)$

$\forall \varsigma \in \{\varsigma \in \mathbb{T}:\alpha \le \varsigma <\beta \}$, $\left|u\left(\varsigma \right)\right|\le v\left(\varsigma \right)\Rightarrow |{\int }_{\alpha }^{\beta }u\left(\varsigma \right)\Delta \varsigma |\le {\int }_{\alpha }^{\beta }v\left(\varsigma \right)\Delta \varsigma$.

Lemma 3

([26]). For two Banach spaces $\mathfrak{X}$, $\mathfrak{Z}$, and nonempty bounded open subset $\Omega \subset \mathfrak{X}$, define some operators $\mathcal{L}:\mathfrak{X}\to \mathfrak{Z}$, $\mathcal{N}:\mathfrak{X}\times [0,1]\to \mathfrak{Z}$, $\mathcal{Q}:\mathfrak{Z}\to \mathfrak{Z}$ and $\mathcal{J}:\mathfrak{Z}\to \mathfrak{Z}$. Assume that $\mathcal{L}$ is zero index Fredholm type, $\mathcal{N}$ is $\mathcal{L}$-compact, $\mathcal{Q}$ is projected, and $\mathcal{J}$ is homotopy. Further suppose that

$\left(\mathrm{a}\right)$

For all $0<\eta <1$, if x is a solution of $\mathcal{L}x=\eta \mathcal{N}(x,\eta )$, then $x\notin \partial \Omega \cap \mathrm{Dom}\left(\mathcal{L}\right)$;

$\left(\mathrm{b}\right)$

If $x\in \partial \Omega \cap \mathrm{Ker}\left(\mathcal{L}\right)$, then $\mathcal{Q}\mathcal{N}(x,0)x\ne 0$;

$\left(\mathrm{c}\right)$

$deg\left(\mathcal{J}\mathcal{Q}\mathcal{N}\right(x,0),\Omega \cap \mathrm{Ker}(\mathcal{L}),0)\ne 0$.

Then there has at least a $x^{\ast }\in \overline{\Omega }\cap \mathrm{Dom}\left(\mathcal{L}\right)$ meeting with $\mathcal{L}{x}^{\ast }=\mathcal{N}(x^{\ast },1)$.

Lemma 4

([11]). Let $\alpha ,\beta ,\gamma ,\vartheta >0$, $F\left(x\right)=\alpha e^{(1+\vartheta )x}-\beta e^x+\gamma$. Assume that $\vartheta {\alpha }^{-\frac{1}{\vartheta }}{\left(\frac{\beta }{1+\vartheta }\right)}^{\frac{1+\vartheta }{\vartheta }}>\gamma$. Then we have the followings

$\left(1\right)$

There exists a unique $x_0=\frac{1}{\vartheta }ln\left[\frac{\beta }{\alpha (1+\vartheta )}\right]\in \mathbb{R}$ such that

$$F_{min}=F\left(x_0\right)=-\vartheta {\alpha }^{-\frac{1}{\vartheta }}{\left(\frac{\beta }{1+\vartheta }\right)}^{\frac{1+\vartheta }{\vartheta }}+\gamma <0,$$

and for all $x_1<x_2$, $F\left(x_1\right)>F\left(x_2\right)$ when $x_1,x_2\in (-\infty ,x_0]$, $F\left(x_1\right)<F\left(x_2\right)$ when $x_1,x_2\in [x_0,+\infty )$.

$\left(2\right)$

There are only two $x_1^{\ast },x_2^{\ast }\in \mathbb{R}$ with $x_1^{\ast }<x_2^{\ast }$ such that $F\left(x_1^{\ast }\right)=F\left(x_2^{\ast }\right)=0$.

Now let’s introduce some symbols.

$$min\{[0,+\infty )\cap \mathbb{T}\}\triangleq \kappa ,I_{ɷ}=[\kappa ,\kappa +ɷ]\cap \mathbb{T},\underset{\varsigma \in I_{ɷ}}{sup}w\left(\varsigma \right)\triangleq \overline{w},$$

$$\underset{\varsigma \in I_{ɷ}}{inf}w\left(\varsigma \right)\triangleq \underline{w},\frac{1}{ɷ}{\int }_{I_{ɷ}}w\left(s\right)\Delta s=\frac{1}{ɷ}{\int }_{\kappa }^{\kappa +ɷ}w\left(s\right)\Delta s\triangleq \widehat{w},$$

where $w(\varsigma +ɷ)=w\left(\varsigma \right)$, $\forall w\in C_{rd}(\mathbb{T},\mathbb{R})$. The whole paper needs the following basic assumptions.

$\left({\mathrm{H}}_1\right)$

Assume that $0<r_1\left(\varsigma \right)$, $r_2\left(\varsigma \right)$, ${\xi }_1\left(\varsigma \right)$, ${\xi }_2\left(\varsigma \right)$, ${\varphi }_1\left(\varsigma \right)$, ${\varphi }_2\left(\varsigma \right)$, $c_{11}\left(\varsigma \right)$, $c_{22}\left(\varsigma \right)$, $c_{12}\left(\varsigma \right)$, $c_{21}\left(\varsigma \right)$, ${\phi }_1\left(\varsigma \right)$, ${\phi }_2\left(\varsigma \right)$, $k_1\left(\varsigma \right)$, $k_2\left(\varsigma \right)\in C_{rd}(\mathbb{T},\mathbb{R})$ are all $ɷ$-periodic, and satisfy ${\int }_{-\xi }^0{k}_1\left(s\right)\Delta s<\infty$, ${\int }_{-\xi }^0{k}_1\left(s\right)\Delta s<\infty$, where $\xi =max\{\overline{{\xi }_1},\overline{{\xi }_2}\}$.

3. Existence of Periodic Solution on Time Scales

In the portion, we shall put to use Lemma 3 to argue that system (1) has a periodic solution. Set $\mathfrak{X}=\mathfrak{Z}=W_1\oplus W_2$, where $W_1=\{w\left(\varsigma \right)\equiv {(C_1,C_2)}^T\in {\mathbb{R}}^2\}$,

$$\begin{array}{c}\hfill W_2=\left\{w\left(\varsigma \right)={(w_1\left(\varsigma \right),w_2\left(\varsigma \right))}^T:w_j\left(\varsigma \right)\in C_{rd}(\mathbb{T},\mathbb{R}),w_j(\varsigma +ɷ)=w_j\left(\varsigma \right),j=1,2\right\},\end{array}$$

equipped with the norm

$$\begin{array}{c}\hfill \parallel w\parallel =\underset{1\le j\le 2}{max}\underset{\varsigma \in I_{ɷ}}{sup}\left|w_j\left(\varsigma \right)\right|,\forall w={(w_1,w_2)}^T\in \mathfrak{X}=\mathfrak{Z}.\end{array}$$

(6)

In the manner of Ref. [27], it is easy to prove Lemmas 5–8. So, we omit their proofs.

Lemma 5.

$\mathfrak{X}=\mathfrak{Z}$ is the Banach space with the norm $\left|\right|·\left|\right|$ defined as (6).

Lemma 6.

$\mathcal{L}:\mathfrak{X}\to \mathfrak{Z}$ defined by

$$\mathcal{L}w\left(\varsigma \right)=w^{\Delta }\left(\varsigma \right)={(w_1^{\Delta }\left(\varsigma \right),w_2^{\Delta }\left(\varsigma \right))}^T,\forall w\left(\varsigma \right)={(w_1\left(\varsigma \right),w_2\left(\varsigma \right))}^T\in \mathfrak{X},$$

then $\mathcal{L}$ is zero index Fredholm type.

Lemma 7.

For all $w={(\mathcal{U},\mathcal{V})}^T\in \mathfrak{X}=\mathfrak{Z}$, $\mathcal{N}(w,\eta ):\mathfrak{X}\times [0,1]\to \mathfrak{Z}$, $\mathcal{P}:\mathfrak{X}\to \mathfrak{Z}$ and $\mathcal{Q}:$$\mathfrak{Z}$$\to \mathfrak{Z}$ are defined by

$$\begin{array}{c}\hfill \mathcal{N}(w,\eta )=\left(\begin{array}{c}{r}_1\left(\varsigma \right)-c_{11}\left(\varsigma \right)e^{{\theta }_1\mathcal{U}\left(\varsigma \right)}-\eta c_{12}\left(\varsigma \right){\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}ds-{\varphi }_1\left(\varsigma \right)e^{-\mathcal{U}\left(\varsigma \right)}\\ r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right)e^{{\theta }_2\mathcal{V}\left(\varsigma \right)}-\eta c_{21}\left(\varsigma \right){\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}ds-{\varphi }_2\left(\varsigma \right)e^{-\mathcal{V}\left(\varsigma \right)}\end{array}\right),\end{array}$$

$$\mathcal{P}w=\mathcal{Q}w={(\widehat{\mathcal{U}},\widehat{\mathcal{V}})}^T={\left(\frac{1}{ɷ}{\int }_{\kappa }^{\kappa +ɷ}\mathcal{U}\left(\varsigma \right)\Delta \varsigma ,\frac{1}{ɷ}{\int }_{\kappa }^{\kappa +ɷ}\mathcal{V}\left(\varsigma \right)\Delta \varsigma \right)}^T,$$

$$\mathrm{Ker}\left(\mathcal{L}\right)=\{w=(\mathcal{U},\mathcal{V})\in \mathfrak{X}:(\mathcal{U},\mathcal{V})=(C_1,C_2),\varsigma \in \mathbb{T}\},$$

$$\mathrm{Im}\left(\mathcal{L}\right)=\{w=(\mathcal{U},\mathcal{V})\in \mathfrak{Z}:(\widehat{\mathcal{U}},\widehat{\mathcal{V}})=(0,0)\},$$

$$\left({\mathcal{L}}^{-1}{|}_{\mathcal{P}}\right)\left(w\left(\varsigma \right)\right)=\left(\begin{array}{c}{\int }_{\kappa }^{\varsigma }\mathcal{U}\left(s\right)\Delta s-\frac{1}{ɷ}{\int }_{\kappa }^{\kappa +ɷ}{\int }_{\kappa }^{\varsigma }\mathcal{U}\left(s\right)\Delta s\Delta \varsigma \\ {\int }_{\kappa }^{\varsigma }\mathcal{V}\left(s\right)\Delta s-\frac{1}{ɷ}{\int }_{\kappa }^{\kappa +ɷ}{\int }_{\kappa }^{\varsigma }\mathcal{V}\left(s\right)\Delta s\Delta \varsigma \end{array}\right).$$

Then $\mathcal{N}(w,\eta )$ is $\mathcal{L}$-compact on $\overline{\Omega }\times [0,1]$.

Lemma 8.

Let $\mathbb{T}$ be an ϖ-periodic time scale. Suppose $\psi :\mathbb{T}\to \mathbb{R}$ be an ϖ-periodic function which is rd-continuous, then

$$\begin{array}{c}\hfill 0\le \underset{s\in I_{ɷ}}{sup}\psi \left(s\right)-\underset{s\in I_{ɷ}}{inf}\psi \left(s\right)\le \frac{1}{2}{\int }_{\kappa }^{\kappa +ɷ}\left|{\psi }^{\Delta }\left(s\right)\right|\Delta s.\end{array}$$

$\left({\mathrm{H}}_2\right)$

Suppose that all inequalities hold as follows:

$${\theta }_1{\left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)}^{-\frac{1}{{\theta }_1}}{\left(\frac{\overline{r_1}}{1+{\theta }_1}\right)}^{\frac{1+{\theta }_1}{{\theta }_1}}>\underline{{\varphi }_1},{\theta }_2{\left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)}^{-\frac{1}{{\theta }_2}}{\left(\frac{\overline{r_2}}{1+{\theta }_2}\right)}^{\frac{1+{\theta }_2}{{\theta }_2}}>\underline{{\varphi }_2}.$$

Theorem 1.

Assume that $\left({\mathrm{H}}_1\right)$ and $\left({\mathrm{H}}_2\right)$ are true, then there has at least an ϖ-periodic solution ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ satisfies model (1) on periodic time scale $\mathbb{T}$ such that $L^{-}-ɷ\overline{r_1}<\tilde{\mathcal{U}}\left(\varsigma \right)<L^{+}$ and $M^{-}-ɷ\overline{r_2}<\tilde{\mathcal{V}}\left(\varsigma \right)<M^{+}$, where $L^{\pm }$ and $M^{\pm }$ can be solved by the below equations

$$\begin{array}{cc}\hfill f\left(L^{\pm }\right)=& \left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)e^{(1+{\theta }_1)L^{\pm }}-\overline{r_1}{e}^{L^{\pm }}+\underline{{\varphi }_1}=0,\hfill \\ \hfill g\left(M^{\pm }\right)=& \left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)e^{(1+{\theta }_2)M^{\pm }}-\overline{r_2}{e}^{M^{\pm }}+\underline{{\varphi }_2}=0.\hfill \end{array}$$

Proof. 

Let $\mathfrak{X}=\mathfrak{Z}$ be same as Lemma 5, and $\mathcal{L},$$\mathcal{N},$$\mathcal{P},$$\mathcal{Q}$ be same as Lemmas 6 and 7. In what follows, we shall prove that model (1) has an $ɷ$-periodic solution based on Lemma 3.

First of all, we find the existence region $\Omega \subset \mathfrak{X}$ of solution. Assume that an $ɷ$-periodic solution $w={(\mathcal{U},\mathcal{V})}^T\in \mathfrak{X}$ solves the operator equation $\mathcal{L}w=\eta \mathcal{N}(w,\eta )$, then we obtain

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{\mathcal{U}}^{\Delta }\left(\varsigma \right)=\eta \left[r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right)e^{{\theta }_1\mathcal{U}\left(\varsigma \right)}-\eta c_{12}\left(\varsigma \right){\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}\Delta s-{\varphi }_1\left(\varsigma \right)e^{-\mathcal{U}\left(\varsigma \right)}\right],\hfill \\ {\mathcal{V}}^{\Delta }\left(\varsigma \right)=\eta \left[r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right)e^{{\theta }_2\mathcal{V}\left(\varsigma \right)}-\eta c_{21}\left(\varsigma \right){\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}\Delta s-{\varphi }_2\left(\varsigma \right)e^{-\mathcal{V}\left(\varsigma \right)}\right].\hfill \end{array}\right.\end{array}$$

(7)

Integrating at both sides of (7) yields

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{\int }_{\kappa }^{\kappa +ɷ}\left[r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right)e^{{\theta }_1\mathcal{U}\left(\varsigma \right)}-\eta c_{12}\left(\varsigma \right){\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}\Delta s-{\varphi }_1\left(\varsigma \right)e^{-\mathcal{U}\left(\varsigma \right)}\right]\Delta \varsigma =0,\hfill \\ {\int }_{\kappa }^{\kappa +ɷ}\left[r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right)e^{{\theta }_2\mathcal{V}\left(\varsigma \right)}-\eta c_{21}\left(\varsigma \right){\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}\Delta s-{\varphi }_2\left(\varsigma \right)e^{-\mathcal{V}\left(\varsigma \right)}\right]\Delta \varsigma =0.\hfill \end{array}\right.\end{array}$$

(8)

In view of periodicity of $\mathcal{U}\left(\varsigma \right)$ and $\mathcal{V}\left(\varsigma \right)$, there exist ${\mu }_1$, ${\mu }_2$, ${\nu }_1$ and ${\nu }_2\in I_{ɷ}$ satisfying $\mathcal{U}\left({\mu }_1\right)=\overline{\mathcal{U}}$, $\mathcal{U}\left({\mu }_2\right)=\underline{\mathcal{U}}$, $\mathcal{V}\left({\nu }_1\right)=\overline{\mathcal{V}}$, $\mathcal{V}\left({\nu }_2\right)=\underline{\mathcal{V}}$. The first equation in (7) and (8) leads

$$\begin{array}{c}\hfill {\int }_{\kappa }^{\kappa +ɷ}\left|{\mathcal{U}}^{\Delta }\left(\varsigma \right)\right|\Delta \varsigma <2ɷ\overline{r_1}.\end{array}$$

(9)

By the first equation of (8) and (9) together with Lemma 8, we have

$$\begin{array}{cc}\hfill ɷ\overline{r_1}\ge & {\int }_{\kappa }^{\kappa +ɷ}{r}_1\left(s\right)\Delta s={\int }_{\kappa }^{\kappa +ɷ}{c}_{11}\left(\varsigma \right)e^{{\theta }_1\mathcal{U}\left(\varsigma \right)}\Delta \varsigma +\eta {\int }_{\kappa }^{\kappa +ɷ}{c}_{12}\left(\varsigma \right)\hfill \\ \hfill & \times \left[{\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)e^{\mathcal{V}(\varsigma +s)}\Delta s\right]\Delta \varsigma +{\int }_{\kappa }^{\kappa +ɷ}{\varphi }_1\left(\varsigma \right)e^{-\mathcal{U}\left(\varsigma \right)}\Delta \varsigma \hfill \\ \hfill >& ɷ\underline{c_{11}}{e}^{{\theta }_1\mathcal{U}\left({\mu }_2\right)}+ɷ\underline{{\varphi }_1}{e}^{-\mathcal{U}\left({\mu }_1\right)}\ge ɷ\underline{c_{11}}{e}^{{\theta }_1[\mathcal{U}\left({\mu }_1\right)-ɷ\overline{r_1}]}+ɷ\underline{{\varphi }_1}{e}^{-\mathcal{U}\left({\mu }_1\right)},\hfill \end{array}$$

which implies that

$$\begin{array}{c}\hfill \left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)e^{(1+{\theta }_1)\mathcal{U}\left({\mu }_1\right)}-\overline{r_1}{e}^{\mathcal{U}\left({\mu }_1\right)}+\underline{{\varphi }_1}<0.\end{array}$$

(10)

We derive from Lemma 4 that there has unique $L_0\in \mathbb{R}$ such that

$$f_{min}=f\left(L_0\right)=-{\theta }_1{\left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)}^{-\frac{1}{{\theta }_1}}{\left(\frac{\overline{r_1}}{1+{\theta }_1}\right)}^{\frac{1+{\theta }_1}{{\theta }_1}}+\underline{{\varphi }_1},f^{\prime }\left(L_0\right)=0,$$

where $L_0=\frac{1}{{\theta }_1}ln\left[\frac{\overline{r_1}}{\left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)(1+{\theta }_1)}\right].$ From $\left({\mathrm{H}}_2\right)$ and Lemma 4, one knows that $f\left(L_0\right)<0$, and $f\left(x\right)=0$ has only two roots $L^{-}$ and $L^{+}$ satisfying

$$\begin{array}{c}\hfill L^{-}<L_0<L^{+},f\left(L^{\pm }\right)=0.\end{array}$$

(11)

By (11) and Lemma 4, the Inequality (10) is solved by

$$\begin{array}{c}\hfill L^{-}<\mathcal{U}\left({\mu }_1\right)<L^{+}.\end{array}$$

(12)

It follows from Lemma 8 and (12) that

$$\begin{array}{c}\hfill L^{-}-ɷ\overline{r_1}<\mathcal{U}\left({\mu }_2\right)\le \mathcal{U}\left({\mu }_1\right)<L^{+}.\end{array}$$

(13)

Similarly, the second equation of (7) gives

$$\begin{array}{cc}\hfill ɷ\overline{r_2}>& {\int }_{\kappa }^{\kappa +ɷ}{r}_2\left(\varsigma \right)\Delta \varsigma =\eta {\int }_{\kappa }^{\kappa +ɷ}{c}_{21}\left(\varsigma \right)\left[{\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)e^{\mathcal{U}(\varsigma +s)}\Delta s\right]\Delta \varsigma \hfill \\ \hfill +& {\int }_{\kappa }^{\kappa +ɷ}{c}_{22}\left(\varsigma \right)e^{{\theta }_2\mathcal{V}\left(\varsigma \right)}\Delta \varsigma +{\int }_{\kappa }^{\kappa +ɷ}{\varphi }_2\left(\varsigma \right)e^{-\mathcal{V}\left(\varsigma \right)}\Delta \varsigma \hfill \\ \hfill >& ɷ\underline{c_{22}}{e}^{{\theta }_2[\mathcal{V}\left({\nu }_1\right)-ɷ\overline{r_2}]}+ɷ\underline{{\varphi }_2}{e}^{-\mathcal{V}\left({\nu }_1\right)},\hfill \end{array}$$

which indicates that

$$\begin{array}{c}\hfill \left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)e^{(1+{\theta }_2)\mathcal{V}\left({\nu }_1\right)}-\overline{r_2}{e}^{\mathcal{V}\left({\nu }_1\right)}+\underline{{\varphi }_2}<0.\end{array}$$

(14)

From Lemma 4, one knows that there has unique $M_0\in \mathbb{R}$ such that

$$g_{min}=g\left(M_0\right)=-{\theta }_2{\left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)}^{-\frac{1}{{\theta }_2}}{\left(\frac{\overline{r_2}}{1+{\theta }_2}\right)}^{\frac{1+{\theta }_2}{{\theta }_2}}+\underline{{\varphi }_2},g^{\prime }\left(M_0\right)=0,$$

where $M_0=\frac{1}{{\theta }_2}ln\left[\frac{\overline{r_2}}{\left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)(1+{\theta }_2)}\right].$ It follows from $\left({\mathrm{H}}_2\right)$ and Lemma 4 that $g\left(M_0\right)<0$, and $g\left(x\right)=0$ has two roots $M^{-}$ and $M^{+}$ satisfying

$$\begin{array}{c}\hfill M^{-}<M_0<M^{+},g\left(M^{\pm }\right)=0.\end{array}$$

(15)

By (15) and Lemma 4, the solution of Inequality (14) is read as

$$\begin{array}{c}\hfill M^{-}<\mathcal{V}\left({\nu }_1\right)<M^{+}.\end{array}$$

(16)

From Lemma 8 and (16), we have

$$\begin{array}{c}\hfill M^{-}-ɷ\overline{r_2}<\mathcal{V}\left({\nu }_2\right)\le \mathcal{V}\left({\nu }_1\right)<M^{+}.\end{array}$$

(17)

Based on (13) and (17), we take

$$\begin{array}{c}\hfill \Omega =\left\{{(\mathcal{U}\left(\varsigma \right),\mathcal{V}\left(\varsigma \right))}^T\in \mathfrak{X}:L^{-}-ɷ\overline{r_1}<\mathcal{U}\left(\varsigma \right)<L^{+},M^{-}-ɷ\overline{r_2}<\mathcal{V}\left(\varsigma \right)<M^{+}\right\}.\end{array}$$

Obviously, $\Omega \subset \mathfrak{X}$ meets with Lemma 3 $\left(\mathrm{a}\right)$.

We next adopt the reduction to absurdity to prove that Lemma 3 $\left(\mathrm{b}\right)$ hold, i.e., $w\in \partial \Omega \cap \mathrm{Ker}\left(\mathcal{L}\right)=\partial \Omega \cap {\mathbb{R}}^2$ implies $\mathcal{Q}\mathcal{N}(w,0)\ne (0,0)$. Indeed, assume that the conclusion is opposite, then there has a constant vector $W^{\ast }=(U^{\ast },V^{\ast })$ with $W^{\ast }\in \partial \Omega \cap {\mathbb{R}}^2$ that fulfills

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{\int }_{\kappa }^{\kappa +ɷ}\left[r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right)e^{{\theta }_1{U}^{\ast }}-c_{12}\left(\varsigma \right)e^{V^{\ast }}{\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_1\left(s\right)\Delta s-{\varphi }_1\left(\varsigma \right)e^{-U^{\ast }}\right]\Delta \varsigma =0,\hfill \\ {\int }_{\kappa }^{\kappa +ɷ}\left[r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right)e^{{\theta }_2{U}^{\ast }}-c_{21}\left(\varsigma \right)e^{U^{\ast }}{\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)\Delta s-{\varphi }_2\left(\varsigma \right)e^{-V^{\ast }}\right]\Delta \varsigma =0.\hfill \end{array}\right.\end{array}$$

(18)

A discussion for (18) analogue to (8)–(17) yields $W^{\ast }=(U^{\ast },V^{\ast })\in \Omega \cap {\mathbb{R}}^2$, this is contrary to $W^{\ast }=(U^{\ast },V^{\ast })\in \partial \Omega \cap {\mathbb{R}}^2$. So the Lemma 3 $\left(\mathrm{b}\right)$ holds.

We choose the identity operator $\mathcal{J}=\mathcal{I}$ and calculate it directly

$$\begin{array}{c}\hfill deg\left\{\mathcal{J}\mathcal{Q}\mathcal{N}(w,0),\Omega \cap \mathrm{Ker}\left(\mathcal{J}\right),{(0,0)}^T\right\}\ne 0,\end{array}$$

which means that Lemma 3 $\left(\mathrm{c}\right)$ holds. Thus, one concludes from Lemma 3 that there has at least an $ɷ$-periodic function ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ satisfying model (1). The proof is completed. □

4. Global Asymptotic Stability

The portion concentrates on the global asymptotic stability of model (1). To this end, we need the following definitions.

Definition 5.

According to Lyapunov stability theory, we call that the ϖ-periodic positive solution $\tilde{w}\left(\varsigma \right)=(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))$ of (1) is globally asymptotically stable on a periodic time scale $\mathbb{T}$ iff, for any positive solution $w\left(\varsigma \right)=\left(\mathcal{U}\right(\varsigma ),\mathcal{V}(\varsigma \left)\right)$ of (1), and for any large real number $R>0$, when $\parallel \tilde{w}\left(0\right)-w\left(0\right)\parallel <R$, we have ${lim}_{\varsigma \to +\infty }\parallel \tilde{w}\left(\varsigma \right)-w\left(\varsigma \right)\parallel =0$.

Definition 6

([28]). Let $U_{\varsigma }$ be a neighborhood of ς, $\forall \varsigma \in \mathbb{T}$, and $V\in C_{rd}(\mathbb{T}\times {\mathbb{R}}^n,{\mathbb{R}}^{+})$, $D^{+}{V}^{\Delta }(\varsigma ,x\left(\varsigma \right))$ is called the Dini derivative of V iff, $\forall ϵ>0$, there has a right neighborhood $U_{ϵ}\cap U_{\varsigma }$ of ς satisfying

$$\frac{V\left(\sigma \right(\varsigma ),x(\sigma \left(\varsigma \right)\left)\right)-V(s,x(s\left)\right)}{\sigma \left(\varsigma \right)-s}<D^{+}{V}^{\Delta }(\varsigma ,x\left(\varsigma \right))+ϵ,\forall s\in U_{ϵ},s>\varsigma .$$

If $V(\varsigma ,x(\varsigma \left)\right)$ is continuous at right-scattered ς, then

$$D^{+}{V}^{\Delta }(\varsigma ,x\left(\varsigma \right))=\frac{V\left(\sigma \right(\varsigma ),x(\sigma \left(\varsigma \right)\left)\right)-V(\varsigma ,x(\varsigma \left)\right)}{\sigma \left(\varsigma \right)-\varsigma }.$$

By Theorem 1, we know that there has an $ɷ$-periodic function ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T\in \Omega$ satisfying model (1). Let $u\left(\varsigma \right)=e^{\mathcal{U}\left(\varsigma \right)}$, $v\left(\varsigma \right)=e^{\mathcal{V}\left(\varsigma \right)}$, then one has ${\mathcal{U}}^{\Delta }\left(\varsigma \right)={(lnu\left(\varsigma \right))}^{\Delta }$ and ${\mathcal{V}}^{\Delta }\left(\varsigma \right)={(lnv\left(\varsigma \right))}^{\Delta }$. Thus the system (1) becomes

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{(lnu\left(\varsigma \right))}^{\Delta }=r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right){\left[u\left(\varsigma \right)\right]}^{{\theta }_1}-c_{12}\left(\varsigma \right){\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)v(\varsigma +s)\Delta s-\frac{{\varphi }_1\left(\varsigma \right)}{u\left(\varsigma \right)},\varsigma \in \mathbb{T},\hfill \\ {(lnv\left(\varsigma \right))}^{\Delta }=r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right){\left[v\left(\varsigma \right)\right]}^{{\theta }_2}-c_{21}\left(\varsigma \right){\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)u(\varsigma +s)\Delta s-\frac{{\varphi }_2\left(\varsigma \right)}{v\left(\varsigma \right)},\varsigma \in \mathbb{T}.\hfill \end{array}\right.\end{array}$$

(19)

Similarly, there exists an $ɷ$-periodic positive function ${(\tilde{u}\left(\varsigma \right),\tilde{v}\left(\varsigma \right))}^T\in \tilde{\Omega }$ that solves system (19), here

$$\begin{array}{c}\hfill \tilde{\Omega }=\left\{{(u\left(\varsigma \right),v\left(\varsigma \right))}^T:e^{L^{-}-ɷ\overline{r_1}}<u\left(\varsigma \right)<e^{L^{+}},e^{M^{-}-ɷ\overline{r_2}}<v\left(\varsigma \right)<e^{M^{+}}\right\}.\end{array}$$

Choose the constants $\rho >0$ and $\theta >0$ such that

$$0<\rho <min\left\{e^{L^{-}-ɷ\overline{r_1}},e^{M^{-}-ɷ\overline{r_2}}\right\},\theta \ge max\left\{1,{\theta }_1,{\theta }_2\right\}.$$

We further assume that

$\left({\mathrm{H}}_3\right)$

The followings are true:

$$-{\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_1}+{\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\Delta s<0,$$

$$-{\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_2}+{\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\Delta s<0.$$

Taking variable substitution $u\left(\varsigma \right)={\left(\rho X\left(\varsigma \right)\right)}^{\frac{1}{\theta }}$ and $v\left(\varsigma \right)={\left(\rho Y\left(\varsigma \right)\right)}^{\frac{1}{\theta }}$, then we have

$${(lnu\left(\varsigma \right))}^{\Delta }={\left[\frac{1}{\theta }ln\left(\rho X\left(\varsigma \right)\right)\right]}^{\Delta }=\frac{1}{\theta }{\left[ln\rho +lnX\left(\varsigma \right)\right]}^{\Delta }=\frac{1}{\theta }{(lnX\left(\varsigma \right))}^{\Delta },$$

and

$${(lnv\left(\varsigma \right))}^{\Delta }={\left[\frac{1}{\theta }ln\left(\rho Y\left(\varsigma \right)\right)\right]}^{\Delta }=\frac{1}{\theta }{\left[ln\rho +lnY\left(\varsigma \right)\right]}^{\Delta }=\frac{1}{\theta }{(lnY\left(\varsigma \right))}^{\Delta }.$$

Consequently, system (19) changes into

$$\begin{array}{c}\hfill \left\{\begin{array}{c}{(lnX\left(\varsigma \right))}^{\Delta }=\theta [r_1\left(\varsigma \right)-{\rho }^{\frac{{\theta }_1}{\theta }}{c}_{11}\left(\varsigma \right)X^{\frac{{\theta }_1}{\theta }}\left(\varsigma \right)-{\rho }^{\frac{1}{\theta }}{c}_{12}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)Y^{\frac{1}{\theta }}(\varsigma +s)\Delta s-{\rho }^{-\frac{1}{\theta }}{\varphi }_1\left(\varsigma \right)X^{-\frac{1}{\theta }}\left(\varsigma \right)],\hfill \\ {(lnY\left(\varsigma \right))}^{\Delta }=\theta [r_2\left(\varsigma \right)-{\rho }^{\frac{{\theta }_2}{\theta }}{c}_{22}\left(\varsigma \right)Y^{\frac{{\theta }_2}{\theta }}\left(\varsigma \right)-{\rho }^{\frac{1}{\theta }}{c}_{21}\left(\varsigma \right)\hfill \\ \times {\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)X^{\frac{1}{\theta }}(\varsigma +s)\Delta s-{\rho }^{-\frac{1}{\theta }}{\varphi }_2\left(\varsigma \right)Y^{-\frac{1}{\theta }}\left(\varsigma \right)].\hfill \end{array}\right.\end{array}$$

(20)

Clearly, system (20) has an $ɷ$-periodic positive solution ${(\tilde{X}\left(\varsigma \right),\tilde{Y}\left(\varsigma \right))}^T={(\frac{1}{\rho }{\tilde{u}}^{\theta }\left(\varsigma \right),\frac{1}{\rho }{\tilde{v}}^{\theta }\left(\varsigma \right))}^T$ within ${\tilde{\Omega }}^{\prime }$, where

$$\begin{array}{c}\hfill {\tilde{\Omega }}^{\prime }=\left\{{(X\left(\varsigma \right),Y\left(\varsigma \right))}^T:\frac{1}{\rho }{e}^{\theta (L^{-}-ɷ\overline{r_1})}<X\left(\varsigma \right)<\frac{1}{\rho }{e}^{\theta L^{+}},\frac{1}{\rho }{e}^{\theta (M^{-}-ɷ\overline{r_2})}<Y\left(\varsigma \right)<\frac{1}{\rho }{e}^{\theta M^{+}}\right\}.\end{array}$$

From Theorem 1 and ${\tilde{\Omega }}^{\prime }$, we have

$$\begin{array}{c}\hfill 1<\frac{1}{\rho }{e}^{\theta (L^{-}-ɷ\overline{r_1})}<\tilde{X}\left(\varsigma \right)<\frac{1}{\rho }{e}^{\theta L^{+}},1<\frac{1}{\rho }{e}^{\theta (M^{-}-ɷ\overline{r_2})}<\tilde{Y}\left(\varsigma \right)<\frac{1}{\rho }{e}^{\theta M^{+}}.\end{array}$$

(21)

Theorem 2.

Assume that $\left({\mathrm{H}}_1\right)$–$\left({\mathrm{H}}_3\right)$ hold, then a unique ϖ-periodic solution ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ of (1) is globally asymptotically stable.

Proof. 

Assume that the $ɷ$-periodic function ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ satisfying (1) is globally asymptotically stable, then ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ is attractive, that is, for each function ${(\mathcal{U}\left(\varsigma \right),\mathcal{V}\left(\varsigma \right))}^T$ satisfying (1), we have $\underset{\varsigma \to +\infty }{lim}[\mathcal{U}\left(\varsigma \right)-\tilde{\mathcal{U}}\left(\varsigma \right)]=0$, $\underset{\varsigma \to +\infty }{lim}[\mathcal{V}\left(\varsigma \right)-\tilde{\mathcal{V}}\left(\varsigma \right)]=0$. If the system (1) has another $ɷ$-periodic solution ${({\mathcal{U}}^{\ast }\left(\varsigma \right),{\mathcal{V}}^{\ast }\left(\varsigma \right))}^T\in \Omega$ with ${({\mathcal{U}}^{\ast }\left(\varsigma \right),{\mathcal{V}}^{\ast }\left(\varsigma \right))}^T\ne {(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$, without loss of generality, assume that ${\mathcal{U}}^{\ast }\left(\varsigma \right)\ne \tilde{\mathcal{U}}\left(\varsigma \right)$, then we obtain $0<|\tilde{\mathcal{U}}\left(\varsigma \right)-{\mathcal{U}}^{\ast }\left(\varsigma \right)|\le |\tilde{\mathcal{U}}\left(\varsigma \right)-\mathcal{U}\left(\varsigma \right)|+|\mathcal{U}\left(\varsigma \right)-{\mathcal{U}}^{\ast }\left(\varsigma \right)|\to 0$, as $\varsigma \to +\infty$, which is an obvious fallacy. Thus, we prove that the $ɷ$-periodic function ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ satisfying model (1) is unique provided that ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ is globally asymptotically stable. In addition, since the global asymptotical stability of $ɷ$-periodic function ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T\in \Omega$ satisfying (1) and ${(\tilde{X}\left(\varsigma \right),\tilde{Y}\left(\varsigma \right))}^T$ satisfying (20) is equivalent, we just need to show that the $ɷ$-periodic function ${(\tilde{X}\left(\varsigma \right),\tilde{Y}\left(\varsigma \right))}^T$ satisfying (20) is globally asymptotically stable. Indeed, it follows from $\left({\mathrm{H}}_1\right)$, $\left({\mathrm{H}}_2\right)$ and Theorem 1 that system (20) has an positive $ɷ$-periodic solution ${(\tilde{X}\left(\varsigma \right),\tilde{Y}\left(\varsigma \right))}^T$. For each positive solution ${(X\left(\varsigma \right),Y\left(\varsigma \right))}^T$ of (20), build a Lyapunov functional $V\left(\varsigma \right)=V_1\left(\varsigma \right)+V_2\left(\varsigma \right),$ here

$$\begin{array}{c}\hfill V_1\left(\varsigma \right)=|lnX\left(\varsigma \right)-ln\tilde{X}\left(\varsigma \right)|+|lnY\left(\varsigma \right)-ln\tilde{Y}\left(\varsigma \right)|,\end{array}$$

(22)

$$\begin{array}{cc}\hfill V_2\left(\varsigma \right)=& \theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\left[{\int }_{\varsigma +s}^{\varsigma }|Y^{\frac{1}{\theta }}\left(\zeta \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\zeta \right)|\Delta \zeta \right]\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\left[{\int }_{\varsigma +s}^{\varsigma }|X^{\frac{1}{\theta }}\left(\zeta \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\zeta \right)|\Delta \zeta \right]\Delta s.\hfill \end{array}$$

(23)

Obviously, $V\left(0\right)<+\infty$ and $V\left(\varsigma \right)\ge V_1\left(\varsigma \right)$. By (21), a direct $\Delta$-derivation along (20) gives

$$\begin{array}{cc}\hfill & D^{+}\left(\right|lnX\left(\varsigma \right)-ln\tilde{X}\left(\varsigma \right){\left|\right)}^{\Delta }\le -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s+\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X^{\frac{-1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{-1}{\theta }}\left(\varsigma \right)|\hfill \\ \hfill \le & -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)|X^{\frac{-1}{\theta }}\left(\varsigma \right){\tilde{X}}^{\frac{-1}{\theta }}\left(\varsigma \right)\hfill \\ \hfill \le & -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)|,\hfill \end{array}$$

(24)

$$\begin{array}{cc}\hfill & D^{+}\left(\right|lnY\left(\varsigma \right)-ln\tilde{Y}\left(\varsigma \right){\left|\right)}^{\Delta }\le -\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s+\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y^{\frac{-1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{-1}{\theta }}\left(\varsigma \right)|\hfill \\ \hfill \le & -\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)|Y^{\frac{-1}{\theta }}\left(\varsigma \right){\tilde{Y}}^{\frac{-1}{\theta }}\left(\varsigma \right)\hfill \\ \hfill \le & -\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)|,\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill & D^{+}{\left({\int }_{-\xi }^0{k}_1\left(s\right)\left[{\int }_{\varsigma +s}^{\varsigma }|Y^{\frac{1}{\theta }}\left(\zeta \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\zeta \right)|\Delta \zeta \right]\Delta s\right)}^{\Delta }={\int }_{-\xi }^0{k}_1\left(s\right)\Delta s·\left|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|\hfill \\ \hfill & -{\int }_{-\xi }^0{k}_1\left(s\right)\left|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)\right|\Delta s,\hfill \end{array}$$

(26)

and

$$\begin{array}{cc}\hfill & D^{+}{\left({\int }_{-\xi }^0{k}_2\left(s\right)\left[{\int }_{\varsigma +s}^{\varsigma }|X^{\frac{1}{\theta }}\left(\zeta \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\zeta \right)|\Delta \zeta \right]\Delta s\right)}^{\Delta }={\int }_{-\xi }^0{k}_2\left(s\right)\Delta s·\left|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|\hfill \\ \hfill & -{\int }_{-\xi }^0{k}_2\left(s\right)\left|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)\right|\Delta s.\hfill \end{array}$$

(27)

Since when $\alpha ,\beta >0$ and $x\ge 1$ $q\left(x\right)=|{\alpha }^x-{\beta }^x|$ is monotonically increasing, and $0<\frac{{\theta }_1}{\theta },\frac{{\theta }_2}{\theta },\frac{1}{\theta }\le 1$, it follows from (21), (24)–(27) and $\left({\mathrm{H}}_3\right)$ that

$$\begin{array}{cc}\hfill & D^{+}{V}^{\Delta }\left(\varsigma \right)\le -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s\hfill \end{array}$$

$$\begin{array}{cc}\hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)|-\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)|\Delta s+\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\Delta s·\left|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|-\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\left|Y^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{Y}}^{\frac{1}{\theta }}(\varsigma +s)\right|\Delta s\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\Delta s·\left|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|-\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\left|X^{\frac{1}{\theta }}(\varsigma +s)-{\tilde{X}}^{\frac{1}{\theta }}(\varsigma +s)\right|\Delta s\hfill \\ \hfill =& -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|+\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)|-\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\Delta s·\left|Y^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{Y}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\Delta s·\left|X^{\frac{1}{\theta }}\left(\varsigma \right)-{\tilde{X}}^{\frac{1}{\theta }}\left(\varsigma \right)\right|\hfill \\ \hfill \le & -\theta {\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|+\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_1}|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)|-\theta {\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|\hfill \\ \hfill & +\theta {\rho }^{\frac{-1}{\theta }}\overline{{\varphi }_2}|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)|+\theta {\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\Delta s·\left|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)\right|\hfill \\ \hfill & +\theta {\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\Delta s·\left|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)\right|\hfill \\ \hfill =& \theta \left[-{\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_1}+{\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)\Delta s\right]\left|X\left(\varsigma \right)-\tilde{X}\left(\varsigma \right)\right|\hfill \\ \hfill & +\theta \left[-{\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_2}+{\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)\Delta s\right]\left|Y\left(\varsigma \right)-\tilde{Y}\left(\varsigma \right)\right|<0.\hfill \end{array}$$

(28)

Thus, from (22), (23) and (28), one concludes that $V\left(\varsigma \right)$ is positive definite and $D^{+}{V}^{\Delta }\left(\varsigma \right)<0$, $\forall \varsigma \ge 0$. Therefore, one draws a conclusion that the $ɷ$-periodic solution ${(\tilde{X}\left(\varsigma \right),\tilde{Y}\left(\varsigma \right))}^T$ of system (20) has global asymptotic stability based on Lyapunov stability theory. The proof is completed. □

5. Numerical Simulation

The portion considers the following nonlinear Ayala-Gilpin competitive ecosystem having distributed lags on time scale $\mathbb{T}=\mathbb{R}$

$$\begin{array}{c}\hfill \left\{\begin{array}{c}\frac{d\mathcal{U}\left(\varsigma \right)}{d\varsigma }=\mathcal{U}\left(\varsigma \right)\left[r_1\left(\varsigma \right)-c_{11}\left(\varsigma \right){\mathcal{U}}^{{\theta }_1}\left(\varsigma \right)-c_{12}\left(\varsigma \right){\int }_{-{\xi }_1\left(\varsigma \right)}^0{k}_1\left(s\right)\mathcal{V}(\varsigma +s)ds\right]-{\varphi }_1\left(\varsigma \right),\hfill \\ \frac{d\mathcal{V}\left(\varsigma \right)}{d\varsigma }=\mathcal{V}\left(\varsigma \right)\left[r_2\left(\varsigma \right)-c_{22}\left(\varsigma \right){\mathcal{V}}^{{\theta }_2}\left(\varsigma \right)-c_{21}\left(\varsigma \right){\int }_{-{\xi }_2\left(\varsigma \right)}^0{k}_2\left(s\right)\mathcal{U}(\varsigma +s)ds\right]-{\varphi }_2\left(\varsigma \right),\hfill \\ \mathcal{U}\left(\varsigma \right)={\phi }_1\left(\varsigma \right),\mathcal{V}\left(\varsigma \right)={\phi }_2\left(\varsigma \right),\varsigma \in (-\xi ,0],\hfill \end{array}\right.\end{array}$$

(29)

where $r_1\left(\varsigma \right)=8+2cos\left(3\varsigma \right)$, $r_2\left(\varsigma \right)=6+sin\left(2\varsigma \right)$, $c_{11}\left(\varsigma \right)=5+2sin\left(\varsigma \right)$, $c_{22}\left(\varsigma \right)=3+cos\left(2\varsigma \right)$, $c_{12}\left(\varsigma \right)=\frac{6+sin\left(2\varsigma \right)}{10}$, $c_{21}\left(\varsigma \right)=\frac{3+cos\left(3\varsigma \right)}{10}$, ${\xi }_1\left(\varsigma \right)=k_1\left(\varsigma \right)=\frac{2+sin\left(\varsigma \right)}{5}$, ${\xi }_2\left(t\right)=k_2\left(\varsigma \right)=\frac{2+cos\left(\varsigma \right)}{4}$, ${\varphi }_1\left(\varsigma \right)=\frac{3+cos\left(2\varsigma \right)}{7}$, ${\varphi }_2\left(\varsigma \right)=\frac{4+sin\left(\varsigma \right)}{7}$, ${\theta }_1=\frac{1}{2}$, ${\theta }_2=\frac{1}{\sqrt{2}}$. Take the initial functions ${\phi }_1\left(\varsigma \right)=7+sin\left(\varsigma \right)$, ${\phi }_2\left(\varsigma \right)=\frac{2+cos\left(\varsigma \right)}{7}$, $\varsigma \in (-\xi ,0]=(-\frac{3}{4},0]$, here $\xi =max\{{\overline{\xi }}_1,{\overline{\xi }}_2\}=max\{\frac{3}{5},\frac{3}{4}\}=\frac{3}{4}$.

Obviously, $r_1\left(\varsigma \right)$, $r_2\left(\varsigma \right)$, $c_{11}\left(\varsigma \right)$, $c_{22}\left(\varsigma \right)$, $c_{12}\left(\varsigma \right)$, $c_{21}\left(\varsigma \right)$, ${\xi }_1\left(\varsigma \right)$, ${\xi }_2\left(\varsigma \right)$, ${\phi }_1\left(\varsigma \right)$, ${\phi }_2\left(\varsigma \right)$, $k_1\left(\varsigma \right)$, $k_2\left(\varsigma \right)$, ${\varphi }_1\left(\varsigma \right)$ and ${\varphi }_2\left(\varsigma \right)$ are all positive periodic functions with period $ɷ=2\pi$. So the conditions $\left({\mathrm{H}}_1\right)$ holds. A direct computation gives $\overline{r_1}=10$, $\underline{r_1}=6$, $\overline{r_2}=7$, $\underline{r_2}=5$, $\overline{c_{11}}=7$, $\underline{c_{11}}=3$, $\overline{c_{22}}=4$, $\underline{c_{22}}=2$, $\overline{c_{12}}=\frac{7}{10}$, $\underline{c_{12}}=\frac{1}{2}$, $\overline{c_{21}}=\frac{2}{5}$, $\underline{c_{21}}=\frac{1}{5}$, $\overline{{\varphi }_1}=\frac{4}{7}$, $\underline{{\varphi }_1}=\frac{2}{7}$, $\overline{{\varphi }_2}=\frac{5}{7}$, $\underline{{\varphi }_2}=\frac{3}{7}$, $\xi =\frac{3}{4}$, ${\int }_{-\xi }^0{k}_1\left(s\right)ds\approx 0.2463$, ${\int }_{-\xi }^0{k}_2\left(s\right)ds\approx 0.5454$. To solve the following algebraic equation

$$f\left(L^{\pm }\right)=\left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)e^{(1+{\theta }_1)L^{\pm }}-\overline{r_1}{e}^{L^{\pm }}+\underline{{\varphi }_1}=0,$$

we find that $L^{-}\approx -3.5553$, $L^{+}\approx 65.2398$. From the algebraic equation

$$g\left(M^{\pm }\right)=\left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)e^{(1+{\theta }_2)M^{\pm }}-\overline{r_2}{e}^{M^{\pm }}+\underline{{\varphi }_2}=0,$$

we have $M^{-}\approx -2.7932$, $M^{+}\approx 45.7540$. Thus we get

$$\begin{array}{cc}\hfill \Omega =\{& {(\mathcal{U}\left(\varsigma \right),\mathcal{V}\left(\varsigma \right))}^T:1.4738\times {10}^{-29}<\mathcal{U}\left(\varsigma \right)<2.1542\times {10}^{28},\hfill \\ \hfill & 4.8491\times {10}^{-21}<\mathcal{V}\left(\varsigma \right)<7.4252\times {10}^{19}\}.\hfill \end{array}$$

The condition $\left({\mathrm{H}}_2\right)$ is verified as

$${\theta }_1{\left(\underline{c_{11}}{e}^{-ɷ{\theta }_1\overline{r_1}}\right)}^{-\frac{1}{{\theta }_1}}{\left(\frac{\overline{r_1}}{1+{\theta }_1}\right)}^{\frac{1+{\theta }_1}{{\theta }_1}}\approx 3.1914\times {10}^{28}>\underline{{\varphi }_1}=\frac{2}{7},$$

$${\theta }_2{\left(\underline{c_{22}}{e}^{-ɷ{\theta }_2\overline{r_2}}\right)}^{-\frac{1}{{\theta }_2}}{\left(\frac{\overline{r_2}}{1+{\theta }_2}\right)}^{\frac{1+{\theta }_2}{{\theta }_2}}\approx 1.0105\times {10}^{20}>\underline{{\varphi }_2}=\frac{3}{7}.$$

By now, $\left({\mathrm{H}}_1\right)$ and $\left({\mathrm{H}}_2\right)$ have been verified. From Theorem 1, one knows that (29) exists at least an $2\pi$-periodic positive solutions ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T\in \Omega$.

Next, we prove periodic positive solution ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ to be globally asymptotically stable. Indeed, take $\rho =1.5\times {10}^{-29}$, $\theta =120$, we get

$$-{\rho }^{\frac{{\theta }_1}{\theta }}\underline{c_{11}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_1}+{\rho }^{\frac{1}{\theta }}\overline{c_{21}}{\int }_{-\xi }^0{k}_2\left(s\right)ds\approx -1.1551<0,$$

and

$$-{\rho }^{\frac{{\theta }_2}{\theta }}\underline{c_{22}}+{\rho }^{-\frac{1}{\theta }}\overline{{\varphi }_2}+{\rho }^{\frac{1}{\theta }}\overline{c_{12}}{\int }_{-\xi }^0{k}_1\left(s\right)ds\approx -0.0104<0.$$

Thus the condition $\left({\mathrm{H}}_3\right)$ holds. From Theorem 2, one knows that the periodic solution ${(\tilde{\mathcal{U}}\left(\varsigma \right),\tilde{\mathcal{V}}\left(\varsigma \right))}^T$ is globally asymptotically stable. We have carried out numerical simulation for example (29) as shown in Figure 1.

6. Summaries and Outlooks

The Ayala-Gilpin dfferential equation model is one of the successful patterns of applying mathematical theories and methods to study ecosystems. Its dynamic behavior has received significant attention and study from mathematicians and ecologists. In this work, we mainly investigate the existence and global asymptotic stability of periodic solutions for a class of nonlinear distributed-lag Ayala-Gilpin vying system (1) in the sense of time scales. By making use of Mawhin’s coincidence degree theorem, we first gain some sufficient criteria for the existence of periodic solutions of model (1). Next, we construct an appropriate Lyapunov functional to demonstrate that model (1) is globally asymptotically stable. Subsequently, we examined the validity and applicability of our essential findings through theoretical analysis and numerical simulation of an example. Our conclusion reveals the existence of periodic oscillations in AG-ecosystem under certain conditions from a mathematical perspective. However, the long-term behavior of AG-ecosystem is globally asymptotically stable. As stated in studies [10,11], an ecosystem may have multiple stable states for different initial values, which means that the ecosystem has multiple positive solutions and is stable in their respective regions. Therefore, we can study the existence and local stability of multiple positive periodic solutions for model (1) in future works. In addition, some scholars have found that using fractional or partial differential models to study some practical problems is more accurate than using integer order differential models. Awaken by some recent papers [29,30,31,32,33,34,35,36,37,38,39,40,41], we plan to apply fractional calculus and PDE theory to further study the AG-ecosystem in the future, in order to explore more dynamic characteristics.