Dynamical Behavior and Stability Analysis in a Hybrid Epidemiological-Economic Model with Incubation

and Applied Analysis 3 model system are performed. Conditions for existence of interior equilibrium of model system are studied. Dynamical behavior of model system without incubation around the interior equilibrium is investigated due to variation of economic interest, and state feedback controllers are designed to stabilize model system around the desired interior equilibria. Furthermore, local stability analysis of model system with incubation is analyzed due to variation of time delay; directions of Hopf bifurcation and stability of the bifurcating periodic solutions are also studied. Numerical simulations are made in Section 5, which are utilized to support the theoretical findings obtained in this paper. Finally, this paper ends with a conclusion. 2. Model Formulation In 1954, Gordon [22] proposed the economic theory of a common-property resource, which studies the effect of harvest effort on ecosystem from an economic perspective. In [22], an algebraic equation is proposed to investigate the economic interest of yield of the harvest effort, which takes form as follows: Net Economic Revenue (NER) = Total Revenue (TR) − Total Cost (TC) . (2) Associated with model (1), an algebraic equation, which considers the economic interest V of the harvest effort on recovered host individuals in epidemiological system, that is, R(t), is established as follows: E (t) (wR (t) − c) = V, (3) where E(t) represents the harvest effort on recovered host individuals at time t. V represents the economic interest of harvest effort on the recovered host individuals. w and c represent unit price of harvested population and cost of harvest effort, respectively. Based on (1) and (3), a delayed hybrid model which consists of three differential equations and an algebraic equation can be established as follows: ̇ S (t) = r (1 − S (t) k ) − βS (t) I (t − τ) , ̇ I (t) = βS (t) I (t − τ) − μ 1 I (t) − mI (t) , ̇ R (t) = mI (t) − μ 2 R (t) − E (t) R (t) , 0 = E (t) (wR (t) − c) − V, (4) where S(t), I(t), R(t), E(t), and other parameters share the same interpretations mentioned in (1) and (3), and initial conditions ψ = (ψ 1 , ψ 2 , ψ 3 , ψ 4 ) for model system (4) are defined in the Banach space: {ψ ∈ C + ([−τ, 0] ,R 4 + ) | ψ 1 (θ) = S (θ) , ψ 2 (θ) = I (θ) , ψ 3 (θ) = R (θ) , ψ 4 (θ) = E (θ) } , (5) where R + = {(S, I, R, E) ∈ R : S ≥ 0, I ≥ 0, R ≥ 0, E ≥ 0}. It is also assumed that ψ i (0) > 0 (i = 1, 2, 3, 4) for a biological reason. Model system (4) can be expressed in the following form: Ξ (t) ̇ X (t) = F (X (t)) , (6) where X (t) = (S (t) , I (t) , R (t) , E (t)) T


