Table of Contents Author Guidelines Submit a Manuscript
International Journal of Stochastic Analysis
Volume 2013 (2013), Article ID 431257, 4 pages
Research Article

Stability Analysis of a Stochastic SIR Epidemic Model with Specific Nonlinear Incidence Rate

1Department of Mathematics and Computer Science, Faculty of Sciences Ben M'sik, Hassan II University, P.O. Box 7955, Sidi Othman, Casablanca, Morocco
2Centre Régional des Métiers de l'Education et de la Formation (CRMEF), Derb Ghalef, Casablanca, Morocco

Received 30 May 2013; Accepted 21 August 2013

Academic Editor: Josefa Linares-Pérez

Copyright © 2013 Jihad Adnani et al. 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.


We investigate a stochastic SIR epidemic model with specific nonlinear incidence rate. The stochastic model is derived from the deterministic epidemic model by introducing random perturbations around the endemic equilibrium state. The effect of random perturbations on the stability behavior of endemic equilibrium is discussed. Finally, numerical simulations are presented to illustrate our theoretical results.

1. Introduction

Many mathematical models have been developed in order to understand disease transmissions and behavior of epidemics. One of the earliest of these models was used by Kermack and Mckendrick [1], by considering the total population into three classes, namely, susceptible individuals, infected individuals, and recovered individuals which is known to us as SIR epidemic model. This SIR epidemic model is very important in today's analysis of diseases.

The disease transmission process is unknown in detail. However, several authors proposed different forms of incidences rate in order to model this disease transmission process. In this paper, we consider the following model with specific nonlinear incidence rate: where is the recruitment rate of the population, is the natural death rate of the population, is the death rate due to disease, is the recovery rate of the infective individuals, is the infection coefficient, and is the incidence rate, where are constants. It is very important to note that this incidence rate becomes the bilinear incidence rate if , the saturated incidence rate if or , the modified saturated incidence rate proposed in [2, 3] when , and Crowley-Martin functional response presented in [46] if .

On the other hand, environmental fluctuations have great influence on all aspects of real life. The aim of this work is to study the effect of these environmental fluctuations on the model (1). We assume that the stochastic perturbations are of white noise type and that they are proportional to the distances of and , respectively. Then, the system (1) will be extended to the following system of stochastic differential equation: where , are the positive points of equilibrium for the corresponding deterministic system (1), are independent standard Brownian motions, and represent the intensities of , respectively.

The rest of paper is organized as follows. In the next section, we present the stability analysis of our stochastic model (2). In Section 3, we present the numerical simulation to illustrate our result. The conclusion of our paper is in Section 4.

2. Stability Analysis of Stochastic Model

Clearly, the system (1) has a basic reproduction number given by Using the results presented by Hattaf et al. in [7], it is easy to show that if , the system (1) has just one disease-free equilibrium which is globally asymptotically stable; otherwise, if , the disease-free equilibrium is still present and is unstable, but there is also a unique positive endemic equilibrium , where , with , and . This endemic equilibrium is globally asymptotically stable.

The system (2) has the same equilibria as the system (1). We assume that , and we discuss the stability of the endemic equilibrium of (2). The stochastic system (2) can be centered at its interior endemic equilibrium by the changes of the variables as follows: Hence, the linearized version corresponding to the stochastic model (2) around is given by the following form: where and the superscript “” represents transposition.

Clearly the endemic equilibrium corresponds to the trivial solution in (5).

Let be the family of nonnegative functions defined on such that they are continuously differentiable with respect to and twice with respect to . From [8], we define the differential operator for a function by where and .

According to [8], we have the following theorem.

Theorem 1. Suppose that there exists a function satisfying the following inequalities: where ,  , and are positive constants. Then, the trivial solution of (5) is exponentially -stable for . Moreover, if , then the trivial solution is also called asymptotically mean square stable and it is globally asymptotically stable in probability.

From Theorem 1, we get the conditions for stochastic asymptotic stability of trivial solution of (5) which are given by the following theorem.

