Research Article | Open Access

# Bifurcation of a Delayed SEIS Epidemic Model with a Changing Delitescence and Nonlinear Incidence Rate

**Academic Editor:**Lu-Xing Yang

#### Abstract

This paper is concerned with a delayed SEIS (Susceptible-Exposed-Infectious-Susceptible) epidemic model with a changing delitescence and nonlinear incidence rate. First of all, local stability of the endemic equilibrium and the existence of a Hopf bifurcation are studied by choosing the time delay as the bifurcation parameter. Directly afterwards, properties of the Hopf bifurcation are determined based on the normal form theory and the center manifold theorem. At last, numerical simulations are carried out to illustrate the obtained theoretical results.

#### 1. Introduction

The outbreak of infectious diseases had not only caused the loss of billions of lives but also badly damaged the social economy in a short time, which brought much pain to human society [1]. Thus, it has been an increasingly urgent issue to understand how to prevent or slow down the transmission of infectious diseases. To this end, many mathematical models have been proposed for describing the spread process of infectious diseases [2–10]. However, all the epidemic models above do not consider the change of delitescence of the infectious diseases. Considering that the diversity of the delitescence period in each infected individual who is infected with disease virus is mainly due to the variation of the virus and the distinct constitution of different people for some disease, such as H1N1 disease, Wang proposed the following SEIS epidemic model with a changing delitescence and a nonlinear incidence rate [11]:where , , and denote the numbers of the susceptible, exposed, and infectious populations at time , respectively. is the recruitment rate of the susceptible population; is the natural death rate of the population; is the death rate due to the disease of the infected population; is the rate at which the exposed population becomes infectious; is the rate at which the infected population returns to the susceptible population because of the treatment; is the rate at which the infected population becomes the exposed one; and is the rate at which the infected population becomes infectious directly. is the nonlinear incidence rate, where measures the infection force of the disease and measures the inhibition effect from the behavioral change of the susceptible population. Wang investigated global stability of system (1).

In fact, many infectious diseases have different kinds of delays during their spreading process in the population, such as latent period delay [9, 12–16], immunity period delay [17, 18], and infection period delay [19]. The time delay may induce Hopf bifurcation and periodic solutions. The occurrence of a Hopf bifurcation means that the state of the epidemic disease prevalence changes from an equilibrium to a limit cycle. Therefore, the time delay can influence the dynamics of infectious diseases. So it is necessary and useful to investigate system (1) with time delay. Based on this fact and taking the period used to cure the infectious population, we consider the following delayed epidemic system:where is the time delay due to the period that is used to cure the infectious population. That is, we assume that all the infectious populations will survive after time . The initial conditions for system (2) are where .

The outline of this paper is as follows. In the next section, stability of the endemic equilibrium is analyzed and the critical value of the time delay at which a Hopf bifurcation occurs is obtained. In Section 3, direction and stability of the Hopf bifurcation are investigated. In Section 4, the obtained theoretical results are verified by some numerical simulations. Finally, this work is summarized in Section 5.

#### 2. Stability of the Endemic Equilibrium and Existence of Hopf Bifurcation

By a direct computation, we know that if (I) and , (II) and , (III) , and , or (IV) and , then system (2) has a unique endemic equilibrium , whereand is the unique positive root of the following equation:where

Let , , . We can rewrite system (2) as the following form:whereThen we obtain the linearized system of system (2)The characteristic equation iswhereWhen , (10) reduces to

Routh-Hurwitz criterion implies that is locally asymptotically stable without delay if condition holds.

() , .

For , substituting () into (10), we obtainThenwhereLet ; thenwhere . According to the analysis about the distribution of roots of (16) in Song et al. [20], we have the following result.

Lemma 1. *For the polynomial equation (16),*(1)*if , then (16) has at least one positive root;*(2)*if and , then (16) has no positive roots;*(3)*if and , then (16) has positive roots if and only if and .**Next, we assume that the coefficients in (16) satisfy the following condition.** (i) or (ii) , , , and .**Thus, (14) has at least one positive root such that (10) has a pair of purely imaginary roots . The corresponding critical value can be obtained from (13)Taking derivative with respect to on both sides of (10), we obtainFurther, we have**Thus, if the condition : holds, then . Then, based on the Hopf bifurcation theorem in [21], we have the following.*

Theorem 2. *For system (2), if the conditions hold, then the endemic equilibrium of system (2) is asymptotically stable for and system (2) undergoes a Hopf bifurcation at the endemic equilibrium when , where is defined in (17).*

#### 3. Direction and Stability of the Hopf Bifurcation

