Abstract

In this paper, we establish two multi-dimensional discrete Halanay-type inequalities. As an application, we consider a malaria transmission model with two delays and employ the two discrete Halanay-type inequalities to establish the global (exponential) asymptotical stability of the disease free equilibrium of the model. It is further shown that the disease free equilibrium of the delayed model is globally asymptotically stable when the basic reproduction number of the model is less than 1.

Keywords

1 Introduction

The Halanay inequality was first introduced in [1] and many generalizations of this inequality have been obtained due to their significance in the analysis of delayed dynamical systems [2–9], and especially in proving the global exponential stability of the equilibria of mathematical models proposed in neural networks and biology with time delays [8–11]. The original form of the Halanay inequality and some of its generalizations are for continuous dynamical systems with only one unknown scalar function involved in the system [1, 6, 9], i.e., the one-dimensional case. In [9], Halanay inequality was generalized to a multi-dimensional form to deal with the stability of the equilibria of continuous dynamical systems in neural networks with impulses.

Along with the development of Halanay-type inequalities for continuous-time dynamical systems, discrete-time Halanay inequalities have also been established to handle the stability of discrete dynamical systems with delays [12–20]. For example, the following result was obtained in [15].

where\(\lambda_{0}\)is the smallest root in the interval\((0,1)\)of the equation

$$ \lambda^{r+1}+(a-1)\lambda^{r}-b=0. $$

(1.3)

Based on Theorem 1.1, different generalizations of discrete Halanay-type inequalities have been developed [12, 13, 16, 19]. However, the generalizations are mostly dealing with one-dimensional functions, such as the sequence \(\{x(n)\}\) in Theorem 1.1; while generalizations of discrete Halanay-type inequalities for multi-dimensional sequences \(\{X(n)\}=\{(x_{1}(n), x_{2}(n), \ldots, x_{m}(n))^{T}\}\) are hardly developed. Motivated by this, the first aim of our present work is to develop certain multi-dimensional discrete Halanay-type inequalities.

On the other hand, the global stability of the equilibria of epidemic models is a key research topic in the quantitative analysis of the transmission of infectious diseases. Much research on this topic has been done on compartmental models of infectious diseases such as influenza, malaria, dengue, cholera, etc. [21–25]. Usually, a compartmental infectious disease model has a disease free equilibrium, which is globally asymptotically stable when the basic reproduction number \(R_{0}\) of the model is less than 1, and it has a positive equilibrium which is globally asymptotically stable when \(R_{0}>1\) [26]. Time delays are frequently involved in these infectious disease models [21–23, 27] and their presence poses great difficulty in stability analysis, especially in analyzing the global asymptotical stability of the equilibria of these models. For infectious disease models with two or more delays involved, the local asymptotical stability of the equilibrium can be obtained through the analysis of the eigenvalues of the Jacobian matrix of the linearized model at the corresponding equilibrium, but the global asymptotical stability of the equilibria of some of these infectious disease models with several delays remains unsolved at present [22, 23]. Hence, another aim of this work is to establish the global asymptotical stability of the disease free equilibrium of a malaria transmission model with two delays. We shall tackle this by applying the multi-dimensional discrete Halanay inequalities developed in this paper. Moreover, we shall show that the disease free equilibrium of the model is globally asymptotically stable when the basic reproduction number \(R_{0}\) is less than 1, which is well consistent with the threshold property of the basic reproduction number.

The paper is organized as follows. In Section 2, we present two generalizations of Theorem 1.1 to multi-dimensional case. Applications of these two generalizations are given in Section 3 to obtain the global asymptotical stability of the disease free equilibrium of a malaria transmission model with two delays. Finally, some concluding remarks are given in Section 4.

2 Generalized discrete Halanay-type inequalities

In this section, we shall establish two generalizations of Theorem 1.1 to multi-dimensional case. First, we introduce some notations as follows. Throughout, we denote

Consequently, there exists \(\lambda_{i}\in(0, 1)\) such that \(F_{i}(\lambda_{i})=0\), i.e., \(\lambda_{i}\) is a root of equation (2.5).