Theorem 2. Assume that , , and hold. Then, the trivial solution of (5) is asymptotically mean square stable.

Proof. We consider the following Lyapunov function: where and are nonnegative constants that will be chosen later. It is easy to verify that inequality (8) holds true with .
By applying the operator on , we get If we choose , then where From the assumptions of the theorem, we deduce that and . Hence, is a symmetric positive definite matrix. Let denote the minimum of its two positive eigenvalues and ; then, we can easily get According to Theorem 1, we conclude that the trivial solution of system (5) is globally asymptotically stable.

3. Numerical Simulations

In this section, we present the numerical simulations to illustrate our theoretical results.

We use the following parameter values: , , , , , , , , , and . In this case, we have , , and . By applying Theorem 2, we deduce that the endemic equilibrium is globally asymptotically stable. Figure 1 demonstrates the above analysis.

Figure 1: Deterministic and stochastic trajectories of models (1) and (2) with parameter values , , , , , , , , , and .

4. Conclusion

The purpose of this work is to study the effects of the environmental fluctuations on dynamical behavior of a deterministic SIR epidemic model with specific nonlinear incidence rate by considering the white noise perturbation around the endemic equilibrium state. We have shown that our stochastic model is globally asymptotically stable in probability when the intensities of white noise are less than certain threshold of parameters. However, if these intensities of white noise are zero, which meant that there is no environmental stochastic perturbation, then the conditions of Theorem 2 are reduced to the condition , which gives a nonlinear stability condition for the deterministic model (1).

From our analytical and numerical results, we conclude that the main factor that affects the stability of the stochastic model is the intensities of white noise. In addition, our main results extend the corresponding results in paper [3] and those in [9] when the value of the parameter is equal to one into the stochastic model [9].


  1. W. O. Kermack and A. G. McKendrick, “Contributions to the mathematical theory of epidemics,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 115, pp. 700–721, 1927. View at Publisher · View at Google Scholar · View at Scopus
  2. A. Kaddar, “On the dynamics of a delayed sir epidemic model with a modified saturated incidence rate,” Electronic Journal of Differential Equations, vol. 2009, no. 133, pp. 1–7, 2009. View at Google Scholar · View at Scopus
  3. P. Das, D. Mukherjee, and Y. H. Hsieh, “An S-I epidemic model with saturation incidence: discrete and stochastic version,” International Journal of Nonlinear Analysis and Applications, pp. 1–9, 2011. View at Google Scholar
  4. P. H. Crowley and E. K. Martin, “Functional responses and interference within and between year classes of a dragony population,” Journal of the North American Benthological Society, vol. 8, pp. 211–221, 1989. View at Google Scholar
  5. X. Zhou and J. Cui, “Global stability of the viral dynamics with crowley-martin functional response,” Bulletin of the Korean Mathematical Society, vol. 48, no. 3, pp. 555–574, 2011. View at Publisher · View at Google Scholar · View at Scopus
  6. X. Q. Liu, S. M. Zhong, B. D. Tian, and F. X. Zheng, “Asymptotic properties of a stochastic predator-prey model with Crowley-Martin functional response,” Journal of Applied Mathematics and Computing, 2013. View at Publisher · View at Google Scholar
  7. K. Hattaf, A. A. Lashari, Y. Louartassi, and N. Yousfi, “A delayed SIR epidemic model with general incidence rate,” Electronic Journal of Qualitative Theory of Differential Equations, vol. 3, pp. 1–9, 2013. View at Google Scholar
  8. V. N. Afanas'ev, V. B. Kolmanowskii, and V. R. Nosov, Mathematical Theory of Control Systems Design, Kluwer Academic, Dordrecht, Netherlands, 1996.
  9. P. Das, D. Mukherjee, and A. K. Sarkar, “Study of an S-I epidemic model with nonlinear incidence rate: discrete and stochastic version,” Applied Mathematics and Computation, vol. 218, no. 6, pp. 2509–2515, 2011. View at Publisher · View at Google Scholar · View at Scopus