# Stability and Threshold of a Stochastic SIRS Epidemic Model with Vertical Transmission and Transfer from Infectious to Susceptible Individuals.

1. IntroductionSince the earlier works of Kermack and McKendrick (see, e.g., [1, 2]), mathematical models have proven their adequacy to describe and analyze the propagation and control of infectious diseases such as Diphtheria, Polio, Tuberculosis, Cholera, and Toxoplasmosis. Various epidemic models of population dynamics have been proposed (see, e.g., [3-6]). Generally, the contact between individuals leads to the dissemination of a disease. The vertical transmission is also possible, which means the direct transmission of a disease from the mother to a fetus during pregnancy. This hypothesis was not taken into consideration in the basic SIRS epidemic model. To be more realistic, our model considers this new hypothesis. In addition, for some bacterial agent infections, recovery cannot produce immunity for a long time. Infected individuals may recover after some treatments and therapies and then go back directly to the susceptible compartment [7]. Considering those two new hypotheses, we obtain the following SIRS epidemic model:

[mathematical expression not reproducible], (1)

where S(t), I(t), and R(t) denote the numbers of susceptible, infected, and recovered individuals at time t, respectively. A is the recruitment rate of susceptible corresponding to immigration. b and [mu] are the birth rate and natural death rate, respectively. It is assumed that [mu] > b[8]. [[gamma].sub.1] is the transfer rate from the infected class to the susceptible class and [[gamma].sub.2] is the transfer rate from the infected class to the recovered class. k is the rate of individuals recovering and returning to S from R. a is the disease-related death rate. q is the vertical transmission coefficient, with 0 < q < 1 and p = 1 - q. [beta] := [??]c represents the transmission rate, which is the product of the contact rate c and the probability of transmission per contact [??]. The parameters involved in system (1) are all positive constants. According to the theory in [7, 8], the basic reproduction number of system (1) is [R.sub.0] = [beta]A/([mu] - b)([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb). If [R.sub.0] [less than or equal to] 1, the deterministic model (1) has only the disease-free equilibrium [P.sup.0](A/([mu] - b), 0, 0), which is globally asymptotically stable. If [R.sub.0] > 1, [P.sup.0] becomes unstable and there exists an endemic equilibrium [P.sup.*]([S.sup.*], [I.sup.*], [Q.sub.*]) such that

[mathematical expression not reproducible]. (2)

The incidence rate is the number of new infected situations per population in a given time phase. In many previous epidemic models, the bilinear incidence rate is frequently used (see, e.g., [9-11]). However, there exist many forms of nonlinear incidence rate and each form presents some advantages as the following examples:

(1) f(S, I) = ([beta] - ([[beta].sub.1]I/(V + I))SI), ([beta] > [[beta].sub.1], V > 0) [12].

(2) f(S, I) = [beta]SI/(1 + k[I.sup.p]), (k > 0, 0 [less than or equal to] p [less than or equal to] 1) [13].

(3) f(S, I) = [beta]SI/(1 + [k.sub.1]I + [k.sub.2]S), ([k.sub.1] > 0, [k.sub.2] [greater than or equal to] 0) [14].

(4) f(S, I) = [beta][S.sup.u]I/(1 + [S.sup.v]), (u [greater than or equal to] V) [15].

There have been many mathematical models [16-21] committed to studying the impacts of nonlinear transmission on the propagation of a disease. The aim of our work is to generalize the results found by those authors. According to the reality, we shall consider our model with the following general functional response:

f (S, I) = [beta]SI/[1 + [[alpha].sub.1]S + [[alpha].sub.2]I + [[alpha].sub.3]SI, (3)

where [[alpha].sub.1], [[alpha].sub.2], [[alpha].sub.3] [greater than or equal to] 0. It is necessary to mention that 1 + [[alpha].sub.1]S + [[alpha].sub.2]I + [[alpha].sub.3]SI is a general form which represents mutual interference between S and I. In particular cases,

(1) when [[alpha].sub.1] = [[alpha].sub.2] = [[alpha].sub.3] = 0, f(S, I) becomes a bilinear mass-action function response (namely, type I Holling functional response) [22];

(2) when [[alpha].sub.2] = [[alpha].sub.3] = 0, f(S, I) becomes a saturated incidence rate (namely, Holling type II functional response) [20];

(3) when [[alpha].sub.3] = 0, f(S, I) becomes a Beddington-DeAngelis functional response (namely, modified type II Holling functional response) [23];

(4) when [[alpha].sub.3] = [[alpha].sub.1][[alpha].sub.2], f(S, I) becomes a Crowley-Martin functional response introduced in [24].

By all what we have introduced, we consider the following system:

[mathematical expression not reproducible], (4)

where [psi](S, I) = 1 + [[alpha].sub.1]S + [[alpha].sub.2]I + [[alpha].sub.3]SI. The basic reproduction number of system (4) can be presented as follows:

[R.sub.0] = [beta]A/[[([mu] - b) + [[alpha].sub.1]A] ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb)]. (5)

The deterministic model constructed above can be improved by taking into account the unpredictable biological condition. In the real world, biological phenomena are often affected by the environmental noise, and the nature of epidemic growth is inherently random due to the unpredictability of person-person contacts in the case of horizontal transmission, or in the other case of vertical transmission (mother-fetus). That is to say, the parameter [beta] involved in system (4) is not absolute constant and may fluctuate around some average values. Both from a biological and mathematical perspectives, the main of our study is to investigate how the stochastic character of human transmission affects the spread of a disease through studying the dynamical behavior of our stochastic model. We assume that the contact rate [beta] is perturbed by Gaussian white noise, which is presented by [beta] + [sigma][??](t), where B(t) is a standard Brownian motion with intensity [sigma] > 0. Then we obtain the following SIRS epidemic model with perturbation stochastic and general functional response:

[mathematical expression not reproducible]. (6)

Many authors have been introduced random effects into population systems by different techniques (see, e.g., [25-34]). Furthermore, in the study of the dynamical behavior of the epidemic models, we are interested in two situations. One is when the disease goes to extinction, the other is when the disease prevails. Thus many authors have studied this interesting topic. For example, Ji and Jiang [35] investigated the threshold of SIR epidemic model with stochastic perturbation. Zhao and Jiang [36] investigated the threshold of a stochastic SIRS epidemic model with saturated incidence [37]. Then, they considered the threshold of a stochastic SIS epidemic model with vaccination [38]. Also, they studied the threshold of a stochastic SIRS epidemic model with varying population size [39]. Zhao and Yuan [22] studied the threshold behavior of a stochastic SIVR epidemic model with standard incidence and imperfect vaccine. In case that the disease goes extinct, they showed that the disease-free equilibrium is almost surely stable by using the nonnegative semimartingale convergence theorem. In this work, we consider a stochastic SIRS epidemic model with general incidence rate (3). This generalization is the main difficulty to be overcome in establishing the threshold of SIRS epidemic model (6), which has never been examined in the previously studies. In addition, we prove the global stability of SIRS epidemic model (4) with stochastic perturbations, which has not been proven in the previous papers.

This paper is organized as follows. In Section 2, we present some preliminaries which will be used in our following analysis. In Section 3, we show that there is a unique global positive solution of system (6). In Section 4, we give the conditions for the moment exponential stability of the equilibrium [P.sup.0] of stochastic system (6). In Sections 5 and 6, we establish sufficient conditions for persistence and extinction of disease, respectively. In Section 7, we present some numerical simulations to illustrate our main results. The paper ends with a brief discussion.

2. Preliminaries

Throughout this paper, we let ([OMEGA], F, P) be a complete probability space with a filtration [{[F.sub.t]}.sub.t[greater than or equal to]0] satisfying the usual conditions (i.e., [{[F.sub.t]}.sub.t[greater than or equal to]0] is increasing and right continuous while [F.sub.0] contains all P-null sets). B(t) is defined on this complete probability space. We also let [R.sup.d.sub.+] = {x [member of] [R.sup.d] : [x.sub.i] > 0, 1 [less than or equal to] i [less than or equal to] d}.

In general, we consider the d-dimensional stochastic differential equation:

dx (t) = f(x (t), t) dt + g(x (t), t) dB (t) [for all]t [greater than or equal to] [t.sub.0], (7)

with initial value x(0) = [x.sub.0] [member of] [R.sup.d]. B(t) denotes n- dimensional standard Brownian motion defined on the complete probability space ([OMEGA], F, [{[F.sub.t]}.sub.t[greater than or equal to]0], P). [C.sup.2,1]([R.sup.d] x [[t.sub.0], [infinity]];[R.sub.+]) is the family of all nonnegative functions V(x, t) defined on [R.sup.d] x [[t.sub.0],[infinity]] such that they are continuously twice differentiable in x and once in t. The differential operator L of (7) is defined by the following [40]:

[mathematical expression not reproducible]. (8)

If L acts on a function V [member of] [C.sup.2,1] ([R.sup.d] x [[t.sub.0],[infinity]]; [R.sub.+]), then

LV (x, t) = [V.sub.t] (x, t) + [V.sub.x] (x, t) f (x, t) + 1/2 trace [[g.sup.T] (x, t) [V.sub.xx] (x, t) g (x, t)]. (9)

By Ito's formula, if x(t) [member of] [R.sup.d], we have

dV (x (t), t) = LV(x (t), t) dt + [V.sub.x] (x (t), t) g (x (t), t) dB (t). (10)

Next, we shall present the definition of hth moment exponentially stability (see [41]).

Definition 1. The equilibrium x = 0 of the system (7) is said to be hth moment exponentially stable, if there is a pair of positive constants [C.sub.1] and [C.sub.2] such that, for all [x.sub.0] [member of] [R.sup.d],

[mathematical expression not reproducible]. (11)

3. Existence and Uniqueness of the Positive Solution

Since S(t), I(t), and R(t) represent the number of the susceptible, the infected, and the recovered individuals at time t, respectively, they should be nonnegative. So, the first step of our study is to prove that system (6) has a unique global positive solution. We define a bounded set [DELTA] as follows:

[DELTA] := {x = ([x.sub.1], [x.sub.2], [x.sub.3]) : [x.sub.1] > 0, [x.sub.2] > 0, [x.sub.3] > 0, [x.sub.1] + [x.sub.2] + [x.sub.3] < A/[[mu] - b] a.s.}. (12)

Theorem 2. For any initial value (S(0), I(0), R(0)) [member of] [DELTA], there exists a unique positive solution (S(t), I(t), R(t)) of system (6) on t [greater than or equal to] 0, and the solution will remain in [R.sup.3.sub.+] with probability one. That is to say, the solution (S(t), I(t), R(t)) [member of] [R.sup.3.sub.+] for all t [greater than or equal to] 0 almost surely.

Proof. Let (S(0), I(0), R(0)) [member of] [DELTA] and N(t) = S(t) + I(t) + R(t). It is easy to check that

dN(t) = [A- ([mu] - b)N(t) - aI] dt. (13)

Then, if (S(s), I(s), R(s)) [member of] [R.sup.3.sub.+] for all 0 [less than or equal to] s [less than or equal to] t a.s. we get

dN(s) < [A- ([mu] - b)N(s)] d sa.s. (14)

Now, by integration we obtain

N(s) < A/[[mu] - b] + (N(0) - A/[[mu] - b]) [[epsilon].sup.-[mu]s] [for all]s [member of] [0, t] a.s. (15)

Then N(s) < A/([mu] - b), and

S (s), I (s), R (s) [member of] (0, A/[[mu] - b]) [for all]s [member of] [0, t] a.s. (16)

Since the coefficients of system (6) are locally Lipschitz continuous, then for any initial value (S(0), I(0), R(0)) [member of] [R.sup.3.sub.+] there is a unique local solution (S(t), I(t), R(t)) on [0, [[tau].sub.[epsilon]]), where [[tau].sub.[epsilon]] is the explosion time. To show that the solution is global, we only need to prove that [[tau].sub.[epsilon]] = [infinity] a.s.

Let [[epsilon].sub.0] > 0 such that S(0), I(0), R(0) > [[epsilon].sub.0]. For each integer [epsilon] [less than or equal to] [[epsilon].sub.0], we define the following stopping times:

[mathematical expression not reproducible]. (17)

Consider the [C.sup.2]-function [V.sub.1] defined for X = (S,I,R) [member of] [R.sup.3.sub.+] by

[V.sub.1] (X) = -ln (([mu] - b)S/A) - ln(([mu] - b)I/A) - ln(([mu] - b)R/A). (18)

Making use Ito's formula to [V.sub.1], we obtain for all t [greater than or equal to] 0 and s [member of] [0, t [conjunction] [[tau].sub.[epsilon]]]

[mathematical expression not reproducible]. (19)

For all s [member of] [0, t [conjunction] [[tau].sup.[epsilon]]], we have

I(s)/[[psi] (S (s), I (s))] [less than or equal to] A/[([mu] - b) + [[alpha].sub.1]A], S(s)/[[psi] (S (s), I (s))] [less than or equal to] A/[([mu] - b) + [[alpha].sub.1]A]. (20)

Therefore

d[V.sub.1] (X(s)) [less than or equal to] M + [[sigma] (I (s) - S(s))]/ [[psi] (S (s), I (s))] dB (s) a.s., (21)

where M = 3[mu] + [[gamma].sub.1] + [[gamma].sub.2] + a + k + [beta]A/(([mu]- b) + [[alpha].sub.1]A) + [[sigma].sup.2][A.sup.2]/[(([mu]- b) + [[alpha].sub.1]A).sup.2]. Integrating both sides of (21) from 0 to t [conjunction] [[tau].sub.[epsilon]], and after taking the expectation on both sides, we obtain that

E[V.sup.1] (X (t [conjunction] [[tau].sub.[epsilon]])) [less than or equal to] [V.sub.1] (X(0)) + ME (t [conjunction] [[tau].sub.[epsilon]]) [less than or equal to] [V.sub.1] (X(0)) + Mt. (22)

Since [V.sub.1] (X(t [conjunction] [[tau].sub.[epsilon]])) > 0, then

[mathematical expression not reproducible]. (23)

For [[tau].sub.[epsilon]], there is some component of X([[tau].sub.[epsilon]]) equal to [epsilon]. Therefore [V.sub.1](X([[tau].sub.[epsilon]])) [greater than or equal to] - ln (([mu] - b)[epsilon]/A).

Thus

[mathematical expression not reproducible]. (24)

Hence, from (24) we conclude

P([[tau].sub.[epsilon]] [less than or equal to] t) [less than or equal to] - [[V.sub.1] (X(0)) + Mt]/[ln (([mu] - b) [epsilon]/A)]. (25)

Extending [epsilon] to 0, we obtain for all t > 0, P([tau] [less than or equal to] t) = 0. Hence P([tau] = [infinity]) = 1. Thus, [tau] = [[tau].sub.[epsilon]] = [infinity] a.s. which completes the proof of the theorem.

From Theorem 2 and (16) we can conclude the following corollary.

Corollary 3. The set [DELTA] is almost surely positively invariant; that is, if (S(0), I(0), R(0)) [member of] [DELTA], then P((S(t), I(t), R(t)) [member of] [DELTA]) = 1 for all t [greater than or equal to] 0.

4. Moment Exponential Stability

In order to obtain the conditions of moment exponential stability, we will use the following theorem (for the proof of this theorem we refer the reader to [41]).

Theorem 4. Suppose there exists a function V(t, x) [member of] [C.sup.1,2](Rx [R.sup.n]) satisfying the following inequalities:

[K.sub.1] [[absolute value of (x)].sup.h] [less than or equal to] V(t, x) [less than or equal to] [K.sub.2] [[absolute value of (x)].sup.h], (26)

LV (t, x) [less than or equal to] - [K.sub.3] [[absolute value of (x)].sup.h], [K.sub.i] > 0, h > 0. (27)

Then the equilibrium of the system (7) is hth moment exponentially stable. When h = 2, it is usually said to be exponentially stable in mean square and the equilibrium x = 0 is globally asymptotically stable.

From Young's inequality, we have the following inequalities.

Lemma 5. Let h [greater than or equal to] 2 and [epsilon], x, y > 0. Then

[x.sup.h-1] y [less than or equal to] [(h - 1) [epsilon]]/h [x.sup.h] + 1/h[[epsilon].sup.h-1] [y.sup.h], [x.sup.h-2][y.sup.2] [less than or equal to] [(h - 2) [epsilon]]/h [x.sup.h] + 1/h[[epsilon].sup.(h-2)/2] [y.sup.h]. (28)

From Theorem 4, we get the sufficient conditions of the moment exponential stability, which are given by the following theorem.

Theorem 6. Let h [greater than or equal to] 2. If the conditions [R.sub.0] < 1 and

[[sigma].sup.2] < 2/[h - 1] [(A/[([mu] - b) + [[alpha].sub.1]A]).sup.-2] x [([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb) - [beta]A/[([mu] - b) + [[alpha].sub.1]A]] (29)

hold, the disease-free equilibrium [P.sup.0] of system (6) is hth moment exponentially stable in [DELTA].

Proof. Let h [greater than or equal to] 2. Considering the following Lyapunov function,

[V.sub.2] = [[lambda].sub.1] [(A/[[mu] - b] - S).sup.h] + 1/h [I.sup.h] + [[lambda].sub.2] [R.sup.h], (30)

where [[lambda].sub.i], i = 1, 2 are real positive constants to be determined later. It is easy to check that inequalities (26) are true. Then we compute

[mathematical expression not reproducible]. (31)

In [DELTA], we have

[mathematical expression not reproducible]. (32)

By using Lemma 5, we get

[mathematical expression not reproducible]. (33)

Then

[mathematical expression not reproducible]. (34)

We chose [epsilon] to be sufficiently small such that the coefficients of [(A/([mu]-b)-S).sup.h] and [R.sup.h] are negative, and as ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb - [beta]A/(([mu]-b) + [[alpha].sub.1]A)) - ((h - 1)/2) ([A.sup.2][[sigma].sup.2]/[[([mu]-b)+[[alpha].sub.1]A].sup.2]) > 0, we can choose [[lambda].sub.1] and [[lambda].sub.2] being positive such that the coefficient of [I.sup.h] is negative. This ends the proof.

Under the Theorems 4 and 6, we have in the case h = 2 the following corollary.

Corollary 7. If the conditions [R.sup.0] < 1 and [[sigma].sub.2] < 2(A/[(([mu]-b)+ [[alpha].sub.1]A)).sup.-2] [([mu] + [[gamma].sub.1] +[[gamma].sub.2] + a - qb) - [beta]A/(([mu] - b) + [[alpha].sub.1]A)] hold, the disease-free [P.sup.0] of system (6) is globally asymptotically stable in [DELTA].

Remark 8. Biologically, the global stability of free disease equilibrium may be interpreted as the inescapable fate of the epidemic extinction, regardless of its initial situation. Furthermore, the hth moment exponential stability means that the number of infected individuals tends to the state of extinction exponentially fast.

5. Persistence

In this section, our main concern is to determine the conditions for the spread and persistence of an infectious disease. Given that the value of the deterministic threshold [R.sub.0] characterizes the dynamical behaviors of system (4) and guarantees persistence or extinction of the disease. Similarly, we define the threshold of our stochastic SIRS epidemic model (6) as follows:

[R.sup.s.sub.0 = [R.sub.0] [1 - [[sigma].sup.2]A/[2[beta] [([mu] - b) + [[alpha].sub.1]A]]]. (35)

To avoid any ambiguity, we define <x(t)> = (1/t) [[integral].sup.t.sub.0] x(s) ds.

Theorem 9. If [R.sup.s.sub.0] > 1, then the solution (S(t), I(t), R(t)) of system (6) with initial value (S(0), I(0), R(0)) [member of] [DELTA] is persistent in the mean. Moreover, we have

[mathematical expression not reproducible], (36)

where

[I.sub.*] = ([R.sup.s.sub.0] - 1) ([mu] + [r.sub.2] + [delta] + [gamma])[[beta].sup.-1] ([mu] - b) ([mu] + k) (([mu] - b) + [[alpha].sub.1]A)]/[(k + b) ([mu] - b) ([mu] - b + a) + [beta]A([mu] - b) + A([mu] + k) ([[alpha].sub.2] ([mu] - b) + [[alpha].sub.3]A)] > 0. (37)

The following result is needed for the proof of Theorem 9 (see [35]).

Lemma 10. Let h [member of] C([0,[infinity]) x [OMEGA], (0, [infinity])) and H [member of] C([0,[infinity]) x [OMEGA], R). If there exist positive constants [[lambda].sub.0] and [lambda] such that

ln h (t) [greater than or equal to] [[lambda].sub.0]t - [lambda] [[integral].sup.t.sub.0] h (s)d s + H(t) a.s. (38)

for all t [greater than or equal to] 0, and [lim.sub.t[right arrow][infinity]](H(t)/t) = 0 a.s., then

[mathematical expression not reproducible]. (39)

Proof of Theorem 9. We consider a function [V.sub.3] defined by

[V.sub.3] = [[omega].sub.1] (S + I + R) + [[omega].sub.2]S + ln I, (40)

where [[omega].sub.1] and [[omega].sub.2] are real positive constants to be determined later.

Using Ito's formula, we have

[mathematical expression not reproducible]. (41)

Since (S(t), I(t), R(t)) [member of] [DELTA], we get

[mathematical expression not reproducible]. (42)

Inject those two inequalities in expression of d[V.sub.3], we get

[mathematical expression not reproducible]. (43)

Then

[mathematical expression not reproducible]. (44)

In order to eliminate (A/[mu]-S) and R from the last inequality, we choose [[omega].sub.1] and [[omega].sub.2] as follows:

[[omega].sub.1] = [[beta] (b + k)]/[([mu] + k) [([mu] - b) + [[alpha].sub.1]A]], [[omega].sub.2] = [[beta] ([mu] - b)]/[([mu] + k) [([mu] - b) + [[alpha].sub.1]A]]. (45)

By integration, we get

[mathematical expression not reproducible]. (46)

Hence,

[mathematical expression not reproducible], (47)

where

[mathematical expression not reproducible]. (48)

Thus, the strong law of large number for martingales implies that [lim.sub.t[right arrow][infinity]] (G(t)/t) = 0 a.s.

By using Lemma 10, we have

[mathematical expression not reproducible]. (49)

Next, the third equation of system (6) gives

[R (t) - R(0)]/t = [[gamma].sub.2] <I (t)> + ([mu] + k) <R (t)> . (50)

Then

[mathematical expression not reproducible]. (51)

Finally, it follow from system (6) that

d (S + I + R) = [([mu] - b)(A/[[mu] - b] - S) - ([mu] - b + a) I - ([mu] - b)R] dt. (52)

Then

[mathematical expression not reproducible]. (53)

6. Extinction

In this section, we investigate the conditions for the extinction of the disease.

Theorem 11. Let (S(t), I(t), R(t)) be the solution of system (6) with initial value (S(0), I(0), R(0)) in [DELTA].

Assume that

(a) [[sigma].sup.2] > [[beta].sup.2]/[2 ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb)],