It is also noted that \(F_{i}(\lambda)=0\) is the characteristic equation of (2.10) when \(0<\lambda<1\). Hence, \(y_{i}(n)=K\lambda_{i}^{n}\), \(n\ge0\) is a solution of (2.10), where K is a constant. In particular, we set \(K=\max_{s\in\mathbb{S}_{1}}\{0, x_{i}^{(s)}\}\) such that (2.8) is fulfilled for \(s\in\mathbb{S}_{1}\). Thus, it follows from (2.11) that

Proof

Set \(k_{i}=\max_{s\in\mathbb{S}_{1}} \{0, x_{i}^{(s)}\}\), \(i\in\mathbb{I}\) and \(\lambda=\max\{\lambda_{1}, \lambda_{2}, \ldots, \lambda_{m}\}\), where \(\lambda_{i}\) has the same meaning as in Theorem 2.1. It is immediate from Theorem 2.1 that \(X(n)\le\mathcal{K}\lambda^{n}\), \(n\ge0\). □

Proof

Since \(C=-(A+B)\) is a nonsingular M-matrix, from (2) of Lemma 2.4 there exists \(z>0\) such that \(-(A+B)z>0\), or equivalently,

$$ (A+B)z< 0. $$

(2.18)

Noting that \(A\lambda^{r}+\lambda^{r}E+B-\lambda^{r+1}E=A+B\) when \(\lambda=1\), it follows from (2.18) that (2.16) has a solution \(z>0\) when \(\lambda=1\). Consequently, (2.16) has a solution \(z>0\) when \(0<\lambda<1\) due to the continuity.

where \(\tau(n)=(\tau_{ij})_{m\times m}\), \(0\le\tau_{ij}\le\tau\), \(i,j\in\mathbb{I}\), τ is a nonnegative integer and the constant \(\mu>0\), the vector \(z>0\) are determined by

$$\bigl[\exp(\mu) \bigl(P+Q\exp(\mu\tau)\bigr)-E \bigr]z\le0, $$

then \(U(n)\le z\exp(-\mu(n-n_{0}))\), \(n\ge n_{0}\).

Comparing (2.23) with (2.1), we see that \(P=A+E\) and \(Q=B\). If we set \(C=-(A+B)\), \(D=E\), and \(M=E-(-(A+B))\), then \(D^{-1}M=M=E+A+B=P+Q\). Hence, the condition \(\rho(P+Q)<1\) is equivalent to the condition that \(-(A+B)\) is a nonsingular M-matrix (according to (3) of Lemma 2.4).

The next theorem shows that the estimation of \(\{X(n)\}\) similar to (2.17) still holds without the assumption (2.15), which imposes a condition on the initial values of (2.1).

it is easy to verify that (2.26) holds for \(\omega=\omega_{1}\omega_{2}\). Noting Remark 2.5, \(\omega z>0\) and \(\lambda\in(0, 1)\) also fulfill (2.16). Thus, (2.25) is obtained as in the proof of Theorem 2.6 by setting \(u_{i}(n)\) in (2.21) as \(u_{i}(n)=(1+\epsilon)\omega z_{i}\lambda^{n}\). □

Remark 2.10

The results of Theorems 2.1, 2.6, and 2.9 can be applied to obtain the global exponential stability of the zero solution of (2.1). As such, we can make use of these theorems to analyze the global exponential stability [27, 30] of the equilibrium of a discrete dynamical system with time delays. It is well known that the equilibrium is globally asymptotically stable if it is globally exponentially stable [30]. Hence, the global asymptotical stability of the equilibria of some discrete dynamical systems with time delays may be established via the results obtained in this section. We shall give such applications in the next section.

3 Global asymptotical stability of equilibrium

In this section, we shall obtain the global stability of the disease free equilibrium of an infectious disease model with two delays that describes malaria transmission. Our proofs employ the generalizations of Halanay-type inequalities obtained in Section 2.

In [22], the following delayed Ross-Macdonald model was established to describe the malaria transmission between human and mosquito:

In model (3.1), \(x(y)\) is the ratio of the number of infected human (mosquito) to the total number of human (mosquito) population at time t (days). The total number of human (mosquito) population is H\((M)\) and H\((M)\) is a constant. Define \(m=M/H\). Moreover, γ is the average recovery rate of human from malaria infection, μ is the natural death rate of mosquito, a is the average number of bites of a mosquito per human per day, b is the rate of a susceptible human becoming infectious after being bitten by an infected mosquito, c is the rate of a susceptible mosquito becoming infectious after it bites an infected human, and \(\tau_{1}\) and \(\tau_{2}\) are the average incubation times of the parasites in the body of human and mosquito, respectively. For more details of (3.1), one can refer to [22].

