The Scientific World Journal

Volume 2013 (2013), Article ID 209256, 5 pages

http://dx.doi.org/10.1155/2013/209256

Research Article

Dynamic Behavior for an SIRS Model with Nonlinear Incidence Rate and Treatment

Junhong Li and Ning Cui

Department of Mathematics and Sciences, Hebei Institute of Architecture & Civil Engineering, Zhangjiakou, Hebei 075000, China

Received 23 September 2013; Accepted 6 November 2013

Academic Editors: W. Jarczyk and L. Kocinac

Copyright © 2013 Junhong Li and Ning Cui. 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.

Abstract.

This paper considers an SIRS model with nonlinear incidence rate and treatment. It is assumed that susceptible and infectious individuals have constant immigration rates. We investigate the existence of equilibrium and prove the global asymptotical stable results of the endemic equilibrium. We then obtained that the model undergoes a Hopf bifurcation and existences a limit cycle. Some numerical simulations are given to illustrate the analytical results.

1. Introduction

Treatment is an important and effective method to prevent and control the spread of various infectious diseases, such as measles, tuberculosis, and flu [1–4]. In classical epidemic models, the treatment rate of the infective is assumed to be proportional to the number of the infective individuals [5]. This is unsatisfactory because the resources for treatment should be quite large. In fact, every community should have a suitable capacity for treatment. If it is too large, the community pays for unnecessary cost. If it is small, the community has the risk of the outbreak of a disease. Thus, it is realistic to maintain a suitable capacity of disease treatment. Wang and Ruan [6] considered an SIR model in which the capacity for the treatment of a disease in a community is a constant. Namely, they used the following function:

which was used by [7]. This seems more reasonable when we consider the limitation of the treatment resource of a community.

There are many reasons for using nonlinear incidence rate, and various forms of nonlinear incidence rates have been proposed recently. Liu et al. [8] proposed a nonlinear saturated incidence function to model the effect of behavioral changes to certain communicable disease, where describes the infection force of the disease and measures the inhibition effect from the behavioral change of the susceptible individuals when the number of infectious individuals increases. The case when , was used by [9]. We assume the population can be partitioned into three compartments: susceptible, infective, and recovered. Let , , and denote the numbers of susceptible, infective, and recovered individuals, respectively. Motivated by the works [6, 7, 9], we will formulate an SIRS model with nonlinear incidence rate and constant immigration rates for susceptible and infectious individuals [10]. Namely, we consider the following SIRS model:

where is the rate of natural death, is the rate for recovery, is the proportionality constant, is the rate at which recovered individuals lose immunity and return to susceptible class, is the parameter measures of the psychological or inhibitory effect, and , are constant recruitments of susceptible and infective individuals, respectively. It is assumed that all the parameters are positive constants. It is easy to see that the total population size implies . Since tends to a constant as tends to infinity, we assume that the population is in equilibrium and investigate the behavior of (2) on the plane . Let be a solution of (2) with initial conditions , , and . This solution will satisfy , , and for all since if  , if , and if . Consequently, is positively invariant for (2). Thus, we restrict our attention to the following reduced model:

where is the treatment constant.

From the epidemiological interpretation, our discussion on (3) will be restricted in the following bounded domain:

which is a positively invariant set for (3).

The paper is organized as follows. In the next section, we investigate the existence and stability of equilibrium for (3). In Section 3, we study the Hopf bifurcation and limit cycle. Some numerical simulations are given to illustrate the analytical results. Section 4 is a brief discussion.

2. Existence and Stability of Equilibrium

In this section, we first consider the existence of equilibrium of (3) and their global stability. In order to find endemic equilibrium of (3), we substitute

into

to obtain the cubic equation

Let , , and be three roots of (7). Then, we get

When , we can see that

then there is a unique positive root of (7). Direct calculations show that

From biological considerations, it is easy to see the positive root

Based on the above analysis, we obtain the following theorem.

Theorem 1. There is a unique endemic equilibrium of (3) if .

Theorem 2. The endemic equilibrium of (3) is locally asymptotically stable if .

Proof. The Jacobian matrix of (3) at takes the form

It is easy to obtain and when . This completes the proof.

Theorem 3. The endemic equilibrium of (3) is globally asymptotically stable if .

Proof. Taking Dulac function

we obtain

where is the vector field of (3). Obviously,

Then by Dulac’s criteria, the system (3) admits no limit cycles or separatrix cycles. The global stability of follows from the Poincare-Bendixson Theorem. This completes the proof.