Let , ; then is the Hopf bifurcation value of system (2). Rescaling the time delay , then system (2) can be transformed into an FDE in as follows:wherewhere , , and are defined by Appendix A.

By the Riesz representation theorem, there exists a matrix function , , whose components are of bounded variation, such thatIn fact, we chooseFor , we defineThen system (20) is equivalent toFor , the adjoint operator of is defined asand a bilinear inner product is defined bywhere .

Let be the eigenvector of belonging to and be the eigenvector of belonging to . By a direct computation, we can getFrom (27), we can getThen we choosesuch that .

Next, we can obtain the coefficients , , , and by using the method introduced in [21] and a computation process similar to that in [22–24]. The expressions of , , , and are defined by Appendix B.

Then, we can get the following coefficients which determine the properties of the Hopf bifurcation:In conclusion, we have the following results.

Theorem 3. *For system (2), if , then the Hopf bifurcation is supercritical (subcritical). If , then the bifurcating periodic solutions are stable (unstable). If , then the bifurcating periodic solutions increase (decrease).*

#### 4. Numerical Simulations

In order to verify the efficiency of the obtained results in the paper, we carry out some numerical simulations in this section. By extracting some values from [11] and considering the conditions for the existence of the Hopf bifurcation, we consider the special case of system (2) with the parameters , , , , , , , and , Then, system (2) becomes the following form:from which we can obtain the unique positive root and then we get the unique endemic equilibrium . Then, we can obtain , , and . Thus, based on Theorem 2, we know that the endemic equilibrium is locally asymptotically stable when , which can be illustrated by Figures 1 and 2. In this case, the disease can be controlled easily. Once the value of the delay passes through the critical value , then the endemic equilibrium loses its stability and a Hopf bifurcation occurs, and a family of periodic solutions bifurcate from the endemic equilibrium . This property can be shown as in Figures 3 and 4. In this case, the disease will be out of control.

In addition, according to (31), we get , , , and . Therefore, we can conclude that the Hopf bifurcation is supercritical and the bifurcating periodic solutions are stable and increase. Since the bifurcating periodic solutions are stable, it can be concluded that the populations in system (32) can coexist from the view of ecology. Based on this fact, we can conclude that the time delay is harmful for system (32).

#### 5. Conclusions

We generalize a delayed SEIS (Susceptible-Exposed-Infectious-Susceptible) epidemic model with a changing delitescence and nonlinear incidence rate in this paper by introducing the time delay due to the period that is used to cure the infectious population into the SEIS model considered in the literature [11]. Compared with the literature [11], we mainly consider the effect of the time delay on the model.

The main results are given in terms of local stability and Hopf bifurcation. Stability of the endemic equilibrium is investigated by analyzing the corresponding characteristic equation. By choosing the time delay as a bifurcation parameter, sufficient conditions have been established for local existence of Hopf bifurcation at the endemic equilibrium. Then, with the help of the normal form theory and the center manifold theorem due to Hassard et al. [21], direction and stability of the Hopf bifurcation are determined. Finally, through numerical simulations, it can be concluded that the period used to cure the infectious population plays an important role in the disease spreading and the disease may be controlled by shortening the period used to cure the infectious population.

#### Appendix

#### A. The Expressions of , , and

#### B. The Expressions of , , , and

#### Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported by Natural Science Foundation of the Higher Education Institutions of Anhui Province (KJ2015A144).

#### References

