Research Article | Open Access

Zhen Wang, Xinhe Wang, "Stability and Hopf Bifurcation Analysis of a Fractional-Order Epidemic Model with Time Delay", *Mathematical Problems in Engineering*, vol. 2018, Article ID 2308245, 8 pages, 2018. https://doi.org/10.1155/2018/2308245

# Stability and Hopf Bifurcation Analysis of a Fractional-Order Epidemic Model with Time Delay

**Academic Editor:**Gisele Mophou

#### Abstract

A fractional-order epidemic model with time delay is considered. Firstly, stability of the disease-free equilibrium point and endemic equilibrium point is studied. Then, by choosing the time delay as a bifurcation parameter, the existence of Hopf bifurcation is studied. Finally, numerical simulations are given to illustrate the effectiveness and feasibility of theoretical results.

#### 1. Introduction

Mathematical model plays an important role in describing the dynamics of biological system [1–3]. The dynamics of the epidemic models have received much attention during the recent years, and to explain the disease spreading and control strategies a series of epidemic models [4–7] was proposed. A stochastic SIRS epidemic model was formulated in [8]; it investigated the effect of stochastic environmental variability on interpandemic transmission dynamics of influenza A. In [9], an age-structured SEIR epidemic model was considered. The authors investigated an SEIR model with varying population size and vaccination strategy in [10], and different threshold parameters were obtained to govern the disease eradication. Many models in biological mathematics involve some time delays. In biological dynamics, time delay was widely applied to reflect some biological facts, such as immunity period [11] and latent period of the disease [12]. An epidemic model with time delay was proposed in [13], and the model is shown as follows:where , represent the number of susceptible and infected population. represents the intrinsic birth rate constant, represents carrying capacity of susceptible population, represents the force of infection or the rate of transmission, represents immigration coefficient, represents death coefficient of , and is the latent period of the disease.

Fractional calculus is a generalization of classical differentiation and integration to arbitrary (noninteger) order [14]. In the past decades, fractional-order calculus garnered considerable attention and it was applied to various fields [15–21]. Recently, many investigators started to study the fractional-order biological models [22–24]. The main reason is that fractional-order models are naturally related to systems with memory which exists in most biological systems [25, 26]. In [27], the authors introduced a fractional-order prey-predator model and deal with the biological behaviors of the model. A fractional-order SIS model with variable population size is considered in [28], and the stability of equilibrium points is studied. A fractional-order model of two-species facultative mutualism with harvesting was presented in [29], and stability of the model was analyzed. In [30], the authors introduced a fractional-order epidemic model with vaccination; it shows that the stability region of the model is related to threshold-value and value of the fractional-order . A delayed fractional-order differential model of HIV infection of CD4^{+} was investigated in [31]. In [32], a fractional-order prey-predator model with time delay and Monod-Haldane function was studied.

In this paper, a fractional-order epidemic model with time delay is studied. We investigate stability and bifurcation of the model with respect to basic reproduction number , fractional-order and time delay . We provide theoretical analysis, using the eigenvalues method and linearization techniques and bifurcation method. The model is depicted as follows:where , , and is in the sense of Caputo fractional derivatives. , , and .

The corresponding linearized system of (2) at any equilibrium point is defined asTaking Laplace transform [33] on both sides of (3), one obtains the characteristic matrix as follows:The properties of eigenvalues of characteristic equation indicate the stability of system (2).

The rest of the paper is organized as follows. In Section 2, some necessary definitions and notions are presented. In Section 3, stability and Hopf bifurcation of the equilibrium point are analyzed. Numerical simulations are given in Section 4 and some conclusions are given in Section 5.

#### 2. Preliminary

There are three main definitions of fractional-order differential, that is, Riemann-Liouville, Grünwald-Letnikov, and Caputo’s definitions. This paper is based on Caputo’s definition.

*Definition 1 (see [34]). *The Caputo fractional derivative with order of a continuous function is defined as follows:where .

Lemma 2 (see [35]). *Considering the fractional differential system with the Caputo derivative,where , , and . The characteristic equation of system (6) is . If the real parts of all the eigenvalues of are negative, then the zero solution to system (6) is locally asymptotically stable.*

