## Advanced Nonlinear Dynamics of Population Biology and Epidemiology

View this Special IssueResearch Article | Open Access

# Hopf Bifurcation of a Delayed Epidemic Model with Information Variable and Limited Medical Resources

**Academic Editor:**Kaifa Wang

#### Abstract

We consider SIR epidemic model in which population growth is subject to logistic growth in absence of disease. We get the condition for Hopf bifurcation of a delayed epidemic model with information variable and limited medical resources. By analyzing the corresponding characteristic equations, the local stability of an endemic equilibrium and a disease-free equilibrium is discussed. If the basic reproduction ratio , we discuss the global asymptotical stability of the disease-free equilibrium by constructing a Lyapunov functional. If , we obtain sufficient conditions under which the endemic equilibrium of system is locally asymptotically stable. And we also have discussed the stability and direction of Hopf bifurcations. Numerical simulations are carried out to explain the mathematical conclusions.

#### 1. Introduction

From an epidemiological viewpoint, it is important to investigate the global dynamics of the disease transmission. In the literature, many authors have researched various epidemic models [1, 2], in which the stability analyses have been carried out extensively. In the recent years, based on SIR epidemic model, in order to investigate the spread of an infectious disease transmitted by a vector, Wang et al. [3] have considered the asymptotic behavior of the following delayed SIR epidemic model:

Since nonlinearity in the incidence rates has been observed in disease transmission dynamics, it has been suggested that the standard bilinear incidence rate will be modified into a nonlinear incidence rate by many authors [4–6]. In [7], incidence rate in (1) was replaced by a nonlinear incidence rate of the form with the following system:

In order to control the spread of epidemic, we consider the new variable : called information variable which summarizes information about the current state of the disease, that is, depending on current values of state variables, and also summarizes information about past values of state variables. Many authors have used this variable in their models (see, e.g., [8–10]).

In this paper, we consider the information variable , nonlinear incidence rate of the form , and limited medical resources . The model can be described by the following system of equations: where , , and , , denote the numbers of susceptible, infective, and recovered individuals at time , respectively. is the intrinsic growth rate of susceptibles, is the carrying capacity of susceptibles, is the saturation factor that measures the inhibitory effect, is the transmission or contact rate, , are the natural death rate of the infective and recovered individuals, is the natural recovery rate, is the disease-related mortality, is the maximal medical resources supplied per unit time, and is half-saturation constant. , , , , , , , , are all positive.

We further assume that the function is continuous on and continuously differentiable on satisfying the following hypotheses:(1) is strictly monotone increasing on with ;(2) is monotone increasing on with .

The organization of this paper is as follows. In Section 2, we explore the existence of disease-free equilibria point and the unique existence of the endemic equilibrium point. In Section 3, we analyze the stability of the disease-free equilibria. In Section 4, we obtain sufficient conditions under which the endemic equilibrium of system is locally asymptotically stable. In Section 5, we also have discussed the stability and direction of Hopf bifurcations. A numerical analysis and a simple discussion are given to conclude this paper in Section 6.

#### 2. The Existence of Equilibria

The nonlinear integrodifferential system (4) can be transformed into the following set of nonlinear ordinary differential questions:

Since the dynamical behavior of the last equation of the system (5), that is, the dynamics of , depends only the dynamics of , we do not consider that equation in our discussion. Here we will study the following nonlinear ordination differential equations: For simplicity, we nondimensionalize system (6) by defining

We note that also satisfies the hypotheses (1) and (2). Dropping the for convenience of readers, system (6) can be written in the following form: The basic reproduction is .

Theorem 1. *(1) The system (8) has a trivial equilibrium and the disease-free equilibrium .**(2) If , the system (8) has one endemic equilibrium except the disease-free equilibria and .*

* Proof. *(1) Let ; we have , or ; it is not easy to find that the system has a trivial equilibrium and the disease-free equilibria and .