or (b) [R.sup.s.sub.0] < 1,

[[sigma].sup.2] [less than or equal to] [[beta] [([mu] - b) + [[alpha].sub.1]A]]/A. (54)

Then

[mathematical expression not reproducible], (55)

[mathematical expression not reproducible]. (56)

Namely, I(t) tend to zero exponentially a.s.; that is, the disease dies out with probability 1.

Proof. It follows from Ito's formula that

d lnI = [[beta]S/[[psi] (S, I)] - ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb) - [[sigma].sup.2]/2 [S.sup.2]/[[[psi].sup.2] (S, I)]] dt + [sigma]S/[[psi] (S, I)] dB. (57)

Integrating this from 0 to t and dividing by t on both sides, we have

[ln I (t)]/t = 1/t [[integral].sup.t.sub.0] [[beta]S(s)/[[psi] (S (s), I (s))] - ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb) - [[sigma].sup.2]/2 S2(s)/[[[psi].sup.2] (S (s), I (s))]] ds + [ln I (0)]/t + M(t)/t, (58)

where M(t) = [[integral].sup.t.sub.0] ([beta]S(s)/[psi](S(s), I(s))) dB(s). By the large number theorem for martingale, we have

[mathematical expression not reproducible]. (59)

If the condition (a) is satisfied, (58) becomes

