Processing math: 100%
Research article Special Issues

Non-smooth dynamics of a SIR model with nonlinear state-dependent impulsive control


  • Received: 17 July 2023 Revised: 14 September 2023 Accepted: 20 September 2023 Published: 09 October 2023
  • The classic SIR model is often used to evaluate the effectiveness of controlling infectious diseases. Moreover, when adopting strategies such as isolation and vaccination based on changes in the size of susceptible populations and other states, it is necessary to develop a non-smooth SIR infectious disease model. To do this, we first add a non-linear term to the classical SIR model to describe the impact of limited medical resources or treatment capacity on infectious disease transmission, and then involve the state-dependent impulsive feedback control, which is determined by the convex combinations of the size of the susceptible population and its growth rates, into the model. Further, the analytical methods have been developed to address the existence of non-trivial periodic solutions, the existence and stability of a disease-free periodic solution (DFPS) and its bifurcation. Based on the properties of the established Poincaré map, we conclude that DFPS exists, which is stable under certain conditions. In particular, we show that the non-trivial order-1 periodic solutions may exist and a non-trivial order-k (k1) periodic solution in some special cases may not exist. Moreover, the transcritical bifurcations around the DFPS with respect to the parameters p and AT have been investigated by employing the bifurcation theorems of discrete maps.

    Citation: Chenxi Huang, Qianqian Zhang, Sanyi Tang. Non-smooth dynamics of a SIR model with nonlinear state-dependent impulsive control[J]. Mathematical Biosciences and Engineering, 2023, 20(10): 18861-18887. doi: 10.3934/mbe.2023835

    Related Papers:

    [1] Dali Makharadze, Alexander Meskhi, Maria Alessandra Ragusa . Regularity results in grand variable exponent Morrey spaces and applications. Electronic Research Archive, 2025, 33(5): 2800-2814. doi: 10.3934/era.2025123
    [2] Kwok-Pun Ho . Martingale transforms on Banach function spaces. Electronic Research Archive, 2022, 30(6): 2247-2262. doi: 10.3934/era.2022114
    [3] Yaning Li, Mengjun Wang . Well-posedness and blow-up results for a time-space fractional diffusion-wave equation. Electronic Research Archive, 2024, 32(5): 3522-3542. doi: 10.3934/era.2024162
    [4] Peng Gao, Pengyu Chen . Blowup and MLUH stability of time-space fractional reaction-diffusion equations. Electronic Research Archive, 2022, 30(9): 3351-3361. doi: 10.3934/era.2022170
    [5] María Guadalupe Morales, Zuzana Došlá, Francisco J. Mendoza . Riemann-Liouville derivative over the space of integrable distributions. Electronic Research Archive, 2020, 28(2): 567-587. doi: 10.3934/era.2020030
    [6] R. F. Snider . Eigenvalues and eigenvectors for a hermitian gaussian operator: Role of the Schrödinger-Robertson uncertainty relation. Electronic Research Archive, 2023, 31(9): 5541-5558. doi: 10.3934/era.2023281
    [7] Yaning Li, Yuting Yang . The critical exponents for a semilinear fractional pseudo-parabolic equation with nonlinear memory in a bounded domain. Electronic Research Archive, 2023, 31(5): 2555-2567. doi: 10.3934/era.2023129
    [8] Mehmet Ali Özarslan, Ceren Ustaoğlu . Extended incomplete Riemann-Liouville fractional integral operators and related special functions. Electronic Research Archive, 2022, 30(5): 1723-1747. doi: 10.3934/era.2022087
    [9] Harish Bhatt . Second-order time integrators with the Fourier spectral method in application to multidimensional space-fractional FitzHugh-Nagumo model. Electronic Research Archive, 2023, 31(12): 7284-7306. doi: 10.3934/era.2023369
    [10] Humberto Rafeiro, Joel E. Restrepo . Revisiting Taibleson's theorem. Electronic Research Archive, 2022, 30(2): 565-573. doi: 10.3934/era.2022029
  • The classic SIR model is often used to evaluate the effectiveness of controlling infectious diseases. Moreover, when adopting strategies such as isolation and vaccination based on changes in the size of susceptible populations and other states, it is necessary to develop a non-smooth SIR infectious disease model. To do this, we first add a non-linear term to the classical SIR model to describe the impact of limited medical resources or treatment capacity on infectious disease transmission, and then involve the state-dependent impulsive feedback control, which is determined by the convex combinations of the size of the susceptible population and its growth rates, into the model. Further, the analytical methods have been developed to address the existence of non-trivial periodic solutions, the existence and stability of a disease-free periodic solution (DFPS) and its bifurcation. Based on the properties of the established Poincaré map, we conclude that DFPS exists, which is stable under certain conditions. In particular, we show that the non-trivial order-1 periodic solutions may exist and a non-trivial order-k (k1) periodic solution in some special cases may not exist. Moreover, the transcritical bifurcations around the DFPS with respect to the parameters p and AT have been investigated by employing the bifurcation theorems of discrete maps.



    In this paper, we consider the Schrödinger operators

    L=+V(x),xRn,n3,

    where Δ=ni=122xi and V(x) is a nonnegative potential belonging to the reverse Hölder class RHq for some qn2. Assume that f is a nonnegative locally Lq(Rn) integrable function on Rn, then we say that f belongs to RHq (1<q) if there exists a positive constant C such that the reverse Hölder's inequality

    (1|B(x,r)|B(x,r)|f(y)|qdy)1qC|B(x,r)|B(x,r)|f(y)|dy

    holds for x in Rn, where B(x,r) denotes the ball centered at x with radius r< [1]. For example, the nonnegative polynomial VRH, in particular, |x|2RH.

    Let the potential VRHq with qn2, and the critical radius function ρ(x) is defined as

    ρ(x)=supr>0{r:1rn2B(x,r)V(y)dy1},xRn. (1.1)

    We also write ρ(x)=1mV(x),xRn. Clearly, 0<mV(x)< when V0, and mV(x)=1 when V=1. For the harmonic oscillator operator (Hermite operator) H=Δ+|x|2, we have mV(x)(1+|x|).

    Thanks to the heat diffusion semigroup etL for enough good function f, the negative powers Lα2(α>0) related to the Schrödinger operators L can be written as

    Iαf(x)=Lα2f(x)=0etLf(x)tα21dt,0<α<n. (1.2)

    Applying Lemma 3.3 in [2] for enough good function f holds that

    Iαf(x)=RnKα(x,y)f(y)dy,0<α<n,

    and the kernel Kα(x,y) satisfies the following inequality

    Kα(x,y)Ck(1+|xy|(mV(x)+mV(y)))k1|xy|nα. (1.3)

    Moreover, we have Kα(x,y)C|xy|nα,0<α<n.

    Shen [1] obtained Lp estimates of the Schrödinger type operators when the potential VRHq with qn2. For Schrödinger operators L=Δ+V with VRHq for some qn2, Harboure et al. [3] established the necessary and sufficient conditions to ensure that the operators Lα2(α>0) are bounded from weighted strong and weak Lp spaces into suitable weighted BMOL(w) space and Lipschitz spaces when pnα. Bongioanni Harboure and Salinas proved that the fractional integral operator Lα/2 is bounded form Lp,(w) into BMOβL(w) under suitable conditions for weighted w [4]. For more backgrounds and recent progress, we refer to [5,6,7] and references therein.

    Ramseyer, Salinas and Viviani in [8] studied the fractional integral operator and obtained the boundedness from strong and weak Lp() spaces into the suitable Lipschitz spaces under some conditions on p(). In this article, our main interest lies in considering the properties of fractional integrals operator Lα2(α>0), related to L=Δ+V with VRHq for some qn2 in variable exponential spaces.

    We now introduce some basic properties of variable exponent Lebsegue spaces, which are used frequently later on.

    Let p():Ω[1,) be a measurable function. For a measurable function f on Rn, the variable exponent Lebesgue space Lp()(Ω) is defined by

    Lp()(Ω)={f:Ω|f(x)s|p(x)dx<},

    where s is a positive constant. Then Lp()(Ω) is a Banach space equipped with the follow norm

    fLp()(Ω):=inf{s>0:Ω|f(x)s|p(x)dx1}.

    We denote

    p:=essinfxΩp(x) and p+:=esssupxΩp(x).

    Let P(Rn) denote the set of all measurable functions p on Rn that take value in [1,), such that 1<p(Rn)p()p+(Rn)<.

    Assume that p is a real value measurable function p on Rn. We say that p is locally log-Hölder continuous if there exists a constant C such that

    |p(x)p(y)|Clog(e+1/|xy|),x,yRn,

    and we say p is log-Hölder continuous at infinity if there exists a positive constant C such that

    |p(x)p()|Clog(e+|x|),xRn,

    where p():=lim|x|p(x)R.

    The notation Plog(Rn) denotes all measurable functions p in P(Rn), which states p is locally log-Hölder continuous and log-Hölder continuous at infinity. Moreover, we have that p()Plog(Rn), which implies that p()Plog(Rn).

    Definition 1.1. [8] Assume that p() is an exponent function on Rn. We say that a measurable function f belongs to Lp(),(Rn), if there exists a constant C such that for t>0,

    Rntp(x)χ{|f|>t}(x)dxC.

    It is easy to check that Lp(),(Rn) is a quasi-norm space equipped with the following quasi-norm

    fp(),=inf{s>0:supt>0Rn(ts)p(x)χ{|f|>t}(x)dx1}.

    Next, we define LipLα,p() spaces related to the nonnegative potential V.

    Definition 1.2. Let p() be an exponent function with 1<pp+< and 0<α<n. We say that a locally integrable function fLipLα,p()(Rn) if there exist constants C1,C2 such that for every ball BRn,

    1|B|αnχBp()B|f(x)mBf|dxC1, (1.4)

    and for Rρ(x),

    1|B|αnχBp()B|f(x)|dxC2, (1.5)

    where mBf=1|B|Bf. The norm of space LipLα,p()(Rn) is defined as the maximum value of two infimum of constants C1 and C2 in (1.4) and (1.5).

    Remark 1.1. LipLα,p()(Rn)Lα,p()(Rn) is introduced in [8]. In particular, when p()=C for some constant, then LipLα,p()(Rn) is the usual weighted BMO space BMOβL(w), with w=1 and β=αnp [4].

    Remark 1.2. It is easy to see that for some ball B, the inequality (1.5) leads to inequality (1.4) holding, and the average mBf in (1.4) can be replaced by a constant c in following sense

    12fLipLα,p()supBRninfcR1|B|αnχBp()B|f(x)c|dxfLipLα,p().

    In 2013, Ramseyer et al. in [8] studied the Lipschitz-type smoothness of fractional integral operators Iα on variable exponent spaces when p+>αn. Hence, when p+>αn, it will be an interesting problem to see whether or not we can establish the boundedness of fractional integral operators Lα2(α>0) related to Schrödinger operators from Lebesgue spaces Lp() into Lipschitz-type spaces with variable exponents. The main aim of this article is to answer the problem above.

    We now state our results as the following two theorems.

    Theorem 1.3. Let potential VRHq for some qn/2 and p()Plog(Rn). Assume that 1<pp+<n(αδ0)+ where δ0=min{1,2n/q}, then the fractional integral operator Iα defined in (1.2) is bounded from Lp()(Rn) into LipLα,p()(Rn).

    Theorem 1.4. Let the potential VRHq with qn/2 and p()Plog(Rn). Assume that 1<pp+<n(αδ0)+ where δ0=min{1,2n/q}. If there exists a positive number r0 such that p(x)p when |x|>r0, then the fractional integral operator Iα defined in (1.2) is bounded from Lp(),(Rn) into LipLα,p()(Rn).

    To prove Theorem 1.3, we first need to decompose Rn into the union of some disjoint ball B(xk,ρ(xk))(k1) according to the critical radius function ρ(x) defined in (1.1). According to Lemma 2.6, we establish the necessary and sufficient conditions to ensure fLipLα,p()(Rn). In order to prove Theorem 1.3, by applying Corollary 1 and Remark 1.2, we only need to prove that the following two conditions hold:

    (ⅰ) For every ball B=B(x0,r) with r<ρ(x0), then

    B|Iαf(x)c|dxC|B|αnχBp()fp();

    (ⅱ) For any x0Rn, then

    B(x0,ρ(x0))Iα(|f|)(x)dxC|B(x0,ρ(x0))|αnχB(x0,ρ(x0))p()fp().

    In order to check the conditions (ⅰ) and (ⅱ) above, we need to find the accurate estimate of kernel Kα(x,y) of fractional integral operator Iα (see Lemmas 2.8 and 2.9, then use them to obtain the proof of this theorem; the proof of the Theorem 1.4 proceeds identically).

    The paper is organized as follows. In Section 2, we give some important lemmas. In Section 3, we are devoted to proving Theorems 1.3 and 1.4.

    Throughout this article, C always means a positive constant independent of the main parameters, which may not be the same in each occurrence. B(x,r)={yRn:|xy|<r}, Bk=B(x0,2kR) and χBk are the characteristic functions of the set Bk for kZ. |S| denotes the Lebesgue measure of S. fg means C1gfCg.

    In this section, we give several useful lemmas that are used frequently later on.

    Lemma 2.1. [9] Assume that the exponent function p()P(Rn). If fLp()(Rn) and gLp()(Rn), then

    Rn|f(x)g(x)|dxrpfLp()(Rn)gLp()(Rn),

    where rp=1+1/p1/p+.

    Lemma 2.2. [8] Assume that p()Plog(Rn) and 1<pp+<, and p(x)p() when |x|>r0>1. For every ball B and fLp(), we have

    B|f(x)|dxCfLp(),χBLp(),

    where the constant C only depends on r0.

    Fo the following lemma see Corollary 4.5.9 in [10].

    Lemma 2.3. Let p()Plog(Rn), then for every ball BRn we have

    χBp()|B|1p(x),if|B|2n,xB,

    and

    χBp()|B|1p(),if|B|1.

    Lemma 2.4. Assume that p()Plog(Rn), then for all balls B and all measurable subsets S:=B(x0,r0)B:=B(x1,r1) we have

    χSp()χBp()C(|S||B|)11p,   χBp()χSp()C(|B||S|)11p+. (2.1)

    Proof. We only prove the first inequality in (2.1), and the second inequality in (2.1) proceeds identically. We consider three cases below by applying Lemma 2.3, and it holds that

    1) if |S|<1<|B|, then χSp()χBp()|S|1p(xS)|B|1p()(|S||B|)1(p)+=(|S||B|)11p;

    2) if 1|S|<|B|, then χSp()χBp()|S|1p()|B|1p()(|S||B|)1(p)+=(|S||B|)11p;

    3) if |S|<|B|<1, then χSp()χBp()|S|1p(xS)|B|1p(xS)|B|1p(xS)1p(xB)C(|S||B|)1(p)+=C(|S||B|)11p, where xSS and xBB.

    Indeed, since |xBxS|2r1, by using the local-Hölder continuity of p(x) we have

    |1p(xS)1p(xB)|log1r1log1r1log(e+1|xSxB|)log1r1log(e+12r1)C.

    We end the proof of this lemma.

    Remark 2.1. Thanks to the second inequality in (2.1), it is easy to prove that

    χ2Bp()CχBp().

    Lemma 2.5. [1] Suppose that the potential VBq with qn/2, then there exists positive constants C and k0 such that

    1) ρ(x)ρ(y) when |xy|Cρ(x);

    2) C1ρ(x)(1+|xy|ρ(x))k0ρ(y)Cρ(x)(1+|xy|ρ(x))k0/(k0+1).

    Lemma 2.6. [11] There exists a sequence of points {xk}k=1 in Rn such that Bk:=B(xk,ρ(xk)) satisfies

    1) Rn=kBk,

    2) For every k1, then there exists N1 such that card {j:4Bj4Bk}N.

    Lemma 2.7. Assume that p()P(Rn) and 0<α<n. Let sequence {xk}k=1 satisfy the propositions of Lemma 2.6. Then a function fLipLα,p()(Rn) if and only if f satisfies (1.4) for every ball, and

    1|B(xk,ρ(xk))|αnχB(xk,ρ(xk))p()B(xk,ρ(xk))|f(x)|dxC,forallk1. (2.2)

    Proof. Let B:=B(x,R) denote a ball with center x and radius R>ρ(x). Noting that f satisfies (1.4), and thanks to Lemma 2.6 we obtain that the set G={k:BBk} is finite.

    Applying Lemma 2.5, if zBkB, we get

    ρ(xk)Cρ(z)(1+|xkz|ρ(xk))k0C2k0ρ(z)C2k0ρ(x)(1+|xz|ρ(x))k0k0+1C2k0ρ(x)(1+Rρ(x))C2k0R.

    Thus, for every kG, we have BkCB.

    Thanks to Lemmas 2.4 and 2.6, it holds that

    B|f(x)|dx=BkBk|f(x)|dx=kG(BBk)|f(x)|dxkGBBk|f(x)|dxkGBk|f(x)|dxCkG|Bk|αnχBkp()C|B|αnχBp().

    The proof of this lemma is completed.

    Corollary 1. Assume that p()P(Rn) and 0<α<n, then a measurable function fLipLα,p() if and only if f satisfies (1.4) for every ball B(x,R) with radius R<ρ(x) and

    1|B(x,ρ(x))|αnχB(x,ρ(x))p()B(x,ρ(x))|f(x)|dxC. (2.3)

    Let kt(x,y) denote the kernel of heat semigroup etL associated to L, and Kα(x,y) be the kernel of fractional integral operator Iα, then it holds that

    Kα(x,y)=0kt(x,y)tα2dt. (2.4)

    Some estimates of kt are presented below.

    Lemma 2.8. [12] There exists a constant C such that for N>0,

    kt(x,y)Ctn/2e|xy|2Ct(1+tρ(x)+tρ(y))N,x,yRn.

    Lemma 2.9. [13] Let 0<δ<min(1,2nq). If |xx0|<t, then for N>0 the kernel kt(x,y) defined in (2.4) satisfies

    |kt(x,y)kt(x0,y)|C(|xx0|t)δtn/2e|xy|2Ct(1+tρ(x)+tρ(y))N,

    for all x,y and x0 in Rn.

    In this section, we are devoted to the proof of Theorems 1.3 and 1.4. To prove Theorem 1.3, thanks to Corollary 1 and Remark 1.2, we only need to prove that the following two conditions hold:

    (ⅰ) For every ball B=B(x0,r) with r<ρ(x0), then

    B|Iαf(x)c|dxC|B|αnχBp()fp();

    (ⅱ) For any x0Rn, then

    B(x0,ρ(x0))Iα(|f|)(x)dxC|B(x0,ρ(x0))|αnχB(x0,ρ(x0))p()fp().

    We now begin to check that these conditions hold. First, we prove (ⅱ).

    Assume that B=B(x0,R) and R=ρ(x0). We write f=f1+f2, where f1=fχ2B and f2=fχRn2B. Hence, by the inequality (1.3), we have

    BIα(|f1|)(x)dx=BIα(|fχ2B|)(x)dxCB2B|f(y)||xy|nαdydx.

    Applying Tonelli theorem, Lemma 2.1 and Remark 1.2, we get the following estimate

    BIα(|f1|)(x)dxC2B|f(y)|Bdx|xy|nαdyCRα2B|f(y)|dyC|B|αnχBp()fp(). (3.1)

    To deal with f2, let xB and we split Iαf2 as follows:

    Iαf2(x)=R20etLf2(x)tα21dt+R2etLf2(x)tα21dt:=I1+I2.

    For I1, if xB and yRn2B, we note that |x0y|<|x0x|+|xy|<C|xy|. By Lemma 2.8, it holds that

    I1=|R20Rn2Bkt(x,y)f(y)dytα21dt|CR20Rn2Btn2e|xy|2t|f(y)|dytα21dtCR20tn+α21Rn2B(t|xy|2)M/2|f(y)|dydtCR20tMn+α21dtRn2B|f(y)||x0y|Mdy,

    where the constant C only depends the constant M.

    Applying Lemma 2.1 to the last integral, we get

    Rn2B|f(y)||x0y|Mdy=i=12i+1B2iB|f(y)||x0y|Mdyi=1(2iR)M2i+1B|f(y)|dyCi=1(2iR)Mχ2i+1Bp()fp().

    By using Lemma 2.4, we arrive at the inequality

    Rn2B|f(y)||x0y|MdyCi=1(R)M(2i)nnp+MχBp()fp()CRMfp()χBp(). (3.2)

    Here, the series above converges when M>nnp+. Hence, for such M,

    |R20etLf2(x)tα21dt|CRMfp()χBp()R20tMn+α21dtC|B|αn1fp()χBp().

    For I2, thanks to Lemma 2.8, we may choose M as above and NM, then it holds that

    |R2etLf2(x)tα22dt|=|R2Rn2Bkt(x,y)f(y)dytα22dt|CR2Rn2BtαnN22ρ(x)Ne|xy|2t|f(y)|dydtCρ(x)NR2tαnN22Rn2B(t|xy|2)M/2|f(y)|dydt.

    As xB, thanks to Lemma 2.5, ρ(x)ρ(x0)=R. Hence we have

    |R2etLf2(x)tα21dt|CRNR2tM+αnN21dtRn2B|f(y)||x0y|Mdy.

    Since M+αnN<0, the integral above for variable t converges, and by applying inequality (3.2) we have

    |R2etLf2(x)tα21dt|C|B|αn1fp()χBp(),

    thus we have proved (ⅱ).

    We now begin to prove that the condition (ⅰ) holds. Let B=B(x0,r) and r<ρ(x0). We set f=f1+f2 with f1=fχ2B and f2=fχRn2B. We write

    cr=r2etLf2(x0)tα21dt. (3.3)

    Thanks to (3.1), it holds that

    B|Iα(f(x))cr|BIα(|f1|)(x)dx+B|Iα(f2)(x)cr|dxC|B|αn1χBp()fp()+B|Iα(f2)(x)cr|dx.

    Let xB and we split Iαf2(x) as follows:

    Iαf2(x)=r20etLf2(x)tα21dt+r2etLf2(x)tα21dt:=I3+I4.

    For I3, by the same argument it holds that

    I3=|r20etLf2(x)tα21dt|C|B|αn1fp()χBp().

    For I4, by Lemma 2.9 and (3.3), it follows that

    |r2etLf2(x)tα21dtcr|r2Rn2B|kt(x,y)kt(x0,y)||f(y)|dytα21dtCδr2Rn2B(|xx0|t)δtn/2e|xy|2Ct|f(y)|dytα21dtCδrδRn2B|f(y)|r2t(nα+δ)/2e|xy|2Ctdttdy.

    Let s=|xy|2t, then we obtain the following estimate

    |r2etLf2(x)tα21dtcr|CδrδRn2B|f(y)||xy|nα+δdy0snα+δ2esCdss.

    Notice that the integral above for variable s is finite, thus we only need to compute the integral above for variable y. Thanks to inequality (3.2), it follows that

    |r2etLf2(x)tα21dtcr|CδrδRn2B|f(y)||xy|nα+δdyCi=1Rαn(2i)αnp+δχBp()fp()C|B|αnnfp()χBp(),

    so (ⅰ) is proved.

    Remark 3.1. By the same argument as the proof of Theorem 1.3, thanks to Lemma 2.2 we immediately obtained that the conclusions of Theorem 1.4 hold.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    Ping Li is partially supported by NSFC (No. 12371136). The authors would like to thank the anonymous referees for carefully reading the manuscript and providing valuable suggestions, which substantially helped in improving the quality of this paper. We also thank Professor Meng Qu for his useful discussions.

    The authors declare there are no conflicts of interest.



    [1] B. Fatima, M. Yavuz, M. ur Rahman, and F.S. Al-Duais. Modeling the epidemic trend of middle eastern respiratory syndrome coronavirus with optimal control. Math. Biosci. Eng., 20 (2023), 11847–11874. http://doi.org/10.3934/mbe.2023527 doi: 10.3934/mbe.2023527
    [2] F. Evirgen, E. Uçar, S. Uçar, N. ¨Ozdemir, Modelling influenza a disease dynamics under Caputo–Fabrizio fractional derivative with distinct contact rates, Math. Mod. Numer. Simul. Appl., 3 (2023), 58–73. https://doi.org/10.53391/mmnsa.1274004 doi: 10.53391/mmnsa.1274004
    [3] H. Joshi, M. Yavuz, S. Townley, B. K. Jha, Stability analysis of a non-singular fractional-order covid-19 model with nonlinear incidence and treatment rate, Phys. Scr., 98 (2023), 045216. https://doi.org/10.1088/1402-4896/acbe7a doi: 10.1088/1402-4896/acbe7a
    [4] A. G. C. Pérez, D. A. Oluyori, A model for COVID-19 and bacterial pneumonia coinfection with community- and hospital-acquired infections, Math. Mod. Numer. Simul. Appl., 2 (2022), 197–210. https://doi.org/10.53391/mmnsa.2022.016 doi: 10.53391/mmnsa.2022.016
    [5] A. O. Atede, A. Omame, S. C. Inyama, A fractional order vaccination model for COVID-19 incorporating environmental transmission: a case study using Nigerian data, Bull. Math. Biol., 1 (2023), 78–110. https://doi.org/10.59292/bulletinbiomath.2023005 doi: 10.59292/bulletinbiomath.2023005
    [6] J. A. Cui, X. X. Mu, H. Wan, Saturation recovery leads to multiple endemic equilibria and backward bifurcation, J. Theoret. Biol., 254 (2008), 275–283. https://doi.org/10.1016/j.jtbi.2008.05.015 doi: 10.1016/j.jtbi.2008.05.015
    [7] H. Wan, J. A. Cui, Rich dynamics of an epidemic model with saturation recovery, J. Appl. Math., 2013 (2013), 314958. https://doi.org/10.1155/2013/314958 doi: 10.1155/2013/314958
    [8] N. C. Grassly, C. Fraser, Mathematical models of infectious disease transmission, Nat. Rev. Microbiol., 6 (2008), 477–487. https://doi.org/10.1038/nrmicro1845 doi: 10.1038/nrmicro1845
    [9] G. R. Jiang, Q. G. Yang, Periodic solutions and bifurcation in an SIS epidemic model with birth pulses, Math. Comput. Model., 50 (2009), 498–508. https://doi.org/10.1016/j.mcm.2009.04.021 doi: 10.1016/j.mcm.2009.04.021
    [10] J. Yang, S. Y. Tang, Holling type Ⅱ predator-prey model with nonlinear pulse as state-dependent feedback control, J. Comput. Appl. Math., 291 (2016), 225–241. https://doi.org/10.1016/j.cam.2015.01.017 doi: 10.1016/j.cam.2015.01.017
    [11] Q. Q. Zhang, S. Y. Tang, X. F. Zou, Rich dynamics of a predator-prey system with state-dependent impulsive controls switching between two means, J. Differ. Equation, 364 (2023), 336–377. https://doi.org/10.1016/j.jde.2023.03.030 doi: 10.1016/j.jde.2023.03.030
    [12] Q. Q. Zhang, B. Tang, S. Y. Tang, Vaccination threshold size and backward bifurcation of SIR model with state-dependent pulse control, J. Theoret. Biol., 455 (2018), 75–85. https://doi.org/10.1016/j.jtbi.2018.07.010 doi: 10.1016/j.jtbi.2018.07.010
    [13] T. Y. Cheng, S. Y. Tang, R. A. Cheke, Threshold dynamics and bifurcation of a state-dependent feedback nonlinear control Susceptible-Infected-Recovered model, J. Comput. Dyn., 14 (2019), 1–14. https://doi.org/10.1115/1.4043001 doi: 10.1115/1.4043001
    [14] S. Y. Tang, Y. N. Xiao, D. Clancy, New modelling approach concerning integrated disease control and cost-effectivity, Nonlinear Anal., 63 (2005), 439–471. https://doi.org/10.1016/j.na.2005.05.029 doi: 10.1016/j.na.2005.05.029
    [15] S. Y. Tang, Y. N. Xiao, L. S. Chen, R. A. Cheke, Integrated pest management models and their dynamical behaviour, Bull. Math. Biol., 67 (2005), 115–135. https://doi.org/10.1016/j.bulm.2004.06.005 doi: 10.1016/j.bulm.2004.06.005
    [16] L. F. Nie, Z. D. Teng, B. Z. Guo, A state dependent pulse control strategy for a SIRS epidemic system, Bull. Math. Biol., 75 (2013), 1697–1715. https://doi.org/10.1007/s11538-013-9865-y doi: 10.1007/s11538-013-9865-y
    [17] S. Y. Tang, B. Tang, A. L. Wang, Y. N. Xiao, Holling Ⅱ predator-prey impulsive semi-dynamic model with complex poincaré map, Nonlinear Dynam., 81 (2015), 1575–1596. https://doi.org/10.1007/s11071-015-2092-3 doi: 10.1007/s11071-015-2092-3
    [18] S. Y. Tang, W. H. Pang, On the continuity of the function describing the times of meeting impulsive set and its application, Math. Biosci. Eng., 14 (2017), 1399–1406. http://doi.org/10.3934/mbe.2017072 doi: 10.3934/mbe.2017072
    [19] S. Y. Tang, C. T. Li, B. Tang, X. Wang, Global dynamics of a nonlinear state-dependent feedback control ecological model with a multiple-hump discrete map, Commun. Nonlinear Sci. Numer. Simul., 79 (2019), 104900. https://doi.org/10.1016/j.cnsns.2019.104900 doi: 10.1016/j.cnsns.2019.104900
    [20] Q. Q. Zhang, B. Tang, T. Y. Cheng, S.Y. Tang, Bifurcation analysis of a generalized impulsive kolmogorov model with applications to pest and disease control, SIAM J. Appl. Math., 80 (2020), 1796–1819. https://doi.org/10.1137/19M1279320 doi: 10.1137/19M1279320
    [21] W. Li, T. H. Zhang, Y. F. Wang, H. D. Cheng, Dynamic analysis of a plankton-herbivore state-dependent impulsive model with action threshold depending on the density and its changing rate, Nonlinear Dynam., 107 (2022), 2951–2963. https://doi.org/10.1007/s11071-021-07022-w doi: 10.1007/s11071-021-07022-w
    [22] Q. Q. Zhang, S. Y. Tang, Bifurcation analysis of an ecological model with nonlinear state-dependent feedback control by poincaré map defined in phase set, Commun. Nonlinear Sci. Numer. Simul., 108 (2022), 106212. https://doi.org/10.1016/j.cnsns.2021.106212 doi: 10.1016/j.cnsns.2021.106212
    [23] Y. Tian, Y. Gao, K. B. Sun, A fishery predator-prey model with anti-predator behavior and complex dynamics induced by weighted fishing strategies, Math. Biosci. Eng., 20 (2023), 1558–1579. http://doi.org/10.3934/mbe.2023071 doi: 10.3934/mbe.2023071
    [24] Y. Tian, Y. Gao, K. B. Sun, Qualitative analysis of exponential power rate fishery model and complex dynamics guided by a discontinuous weighted fishing strategy, Commun. Nonlinear Sci. Numer. Simul., 118 (2023), 107011. https://doi.org/10.1016/j.cnsns.2022.107011 doi: 10.1016/j.cnsns.2022.107011
    [25] A. d'Onofrio, Stability properties of pulse vaccination strategy in SEIR epidemic model, Math. Biosci., 179 (2002), 57–72. https://doi.org/10.1016/S0025-5564(02)00095-0 doi: 10.1016/S0025-5564(02)00095-0
    [26] A. d'Onofrio, On pulse vaccination strategy in the SIR epidemic model with vertical transmission, Appl. Math. Lett., 18 (2005), 729–732. https://doi.org/10.1016/j.aml.2004.05.012 doi: 10.1016/j.aml.2004.05.012
    [27] P. Cull, Global stability of population models, Bull. Math. Biol., 43 (1981), 47–58. https://doi.org/10.1016/S0092-8240(81)80005-5 doi: 10.1016/S0092-8240(81)80005-5
    [28] N. Ferguson, D. Cummings, S. Cauchemez, C. Fraser, S. Riley, A. Meeyai, et al., Strategies for containing an emerging influenza pandemic in Southeast Asia, Nature, 437 (2005), 209–214. https://doi.org/10.1038/nature04017 doi: 10.1038/nature04017
    [29] B. Shulgin, L. Stone, Z. Agur, Pulse vaccination strategy in the SIR epidemic model, Bull. Math. Biol., 60 (1998), 1123–1148. https://doi.org/10.1006/S0092-8240(98)90005-2 doi: 10.1006/S0092-8240(98)90005-2
    [30] Z. Agur, L. Cojocaru, G. Mazor, R. M. Anderson, Y. L. Danon, Pulse mass measles vaccination across age cohorts, Proc. Pakistan Acad. Sci., 90 (1993), 11698–11702. https://doi.org/10.1073/pnas.90.24.11698 doi: 10.1073/pnas.90.24.11698
    [31] F. Albrecht, H. Gatzke, A. Haddad, N. Wax, The dynamics of two interacting populations, J. Math. Anal. Appl., 46 (1974), 658–670. https://doi.org/10.1016/0016-0032(74)90039-8 doi: 10.1016/0016-0032(74)90039-8
    [32] M. E. Fisher, B. S. Goh, T. L. Vincent, Some stability conditions for discrete-time single species models, Bull. Math. Biol., 41 (1979), 861–875. https://doi.org/10.1007/BF02462383 doi: 10.1007/BF02462383
    [33] S. K. Kaul, Stability and asymptotic stability in impulsive semidynamical systems, J. Appl. Math. Stoch. Anal., 7 (1994), 509–523.
    [34] S. K. Kaul, On impulsive semidynamical systems Ⅲ: Lyapunov stability, in Recent Trends in Differential Equations, World Scientific, (1992), 335–345.
    [35] D. D. Bainov, P. S. Simeonov, Impulsive differential equations: periodic solutions and applications, CRC Press, New York, 1993.
    [36] A. Lakmeche, Bifurcation of non trivial periodic solutions of impulsive differential equations arising chemotherapeutic treatment, Dynam. Contin. Discrete Impuls., 7 (2000), 265–287.
    [37] D. Singer, Stable orbits and bifurcation of maps of the interval, SIAM J. Appl. Math., 35 (1978), 260–267. https://doi.org/10.1137/0135020 doi: 10.1137/0135020
    [38] J. E. Marsden, M. McCracken, The hopf bifurcation and its applications, Springer-Verlag, New York, 2012.
    [39] J. M. Grandmont, Periodic and aperiodic behaviour in discrete One-Dimensional dynamical systems, Princeton University Press, Princeton, 1992.
    [40] S. Kaul, On impulsive semidynamical systems, J. Math. Anal. Appl., 150 (1990), 120–128.
    [41] E. M. Bonotto, M. Federson, Limit sets and the Poincaré-Bendixson Theorem in impulsive semidynamical systems, J. Differ. Equation, 244 (2008), 2334–2349. https://doi.org/10.1016/j.jde.2008.02.007 doi: 10.1016/j.jde.2008.02.007
    [42] P. S. Simeonov, D. D. Bainov, Orbital stability of the periodic solutions of autonomous systems with impulse effect, Int. J. Syst. Sci., 19 (1988), 2561–2585. https://doi.org/10.2977/prims/1195173347 doi: 10.2977/prims/1195173347
    [43] J. M. Grandmont, Nonlinear difference equations, bifurcations and chaos: An introduction, Res. Econ., 62 (2008), 122–177. https://doi.org/10.1016/j.rie.2008.06.003 doi: 10.1016/j.rie.2008.06.003
    [44] H. Zhou, S. Y. Tang, Complex dynamics and sliding bifurcations of the Filippov Lorenz–Chen system, Int. J. Bifur. Chaos, 32 (2022), 2250182. https://doi.org/10.1142/S0218127422501826 doi: 10.1142/S0218127422501826
    [45] H. Zhou, S. Y. Tang, Bifurcation dynamics on the sliding vector field of a Filippov ecological system, Appl. Math. Comput., 424 (2022), 127052. https://doi.org/10.1016/j.amc.2022.127052 doi: 10.1016/j.amc.2022.127052
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1690) PDF downloads(87) Cited by(1)

Figures and Tables

Figures(5)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog