Hopf bifurcation analysis of a delayed ecoepidemiological model with nonlinear incidence rate and Holling type II functional response is investigated. By analyzing the corresponding characteristic equations, the conditions for the stability and existence of Hopf bifurcation for the system are obtained. In addition, a hybrid control strategy is proposed to postpone the onset of an inherent bifurcation of the system. By utilizing normal form method and center manifold theorem, the explicit formulas that determine the direction of Hopf bifurcation and the stability of bifurcating period solutions of the controlled system are derived. Finally, some numerical simulation examples confirm that the hybrid controller is efficient in controlling Hopf bifurcation.
National Natural Science Foundation of China1147211611472115Postgraduate Research & Practice Innovation Program of Jiangsu ProvinceKYCX17_1784Qinglan Project of JiangsuResearch Grant Council of Hong Kong112154151. Introduction
In the natural world, transmissible diseases in ecological environment cannot be ignored. Since the pioneering study of Anderson and May [1], great and interesting predator-prey models with disease were discussed by researchers recently [2–11]. In [2], Liu and Wang considered a predator-prey model with disease in the prey; the Bogdanov-Takens bifurcation and the Hopf bifurcation were analyzed. However, referring to diseases that are transmissible in different populations, Guo et al. [12] studied an ecoepidemiological model with disease spreading within the predator population as follows:(1)x˙t=xtr-a11xt-a12xtSt1+mxt,S˙t=a21xtSt1+mxt-r1St-βStIt,I˙t=βStIt-r2It,where x(t), S(t), and I(t) denote the densities of the prey, the susceptible predator, and the infected predator population at time t, respectively. The parameters r, a11, a12, a21, m, r1, r2, and β in model (1) are all positive constants in which r is the intrinsic growth rate of prey and r/a11 represents the carrying capacity of prey; only the susceptible predators have the ability to capture the prey with capturing rate a12; a21/a12 is the conversion rate of the susceptible predators; m is the half-capturing saturation constant, r1 is the natural death rate of the susceptible predator, r2 is the natural and disease-related mortality rate of the infected predator, and β>0 is called the disease transmission coefficient. Guo et al. discussed the sufficient conditions for the Hopf bifurcation analysis of model (1).
The term βSI in model (1) is called the bilinear incidence rate. In the ecoepidemiological model, it is generally assumed that the average perinfected individual is effectively connected to the other members of the βN population at the same time (N represents the total scale of population), but the activity ability about them at the same time is always limited. Therefore, as the scale of population increases infinitely, the contact rate does not increase, but it gradually tends to a saturated state. This saturation contact rate is also known as the nonlinear incidence rate. Capasso and Serio [13] considered the cholera epidemic spread in Bari in 1973. They introduced a nonlinear incidence rate βIS/(1+αI). This incidence rate seems more reasonable than the bilinear incidence rate βSI, because the nonlinear incidence rate is faster than the linear growth. For instance, before the outbreak of an infectious disease, there would be many contacts with the infected individuals. It includes the behavioral change and crowding effect of the infective individuals and prevents the unboundedness of the contact rate by choosing suitable parameters.
Motivated by the works of Guo et al. [12] and Capasso and Serio [13] and based on the influence about the time delay, in this paper, we consider the following ecoepidemiological model with nonlinear incidence rate, time delay, and Holling type II functional response:(2)x˙t=xtr-a11xt-a12xtSt1+mxt,S˙t=a21xt-τSt-τ1+mxt-τ-r1St-βStIt1+αIt,I˙t=βStIt1+αIt-r2It,where the parameters a11, a12, a21, r, r1, r2, and β are defined in system (1). In this model, only the susceptible predators have ability to capture prey with Holling type II functional response x/1+mx and τ>0 is the time delay due to the gestation of the susceptible predator.
The initial conditions for system (2) take the form of(3)x0>0,S0>0,I0>0.According to the fundamental theory of functional differential equations [14], system (2) has a unique solution (x(t),S(t),I(t)) that satisfies the initial conditions (3). It is easy to show that all solutions of system (2) with initial conditions (3) are defined on [0,+∞) and remain positive for all t≥0.
In recent years, bifurcation control has been extensively concerned by researchers from various disciplines. The aim of bifurcation control is to design a controller to modify the bifurcation properties of a given nonlinear system, thereby achieving some desirable dynamical behaviors. From the control theory point of view, many effective control methods have been proposed, such as the state feedback control [15, 16] and hybrid control strategy [17–24]. Especially, the hybrid control has also been widely used recently. Cheng and Cao [20] considered Hopf bifurcation control for a complex network model with time delays, and they presented a hybrid control strategy to control the model. Kiani et al. [21] used the hybrid control method for a three-pole active magnetic bearing (AMB), and the method showed that the power usage decreased in the hybrid control method comparing to a simple linear control. Peng et al. [24] studied the Hopf bifurcation control for a Lotka-Volterra predator-prey model with two delays by using a hybrid control strategy.
From the viewpoint of an ecological model, the corresponding complex bifurcation behavior means that the system changes from a stable state to an unstable one. It even causes the system to explode, which may be harmful to the ecological balance. Based on this point, a hybrid control strategy by combining the state feedback control and perturbation parameter is used in order to postpone the onset of an inherent bifurcation and enlarge the stable range in model (2).
The rest of this paper is organized as follows. In Section 2, the local stability of the positive equilibrium and the existence of Hopf bifurcation for system (2) are discussed. In Section 3, we propose a hybrid control strategy in which the state feedback and parameter perturbation are combined into system (2) and it is used to control the Hopf bifurcation. The formulas for determining the direction of Hopf bifurcation and the stability of bifurcating period solutions of the controlled system are derived. In Section 4, some numerical simulation examples are carried out to illustrate the validity of the main results. A brief conclusion is given in the last section to conclude this work.
2. A Delayed Ecoepidemiological Model without Control
It is easy to see that system (2) has a unique positive equilibrium E∗(x∗,S∗,I∗), where (4)x∗=-a21rα-rr1αm-rmβ+a11r1α+a11β+Δ12-a11a21α+a11r1mα+a11mβ,S∗=-r21+mx∗x∗a21α-r1mα-mβ-r1α-β,I∗=βS∗-r2r2α,with (5)Δ1=a21rα-rr1mα-rmβ+a11rα+a11β2-4-a11a21α+a11r1mα+a11mβa12r2-rβ-rr1αif the following condition holds: (6)H1:a11r1mα+a11mβ>a11a21α,a12r2<rβ+rr1α.
In this part, we shall study the local stability of linearized system at the positive equilibrium and the existence of Hopf bifurcations for system (2).
Consider the linearized system of system (2) at the positive equilibrium E∗(x∗,S∗,I∗),(7)x˙t=a~11xt+a~12St,S˙t=a~22St+a~23It+b~21xt-τ+b~22St-τ,I˙t=a~32St+a~33It,where (8)a~11=r-2a11x∗-a12S∗1+mx∗2,a~12=-a12x∗1+mx∗,a~22=-r1-βI∗1+αI∗,a~23=-βS∗1+αI∗2,a~32=βI∗1+αI∗,a~33=βS∗1+αI∗2-r2,b~21=a21S∗1+mx∗2,b~22=a21x∗1+mx∗.The characteristic equation of the linearized system (7) is(9)λ3+A2λ2+A1λ+A0+B2λ2+B1λ+B0e-λτ=0,where (10)A0=a~11a~23a~32-a~11a~22a~33,A1=a~11a~22+a~22a~33+a~11a~33-a~23a~32,A2=-a~11+a~22+a~33,B0=a~12a~33b~21-a~11a~33b~22,B1=a~11b~22+a~33b~22-a~12b~21,B2=-b~22.
For τ=0, (9) reduces to(11)λ3+A2+B2λ2+A1+B1λ+A0+B0=0.
It is not difficult to verify that A0+B0>0,A2+B2>0. According to the Routh-Hurwitz criteria, the necessary and sufficient condition for all roots of (11) to have a negative real part is given in the following form: (12)H2:A2+B2A1+B1-A0+B0>0.Namely, the equilibrium E∗(x∗,S∗,I∗) is locally asymptotically stable when the condition H2 is satisfied.
For τ≠0, substituting λ=iω(ω>0) into (9) and separating real and imaginary parts, we obtain(13)-ω3+A1ω=-B2ω2+B0sinωτ-B1ωcosωτ,-A2ω2+A0=B2ω2-B0cosωτ-B1ωsinωτ.Squaring and adding the two equations of (13), it follows that(14)-ω3+A1ω=-B2ω2+B0sinωτ-B1ωcosωτ,-A2ω2+A0=B2ω2-B0cosωτ-B1ωsinωτ,where e20=A02-B02, e21=A12-B12-2A0A2+2B0B2, and e22=A22-2A1-B22.
Let z=ω2. Equation (14) can be written as (15)z3+e22z2+e21z+e20=0.Denote(16)h1z=z3+e22z2+e21z+e20.Since h1(0)=e20, limz→+∞h1(z)=+∞, and from (16), we have(17)h1′z=3z2+2e22z+e21.After discussion about the roots of (17) that are similar to those in [25], we have the following lemma.
Lemma 1.
For the polynomial equation (15), we have the following results.
(1) If (H21) e20≥0 and Δ=e222-3e21≤0 holds, then (15) has no positive root.
(2) If (H22) e20≥0, Δ=e222-3e21>0, z∗=-e21+Δ/3>0, h1(z∗)≤0, or (H23) e20<0 holds, then (15) has positive roots.
Suppose that (15) has positive roots. Without loss of generality, we assume that it has three positive roots, denoted by z1, z2, and z3. Then (14) has three positive roots ωk=k,k=1,2,3. The corresponding critical value of time delay τk(j) is(18)τkj=1ωkarccos-B1+A2B2ωk4+A1B1-A2B0-A0B2ωk2+A0B0-B22ωk4+B12-2B0B2ωk2+B02+2πjωk,where j=0,1,2….
Thus ±ωk is a pair of purely imaginary roots of (9) with τ=τk(j), and let τ0=mink∈1,2,3τk(0),ω0=ωk0.
According to the Hopf Bifurcation Theorem [26], we need to verify the transversality condition. Differentiating the two sides of (9) with respect to τ, we obtain(19)dλdτ-1=-3λ2+2A2λ+A1λλ3+A2λ2+A1λ+A0+2B2λ+B1λB2λ2+B1λ+B0-τλ,Then(20)signdReλdτλ=iω0=signRedλdτ-1λ=iω0=sign3ω04+2A22-2A1ω02+A12-2A0A2ω06+-2A1+A22ω04+A12-2A0A2ω02+A02-2B22ω02+B12-2B0B2B22ω04+B12-2B0B2ω02+B02.We derive from (13) that (21)ω06+-2A1+A22ω04+A12-2A0A2+2B0B2+A02=B22ω04+B12-2B0B2ω02+B02.Then, it follows that (22)dReλdτλ=iω0=sign3ω04+2A22-2A1-B22ω02+A12-2A0A2-B12+2B0B2B22ω04+B12-2B0B2ω02+B02=signh1′z0B22ω04+B12-2B0B2ω02+B02.Therefore, d(Reλ)/dτλ=iω0≠0 if the following condition holds: (23)H24:h1′ω02≠0.According to the analysis above, we have the following results.
Theorem 2.
For system (2),
(1) If (H21) holds, then the positive equilibrium E∗(x∗,S∗,I∗) is asymptotically stable for all τ≥0.
(2) If (H22) or (H23) and (H24) hold, then the positive equilibrium E∗(x∗,S∗,I∗) is asymptotically stable for all τ∈[0,τ0) and unstable for τ>τ0. Furthermore, system (2) undergoes a Hopf bifurcation at the positive equilibrium E∗(x∗,S∗,I∗) when τ=τ0.
3. A Delayed Ecoepidemiological Model with Hybrid Control
In this part, a hybrid control strategy is proposed, in which the state feedback and parameter perturbation are combined in an effort to postpone the occurrence of Hopf bifurcation in system (2). Here, a controlled model as follows is considered:(24)x˙t=pxtr-a11xt-a12xtSt1+mxt+qxt,S˙t=pa21xt-τSt-τ1+mxt-τ-r1St-βStIt1+αIt+qSt,I˙t=pβStIt1+αIt-r2It,where p>0 and q∈R are control parameters. The parameters a11, a12, a21, m, r, r1, r2, β, and τ are defined in system (2), qx(t) and qS(t) affect the densities of prey and susceptible predator at time t, respectively, and q>0 denotes increase in the quantity, while q<0 otherwise.
Similar to the discussion in Section 2, model (24) has a unique positive equilibrium E∘(x∘,S∘,I∘), where (25)x∘=-C2+Δ22C1,S∘=-pr21+mx∘x∘αpa21-αpr1m-pβm+αqm-αpr1-pβ+αq,I∘=βS∘-r2r2α,with (26)C1=-a11pαpa21-αpr1m-pβm+αqm,C2=αp2ra21-αp2rr1m-p2βrm+pqrαm+αp2r1a11+p2a11β-αpqa11+αpqa21-αpqr1m-pqβm+αq2m,Δ2=C22-4C1-αp2rr1+αpqr-p2βr-αpqr1-pqβ+αq2+p2a12r2,if the following condition holds: (27)H3:-a11pαpa21-αpr1m-pβm+αqm>0,-αp2rr1+αpqr-p2βr-αpqr1-pqβ+αq2+p2a12r2<0.Let x¯(t)=x(t)-x∘, S¯(t)=S(t)-S∘, and I¯(t)=I(t)-I∘ and still denote x¯(t),S¯(t),I¯(t), respectively. Using Taylor expansion to expand system (24) at the positive equilibrium E∘(x∘,S∘,I∘), we have(28)x˙t=a11′xt+a12′St+f1′,S˙t=a22′St+a23′It+b21′xt-τ+b22′St-τ+f2′,I˙t=a32′St+a33′It+f3′,where (29)a11′=pr-2pa11x∘-a12pS∘1+mx∘2+q,a12′=-a12px∘1+mx∘,a22′=-pr1-pβI∘1+αI∘+q,a23′=-pβS∘1+αI∘,a32′=pβI∘1+αI∘,a33′=pβS∘1+αI∘2-r2p,b21′=pa21S∘1+mx∘∘,b22′=pa21x∘1+mx∘,f1′=a13′x2+a14′xS+a15′x2S+a16′x3,f2′=a24′x2t-τ+a25′xt-τSt-τ+a26′x2t-τSt-τ+a27′x3t-τ+a28′SI+a29′I2+a30′I2S+a31′I3,f3′=a34′SI+a35′I2+a36′I2S+a37′I3,with (30)a13′=-a11p+pa12mS∘1+mx∘3,a14′=-pa121+mx∘2,a15′=pa12m1+mx∘3,a16′=-pa12m2S∘1+mx∘4,a24′=-pa21mS∘1+mx∘3,a25′=pa211+mx∘2,a26′=-pa21m1+mx∘3,a27′=pa21m2S∘1+mx∘4,a28′=-pβ1+αI∘2,a29′=pαβS∘1+αI∘3,a30′=pαβ1+αI∘3,a31′=-pα2βS∘1+αI∘4,a34′=pβ1+αI∘2,a35′=-pαβS∘1+αI∘3,a36′=-pαβ1+αI∘3,a37′=pα2βS∘1+αI∘4.Then we obtain the linearized system of system (24) as follows:(31)x˙t=a11′xt+a12′St,S˙t=a22′St+a23′It+b21′xt-τ+b22′St-τ,I˙t=a32′St+a33′It.Therefore, the corresponding characteristic equation of system (31) is given by(32)λ3+A2′λ2+A1′λ+A0′+B2′λ2+B1′λ+B0′e-λτ′=0,where (33)A0′=a11′a23′a32′-a11′a22′a33′,A1′=a11′a22′+a22′a33′+a11′a33′-a23′a32′,A2′=-a11′+a22′+a33′,B0′=a12′a33′b21′-a11′a33′b22′,B1′=a11′b22′+a33′b22′-a12′b21′,B2′=-b22′.Obviously, the characteristic equation of system (24) is similar to (9). Therefore, the analysis method is very similar to Section 2; we will omit the local stability and Hopf bifurcation analysis of system (24). We obtain the corresponding critical value of time delay (τk′)(j) as(34)τk′j=1ωk′arccos-B1′+A2′B2′ωk′4+A1′B1′-A2′B0′-A0′B2′ωk′2+A0′B0′-B2′2ωk′4+B1′2-2B0′B2′ωk′2+B0′2+2πjωk′,where k=1,2,3;j=0,1,2,…,ωk′ is a positive root of(35)ω′6+e32ω′4+e31ω′2+e30=0,with (36)e30=A0′2-B0′2,e31=A1′2-B1′2-2A0′A2′+2B0′B2′,e32=A2′2-2A1′-B2′2.Thus ±ωk′ is a pair of purely imaginary roots of (32) with τ′=(τk′)(j), and let τ0′=mink∈1,2,3(τk′)(0),ω0′=ωk0′.
In the following, we will investigate the direction of Hopf bifurcation and the stability of bifurcating periodic solutions of the controlled system (24) at τ0′. The theoretical approach we will apply is based on the normal form theory and center manifold theorem [26].
Let τ=τ0′+μ,μ∈R, t=sτ, x(sτ)=x^(s), S(sτ)=S^(s), and I(sτ)=I^(s). Then μ=0 is the Hopf bifurcation value of the controlled system (24). Denote x=x^, S=S^, I=I^, and t=s, then system (24) can be written as a functional differential equation (FDE) in C=C([-1,0],R3):(37)u′t=Lμut+Fμ,ut,where u(t)=(x(t),S(t),I(t))T∈C, ut(θ)=u(t+θ)=(x(t+θ),S(t+θ),I(t+θ))T∈C, and Lμ:C→R3, F:R×C→R3 are given by(38)Lμϕ=τ0′+μA′ϕ0+B′ϕ-1,Fμ,ϕ=τ0′+μF1,F2,F3T,where (39)ϕθ=ϕ1θ,ϕ2θ,ϕ3θT∈C,A′=a11′a12′00a22′a23′0a32′a33′,B′=000b21′b22′0000,F1=a13′ϕ120+a14′ϕ10ϕ20+a15′ϕ120ϕ20+a16′ϕ130,F2=a24′ϕ12-1+a25′ϕ1-1ϕ2-1+a26′ϕ12-1ϕ2-1+a27′ϕ13-1+a28′ϕ20ϕ30+a29′ϕ320+a30′ϕ20ϕ320+a31′ϕ330,F3=a34′ϕ20ϕ30+a35′ϕ320+a36′ϕ20ϕ320+a37′ϕ330.
Hence, by the Riesz representation theorem, there exists a 3×3 matrix function η(θ,μ) of bounded variation for θ∈[-1,0], such that(40)Lμϕ=∫-10dηθ,μϕθ,for ϕ∈C.In fact, we can choose(41)ηθ,μ=τ0′+μA′δθ-τ0′+μB′δθ+1,where δ(θ) is the Dirac function.
For ϕ∈C([-1,0],R3), define(42)Aμϕ=dϕθdθ,-1≤θ<0,∫-10dηθ,μϕθ,θ=0,Rμϕ=0,-1≤θ<0,Fμ,ϕ,θ=0.Then (37) can be transformed into the following operator equation:(43)ut′=Aμut+Rμut.The adjoint operator A∗ of A(0) is defined by(44)A∗φs=-dφsds,s∈0,1,∫-10dηTt,0φ-t,s=0.For ϕ∈C([-1,0],R3) and φ∈C([-1,0],(R3)∗), define the bilinear form:(45)φs,ϕs=φ¯0ϕ0-∫-10∫ξ=0θφ¯ξ-θdηθϕξdξ,where η(θ)=η(θ,0),A=A(0) and A∗ are adjoint operators.
Referring to the previous discussion, we know that ±iω0′τ0′ are the eigenvalues of A(0); thus they are also the eigenvalues of A∗. Suppose that q(θ)=(1,q1,q2)Teiω0′τ0′θ is the eigenvector of A(0) corresponding to iω0′τ0′ and q∗(s)=D(1,q1∗,q2∗)eiω0′τ0′s is the eigenvector of A∗ corresponding to -iω0′τ0′. By direct computation, we obtain(46)q1=iω0′-a11′a12′,q2=a32′iω0′-a11′a12′iω0′-a33′,q1∗=-iω′-a11′b21′eiω0′τ0′,q2∗=a23′iω0′+a11′b21′eiω0′τ0′iω0′+a33′,q∗(s),q(θ)=1, and q∗(s),q¯(θ)=0, where(47)D¯=11+q1q¯1∗+q2q¯2∗+e-iω0′τ0′τ0′q1∗b21′+b22′q1.
Next, we can obtain the coefficients used in determining the direction of Hopf bifurcation and the stability of the bifurcation periodic solutions by the algorithms given in [26].(48)g20=2D¯τ0′a13′+a14′q1+q¯1∗a24′e-2iω0′τ0′+a25′q1e-2iω0′τ0′+a28′q1q2+a29′q22+q¯2∗a34′q1q2+a35′q22,g11=D¯τ0′2a13′+a14′q1+q¯1+q¯1∗2a24′+a25′q1+q¯1+a28′q1q¯2+q¯1q2+2a29′q2q¯2+q¯2∗a34′q1q¯2+q¯1q2+2a35′q2q¯2,g02=2D¯τ0′a13′+a14′q¯1+q¯1∗a24′e2iω0′τ0′+a25′q¯1e2iω0′τ0′+a28′q¯1q¯2+a29′q¯22+q¯2∗a34′q¯1q¯2+a35′q¯22,g21=2D¯τ0′a13′2W1110+W2010+a14′W1120+12W2020+12q¯1W2010+q1W1110+a15′2q1+q¯1+3a16′+q¯1∗a24′2W111-1e-iω0′τ0′+W201-1e-iω0′τ0′+a25′W112-1e-iω0′τ0′+12W202-1eiω0′τ0′+12W201-1q¯1eiω0′τ0′+W111-1q1e-iω0′τ0′+a26′2q1e-iω0′τ0′+q¯1e-iω0′τ0′+3a27′e-iω0′τ0′+a28′q1W1130+12q¯1W2030+12q¯2W2020+q2W1120+a29′2q2W1130+q¯2W2030+a30′2q1q2q¯2+q¯1q22+3a31′q22q¯2+q¯2∗a34′q1W1130+12q¯1W2030+12q¯2W2020+q2W1120+a35′2q2W1130+q¯2W2030+a36′2q1q2q¯2+q¯1q22+2a37′q22q¯2.However,(49)W20θ=ig20ω0′τ0′q0eiω0′τ0′θ+ig¯023ω0′τ0′q¯0e-iω10′τ0′θ+E1e2iω10′τ10′θ,W11θ=-ig11ω0′τ0′q0eiω0′τ0θ+ig¯11ω10′τ10′q¯0e-iω0′τ0′θ+E2,where E1=(E1(1),E1(2),E1(3))T∈R3 and E2=(E2(1),E2(2),E2(3))T∈R3 are also constant vectors and they can be determined, respectively, by (50)2iω0′-a11′-a12′0-b21′e-iω0′τ0′2iω10′-a22′-b22′e-iω0′τ0′-a23′0-a32′2iω10′-a33′E1=2H1H2H3,-a11′-a12′0-b21′-a22′-b22′-a23′0-a32′-a33′E2=P1P2P3,with (51)H1=a13′+a14′q1,H2=a24′e-2iω0′τ0′+a25′q1e-2iω0′τ0′+a28′q1q2+a29′q22,H3=a34′q1q2+a35′q22,P1=2a13′+a14′q1+q¯1,P2=2a24′+a25′q1+q¯1+a28′q1q¯2+q¯1q2+2a29′q2q¯2,P3=a34′q1q¯2+q¯1q2+2a35′q2q¯2.Therefore, we can determine g21 and derive the expressions(52)c10=i2ω0′τ0′g20g11-2g112-g0223+g212,μ2=-Rec10Reλ′τ0′,β2=2Rec10,T2=-Imc10+μ2Imλ′τ0ω0′τ0′,which describe the properties of bifurcation period solutions at τ=τ0′ on the center manifold. From the discussion above, we have the following result.
Theorem 3.
For system (24), the direction of Hopf bifurcation is determined by the sign of μ2: if μ2>0(μ2<0), then the Hopf bifurcation is supercritical (subcritical). The stability of the bifurcating periodic solutions is determined by the sign of β2: if β2<0(β2>0), the bifurcating periodic solutions are stable (unstable). The period of the bifurcating periodic solutions is determined by the sign of T2: if T2>0(T2<0), the bifurcating periodic solutions increase (decrease).
4. Numerical Examples
In this section, we present some numerical examples by using Matlab to verify the analytical predictions obtained in the previous sections. The hybrid control strategy to gain control of the Hopf bifurcation in model (2) is applied.
Let a11=0.125, a12=1.8, a21=1.35, r=0.55, r1=0.18, r2=0.25, m=0.01, α=0.0005, and β=0.96. Then, we have the following particular example of system (2):(53)x˙t=xt0.55-0.125xt-1.8xtSt1+0.01xt,S˙t=1.35xt-τSt-τ1+0.01xt-τ-0.18St-0.96StIt1+0.0005It,I˙t=0.96StIt1+0.0005It-0.25It.
It is easy to show that if (H1) holds, system (53) has a unique coexistence equilibrium E∗(0.6723,0.2605,0.7519). For τ=0, (H2) is satisfied and then the equilibrium is locally asymptotically stable. For τ≠0, we obtain ω0=0.7176, τ0=0.1718, and h1′=0.3476≠0; that is, the transversal condition is satisfied. From Theorem 2, the coexistence equilibrium E∗(0.6723,0.2605,0.7519) is asymptotically stable for τ∈(0,τ0]. For τ=0.1<τ0, which can be shown in Figure 1, the positive equilibrium E∗(0.6723,0.2605,0.7519) is unstable for τ>τ0. For τ=0.3>τ0, this property can be illustrated in Figure 2.
Behavior and phase portrait of the uncontrolled system (2) with τ=0.1<τ0. The positive equilibrium E∗ is asymptotically stable.
Behavior and phase portrait of the uncontrolled system (2) with τ=0.3>τ0. Hopf bifurcation occurs from the positive equilibrium E∗.
Next, we choose appropriate values of p,q to control system (2). Let us consider the following system with hybrid control strategy:(54)x˙t=0.3xt0.55-0.125xt-1.8xtSt1+0.01xt+0.05xt,S˙t=0.31.35xt-τSt-τ1+0.01xt-τ-0.18St-0.96StIt1+0.0005It+0.05St,I˙t=0.30.96StIt1+0.0005It-0.25It.
It is not difficult to verify that if (H3) holds, we obtain the positive equilibrium E∘(2.0434,0.2606,1.4080). From the analysis in Section 3, we obtain ω0′=0.3192 and τ0′=0.9075. By choosing τ′=0.3 and τ′=1.3, the dynamical behavior of this controlled model (54) is illustrated in Figures 3 and 4. That is, for τ′=0.3<τ0′, the positive equilibrium E∘ is asymptotically stable. However, when the time delay τ′ passes through the critical value τ0′, the positive equilibrium E∘ will lose its stability, a Hopf bifurcation occurs, and a family of periodic solutions bifurcate from the positive equilibrium E∘.
Behavior and phase portrait of the controlled system (24) with τ′=0.3<τ0′. The positive equilibrium E∘ is asymptotically stable.
Behavior and phase portrait of the controlled system (24) with τ′=1.3>τ0′. Hopf bifurcation occurs from the positive equilibrium E∘.
Comparing Figures 3 and 4 with Figures 1 and 2, we can easily make the Hopf bifurcation of the uncontrolled system (2) disappear. It is shown that the onset of Hopf bifurcation is delayed when the hybrid controller has been incorporated into the model and the critical value of delay increases from τ0=0.1718 to τ0′=0.9075. After computation of (52), we obtain c1(0)=-8.5035+6.0113i, T2=68.0768, μ2=368.1169, and β2=-17.0070. From Theorem 3, the Hopf bifurcation is supercritical, the bifurcation period solutions are stable, and the bifurcating periodic solutions increase.
These numerical simulation results illustrate excellent validations of the new theoretical analysis presented in this paper. Because the bifurcation periodic solutions are stable, the species in model (24) imply coexistence in an oscillatory mode from the viewpoint of biology.
5. Conclusions
In this paper, we have incorporated nonlinear incidence rate and time delay into an ecoepidemiological model. By analyzing the associated characteristic equation, its local stability and the existence of Hopf bifurcation with respect to time delay are established. To postpone the onset of the Hopf bifurcation, we use a hybrid controller in system (2). It has been shown that the critical value of the delay increases from τ0=0.1718 to τ0′=0.9075. By the normal form theory and the center manifold theorem, we analyze the stability and direction of the bifurcating periodic solutions. The hybrid control strategy is closely related to the continuous survival of the population. The unstable equilibrium point of a uncontrolled system becomes an asymptotically stable equilibrium point in the controlled system after using hybrid control by combining the state feedback control and perturbation parameter. That is, the number of predators and preys eventually implies stability and coexistence. The numerical results are in accord with theoretical analysis. It has certain ecological significance and provides a theoretical basis for the continuous survival of the population in nature.
In addition, the model with the disease spreading in predator population is investigated in this paper. If diseasees spreading in prey and predator population coexist together, the influence will be analyzed in a future work.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Nature Science foundation of China (11472116, 11472115), Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX17_1784), Qinglan Project of Jiangsu, and the Research Grant Council of Hong Kong (11215415).
AndersonR. M.MayR. M.Regulation and Stability of Host-Parasite Population Interactions: I. Regulatory Processes197847121910.2307/3933LiuX.WangC.Bifurcation of a predator-prey model with disease in the prey201062484185010.1007/s11071-010-9766-7MR2745944Zbl1215.370532-s2.0-78049453226ZhouX.CuiJ.Stability and Hopf bifurcation analysis of an eco-epidemiological model with delay201034791654168010.1016/j.jfranklin.2010.08.001MR2729478Zbl1208.341362-s2.0-77958089819ShiX.CuiJ.ZhouX.Stability and Hopf bifurcation analysis of an eco-epidemic model with a stage structure20117441088110610.1016/j.na.2010.09.038MR2746791HuG.-P.LiX.-L.Stability and Hopf bifurcation for a delayed predator-prey model with disease in the prey201245322923710.1016/j.chaos.2011.11.011MR2881651GaoX.PanQ.HeM.KangY.A predator-prey model with diseases in both prey and predator2013392235898590610.1016/j.physa.2013.07.077MR31238882-s2.0-84884673143CaiY.WangW.Stability and Hopf bifurcation of the stationary solutions to an epidemic model with cross-diffusion20157081906192010.1016/j.camwa.2015.08.003MR3400496WangX.PengM.LiuX.Stability and Hopf bifurcation analysis of a ratio-dependent predator-prey model with two time delays and Holling type III functional response201526849650810.1016/j.amc.2015.06.108MR3399436ZhangX.HuangY.WengP.Stability and bifurcation of a predator-prey model with disease in the prey and temporal-spatial nonlocal effect201629046748610.1016/j.amc.2016.06.050MR35234432-s2.0-84978891152SaifuddinM.BiswasS.SamantaS.SarkarS.ChattopadhyayJ.Complex dynamics of an eco-epidemiological model with different competition coefficients and weak Allee in the predator20169127028510.1016/j.chaos.2016.06.009MR3551708HuZ.TengZ.ZhangT.ZhouQ.ChenX.Globally asymptotically stable analysis in a discrete time eco-epidemiological system201799203110.1016/j.chaos.2017.03.042MR3650006GuoZ. K.LiW. L.ChengL. H.LiZ. Z.Eco-epidemiological model with epidemic and response function in the predator2009453117121MR2560354CapassoV.SerioG.A generalization of the Kermack-McKendrick deterministic epidemic model1978421-2436110.1016/0025-5564(78)90006-8MR529097HaleJ. K.1977New York, NY, USASpringerMR0508721Zbl0352.34001JiangG.LuQ.QianL.Complex dynamics of a Holling type II prey-predator system with state feedback control200731244846110.1016/j.chaos.2005.09.077MR2259769LiM.ChenB.WangG.LiZ.LiuW.Hopf bifurcation and nonlinear state feedback control for a modified Lotka-Volterra differential algebraic predator-prey systemProceedings of the 5th International Conference on Intelligent Control and Information Processing, ICICIP 2014August 201423323810.1109/ICICIP.2014.70103452-s2.0-84925330090LuoX. S.ChenG. R.WangB. H.FangJ. Q.Hybrid control of period-doubling bifurcation and chaos in discrete nonlinear dynamical systems200318477578310.1016/s0960-0779(03)00028-62-s2.0-0038239703Zbl1073.37512LiuZ.ChungK. W.Hybrid control of bifurcation in continuous nonlinear dynamical systems200515123895390310.1142/S0218127405014374MR2205624Zbl1093.935452-s2.0-33644691604ChengZ.CaoJ.Hopf bifurcation control for delayed complex networks2007344684685710.1016/j.jfranklin.2006.10.007MR2344914Zbl1126.370332-s2.0-34547420752ChengZ. S.CaoJ. D.Hybrid control of Hopf bifurcation in complex networks with delays201413116417010.1016/j.neucom.2013.10.0282-s2.0-84894098406KianiM.SalariehH.AlastyA.DarbandiS. M.Hybrid control of a three-pole active magnetic bearing20163928412-s2.0-8497992998710.1016/j.mechatronics.2016.07.004CaiP.YuanZ. Z.Hopf bifurcation and chaos control in a new chaotic system via hybrid control strategy2017551647010.1016/j.cjph.2016.12.002MR36007372-s2.0-85010669601AlfiA.KalatA. A.FarrokhnejadF.Hybrid control strategy applied to chaos synchronization: new control design and stability analysis201721810.1007/s40435-017-0343-5PengM.ZhangZ.WangX.Hybrid control of Hopf bifurcation in a Lotka-Volterra predator-prey model with two delays201738710.1186/s13662-017-1434-5MR3738699SongY.WeiJ.Bifurcation analysis for Chen's system with delayed feedback and its application to control of chaos2004221759110.1016/j.chaos.2003.12.075MR2057549Zbl1112.373032-s2.0-1842581482HassardB. D.KazarinoffN. D.WanY.-H.1981Cambridge, UKCambridge University PressMR603442