[mathematical expression not reproducible]. (60)

Taking the limit superior of both sides, we obtain the desired assertion (55).

If the condition (b) is satisfied, then

[mathematical expression not reproducible]. (61)

Taking the limit superior of both sides, we obtain the desired assertion (56). This finishes the proof.

Remark 12. From Theorem 11, we concluded that when the noise is sufficiently small, the value of [R.sup.s.sub.0] which is below 1 will lead to the extinction of the disease. Therefore, we consider the value [R.sup.s.sub.0] as the threshold of stochastic system (6).

Remark 13. If the condition (a) holds, it easy to prove that [R.sup.s.sub.0] < 1; and if the condition (b) holds, we directly have [R.sup.s.sub.0] < 1. Therefore an interesting open problem is whether we can establish the extinction of the disease only when [R.sup.s.sub.0] < 1.

7. Numerical Simulations

In this section, we present the numerical simulations to support the above theoretical results, illustrating extinction and persistence in mean of the disease. We mainly use Milstein's higher-order method to discrete the system (6). Moreover, we numerically simulate the solution of a corresponding deterministic system for comparison.

Example 1. For the deterministic model (4) and its stochastic version (6), the parameters are taken as follows:

A = 0.4, [beta] = 0.2, [mu] = 0.3 [[gamma].sub.1] = 0.01, [[gamma].sub.2] = 0.03, b = 0.2, q = 0.1, a = 0.1, k = 0.2, [[alpha].sub.1] = 0.05, [[alpha].sub.2] = 0.02, [[alpha].sub.3] = 0.01. (62)

