Discrete Dynamics in Nature and Society

Volume 2016 (2016), Article ID 3458965, 6 pages

http://dx.doi.org/10.1155/2016/3458965

## Modelling the Impact of Media in Controlling the Diseases with a Piecewise Transmission Rate

Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China

Received 23 October 2015; Accepted 5 January 2016

Academic Editor: Amit Chakraborty

Copyright © 2016 Maoxing Liu 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.

#### Abstract

An epidemic model with media is proposed to describe the spread of infectious diseases in a given region. A piecewise continuous transmission rate is introduced to describe that the media has its effect when the number of the infected exceeds a certain critical level. Furthermore, it is assumed that the impact of the media on the contact transmission is described by an exponential function. Stability analysis of the model shows that the disease-free equilibrium is globally asymptotically stable if the basic reproduction number is less than unity. On the other hand, when the basic reproduction number is greater than unity, a unique endemic equilibrium exists, which is also globally asymptotically stable. Our analysis implies that media coverage plays an important role in controlling the spread of the disease.

#### 1. Introduction

When an infectious disease outbreaks, the spread and the control of the disease will be reported by the media including television programs, newspapers, and online social networks. Many examples are the massive reports and daily updates in the public media on the number of the infections and deaths, which had important impacts on the diseases control [1, 2]. It is shown that media coverage plays an important role in the spread and control of the infectious disease, such as the SARS in 2003 [3], the H1N1 in 2009 influenza epidemic [4, 5], and the Ebola in 2014 in Africa [6].

Recently, such impact on disease spreading and controlling has been investigated by mathematical modeling approach [5, 7–20]. In these research works some focused their attention on the incidence rate, and the recent survey identified three typical terms in [5]. Liu et al. [7] described the impact of media coverage using the transmission coefficient , and this impact leads to the change of avoidance and contact patterns at both individual and community levels. Cui et al. [8–11] developed a compartment model using incidence rate with to investigate the impact of media coverage on the transmission, and stability analysis of the models has shown that Hopf bifurcation can occur. Tchuenche et al. [14] used an exponentially decreasing function to reveal the force of media and showed the potential short-term beneficial effect of awareness programs.

More recently nonsmooth media functions [21–25] have been studied. Xiao et al. [21] introduced a segmented function to describe the media impact , here is also the strength of the media effect, and is a threshold that people take the controlling measures or not. A Filippov epidemic model was proposed to describe the real characteristics of media impact in the spread of infectious diseases by incorporating a piecewise continuous transmission rate in [24], and mathematical analysis with regard to the local and global stability of equilibria and local sliding bifurcations are performed.

In fact, during infectious disease outbreaks, individuals may reduce their activities after receiving information about the risk of infection. For example, people will reduce the time that they go out, students will not attend school, and so on, and such information on the ongoing epidemics may impact the dynamics itself. In fact, at the initial stages of the prevalence of disease, most people and public mass media are unaware of the disease; thus, the individuals will not do any protective measures. Only when the number of infectious individuals reaches and exceeds a certain level, the individuals will take precautionary measures against the diseases. Based on the above facts, in this paper we also focus on the incidence rate. Here we introduce to show this fact, where is the threshold value. It is assumed that the impact of media is described by an exponential decreasing factor and the population obeys the logistic growth.

The rest of this paper is organized as follows: in the next section, a mathematical model is proposed in order to reveal the effect of media; then the existence of equilibria is given in Section 3. In Section 4 the local and global stability of the disease-free and the unique endemic equilibria are analyzed. Furthermore, in Section 5 some numerical simulations and discussions are given in the last section.

#### 2. Mathematical Model with Media

In this model the population is divided into three types: the susceptible, the infective and the recovered. Let , , and denote the number of susceptible, infective, and recovered individuals at time , respectively. It is assumed that the growth of the susceptible population obeys the logistic growth, the intrinsic growth rate of the population is , and the carrying capacity for the population is . The interactions between susceptible and infective individuals are assumed to be bilinear, and is the contact rate of susceptible with infective individuals, is the recovery rate of the infective individuals, and is the natural death rate of the population. Consider that the cumulative density of media about the disease is . The growth rate of cumulative density of media, , is assumed to be proportional to the number of infective individuals in the population, and is the depletion rate of cumulative density of media. Thus, we have the following model:

In this model we assume that the impact of media is described by an exponential decreasing factor as , where the parameter is to reflect the impact of media. When , the transmission rate is the constant . When , as one can see if is relatively small, the transmission rate is close to the constant . Thus, the model can be rewritten as follows:withwhere the function means that when infective individuals reach and exceed the certain level , the media has its effects.