In model (3.2), \(n\ge0\) and \(x(n)\) (\(y(n)\)) is the ratio of the number of infected humans (mosquitos) to the total number of humans (mosquitos) of the population on the nth day. Moreover, \(r_{1}\) and \(r_{2}\) are now positive integers to reflect the incubation time of the parasites in the body of human and mosquito, respectively. The implications of the other parameters are the same as in (3.1) and a, b, c, m, γ, and μ are all positive constants.

In view of the practical background of model (3.2), it is natural to assume that

$$ 0< x^{(s)}< 1,\qquad 0< y^{(s)}< 1,\quad s=-r,-r+1, \dots,0. $$

(3.3)

Remark 3.1

The basic reproduction number\(R_{0}\) of model (3.2) (refer to [22, 23]) is defined as

$$ R_{0}=\frac{a^{2}bcm\alpha\beta}{\gamma\mu}, $$

(3.4)

where \(\alpha=\exp(-\gamma r_{1})\) and \(\beta=\exp(-\mu r_{2})\). If we denote \(R_{01}=abm\alpha/\gamma\) and \(R_{02}=ac\beta/\mu\), then \(R_{0}=R_{01}R_{02}\).

Remark 3.2

Since γ is the recovery rate of human from malaria infection, it follows that \(1/\gamma\) is the average recovery time of human from the infection. As it takes 7-30 days for a human to recover from the infection [22, 23], it is reasonable to assume that

$$ 0< \gamma< 1. $$

(3.5)

Moreover, since the natural lifespan of mosquito is 14-60 days [22, 23], it also makes sense to assume that

Similarly, we can prove that \(0< y(n)<1\) for \(0< n\le r\). By assuming that \(0< x(n)<1\) and \(0< y(n)<1\) for \((k-1)r< n\le kr\), using a similar technique we can show that \(0< x(n)<1\) and \(0< y(n)<1\) for \(kr< n\le (k+1)r\). Hence, it is shown by induction that \(0< x(n)<1\) and \(0< y(n)<1\) for \(n\ge0\). □

Remark 3.4

Lemma 3.3 gives sufficient conditions to guarantee that the solutions of (3.2) with (3.3) are positive and bounded. As pointed out in [31], sometimes these sufficient conditions may be quite rigorous. Indeed, noting Remark 3.1 we see that condition (3.7) implies \(R_{0}<1\), but in fact numerical simulations show that solutions of (3.2) with (3.3) can be positive and bounded even when \(R_{0}>1\). In view of this, we shall introduce a set Ω of parameters that ensures the positivity and boundedness of the solutions of (3.2) with (3.3), let

for \(n\ge0\), where \(k_{1}=\max_{s\in\mathbb{S}_{1}}\{x^{(s)}\}\) and \(k_{2}=\max_{s\in\mathbb{S}_{1}}\{y^{(s)}\}\).

In view of (3.12), the zero solution of (3.2), which is the disease free equilibrium \(E_{0}\), is globally exponentially stable. Noting Remark 2.10, it follows that \(E_{0}\) is also globally asymptotically stable. □

Remark 3.7

When the conditions of Theorem 3.6 are satisfied, noting (3.12) the exponential convergence rate of the solutions of (3.2) with (3.3) to the disease free equilibrium is \(\min\{-\ln\lambda_{1}, -\ln\lambda _{2}\}\).

Remark 3.8

Condition (3.7) implies that \(R_{0}=R_{01}R_{02}<1\), but \(R_{0}<1\) may not imply (3.7). In the literature of compartmental infectious disease models, usually \(E_{0}\) is globally asymptotically stable when \(R_{0}<1\). This expected result is not obtained when Theorem 2.1 is applied to model (3.2). As such, in the sequel we shall apply Theorem 2.9 to model (3.2) to see whether \(E_{0}\) is globally asymptotically stable when \(R_{0}<1\).