For system (4), we get [R.sub.0] = 1.5873 > 1; thus it admits a unique endemic equilibrium [P.sup.*] which is globally stable for any initial values (S(0), I(0), R(0)).

For the corresponding stochastic model (6), we choose [sigma] = 0.5; then we have

[[beta].sup.2]/[2 ([mu] + [[gamma].sub.1] + [[gamma].sub.2] + a - qb)] - [[sigma].sup.2] = -0.2024 < 0. (63)

Then condition (a) of Theorem 11 is satisfied. We conclude that, for any initial value, I(t) obeys

[mathematical expression not reproducible]. (64)

That is, I(t) will tend to zero exponentially with probability one (see Figure 1(a)).

Example 2. We choose the parameter values of our stochastic system (6) as follows:

A = 0.2, [beta] = 0.52, [mu] = 0.3 [[gamma].sub.1] = 0.01, [[gamma].sub.2] = 0.03, b = 0.2, q = 0.01, a = 0.1, k = 0.2, [[alpha].sub.1] = 0.05, [[alpha].sub.2] = 0.02, [[alpha].sub.3] = 0.01, [sigma] = 0.1. (65)

In this case, we have

[mathematical expression not reproducible]. (66)

Then, condition (b) of Theorem 11 is satisfied. We conclude that, for any initial value, I(t) obeys

