Abstract

We discuss multigroup SIRS (susceptible, infectious, and recovered) epidemic models with random perturbations. We carry out a detailed analysis on the asymptotic behavior of the stochastic model; when reproduction number , we deduce the globally asymptotic stability of the endemic equilibrium by measuring the difference between the solution and the endemic equilibrium of the deterministic model in time average. Numerical methods are employed to illustrate the dynamic behavior of the model and simulate the system of equations developed. The effect of the rate of immunity loss on susceptible and recovered individuals is also analyzed in the deterministic model.

1. Introduction

To curb the spread and impact of viruses, it is important to study their feature, propagating methods, means, and limitation. Many studies about the outbreak and spread of disease have been done by means of establishing epidemic models. These researches provided some useful and valid reference for the characteristic of disease transmission. Based on these results of theoretical analysis, one can predict the future course of an outbreak and evaluate strategies to control an epidemic. In 1927, Kermack and McKendrick created an SIR model in which they considered a fixed population with only three compartments of three classes: susceptible , infected , and removed . represents the number of individuals not yet infected with the disease at time , or those susceptible to the disease. denotes the number of individuals who have been infected with the disease and are capable of spreading the disease to those in the susceptible class. is the compartment used for those individuals who have been infected and then removed from the disease, either due to immunization or due to death. Those in this class are not able to be infected again or to transmit the infection to others. A single group SIRS model is an extension of the SIR model. The SIRS model for infections that do not confer permanent immunity, that is, an infection that does not leave a long lasting immunity: thus individuals that have recovered return to being susceptible again, moving back into the compartment. The only difference between SIR and SIRS is that SIRS model allows members of the recovered class to be free of infection and rejoin the susceptible class.

Considering different contact patterns, a distinct number of sexual partners, or different geography, and so forth, it is more appropriate to divide individual hosts into groups in modeling epidemic disease [1]. Therefore, it is reasonable to propose multigroup models to describe the transmission dynamics of viruses in heterogeneous host populations on epidemic models. In fact, there are already many scholars focusing their study on various forms of multigroup epidemic models (see [27]). They have also proved the global stability of the unique endemic equilibrium through Lyapunov function, which is one of the main mathematical challenges in analyzing multigroup models. Recently, Guo et al. have paid more attention to multigroup models [810]. They first proposed a graph-theoretic approach to the method of global Lyapunov functions and used it to establish the global stability of the interior equilibrium for more general models [8]. Muroya et al. considered a class of -group () SIRS epidemic models described by the following system of equations [11]: Nakata et al. [12] and Enatsu et al. [13] proposed an idea to extend Lyapunov functional techniques in McCluskey [14] for SIR epidemic models to SIRS epidemic models. By extending well-known Lyapunov function techniques, Muroya et al. [11] succeeded to prove the global stability of system (1) without use of the grouping technique by graph theory 4 in Guo et al. [10]. Note that because of environmental noises, the deterministic approach has some limitations in the mathematical modeling transmission of an infectious disease, several authors began to consider the effect of white noise in epidemic models [1518]. Beretta et al. proved the stability of epidemic model with stochastic time delays influenced by probability under certain conditions [19]. Such type of stochastic perturbations firstly was proposed in [19, 20] and later was successfully used in many other papers for many other different systems (see, e.g., [2127]). Yuan et al. in [28] and Yu et al. in [29] all investigated epidemic models with fluctuations around the positive equilibrium and they proved locally stochastically asymptotic stability of the positive equilibrium. Ji et al. also discuss a multigroup SIR model with stochastic perturbation and deduce the globally asymptotic stability of the disease-free equilibrium when , which means the disease will die out. When , they derive the disease will prevail, which is measured through the difference between the solution and the endemic equilibrium of the deterministic model in time average [1]. Imhof and Walcher [30] considered a stochastic chemostat model and they proved that the stochastic model led to extinction even though the deterministic counterpart predicts persistence. In our previous work, we considered an SEIR epidemic model with constant immigration and random fluctuation around the endemic equilibrium, we carried out a detailed analysis on the asymptotic behavior of the stochastic model [31], and we also investigated a two-group epidemic model with distributed delays and random perturbation [32]. In addition, we investigated multigroup SEIQR models with and without random perturbation in computer network [33]. In the present paper, based on system (1), to examine the influence of white noise on system (1) we also consider a stochastic version of the SIRS model by perturbing the deterministic system (1) by a white noise and assume that the perturbations are around the positive endemic equilibrium of epidemic models.

This paper is organized as follows. We begin in Section 2 with necessary background with respect to the deterministic multigroup SIRS model. We establish the global dynamics determined by the basic reproduction number and introduce some results of graph theory used by Guo et al. We also review the important results in Theorem 2 on deterministic model (1) by Muroya et al. In Section 3 we derive the stochastic version from the deterministic model (1). In Section 4 we analyse the asymptotic behavior of the stochastic model by means of the method of Lyapunov functions and the theories of stochastic differential equation in Theorem 5. Numerical methods are employed to simulate the dynamic behavior of the model, and the effect of the rate of immunity loss on the recovered is also analyzed in the deterministic models and the corresponding stochastic models in Section 5. Finally, we give the conclusion of our paper in Section 6.