Theorem 4. There will be one or two endemic equilibria of the system (3) if .

Proof. Based on the analysis of Theorem 1, when , we obtain the roots of (7) satisfying

then there exist one positive real root and two conjugate complex roots with negative real parts or two positive real roots and one negative root. Thus, there will be one or two endemic equilibria in (3). This completes the proof.

3. Hopf Bifurcation

In this section, we study the Hopf bifurcation and limit cycle of the system (3). For simplicity of computation, we consider the following system which is equivalent to (3):

Let , to translate to . Then, (17) becomes

where and represent the higher order terms and

To obtain the Hopf bifurcation, we fix parameters such that , which is equivalent to . Let

then (18) is reduced to

where

Let

we obtain the normal form of the Hopf bifurcation:

where

Set the Lyapunov number by

which can be reduced to

So we have the following Hopf bifurcation results.

Theorem 5. There exist Hopf bifurcation and limit cycle in the system (17), when

To illustrate the theorem, let us consider the following parameters.

(see [11]), (see [12]), , , , and .

We have , and the equilibria , exist (see Theorem 4). The equilibrium is unstable saddle. The parameter values satisfy conditions (29) of Theorem 5 and . Therefore, (17) has an unstable periodic orbit which encircles . Its phase portrait is illustrated in Figure 1. The time series of the infective and recovered individuals are given in Figures 2 and 3, respectively.

Figure 1: The phase portraits of (17).

Figure 2: Time series of infective individuals.

Figure 3: Time series of recovered individuals.

4. Conclusion

In this paper, we discuss an SIRS epidemic model with nonlinear incidence rate and treatment. It is assumed that susceptible and infectious individuals have constant immigration rates. We investigate the existence and stability of equilibria of (3) and study the Hopf bifurcation and limit cycle. Some numerical simulations are given to illustrate the analytical results. Without the treatment and recruitment of infectious, (2) becomes the SIRS model (see [9]).

Acknowledgments

The authors are very grateful to the reviewers for their valuable comments and suggestions. This work was supported by the Youth Science Foundations of Education Department of Hebei Province (nos. 2010233 and 2011236).

References

Z. Feng and H. R. Thieme, “Recurrent outbreaks of childhood diseases revisited: the impact of isolation,” Mathematical Biosciences, vol. 128, no. 1-2, pp. 93–130, 1995.

J. M. Hyman and J. Li, “Modeling the effectiveness of isolation strategies in preventing STD epidemics,” SIAM Journal on Applied Mathematics, vol. 58, no. 3, pp. 912–925, 1998.

L.-I. Wu and Z. Feng, “Homoclinic bifurcation in an SIQR model for childhood diseases,” Journal of Differential Equations, vol. 168, no. 1, pp. 150–167, 2000.

Z. Qiu and Z. Feng, “Transmission dynamics of an influenza model with vaccination and antiviral treatment,” Bulletin of Mathematical Biology, vol. 72, no. 1, pp. 1–33, 2010.

R. M. Anderson and R. M. May, Infectious Diseases of Humans, Oxford University Press, London, UK, 1991.

W. Wang and S. Ruan, “Bifurcations in an epidemic model with constant removal rate of the infectives,” Journal of Mathematical Analysis and Applications, vol. 291, no. 2, pp. 775–793, 2004.

Z. Bai and Y. Zhou, “Existence of two periodic solutions for a non-autonomous SIR epidemic model,” Applied Mathematical Modelling, vol. 35, no. 1, pp. 382–391, 2011.

W.-M. Liu, S. A. Levin, and Y. Iwasa, “Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models,” Journal of Mathematical Biology, vol. 23, no. 2, pp. 187–204, 1985.

D. Xiao and S. Ruan, “Global analysis of an epidemic model with nonmonotone incidence rate,” Mathematical Biosciences, vol. 208, no. 2, pp. 419–429, 2007.

G. Li, W. Wang, and Z. Jin, “Global stability of an SEIR epidemic model with constant immigration,” Chaos, Solitons and Fractals, vol. 30, no. 4, pp. 1012–1019, 2006.

L. Zhou and M. Fan, “Dynamics of an SIR epidemic model with limited medical resources revisited,” Nonlinear Analysis: Real World Applications, vol. 13, no. 1, pp. 312–324, 2012.

X.-Z. Li, W.-S. Li, and M. Ghosh, “Stability and bifurcation of an SIR epidemic model with nonlinear incidence and treatment,” Applied Mathematics and Computation, vol. 210, no. 1, pp. 141–150, 2009.