[mathematical expression not reproducible]. (67)

That is, I(t) will tend to zero exponentially with probability one (see Figure 1(b)).

Example 3. We choose the parameter values of our stochastic system (4) as follows:

A = 0.2, [beta] = 0.8, [mu] = 0.3 [[gamma].sub.1] = 0.1, [[gamma].sub.2] = 0.5, b = 0.2, q = 0.01, a = 0.1, k = 0.2, [[alpha].sub.1] = 0.05, [[alpha].sub.2] = 0.02, [[alpha].sub.3] = 0.01, [sigma] = 0.1. (68)

In this case, we have

[mathematical expression not reproducible]. (69)

Then, according to Theorem 9, the solution of a stochastic system (6) is persistent in mean (see Figures 2 and 3).

8. Conclusion and Future Directions

In this paper, we considered a general stochastic SIRS epidemic model with vertical transmission and transfer from infectious to susceptible. Firstly, we proved that the solution of the stochastic system (6) is positive and bounded. Then, we obtained sufficient conditions for the stochastic stability of disease-free equilibrium by using a suitable Lyapunov function and other techniques of stochastic analysis. Furthermore, we showed that the disease persists when the basic reproduction number [R.sup.s.sub.0] > 1. Finally, we given tow sufficient conditions for extinction of disease with probability 1.