(2) If , from the third question of (8), we have ; from the second question of (8), we have

Then substituting them into the first question of (8) yields
Let . By hypothesis (2), we obtain

Since is strictly monotone decreasing function on , it suffices to show that holds for sufficiently large. From (1), is either unbounded above or bounded above on .

First, we suppose that is unbounded above. Then there exists an such that , from which we have for all . Second we suppose that is bounded above. Then, from (2), is unbounded above on ; that is, there exists an such that . This yields for all . Therefore, for the both cases, there exists a unique endemic such that . By the second and third equations of (8), there exists a unique endemic equilibrium of system (8) if .

Second, we assume ; then it is obvious that system (8) has no equilibria. Hence the proof is complete.

#### 3. The Stability Analysis of Disease-Free Equilibrium Point

In this section, we will examine the local stability of the equilibria by analyzing the eigenvalues of the Jacobian matrices of (8) at the equilibria and using Routh-Hurwitz criterion.

Let be the arbitrarily equilibrium point of system (8); then the Jacobian matrix of (8) at is

Then the characteristic equation of the system (8) at equilibrium is

Theorem 2. *The trivial equilibrium of system (8) is always unstable.*

* Proof. *The characteristic equation (13) at becomes as follows:
Since (14) has a positive root , is unstable.

Theorem 3. *The disease-free equilibrium of system (8) is locally asymptotically stable if and it is unstable if .*

* Proof. *For , the characteristic equation (13) at becomes as follows:
It is clear that both and are all the negative root of (15). Then the other root of (15) is determined as the following equation:
For the case , we suppose on the contrary that is not locally asymptotically stable; that is, . Then, there exists a root , such that . However, from (16), we obtain

which is a contradiction. Hence, if , the disease-free equilibrium of system (8) is locally asymptotically stable.

Now, we put

For the case , we have and ; then has at least one positive root. Hence, is unstable if and only if . The proof is complete.

#### 4. The Stability Analysis of the Endemic Equilibrium Point

Theorem 4. *If , , and , then the positive equilibrium of system (8) is locally asymptotically stable.*

* Proof. *The characteristic equation of (13) at becomes as follows:
The above equation can be rewritten as
where , , and

Let .

Then if , (20) becomes ; that is,
where

Let ; we have , , and

By using the Routh-Hurwitz theorem, has negative real part for . So the positive equilibrium is locally asymptotically stable.

In the following, we investigate the existence of purely imaginary roots to (19). Equation (19) takes the form of a third-degree exponential polynomial in , with all the coefficients of and depending on . Beretta and Kuang [11] established a geometrical criterion which gives the existence of purely imaginary root of a characteristic equation with delay dependent coefficients.

Now we let be a root of (20) from which we have that

Hence, we have that

From (26), it follows thatBy the definitions of , is as in (20), and applying the property (1), (27a) and (27b) can be written as

which yields .

Assume that is the set where is a positive root of From

we have , where , , , and, for , is not defined. Then, for all in , satisfied .

Let ; then we have that

Assume that has only one positive real root; we denote by this positive real root. Thus, (29) has only one positive real root . And the critical values of and are impossible to solve explicitly, so we will use the procedure described in Beretta and Kuang [11] and Song et al. [12]. According to this procedure, we define such that and are given by the right-hand sides of (27a) and (27b), respectively, with given by (19).

And the relation between the argument and in (28) for must be Hence we can define the maps given by where a positive root of (31) exists in . Let us introduce the functions which are continuous and differentiable in . Thus, we give the following theorem which is due to Beretta and Kuang [11].

Theorem 5. *Assume that is a positive root of (19) defined for , , and, at some , for some . Then a pair of simple conjugate pure imaginary roots exists at which crosses the imaginary axis from left to right if
*

Applying Theorems 2 and 3 and the Hopf bifurcation theorem for functional differential equation [13], we can conjecture the existence of a Hopf bifurcation as stated in Theorem 6.

