#### Abstract

We investigate the complex dynamics of a SIRS epidemic model incorporating media coverage with random perturbation. We first deal with the boundedness and the stability of the disease—free and endemic equilibria of the deterministic model. And for the corresponding stochastic epidemic model, we prove that the endemic equilibrium of the stochastic model is asymptotically stable in the large. Furthermore, we perform some numerical examples to validate the analytical finding, and find that if the conditions of stochastic stability are not satisfied, the solution for the stochastic model will oscillate strongly around the endemic equilibrium.

#### 1. Introduction

Epidemiology is the study of the spread of diseases with the objective of tracing factors that are responsible for or contribute to their occurrence. Mathematical modeling has become an important tool in analyzing the epidemiological characteristics of infectious diseases and can provide useful control measures. Various models have been used to study different aspects of diseases spreading [1–11].

Let be the number of susceptible individuals, the number of infective individuals, and the number of removed individuals at time , respectively. A general SIRS epidemic model can be formulated as where is the recruitment rate of the population, is the natural death rate of the population, is the natural recovery rate of the infective individuals, is the rate at which recovered individuals lose immunity and return to the susceptible class, and is the disease-induced death rate. The transmission of the infection is governed by the incidence rate , and is called the infection force.

In modelling of communicable diseases, the incidence rate has been considered to play a key role in ensuring that the models indeed give reasonable qualitative description of the transmission dynamics of the diseases. Some factors, such as media coverage, density of population, and life style, may affect the incidence rate directly or indirectly [12–18]. It is worthy to note that, during the spreading of severe acute respiratory syndrome (SARS) from 2002 to 2004 and the outbreak of influenza A (H1N1) in 2009, media coverage plays an important role in helping both the government authority make interventions to contain the disease and people response to the disease [12, 18]. And a number of mathematical models have been formulated to describe the impact of media coverage on the transmission dynamics of infectious diseases. Especially, Liu and Cui [15], Tchuenche et al. [17], and Sun et al. [16] incorporated a nonlinear function of the number of infective (2) in their transmission term to investigate the effects of media coverage on the transmission dynamics: where is the contact rate before media alert; the terms measure the effect of reduction of the contact rate when infectious individuals are reported in the media. Because the coverage report cannot prevent disease from spreading completely we have . The half-saturation constant reflects the impact of media coverage on the contact transmission. The function is a continuous bounded function which takes into account disease saturation or psychological effects. Then model (1) becomes where all the parameters are nonnegative and have the same definitions as before.

For model (3), the basic reproduction number is the threshold of the system for an epidemic to occur. Model (3) has a the disease-free equilibrium which exists for all parameter values. And the endemic equilibrium of model (3) satisfies which yields where When , we know that , ; hence, model (3) has a unique endemic equilibrium . These results of model (3) were studied in [15].

On the other hand, if the environment is randomly varying, the population is subject to a continuous spectrum of disturbances [19, 20]. That is to say, population systems are often subject to environmental noise; that is, due to environmental fluctuations, parameters involved in epidemic models are not absolute constants, and they may fluctuate around some average values. Therefore, many stochastic models for the populations have been developed and studied [21–40]. But, to our knowledge, the research on the dynamics of SIRS epidemic model incorporating media coverage with random perturbation seems rare.

In this paper, our basic approach is analogous to that of Beretta et al. [24]. They assumed that stochastic perturbations were of white noise type, which were directly proportional to distances , , and from values of , , and , influenced the , , and , respectively. By this method, we formulate our stochastic differential equation corresponding to model (3) as follows: where , , and are real constants and known as the intensity of environmental fluctuations; , , and are independent standard Brownian motions.

The aim of this paper is to consider the stochastic dynamics of model (8). The paper is organized as follows. In Section 2, we carry out the analysis of the dynamical properties of stochastic model (8). And in Section 3, we give some numerical examples and make a comparative analysis of the stability of the model with deterministic and stochastic environments and have some discussions.

#### 2. Mathematical Properties of the Deterministic Model (3)

The following result shows that the solutions for model (3) are bounded and, hence, lie in a compact set and are continuable for all positive time.

Lemma 1. *The plane is an invariant manifold of model (3), which is attracting in the first octant.*

*Proof. *Summing up the three equations in (3) and denoting , we have
Hence, by integration, we check
Hence,
which implies the conclusion.

Therefore, from biological consideration, we study model (3) in the closed set

Proposition 2 is proved in [15] and is here just recalled.

Proposition 2. *(i) The disease-free equilibrium is globally asymptotically stable if and unstable if in the set .**(ii) The endemic equilibrium of model (3) is locally asymptotically stable if in the set .*

Next, we present the following theorem which gives condition for the global asymptotical stability of the endemic equilibrium of model (3).

Theorem 3. *If , the endemic equilibrium of model (3) is globally asymptotically stable in the set .*

*Proof. *By summing all the equations of model (3), we find that the total population size verify the following equation:
where .