Besides the white noise perturbation, epidemic models may be perturbed by telegraph noise (or burst noise) which can lead to switch the system from an environmental regime to another [42]. The telegraph noise can be illustrated as a switching between two or more subregimes of different environments [43]. Those switching regimes are often memory less and the waiting time for the next switching follows the exponential distribution [44]. As a matter of fact, the effect of the telegraph noise on the population dynamics has received great attention recently. For example, Li et al. [45] investigated the threshold dynamics and ergodicity of an SIRS epidemic model with Markovian switching. Liu and Zhu [46] studied the stability and ergodicity of a budworm growth model with random perturbations. They showed that both white noises and regime-switching can change the stability of the model greatly. Liu et al. [47] established an interesting results on the dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays. We seek in our future works to see the impact of the telegraph noise on the dynamics of complicated population systems.

https://doi.org/10.1155/2018/7570296

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

[1] W. O. Kermack and A. G. McKendrick, "A contribution to the mathematical theory of epidemics," Proceedings of the Royal Society A Mathematical, Physical and Engineering Sciences, vol. 115, no. 772, pp. 700-721, 1927.

[2] W. O. Kermack and A. G. McKendrick, "Contributions to the mathematical theory of epidemics-I," Bulletin of Mathematical Biology, vol. 53, no. 1-2, pp. 33-55, 1991.

[3] M. Zhien, Z. Yicang, and W. Jianhong, Modeling and Dynamics of Infectious Disease, Higher Education Press, 2009.

[4] R. Anderson and R. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, 1991.

[5] R. Anderson and R. May, Population Biological of Infectious Disease, Springer, Heidelberg, Germany, 1982.

[6] V. Capasso and G. Serio, "A generalization of the Kermack-McKENdrick deterministic epidemic model," Mathematical Biosciences, vol. 42, no. 1-2, pp. 43-61, 1978.

[7] T. Li, F. Zhang, H. Liu, and Y. Chen, "Threshold dynamics of an SIRS model with nonlinear incidence rate and transfer from infectious to susceptible," Applied Mathematics Letters, vol. 70, pp. 52-57, 2017.

[8] C. Zhu, G. Zeng, and Y. Sun, "The threshold of a stochastic SIRS model with vertical transmission and saturated incidence," Discrete Dynamics in Nature and Society, vol. 2017, Article ID 5620301, 9 pages, 2017.

[9] L. J. Allen, "An introduction to stochastic epidemic models," in Mathematical Epidemiology, vol. 1945 of Lecture Notes in Math., pp. 81-130, Springer, Berlin, Germany, 2008.

[10] J. Li and Z. Ma, "Qualitative analyses of SIS epidemic model with vaccination and varying total population size," Mathematical and Computer Modelling, vol. 35, no. 11-12, pp. 1235-1243, 2002.

[11] J. Li and Z. Ma, "Stability analysis for SIS epidemic models with vaccination and constant population size," Discrete and Continuous Dynamical Systems--Series B, vol. 4, no. 3, pp. 635- 642, 2004.

[12] Y. Liu and J.-A. Cui, "The impact of media coverage on the dynamics of infectious disease," International Journal of Biomathematics, vol. 1, no. 1, pp. 65-74, 2008.

[13] A. Lahrouz, L. Omari, and D. Kiouach, "Global analysis of a deterministic and stochastic nonlinear SIRS epidemic model," Nonlinear Analysis: Modelling and Control, vol. 16, no. 1, pp. 59-76, 2011.

[14] A. Kaddar, "On the dynamics of a delayed SIR epidemic model with a modified saturated incidence rate," Electronic Journal of Differential Equations, no. 133, pp. 1-7, 2009.

[15] X. Wang and S. Liu, "An epidemic model with different distributed latencies and nonlinear incidence rate," Applied Mathematics and Computation, vol. 241, pp. 259-266, 2014.

[16] 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, 1986.

[17] A. Korobeinikov, "Lyapunov functions and global stability for SIR and SIRS epidemiological models with non-linear transmission," Bulletin of Mathematical Biology, vol. 68, no. 3, pp. 615-626, 2006.

[18] A. Korobeinikov, "Global properties of infectious disease models with nonlinear incidence," Bulletin of Mathematical Biology, vol. 69, no. 6, pp. 1871-1886, 2007.

[19] A. Korobeinikov and P. K. Maini, "Non-linear incidence and stability of infectious disease models," Mathematical Medicine and Biology, vol. 22, no. 2, pp. 113-128, 2005.

[20] X. Liu and L. Yang, "Stability analysis of an SEIQV epidemic model with saturated incidence rate," Nonlinear Analysis: Real World Applications, vol. 13, no. 6, pp. 2671-2679, 2012.

[21] S. Ruan and W. Wang, "Dynamical behavior of an epidemic model with a nonlinear incidence rate," Journal of Differential Equations, vol. 188, no. 1, pp. 135-163, 2003.

[22] D. Zhao and S. Yuan, "Persistence and stability of the disease-free equilibrium in a stochastic epidemic model with imperfect vaccine," Advances in Difference Equations, vol. 2016, no. 280, 2016.

[23] J. R. Beddington, "Mutual interference between parasites or predators and its effect on searching efficiency," Journal of Animal Ecology, vol. 44, pp. 331-340, 1975.

[24] P. Crowley and E. Martin, "Functional responses and interference within and between year classes of a dragonfly population," Journal of the North American Benthological Society, vol. 8, pp. 211-221, 1989.

[25] B. Du, M. Hu, and X. Lian, "Dynamical behavior for a stochastic predator-prey model with HV type functional response," Bulletin of the Malaysian Mathematical Sciences Society, vol. 40, no. 1, pp. 487-503, 2017.

[26] C. Y. Ji, D. Q. Jiang, and N. Z. Shi, "Multigroup SIR epidemic model with stochastic perturbation," Physica A: Statistical Mechanics and Its Applications, vol. 390, no. 10, pp. 1747-1762, 2011.

[27] Y. Zhou, W. Zhang, and S. Yuan, "Survival and stationary distribution of a SIR epidemic model with stochastic perturbations," Applied Mathematics and Computation, vol. 244, pp. 118-131, 2014.

[28] F. Rao, W. Wang, and Z. Li, "Stability analysis of an epidemic model with diffusion and stochastic perturbation," Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 6, pp. 2551-2563, 2012.

[29] Z. Chang, X. Meng, and X. Lu, "Analysis of a novel stochastic SIRS epidemic model with two different saturated incidence rates," Physica A: Statistical Mechanics and Its Applications, vol. 472, pp. 103-116, 2017.

[30] E. Tornatore, S. M. Buccellato, and P. Vetro, "Stability of a stochastic SIR system," Physica A: Statistical Mechanics and Its Applications, vol. 354, no. 1-4, pp. 111-126, 2005.

[31] A. Gray, D. Greenhalgh, L. Hu, X. Mao, and J. Pan, "A stochastic differential equation SIS epidemic model," SIAM Journal on Applied Mathematics, vol. 71, no. 3, pp. 876-902, 2011.

[32] P. J. Witbooi, "Stability of an SEIR epidemic model with independent stochastic perturbations," Physica A: Statistical Mechanics and Its Applications, vol. 392, no. 20, pp. 4928-4936, 2013.

[33] Q. Liu, D. Jiang, N. Shi, T. Hayat, and A. Alsaedi, "Nontrivial periodic solution of a stochastic non-autonomous SISV epidemic model," Physica A: Statistical Mechanics and Its Applications, vol. 462, pp. 837-845, 2016.

[34] F. Wei and J. Liu, "Long-time behavior of a stochastic epidemic model with varying population size," Physica A: Statistical Mechanics and Its Applications, vol. 470, pp. 146-153, 2017.

[35] C. Ji and D. Jiang, "Threshold behaviour of a stochastic SIR model," Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems, vol. 38, no. 21-22, pp. 5067-5079, 2014.

[36] Y. Zhao and D. Jiang, "The threshold of a stochastic SIRS epidemic model with saturated incidence," Applied Mathematics Letters, vol. 34, pp. 90-93, 2014.

[37] Y. Zhao and D. Jiang, "The threshold of a stochastic SIS epidemic model with vaccination," Applied Mathematics and Computation, vol. 243, pp. 718-727, 2014.

[38] Y. Zhao and D. Jiang, "Dynamics of stochastically perturbed SIS epidemic model with vaccination," Abstract and Applied Analysis, vol. 2013, Article ID 517439, 12 pages, 2013.

[39] Y. Zhao, D. Jiang, X. Mao, and A. Gray, "The threshold of a stochastic sirs epidemic model in a population with varying size," Discrete and Continuous Dynamical Systems--Series B, vol. 20, no. 4, pp. 1289-1307, 2015.

[40] X. Mao, Stochastic Differential Equations and Applications, Horwood, Chichester, UK, 2nd edition, 1997.

[41] V. N. Afanas'ev, V. B. Kolmanovskii, and V. R. Nosov, Mathematical Theory of Control Systems Design, Springer, Dorderecht, The Netherlands, 1996.

[42] Q. Liu, D. Jiang, and N. Shi, "Threshold behavior in a stochastic SIQR epidemic model with standard incidence and regime switching," Applied Mathematics and Computation, vol. 316, pp. 310-325, 2018.

[43] D. Greenhalgh, Y. Liang, and X. Mao, "Modelling the effect of telegraph noise in the SIRS epidemic model using Markovian switching," Physica A: Statistical Mechanics and Its Applications, vol. 462, pp. 684-704, 2016.

[44] Q. Luo and X. Mao, "Stochastic population dynamics under regime switching," Journal of Mathematical Analysis and Applications, vol. 334, no. 1, pp. 69-84, 2007.

[45] D. Li, S. Liu, and J. Cui, "Threshold dynamics and ergodicity of an SIRS epidemic model with Markovian switching," Journal of Differential Equations, vol. 263, no. 12, pp. 8873-8915, 2017.

[46] M. Liu and Y. Zhu, "Stability of a budworm growth model with random perturbations," Applied Mathematics Letters, vol. 79, pp. 13-19, 2018.

[47] M. Liu, X. He, and J. Yu, "Dynamics of a stochastic regime-switching predator-prey model with harvesting and distributed delays," Nonlinear Analysis: Hybrid Systems, vol. 28, pp. 87-104, 2018.

Driss Kiouach (iD) and Yassine Sabbar

MSTI Team, High School of Technology, Ibn Zohr University, Agadir, Morocco

Correspondence should be addressed to Driss Kiouach; d.kiouach@uiz.ac.ma

Received 22 December 2017; Revised 14 March 2018; Accepted 25 March 2018; Published 6 May 2018

Academic Editor: Zhengqiu Zhang

Caption: Figure 1: Computing simulation of the paths I(t) for the SIRS epidemic model (6) and I(t) of the corresponding deterministic system (4).

Caption: Figure 2: Computing simulation of the paths (S(t), I(t)) for the SIRS epidemic model (6).

Caption: Figure 3: Computing simulation of the paths R(t) for the SIRS epidemic model (6).

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Kiouach, Driss; Sabbar, Yassine |

Publication: | Discrete Dynamics in Nature and Society |

Date: | Jan 1, 2018 |

Words: | 6349 |

Previous Article: | The Car Following Model with Relative Speed in Front on the Three-Lane Road. |

Next Article: | Nash Equilibrium Strategy for a DC Pension Plan with State-Dependent Risk Aversion: A Multiperiod Mean-Variance Framework. |

Topics: |