Abstract

In this paper, we study the stability of positive steady states in a delayed competition system on a weighted network, which does not satisfy the comparison principle appealing to classical competitive systems. By introducing some auxiliary equations and constructing proper contracting rectangles, we present some sufficient conditions on the stability of the unique positive steady state. Moreover, some numerical examples are given to explore the complex dynamics of this nonmonotone model, which implies the nontrivial roles of weights and time delays.

1. Introduction

Competition is one of the most universal phenomena in the natural world due to the limit resources including water and sunshine. When the amount of individuals is concerned in population dynamics, some differential systems modeling both interspecific and intraspecific competitions have been proposed and studied. For example, if only two competitive species are concerned in population dynamics, then one famous model isin which all the parameters are positive, and present the densities of two competitive species; we may refer to Murray [1], section 3.5, for the dynamics and biology background of (1). When multispecies and delayed effect [24] were considered, Faria [5] and Smith [6] studied the following delayed system:in which , all the parameters are positive, and formulates the history delayed effect and satisfies

In some cases, the spatial distribution of individuals must be considered. For example, to control the disease spreading, the spatial distribution of the infected is important. To model the spatial distribution of individuals, many reaction-diffusion systems have been established and studied; see [710]. In particular, we refer to He and Ni [11] and references cited therein for some results on the stability of different steady states in the corresponding reaction-diffusion system of (1). Moreover, when the spatial habitat is , Weinberger et al. [12] studied the invasion dynamics between the invader and the native in the corresponding diffusion model of (1). For the corresponding diffusion model of (2), Martin and Smith [13] proved the stability of steady states in the bounded domain. In the classical reaction-diffusion systems, it is necessary to assume that the density of individuals is continuous in both spatial and temporal variables. However, the spatial continuity is difficult to be satisfied due to the limitation of data in some cases. To overcome the deficiency, one recipe is to study the ordinary differential equations on a weighted network by considering finite patches as habitat of individuals; we may refer to Beretta et al. [14], Faria [5], Gao and Ruan [15], Liao and Lou [16], Mai et al. [17], Sun and Mai [18], Xu and Chen [19], and Zhang and Wang [20] for patch/network models with time delay.

Recently, Liu et al. [21] considered the competition dynamics of two competitors on a network or patches ( and ); their model can be written asin which describes the movement between the -th and -th patches of the -th species, which is called the weight between patches . In [18], Sun and Mai studied the following system on a network:in which are time delays. A natural question is how to further study the dynamics of multicompetitive species with delayed effect on networks. In this paper, we study the following competition system with time delays on the network:in which all the parameters are positive, , denotes the population density of the -th species on the -th patch, and satisfies

That is, (6) is a coupled system of functional differential equations. Evidently, there are patches in this model, the same species on different patches may move from one to another due to , and different species on the same patch compete each other. So, (6) is a competitive system on a weighted network or patch environment.

To compare the literature results with the complex phenomena in the real world, the stability of mathematical models is important since the disturbance is inevitable. Moreover, the controllability of many mathematical models in finite time is also important, which partly depends on the stability conditions of some states; see a number of examples about switched nonlinear systems [2224]. Liu et al. [21] and Sun and Mai [18] studied the stability of different steady states of (4) and (5) by using the comparison principle and other techniques, which also include some applications in population dynamics. However, because of the intraspecific delayed effect, (6) does not satisfy the comparison principle of classical competitive systems, and the recipe in [21] does not work directly. Moreover, since there are functional differential equations, the analysis of eigenvalues would be complex if we analyse the stability by the characteristic equation at the positive steady state, which may involve transcendental equations. Therefore, to study the dynamics of (6), some other techniques are necessary.

Evidently, (6) may have several different steady states under proper conditions ((1) may have four steady states [1]), which reflect different biological sense in population dynamics. Coexistence of competitive species has been widely studied in many cases due to its important biological backgrounds on the biodiversity [25, 26]. In this paper, by combining the technique in [21] with the contracting rectangles in [6], we study the stability of the positive steady state by giving some sufficient conditions. Moreover, we also show some numerical results to support our theoretical conclusions and present the complex dynamics of this system when our stability conditions do not hold, which implies the nontrivial role of time delays and partial degenerate weights.

The rest of this paper is organized as follows. In section 2, we present some preliminaries on the stability via contracting rectangles. In Section 3, we show our main results on the stability of the unique positive steady state. Some numerical examples will be given in Section 4. Furthermore, we make a discussion in Section 5.

2. Preliminaries

In this paper, we use the standard partial ordering in . That is, ifthen

Moreover, implies , but , and implies . , denote the set of uniform continuous vector functions from to .

In this part, we recall the stability of the steady state by contracting rectangles [6]. Consider the functional differential systemin which is the delay, , and satisfies the following assumptions:(H1) There exists with such that where denotes the constant-valued function in .(H2) There exists a one-parameter family of ordered intervals given bysuch that , and for ,where and are continuous in .(H3) If and for , thenis Lipschitz continuous in the sense of supremum norm for each . (H4) is a strict contracting rectangle, namely, let

Then, for any and , we have

By contracting rectangles, we have the following conclusion [6].

Lemma 1. Assume that (H1)–(H4) hold. If for some with , then

3. Main Results

In this section, we study the long-time behavior of (6). We denoteand set

In population dynamics, implies the existence of the instantaneous self-limitation effect.