Theorem 6 (a conjecture). *For system (2), there exists , such that the equilibrium is asymptotically stable for , and it becomes unstable for staying in some tight neighborhood of , with a Hopf bifurcation occurring when .*

#### 5. Stability and Direction of Hopf Bifurcations

In this section, we will study the direction of the Hopf bifurcation and stability of bifurcating periodic solutions by using the normal theory and center manifold theorem due to Hassard et al. [14]. Letting , , , , , and dropping the bars for simplification of notations, system (8) becomes functional differential equations in as Then system (36) is equivalent to where and , are given, respectively, by By the Riesz representation theorem, there exists a function of bounded variation for , such that In fact, we can choose where denote the Dirac delta function. For , define

Then system (37) is equivalent to where for .

For , define and a bilinear inner product where . Then and are adjoint operators. By the discussion in (20), we know that are eigenvalues of . Hence, they are also eigenvalues of . We first need to compute the eigenvectors of and corresponding to and , respectively.

Suppose that is the eigenvectors of corresponding to ; then . Then, from the definition of and (38), (40), (41), and , we have We obtain Similarly, we can obtain the eigenvector of corresponding to , where In order to assure that , , we need to determine the value of . By (45), we have Therefore, we can choose as We use the way of [14] and similarly way of [3]; we obtain that the coefficients are where

Besides, , are satisfied with the following equation:

Thus, we can determine and from (52). Furthermore, we can compute by (51). Thus we can compute the following values: By the result of Hassard et al. [14], we have the following.

Theorem 7. *In (54), the following results hold:*(i)*the sign of determines the directions of the Hopf bifurcation: if , then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for ;*(ii)*the sign of determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if ;*(iii)*the sign of determines the period of the bifurcating periodic solutions: the period is increasing (decreasing) if .*

#### 6. Numerical Simulations

To demonstrate the theoretical results obtained from this paper, letting , we will give some numerical simulations. We consider the hypothetical set of parameter values as follows.(1)Consider ; ; ; ; ; ; ; . By directly computing, we obtain . According to Theorem 4, we know that the disease-free equilibrium of system (8) is locally asymptotically stable for this case (see Figures 1(a)–1(d)).(2)Consider ; ; ; ; ; ; ; ; . By directly computing, we obtain . According to Theorem 4, we know that the disease-free equilibrium of system (8) is locally asymptotically stable for this case (see Figures 2(a)–2(d)).(3)Consider ; ; ; ; ; ; ; , . By directly computing, we obtain . According to Theorem 6, we know that the disease-free equilibrium of system (8) is unstable for this case (see Figures 3(a)–3(d)).

**(a) Time series of of the system**

**(b) Time series of of the system**

**(c) Time series of of the system**

**(d)**

**(a) Time series of of the system**

**(b) Time series of of the system**

**(c) Time series of of the system**

**(d)**

**(a) Time series of of the system**

**(b) Time series of of the system**

**(c) Time series of of the system**

**(d)**

#### 7. Conclusion

In this paper, we formulate and analyze a new delayed epidemic model with information variable and limited medical resources, the conditions for Hopf bifurcation to occur are derived. By analyzing the model, we have found the existence of disease-free equilibria and and have a unique positive equilibrium . The basic reproduction number changes the stability of the disease-free equilibrium. When , we discuss the stability of the disease-free equilibrium by analyzing the corresponding characteristic equations and constructing a Lyapunov functional, respectively. The conclusion reveals that the disease dies out and when , we also get the sufficient criteria of stability switch at the positive equilibrium. Using the time delay (i.e., incubation time) as a bifurcation parameter, the local stability of the endemic equilibrium is investigated, and the conditions for Hopf bifurcation to occur are derived. Using the normal form theory and the center manifold theorem introduced by Hassard et al., we have studied the direction and stability of the bifurcating periodic solutions. Our theoretical results show that the time delay must be responsible for the observed regular cycles of disease incidence.