2. Global Stability of Deterministic Multigroup SIRS Models

First, let us review some theories and results on deterministic multigroup SIRS models; we summarize the parameters in the model (1) by the following list:: the recruitment rate of the population,: transmission coefficient between compartments and ,, , : the natural death rates of susceptible, infected, and recovered individuals in city , respectively,: the rate of immunity loss in the th group,: the natural recovery rate of the infected individuals in city .

We assume and the rest of the parameters are nonnegative for all . In particular, if there is no transmission of the disease between compartments and .

By the biological meanings, we may assume that For each , adding the three equations in (1) gives . Hence, . Therefore, omega limit sets of system (1) are contained in the following bounded region in the non-negative cone of : It can be verified that region is positively invariant and system (1) always has the disease-free equilibrium , where is the equilibrium of the population in the absence of disease . Let denote the interior of . An endemic equilibrium belongs to satisfying the equilibrium equations Let denote the spectral radius of the matrix The parameter is referred to as the basic reproduction number. Its biological significance is that if , the disease dies out while if , the disease becomes endemic [9]. In the following theorem, we show that the multi-group model (1) has at least one endemic equilibrium when , and is globally stable [9].

The matrix denotes the contact matrix. Associated with , one can construct a directed graph whose vertex represents the th group, . A directed edge exists from vertex to vertex if and only if . Throughout the paper, we assume that is irreducible. This is equivalent to being strongly connected. Biologically, this is the same as assuming that any two groups and have a direct or indirect route of transmission. More specifically, individuals in can infect ones in directly or indirectly. Now consider the linear system where and , , , . Let denote the directed graph associated with matrix (and ), and let denote the cofactor of the entry of .

We have the following fundamental lemma[8].

Lemma 1 (Kirchhoff’s Matrix-Tree Theorem). Assume that is irreducible and . Then the following results hold.(1)The solution space of system (9) has dimension 1, with a basis .(2)For , where is the set of all directed spanning subtrees of that are rooted at vertex , is the weight of a directed tree T, and denotes the set of directed arcs in a directed tree T.

Y. Muroya et al. proved the following theorem.

Theorem 2. Assume that is irreducible and inequality (2) holds. If and , then system (1) has at least one endemic equilibrium in , and is globally asymptotically stable.

3. Multigroup Stochastic SIRS Model

Under assumptions that , , and is irreducible in Theorem 2, we know from Section 2 that there exists a positive endemic equilibrium in . Furthermore, we assume stochastic perturbations are of white noise type, which are directly proportional to distances , and from values of , influence on the , and , respectively. So system (1) results in where , and are independent standard Brownian motions and represent the intensities of , respectively. Obviously, stochastic system (12) has the same equilibrium points as system (1). Hence, when the stochastic model (12) satisfies the conditions of Theorem 2, we can ensure the existence of the positive equilibrium point of the stochastic model (12). Note that because of the influence of the white noise, the components of solutions to (12) can also be negative at some time. But because it is impossible for populations to be less than zero physically, once the negative sections appear on components of solutions of (12), we only need to consider the nonnegative section of the components of solutions. In the next section, we will investigate asymptotic stability of the equilibrium of system (12). To this end, we will construct a class of different Lyapunov functions to achieve our proof under certain conditions.

4. Stochastic Stability of the Endemic Equilibrium

Let us now proceed to discuss asymptotic stability of the system (12). In this paper, unless otherwise specified, let be a complete probability space with a filtration satisfying the usual conditions (i.e., it is increasing and right continuous while contains all -null sets). Let be the Brownian motions defined on this probability space. If , then the stochastic system (12) can be centered at its endemic equilibrium ; by the change of variables we obtain the following system:

It is clear that the stability of equilibrium of system (12) is equivalent to the stability of zero solution of system (14). Considering the -dimensional stochastic differential equation If the assumptions of the existence and uniqueness theorem are satisfied, then for any given initial value , (15) has a unique global solution denoted by . For the purpose of stability, we assume in this section and for all . So (15) admits a solution which is called the trivial solution or equilibrium position. We let denote the family of all nonnegative functions on which are continuously twice differentiable in and once in . Define the differential operator associated with (15) by

If acts on a function , then

Definition 3. (1) The trivial solution of (15) is said to be stochastically stable or stable in probability if for every pair of and , there exists a such that whenever . Otherwise, it is said to be stochastically unstable.
(2) The trivial solution is said to be stochastically asymptotically stable if it is stochastically stable and for every , there exists a such that whenever .
(3) The trivial solution is said to be stochastically asymptotically stable in the large if it is stochastically stable and for all

Before proving the main theorem we put forward a lemma in [34].