we have \(\gamma>0\) and \(\det(-(A+B))=\gamma\mu-a^{2}bcm\alpha\beta=\gamma\mu(1-R_{0})>0\) due to (3.13), i.e., the successive principal minors of \(-(A+B)\) are both positive. Therefore, according to (1) of Lemma 2.4, \(-(A+B)\) is a nonsingular M-matrix. The conditions of Theorem 2.9 are all satisfied, hence the disease free equilibrium \(E_{0}\) of (3.2) with (3.3) is globally exponentially stable and consequently it is globally asymptotically stable in Ω. □

We shall now give an example to illustrate Theorem 3.9 and Remark 3.4.

By direct computation, we have \(R_{01}=0.2624\), \(R_{02}=1.0570\), and \(R_{0}=0.2774\). According to Theorem 3.9, the disease free equilibrium \(E_{0}\) of (3.14) is globally asymptotically stable since (3.5), (3.6), and (3.13) are satisfied. Indeed, the dynamics of Example 3.10 is depicted in Figure 1 and we observe that both \(x(n)\) and \(y(n)\) tend to 0 as \(n\to+\infty\), i.e., the disease free equilibrium \(E_{0}\) is globally asymptotically stable.

Figure 1

The dynamics of Example 3.10. The solid line is the variation of the infected human population density and the dashed line is the variation of the infected mosquito population.

In this example, we note that (3.7) is not satisfied since \(R_{02}>1\). However, it is observed from Figure 1 that the solution of (3.14) is positive and bounded. This illustrates Remark 3.4 and shows that even if (3.7) does not hold, the set Ω (see (3.8)) may not be empty.

4 Concluding remarks

Halanay-type inequalities, whether continuous or discrete, have been widely applied to obtain the global exponential (asymptotical) stability of the equilibria of dynamical systems with several delays, especially dynamical systems of neural networks. In this paper, we have derived two generalizations of multi-dimensional discrete Halanay-type inequalities. Further, the generalizations are applied to a discrete malaria transmission model with two delays. We have shown that the disease free equilibrium is globally asymptotically stable when the basic reproduction number \(R_{0}\) is less than 1, which is well consistent with the threshold property of the basic reproduction number.

The global asymptotical stability of the equilibria of infectious disease models with time delays is usually obtained via the construction of suitable Lyapunov functionals together with Razumikhin-type theorem and/or LaSalle invariant sets. However, a suitable Lyapunov functional is somewhat difficult to construct for a delayed dynamical system. Hence, it is reasonable to try other methods to obtain the global asymptotical stability of equilibria of dynamical systems with time delays. From our present work, we have observed that it is direct and simple to obtain the global asymptotical stability of equilibria of dynamical systems with time delays via Halanay-type inequalities.

When using Halanay-type inequalities established in this paper to obtain the global asymptotical stability of the equilibrium of a discrete dynamical system with time delays, the positivity of the solutions should be initially guaranteed. This can be proved for many kinds of compartmental infectious disease models. Hence, it is direct to obtain the global asymptotical stability of the disease free equilibria of dynamical systems with time delays by applying Halanay-type inequalities. It is well known that, in order to obtain the global asymptotical stability of the positive equilibrium \((x^{*}, y^{*})\) of a dynamical system with time delays, the change of variables \(\bar{x}=x-x^{*}\) and \(\bar{y}=y-y^{*}\) is usually applied to transfer the global asymptotical stability of the positive equilibrium to the global asymptotical stability of the zero solution of the system with respect to x̄ and ȳ. However, after the change of variables, the positivity of x̄ and ȳ cannot be guaranteed. Hence, if one intends to employ Halanay-type inequalities to obtain the global asymptotical stability of the positive equilibrium, new techniques are needed to deal with x̄ and ȳ. This remains as future work.

Declarations

Acknowledgements

The research of Chunqing Wu is sponsored by the Jiangsu Overseas Research & Training Program. The authors would like to thank the referees for their comments which help to improve the paper.

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Competing interests

None of the authors have any competing interests to the paper.

Authors’ contributions

All authors contributed equally. All authors read and approved the final manuscript.

Authors’ Affiliations

(1)

School of Mathematics and Physics, Changzhou University

(2)

School of Electrical and Electronic Engineering, Nanyang Technological University