It is convenient to choose the variable instead of . That is, consider the following model:
changing the variables such that , , and , where , so model (14) becomes as follows:
Consider the function
where and are positive constants which will be chosen later. Then the derivative of along the solution for model (15) is given by
Let us choose and such that
then Thus, we have
By applying the Lyapunov-LaSalle asymptotic stability theorem [41, 42], the endemic equilibrium of model (3) is globally asymptotically stable. This completes the proof.

*Example 4. *We now use the parameter values
and show the stability of the endemic equilibrium of model (3). Model (3) becomes
Note that
From Theorem 3, one can therefore conclude that, for any initial values , the endemic equilibrium of model (21) is globally stable (see Figure 1).

#### 3. Stochastic Stability of the Endemic Equilibrium of Model (8)

Throughout this paper, let be a complete probability space with a filtration satisfying the usual conditions (i.e., it is right continuous and increasing while contains all -null sets).

Considering the general -dimensional stochastic differential equation on with initial value , the solution is denoted by . Assume that and for all , so (23) has the solution . This solution is called the trivial solution.

*Definition 5 ( see [43]). *The trivial solution of (23) is said to be as follows:(i)stable in probability if for all ,
(ii)asymptotically stable if it is stable in probability and, moreover,
(iii)asymptotically stable in the large if it is stable in probability and, moreover, for all Define the differential operator associated to (23) by
If acts on a function , then
where means transposition. For more definitions of stability we refer to [43].

In the following, we will give the result of the asymptotical stability in the large of the endemic equilibrium of model (8).

If , stochastic model (8) can center at its endemic equilibrium . By the change of variables we obtain the following system: where

It is easy to see that the stability of the endemic equilibrium of model (8) is equivalent to the stability of the trivial solution for model (30).

Before proving the stochastic stability of the trivial solution for model (30), we put forward a Lemma in [44].

Lemma 6 (see [44]). *Suppose that there exists a function satisfying the following inequalities:
**
where and is positive constant. Then the trivial solution for model (30) is exponentially -stable for all time . When , it is usually said to be exponentially stable in mean square and the trivial solution is asymptotically stable in the large.*

From the above Lemma, we obtain the following theorem.

Theorem 7. *Assume that . If the following conditions are satisfied:
**
where
**
then the trivial solution of model (30) is asymptotically stable in the large. And the endemic point of model (8) is asymptotically stable in the large.*

*Proof. *We define the Lyapunov function as follows:
where , and are real positive constants to be chosen later. It is easy to check that inequalities (32) are true.

Furthermore, by the Itô formula, we have
Then we have
Choose and
then
Moreover, using Cauchy inequality to , , and , we can obtain
Substituting (39) and (40) into (37), yields
where
Let us choose such that

On the other hand, the conditions in (33) are satisfied, so , , and are positive constants. Let ; then . From (41), one sees that
According to Lemma 6, we therefore conclude that the trivial solution of model (30) is asymptotically stable in the large. We therefore have the assertion.

Next, for further studying the effects of noise on the dynamics of model (8), we give some numerical examples to illustrate the dynamical behavior of stochastic model (8) by using the Milstein method mentioned in Higham [45]. In this way, model (8) can be rewritten as the following discretization equations: where , , and , , are the Gaussian random variables .

The parameters of model (8) are fixed as (20). Then model (8) has the endemic point . And model (8) becomes Choosing and noting that

It is easy to see that all the conditions of Theorem 7 are satisfied, and we can therefore conclude that the endemic point of model (47) is asymptotically stable in the large. The numerical examples shown in Figure 2(b) clearly support these results. To further illustrate the effect of the noise intensity on model (47), we keep all the parameters in (20) unchanged but increase to . In this case, we can therefore conclude, by Theorem 7, that for any initial value , the endemic point of model (47) is asymptotically stable in the large (see Figure 2(b)).

(a) |

(b) |

In the above case, if we adopt and keep the other parameters unchanged, in this case, model (47) has the endemic point . And it is easy to compute Therefore, the conditions of Theorem 7 are not satisfied, and the solution of model (47) will oscillate strongly around the endemic point , which is not asymptotically stable in the large (see Figure 3).

#### 4. Conclusions and Discussions

In this paper, by using the theory of stochastic differential equation, we investigate the dynamics of a SIRS epidemic model incorporating media coverage with random perturbation. The value of this study lies in two aspects. First, it presents some relevant properties of the deterministic model (3), including boundedness and the stability of the disease-free and endemic points. Second, it verifies the stochastic stability in the large of the endemic equilibrium for the stochastic model (8).

From the theoretical and numerical results, we can know that, when the noise density is not large, the stochastic model (8) preserves the property of the stability of the deterministic model (3). To a great extent, we can ignore the noise and use the deterministic model (3) to describe the population dynamics. However, when the noise is sufficiently large, it can force the population to become largely fluctuating. In this case, we can not use deterministic model (3) but instead stochastic model (8) to describe the population dynamics. Needless to say, both deterministic and stochastic epidemic models have their important roles.

#### Acknowledgments

The author would like to thank the editors and referees for their helpful comments and suggestions. This research was supported by the National Science Foundation of China (61272018 & 61373005) and Zhejiang Provincial Natural Science Foundation (R1110261).