Lemma 3 (see [33]). *Considering the fractional delayed differential system with the Caputo derivative,where , , , and . The characteristic equation of the system (7) is . If all the roots of the characteristic equation have negative real parts, then the zero solution of system (7) is locally asymptotically stable.*

#### 3. Main Results

##### 3.1. Basic Production Number and the Existence of the Equilibrium Point

Following from [16], system (2) has a disease-free equilibrium point and the basic reproduction number for the model is . Endemic equilibrium point is . Obviously, . Then we know that model (2) has an endemic equilibrium point when .

##### 3.2. Stability of the Disease-Free Equilibrium Point

Theorem 4. *The disease-free equilibrium point of system (2) is locally asymptotically stable if .*

*Proof. *The characteristic matrix of system (3) evaluated at the equilibrium point isand the characteristic equation isLet ; we can rewrite (9) asClearly, , . When , we get . According to Lemma 2, the disease-free equilibrium point is locally asymptotically stable. This completes the proof.

##### 3.3. Stability of the Endemic Equilibrium Point

The characteristic matrix of system (3) evaluated at the equilibrium point isfrom which we have the characteristic equationwhere

Theorem 5. *When , the endemic equilibrium point of system (2) is locally asymptotically stable if .*

*Proof. *Let ; we can rewrite (12) asIf , one obtains . Obviously, the two roots of (14) are negative. According to Lemma 2, the endemic equilibrium point is locally asymptotically stable. This completes the proof.

When , (12) can be rewritten as whereAssume that (15) has a pair of pure imaginary roots , and then substitute into equation (15); one obtainsand then separating the real and imaginary parts of (17) one hasSquaring and adding the two equations in (18), we obtainDenote , where ; therefore (19) has one positive root at least. If , , , are the roots of , we assume is positive. Substituting into (18), one obtainsSquaring and adding the two equations in (20), one obtainsFrom (21), can be obtained

Theorem 6. *When and , the endemic equilibrium point of system (2) is locally asymptotically stable if and unstable if , where .*

*Proof. *Let and ; (12) reduces toDenote ; one hasDifferentiating both sides of (24) with respect to , we obtainDifferentiating both sides of (24) with respect to yieldsFrom (26), one obtainsWith , the above equality becomesWhen , one gets . Hence, the endemic equilibrium point of system (2) is locally asymptotically stable if and unstable if . This completes the proof.

*Remark 7. *It is worth noting that there will be some future directions to apply our main results to more complex ones like models with time varying delay [36] and models with perturbed parameters [37] or to study the Hopf bifurcation of models with discrete and distributed delays [38].

#### 4. Numerical Simulations

In this section, we give some numerical simulations for system (2) by using the method mentioned in [39, 40].

In Figure 1, we select parameters as , , , , , and , with initial conditions , . After calculation, one obtains disease-free equilibrium point and . In (a), we take , and in (b) we take . According to Theorem 4, the disease-free equilibrium point of system (2) is locally asymptotically stable when .

**(a)**

**(b)**

In Figure 2, the selected parameters are , , , , , , and , with initial conditions , After calculation, one obtains endemic equilibrium and . According to Theorem 5, if and , the endemic equilibrium point of system (2) is locally asymptotically stable. The numerical simulation results are shown in Figure 2.

In Figures 3(a) and 3(b), we plotted the effect measure of immigration coefficient on susceptible and infected populations. The selected parameters are same as Figure 2 with initial conditions , . Values of are shown in the legend. From Figures 3(a) and 3(b), we observe that the number of susceptible individuals increases as increases at the beginning but is finally stable at the same fixed value. The number of infected individuals increases as increases. It shows that after the endemic formation, the number of the susceptible individuals increases as the number of floating population increases in the short term, but in the long run the number of susceptible individuals is the same, and only the number of infected individuals increases.

**(a)**

**(b)**