- T. Zhang, X. Meng, T. Zhang, and Y. Song, “Global dynamics for a new high-dimensional {SIR} model with distributed delay,”
*Applied Mathematics and Computation*, vol. 218, no. 24, pp. 11806–11819, 2012. View at: Publisher Site | Google Scholar | MathSciNet - Y. Xiao, T. Zhao, and S. Tang, “Dynamics of an infectious diseases with media/psychology induced non-smooth incidence,”
*Mathematical Biosciences and Engineering. MBE*, vol. 10, no. 2, pp. 445–461, 2013. View at: Publisher Site | Google Scholar | MathSciNet - W. Qin, S. Tang, and R. A. Cheke, “Nonlinear pulse vaccination in an SIR epidemic model with resource limitation,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 670263, 13 pages, 2013. View at: Publisher Site | Google Scholar - T. Zhang, R. Kang, K. Wang, and J. Liu, “Global dynamics of an {SEIR} epidemic model with discontinuous treatment,”
*Advances in Difference Equations*, vol. 2015: Article 361, 16 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - X. Wang, “An SIRS Epidemic Model with Vital Dynamics and a Ratio-Dependent Saturation Incidence Rate,”
*Discrete Dynamics in Nature and Society*, vol. 2015, Article ID 720682, 9 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - Y. Tian and X. Liu, “Global dynamics of a virus dynamical model with general incidence rate and cure rate,”
*Nonlinear Analysis. Real World Applications*, vol. 16, pp. 17–26, 2014. View at: Publisher Site | Google Scholar | MathSciNet - D. Zhao and S. Yuan, “Persistence and stability of the disease-free equilibrium in a stochastic epidemic model with imperfect vaccine,”
*Advances in Difference Equations*, vol. 2016: Article ID 280, 14 pages, 2016. View at: Publisher Site | Google Scholar | MathSciNet - L. Liu, “A delayed {SIR} model with general nonlinear incidence rate,”
*Advances in Difference Equations*, vol. 2015: Article 329, 11 pages, 2015. View at: Publisher Site | Google Scholar | MathSciNet - Y. N. Xiao and L. S. Chen, “An SIS epidemic model with stage structure and a delay,”
*Acta Mathematicae Applicatae Sinica*, vol. 18, pp. 607–618, 2002. View at: Google Scholar | MathSciNet - Q. Liu and P. Van Mieghem, “Evaluation of an analytic, approximate formula for the time-varying {SIS} prevalence in different networks,”
*Physica A. Statistical Mechanics and its Applications*, vol. 471, pp. 325–336, 2017. View at: Publisher Site | Google Scholar | MathSciNet - J. Wang, “Analysis of an SEIS Epidemic Model with a Changing Delitescence,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 318150, 10 pages, 2012. View at: Publisher Site | Google Scholar - 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 | MathSciNet - F.-F. Zhang, Z. Jin, and G.-Q. Sun, “Bifurcation analysis of a delayed epidemic model,”
*Applied Mathematics and Computation*, vol. 216, no. 3, pp. 753–767, 2010. View at: Publisher Site | Google Scholar | MathSciNet - W. Wang and L. Chen, “Stability and Hopf Bifurcation Analysis of an Epidemic Model by Using the Method of Multiple Scales,”
*Mathematical Problems in Engineering*, vol. 2016: Article 2034136, 8 pages, 2016. View at: Publisher Site | Google Scholar | 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 | MathSciNet - J.-Z. Zhang, Z. Jin, Q.-X. Liu, and Z.-Y. Zhang, “Analysis of a delayed {SIR} model with nonlinear incidence rate,”
*Discrete Dynamics in Nature and Society*, vol. 2008, Article ID 636153, 16 pages, 2008. View at: Publisher Site | Google Scholar | MathSciNet - L. Wen and X. Yang, “Global stability of a delayed {SIRS} model with temporary immunity,”
*Chaos, Solitons & Fractals*, vol. 38, no. 1, pp. 221–226, 2008. View at: Publisher Site | Google Scholar | MathSciNet - Y. N. Kyrychko and K. B. Blyuss, “Global properties of a delayed {SIR} model with temporary immunity and nonlinear incidence rate,”
*Nonlinear Analysis. Real World Applications*, vol. 6, no. 3, pp. 495–507, 2005. View at: Publisher Site | Google Scholar | MathSciNet - P. Yongzhen, S. Li, C. Li, and S. Chen, “The effect of constant and pulse vaccination on an {SIR} epidemic model with infectious period,”
*Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems*, vol. 35, no. 8, pp. 3866–3878, 2011. View at: Publisher Site | Google Scholar | MathSciNet - Y. Song, M. Han, and J. Wei, “Stability and Hopf bifurcation analysis on a simplified {BAM} neural network with delays,”
*Physica D. Nonlinear Phenomena*, vol. 200, no. 3-4, pp. 185–204, 2005. View at: Publisher Site | Google Scholar | MathSciNet - B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan,
*Theory and Applications of Hopf Bifurcation*, Cambridge University Press, Cambridge, UK, 1981. - M. Ferrara, L. Guerrini, and G. M. Bisci, “Center Manifold Reduction and Perturbation Method in a Delayed Model with a Mound-Shaped Cobb-Douglas Production Function,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 738460, 6 pages, 2013. View at: Publisher Site | Google Scholar - X.-Y. Meng, H.-F. Huo, X.-B. Zhang, and H. Xiang, “Stability and Hopf bifurcation in a three-species system with feedback delays,”
*Nonlinear Dynamics*, vol. 64, no. 4, pp. 349–364, 2011. View at: Publisher Site | Google Scholar | MathSciNet - C. Xu, X. Tang, and M. Liao, “Stability and bifurcation analysis of a six-neuron BAM neural network model with discrete delays,”
*Neurocomputing*, vol. 74, no. 5, pp. 689–707, 2011. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2017 Juan Liu. 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.