Introduction
In recent decades, plenty of mathematical models describing the population dynamics of infectious disease have been extensively utilized to understand the transmission mechanism of infectious disease within population ecosystem (see [1][2][3][4] and references therein).Much research efforts have been paid to susceptible-infective-recovered (SIR) vector disease model and corresponding model dynamics (see [5][6][7][8][9][10][11][12] and references therein).Generally, in modelling of communicable disease, the incidence rate (the rate of new infections) is considered to play a vital role in ensuring that the model can provide a reasonable qualitative description of the infectious disease dynamics [3,4].
In order to discuss the spread of an infectious disease transmitted by a vector (e.g., mosquitoes and rats), Takeuchi et al. [7] formulated a delayed SIR epidemic model with a bilinear incidence rate.Beretta et al. [8] considered the global stability of disease free equilibrium and endemic equilibrium of model system; it was shown that the disease free equilibrium is globally stable for any time delay while the endemic equilibrium is not feasible.By constructing a suitable Lyapunov functional, sufficient conditions were derived to guarantee that if the endemic equilibrium is feasible, it is also globally stable for the delay being sufficiently small.Ruan and Wang [13] studied the global dynamics of an SIR model with vital dynamics and nonlinear incidence rate of saturated mass action and global qualitative and bifurcation analyses are carried out.Ma et al. [14] derived an explicit expression of lower bound of the infective individual of solution of model system, which was proposed as an open problem.They therefore gave an estimation of the length of time delay ensuring global asymptotic stability of the endemic equilibrium.Xu and Ma [15] proposed an SIR epidemic 2 Abstract and Applied Analysis model with nonlinear incidence rate and time delay.By analyzing the corresponding characteristic equations, local stability of an endemic equilibrium and a disease free equilibrium are discussed.An SIR model with distributed delay and a general incidence function is studied in McCluskey [9], and the global dynamics for the SIR epidemiological system is analyzed in Zhou and Cui [10].Wang et al. [11] considered the asymptotic behavior of the following SIR vector model: where (), (), and () represent the population density of susceptible, infective, and recovered host individuals at time , respectively.It is assumed that the population growth of susceptible host individuals is governed by the logistic growth with a carrying capacity  > 0 as well as intrinsic birth rate constant  > 0.  > 0 is the average number of constants per infective per unit time and  ≥ 0 denotes the incubation time, and  1 > 0 and  2 > 0 stand for the death rate of infective and recovered host individuals, respectively. > 0 represents the recovery rate of infective host individuals.The local stability of endemic equilibrium is investigated, and conditions for Hopf bifurcation to occur are derived in [11].Along with the line of this research, Enatsu et al. [12] analyze stability of equilibria for a delayed SIR epidemic model, in which population growth is subject to logistic growth in absence of disease and the proposed model with a nonlinear incidence rate satisfying suitable monotonicity conditions.Nowadays, biological resource within ecosystem is commercially harvested and sold with aim of achieving economic interest [16,17].It is well known that harvesting has a strong impact on the dynamic evolution of population and several mathematical models have been established to discuss dynamic effects of harvest effort on population in ecological-epidemiological system, which can be found in [18][19][20][21] and the references therein.The role of harvesting in a predator-prey-parasite system is discussed in [18]; theoretical results show that, using impulsive harvesting effort as control parameter, it is not only possible to control the cyclic behavior of the system populations leading to the persistence of all species but other desired stable equilibrium including disease free can be obtained.A ratio-dependent eco-epidemiological system is proposed in [19] where prey population is subject to harvesting.Positive invariance, boundedness, stability of equilibria, and permanence of system have been established.In [20], an eco-epidemiological model is studied where prey disease is modeled by a susceptible-infective scheme, and the role of harvesting and switching on the dynamics of disease propagation and/or eradication is discussed.An eco-epidemiological model with distributed time delay and impulsive control strategy is investigated in [21]; local stability and complex dynamical behavior are discussed.Under the system of market economy, harvest effort is usually influenced by variation of economic interest of commercial harvesting [16,22].It should be noted that the above mentioned related work [18][19][20][21] only concentrate on the role of harvest effort on population dynamics, while the dynamic effect of economic interest on commercial harvesting and indirect dynamic effect on ecosystem are not considered.The work done in [12] is an extension of [11] with nonlinear incidence rate, while dynamic effect of harvest effort on population dynamics is not considered.
Recently, some hybrid dynamical models are proposed in [23][24][25][26][27][28], which are utilized to discuss the interaction mechanism of harvested ecosystem from an economic perspective.Compared with the traditional mathematical models (differential equations or difference equations) discussing the population dynamics in ecosystem, the hybrid mathematical models proposed in [23][24][25][26][27][28] are made up of differential equations and algebraic equations, where differential equations concentrate on coexistence and interaction mechanism of population and algebraic equations offer a simpler way to study the effect of harvest effort on ecosystem from an economic perspective.Complex dynamical behavior and stability analysis in prey-predator ecosystems with stagestructured population and gestation delay are considered in [23][24][25][26][27][28].In general, differential-algebraic models exhibit more complicated dynamics than ordinary differential models.The differential-algebraic models have been applied widely in power systems, aerospace engineering, chemical processes, social management systems, biological systems, network analysis and oil catalysis, and cracking process (see [29][30][31] and references therein).With the help of differentialalgebraic model for the power systems and bifurcation theory, complex dynamical behaviors of the power systems, especially the bifurcation phenomena that reveal the instability mechanism of power systems have been extensively studied, which can be found in [32][33][34] and the references therein.Furthermore, some applications of differentialalgebraic models in the field of economy, which can be found in [35,36].
It is well known that the recovered host individuals are naturally immune to vector disease [1], and its potential economic interest can be commercially exploited.Furthermore, harvest effort is usually influenced by variation of economic interest of commercial harvesting [16,22] under the system of market economy.Consequently, it is necessary to discuss the coexistence and interaction mechanism of population within harvested epidemiological ecosystem as well as dynamical effect of harvest effort due to variation of economic interest.However, as far as knowledge goes, nobody has explicitly proposed a mathematical model to discuss the dynamic effect of commercial harvest on epidemiological system under the system of market economy.The main objective of this paper is to investigate the transmission mechanism of infectious disease and dynamical effect of commercial harvest on population dynamics, especially the complex dynamical behavior and stability switch due to variation of incubation and commercial harvest economic interest.The organisation of the rest section of this paper is as follows.By introducing commercial harvest effort into model system (1), a hybrid epidemiological-economic model is established in Section 2. Positivity and permanence of solutions of model system are discussed in Section 3. In Section 4, qualitative analyses of model system are performed.Conditions for existence of interior equilibrium of model system are studied.Dynamical behavior of model system without incubation around the interior equilibrium is investigated due to variation of economic interest, and state feedback controllers are designed to stabilize model system around the desired interior equilibria.Furthermore, local stability analysis of model system with incubation is analyzed due to variation of time delay; directions of Hopf bifurcation and stability of the bifurcating periodic solutions are also studied.Numerical simulations are made in Section 5, which are utilized to support the theoretical findings obtained in this paper.Finally, this paper ends with a conclusion.