Figure 4 depicts the Hopf bifurcation of the endemic equilibrium. The parameters are taken as , , , , and , with initial conditions , . After calculation, one obtains , , and . When , is calculated. In (a), we let , and in (b) . (a) and (b) show Hopf bifurcation occurs at . Then one selects different order , we get different time delay, and the results are shown in (c). Figure 4(c) shows that as the value of becomes smaller, the stability domain becomes larger. When , , and , Figure 4(d) shows that the endemic equilibrium point becomes stable.

**(a)**

**(b)**

**(c)**

**(d)**

#### 5. Conclusion

In this paper, a fractional-order epidemic model with time delay is studied and stability and bifurcation of the model are analyzed. The results show that when , the disease-free equilibrium point is locally asymptotically stable for . And we get that when and , the endemic equilibrium point is locally asymptotically stable. According to Theorem 6, when and , the stability of the endemic equilibrium point changes at bifurcation point . Some numerical simulations are given to verify the correctness of the theory, and stability region of model is related to the value of , , and fractional-order .

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare no conflicts of interest.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (nos. 61573008, 61473178), the Natural Science Foundation of Shandong Province (no. ZR201709230160), Post-Doctoral Applied Research Projects of Qingdao (no. 2016115), and SDUST Research Fund (no. 2014TDJH102).

#### References

