International Journal of Modern Nonlinear Theory and Application
Vol.05 No.04(2016), Article ID:72114,14 pages
The Dynamic Behavior of a Discrete Vertical and Horizontal Transmitted Disease Model under Constant Vaccination
Mingshan Li, Xiumin Liu, Xiaoliang Zhou*
School of Mathematics and Statistics, Lingnan Normal University, Zhanjiang, China
Copyright © 2016 by authors and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
Received: October 16, 2016; Accepted: November 15, 2016; Published: November 18, 2016
In this paper, a class of discrete vertical and horizontal transmitted disease model under constant vaccination is researched. Under the hypothesis of population being constant size, the model is transformed into a planar map and its equilibrium points and the corresponding eigenvalues are solved out. By discussing the influence of coefficient parameters on the eigenvalues, the hyperbolicity of equilibrium points is determined. By getting the equations of flows on center manifold, the direction and stability of the transcritical bifurcation and flip bifurcation are discussed.
Vertical and Horizontal Transmission, Vaccination, Center Manifold, Transcritical Bifurcation, Flip Bifurcation
The SIR infections disease model is an important model and has been studied by many authors  -  . The basic and important research subjects for these systems are local and global stability of the disease-free equilibrium and the endemic equilibrium, existence of periodic solutions, persistence and extinction of the disease, etc. In recent years, the study of vaccination, treatment, and associated behavioral changes related to disease transmission has been the subject of intense theoretical analysis      . In 2008, Meng and Chen  considered a class of continuous vertical and horizontal transmitted epidemic model under constant vaccination
where S represents the proportion of individuals susceptible to the disease, who are born (with b) and die (with d) at the same rate b (b = d) and have mean life expectancy. The susceptible become infectious at a bilinear rate, where I is the proportion of infectious individuals and is the contact rate. The infectious recover (i.e. acquire lifelong immunity) at a rate r, so that is the mean infectious period. The constant p, q, 0 < p < 1, 0 < q < 1, and p + q = 1, where p is the proportion of the offspring of infective parents that are susceptible individuals, and q is the proportion of the offspring of infective parents that are infective individuals. In their work, the basic reproductive rate determining the stability of disease-free equilibrium point and endemic equilibrium point was found out and the local and global stability of the equilibrium points have been researched by using Lyapunov function and Dulac function.
Due to a lot of discrete-time models are not trivial analogues of their continuous ones and simple discrete-time models can even exhibit complex behavior (see  ), in this paper, we pay attention to the discrete situation of Equation (1) as follows
where, and represent susceptible, infective and recovered subgroups, n represent a fixed time. Under the hypothesis of population being constant size, the model is transformed into a planar map and its equilibrium points and the corresponding eigenvalues are solved out. By discussing the influence of coefficient parameters on the eigenvalues, we determine the hyperbolicity of equilibrium points. Further, we get the equations of flows on center manifold and discuss the direction and stability of the transcritical bifurcation and flip bifurcation.
2. Hyperbolic and Non-Hyperbolic Cases
In this section, we will discuss the hyperbolic and non-hyperbolic cases in a two parameters space parameter. In view of assumption that population is a constant size, i.e.,
system Equation (2) can be changed into
Rewrite Equation (4) as a planar map F:
It is obvious that this map has a disease-free equilibrium point and an endemic equilibrium point where
Theorem 1. The equilibrium point is non-hyperbolic if and only if lies on the lines:
Otherwise, the equilibrium point is an one of the following types: (See Table 1).
Proof. The Jacobian matrix of map (5) at is:
And its eigenvalues are
From the assumption, we see that. Then non-hyperbolic will be happened in the case. From and, we get that and
lies on. Also, from, we know which means lies on. When (referred to the case), the eigenvalue
satisfies, then the equilibrium point P is a saddle. When
(referred to the case), the eigenvalue satisfie,
so the equilibrium point P is a stable node and meanwhile when (referred to the case), the equilibrium point P is a saddle since. The proof is complete.
Theorem 2. We select s, r as parameters. There does not exist non-hyperbolic case for the equilibrium. But the hyperbolicity can be divided into the following cases (I), (II).
(I) When, there exist six types for hyperbolic equilibrium point Q: (See Table 2).
Table 1. Types of hyperbolic equilibrium point.
Table 2. Types of hyperbolic equilibrium point.
(II) When, there exist four types for hyperbolic equilibrium point Q: (See Table 3).
Proof. Performing a coordinate shift as follows:
and letting denote the transformed F, we translate equilibrium into and discuss equilibrium point of the map. The matrix of linearization of at is
where,. Its eigenvalues are
Table 3. Types of hyperbolic equilibrium.
It is known that is hyperbolic if and only if none of eigenvalues, lies on the unit circle. In the following we discuss the eigenvalues in two case, i.e., and.
When discriminant, then and are both real . Because non-hyperbolicity happens if and only if or. For whether or, we can get. By condition
and, we see that. This is a contradiction with and, so and are impossible. Next, let’s examine and. From
whether or, we can get, By condition we see that, , This is a contra-
diction with, so and are impossible.
When, and are a pair of conjugate complex. Since
Therefore, and lie inside of and the equilibrium point Q is a stable focus referred to the case.
When, the equilibrium point Q Is hyperbolic. If, i.e.
The matrix has a double real eigenvalue. From the constraint condition, it is obvious that. Therefore, equilibrium point Q is stable node in the case of and.
If, i.e., and, the eigenvalue and are different real numbers. We first discuss the case that, i.e., , In this case we have
We have for, On the other hand, there also exists for. In fact, since
We have. Therefore, the equilibrium Q is a stable node as .
For the case, i.e., , we have and
We assume, by condition, we see that
, i.e., and by condition. This is a con-
tradiction with and. So are impossible,
i.e.,. Therefore, we have. Therefore, the equilibrium Q is a stable node as.
Finally, we study the case of,. We have
Then, we have for. Moreover, there also has for. In fact that,
We have. This means that the equilibrium Q is a stable node for.
When discriminant, because non-hyperbolicity happens if and only if or. Similar to the proof in case (I), neither nor is possible.
When, and are a pair of conjugate complex. Since
Therefore, and lie inside of and the equilibrium point Q is a stable node referred to the case.
When, the equilibrium point Q is hyperbolic. If, the matrix has a
double real eigenvalue. From the constraint condition,
it is obvious that. Therefore, equilibrium point Q is stable node in the case of. If, we first discuss the case that, i.e., , In this case we have
We have for, On the other hand, there also exists for
. In fact, since Therefore, we have.
Therefore, the equilibrium Q is a saddle as.
Finally, we study the case of, i.e., We easily prove by same methods as in case (I). This means that the equilibrium Q is a stable node for. The proof is complete.
3. Transcritical Bifurcation of the Model
The following lemmas were be derived from reference  .
Lemma 1. (  , Theorem 2.1.4) The map
satisfies that A is cxc matrix with eigenvalues of modulus one, and B is sxs matrix with eigenvalues of modulus less than one, and
where f and g are () in some neighborhood of the origin. Then there exists a center manifold for equation (7) which can be locally represented as a graph as follows
For sufficiently small. Moreover, the dynamics of equation (4.1) restricted to the center manifold is, for sufficiently small, given by the c-dimensional map
Lemma 2. (  , in page 365) A one-parameter family of () one-dimensional
Having a non-hyperbolic fixed point, i.e.,
Undergoes a transcritical bifurcation at if
Theorem 3. A transcritical bifurcation occurs at the equilibrium when. More concretely, for slightly there are two equilibriums: a stable point P and an unstable negative equilibrium which coalesce at, for slightly there are also two equilibriums: an unstable equilibrium P and a stable positive equilibrium Q. Thus an exchange of stability has occurred at.
Proof. For, we have and. Consider as the bifurcation parameter and write F as to emphasize the dependence on. Performing a coordinate shift as follows,. One can easily see that the matrix is
and it has eigenvectors
Corresponding to and respectively, where T means the transpose of matrices. First, we put the matrix into a diagonal form. Using the eigenvectors (9), we obtain the transformation
which transform system Equation (5) into
Rewrite system (12) in the suspended form with assumption,
Thus, from Lemma 1, the stability of equilibrium near can be determined by studying an one parameter family of map on a center manifold which can be represented as follows,
for sufficiently small v and.
We now want to compute the center manifold and derive the mapping on the center manifold. We assume
near the origin, where means terms of order. By Lemma 1, those coefficients can be determined by the equation
Substituting (16)into (15) and comparing coefficients of and in (15), we get
from which we solve
Therefore, the expression of (15) is approximately determined:
Substituting (17) into (14), we obtain a one dimensional map reduced to the center manifold
It is easy to check that
The condition (19) implies that in the study of the orbit structure near the bifurcation point terms of do not qualitatively affect the nature of the bifurcation, namely they do not affect the geometry of the curves of equilibriums passing through the bifurcation point. Thus, the orbit structure of (18) near is qualitatively the same as the orbit structure near of the map
Map (20) can be viewed as truncated normal form for the transcritical bifurcation (see Lemma 2). The stability of the two branches of equilibriums lying on both sides of are easily verified.
4. Degenerate Flip Bifurcation of the Model
This section is devoted to the analysis for the case. From section 2, we
have for. For this case, degenerate flip
bifurcation happens at the equilibrium point.
Theorem 4. For map (5) when, degenerate flip bifurcation happens at the equilibrium point.
Proof. Performing a coordinate shift as follows
We translate equilibrium into, and letting denote the transformed
Therefore, we discuss equilibrium point of the map. The matrix of linearization of at is
For, considering as the bifurcation parameter and write as to emphasize the dependence on w. Therefore, we have
The matrix have eigenvectors and corresponding
to and. Therefore, by transformation
Therefore, we obtain the inverse of transformation (23)
Therefore can be changed into the maps:
Rewrite system (25) in the suspended form
Equivalently, the suspended system (26) has a two-dimensional center manifold of the form
Near the origin, where means terms of order. By Lemma 1, those coefficients can be determined by the equation
Comparing coefficients of, and in (27), we get
from which we solve
Thus, the expression of (27)is determined, i.e.,
Substituting (30) into the first equation in (26), we obtain a one-dimensional map, where
From (31), we can check that
Thus, the conditions and of Theorem 3.5.1 in  are not satisfied. Therefore, this is a degenerate flip bifurcation.
Due to a lot of discrete-time models are not trivial analogues of their continuous ones and simple discrete-time models can even exhibit complex behavior (see  ), motivated mainly by Meng and Chen  considering a class of continuous vertical and horizontal transmitted epidemic model (1) under constant vaccination, we study a class of discrete vertical and horizontal transmitted disease model (2) under constant vaccination. By detailed studies, we found discrete model (2) has a flip bifurcation which did not occurred for continuous model. However, the result of flip bifurcation in current paper is a degenerate situation, for which the more in-depth research needs to be continued.
This work has been supported by the Innovation and Developing School Project of Department of Education of Guangdong province (Grant No. 2014KZDXM065) and the Key project of Science and Technology Innovation of Guangdong College Students (Grant No. pdjh2016a0301).
Cite this paper
Li, M.S., Liu, X.M. and Zhou, X.L. (2016) The Dynamic Behavior of a Discrete Vertical and Horizontal Transmitted Disease Model under Constant Vaccination. International Journal of Mo- dern Nonlinear Theory and Application, 5, 171-184. http://dx.doi.org/10.4236/ijmnta.2016.54017
- 1. Piyawong, W., Twizell, E.H. and Gumel, A.B. (2003) An Unconditionally Convergent Finite-Difference Scheme for the SIR Model. Applied Mathematics and Computation, 146, 611-625.
- 2. Pourabbas, E., d’Onofrio, A. and Rafanelli, M. (2001) A Method to Estimate the Incidence of Communicable Diseases under Seasonal Fluctuations with Application to Cholera. Applied Mathematics and Computation, 118, 161-174.
- 3. Beretta, E. and Takeuchi, Y. (1997) Convergence Results in SIR Epidemic Model with Varying Population Sizes. Nonlinear Analysis, 28, 1909-1921.
- 4. Meng, X., Chen, L. and Song, Z. (2007) The Global Dynamics Behaviors for a New Delay SEIR Epidemic Disease Model with Vertical Transmission and Pulse Vaccination. Applied Mathematics and Computation, 28, 1259-1271.
- 5. Allen, L.J.S. (1994) Some Discrete-Time SI, SIR and SIS Epidemic Models. Mathematical Biosciences, 124, 83-105.
- 6. Ma, Z. Zhou, Y. Wang, W. and Jin, Z. (2004) Mathematical Modelling and Research of Epidemic Dynamical Systems (in Chinese). Science Press, Beijing.
- 7. Zhou, X., Li, X. and Wang, W.S. (2014) Bifurcations for a Deterministic SIR Epidemic Model in Discrete Time. Advances in Difference Equations, 168.
- 8. Liao, X., Wang, H., Huang, X., Zeng, W. and Zhou, X. (2015) The Dynamic Properties of a Deterministic SIR Epidemic Model in Discrete Time. Applied Mathematics, 6, 1665-1675.
- 9. Meng, X., Chen, L. and Cheng, H. (2007) Two Profitless Delays for the SEIRS Epidemic Disease Model with Nonlinear Incidence and Pulse Vaccination. Applied Mathematics and Computation, 186, 516-529.
- 10. Agur, Z.L., et al. (1993) Pulse Mass Measles Vaccination across Age Cohorts. Proceedings of the National Academy of Sciences of the USA, 90, 11698-11702.
- 11. Shulgin, B., et al. (1998) Pulse Vaccination Strategy in the SIR Epidemic Model. Bulletin of Mathematical Biology, 60, 1-26.
- 12. Meng, X., Chen, L. and Song, Z. (2007) The Global Dynamics Behaviors for a New Delay SEIR Epidemic Disease Model with Vertical Transmission and Pulse Vaccination. Applied Mathematics and Mechanics (English Edition), 28, 1259-1271.
- 13. Meng, X. and Chen, L. (2008) The Dynamics of a New SIR Epidemic Model Concerning Pulse Vaccination Strategy. Applied Mathematics and Computation, 197, 582-597.
- 14. Anderson, R.M. and May, R.M. (1991) Infections Diseases of Humans: Dynamics and Control. Oxford University Press, Oxford.
- 15. Wiggins, S. (1990) Introduction to Applied Nonlinear Dynamical Systems and Chaos. Springer, New York.
- 16. Guckenheimer, J. and Holmes, P. (1983) Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields. Springer, New York.