By the basic theory of functional differential equations [6, 13], we have the following global existence of (6), in which the unknown function is

Lemma 2. Assume that . Ifthen (6) has a global solution such that

Moreover, ifthen (6) has a positive steady state

Once the global existence is obtained, we shall prove the following result on the global asymptotic stability of the positive steady state, which also implies the uniqueness of the positive steady state of (6).

Theorem 1. Assume that , such that (20) holds andFor each , if there exists such that , thenIn the following, we shall prove Theorem 1 by several lemmas. The first is the following result on the positivity of solutions.

Lemma 3. Assume that for some . Ifthen , and

Proof. In the proof, we fix for simplicity. By , it is clear that . We now prove (27) by Lemma 2. From the positivity, we haveFurther considerThen, the comparison principle leads tosince , implies that (29) is a cooperative, irreducible system. By the theory of monotone dynamical systems [6], we haveUsing the comparison principle of cooperative systems, we have (27) for such a fix . We complete the proof.
Consider the initial value problemAfter direct calculation, we have the following conclusion by Lemma 1 (we also refer to Smith [6], Lemma 7.4, and Appendix).

Lemma 4. Assume that (24) holds. Definewith . If there exists such thatand with satisfy (H3) and (H4) in (32), then

Lemma 5. Assume that (24) holds and satisfies Lemma 4. Then, there exist such that

Proof. By (27), we have such thatDue to the positivity of solutions in Lemma 3, the result is true by selecting small enough. The proof is complete.
By what we have done, we now give the proof of Theorem 1.

Proof. Definewith small enough. By direct analysis, we can verify that satisfy (H3) and (H4) for (6) if (24) holds (we present the calculation in Appendix). By Lemmas 1 and 5, we complete the proof.

4. Numerical Simulation

In this part, we simulate the dynamics of (6) when two patches and three species are involved by selecting some parameters, which may show the plentiful dynamics of such a delayed system without monotonicity. Since in (6) do not change the stability in our results, we shall not consider the delayed effect modeling the spatial movement among different patches.

Example 1. We first try to simulate the dynamics of a system satisfying Theorem 1. For the purpose, we consider the following delayed system:with and initial conditionEvidently, (39) satisfies (24) and has a unique positive steady state, of which the approximate value isFrom our main result, we see that the positive steady state is asymptotic stable, which is also illustrated in Figure 1.

Example 2. We now consider the possible extinction of some unknown functions when Theorem 1 does not hold and study the dynamics of the following delayed system:with and initial condition (40). Evidently, this system does not satisfy (24), and its dynamics may be different from our Theorem 1, which is presented in Figure 2.

Example 3. The weights among the same species partly play the role of irreducibility, which leads to the possible synchronism. What will happen if the condition does not hold? We consider the dynamics of the following delayed system:with and initial condition (40), and are positive.
In this system, the weights may be zero, and it does not satisfy our condition. By simulation in Figures 3 and 4, we may observe different persistence-extinction phenomena with partially degenerate weights.

Example 4. Finally, we try to explore the influence of delay size. In the literature, some nontrivial results of the delay size have been reported even in scalar equations [2]. We consider the dynamics of the following delayed system:with , and initial value conditionBy our simulations in Figure 5, we find that large delay may lead to the instability of the positive steady state, while the small delay may be harmless to the stability of the positive steady state.

5. Discussion

In this paper, we studied a delayed model on a weighted network, which may be regarded as a delayed competition model on patches. We established the global stability of the unique positive steady state by constructing contracting rectangles under condition (24). Moreover, we gave some numerical examples to model the dynamics of (6). From Examples 14, we found that the dynamics of this nonmonotone system may be very complex, which includes the stability of different steady states. In particular, if the weight coefficients among some patches are zero, the uniform convergence on different patches may be false even if (24) holds, and we observed this from Example 2. That is, weights among different patches may affect the dynamics of this system, and the role of weights should be further studied.

Moreover, if (24) does not hold, then the stability of the positive steady state may depend on the size of time delay, which is partially motivated by the dynamics of the classical delayed logistic model [27]. From Example 4, we observed different convergence results with different delays. To further show the effect of time delay in such a delayed model is an interesting question. However, with the increase in the number of unknown functions, the analysis of eigenvalues becomes complex, and this is a challenging question. We shall further consider these questions in our future works.

Finally, in many other mathematical models including uncertain discrete-time descriptor systems with multiple time-varying delays and parametric uncertainties [28], nonholonomic systems with time-varying delays and nonlinear disturbances [29], and uncertain chaotic systems with nonlinearity and time delay [30], it is important to estimate the dynamics in finite time. Due to the effect of time-varying delays, this is not a trivial work and deserves further investigation. Because the contracting rectangles could reflect certain uniform convergence, it is possible to consider the finite-time dynamics of these important models by contracting rectangles, at least when the initial data satisfy proper uniform conditions in the delayed interval.

Appendix

In this part, we verify the definition of contracting rectangles in (6) and (32). Since the verification of (32) can be obtained by that of (6) and can be found in [6], we only verify the result on (6).

We first verify the conclusion on in (H4); if , andfor some , thenby selecting small enough since (24) implieswhich also implies the admissibility of in the proof of Theorem 1.

We now verify the conclusion on in (H4); if , andfor some , thenby . Due to the arbitrary of , we finish the verification.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Authors’ Contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.