- 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. View at: Publisher Site | Google Scholar - X. Meng and L. Zhang, “Evolutionary dynamics in a Lotka-Volterra competition model with impulsive periodic disturbance,”
*Mathematical Methods in the Applied Sciences*, vol. 39, no. 2, pp. 177–188, 2016. View at: Publisher Site | Google Scholar - X.-z. Meng, S.-n. Zhao, and W.-y. Zhang, “Adaptive dynamics analysis of a predator-prey model with selective disturbance,”
*Applied Mathematics and Computation*, vol. 266, pp. 946–958, 2015. View at: Publisher Site | Google Scholar | MathSciNet - J.-J. Cheng, Y. Liu, B. Shen, and W.-G. Yuan, “An epidemic model of rumor diffusion in online social networks,”
*The European Physical Journal B*, vol. 86, no. 1, Art. 29, 7 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - Z. Bai and Y. Zhou, “Global dynamics of an {SEIRS} epidemic model with periodic vaccination and seasonal contact rate,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 3, pp. 1060–1068, 2012. View at: Publisher Site | Google Scholar | MathSciNet - A. d'Onofrio, “Stability properties of pulse vaccination strategy in {SEIR} epidemic model,”
*Mathematical Biosciences*, vol. 179, no. 1, pp. 57–72, 2002. View at: Publisher Site | Google Scholar | MathSciNet - X. Zhou and J. Cui, “Analysis of stability and bifurcation for an {SEIR} epidemic model with saturated recovery rate,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 11, pp. 4438–4450, 2011. View at: Publisher Site | Google Scholar | MathSciNet - D. Li, J. Cui, M. Liu, and S. Liu, “The evolutionary dynamics of stochastic epidemic model with nonlinear incidence rate,”
*Bulletin of Mathematical Biology*, vol. 77, no. 9, pp. 1705–1743, 2015. View at: Publisher Site | Google Scholar | MathSciNet - X.-Z. Li, G. Gupur, and G.-T. Zhu, “Threshold and stability results for an age-structured {SEIR} epidemic model,”
*Computers & Mathematics with Applications. An International Journal*, vol. 42, no. 6-7, pp. 883–907, 2001. View at: Publisher Site | Google Scholar | MathSciNet - C. Sun and Y.-H. Hsieh, “Global analysis of an {SEIR} model with varying population size and vaccination,”
*Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems*, vol. 34, no. 10, pp. 2685–2697, 2010. View at: Publisher Site | Google Scholar | MathSciNet - Z. Jiang, W. Ma, and J. Wei, “Global Hopf bifurcation and permanence of a delayed {SEIRS} epidemic model,”
*Mathematics and Computers in Simulation*, vol. 122, pp. 35–54, 2016. View at: Publisher Site | Google Scholar | MathSciNet - T. Zhang, X. Meng, and T. Zhang, “Global analysis for a delayed {SIV} model with direct and environmental transmissions,”
*Journal of Applied Analysis and Computation*, vol. 6, no. 2, pp. 479–491, 2016. View at: Google Scholar | MathSciNet - G. Ranjith Kumar, K. Lakshmi Narayan, and B. Ravindra Reddy, “Stability and Hopf bifurcation analysis of SIR epidemic model with time delay,”
*ARPN Journal of Engineering and Applied Sciences*, vol. 11, no. 3, pp. 1419–1423, 2016. View at: Google Scholar - I. Podlubny,
*Fractional Differential Equations : an Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications*, Academic Press, San Diego, Calif, USA, 1999. - Z. Wang, X. Wang, Y. Li, and X. Huang, “Stability and Hopf Bifurcation of Fractional-Order Complex-Valued Single Neuron Model with Time Delay,”
*International Journal of Bifurcation and Chaos*, vol. 27, no. 13, 1750209, 13 pages, 2017. View at: Publisher Site | Google Scholar | MathSciNet - L. Ferrara and D. Guegan, “Fractional seasonality: Models and Application to Economic Activity in the Euro Area,”
*Aps March Meeting*, vol. 74, no. 1, pp. 2400–2403, 2006. View at: Google Scholar - T. J. Freeborn, “A survey of fractional-order circuit models for biology and biomedicine,”
*IEEE Journal on Emerging and Selected Topics in Circuits and Systems*, vol. 3, no. 3, pp. 416–424, 2013. View at: Publisher Site | Google Scholar - Y. Fan, X. Huang, Z. Wang, and Y. Li, “Nonlinear dynamics and chaos in a simplified memristor-based fractional-order neural network with discontinuous memductance function,”
*Nonlinear Dynamics*, vol. 93, no. 2, pp. 611–627, 2018. View at: Publisher Site | Google Scholar - X. Huang, Y. Fan, J. Jia, Z. Wang, and Y. Li, “Quasi-synchronisation of fractional-order memristor-based neural networks with parameter mismatches,”
*IET Control Theory & Applications*, vol. 11, no. 14, pp. 2317–2327, 2017. View at: Publisher Site | Google Scholar | MathSciNet - Z. Bai, Y. Chen, H. Lian, and S. Sun, “On the existence of blow up solutions for a class of fractional differential equations,”
*Fractional Calculus and Applied Analysis*, vol. 17, no. 4, pp. 1175–1187, 2014. View at: Publisher Site | Google Scholar - Q. Song, X. Dong, Z. Bai, and B. Chen, “Existence for fractional Dirichlet boundary value problem under barrier strip conditions,”
*Journal of Nonlinear Sciences and Applications. JNSA*, vol. 10, no. 7, pp. 3592–3598, 2017. View at: Publisher Site | Google Scholar | MathSciNet - E. Ahmed and A. S. Elgazzar, “On fractional order differential equations model for nonlocal epidemics,”
*Physica A: Statistical Mechanics and its Applications*, vol. 379, no. 2, pp. 607–614, 2007. View at: Publisher Site | Google Scholar - E. Ahmed, A. M. El-Sayed, and H. A. El-Saka, “Equilibrium points, stability and numerical solutions of fractional-order predator-prey and rabies models,”
*Journal of Mathematical Analysis and Applications*, vol. 325, no. 1, pp. 542–553, 2007. View at: Publisher Site | Google Scholar | MathSciNet - I. Ameen and P. Novati, “The solution of fractional order epidemic model by implicit Adams methods,”
*Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems*, vol. 43, pp. 78–84, 2017. View at: Publisher Site | Google Scholar | MathSciNet - Y. Yan and C. Kou, “Stability analysis of a fractional differential model of {HIV} infection of CD4
^{+}T-cells with time delay,”*Mathematics and Computers in Simulation*, vol. 82, no. 9, pp. 1572–1585, 2012. View at: Publisher Site | Google Scholar | MathSciNet - H. A. El-Saka, “Backward bifurcations in fractional-order vaccination models,”
*Journal of the Egyptian Mathematical Society*, vol. 23, no. 1, pp. 49–55, 2015. View at: Publisher Site | Google Scholar | MathSciNet - M. Javidi and N. Nyamoradi, “Dynamic analysis of a fractional order prey-predator interaction with harvesting,”
*Applied Mathematical Modelling: Simulation and Computation for Engineering and Environmental Systems*, vol. 37, no. 20-21, pp. 8946–8956, 2013. View at: Publisher Site | Google Scholar | MathSciNet - H. A. A. El-Saka, “The fractional-order SIS epidemic model with variable population size,”
*Journal of the Egyptian Mathematical Society*, vol. 22, no. 1, pp. 50–54, 2014. View at: Publisher Site | Google Scholar | MathSciNet - N. Supajaidee and S. Moonchai, “Stability analysis of a fractional-order two-species facultative mutualism model with harvesting,”
*Advances in Difference Equations*, Paper No. 372, 13 pages, 2017. View at: Publisher Site | Google Scholar | MathSciNet - D. Rostamy and E. Mottaghi, “Forward and Backward BIFurcation in a Fractional-Order {SIR} Epidemic Model with Vaccination,”
*Iranian Journal of Science & Technology*, vol. 42, no. 2, pp. 663–671, 2018. View at: Publisher Site | Google Scholar | MathSciNet - Z. Liu and P. Lu, “Stability analysis for {HIV} infection of CD4
^{+}T-cells by a fractional differential time-delay model with cure rate,”*Advances in Difference Equations*, 2014:298, 20 pages, 2014. View at: Publisher Site | Google Scholar | MathSciNet - R. Chinnathambi and F. A. Rihan, “Stability of fractional-order prey–predator system with time-delay and Monod–Haldane functional response,”
*Nonlinear Dynamics*, vol. 92, no. 4, pp. 1637–1648, 2018. View at: Publisher Site | Google Scholar - W. Deng, C. Li, and J. Lü, “Stability analysis of linear fractional differential system with multiple time delays,”
*Nonlinear Dynamics*, vol. 48, no. 4, pp. 409–416, 2007. View at: Publisher Site | Google Scholar | MathSciNet - I. Podlubny,
*Fractional Differential Equations*, Academic Press, New York, NY, USA, 1999. - C.-p. Li and Z.-g. Zhao, “Asymptotical stability analysis of linear fractional differential systems,”
*Journal of Shanghai University (English Edition)*, vol. 13, no. 3, pp. 197–206, 2009. View at: Publisher Site | Google Scholar | MathSciNet - S. Jiao, H. Shen, Y. Wei, X. Huang, and Z. Wang, “Further results on dissipativity and stability analysis of Markov jump generalized neural networks with time-varying interval delays,”
*Applied Mathematics and Computation*, vol. 336, pp. 338–350, 2018. View at: Publisher Site | Google Scholar | MathSciNet - J. Wang, K. Liang, X. Huang, Z. Wang, and H. Shen, “Dissipative fault-tolerant control for nonlinear singular perturbed systems with Markov jumping parameters based on slow state feedback,”
*Applied Mathematics and Computation*, vol. 328, pp. 247–262, 2018. View at: Publisher Site | Google Scholar | MathSciNet - L. Li, Z. Wang, Y. Li, H. Shen, and J. Lu, “Hopf bifurcation analysis of a complex-valued neural network model with discrete and distributed delays,”
*Applied Mathematics and Computation*, vol. 330, pp. 152–169, 2018. View at: Publisher Site | Google Scholar | MathSciNet - Z. Wang, “A numerical method for delayed fractional-order differential equations,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 256071, 7 pages, 2013. View at: Publisher Site | Google Scholar | MathSciNet - Z. Wang, X. Huang, and J. Zhou, “A numerical method for delayed fractional-order differential equations: based on G-L definition,”
*Applied Mathematics & Information Sciences*, vol. 7, no. 2L, pp. 525–529, 2013. View at: Publisher Site | Google Scholar | MathSciNet

#### Copyright

Copyright © 2018 Zhen Wang and Xinhe Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.