For the solutions of system (2) with (3), the region of attraction is given by the set: and it attracts all solutions initiating in the interior of the positive orthant. The phase space is split into two parts: and . In region , there is no effect of the media and the transmission rate , while in region , the transmission rate declines to

For convenience, let the vector , and denote Then system (2) with (3) can be rewritten as the following switching system:In the following, system (6) defined in region is the subsystem and that defined in region is the subsystem , and we consider the two subsystems in isolation firstly.

#### 3. Existence of Equilibria

##### 3.1. Equilibria of

Let be zero; one can verify that the origin is a hyperbolic saddle point with eigenvalues , , . The subsystem has one disease-free equilibrium , and the stability of can be obtained in the following section.

It follows from [26] that the reproduction number and we can verify that when , has a unique endemic equilibrium , where

##### 3.2. Equilibria of

The disease-free equilibria of the subsystem are same with the subsystem ; thus, we study the existence of the endemic equilibrium of .

Let be zero; one can see that the endemic equilibrium must satisfy If there exists a positive equilibrium, it is a positive solution of

One can verify that , , and if , . Hence, the two curves and have at least one positive intersection. In order to determine the number of positive intersections we consider the tangency of two curves. If the two curves intersect, it must have , ; that is, From (10) must satisfy the quadratic equation:Let , ; then we have

Let ; if , (11) has no roots and the system has a unique endemic equilibrium. If , (11) has one unique root and the system has one endemic equilibria of multiplicity at least two. If , (11) has two roots and the system has three endemic equilibria. Thus, we have the following theorem.

Theorem 1. *When and subsystem has a unique endemic equilibrium .*

#### 4. Stability of Equilibria

In this section we present the locally and globally asymptotical stability of the equilibria of subsystem and subsystem .

Theorem 2. *The disease-free equilibrium is locally asymptotically stable if and unstable when .*

*Proof. *The characteristic equation corresponding to the subsystem at is the form where is the eigenvalue. We get Thus, is locally asymptotically stable if and unstable when . The proof of the disease-free equilibrium of the subsystem is similar to ; then this proof is omitted.

Theorem 3. *The disease-free equilibrium is globally asymptotically stable if .*

*Proof. *To establish the global stability of the disease-free equilibrium of subsystem , we consider the positive definite function by using Lyapunov’s method. Now differentiating with respect to , we get If , then . Thus, LaSalle’s invariance principle implies that is globally asymptotically stable in . We can use the same positive definite function when we prove the global stability of the disease-free equilibrium of subsystem .

In fact, when curves starting from will converge to that belongs to , and curves starting from also converge to ; therefore, is globally asymptotically stable. In the following we consider the stability of the endemic equilibrium of subsystem and of subsystem . According to Hurwitz criterion, it is easy to get the following proposition.

Proposition 4. *The unique endemic equilibrium of subsystem is locally asymptotically stable if . The unique endemic equilibrium of subsystem is locally asymptotically stable if .*

When the endemic equilibrium belongs to , we get from the expression of , so if , belongs to , which lead to the existence of that belongs to . Let ; then we have the following result.

Theorem 5. *When the endemic equilibrium is globally asymptotically stable; when the endemic equilibrium is globally asymptotically stable.*

*Proof. *Let a Dulac function and we make , , and for the subsystem . Then Hence, we can exclude the existence of limit cycles in region . The same method can be applied in discussion for region .

When , it is easy to get the endemic equilibrium that belongs to and all curves starting from will verge to . We know that all curves starting from will converge to , but also belongs to if , and there is a unique endemic equilibrium and no limit cycles in region ; hence, curves starting from will enter and converge to . Thus, when the endemic equilibrium is globally asymptotically stable. Similarly when the endemic equilibrium belongs to ; hence, curves starting from will enter and converge to the unique endemic equilibrium in region . So when the endemic equilibrium is globally asymptotically stable.

#### 5. Numerical Simulations and Results

To check the analysis in Section 4 about and , we do some numerical simulations in this section.

When , we set the following values of parameters: , , , , , , and , and we get , , and . Next we consider another case by taking , , and , and get , . We, respectively, choose two sets of initial values that belong to and two sets of initial values that belong to . Thus, we get , ; see Figure 1. As shown in Figure 1(a) the trajectories all approach towards no matter whether or . Hence, is globally asymptotically stable when . Whereas when the trajectories approach towards what Figure 1(b) responses. In conclusion, and are globally asymptotically stable if .