Global Analysis of an SEIR Epidemic Model with a Ratio-Dependent Nonlinear Incidence Rate

In this paper, a SEIR model with ratio-dependent transmission rate in the form h h h kS I S I α + 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 0 1 R ≤ , the disease-free equilibrium point of the model is globally asymptotically stable and the disease always dies out; if 0 1 R > , the endemic equilibrium point is globally asymptotically stable and the disease persists.


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 ( ) , g S I .
The transmission function ( ) , g S I plays a key role in determining disease dynamics, see [12] [13].Traditionally, the density-dependent transmission (or the bilinear incidence rate, ( ) , g S I kSI = , k 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 ( ) f I 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 where α is a parameter which measures the psychological or inhibitory effect.
It should be noted that if 1 α = and 1 h l = = , (1) becomes the well-known frequency-dependent transmission rate kSI S I + .
In the case of 1 l = , we can obtain the rate-dependent transmission rate ( ) , g S I as the following form: which indicates that the transmission rate of disease is approximately governed by kI if I S is small (e.g., at the beginning of disease's spreading) or it is approximately governed by 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 0 1 R ≤ , the disease-free equilibrium point is globally asymptotically stable, while if 0 1 R > , the disease-free equilibrium point is un- stable 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.

Model Formulation
The whole population ( ) N t is divided into four subclasses based on disease status: the susceptible population, the exposed population, the infected population and removed population, denoted by ( ) ( ) ( ) ( ) . 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: where A is the recruitment rate of the population, µ is the natural death rate, q is the constant rate such that the exposed individuals become infective, γ is the constant rate for recovery, δ is the disease inducing death.Since R does not appear in the first three equations of system (3), thus (3) reduces to the following three-dimensional system (4) ≥ , the non-negativity of the initial value of system (4) in 3 R + is guaranteed, where It follows from system (4) that: 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 0 , 0, 0 for all parameter values.Let ( ) T , , = X E I S , system (4) can be written as (see [17]) And ( ) ) .
The jacobian matrices of ( ) X  and ( ) at the disease-free equilibrium point 0 P are, respectively, ( ) ( ) where 0 , 0 0 So the regeneration matrix of system ( 4) is the spectral radius of Hence the basic reproduction number 0 R (see [18]) is Without difficulty we can get unique endemic equilibrium point state ( ) and the endemic equilibrium point is written in the following form: ( ) Through the above analysis, system (4) has no endemic equilibrium point for 0 1 R ≤ , and from ( 19) we know that system (4) has a unique endemic equilibrium point if 0 1 R > .

Global Stability of the Equilibrium Point
In this section, we first consider the global stability of model ( 4) at the disease-free equilibrium point 0 P .
Theorem 1: If 0 1 R ≤ , the disease-free equilibrium point 0 P is globally asymp- Proof.The characteristic equation of system (4) at 0 P is give by It is clear that λ µ = − is one root of (20).The other roots of ( 20) are determined by the following equation has two roots with negative real parts, therefore, the disease-free equilibrium point 0 P is locally asymptoti- cally stable.
To complete the proof, we construct the following Lyapunov function ( ) The time derivative of 1 L along the solution of ( .
When 0 1 R ≤ , we have By the LaSalle's invariance principle, see [19], we conclude that 0 P is globally asymptotically stable if 0 1 R ≤ .
When 0 1 R > , the Jacobian matrix of model ( 4) evaluated at 0 P is which has an eigenvalue 0 µ − < .Denoted by ( ) ( ) we find that 0 trB < ( trB represents the trace of matrix B, which is the sum of the elements of the main diagonal of the matrix B).When 0 1 R > , and 0 detB < ( detB is the determinant of matrix B), the matrix B must have a positive eigenvalue, thus the disease-free equilibrium point 0 P is unstable whenever 0 1 R > .
This completes the proof.
For the stability of endemic equilibrium point is the solution to the system (4), the upper right derivative of 2 L can be esti- mated: for the other situations, the discussion is similar.
Firstly there is ( ) Integrate from 0 t to t on both sides of (28), we have By Barbalat Lemma in [20], there is ( ) and the unique endemic equilibrium point * P of model ( 4) is globally asymptotically stable.This completes the proof.

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 α + h h h kS I S I , and the regeneration matrix is used to obtain the basic reproductive number 0 R ; the existence of equilibrium is obtained by direct calcula- tion.
By constructing the proper Lyapunov functions, we prove that if 0 1 ≤ R , there exists only the disease-free equilibrium point which is globally asymptotically stable, and if 0 1 > R , there is a unique endemic equilibrium point and this endemic equilibrium point is globally asymptotically stable.

2 0
L P ≠ .Remember that * P set Γ has a boundary, , S E and I must have boundaries, and their derivatives are bounded.It means that 2L is uniformly continuous.