Lastly, a numerical simulation provided that when is less than 1, the disease-free equilibrium is stable and while is more than 1, the disease-free equilibrium is unstable; that is, the endemic equilibrium exists (see Figure 1). Further, for , there will exist , such that the endemic equilibrium is asymptotically stable for (see Figure 2) and becomes unstable for staying in some right neighborhood of , with a Hopf bifurcation occurring when . If , the endemic equilibrium is unstable (see Figure 3).

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors would like to thank the anonymous referees for their careful reading of the original paper and their many valuable comments and suggestions that greatly improved the presentation of this work. This work is supported by the Natural Science Foundation of Shanxi province (2013011002-2).

#### References

- J. W. Jia and Q. Y. Li, “Qualitative analysis of an SIR epidemic model with stage structure,”
*Applied Mathematics and Computation*, vol. 193, no. 1, pp. 106–115, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J.-Q. Li, Z.-E. Ma, and J. Zhang, “Global analysis of some epidemic models with general contact rate and constant immigration,”
*Applied Mathematics and Mechanics*, vol. 25, no. 4, pp. 396–404, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. J. Wang, J. Z. Zhang, and Z. Jin, “Analysis of an SIR model with bilinear incidence rate,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 4, pp. 2390–2402, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - G. Huang and Y. Takeuchi, “Global analysis on delay epidemiological dynamic models with nonlinear incidence,”
*Journal of Mathematical Biology*, vol. 63, no. 1, pp. 125–139, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. Korobeinikov, “Global properties of infectious disease models with nonlinear incidence,”
*Bulletin of Mathematical Biology*, vol. 69, no. 6, pp. 1871–1886, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - R. Xu and Z. E. Ma, “Global stability of a SIR epidemic model with nonlinear incidence rate and time delay,”
*Nonlinear Analysis: Real World Applications*, vol. 10, no. 5, pp. 3175–3189, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. Enatsu, E. Messina, Y. Muroya, Y. Nakata, E. Russo, and A. Vecchio, “Stability analysis of delayed SIR epidemic models with a class of nonlinear incidence rates,”
*Applied Mathematics and Computation*, vol. 218, no. 9, pp. 5327–5336, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. d'Onofrio, P. Manfredi, and P. Manfredi, “Bifurcation thresholds in an SIR model with information-dependent vaccination,”
*Mathematical Modelling of Natural Phenomena*, vol. 2, no. 1, pp. 26–43, 2007. View at: Publisher Site | Google Scholar | MathSciNet - A. d'Onofrio, P. Manfredi, and E. Salinelli, “Vaccinating behaviour, information, and the dynamics of SIR vaccine preventable diseases,”
*Theoretical Population Biology*, vol. 71, no. 3, pp. 301–317, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH - B. Buonomo, A. d'Onofrio, and D. Lacitignola, “Global stability of an SIR epidemic model with information dependent vaccination,”
*Mathematical Biosciences*, vol. 216, no. 1, pp. 9–16, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - E. Beretta and Y. Kuang, “Geometric stability switch criteria in delay differential systems with delay dependent parameters,”
*SIAM Journal on Mathematical Analysis*, vol. 33, no. 5, pp. 1144–1165, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Y. Song, S. L. Wang, and J. Dong, “Stability properties and Hopf bifurcation of a delayed viral infection model with lytic immune response,”
*Journal of Mathematical Analysis and Applications*, vol. 373, no. 2, pp. 345–355, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Hale and S. M. V. Lunel,
*Introduction to the Theory of Functional Differential Equations Methods and Applications*, Spring, 1993. - B. Hassard, N. Kazarinoff, and Y. Wan,
*Theory and Applications of Hopf Bifurcation*, Cambridge University Press, Cambridge, UK, 1981. View at: MathSciNet

#### Copyright

Copyright © 2014 Caijuan Yan and Jianwen Jia. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.