Loading [MathJax]/jax/output/SVG/jax.js
Research article

A new approach to persistence and periodicity of logistic systems with jumps

  • Received: 02 July 2021 Accepted: 15 August 2021 Published: 25 August 2021
  • MSC : 34A37, 34D20

  • This paper considers a class of logistic type differential system with jumps. Based on discontinuous control theory, a new approach is developed to guarantee the persistence and existence of a unique globally attractive positive periodic solution. The development results of this paper emphasize the effects of jumps on system, which are different from the existing ones in the literature. Two examples and their simulations are given to illustrate the effectiveness of the proposed results.

    Citation: Kegang Zhao. A new approach to persistence and periodicity of logistic systems with jumps[J]. AIMS Mathematics, 2021, 6(11): 12245-12259. doi: 10.3934/math.2021709

    Related Papers:

    [1] Yuanfu Shao . Dynamics and optimal harvesting of a stochastic predator-prey system with regime switching, S-type distributed time delays and Lévy jumps. AIMS Mathematics, 2022, 7(3): 4068-4093. doi: 10.3934/math.2022225
    [2] Yassine Sabbar, Aeshah A. Raezah . Modeling mosquito-borne disease dynamics via stochastic differential equations and generalized tempered stable distribution. AIMS Mathematics, 2024, 9(8): 22454-22485. doi: 10.3934/math.20241092
    [3] Nana Tao, Chunxiao Ding . Global attractivity for uncertain differential systems. AIMS Mathematics, 2022, 7(2): 2142-2159. doi: 10.3934/math.2022122
    [4] Taras Lukashiv, Igor V. Malyk, Maryna Chepeleva, Petr V. Nazarov . Stability of stochastic dynamic systems of a random structure with Markov switching in the presence of concentration points. AIMS Mathematics, 2023, 8(10): 24418-24433. doi: 10.3934/math.20231245
    [5] Yaw Chang, Wei Feng, Michael Freeze, Xin Lu, Charles Smith . Elimination, permanence, and exclusion in a competition model under Allee effects. AIMS Mathematics, 2023, 8(4): 7787-7805. doi: 10.3934/math.2023391
    [6] Xiaodong Wang, Kai Wang, Zhidong Teng . Global dynamics and density function in a class of stochastic SVI epidemic models with Lévy jumps and nonlinear incidence. AIMS Mathematics, 2023, 8(2): 2829-2855. doi: 10.3934/math.2023148
    [7] Hong Qiu, Yanzhang Huo, Tianhui Ma . Dynamical analysis of a stochastic hybrid predator-prey model with Beddington-DeAngelis functional response and Lévy jumps. AIMS Mathematics, 2022, 7(8): 14492-14512. doi: 10.3934/math.2022799
    [8] Yassine Sabbar, Aeshah A. Raezah, Mohammed Moumni . Enhancing epidemic modeling: exploring heavy-tailed dynamics with the generalized tempered stable distribution. AIMS Mathematics, 2024, 9(10): 29496-29528. doi: 10.3934/math.20241429
    [9] Chunwei Wang, Jiaen Xu, Naidan Deng, Shujing Wang . Two-sided jumps risk model with proportional investment and random observation periods. AIMS Mathematics, 2023, 8(9): 22301-22318. doi: 10.3934/math.20231137
    [10] Wanjuan Du . A further study on an attraction-repulsion chemotaxis system with logistic source. AIMS Mathematics, 2024, 9(7): 16924-16930. doi: 10.3934/math.2024822
  • This paper considers a class of logistic type differential system with jumps. Based on discontinuous control theory, a new approach is developed to guarantee the persistence and existence of a unique globally attractive positive periodic solution. The development results of this paper emphasize the effects of jumps on system, which are different from the existing ones in the literature. Two examples and their simulations are given to illustrate the effectiveness of the proposed results.



    For modelling the dynamics of some biological populations, various logistic type differential systems have been proposed and studied in the past several years, see ([1,2,3,4,5,6,7,8,9,10,11]). The classical nonautonomous logistic differential system can be described by

    N(t)=N(t)r(t)[1N(t)K(t)], (1.1)

    where N is the density of population at time t; r and K admit positive upper and lower bounds which models the growth rate and the environmental carrying capacity or saturation level at time t, respectively. Since most populations are affected by the outside environment, it is necessary to study the models of the population dynamics with harvesting, and the following system and its various generalized forms are considered

    N(t)=N(t)r(t)[1N(t)K(t)]E(t,N(t)), (1.2)

    where function E is a harvesting strategy for the population. In 1959, Holling ([12]) proposed three basic types of the harvesting term:

    Type I (linear): E(t,N(t))=α(t)N+β(t),

    Type II (cyrtoid): E(t,N(t))=α(t)NN+β(t),

    Type III (sigmoid): E(t,N(t))=α(t)N2γ2(t)+β(t)N2,

    where α,β and γ are some scalar functions with positive upper and lower bounds. We refer to the monographs ([13,14,15]) for the discussions of three types harvesting. In general, population models (1.1) or models (1.2) with the above three basic types harvesting can be uniformly described by the following differential system

    N(t)=N(t)f(t,N(t)), (1.3)

    where f is a scalar function of t and N.

    Recently, many research works have paid much attention to the population models with jumps which includes impulsive harvesting and planting, since the discontinuous models governed by impulsive differential systems are more feasible and beneficial than the continuous ones at times, see ([16,17,18,19,20,21,22,23,24,25,26,27,28,29]). For instance, it is often the case that planting and harvesting of the species are intermittent or occur irregularly. Besides, continuous changes such as temperature or rainfall in environment parameters can create some discontinuous outbreaks in biological populations. For fishery management it is unreasonable to assume that fisherman to fish the whole day, and in fact they only fish for some time, and moreover, the seasons and weather variations will also affect the fishing. It has been shown ([17,29]) that the continuous harvesting policy is superior to the impulsive harvesting policy, however, the latter is more beneficial in realistic operation. Hence, it is significant to consider jumps in the investigation of population models.

    Motivated by the above discussions, the main objective of this paper is to study the logistic system (1.3) with jumps

    {N(t)=N(t)f(t,N(t)),  ttk,ΔN|t=tk=N(tk)N(tk)=Ik(N(tk)),  kZ+. (1.4)

    It shows that there exist jumps when human activities are considered continuously on population model. Especially, it reflects the combination of continuous harvesting and impulsive harvesting. Based on impulsive control theory, we shall investigate the effects of impulsive harvesting and stocking, and establish conditions for the persistence and existence of a unique globally attractive positive periodic solution of system (1.4). The development results are different from the existing ones in the literature and the relation between dynamics of population models and jumps will be emphasized in this paper.

    The remainder of this paper is organized as follows. In Section 2, some necessary definitions and preliminary results are presented. In Section 3, some new criteria for persistence and periodicity are presented. In Section 4, simulations are given to illustrate the effectiveness of the main results. Finally, conclusions are drawn in Section 5.

    Notations. Let R denote the set of real numbers, R+ the set of nonnegative real numbers, R the set of nonpositive real numbers and Z+ the set of positive integers. For any interval JR, set SRk(1kN),C(J,S)={φ:JS is continuous} and PC(J,S)={φ:JS is continuous everywhere except at finite number of points t, at which φ(t+), φ(t) exist and φ(t+)=φ(t)}. Σ+={c(t): cPC(R+,R+) and for any interval (α,β)R+,c(t)0. }. Σ={c(t): cPC(R+,R) and for any interval (α,β)R+,c(t)0. }. Given a continuous function f which is defined on ΛR, we set fIinfsΛf(s),fSsupsΛf(s). The jump times tk,kZ+, satisfy 0t0<t1<<tk+ as k+.

    System (1.4) may be rewritten as:

    {x(t)=xf(t,x),  tt0, ttk,Δx|t=tk=x(tk)x(tk)=Ik(x(tk)),  kZ+,x(t0)=x0, (2.1)

    where x0R+, fPC([t0,)×R+,R). The numbers x(tk) and x(tk) denote the population densities of the species before and after jumps at the moments tk, respectively. IkC(R+,R), which characterize the magnitude of the jumps on the species at the moments tk and satisfy Ik(s)+s>0 for any sR+,kZ+; In particular, when Ik>0, the perturbation stands for planting of the species, while Ik<0 stands for harvesting. We assume that system (2.1) satisfies some fundamental conditions which guarantee the global existence and uniqueness of the solutions on [t0,), see ([17,30]). In the following, denote by x(t)=x(t,t0,x0) the solution of system (2.1) with initial value (t0,x0).

    Definition 2.1. (see [17]) System (2.1) is said to be persistent, if there exist constants M>0 and m>0 such that each positive solution x(t)=x(t,t0,x0) of system (2.1) satisfies

    mlim inft+x(t)lim supt+x(t)M.

    Definition 2.2. (see [17]) A map x:R+R+ is said to be a ω-periodic solution of system (2.1), if

    (i) x(t) is a piecewise continuous map with first-class discontinuity points and satisfies (2.1);

    (ii) x(t) satisfies x(t+ω)=x(t),ttk and x(tk+ω+)=x(t+k),kZ+.

    Definition 2.3. (see [17]) Assume that x(t)=x(t,t0,x0) be a positive periodic solution of system (2.1). Then x is said to be globally attractive, if for any positive solution x(t)=x(t,t0,x0) of system (2.1), it holds that

    limt+x(t)x(t)=0.

    Lemma 3.1. The set R+ is the positively invariant set of system (2.1).

    Proof. Note that Ik(s)+s>0 for any sR+,kZ+. The proof of Lemma 3.1 is obvious.

    Theorem 3.1. Assume that there exist constants q>1,λ>1 and M>m>0 such that

    (i) 1qqIk(s)sλ1,  s>0;

    (ii) f(t,M)Σ and supkZ+tktk1f(s,M)ds<lnλ;

    (iii) f(t,m)Σ+ and infkZ+tktk1f(s,m)ds>lnq;

    (iv) [f(t,x)f(t,y)]sgn(xy)0,  (t,x,y)R+×R+×R+.

    Then the set

    Δ={xR+:mq3xMλ3}

    is the ultimately bounded set of system (2.1), that is, system (2.1) is persistent.

    Proof. Let x(t)=x(t,t0,x0) be the solution of system (2.1) with initial value (t0,x0), where (t0,x0)R+×R+. First, {{we show that} there exists tt0 such that x(t)M. If x0(0,M], then the conclusion is obvious. Or else, assume that x(t)>M for all tt0. Then it follows from (iv) and system (2.1) that x(t)=xf(t,x)xf(t,M),t[tk1,tk),kZ+, which, together with (i) and (ii), yields

    x(tk+1)λx(tk+1)λk+1x(t0)exp(tk+1t0f(s,M)ds)=x(t0)exp((k+1)lnλ+tk+1t0f(s,M)ds)x(t0)exp((k+1)[lnλ+μ])0,  as k,

    where μsupkZ+tktk1f(s,M)ds. This is a contradiction with the above assumption that x(t)>M for all tt0 and thus there exists tt0 such that x(t)M.

    Now we show that x(t)Mλ3,tt. Suppose not, then there exists ¯t>t such that x(¯t+)Mλ3 and x(¯t)Mλ3. Since x(t)M, there exists t_[t,¯t] such that x(t_+)M and x(t_)M. Moreover, t_<¯t. In fact, if t_=¯t, then λMλx(t_)x(t_+)Mλ3, which contradicts λ>1. Thus we get Mx(t)Mλ3,t[t_,¯t]. Then there are three cases: (a) If there is no jump on [t_,¯t], then x(t)=xf(t,x)xf(t,M)0,t[t_,¯t], which implies that Mλ2x(¯t)x(t_+)λM. This is a contradiction. (b) If there exists a jump on [t_,¯t], assume that t_<tσ<¯t, where σZ+. Then it follows from f(t,M)Σ that

    Mλ2x(¯t)λx(t_+)exp(¯tt_f(s,M)ds)λ2Mexp(¯tt_f(s,M)ds)<λ2M,

    which is also a contradiction. (c) If there exist some jumps on [t_,¯t], assume that t_<tσ<tσ+1<<tσ+l<¯t, where σ,lZ+. Then it can be deduced that

    x(tσ+l)λlx(t_+)exp(tσ+lt_f(s,M)ds),

    which lead to

    Mλ2x(¯t)x(t+σ+l)exp(¯ttσ+lf(s,M)ds)λl+2Mexp(¯tt_f(s,M)ds)λ2Mexp(llnλ+tσ+ltσf(s,M)ds)λ2Mexp(l[lnλ+μ])<λ2M.

    Obviously, this is also a contradiction and thus all cases (a-c) are impossible. Hence, it holds that x(t)Mλ3,tt.

    Next we show that there exists tt such that x(t)>m. Suppose not, then x(t)m,tt. It follows from (iv) and (2.1) that x(t)=xf(t,x)xf(t,m),t[tk1,tk)[t,),kZ+. Assume that t[tl1,tl), for some lZ+, then it can be deduced from (i) and (iii) that

    x(tl+k)1qx(tl+k)1qk+1x(t)exp(tl+ktf(s,m)ds)1qx(t)exp(klnq+tl+ktlf(s,m)ds)1qx(t)exp(k[ηlnq]),  as k,

    where ηinfkZ+tktk1f(s,m)ds. This is a contradiction with x(t)m,tt and thus there exists tt such that x(t)>m.

    Furthermore, we show that x(t)mq3,tt. Suppose not, then there exists ˆtt such that x(ˆt+)mq3 and x(ˆt)mq3. Note that x(t)>m and q>1, there exists ˜t[t,ˆt) such that x(˜t+)m and x(˜t)m. Thus we get that mq3x(t)m,t[˜t,ˆt]. Then there are also three cases: (d) If there is no jump on [˜t,ˆt], then x(t)=xf(t,x)xf(t,m)0,t[˜t,ˆt], which implies that mqx(˜t+)x(ˆt)mq2. This is a contradiction. (e) If there exists a jump on [˜t,ˆt], assume that ˜t<tρ<ˆt, where ρZ+. Then it follows from f(t,m)Σ+ that

    mq2x(ˆt)1qx(˜t+)exp(ˆt˜tf(s,m)ds)mq2exp(ˆt˜tf(s,m)ds)>mq2,

    which is a contradiction. (f) If there exist some jumps on [˜t,ˆt], assume that ˜t<tρ<tρ+1<<tρ+s<ˆt, where ρ,sZ+. Then it can be deduced that

    x(tρ+s)1qsx(˜t+)exp(tρ+s˜tf(s,m)ds),

    which lead to

    mq2x(ˆt)x(t+ρ+s)exp(ˆttρ+sf(s,m)ds)1qs+1x(˜t+)exp(ˆt˜tf(s,m)ds)1qx(˜t+)exp(slnq+ˆt˜tf(s,m)ds)mq2exp(s[ηlnq])>mq2.

    Obviously, it is also a contradiction and thus all cases (d-f) are impossible. The proof of Theorem 3.1 is there completed.

    Remark 3.1. Note that condition (iv) in Theorem 3.1 can be replaced by the following stronger one: f(t,x)x0, (t,x)R+×R+, which can be checked more easily in practical problems. In addition, from Theorem 3.1 one may note that it is possible that Ik0 or Ik0. Thus the development results can be applied to logistic systems with impulsive planting or/and impulsive harvesting. In particular, when there is no jump, we have the following two Corollaries.

    Corollary 3.1. Assume that there exist constants M>m>0 such that

    (i) f(t,M)0,tt0, and +t0f(s,M)ds;

    (ii) f(t,m)0,tt0, and +t0f(s,m)ds+;

    (iii) [f(t,x)f(t,y)]sgn(xy)0,  (t,x,y)R+×R+×R+.

    Then the set Δ={xR+:mxM} is the ultimately bounded set of system (2.1) without jumps, that is, system (2.1) without jumps is persistent.

    Corollary 3.2. Assume that there exist constants M>m>0 such that f(t,M)0 and f(t,m)0,tt0, then Δ={xR+:mxM} is the invariant set of system (2.1) without jumps.

    Remark 3.2. Corollaries 3.1 and 3.2 can be easily derived by the proof process of Theorem 3.1.

    Consider the following logistic differential system with jumps

    {x(t)=x(t)[r(t)a(t)x(t)],  tt0, ttk,Δx|t=tk=x(tk)x(tk)=Ik(x(tk)),  kZ+,x(t0)=x0, (3.1)

    where x0R+, r,aPC(R+,R+) and admit positive upper and lower bounds which are natural for biological meanings.

    Corollary 3.3. Suppose that there exist constants q>1,λ>1 and M>m>0 such that

    (i) 1qqIk(s)sλ1,  s>0;

    (ii) m<rIaS and M>rSaI;

    (iii) infkZ+{tktk1}>max{lnqrIaSm, lnλaIMrS}.

    Then the set

    Δ={xR+:mq3xMλ3}

    is the ultimately bounded set of system (3.1).

    Corollary 3.4. Suppose that there exist constants q>1,λ>1 such that

    (i) 1qqIk(s)sλ1,  s>0;

    (ii) infkZ+{tktk1}>lnqrI.

    Then system (3.1) is persistent.

    Consider the following logistic differential system with jumps and cyrtoid type harvesting

    {x(t)=x(t)[r(t)a(t)x(t)]α(t)x(t)x(t)+β(t),  tt0, ttk,Δx|t=tk=x(tk)x(tk)=Ik(x(tk)),  kZ+,x(t0)=x0, (3.2)

    where x0R+, r,a,α,βPC(R+,R+) and admit positive upper and lower bounds which are natural for biological meanings.

    Corollary 3.5. Suppose that rIβI>αS and aI(βI)2>αS. Moreover, there exist constants q>1,λ>1 and M>m>0 such that

    (i) 1qqIk(s)sλ1,  s>0;

    (ii) rSaIMαIM+βS<0;

    (iii) rIaSmαSm+βI>0;

    (iii) infkZ+{tktk1}>max{lnqrIaSmαSm+βI, lnλaIM+αIM+βSrS}.

    Then the set

    Δ={xR+:mq3xMλ3}

    is the ultimately bounded set of system (3.2).

    Corollary 3.6. Suppose that rIβI>αS and aI(βI)2>αS. Moreover, there exist constants q>1,λ>1 such that

    (i) 1qqIk(s)sλ1,  s>0;

    (ii) infkZ+{tktk1}>lnqrIαSβI.

    Then system (3.2) is persistent.

    Remark 3.3. Corollaries 3.3 and 3.5 can be directly derived by Theorem 3.1. For Corollaries 3.4 and 3.6, {{one only needs} choose small enough m>0 and large enough M>0 such that (ii) and (iii) in the corresponding Corollary hold.

    In the following, we shall investigate the periodic solution problem of system (2.1).

    Theorem 3.2. Assume that there exist constants q>1,m>0,ω>0 and θZ+ such that

    (i) Ik(s)=(ρk1)s,sR+, where ρk1q, kZ+;

    (ii) f(t,m)Σ+ and infkZ+tktk1f(s,m)ds>lnq;

    (iii) f(t+ω,)=f(t,), tk+ω=tk+θ and ρk+θ=ρk,kZ+;

    (iv) f(t,x)f(t,y)xyp(t), where pPC(R+,R+) and satisfies t0p(s)ds=+.

    Then system (2.1) has a unique positive ω-periodic solution, which is globally attractive.

    Proof. Let y(t)=lnx(t), then system (2.1) may be rewritten as

    {y(t)=f(t,ey(t)),  tt0, ttk,y(tk)y(tk)=lnρk,  kZ+,y(t0)=lnx0R. (3.3)

    Obviously, the investigation of the positive periodic solution problem for (2.1) is equal to investigate the periodic solution problem for system (3.3).

    First, by condition (iii) it is easy to show that system (3.3) has an ω-periodic solution if there exists a y0R such that y(t0+ω,t0,y0)=y0, where y(t,t0,y0) is the solution of system (3.3) through (t0,y0). In fact, one may establish a solution as follows:

    ˜y(t)={y(t),t[t0,t0+ω],y(tnω),t[t0+nω,t0+(n+1)ω].

    Obviously, ˜y is ω-periodic. Next we show that ˜y is a solution of system (3.3). For any t[t0+ nω,t0+(n+1)ω], if ttk,kZ+, then

    ˜y(t)=y(tnω)=f(tnω,ey(tnω))=f(t,ey(tnω))=f(t,e˜y(t));

    if t=tm for some mZ+, then it follows from (iii) that tmnω=tmnθ, which yields that

    ˜y(tm)=y(tmnθ)=y(tmnθ)+lnρmnθ=y(tmnω)+lnρm=˜y(tm)+lnρm.

    By the above discussion, ˜y is a solution of system (3.3). Moreover, by the existence-uniqueness theorem, ˜yy,tt0. That is, y is an ω-periodic solution of system (3.3).

    Let z(t)=z(t,t0,z0) and h(t)=h(t,t0,h0) be two solutions of system (3.3) through (t0,z0) and (t0,h0), respectively, where z0h0. Define Γ(t)=|h(t)z(t)|,tt0. Then

    Γ(tk)=|h(tk)z(tk)|=|h(tk)z(tk)|=Γ(tk), kZ+. (3.4)

    In addition, it can be deduced that

    Γ(t)=[f(t,eh(t))f(t,ez(t))]sgn(hz)=f(t,eh(t))f(t,ez(t))eh(t)ez(t)|eh(t)ez(t)|p(t)eξ(t)Γ(t),tt0,

    where ξ(t) is a real value between h(t) and z(t). By conditions (i), (ii) and Theorem 3.1, we know that there exists T0t0 such that eξ(t)mq3,tT0, which yields that

    Γ(t)p(t)mq3Γ(t),tT0. (3.5)

    For given T0t0, define

    G(T0)|h(T0)z(T0)||h(t0)z(t0)|. (3.6)

    Considering (3.4), (3.5) and (3.6), we have

    Γ(t)=|h(t)z(t)|G(T0)|h(t0)z(t0)|exp(mq3tT0p(s)ds),tT0, (3.7)

    which implies that there exists T1T0 such that

    |h(t)z(t)|12|h(t0)z(t0)|,tT1. (3.8)

    Define an operator

    F:r0r(t0+ω,t0,r0),

    where r(t,t0,r0) is the solution of system (3.3) through (t0,r0). Obviously, operator F maps the set R into itself. Moreover,

    Fk(r0)=r(t0+kω,t0,r0),kZ+.

    Let k large enough such that t0+kωT1, then it follows from (3.8) that

    |Fk(r0)Fk(s0)|12|r0s0|,

    where s0R. Thus, operator F is a contraction mapping in Banach space. Using Banach fixed point theorem, there exists a unique r0R such that F(r0)=r0. Hence, system (3.3) has a ω-periodic solution r(t)=r(t,t0,r0). That is, system (2.1) has a positive ω-periodic solution x(t)=er(t).

    Next we show that r(t)=r(t,t0,r0) is the unique ω-periodic solution of system (3.3) and all other solutions converge to it. Suppose that r(t)=r(t,t0,r0) is any another solution of system (3.3) through (t0,r0), then it follows from (3.7) that

    |r(t)r(t)|G(T0)|r0r0|exp(mq3tT0p(s)ds),tT0,

    where

    G(T0)|r(T0)r(T0)||r0r0|.

    It is obvious that |r(t)r(t)|0, as t, which implies that system (2.1) has a unique positive ω-periodic solution, which is globally attractive. The proof of Theorem 3.1 is completed.

    Remark 3.4. One may observe from Theorem 3.2 that, to investigate the periodic problem of system (2.1), there is no any restriction on the upper bound of jump constant ρk, that is, ρk may be large enough if it satisfies the periodic condition ρk+θ=ρk. Moreover, condition ρk1q in Theorem 3.2 can be replaced by min1iθρi1q.

    For system (3.1) and (3.2), we have

    Corollary 3.7. Suppose that there exist constants q>1,ω>0 and θZ+ such that

    (i) Ik(s)=(ρk1)s,sR+, where min1iθρi1q;

    (ii) min1iθ{titi1}>lnqrI;

    (iii) r(t+ω)=r(t),a(t+ω)=a(t), tk+ω=tk+θ and ρk+θ=ρk,kZ+.

    Then system (3.1) has a unique positive ω-periodic solution, which is globally attractive.

    Corollary 3.8. Suppose that rIβI>αS and aI(βI)2>αS. Moreover, there exist constants q>1,ω>0 and θZ+ such that

    (i) Ik(s)=(ρk1)s,sR+, where min1iθρi1q;

    (ii) min1iθ{titi1}>lnqrIαSβI;

    (iii) r(t+ω)=r(t),a(t+ω)=a(t), tk+ω=tk+θ and ρk+θ=ρk,kZ+.

    Then system (3.2) has a unique positive ω-periodic solution, which is globally attractive.

    Remark 3.5. The research thought in the paper is completely new and can be extended to the investigation of jump for delay logistic differential system.

    In this section, two examples and their simulations are presented to show the effectiveness of our obtained results.

    Example 4.1. Consider the logistic type differential equations with jumps:

    {x(t)=x[ln(7+2t1+t2)1+t21+t+t2ln(1+x2)],t0,ttk,x(t+k)=γx(tk),kZ+,x(0)=x0R+, (4.1)

    where γ>0 is a given constant. For system (4.1), we have

    Property 4.1. Case γ<1. System (4.1) is persistent, if infkZ+{tktk1}>logγ7.

    Proof. Let

    f(t,x)=ln(7+2t1+t2)1+t21+t+t2ln(1+x2),(t,x)R+×R+.

    Then it is easy to check that

    ln7ln(1+x2)f(t,x)ln923ln(1+x2).

    Since infkZ+{tktk1}>logγ7, there exists m>0 small enough such that

    infkZ+{tktk1}>lnγln71+m2>0. (4.2)

    For given m, let λ=2, then there exists M>m large enough such that

    0<ln223ln(1+M2)ln9<lnγln71+m2. (4.3)

    Choose q=1γ and considering (4.2), (4.3), it is easy to check that all conditions in Theorem 3.1 hold and thus system (4.1) is persistent for the case γ<1.

    Property 4.2. Case γ>1. System (4.1) is persistent, if infkZ+{tktk1}>0.

    Proof. Choose q(1,γ] and m>0 small enough such that

    0<lnqln71+m2<η, (4.4)

    where ηinfkZ+{tktk1}>0. Let λ=γ>1 and then choose M>m large enough such that

    lnλ23ln(1+M2)ln9<η. (4.5)

    By (4.4) and (4.5), it is easy to check that all conditions in Theorem 3.1 hold and thus system (4.1) is persistent for the case γ>1.

    Property 4.3. Case γ=1. System (4.1) is persistent.

    Proof. In view of Corollary 3.1 and the analysis of Properties 4.1 and 4.2, the above property is obvious.

    Remark 4.1. In particular, if let γ=0.5, then by Property 4.1 system (4.1) is persistent, if infkZ+{tktk1}>0.36. For example, when tk=0.4k, Figure 1(a) shows that system (4.1) is persistent. However, when tk=0.3k, Property 4.1 is invalid. In this case, it is interesting to see that system (4.1) will become extinct, which is shown in Figure 1(b). This partially reflects the advantage of our development results. In addition, if let γ=2 and 8, then by Property 4.2 system (4.1) is persistent, if infkZ+{tktk1}>0. Figure 1(c) shows the case that γ=2,tk=0.2k, Figure 1(d) shows the case that γ=8,tk=k.

    Figure 1.  (a) State trajectory of system (4.1) with γ=0.5 and tk=0.4k; (b) State trajectory of system (4.1) with γ=0.5 and tk=0.3k; (c) State trajectory of system (4.1) with γ=2 and tk=0.2k; (d) State trajectory of system (4.1) with γ=8 and tk=k.

    Example 4.2. Consider the logistic type differential equations with jumps:

    {x(t)=x[r+0.1sin2πωt(1+0.2cos2πωt)x],t0,ttk,x(t+k)=γkx(tk),kZ+,x(0)=x0R+, (4.6)

    where tk=kωθ;r>0.1,ω>0 and θZ+ are some given constants; γk>0 satisfying γk+θ=γk,kZ+. For system (4.6), we have

    Property 4.4. System (4.6) has a unique positive ω-periodic solution, which is globally attractive, if

    min1iθγi<1  and  1min1iθγi<exp[ωθ(r0.1)].

    Property 4.5. System (4.6) has a unique positive ω-periodic solution, which is globally attractive, if min1iθγi>1.

    Proof. By Corollary 3.6, then above two Properties can be easily derived.

    Remark 4.2. In particular, if let γk=0.8,r=0.4 and ω=θ, then by Property 4.4 system (4.6) has a unique globally attractive positive ω-periodic solution. When ω=2 and 5, the corresponding simulations are given in Figure 2(a) and 2(b), respectively. Under the same conditions, if let γk=2, then by Property 4.5 system (4.6) has a unique globally attractive positive ω-periodic solution, which are shown in Figure 2(c) and 2(d) for ω=1 and 10, respectively.

    Figure 2.  (a) State trajectories of system (4.6) with γk=0.8 and ω=2; (b) State trajectories of system (4.6) with γk=0.8 and ω=5; (c) State trajectories of system (4.6) with γk=2 and ω=1; (d) State trajectories of system (4.6) with γk=2 and ω=10.

    In this paper, we investigated a class of logistic type differential system with jumps. Based on impulsive control theory, some new sufficient condition ensuring the permanence and existence of a unique globally attractive positive periodic solution were derived. The developed method is different from the usual methods in other literatures. Two numerical examples were given to illustrate the effectiveness and advantages of the results. The research thought in the paper can be extended to the investigation of jump for impulsive logistic differential system with time delays. In the near future, we shall do some further research on this topic.

    The author declares no conflict of interest.



    [1] K. Gopalsamy, Stability and oscillations in delay differential equations of population dynamics, Springer Science & Business Media, 1992.
    [2] G. Consolini, M. Materassi, A stretched logistic equation for pandemic spreading, Chaos Soliton. Fract., 140 (2020), 110113. doi: 10.1016/j.chaos.2020.110113
    [3] C. Burgos, J. C. Cortés, E. López-Navarro, R. J. Villanueva, Probabilistic analysis of linear-quadratic logistic-type models with hybrid uncertainties via probability density functions, AIMS Mathematics, 6 (2021), 4938–4957. doi: 10.3934/math.2021290
    [4] G. Ascione, Abstract Cauchy problems for generalized fractional calculus, Nonlinear Anal., 209 (2021), 112339. doi: 10.1016/j.na.2021.112339
    [5] W. Cintra, C. Morales-Rodrigo, A. Suárez, Refuge versus dispersion in the logistic equation, J. Differ. Equations, 262 (2017), 5606–5634. doi: 10.1016/j.jde.2017.02.012
    [6] I. Area, J. Losada, J. J. Nieto, A note on the fractional logistic equation, Physica A, 444 (2016), 182–187. doi: 10.1016/j.physa.2015.10.037
    [7] M. Delgado, M. Molina-Becerra, A. Suárez, A logistic type equation in RN with a nonlocal reaction term via bifurcation method, J. Math. Anal. Appl., 493 (2021), 124532. doi: 10.1016/j.jmaa.2020.124532
    [8] M. Liu, K. Wang, On a stochastic logistic equation with impulsive perturbations, Comput. Math. Appl., 63 (2012), 871–886.
    [9] P. Korman, On separated solutions of logistic population equation with harvesting, J. Math. Anal. Appl., 455 (2017), 2024–2029. doi: 10.1016/j.jmaa.2017.06.029
    [10] M. D'Ovidio, Non-local logistic equations from the probability viewpoint, arXiv: 2105.00779.
    [11] M. D'Ovidio, Solutions of fractional logistic equations by Euler's numbers, Physica A, 506 (2018), 1081–1092. doi: 10.1016/j.physa.2018.05.030
    [12] C. Holling, Some characteristics of simple types of predation and parasitism, Can. Entomol., 91 (1959), 385–398. doi: 10.4039/Ent91385-7
    [13] A. Hasting, Population biology: Concepts and models, New York: Springer, 1997.
    [14] M. Kot, Elements of mathematical ecology, Cambridge University Press, 2001.
    [15] F. Brauer, C. Castilo-Chavez, Mathematical models in population biology and epidemiology, New York: Springer, 2001.
    [16] X. Yang, X. Li, Q. Xi, P. Duan, Review of stability and stabilization for impulsive delayed systems, Math. Biosci. Eng., 15 (2018), 1495–1515. doi: 10.3934/mbe.2018069
    [17] D. Bainov, P. Simeonov, Impulsive differential equations: Periodic solutions and applications, CRC Press, 1993.
    [18] Z. Li, X. Shu, F. Xu, The existence of upper and lower solutions to second order random impulsive differential equation with boundary value problem, AIMS Mathematics, 5 (2020), 6189–6210. doi: 10.3934/math.2020398
    [19] X. Li, J. Shen, R. Rakkiyappan, Persistent impulsive effects on stability of functional differential equations with finite or infinite delay, Appl. Math. Comput., 329 (2018), 14–22.
    [20] G. Ballinger, X. Liu, Permanence of population growth models with impulsive effects, Math. Comput. Model., 26 (1997), 59–72.
    [21] I. Stamova, G. T. Stamov, Applied impulsive mathematical models, Switzerland: Springer, 2016.
    [22] Z. Teng, L. Nie, X. Fang, The periodic solutions for general periodic impulsive population systems of functional differential equations and its applications, Comput. Math. Appl., 61 (2011), 2690–2703. doi: 10.1016/j.camwa.2011.03.023
    [23] D. Yang, X. Li, Z. Liu, J. Cao, Persistence of nonautonomous logistic system with time-varying delays and impulsive perturbations, Nonlinear Anal. Model., 25 (2020), 564–579.
    [24] A. N. Churilov, Orbital stability of periodic solutions of an impulsive system with a linear continuous-time part, AIMS Mathematics, 5 (2020), 96–110. doi: 10.3934/math.2020007
    [25] G. Stamov, E. Gospodinova, I. Stamova, Practical exponential stability with respect to h-manifolds of discontinuous delayed Cohen-Grossberg neural networks with variable impulsive perturbations, Math. Model. Control, 1 (2021), 26–34. doi: 10.3934/mmc.2021003
    [26] P. Li, X. Li, J. Cao, Input-to-state stability of nonlinear switched systems via Lyapunov method involving indefinite derivative, Complexity, 2018 (2018), 8701219.
    [27] M. Liu, S. Li, X. Li, L. Jin, C. Yi, Intelligent controllers for multirobot competitive and dynamic tracking, Complexity, 2018 (2018), 4573631.
    [28] T. Wei, X. Xie, X. Li, Input-to-state stability of delayed reaction-diffusion neural networks with multiple impulses, AIMS Mathematics, 6 (2021), 5786–5800. doi: 10.3934/math.2021342
    [29] X. Zhang, Z. Shuai, K. Wang, Optimal impulsive harvesting policy for single population, Nonlinear Anal. Real, 4 (2003), 639–651. doi: 10.1016/S1468-1218(02)00084-6
    [30] V. Lakshmikantham, D. D. Bainov, P. S. Simeonov, Theory of impulsive differential equations, World Scientific, 1989.
  • Reader Comments
  • © 2021 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(2341) PDF downloads(78) Cited by(0)

Figures and Tables

Figures(2)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog