key: cord-0006422-x4dgseeu authors: Fan, Xiaolin; Wang, Lei; Teng, Zhidong title: Global dynamics for a class of discrete SEIRS epidemic models with general nonlinear incidence date: 2016-05-06 journal: Adv Differ Equ DOI: 10.1186/s13662-016-0846-y sha: 73f483123ea6653fbd59db34fed5f432300b756d doc_id: 6422 cord_uid: x4dgseeu In this paper, a class of discrete SEIRS epidemic models with general nonlinear incidence is investigated. Particularly, a discrete SEIRS epidemic model with standard incidence is also considered. The positivity and boundedness of solutions with positive initial conditions are obtained. It is shown that if the basic reproduction number [Formula: see text] , then disease-free equilibrium is globally attractive, and if [Formula: see text] , then the disease is permanent. When the model degenerates into SEIR model, it is proved that if [Formula: see text] , then the model has a unique endemic equilibrium, which is globally attractive. Furthermore, the numerical examples verify an important open problem that when [Formula: see text] , the endemic equilibrium of general SEIRS models is also globally attractive. As is well known, many infectious diseases possess a latent period, such as Hepatitis, HIV, SARS, Ebola, MERS, etc. When a susceptible individual is infected at the beginning, the disease incubates inside the susceptible for a period of time, then the susceptible becomes an exposed individual before becoming infectious. For such infectious diseases, the resulting model is SEIR (susceptible S, exposed E, infectious I, removed R) epidemic type. The study on SEIR-type epidemic dynamical models is a very important subject in the mathematical theory of epidemiology, and in the last two decades there have been a number of researches on modeling, theoretical analysis, and applications. Continuous SEIR-type epidemic models described by the differential equations have been widely studied. Many important and interesting results can be found in [-] and the references therein. As we all know, it is very difficult to accurately solve a nonlinear differential equation with a given initial condition. Therefore, for many practical requirements, such as numerical calculation, it is often necessary to discretize a continuous model to obtain the corresponding discrete model. At the present time, there are various discretization methods to discretize a continuous model, including the standard methods, such as Euler method, Runge-Kutta method, and some other standard finite difference schemes, and the nonstandard finite difference (NSFD) scheme, which is originally developed by Mickens [-] . In recent years, discrete epidemic models have been widely studied. The basic and important research subjects for these models are the computing of the thresholds values and basic reproduction numbers, the local and global stability of disease-free equilibrium and the endemic equilibrium, the persistence, permanence, and extinction of the disease, and bifurcations and chaos phenomena of the models when some parameters of the models vary, and so on. Many important and interesting results can be found in [-] and the references therein. Particularly, we see that in [, , , -, , ] discrete SI-type epidemic models are investigated, and in [, , , , , ] discrete SIR-type epidemic models are discussed. However, we see that up to now there have been fewer research works on discrete SEIand SEIR-type epidemic models, where the disease has a latent period. Cao and Zhou [] formulated and studied a discrete age-structured SEIT epidemic model, and as an application, discussed the tuberculosis transmission in China. In [], the authors applied Micken's discretization method to obtain a discrete SEIR epidemic model. The positivity of solutions and the existence and stability of equilibrium are discussed. The design of a state observer for the model is tackled. Some sufficient conditions to ensure the asymptotic stability of the observer are provided in terms of a matrix inequality. In [], the authors studied a discrete plant virus disease model with roguing and replanting, which is derived from the continuous case by using the backward Euler method. The basic reproduction number R  is obtained. It is showed that the disease-free equilibrium is globally attractive if R  ≤ , and otherwise, the disease is permanent if R  > . In [, ], the authors proposed a class of discrete SEIS epidemic models with bilinear incidence, which is established from the corresponding continuous SEIS epidemic model by applying the well-known backward difference scheme. The positivity of solutions and the permanence of the model are established. Furthermore, using the Lyapunov function method, the authors proved that if the basic reproduction number R  ≤ , then the disease-free equilibrium is globally asymptotically stable, and if R  > , then the endemic equilibrium exists and is globally asymptotically stable. Consider the following continuous SEIRS epidemic model with general nonlinear incidence: Some particular cases for this model have been investigated in [, , ], where the basic reproduction number is calculated, and the dynamical properties, such as the local and global stability of disease-free equilibrium and endemic equilibrium and the extinction and persistence of the disease are established. Motivated by this work, in this paper, we propose the following discrete SEIRS epidemic model with general nonlinear incidence established by using the backward difference scheme to discretize model (): where X(n) = (S(n), E(n), I(n), R(n)). Our purpose in this paper is to investigate the dynamical behaviors of model (). The basic reproduction number R  is defined. We will prove by using the linearization method and Lyapunov function that if R  ≤ , then disease-free equilibrium is globally attractive, and as a result, the disease is also extinct, and by using the theory of persistence for dynamical systems that if R  > , then the disease is permanent. Furthermore, when model () degenerates into the particular case f (S, E, I, R) = f (S, E, I) and σ = , by constructing the suitable discrete type Lyapunov function we also will prove that if R  > , then model () has a unique endemic equilibrium, which is globally attractive. The organization of this paper is as follows. In Section , the model description and some basic properties are given. Section  deals with the global attractivity of disease-free equilibrium of model (). In Section , the criterion on the permanence of the disease for model () is stated and proved. In Section , the criterion on the global attractivity of the endemic equilibrium for model () in the particular case f (S, E, I, R) = f (S, E, I) and σ =  is stated and proved. Furthermore, in Section , some numerical examples are provided to illustrate the validity of main results obtained in this paper and verify the interesting open problem given in Remark .. Lastly, a discussion is given in Section . In model (), S(n), E(n), I(n), and R(n) denote the numbers of susceptible, exposed, infectious, and recovered classes at nth generation, respectively, is the recruitment rate of the susceptible, μ i (i = , , , ) are the death rates of susceptible, exposed, infectious, and recovered individuals, respectively. Particularly, μ  includes the natural death rate and the disease-related death rate of the infectious class. δ is the translation rate from exposed to infectious, γ is the recovery rate of the infectious individuals, and σ is the rate of losing immunity of the recovered; σ >  indicates that the recovered individuals possess the provisional immunity, and σ =  predicates that the recovered individuals acquire permanent immunity. The incidence rate of the infectious is described by a nonlinear function f (S, E, I, R). In this paper, we always assume that the parameters , μ i (i = , , , ), δ, and γ are positive constants, σ is a nonnegative constant, and μ  ≤ min{μ  , μ  , μ  }. We set For a nonlinear incidence f (S, I), we introduce the following assumption. (H * ) h(S) and g(I) are continuously differentiable with respect to S ≥  and I ≥ , respectively, h(S) is increasing for S ≥ , and g(I) The initial condition for model () is given by We have the following result on the positivity and ultimate boundedness of solutions. Proof We can prove this theorem by using an argument similar to that introduced in [], Theorem .. In fact, we only need to prove by induction that, for any integer n ≥ , if (S(n), E(n), I(n), R(n)) exists and S(n) > , E(n) > , I(n) > , and R(n) ≥ , then (S(n + ), E(n + ), I(n + ), R(n + )) also exists, and S(n + ) > , E(n + ) > , I(n + ) > , and R(n + ) > . From model () by calculating we can obtain S(n + ) = a -bE(n + ), and Let y = E(n + ). By the second equation of model () and by () and (), y satisfies the equation from (H) we obtain that (y) is increasing with respect to y ∈ (, y  ). Then, we obtain Therefore, (y) =  has a unique positive solutionȳ ∈ (, y  ). This shows that E(n + ) exists and E(n + ) =ȳ > . By (), when E(n + ) >  exists, then I(n + ) also exists, and I(n + ) > . By the fourth equation of model () we further have that R(n + ) exists and R(n + ) > . Let x = S(n + ). By the first equation of model () it follows that we obtain that (x) =  has a unique positive solutionx. Therefore, S(n + ) exists, and S(n + ) =x > . From the previous discussions we finally obtain that (S(n + ),E(n + ),I(n + ), R(n + )) exists and is positive. Therefore, solution (S(n), E(n), I(n), R(n)) uniquely exists and is positive for all n > . From model () we have When N() ≤ S  , where S  = μ  , we have N(n) ≤ S  for all n > . In a general way, for any N() >  we can obtain lim sup n→∞ N(n) ≤ S  . Therefore, (S(n), E(n), I(n), R(n)) is ultimately bounded. This completes the proof. Remark . From the previous discussion we see that the region is a positive invariable set for model () and absorbs all nonnegative solutions of model (). Therefore, we can assume in the rest of this paper that The basic reproduction number for model () is given by . Particularly, when f (S, E, I, R) = βh(S)g(I) and f (S, E, I, R) = β SI N , R  becomes of the following forms, respectively, . On the existence of equilibria of model (), we have the following result. Proof Any equilibrium (S, E, I, R) of model () satisfies the equations and Thus, By (H), (I) is decreasing for I > , (I * ) = -(μ  +δ)(μ  +γ ) δ < , and If R  ≤ , then lim I→ + (I) ≤ . Hence, (I) =  has no positive roots. This shows that model () has only a disease-free equilibrium P  . If R  > , then lim I→ + (I) > . Hence, (I) =  has a unique positive root I * . This shows that model () has a unique endemic equilibrium P * (S * , E * , I * , R * ), where This completes the proof. We have the following result on the local stability of the disease-free equilibrium and endemic equilibrium. Theorem . When R  < , the disease-free equilibrium P  of model () is locally asymptotically stable, and when R  > , P  is unstable. Proof The linearization system of model () at equilibrium P  is From the second and third equations of system () we have where Since R  < , we easily prove that two eigenvalues λ i (i = , ) of the matrix A satisfy |λ i | > . Therefore, two eigenvalues ρ i (i = , ) of the matrix A - satisfy |ρ i | < . From the first and fourth equations of system () we have Obviously, the matrix B - has eigenvalues ρ i (i = , ) satisfying |ρ i | < . Therefore, equilibrium (, , , ) of system () is asymptotically stable. Consequently, when R  < , the equilibrium P  of model () is locally asymptotically stable. When R  > , we easily prove that two eigenvalues ρ i (i = , ) of the matrix A - are real numbers and |ρ  | <  and |ρ  | > . Hence, the equilibrium (, ) of system () is unstable. This shows that the equilibrium P  is unstable when R  > . Remark . It is unfortunate that we do not establish the local asymptotic stability of endemic equilibrium P * of model (). In fact, the linearization system of model () at endemic equilibrium P * is In order to obtain the local asymptotic stability of endemic equilibrium P * , we only need to prove that all eigenvalues λ of matrix C - satisfy |λ| < . However, it is a pity that here we do not obtain this. Therefore, when R  > , whether the endemic equilibrium P * of model () also is locally asymptotically stable still is an interesting open problem. In this section, we discuss the global attractivity of disease-free equilibrium of model (). We have the following result. Proof The necessity is obvious because when R  > , model () has an endemic equilibrium P * . Now, we prove the sufficiency. When R  ≤ , we can choose a constant p >  such that Proof The necessity is obvious. In fact, if R  ≤ , then by Theorem . the disease-free equilibrium P  is globally attractive. Now, we prove the sufficiency. When R  > , we can choose constants p >  and ε  >  such that We will use the persistence theory of dynamical systems (see [] , Section . in Chapter ) to prove the theorem. Define the sets It is clear that model () restricted to M ∂ has a globally attractive equilibrium P  (S  , , , ). This shows that {P  } in M ∂ is isolated invariable and acyclic. Now, we prove that for the above ε  > , there is η  >  such that when |S -S  | < η  , E < η  , I < η  , and R < η  , we have We can choose an integer n  >  such that |S(n) -S  | < η  , E(n) < η  , I(n) < η  , and R(n) < η  for all n ≥ n  . Consider the Lyapunov function V (n) = pE(n) + I(n). We have that, for n > n  , = p f S(n + ), E(n + ), I(n + ), R(n + ) -(μ  + δ)E(n + ) + δE(n + ) -(μ  + γ )I(n + ) ≥ p ∂f (S  , , , ) ∂I ε  I(n + ) -(μ  + δ)E(n + ) Hence, we finally have lim n→∞ V (n) = ∞, which leads to a contradiction with lim n→∞ V (n) = . It follows that W s (P  ) ∩ X  = ∅. Thus, by the theorems of uniform persistence for dynamical systems given in [] , we obtain that model () is permanent. This completes the proof. Remark . When f (S, E, I, R) = SI N , by Theorem ., if R  = βδ (μ  +γ )(μ  +δ) > , then the disease in model () is permanent. Remark . Theorem . only obtains the permanence of the disease for model (). However, whether we can also prove that an endemic equilibrium P * is globally attractive for model () when R  > ? In the following section, we will give a partial positive answer. We will prove that, for special case σ =  of model (), an endemic equilibrium P * is globally attractive only when R  > . In this section, we consider a particular case of model ( Let (S(n), E(n), I(n)) be any positive solution of system (). Define the functions Calculating V  (n) = V  (E(n + )) -V  (E(n)), we obtain Using the inequality ln(x) ≤ -x for x < , we have . Example . In model (), we take f (S, E, I, R) = βSI +αI+ωS , = ., μ  = ., μ  = ., μ  = ., β = ., δ = ., ω = ., and γ = .. The parameters μ  , α, and σ will be chosen later. By calculating we have the basic reproduction number R  = . > . We further take μ  = ., α = ., and σ = .. Then the endemic equilibrium P * = (., ., ., .). From the numerical simulations (see Figure  ) we obtain that P * may be globally attractive. Example . In model (), we take f (S, E, I, R) = βSI +αI  , = ., μ  = ., μ  = ., μ  = ., β = ., δ = ., and γ = .. The parameters μ  , α, and σ will be chosen later. By calculating we have the basic reproduction number R  = . > . We further take μ  = ., α = ., and σ = .. Then the endemic equilibrium P * = (., ., ., .). By numerical simulations (see Figure  ) we obtain that P * may be globally attractive. disease-free equilibrium and the permanence of the disease are established, that is, if the basic reproduction number R  ≤ , then the disease-free equilibrium is globally attractive, and if R  > , then the disease is permanent. Furthermore, when the model degenerates into SEIR model, it is proved that when R  > , the model has a unique globally attractive endemic equilibrium. Unfortunately, for SEIRS model (), when the basic reproduction number is greater than one, we do not obtain the local asymptotic stability and global attractivity of the endemic equilibrium. But the numerical examples given in Section  show that the endemic equilibrium for general SEIRS model () may be globally attractive. Therefore, it is still an important and interesting open problem how to apply the linearization method to establish the local asymptotic stability of the endemic equilibrium and how to construct the discrete analogue Lyapunov functions to study the global attractivity of the endemic equilibrium for general SEIRS model (). In addition, the dynamical behaviors for the nonautonomous discrete SEIRS epidemic models, discrete SEIRS epidemic models with vaccination, stage-structured discrete SEIRS epidemic models, and delayed discrete SEIRS epidemic models with nonlinear incidence described by the backward difference scheme are rarely considered. Whether similar results on the permanence and extinction of the disease and the global attractivity of the disease-free equilibrium for these models can be obtained is also an interesting open question. On the other hand, corresponding to continuous model (), we also have the following discrete SEIRS epidemic models with general nonlinear incidence described by the for- = γ I(n + ) -(μ  + σ )R(n + ), with the denominator function φ(h) = e μ  h - μ  , and h >  is the time-step size. An important open problem is whether the results obtained in this paper for model () can also be extended to these models. Global dynamics of an SEIRS epidemic model with periodic vaccination and seasonal contact rate Z: Impulsive vaccination of an SEIRS model with time delay and varying total population size A: Lyapunov functions and global properties for SEIR and SEIS epidemic models A: Global properties of infectious disease models with nonlinear incidence P: Non-linear incidence and stability of infectious disease models Permanence and extinction for a nonautonomous SEIRS epidemic model AB: Global asymptotic properties of an SEIRS model with multiple infectious stages Global dynamics of a class of SEIRS epidemic models in a periodic environment On a nonautonomous SEIRS model in epidemiology Application of Nonstandard Finite Difference Scheme Dynamics consistency: a fundamental principle for constructing nonstandard finite difference scheme for differential equation RE: Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition P: The basic reproduction number in some discrete-time epidemic models Spatial patterns in a discrete-time SIS patch model C: Discrete epidemic models The discrete age-structured SEIT model with application to tuberculosis transmission in China C: Discrete-time SIS models with complex dynamics The existence of codimension-two bifurcation in a discrete SIS epidemic model with standard incidence Global stability for a class of discrete SIR epidemic models A: Discrete-time SIS epidemic model in a seasonal environment AA: Disease-induced mortality in density-dependent discrete-time S-I-S epidemic models AA: Periodically forced discrete-time SIS epidemic model with disease induced mortality JMS: Dynamically-consistent non-standard finite difference method for an epidemic model Complex dynamical behaviors in a discrete eco-epidemiological model with disease in prey Stability analysis in a class of discrete SIRS epidemic models I: Stability analysis and observer design for discrete-time SEIR epidemic models Global analysis of discrete-time SI and SIS epidemic models A discrete plant disease model with roguing and repanting Global stability for a discrete epidemic model for disease with immunity and latency spreading in a heterogeneous host population Y: A discrete-time analogue preserving the global stability of a continuous SEIS epidemic model A: Permanence and global stability of a class of discrete epidemic models HL: Persistence in a discrete-time stage-structured epidemic model M: Permanence of a discrete SIRS epidemic model with delays Z: Global dynamics in a class of discrete-time epidemic models with disease courses Global attractivity of a discrete SIRS epidemic model with standard incidence rate Alee effects in a discrete-time SIS epidemic model with infected newborns The Stability of Dynamical Systems. SIAM Dynamical Systems in Population Biology From (H) we easily obtain that when S(n) = S * ,Since g(x) = x - -ln x >  for x >  and x = , we obtain that when E(n) = E * and I(n) = I * , V  (E(n)) >  and V  (I(n)) > . Computing V  (n) = V  (S(n + )) -V  (S(n)), we have V  (n) = S(n + ) -S(n) -From (H) it follows that, for any η between S(n) and S(n + ),We haveTherefore, from () we obtainf S(n + ), E(n + ), I(n + )Therefore,Choose the Lyapunov functionFor convenience of calculations, we denote S = S(n + ), E = E(n + ), and I = I(n + ).- By calculating we have the basic reproduction number R  = . > . We further take μ  = ., α = ., and σ = .. Then the endemic equilibrium P * = (., ., ., .). By numerical simulations (see Figure  ) we obtain that P * may be globally attractive.Example . In model (), we take f (S, E, I, R) = βS  I (+ωS)(+αI  ) , = ., μ  = ., μ  = ., μ  = ., β = ., δ = ., ω = ., and γ = .. The parameters μ  , α, and σ will be chosen later.By calculating we have the basic reproduction number R  = . > . We further take μ  = ., α = ., and σ = .. Then the endemic equilibrium P * = (., ., ., .). By numerical simulations (see Figure  ) we obtain that P * may be globally attractive.All these examples of numerical simulations show that when R  > , no matter sufficiently greater than one or closer to one but still greater than one, we always obtain that the endemic equilibrium P * is globally attractive, which may offer an affirmative conjecture to the open problem given in Remark ., that is, for the general SEIRS model () the endemic equilibrium P * is globally attractive only when R  > . Therefore, in our future work, we expect to obtain the corresponding theoretical results for this open problem. In this paper, we proposed a discrete SEIRS epidemic model () with general nonlinear incidence, which is described by the backward difference scheme. By our discussions presented in this paper, necessary and sufficient conditions for the global attractivity of the The authors declare that they have no competing interests. All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript. 1 College of Mathematics and System Sciences, Xinjiang University, Urumqi, 830046, People's Republic of China.