Model Formulation
In 1954, Gordon [22] proposed the economic theory of a common-property resource, which studies the effect of harvest effort on ecosystem from an economic perspective.In [22], an algebraic equation is proposed to investigate the economic interest of yield of the harvest effort, which takes form as follows: Associated with model ( 1), an algebraic equation, which considers the economic interest V of the harvest effort on recovered host individuals in epidemiological system, that is, (), is established as follows: where () represents the harvest effort on recovered host individuals at time .V represents the economic interest of harvest effort on the recovered host individuals. and  represent unit price of harvested population and cost of harvest effort, respectively.Based on (1) and (3), a delayed hybrid model which consists of three differential equations and an algebraic equation can be established as follows: where (), (), (), (), and other parameters share the same interpretations mentioned in (1) and (3), and initial conditions  = ( 1 ,  2 ,  3 ,  4 ) for model system (4) are defined in the Banach space: where R 4 + = {(, , , ) ∈ R 4 :  ≥ 0,  ≥ 0,  ≥ 0,  ≥ 0}.It is also assumed that   (0) > 0 ( = 1, 2, 3, 4) for a biological reason.
Model system (4) can be expressed in the following form: where Remark 1.The algebraic equation in model system (6) contains no differentiated variables; hence, the leading matrix Ξ() in model system (6) has a corresponding zero row.
Remark 2. The model proposed in [11], which composed of differential equations, only discusses the interaction and coexistence mechanism of susceptible, infective, and recovered host individuals.Compared with the model proposed in [11], algebraic equations are incorporated into the model system (4), which focus on the economic interest of harvesting on recovered host individuals.Hence, the established model not only investigates interaction and coexistence mechanism of population in harvested ecosystem but also studies the dynamical behavior due to the variation of economic interest of commercial harvesting and incubation.

Positivity and Permanence
Theorem 3. Any solutions of model system (4) with initial conditions are positive.
Due to the lemma in [37] and Theorem A.4 in [38], any solution of the model system (4) with positive initial conditions exists uniquely and each component of the solution remains within the interval [0,  0 ) for some  0 > 0. Standard and simple arguments show that solutions of model system (4) always exist and stay positive.Hence, this completes the positivity of the solutions of model system (4).
From a viewpoint of biological and economic interest perspective, persistence of solutions of model system (4) in the case of economic interest V ≥ 0 will be investigated in this section.Some preliminaries are introduced as follows.
Definition 4 (see [39]).Model system (4) is said to be permanent if there exists a compact region Ω 0 ∈ int Ω such that every solution of model system (4) with initial conditions will eventually enter and remain in region Ω 0 .
Definition 5 (see [39]).Consider a metric space  with metric .The distance (, ) of a point  ∈  from a subset  of  is defined by It is further assumed that  is the closure of an open set  0 , and  0 =  0 is nonempty and is the boundary of  0 .Consequently,  0 ∪ 0 = ,  0 ∩ 0 = Ø.We will also suppose that () is a  0 semigroup on  satisfying  () :  0 →  0 , () :  0 →  0 .
According to Theorem 3 and the first and second equation of model system (4), it gives that where   > 0 and   > 0 are independent of corresponding initial values of model system (4). Proof. Denote Next, all conditions in Lemma 6 will be checked.In order to facilitate the proof, we consider the following subsystem of model system (4): where () ≥ 0, () ≥ 0 are continuous on  ∈ [−, 0] and (0) > 0, (0) > 0.
By Definition 5 and model system (16), it is easy to see that  0 and  0 are positively invariant, and conditions (i) and (ii) of Lemma 6 clearly hold.
Since model system (16) possesses two constant solutions in  0 : P0 ∈  1 , P1 ∈  2 with the following form: It follows from simple computation that Moreover, it follows from the second equation of model system (16) that which reveals that all points in  1 approach to P0 ; that is,  1 =   ( P0 ).By using the similar analysis mentioned above, it can be also concluded that all points in  2 approach to P1 ; that is,  2 =   ( P1 ).Based on the above analysis, it shows that invariant sets P0 and P1 are isolated invariant, and { P0 , P1 } is isolated and an acyclic covering.It can be concluded that condition (iii) of Lemma 6 holds.Finally, we will show that   ( P ) ∩  0 = Ø for  = 0, 1.Based on the definition of P0 , it is easy to show that   ( P0 ) ∩  0 = Ø.We will show   ( P1 ) ∩  0 = Ø in the following part.
By the second equation of model system (16), it derives that holds for all  ≥  + .
Consider the following equation: Based on (21) and the comparison principle, it derives that () ≥ () for all  > .
According to the above analysis, all conditions of Lemma 6 hold.By using Lemma 6, it can be obtained that lim inf where   > 0 and   > 0 are independent of the corresponding initial values of model system (4).
Proof .According to Lemmas 8 and 9, it can be obtained that hold for all  > 0, which derive that When the economic interest V = 0, it follows from Theorem 3 and the fourth equation of model system (4) that and Ṙ () = 0. Based on the third equation of model system (4), it can be computed that () = (/)() −  2 .Accordng to (24), it derives that provided that  2 <   .
In the case of V > 0, it derives that () = V/(() − ) based on implicit function theory [41].According to the third equation of model system (4), it can be obtained that which derives that lim inf provided that 0 < V <  2 +   .It follows from Theorem 3 and the third equation of model system (4) that which derives that lim sup Hence, it gives that V/( − ) ≤ () ≤ V/( − ), and it can be rewritten as follows: Based on ( 23), ( 24), (25), and ( 26), it can be concluded that all solutions of model system (4) with initial conditions are persistent in the case of V = 0, and it follows from ( 23), ( 24), ( 28), ( 30) and ( 31) that all solutions of model system (4) with initial conditions are persistent in the case of 0 < V <  2 +   .

Qualitative Analysis of Model System
Dynamical effects of harvest effort and time delay on population dynamics are discussed in this section.It should be noted that the interior equilibrium biologically interprets that susceptible, infective, and recovered host individuals survive as well as harvest on recovered host individuals exists.Bifurcation phenomenon around the interior equilibria can reveal instability mechanism of model system, which are theoretically relevant to infectious disease control and sustainable yield on recovered host individuals in the real world.Consequently, we will mainly concentrate on dynamical behavior and local stability analysis around interior equilibrium of model system (4) in this paper.

Model System without Time Delay.
In this section, dynamical behavior of model system (4) without time delay is investigated, and local stability analysis around the interior equilibrium is discussed due to variation of economic interest of commercial harvesting.Furthermore, state feedback controllers are designed to stabilize model system around the desired interior equilibria in the case of zero economic interest and positive economic interest, respectively.

