Journal of Applied Mathematics and Physics
Vol.05 No.12(2017), Article ID:80899,9 pages
10.4236/jamp.2017.512188
Global Analysis of an SEIR Epidemic Model with a Ratio-Dependent Nonlinear Incidence Rate
Xiaomei Ren, Tiansi Zhang
College of Science, University of Shanghai for Science and Technology, Shanghai, China
Copyright © 2017 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
Received: October 10, 2017; Accepted: December 4, 2017; Published: December 7, 2017
ABSTRACT
In this paper, a SEIR model with ratio-dependent transmission rate in the form is studied and the basic reproduction number which determines the disease’s extinction or continued existence is obtained. By constructing the proper Lyapunov function, we prove that if , the disease-free equilibrium point of the model is globally asymptotically stable and the disease always dies out; if , the endemic equilibrium point is globally asymptotically stable and the disease persists.
Keywords:
SEIR Model, the Ratio-Dependent Transmission Rate, Basic Reproduction Number, Equilibrium, Global Stability
1. Introduction
As we all know, infectious diseases are enemies of human health. For a long time, people have been fighting various infectious disease; and many methods have been used to study the spread of infectious diseases, so as to control and eliminate infectious diseases, see [1] . Amongst them, mathematical model has become an important tool to analyze the epidemiological characteristics of infectious diseases since the pioneer work of Kermack and McKendrick [2] [3] , which provides us useful control measures in [4] . In standard epidemiological models, the incidence rate (the rate of new infections) is bilinear in the infective and susceptible individuals, see [5] [6] [7] . It has been suggested that the diseases transmission process may have a nonlinear incidence rate, see [8] [9] [10] .
In real life, epidemics tend to have an incubation period, as susceptible to infection after contacting with infected people. First of all, carrying virus, the virus is not immediately, but after a period of time, to onset and into the herd of infected people. In the paper, on the basis of the work of Yuan and Li in [11] an SEIR (Susceptible-Exposed-Infected-Removed) epidemic model is considered with a ratio-dependent nonlinear incident rate .
The transmission function plays a key role in determining disease dynamics, see [12] [13] . Traditionally, the density-dependent transmission (or the bilinear incidence rate, , the proportionality constant) and the frequency-dependent transmission (or the standard incidence rate,
) are two extreme forms of disease transmission, which have
been frequently used in well-know epidemic models [14] [15] . For example, Capasso and Serio [16] introduced a saturated transmission rate , where the infectious force is a function of infectious individuals which has been used in many classic disease model. Especially, Yuan and Li [11] studied a rate-dependent nonlinear incident rate with the following form
(1)
where is a parameter which measures the psychological or inhibitory effect. It should be noted that if and , (1) becomes the well-known
frequency-dependent transmission rate .
In the case of , we can obtain the rate-dependent transmission rate as the following form:
(2)
which indicates that the transmission rate of disease is approximately governed
by if is small (e.g., at the beginning of disease’s spreading) or it is approximately governed by if is large (e.g., in the endemic when
almost everyone is infected). Therefore, the ratio-dependent transmission rate (2) indeed takes accounts of the crowding effects and behavior changes during epidemics. In this paper, we mainly focus on a SEIR epidemic model with the ratio-dependent incidence rate (2).
We consider the global properties of this SEIR model and show that if the basic reproduction number , the disease-free equilibrium point is globally asymptotically stable, while if , the disease-free equilibrium point is unstable and the unique endemic equilibrium point is globally asymptotically stable.
The organization of this paper is as follows: in the next section, we present the model and derive the disease-free equilibrium point and the endemic equilibrium point. In Section 3 we analyze the global stability of the equilibrium point. A brief discussion and summarize are given in Section 4.
2. Model Formulation
The whole population is divided into four subclasses based on disease status: the susceptible population, the exposed population, the infected population and removed population, denoted by , respectively, and . We assume that infectious disease can cause additional mortality, then the SEIR model can be modeled by the following set of nonlinear differential Equations (3) deterministically:
(3)
where is the recruitment rate of the population, is the natural death rate, is the constant rate such that the exposed individuals become infective, is the constant rate for recovery, is the disease inducing death. Since does not appear in the first three equations of system (3), thus (3) reduces to the following three-dimensional system (4)
(4)
Because , , , the
non-negativity of the initial value of system (4) in is guaranteed, where
(5)
It follows from system (4) that:
(6)
Since , the feasible region for system (4) is thus a
bounded set :
(7)
The region is positively invariant with respect to systems (4). So, the only solution with the associated initial conditions will be considered inside the region , where the uniqueness of solutions, usual existence, and continuation results are satisfied. Hence, system (4) is considered mathematically and epidemiologically well posed in . Notice that model (4) has a disease-free equilibrium
point for all parameter values. Let , system (4) can
be written as (see [17] )
(8)
And
(9)
(10)
The jacobian matrices of and at the disease-free equilibrium point are, respectively,
(11)
(12)
where
(13)
(14)
So the regeneration matrix of system (4) is
(15)
the spectral radius of is
(16)
Hence the basic reproduction number (see [18] ) is
(17)
Without difficulty we can get unique endemic equilibrium point state of model (4) with
(18)
and the endemic equilibrium point is written in the following form:
(19)
if .
Through the above analysis, system (4) has no endemic equilibrium point for , and from (19) we know that system (4) has a unique endemic equilibrium point if .
3. Global Stability of the Equilibrium Point
In this section, we first consider the global stability of model (4) at the disease-free equilibrium point .
Theorem 1: If , the disease-free equilibrium point is globally asymptotically stable; if , is unstable.
Proof. The characteristic equation of system (4) at is give by
(20)
It is clear that is one root of (20). The other roots of (20) are determined by the following equation
(21)
If , , thus (21) has two roots with negative real parts, therefore, the disease-free equilibrium point is locally asymptotically stable.
To complete the proof, we construct the following Lyapunov function
(22)
The time derivative of along the solution of (4) is
(23)
When , we have . By the LaSalle’s invariance principle, see
[19] , we conclude that is globally asymptotically stable if .
When , the Jacobian matrix of model (4) evaluated at is
(24)
which has an eigenvalue . Denoted by
(25)
we find that ( represents the trace of matrix B, which is the sum of the elements of the main diagonal of the matrix B). When , and ( is the determinant of matrix B), the matrix B must have a positive eigenvalue, thus the disease-free equilibrium point is unstable whenever . This completes the proof.
For the stability of endemic equilibrium point of model (4), we have the following theorem:
Theorem 2: If , the unique endemic equilibrium point of model (4) is globally asymptotically stable.
Proof. Suppose that is any positive solution to system (4). Define a Lyapunov function
(26)
Obviously and when , . Remember that is the solution to the system (4), the upper right derivative of can be estimated:
(27)
and we all know that when , ; , ; , .
In (27), there are 8 kinds of situation for the size of and , and , and , it is enough to analyze the situation of , , , while for the other situations, the discussion is similar.
Firstly there is
(28)
Integrate from to on both sides of (28), we have
(29)
Since the front set has a boundary, and must have boundaries, and their derivatives are bounded. It means that is uniformly continuous.
By Barbalat Lemma in [20] , there is
(30)
and the unique endemic equilibrium point of model (4) is globally asymptotically stable. This completes the proof.
4. Brief Summary
In this paper, we consider the SEIR epidemic model which is different from the classical nonlinear incident rate. We assume that the infectious force is a function of a number ratio of the infective to that of the susceptible which takes the
form , and the regeneration matrix is used to obtain the basic reproductive number ; the existence of equilibrium is obtained by direct calculation.
By constructing the proper Lyapunov functions, we prove that if , there exists only the disease-free equilibrium point which is globally asymptotically stable, and if , there is a unique endemic equilibrium point and this endemic equilibrium point is globally asymptotically stable.
In the future work, we can further consider adding pulse condition to the model studied in this paper.
Cite this paper
Ren, X.M. and Zhang, T.S. (2017) Global Analysis of an SEIR Epidemic Model with a Ratio-Dependent Nonlinear Incidence Rate. Journal of Applied Mathematics and Physics, 5, 2311-2319. https://doi.org/10.4236/jamp.2017.512188
References
- 1. Ma, Z.E., Zhou, Y.C. and Wang, W.D. (2004) Mathematical Modeling and Research on the Dynamics of Infectiuos Diseases. Science Press, Beijing.
- 2. Kermack, W.O. and McKendrick, A.G. (1927) Cotributions to the Mathematical Theory of Epidemics-I. Proceedings of the Royal Society, 115A, 700-721. https://doi.org/10.1098/rspa.1927.0118
- 3. Cai, Y.L., Kang, Y. and Wang, W.M. (2017) A Stochastic SIRS Epidemic Model with Nonlinear Incidence Rate. Applied Mathematics and Computation, 305, 221-240.https://doi.org/10.1016/j.amc.2017.02.003
- 4. Ma, Z.E. and Zhou, Y.C. (2009) Modeling and Dynamics of Infectiuos Diseases. Higher Education Press, Beijing. https://doi.org/10.1142/7223
- 5. Zhao, Y. (2013) Analysis of a Dynamic Model of Host-Parasite Interaction with Delay and Treatment. Chin Quart J of Math, 28, 118-128.
- 6. Wang, Z.W. and Wang, W.D. (2010) Mathematical Analysis of Immune Response of HIV-I Including Delay. Chin Quarrt J of Math, 25, 45-51.
- 7. Zhang, X.Q. and Zhao, Z. (2010) Differential Susceptibility SIR Epidemic Model with Time Delay and Pulse Vaccination. Chin Quart J of Math, 25, 228-235.
- 8. Zhang, J. and Ma, Z.E. (2003) Global Dynamics of an SEIR Epidemic Model with Saturating Contact Rate. Mathematical Biosciences, 185, 15-32. https://doi.org/10.1016/S0025-5564(03)00087-7
- 9. Korobeinikov, A. and Maini, P.K. (2004) A Lyapunov Function and Global Properties for SIR and SEIR Epidemiological Models with Nonlinear Incidence. Mathematical Bioscience and Engineering, 1, 57-60. https://doi.org/10.3934/mbe.2004.1.57
- 10. Liu, W.M., Hethcote, H.W. and Levin, S.A. (1987) Dynamical Behavior of Epidemiological Models with Nonlinear Incidence Rates. Mathematical Biosciences, 25, 359-380.
- 11. Yuan, S.L. and Li, B. (2009) Global Dynamics of an Epidemic Model with a Ratio-Dependent Nonlinear Incidence Rate. Discrete Dynamics in Nature and Society, 609306, 13.
- 12. Bain, C. (1990) Applied Mathematical Ecology. Journal of Epidemiology and Community Health, 44, 254. https://doi.org/10.1136/jech.44.3.254-b
- 13. Korobeinikov, A. and Maini, P.K. (2005) Nonlinear Incidence and Stability of Infectious Disease Models. Mathematical Medicine and Biology, 22, 113-128. https://doi.org/10.1093/imammb/dqi001
- 14. Fred, B. and Carlos, C.C. (2012) Mathematical Models in Population Biology and Epidemiology. 2nd Edition, Springer, Berlin.
- 15. Hethcote, H.W. (2000) The Mathematica of Infectious Diseases. SIAM Review, 42, 599-653. https://doi.org/10.1137/S0036144500371907
- 16. Capasso, V. and Serio, G. (1978) A Generalizaton of the Kermack-Mckendrick Deterministic Epidemic Model. Mathematical Biosciences, 42, 43-61. https://doi.org/10.1016/0025-5564(78)90006-8
- 17. Driessche, P.V.D. and Wathmough, J. (2002) Reproduction Number and Suc-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48. https://doi.org/10.1016/S0025-5564(02)00108-6
- 18. Bai, Z.G. (2013) The Basic Reproduction Number of the Periodic Infectious Disease Models. Chinese Journal of Engineering Msthematics.
- 19. Elaydy, S. (2004) An Introduction to Difference Equations. Springer, New York, 204-255.
- 20. Min, Y.Y. and Liu, Y.G. (2007) Barbalat Lemma and Its Application in System Stability Analysis. Journal of Shandong University, 37, 51-55.