Lemma 4 (see [34]). If there exists a positive-definite decrescent radially unbounded function such that is negative-definite, then the trivial solution of (15) is stochastically asymptotically stable in the large.

From the above lemma, we can obtain the stochastically asymptotical stability of equilibrium as following.

Theorem 5. Assume that is irreducible and inequality (2) holds and , ; then if the following condition is satisfied: the endemic equilibrium of system (12) is stochastically asymptotically stable in the large, where

Proof. Set and , , , .
Note that is the Laplacian matrix of the matrix (see Lemma 1). Since is irreducible, the matrices and are also irreducible. Let denote the cofactor of the entry of . We know that the system has a positive solution , where for , by Lemma 1.
It is easy to see that we only need to prove the zero solution of (14) is stochastically asymptotically stable in the large. Let , and . We define the Lyapunov function as follows: where , , are real positive constants to be chosen later. Then it can be described as the quadratic form where is a symmetric positive-definite matrix. So it is obvious that is positive-definite decrescent radially unbounded. For the sake of simplicity, (24) may be divided into three functions: , where
Using Itô’s formula and (5), we compute Let , so . It follows from and that
which implies
Hence, inequality (28) becomes Similarly, from Itô’s formula, we obtain Moreover, using Cauchy inequality, we can obtain where
Hence, we can calculate Set then we have where We let The proofs above show that if the condition (21) is satisfied, , , and are positive constants. Let ; then . From (37) and (38), one sees that where , and is an infinitesimal of higher order of for . Hence, is negative-definite for . According to Lemma 4, we therefore conclude that the zero solution of (12) is stochastically asymptotically stable in the large. The proof is complete.

5. Numerical Simulation

Numerical methods are employed to solve the system (1) and (12) and to depict the behavior of the susceptible, infectious, and recovered nodes with respect to time. We numerically simulate the solution of system (1) and (12) with . In this case, we have Let Hence we obtain , , , , , . We can calculate In the absence of noise, we simulate the global stability of the endemic equilibrium of deterministic system (1) in Figure 1, and we always choose initial value , , , , and .

Next, we consider the effect of stochastic fluctuations of environment to the endemic equilibrium of the corresponding deterministic system. Given the discretization of system (12) for , and , where time increment , and , , are -distributed independent random variables which can be generated numerically by pseudorandom number generators.

As mentioned above, there is an endemic equilibrium of system (1) when , where Figure 2 corresponds to , , , , , and ; the numerical simulation shows that the endemic equilibrium of stochastic system (12) is global asymptotically stable under the condition (21), while Figure 3 corresponds to , , , , , and . Moreover, comparison of Figures 2 and 3 suggests that fluctuations enhance as the noise level increases. Note that the condition (21) is just a sufficient condition. When this condition is not satisfied, the stochastic system (12) may be or may not be stable. For instance, the intensities of Brownian motions , , , , , and do not meet the condition (21), but we can see from Figure 3 that the stochastic system (12) is still asymptotically stable. If we choose , , , , , and , then the solution of the stochastic system (12) is not asymptotically stable but explodes to infinity at the finite time (see Figure 4).

The relationships between and are also observed and are depicted in Figures 5 and 6. Note that these two curves of Figure 5 show that the number of susceptible individuals sharply declines to near as the number of the recovery individuals increases slightly and then increases with the number of the increasing recovered individuals gently and goes on increasing with the number of the decreasing recovered individuals, and finally both reach the steady state values. For the stochastic version, Figure 6 corresponds to , , , , , and ; oscillation appears under environmental driving forces which actually affect the deterministic curves shown in Figure 5. But Figure 6 still maintains the same total trend as Figure 5 and reaches the same equilibrium point as deterministic version. Simulation result agrees with the real life situation.

Figure 7 and Tables 1 and 2 show that when the rate of immunity loss gradually reduces, decreases, but increases; the lower the rate of immunity loss is, the lower the steady state value of the susceptible is, the higher the steady state value of the recovered is. Thus, it will be of great importance for health management to control the rate of immunity loss to keep it in a lower level; for example, when antibody concentrations of recovered individuals are at a low level or are zero, they can be required to undergo vaccination to achieve the protective antibody levels.

6. Conclusion

This paper presented a mathematical study describing the dynamical behavior of an SIRS epidemic model with stochastic perturbations. Our purpose was based on analyzing this behavior using a stochastic model. When the reproduction number is greater than one, we obtained sufficient conditions for stochastic stability of the endemic equilibrium by using a suitable Lyapunov function and other techniques of stochastic analysis. The investigation of this stochastic model revealed that the stochastic stability of depends on the magnitude of the intensity of noise as well as the parameters involved within the model system. Finally, numerical simulations are given to validate the theoretical results. The proposed model, a more accurate epidemic model, can help us to understand the dynamic behavior of virus. Moreover, the theoretical results may provide some useful guidance for making effective countermeasures on virus propagation.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author was partially supported by Project of Science and Technology of Heilongjiang Province of China no. 12531187.