Singularity Induced Bifurcation
Theorem 11.Model system (4) without time delay has a singularity induced bifurcation around the interior equilibrium, and V = 0 is a bifurcation value.Furthermore, local stability switch occurs as V increases through 0.
Proof.Based on the economic theory of a common-property resource [22], there is a phenomenon of bioeconomic equilibrium in the case of zero harvest economic interest; that is, V = 0.An interior equilibrium can be obtained as follows:  * ( * ,  * ,  * ,  * ), where  * = ( 1 + )/,  * = ( −  1 − )/ 2 ,  * = /, and According to biological interpretation of the interior equilibrium, it follows that  * > 0,  * > 0,  * > 0 and  * > 0. In order to guarantee the existence of interior equilibrium, some inequalities are satisfied: Let V be a bifurcation parameter, () = ((), (), ())  , It can be calculated that trace By virtue of (32), it can be obtained that trace Furthermore, it can be also calculated that It follows from (32) that Based on Section IV(A) in [42], ℎ 3 ((), (), V) can be defined as follows: By simple computing, According to (32), it derives that Based on the above analysis, four items (i-iv) can be obtained as follows.
(i) It is easy to show that   ℎ 2 has a simple zero eigenvalue: and trace( (iii) By virtue of (40), it can be shown that (iv) It is easy to show rank (ℎ 1 ((), (), V)) = 3 and rank(ℎ 2 ((), (), It should be noted that the conditions for singularity induced bifurcation, which is introduced in Section III (A) in [42], consist of three conditions, that is, SI1, SI2, and SI3.According to the above items (i)-(iv), SI1, SI2, and SI3 are all satisfied; hence model (4) without time delay has a singularity induced bifurcation around the interior equilibrium  * and the bifurcation value is V = 0.
Along with the line of the above proof, for model (4) without time delay, it follows from simple computing that It follows from (32) that Inequality (44) satisfies Theorem 3 of [42].According to Theorem 3 of [42], when V increases through 0, one eigenvalue (denoted by  1 ) of model system (4) without time delay moves from C − to C + along the real axis by diverging through infinity; the movement behavior of this eigenvalue influences the stability of model system (4) without time delay.
Since the Jacobian of model system (4) without time delay evaluated around  * takes the following form: according to the leading matrix Ξ() in model system (4) and   * , the characteristic equation of the model system (4) without time delay around By virtue of simple computation, the characteristic equation is as follows: It can be concluded that the rest eigenvalues of model system (4) without time delay (denoted by  2 and  3 ) have negative real parts by using the Routh-Hurwitz criteria [43].It follows from Theorem 3 in [42] that there is only one eigenvalue diverging to infinity as V increases through 0, and the rest eigenvalues are continuous and nonzero and cannot jump from one half open complex plane to another one as V increases through 0. It has been shown that  1 moves from C − to C + along the real axis by diverging through infinity.Therefore,  2 and  3 are continuous and bounded in the C − half plane as V increases through 0 and their movement behaviors have no influence on the stability of model system (4) without time delay around the interior equilibrium  * .
According to Table 1 and the stability theory, it can be concluded that model system (4) without time delay is stable around  * as V < 0 and model system (4) without time delay is unstable around  * as V > 0. Consequently, a stability switch occurs as V increases through 0.
Remark 12.Some preliminaries of singularity induced bifurcation are introduced below.Parameter dependent differential-algebraic hybrid system of the form where (), (), and  have appropriate dimensions.It has been shown recently that there are generically three types of codimension one local bifurcation associated with the differential-algebraic model (48), namely, saddle-node bifurcation, Hopf bifurcation, and singularity induced bifurcation (see [42]).
The singularity induced bifurcation is firstly introduced and analyzed in [42,44].It is a new type of bifurcation and does not occur in usual ordinary differential equation system, which has been characterized for differential-algebraic system, and later improved in [45,46].Roughly speaking, the singularity induced bifurcation refers to a stability change of the differential-algebraic hybrid model (48) owing to some eigenvalues of related linearization ℎ  − ℎ   −1    diverging to infinity when Jacobian   is singular.
One of the important consequences of the singularity induced bifurcation is that it leads to an impulse phenomenon, which may result in the collapse of the differentialalgebraic system (see [45]).More detailed introductions of the singularity induced bifurcation can be found in [42,[44][45][46].Remark 13.It follows from Theorem 11 that there is a phenomenon of singularity induced bifurcation around the interior equilibrium when economic interest increases through zero, which can cause local stability switch of model system (4).As stated in Remark 12, the singularity induced bifurcation can result in impulse phenomenon, which may lead to the collapse of the proposed model.In the harvested epidemiological-economic system, the impulse phenomenon is vividly reflected with the outbreak of infectious disease during a short period in the real world.Under this climate, the infected population will be beyond the carrying capacity of environment, which is disastrous for sustainable development of the harvested ecosystem as well as prosperous yield on recovered host individuals.

State Feedback Controller.
In order to maintain the sustainable yield on recovered host individuals biological resource as well as economic interest of commercial harvesting at an ideal level, some corresponding control strategies should be taken to eliminate the impulse phenomenon caused by singularity induced bifurcation and stabilize model (4) without time delay.In this subsection, state feedback controllers are designed to stabilize model system (4) without time delay around corresponding interior equilibria in the case of V = 0 and V > 0, respectively.
According to the leading matrix Ξ() in model system (4) and   * in (45) (the Jacobian of model system (4) without time delay around the interior equilibrium  * ), it can be calculated that rank (  * , Ξ  * , Ξ 2   * , Ξ 3   * ) = 4.By using Theorem 2-2.1 in [47], it is easy to show that the model system (4) without time delay is locally controllable around the interior equilibrium  * in the case of V = 0. Consequently, a state feedback controller can be applied to stabilize the model system (4) without time delay around  * .By using Theorem 3-1.2 in [47], a state feedback controller () = (() −  * ) ( is a feedback gain and  * is the component of the interior equilibrium  * ) can be applied to stabilize model system (4) without time delay around  * .
Furthermore, the controlled model system (4) without time delay takes the following form: Theorem 14.When economic interest of harvesting is zero, V = 0, if the feedback gain  satisfies the following inequality: then singularity induced bifurcation is eliminated and model system (49) is stable around  * .
Proof.The Jacobian of the model system (49) evaluated at the interior equilibrium  * takes the form According to the leading matrix Ξ() in the model system (4) and J * , the characteristic equation of model system (49) around  * is det(Ξ − J * ) = 0, which can be expressed as follows: where Consequently, if the feedback gain satisfies the above inequality, then model system (49) is stable around  * in the case of zero interest of harvesting. where Based on Routh-Hurwitz criteria [43], (54) has two positive roots if economic interest V satisfies the following inequalities: As analyzed above, there are two interior equilibria (denoted by P * 1 and P * 2 ) when 0 < V * < Ṽ.In this subsection, we only design the controller for the model (4) around the interior equilibrium P * 1 .Some symmetric results about P * 2 can be also obtained, and P * 1 is denoted as P * for simplicity in the following part.Theorem 15.When economic interest of harvesting is positive, 0 < V * < Ṽ, if feedback gain l of controller () = l(() − Ẽ * ) satisfies following inequality: then model system (49) is stable around the interior equilibrium P * ( S * , Ĩ * , R * , Ẽ * ).
Proof.The proof is similar to the proof of Theorem 14 of this paper.
Remark 16.It follows from (55) and Theorem 15 that economic interest of commercial harvesting should be regulated within certain interval V ∈ (0,Ṽ), which guarantees the existence of interior equilibrium in the case of positive economic interest.After applying the state feedback controller into model system (4) without time delay, model system can be stabilized around the corresponding interior equilibrium, respectively.The elimination of the singularity induced bifurcation means the harvested epidemiologicaleconomic system restores to ecological balance and avoidance of infectious disease outbreak.In this subsection, we only investigate dynamical behavior of model system (4) around the interior equilibrium P * 1 .Some symmetric results about the interior equilibrium P * 2 can be also obtained, and P * 1 is denoted as P * for simplicity.According to Jacobian evaluated at the interior equilibrium P * and the leading matrix Ξ() in model system (4), we can obtain the characteristic equation of model system (4) around P * , which can be expressed as follows:

Model
where Now substituting  =  ( is a positive real number) into (58) and separating the real and imaginary parts, two transcendental equations can be obtained as follows: By squaring and adding (60) and (61), it can be calculated that where According to the values of   , ( = 1, 2, 3) and the Routh-Hurwitz criteria [43], a simple assumption that (58) has at least one positive real root  0 is  3 < 0, which derives that  > 3( 1 + ).Hence, under this assumption, (58) will have a pair of purely imaginary roots of the form ± 0 .
In the following part, we determine the direction of motion of  =  0 as  is varied; namely, we determine By differentiating (58) with respect to , it can be obtained that From (62) and the above equation, it can be obtained that Abstract and Applied Analysis 11 According to the values of   , ( = 1, 2, 3) given in (58) of this paper, it is easy to show that Furthermore, if  > 3( 1 + ), then it can be shown that  3 > 0. Hence, it can be concluded that 2 6 0 +  1  4 0 −  3 > 0, which derives sign[d(Re )/d] = *  ,= 0 > 0. Consequently, the transversality condition holds and Hopf bifurcation occurs at  =  0 ,  =  *  .Furthermore, an attracting invariant closed curve bifurcates from interior equilibrium P * when  >  * 0 and ‖ −  * 0 ‖ ≪ 1.
Remark 18. sign[d(Re )/d] = *  > 0 signifies that there exists at least one eigenvalue with positive real part for  =  *  , and the conditions for Hopf bifurcation in [39] are also satisfied yielding the required periodic solution.
In fact, we can choose where  denotes the Dirac delta function.
By virtue of ( 78) and (86), we have where By substituting the corresponding series into (95) and comparing the coefficients, we have Based on the definition of  and ( 97) and (99), it can be obtained that For where ) is a constant vector.Similarly, it follows from (97) and (100) that where  2 = ( (1)  2 ,  (2)  2 ,  (3)  2 ) is a constant vector.Subsequently, values of  1 and  2 should be computed.By using the definition of  and (95), we have where () = (0, ).Based on (95), it derives that in the case of  = 0, which follows that By virtue of (72), it gives that ) .
By virtue of (86), it can be obtained that Then we have According to (107) and (110), we have ) . ( Since  0   is the eigenvalue of (0) and (0) is the corresponding eigenvector, we obtain that where  is identity matrix.

Numerical Simulation
In this section, some numerical simulations are provided to substantiate the theoretical results obtained in Section 4 of this paper.

Numerical Simulation of State Feedback Controller for
Singularity Induced Bifurcation and Local Stability Switch.In this subsection, values of parameters are partially taken from Section 5 of [11] and set in appropriate units, which can be found in Table 2. Numerical simulations are provided to illustrate the effectiveness of the state feedback controllers designed in Section 4.1 in the case of zero economic interest and positive economic interest, respectively.In the case of zero economic interest, model (4) without time delay takes the following form:   By using Theorem 11 of this paper, it can be shown that the model (4) without time delay has a singularity induced bifurcation around the interior equilibrium  * (0.2, 0.08, 0.05, 0.06), and local stability switch occurs as V increases through 0.
Based on the analysis in Section 4.1.2and Theorem 14 of this paper, a state feedback controller () = (() − 0.06) can be applied to stabilize model system (122) around  * , and then the model system (122) with the state feedback controller takes the following form: (123) By using Theorem 14 of this paper, if the feedback gain  satisfies  > 2.2429, then model system (122) is stable around  * and singularity induced bifurcation of the model system (122) is also eliminated.The dynamical responses of model system (122) can be shown in Figure 1.
Furthermore, based on the analysis and inequality (55) in Section 4.1.2,there are two interior equilibria (denoted by P * 1 and P * 2 ) when 0 < V < 0.00668.In the following part, we focus on the case of 0 < V < 0.00668, and the economic interest is set as V * = 0.005 ∈ (0, 0.00668) in appropriate unit, which is arbitrarily selected within the interval (0, 0.00668) and is enough to merit the theoretical analysis obtained in Section 4.1.2.By virtue of the given values of parameters in using Theorem 15 of this paper, if the feedback gain l satisfies l > 3.753, then model system (125) is stable around P * 1 and model system (126) is stable around P * 2 .The dynamical responses of model system (125) and (126) can be shown in Figures 2 and 3, respectively.

Numerical Simulation for Hopf Bifurcation and Local
Stability Switch.In this subsection,values of parameters are partially taken from Section 5 of [11] and set in appropriate units, which can be found in Table 3. Numerical simulations are provided to support the theoretical findings obtained in Section 4.2 of this paper.
Based on the analysis and inequality (55) in Section 4.1.2,there are two interior equilibria (denoted by P * 1 and P * 2 ) when 0 < V < 0.01831.In the following part, we focus on the case of 0 < V < 0.01831, and the economic interest is set as V * = 0.012 ∈ (0, 0.01831) in appropriate unit, which is arbitrarily selected within the interval (0, 0.01831) and is enough to merit the theoretical analysis obtained in Section 4.2.By virtue of the given values of parameters in remains stable for  <  * 0 , and dynamical responses of model system (4) with  = 1.54 are plotted in Figure 4.It should be noted that  = 1.54 in Figure 4 is ran-domly selected in the interval (0, 2.7814), which is enough to merit the above mathematical study.Only the dynamical responses and corresponding phase portrait of model ( 4  are plotted; some symmetric results about P * 2 can be also obtained.As  increases through  * 0 , a periodic solution caused by the phenomenon of Hopf bifurcation occurs; that is, a family of periodic solutions bifurcate from the interior equilibrium P * 1 .Since  2 > 0 and  2 < 0, the Hopf bifurcation is supercritical, the directions of the Hopf bifurcation is  >  * 0 , and these bifurcating periodic solutions from the interior equilibrium P * 1 at  * 0 are stable.Dynamical responses of model ( 4) with  = 3 >  * 0 are plotted in Figure 5. Figures 6 and 7 show a limit cycle corresponding to the periodic solution in Figure 5 in the S-I-R and S-I-E space, respectively.

Conclusion
It is well known that the recovered host individuals are naturally immune to vector disease [1], and its potential economic interest can be commercially exploited.Furthermore, harvest effort is usually influenced by variation of economic interest under market economy.Consequently, it is necessary to discuss the coexistence and interaction mechanism of population within harvested epidemiological ecosystem as well as dynamical effect of harvest effort due to variation of economic interest.
By introducing commercial harvest effort into model proposed in [11], a delayed hybrid mathematical model is established to investigate the dynamical effect of commercial harvesting and incubation time delay on epidemiological economic system, which extends the work done in [11] from a bioeconomic perspective.Positivity and persistence of solutions of model system are discussed in Theorems 3 and 10, respectively.The economic interest of commercial harvesting should be restricted within certain interval that guarantees the existence of interior equilibrium, which can be found in Remark 16.Since the interior equilibrium biologically interprets that susceptible, infective, and recovered host individuals survive as well as harvest on recovered host individuals exist, the bifurcation phenomena around the equilibria can reveal the instability mechanism of model system, which are theoretically relevant to infectious disease control and sustainable yield on recovered host individuals.Consequently, we will mainly concentrate on dynamical behavior and local stability switch around interior equilibrium of model system (4) in this paper.As analyzed in Theorem 11 of this paper, a singularity induced bifurcation occurs which leads to local stability switch in the case of positive economic interest of harvesting.In the perspective of practical viewpoint, a direct damage done by the singularity induced bifurcation to the proposed model is impulse phenomenon, which may lead to outbreak of infectious disease and hamper prosperous harvesting on recovered host individual population resource in the harvested ecosystem, which can be found in Remark 13.With the purpose of maintaining the economic interest at an ideal level, state feedback controllers are designed to stabilize model system around the desirable interior equilibria in the case of zero economic interest and positive economic interest, respectively.The design of the state feedback controller can be found in Theorems 14 and 15 of this paper.The theoretical results and numerical simulations obtained in this paper suggest that incorporating harvest effort on recovered host individuals can not only prevent the stability switch of model system,but also drive model system to stable equilibrium, which will contribute to the persistence and sustainable yield of the harvested ecosystem.
Further attempts are made to understand the dynamical effect of incubation time delay and economic interest on local stability of model system around interior equilibrium.Local stability analysis reveals that incubation delay is responsible for local stability switch of the proposed model, and a family of periodic solutions bifurcate from the interior equilibrium which occurs as incubation delay increases through a critical threshold, which can be found in Theorem 17.The direction and stability of Hopf bifurcation are also discussed in Theorem 19 of this paper, which reveals that Hopf bifurcation is supercritical, the directions of Hopf bifurcation is  >  * 0 , and these bifurcating periodic solutions from the interior equilibrium are stable.The model proposed in [11] does not discuss the harvest effort on economic population.For the model proposed in [11], the threshold value of incubation delay where Hopf bifurcation occurs in [11] is τ0 = 2.0842.However, the harvest effort on recovered host individuals is considered in this paper.As calculated in Section 5.2 of this paper, the Hopf bifurcation occurs at  * 0 = 2.7814 in the case of positive economic interest.It is obvious that  * 0 > τ0 , which implies that the harvesting has a stabilizing impact on the dynamical behavior of population dynamics; cyclic behavior caused by incubation delay can be deferred by introduction of commercial harvesting effort.
It should be noted that some hybrid dynamical models are proposed in [23][24][25][26][27][28], which are utilized to discuss the interaction mechanism of harvested ecosystem from an economic perspective in recent years.Complex dynamical behavior and stability analysis in prey-predator ecosystems with stagestructured population and gestation delay are considered.However, as far as knowledge goes, nobody has explicitly proposed a mathematical model to discuss the dynamic effect of commercial harvest on epidemiological system under the market economy environment.The main objective of this paper is to investigate the transmission mechanism of infectious disease and dynamical effect of commercial harvest on population dynamics, especially the complex dynamical behavior and stability switch due to variation of incubation and commercial harvest economic interest, which makes the work studied in this paper has some new and positive features.

Figure 6 :
Figure 6: A limit cycle corresponding to the periodic solution in Figure 5 in the S-I-R space.

Figure 7 :
Figure 7: Limit cycle corresponding to the periodic solution in Figure 5 in the S-I-E space.

Table 1 :
Signs of real parts of eigen values of model (4) without time delay around interior equilibrium  * .
System with Time Delay.By analyzing corresponding characteristic equation of model system with time delay, local stability analysis around the interior equilibrium due to variation of time delay is discussed.Conditions for existence of Hopf bifurcation are studied.Furthermore, directions of Hopf bifurcation and stability of periodic solutions are investigated.4.2.1.Local Stability and Hopf Bifurcation.As analyzed in the above subsection, in the case of time delay and positive economic interest of harvesting 0 < V * < Ṽ where Ṽ is defined in (55), there are two interior equilibria P * 1 and P * 2 for model system (4) with respect to the positive economic interest V * .

Table 2 :
[49]es of parameters for numerical simulation.By using similar arguments in[49], some properties of bifurcating periodic solutions of model (4) in the center manifold at the critical values are discussed in this paper.Based on the analysis in Section 4.2.2 of this paper, the following theorem can be concluded. 2 determines directions of Hopf bifurcation: if  2 > 0 ( 2 < 0), then Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for  >   ( <   ).