Hindawi Discrete Dynamics in Nature and Society Volume 2017, Article ID 7037245, 15 pages https://doi.org/10.1155/2017/7037245
Research Article Coexistence for an Almost Periodic Predator-Prey Model with Intermittent Predation Driven by Discontinuous Prey Dispersal Yantao Luo, Long Zhang, Zhidong Teng, and Tingting Zheng College of Mathematics and System Science, Xinjiang University, Urumqi 830046, China Correspondence should be addressed to Long Zhang; longzhang
[email protected] Received 26 July 2017; Accepted 23 October 2017; Published 6 December 2017 Academic Editor: Guang Zhang Copyright Β© 2017 Yantao Luo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. An almost periodic predator-prey model with intermittent predation and prey discontinuous dispersal is studied in this paper, which differs from the classical continuous and impulsive dispersal predator-prey models. The intermittent predation behavior of the predator species only happens in the channels between two patches where the discontinuous migration movement of the prey species occurs. Using analytic approaches and comparison theorems of the impulsive differential equations, sufficient criteria on the boundedness, permanence, and coexistence for this system are established. Finally, numerical simulations demonstrate that, for an intermittent predator-prey model, both the intermittent predation and intrinsic growth rates of the prey and predator species can greatly impact the permanence, extinction, and coexistence of the population.
1. Introduction In real ecosystems, since the spatial distribution and dynamics of a population are greatly affected by their spatial heterogeneity and population mobility, dispersal becomes one of the dominant themes in mathematical biology. In fact, animal dispersal movements between patches are extremely prevalent in ecological environments; for example, many types of birds and mammals will migrate from cold regions to warm regions in search of a better habitat or a breeding site [1]. Therefore, to take spatial heterogeneity into account, realistic population models should contain the dispersal process. During the past couple of decades, predator-prey models with diffusion in a patchy environment have attracted significant attention from ecologists, biologists, and biomathematicians. Many important works and monographs about the properties of population dynamics in a spatial idiosyncratic environment, for example, permanence, extinction, and global asymptotic stability of positive periodic solutions, have been written (see [2β14]). Teng and Chen [6] considered a nonautonomous predator-prey Lotka-Volterra type dispersal system with periodic coefficients and distributed delays: π₯1Μ (π‘) = π₯1 (π‘) [π1 (π‘) β π1 (π‘) π₯1 (π‘)
0
β π (π‘) β«
ββ
π12 (π ) π¦ (π‘ + π ) ππ ]
π
+ βπ1π (π‘) [π₯π (π‘) β π₯1 (π‘)] , π=1
π
π₯πΜ (π‘) = π₯π (π‘) [ππ (π‘) β ππ (π‘) π₯π (π‘)] + β πππ (π‘) [π₯π (π‘) π=1
β π₯π (π‘)] , π = 2, 3, . . . , π, 0
π¦Μ (π‘) = π¦ (π‘) [βπ (π‘) + π (π‘) β«
ββ
0
β π (π‘) β«
ββ
π21 (π ) π₯1 (π‘ + π ) ππ
π21 (π ) π₯1 (π‘ + π ) ππ ] , (1)
where π¦ is the population density of the predator species confined to the 1st patch. Criteria for the permanence, extinction, and existence of positive periodic solutions for system (1) were established. In this model, the prey dispersal behavior occurs at every point in time and simultaneously
2
Discrete Dynamics in Nature and Society
between any patches; that is, it is a continuous bidirectional dispersal. However in practice, it is often the case that diffusion occurs in the form of regular pulses. For example, when winter comes, birds will migrate between patches to seek a better habitat, whereas they do not diffuse in other seasons, and the dispersion of foliage seeds occurs at a fixed period of time every year. For another example, in the Pacific Northwest, Larimichthys polyactis cross over deep water during the winter and migrate to the coast during the spring; then, 3β6 months after spawning, they scatter offshore and return to the depths of the sea during late autumn [15]. All these types of migratory behaviors are appropriately assumed to be in the form of pulses in the modeling process. Thus, impulsive diffusion provides a more natural description. Currently, theories of impulsive differential equations [16] have been introduced into population dynamics. A large number of models have been described by impulsive diffusion (see [14, 17β24]) during the past couple of decades. Shao [23] considered the following predator-prey models with impulsive prey diffusion between two patches:
they are extensively killed by humans, and other carnivorous animals can prolong the journey [25]. As another example, the wild goose will fly to the south when winter comes; in this process, they will stop to rest in some places and at certain time periods. In other words, their diffusion behavior is neither continuous at all times nor impulsive at a fixed time, but it is intermittent within some time intervals. Therefore, it is more reasonable to model this kind of population dynamics with intermittent dispersals. Zhang et al. [26] considered the following nonautonomous almost periodic single species model with intermittent dispersals and dispersal delays between two patches: π₯1Μ (π‘) = π₯1 (π‘) [π1 (π‘) β π1 (π‘) π₯1 (π‘)] , π₯2Μ (π‘) = π₯2 (π‘) [π2 (π‘) β π2 (π‘) π₯2 (π‘)] , π‘ β [π2π , π2π+1 ) , β ), π₯1 (π2π+1 ) = π1 π₯1 (π2π+1 β ), π₯2 (π2π+1 ) = π2 π₯2 (π2π+1
π‘ = π2π+1 ,
π₯1Μ (π‘) = π₯1 (π‘) (π1 β π1 π₯1 (π‘) β π1 π¦ (π‘)) , π1 (π‘) β Μπ1 (π‘) π₯1 (π‘)] π₯1Μ (π‘) = π₯1 (π‘) [Μ
π₯2Μ (π‘) = π₯2 (π‘) (π2 β π2 π₯2 (π‘)) ,
π‘ =ΜΈ ππ, Ξπ₯1 (π‘) = π1 (π₯2 (π‘) β π₯1 (π‘)) ,
(3)
+ π·12 (π‘) (π₯2 (π‘ β π1 ) β π₯1 (π‘)) ,
π¦Μ (π‘) = π¦ (π‘) (βπ3 + π3 π₯1 (π‘ β π1 ) β π2 π¦ (π‘ β π2 )) , (2)
π₯2Μ (π‘) = π₯2 (π‘) [Μ π2 (π‘) β Μπ2 (π‘) π₯2 (π‘)] + π·21 (π‘) (π₯1 (π‘ β π2 ) β π₯2 (π‘)) ,
Ξπ₯2 (π‘) = π2 (π₯1 (π‘) β π₯2 (π‘)) ,
π‘ β [π2π+1 , π2π+2 ) ,
Ξπ¦ (π‘) = 0, π‘ = ππ, where the pulse diffusion of the species π₯ occurs in every period π (a positive constant) and ππ is the dispersal rate in the πth patch satisfying 0 < ππ < 1 for π = 1, 2. The system evolves from its initial state without being further affected by diffusion until the next pulse appears. Ξπ₯π (ππ) = π₯π (ππ+ )βπ₯π (ππ), where π₯π (ππ+ ) represents the density of the population of the prey species π₯ in the πth patch immediately after the πth diffusion pulse at time π‘ = ππ; π₯π (ππ) represents the density of the population of the prey species in the πth patch before the πth diffusion pulse at time π‘ = ππ, π = 0, 1, . . .. Criteria for the global attractivity and permanence of system (2) were obtained. Furthermore, migration movements of the population will be influenced by many uncertain factors (such as the landscape and weather). Therefore, the dispersal movement of migratory species will have to be suspended when the environment becomes unavailable. In other words, patches would permit normal movement patterns between patches to occur only during certain time intervals instead of all the time. For instance, in the Canary Islands of Spain, Anas platyrhynchos undergo a spring migration from early March to the end of March and a fall migration from late September to the end of October, departing as late as early November, during which
β ), π₯1 (π2π+2 ) = π·1 π₯1 (π2π+2 β ), π₯2 (π2π+2 ) = π·2 π₯2 (π2π+2
π‘ = π2π+2 , where all the parameters are almost periodic and the dispersal movement happens only in the time interval π‘ β [π2π+1 , π2π+2 ) but not in π‘ β [π2π , π2π+1 ). Here, π β π. Criteria for the existence, uniqueness, and global attractivity of positive almost periodic solution for system (3) were established. Moreover, in a real ecological system, there always exist natural enemies during the migratory process between patches. For example, annually, at the end of July, with the arrival of the dry season, millions of wildebeests, zebras, and other herbivorous wildlife form a migratory army, migrating from the Serengeti National Park, Tanzania, Africa, to Kenyaβs Masai Mara National Nature Reserve to find enough water and food. Along the way, they will be preyed upon by lions, leopards, and so on. Additionally, crocodiles and hippopotamus will wait and ambush the migration species in the Mara River. As the seasons alternate, i.e., when the rainy season comes, the migration movement starts again, and these species will return to the Serengeti National Park, and vice versa [27]. Obviously, the predation behavior for the above situation is intermittent and only happens in the channels where the dispersals of migratory species occur.
Discrete Dynamics in Nature and Society
3
Motivated by the above consideration, in this paper, we introduce an almost periodic predator-prey model with intermittent predation and discontinuous prey dispersal between two patches: π₯1Μ (π‘) = π₯1 (π‘) [π1 (π‘) β π1 (π‘) π₯1 (π‘)] , π₯2Μ (π‘) = π₯2 (π‘) [π2 (π‘) β π2 (π‘) π₯2 (π‘)] , π¦Μ (π‘) = π¦ (π‘) [π3 (π‘) β π3 (π‘) π¦ (π‘)] , π‘ β [π2π , π2π+1 ) , β ), π₯1 (π2π+1 ) = π1 π₯1 (π2π+1 β ), π₯2 (π2π+1 ) = π2 π₯2 (π2π+1 β ), π¦ (π2π+1 ) = π¦ (π2π+1
π‘ = π2π+1 , π1 (π‘) β Μπ1 (π‘) π₯1 (π‘) β π1 (π‘) π¦ (π‘)] π₯1Μ (π‘) = π₯1 (π‘) [Μ + π·12 (π‘) (π₯2 (π‘ β π1 ) β π₯1 (π‘)) ,
(4)
π₯2Μ (π‘) = π₯2 (π‘) [Μ π2 (π‘) β Μπ2 (π‘) π₯2 (π‘) β π2 (π‘) π¦ (π‘)] + π·21 (π‘) (π₯1 (π‘ β π2 ) β π₯2 (π‘)) , π¦Μ (π‘) = π¦ (π‘) β
[Μ π3 (π‘) β Μπ3 (π‘) π¦ (π‘) + πΜ1 (π‘) π₯1 (π‘) + πΜ2 (π‘) π₯2 (π‘)] , π‘ β [π2π+1 , π2π+2 ) , π₯1 (π2π+2 ) =
β π·1 π₯1 (π2π+2 ),
β ), π₯2 (π2π+2 ) = π·2 π₯2 (π2π+2 β ), π¦ (π2π+2 ) = π¦ (π2π+2
π‘ = π2π+2 , where π₯π (π‘) denotes the prey population density in the πth patch (π = 1, 2) and π¦(π‘) represents the predator population density in channels between two patches. When π‘ β [π2π , π2π+1 ) with π β π, the prey species π₯π inhabits the πth patch and does not disperse. At the same time, the predator species π¦ inhabits the channels between two patches with other food sources. When π‘ β π2π+1 , the intrinsic discipline of the species π₯ in each patch changes. The channels between the two patches will open, and the species π₯ disperses bidirectionally from one patch to another; this dispersal movement will continue for the time interval π‘ β [π2π+1 , π2π+2 ). Meanwhile, the predator species π¦ preys on the species π₯ in the channels. When π‘ β π2π+2 , the gate of the channels will close, the species π₯π will stop dispersing and inhabit patch π. At the same time, the predator species π¦ will also stop preying on species π₯. Obviously, the predation behavior only happens in the time interval π‘ β [π2π+1 , π2π+2 ); that is, it is intermittent. Here, ππ (π‘), πΜπ (π‘) (π = 1, 2) represent the intrinsic growth rates of the species π₯ in the πth patch over
the time intervals [π2π , π2π+1 ) and [π2π+1 , π2π+2 ), respectively; π3 (π‘), πΜ3 (π‘) denote the intrinsic growth rates of the species π¦ over the time intervals [π2π , π2π+1 ) and [π2π+1 , π2π+2 ), respectively; ππ (π‘), Μππ (π‘) (π = 1, 2) represent the intercompetition rates of the species π₯ in the πth patch over the time intervals [π2π , π2π+1 ) and [π2π+1 , π2π+2 ), respectively; π3 (π‘), Μπ3 (π‘) denote the intercompetition rates of the species π¦ over the time intervals [π2π , π2π+1 ) and [π2π+1 , π2π+2 ), respectively; ππ represents the survival rates of switching from stage 1 (without dispersal) to stage 2 (dispersal movement); π·π denotes the survival rates of switching from stage 2 to stage 1; π·ππ (π‘) is the dispersal rates from the πth patch to the πth patch during the time interval [π2π+1 , π2π+2 ); and ππ is the time for the population to disperse from patch π to π (π =ΜΈ π, π, π = 1, 2). During the whole process, the predator species π¦ never disperses. In allusion to system (4) above, our main purpose in this paper is to establish a series of criteria on the ultimate boundedness, permanence, and coexistence of the two populations for system (4). The methods used in this paper are motivated by the works on the permanence and extinction for periodic predator-prey systems in patchy environments given by Teng and Chen in [6] and the works on the survival analysis for a periodic predator-prey model given by Zhang et al. in [22]. This paper is organized as follows. In Section 2, some definitions, assumptions, and useful lemmas are introduced. In Section 3, we state and prove the main results. Finally, special examples and numerical simulations are illustrated to demonstrate our theoretical results in Section 4.
2. Preliminaries Let π
and π
2 denote the set of real numbers and the 2dimensional Euclidean linear space, respectively, and ππ be a time sequence, satisfying ππ < ππ+1 , with ππ β β as π β β. Here, π β π. Define πΜ = max{π1 , π2 }, ππΆ([π‘ β πΜ, π‘], π
2 ) = {π : ([π‘ β πΜ, π‘] β π
2 ) | π = (π1 , π2 ) is continuous everywhere π, π‘], and π(ππ+ ), π(ππβ ) exist with π(ππ+ ) = except at π‘ = ππ β [π‘βΜ π(ππ )}. Define π΅ππΆ = {π β ππΆ((ββ, 0], π
2 ] : π is bounded}; the norm of π is defined by βπβ = supπβ[βΜπ,0] |π(π)|. Let π΅ππΆ+ = {π = (π1 , π2 ) β π΅ππΆ : ππ (π) β₯ 0 for all π β [ββ, 0] and ππ (0+ ) > 0 for π = 1, 2}. In this paper, we assume that all solutions of system (4) satisfy the following initial conditions: π₯π (π , π) = ππ (π ) , π β [βΜ π, 0] , π₯π (0+ , π) = ππ (0+ ) ,
(5)
π β [βΜ π, 0] , where π = (π1 , π2 ) β π΅ππΆ+ . It is not hard to prove that the functional of right of system (4) is continuous and satisfies the local Lipschitz condition with respect to π in the space π
Γ π΅ππΆ. Therefore, by the fundamental theory of the impulsive functional differential equations with finite delays [16, 28, 29], system (4) has a unique solution (π₯1 (π‘, π), π₯2 (π‘, π), π¦(π‘)) satisfying the initial conditions (5). Obviously, the solution
4
Discrete Dynamics in Nature and Society
(π₯1 (π‘, π), π₯2 (π‘, π), π¦(π‘)) is positive in its maximal interval of the existence. Before going into details, we first draw some very useful definitions and lemmas. Definition 1 (see [30]). The set of sequences {ππ π | ππ π = ππ+1 β ππ , π, π β π, ππ β π΅} is said to be uniformly almost periodic if, for arbitrary π > 0, there exists a relatively dense set in π
of π-almost periodic common for all of the sequences; here π΅ = {ππ | ππ β π
, ππ < ππ+1 , π β π, limπβΒ±β ππ = Β±β}. Definition 2 (see [30]). Assume that the following conditions hold:
(π2 ) We let π π (π‘) = π‘ β ππ , π π β1 (π‘) be an inverse function of the function π π (π‘) (π = 1, 2). (π»1 ) Functions ππ (π‘), πΜπ (π‘), ππ (π‘), Μππ (π‘), and π·ππ (π‘) are πalmost periodic and bounded continuous functions for all π‘ β π
and πππ β₯ 0, Μπππ β₯ 0, 0 < ππ , π·π β€ 1 (π = 1, 2, 3, π = 1, 2). ππ (π‘)) > 0, π΄(ππ (π‘)) > 0, π΄(Μππ (π‘)) > (π»2 ) π΄(ππ (π‘)) > 0, π΄(Μ 0 (π = 1, 2). (π»3 ) There exists a constant π > 0, such that, for any π‘ β₯ 0,
(1) The set of sequences {πππ } is almost periodic, π, π β π. (2) For any π > 0 there exists a real number πΏ > 0 such that if the points π‘σΈ , π‘σΈ σΈ belong to the same interval of continuity of π(π‘) and satisfy the inequality |π‘σΈ β π‘σΈ σΈ | < πΏ, then |π(π‘σΈ ) β π(π‘σΈ σΈ )| < π. (3) For any π > 0 there exists a relatively dense set πΆ such that if π β πΆ, then |π(π‘ + π) β π(π‘)| < π for all π‘ β π
satisfying the condition |π‘ β ππ | > π, π β π. The elements of πΆ are called π-almost periods.
β«
π‘+π
π‘
π β€ lim inf π₯π (π‘) β€ lim supπ₯π (π‘) β€ π, π‘ββ
π‘ββ
π = 1, 2,
π β€ lim inf π¦ (π‘) β€ lim supπ¦ (π‘) β€ π. π‘ββ
(6)
π‘ββ
Definition 4 (see [24]). Let π : π
Γ π
2 β π
2 , then π is said to belong to class πΏ if (i) π is continuous in (ππ , ππ+1 ) Γ π
2 for each π₯ β π
2 , π β π, limβ
π (π‘, π¦) = π (ππ+1 , π₯) .
(π‘,π¦)β(ππ+1 ,π₯)
1 π‘+π β« π (π ) ππ , π‘ββ π π‘
(9)
(π = 1, 2, π β π) .
σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨σ΅¨ σ΅¨ σ΅¨σ΅¨ β ln β σ΅¨σ΅¨ ππ σ΅¨σ΅¨σ΅¨ β€ ππ , π β [0, π] . σ΅¨σ΅¨ σ΅¨σ΅¨π‘ 0 (π, π = 1, 2, π β π). (π»6 ) There exist constants π > 0, ππ > 0 (π = 1, 2), such that limπ‘ββ inf π΅π (π‘) = limπ‘ββ [ππ (πππ (π‘)) β Μ β₯ π, where π Μ = max{1, [π]} and [π] ππππ (π π β1 (π‘))/π] is the integer part of π (π β π). (π»7 ) The set of sequences {πππ }, π β π, is uniformly almost periodic, and ππ+2 = ππ + π, inf πβπ|ππ+1 β ππ | > 0, where π belongs to the relatively dense set π (π = 1, 2, 3). (π»8 ) There exists a constant πΜ > 0, such that, for any π‘ β₯ 0, π‘+Μ π
π‘
(π1 ) If π(π‘), π‘ β π
, is an almost periodic function, we define
ln βππ > 0
(π»4 ) There exists a constant ππ > 0 (π = 1, 2), such that, for any π‘ β₯ 0,
β«
In this paper, there are some notations and assumptions that shall be used:
β π‘ 0. π3 (π ) + 2 2
(11)
Now, we give some useful lemmas which will be used in the proofs of the main results. If there is no predator π¦ in system (4), we have the following predator-free system:
π΄ (π) = lim
ππ’ = supπ (π‘) , π‘βπ
π₯πΜ (π‘) = π₯π (π‘) [πππ (π‘) β πππ (π‘) π₯π (π‘)] (8)
+ ππππ (π‘) [π₯π (π‘ β ππ ) β π₯π (π‘)] ,
ππ = inf π (π‘) ,
π‘ β [ππ , ππ+1 ) ,
π‘βπ
where π is a positive constant.
β ) , π‘ = ππ+1 , π₯π (ππ+1 ) = βππ π₯π (ππ+1
(12)
Discrete Dynamics in Nature and Society
5 β 2 β V β π
2 , lim(π‘,π’)β(ππ+1 ) π(π‘, π’) = π(ππ+1 , V) exists, and ππ : π
β 2 π
is nondecreasing. Let π(π‘) be the maximal solution of the following vector impulsive differential system:
where πππ (π‘) =
(β1)π + 1 (β1)π+1 + 1 πΜπ (π‘) , ππ (π‘) + 2 2
πππ (π‘) =
(β1)π + 1 (β1)π+1 + 1 Μ ππ (π‘) + ππ (π‘) , 2 2
ππππ (π‘) = βππ =
π+1
(β1)
2
+1
π’Μ = π (π‘, π (π‘, π₯)) (13)
π·ππ (π‘) ,
π‘ β [ππ , ππ+1 ) ,
π’ (π‘) = ππ (π’ (π‘β )) ,
(18)
π’ (π‘0 ) = π’ (π‘0β ) = π’0 π‘ = ππ+1 ,
π
(β1) + 1 (β1)π+1 + 1 ππ + π·π , 2 2
existing on [π‘0 , β). Then π(π‘0 , π₯0 ) β€ π’0 implies that
π =ΜΈ π, π, π = 1, 2.
π (π‘, π₯ (π‘)) β€ π (π‘) , π‘ β₯ π‘0 ,
For system (12), we have the following result.
(19)
Lemma 5 (see [26] Theorem 3.4). Suppose that assumptions (π»1 )β(π»7 ) hold; then system (12) has a unique globally attractive positive π-almost periodic solution π₯β (π‘) = (π₯1β (π‘), π₯2β (π‘)).
where π₯(π‘) = π₯(π‘, π‘0 , π₯0 ) is any solution of system (12) existing on [π‘0 , β). Here, we state that function π(π‘, π’) is quasimonotone nondecreasing in π’; if π’, V β π
2 , π’ β€ V, and π₯π = π¦π for some 1 β€ π β€ 2, then ππ (π‘, π₯) β€ ππ (π‘, π¦).
If there is no prey π₯ in system (4), we have the following prey-free system:
3. Main Results
π¦Μ (π‘) = π¦ (π‘) [π3π (π‘) β π3π (π‘) π¦ (π‘)] , π‘ β [ππ , ππ+1 ) , (14) β ) , π‘ = ππ+1 , π¦ (ππ+1 ) = π¦ (ππ+1
First, in terms of the ultimate boundedness for system (4), we obtain the following result. Theorem 8. Suppose that assumptions (π»1 )β(π»9 ) hold. Then there is a constant π > 0 such that lim sup π₯π (π‘, π) β€ π
where π3π (π‘) = π3π (π‘) =
(β1)π + 1 (β1)π+1 + 1 πΜ3 (π‘) , π3 (π‘) + 2 2 (β1)π + 1 (β1)π+1 + 1 Μ π3 (π‘) + π3 (π‘) . 2 2
π‘ββ
(20)
lim sup π¦ (π‘) β€ π
(15)
For system (14), we have the following result. Lemma 6 (see [26]). Suppose assumptions (π»1 )-(π»2 ) and (π»7 )-(π»8 ) hold; then system (14) has a unique globally attractive positive π-almost periodic solution π¦β (π‘) β π΄ππΆ. For (π‘, π₯) β (ππ , ππ+1 ) Γ π
2 , we define
π‘ββ
for any positive solutions (π₯1 (π‘, π), π₯2 (π‘, π), π¦(π‘, π)) of system (4). Proof. Let (π₯1 (π‘, π), π₯2 (π‘, π), π¦(π‘)) be any positive solution of system (4) satisfying the initial condition (5). From system (4), we have π₯πΜ (π‘) β€ π₯π (π‘) [πππ (π‘) β πππ (π‘) π₯π (π‘)] + ππππ (π‘) [π₯π (π‘ β ππ ) β π₯π (π‘)] ,
π·+ π (π‘, π₯)
π‘ β [ππ , ππ+1 ) ,
1 = lim sup [π (π‘ + β, π₯ + βπ (π‘, π₯)) β π (π‘, π₯)] , + β ββ0
(16)
where π = (π1 , π2 ) is the right-hand side of system (12). We give the following vector comparison results of the impulsive differential equations. Lemma 7 (see [1]). Let π : π
+ Γ
(π = 1, 2) ,
π
+2
+
and π β πΏ. Assume that
π· π (π‘, π₯) β€ π (π‘, π (π‘, π₯))
π‘ β [ππ , ππ+1 ) ,
π (π‘, π₯ (π‘)) β€ ππ (π (π‘, π₯)) ,
π‘ = ππ+1 ,
(17)
where π(π‘, π’) : π
Γ π
2 β π
2 is continuous in (ππ , ππ+1 ) Γ π
2 and quasi-monotone nondecreasing in π’, for
(21)
β ) , π‘ = ππ+1 , π₯π (ππ+1 ) = βππ π₯π (ππ+1
for all π‘ β₯ 0, π, π = 1, 2 (π =ΜΈ π). Then by the Lemmas 5 and 7, we obtain π₯π (π‘) β€ π₯π (π‘) , βπ‘ β₯ 0, π = 1, 2,
(22)
where (π₯1 (π‘), π₯2 (π‘)) is the solution of system (12) with the initial condition (5). Under assumptions (π»1 )β(π»7 ), from Lemma 5 we obtain π₯π (π‘) β π₯β (π‘) (π = 1, 2, π‘ β β), where π₯β (π‘) = (π₯1β (π‘), π₯2β (π‘)) is the globally asymptotically stable positive π-almost periodic solution of system (12). Hence, π₯π (π‘) (π = 1, 2) is bounded on π
+ . Next, we choose a constant π1 = supπ‘βπ
|π₯β (π‘)|, where |π₯β (π‘)| = β2π=1 π₯πβ (π‘). Since π₯π (π‘) β
6
Discrete Dynamics in Nature and Society
π₯β (π‘) (π = 1, 2, π‘ β β), there is a pair of π1 > 0 and π > 0 such that π₯π (π‘) β€ π₯πβ (π‘) + π β π1
βπ‘ β₯ π1 , π = 1, 2.
βπ‘ β₯ π1 , π = 1, 2.
lim sup π¦ (π‘) β€ π2 .
Therefore, (26) holds. Choose a constant π = max{π1 , π2 }; then we can see (24)
lim sup π₯π (π‘, π) β€ π
π‘ββ
Consequently, lim sup π₯π (π‘, π) β€ π1 , π = 1, 2.
π‘ββ
lim sup π¦ (π‘) β€ π2 .
(26)
From systems (4) and (24) we obtain
π‘ β [ππ , ππ+1 ) ,
(27)
β ) , π‘ = ππ+1 , π¦ (ππ+1 ) = π¦ (ππ+1
πΜ2π (π‘) =
(β1)
+1
2
(28)
π2 (π‘) ,
From the comparison theorem of the impulsive differential equations [16, 28, 29], we have π¦(π‘) β€ π§(π‘) for all π‘ β₯ π1 , where π§(π‘) is the solution of the following auxiliary equation: π§Μ (π‘) = π§ (π‘) β
[π3π (π‘) + πΜ1π (π‘) π1 + πΜ2π (π‘) π2 β π3π (π‘) π§ (π‘)] , π‘ β [ππ , ππ+1 ) ,
with initial value π§(π1 ) = π¦(π1 ). According to assumptions (π»1 )-(π»2 ) and (π»7 )-(π»8 ), we also have π‘
(β1)π+1 + 1 (β1)π + 1 { π3 (π ) + [Μ π3 (π ) 2 2
(30)
+ πΜ1π (π‘) π1 + πΜ2π (π‘) π1 ]} ππ > 0. Hence from Lemma 6, system (29) has a unique globally attractive positive π-almost periodic solution π§β (π‘). For any constant π > 0, there is a π2 β₯ π1 such that π§(π‘) < π§β (π‘) + π for all π‘ β₯ π2 . Therefore, we have π¦ (π‘) < π§β (π‘) + π β€
Proof. Let (π₯1 (π‘, π), π₯2 (π‘, π), π¦(π‘)) be any positive solutions of system (4) satisfying the initial condition (5). By system (4), we easily obtain
1 max π§β (π‘) + π β π2 , 2 π‘β[0,π]
(31)
π‘ = ππ+1 .
Using comparison theorem of the impulsive differential equations [16, 28, 29] and Lemma 6, we can easily obtain π¦(π‘) β₯ π¦(π‘), for all π‘ β π
, where π¦(π‘) is the solution of system (14) with condition (5). Under the assumptions (π»1 )-(π»2 ) and (π»7 )-(π»8 ), by Lemma 6, we can obtain π¦(π‘) β π¦β (π‘) as π‘ β β, where π¦β (π‘) is the globally asymptotically stable positive π-almost periodic solution of system (14). Hence, we can easily obtain that, for π0 = (1/2)minπ‘β[0,π] π¦β (π‘), there exists a π1 > 0 such that π¦ (π‘) > π¦β (π‘) β π0 β₯
(29)
β ) , π‘ = ππ+1 , π§ (ππ+1 ) = π§ (ππ+1
β«
Next, on the permanence of system (4), we have the following results.
β ), π¦ (ππ+1 ) = π¦ (ππ+1
π β π.
π‘+π
this completes the proof of Theorem 8.
π‘ β [ππ , ππ+1 ) , (34)
(β1)π+1 + 1 π1 (π‘) , 2 π+1
π‘ββ
π¦Μ (π‘) β₯ π¦ (π‘) [π3π (π‘) β π3π (π‘) π¦ (π‘)] ,
for all π‘ β₯ π1 , where πΜ1π (π‘) =
(33)
Theorem 9. Suppose that assumptions (π»1 )-(π»2 ) and (π»7 )(π»8 ) hold. Then the predator species π¦ of system (4) is permanent.
π¦Μ (π‘) β€ π¦ (π‘) β
[π3π (π‘) + πΜ1π (π‘) π1 + πΜ2π (π‘) π2 β π3π (π‘) π¦ (π‘)] ,
(π = 1, 2) ,
lim sup π¦ (π‘) β€ π;
(25)
Further, we prove that there is a constant π2 > 0 such that π‘ββ
(32)
π‘ββ
(23)
Hence, by (22) we have π₯π (π‘) β€ π1
for all π‘ β₯ π2 . Consequently,
1 min π¦β (π‘) β π2 , 2 π‘β[0,π]
(35)
for all π‘ β₯ π1 . Together with Theorem 8, we have that the predator species π¦ of system (4) is permanent. This completes the proof of Theorem 9. Remark 10. Set π=β«
π‘+π
π‘
[
(β1)π+1 + 1 (β1)π + 1 πΜ3 (π )] ππ . π3 (π ) + 2 2
(36)
In system (4), based on the assumptions and the actual biological meanings of the parameters π3 (π‘), π3 (π‘) and πΜ3 (π‘), Μπ3 (π‘), we can see that the conditions in Theorem 9 are easily satisfied. Additionally, the constant π represents the minimal total growth rate of the predator species π¦. If π > 0, which means that the predator species π¦ has another food resource. As a result, Theorem 9 implies that if π > 0, then the predator species π¦ will be permanent regardless of whether the prey species π₯ exists or not.
Discrete Dynamics in Nature and Society
7
Next, on the permanence of the prey species π₯ of system (4), we have the following result. Theorem 11. Suppose that assumptions (π»1 )β(π»7 ) and the following inequality π‘+π
β«
π‘
(β1)π + 1 ππ (π ) 2 π+1
+
(β1) 2
+
β
π‘β€ππ+1 β€π‘+π
(Μ ππ (π ) β ππ (π ) π¦β (π )) ππ
(37)
π‘
+
β π‘ 0.
ππ (ππ+1 ) =
β βππ ππ (ππ+1 ),
(39)
π‘ = ππ+1 ,
σ΅¨σ΅¨ σ΅¨ π σ΅¨σ΅¨π’πΎ (π‘, π‘0 , π’0 ) β π’πΎβ (π‘)σ΅¨σ΅¨σ΅¨ < 1 σ΅¨ σ΅¨ 2
(β1)π+1 + 1 π1 (π‘) , 2
(β1)π+1 + 1 π2π (π‘) = π2 (π‘) . 2
for all π‘ β₯ π‘0 + π‘3 , where π’πΎ (π‘, π‘0 , π’0 ) is the solution of system (42) with the initial condition π’πΎ (π‘0 ) = π’0 . According to the theorem of the continuity of solutions with respect to parameters of the impulsive differential equations [17], there exist πΎ0 β (0, π1 ) and πΎ0 < ππβ1 (π‘) β πΎ0 such that
π‘+π
π‘
+ +
β π‘ 0, solution π’πΎβ (π‘). For above π1 > 0 and π 2 Μ > 0, such that, for any π‘0 β₯ 0 and there is a π‘3 = π‘(π1 , π) Μβ1 , π], Μ π’0 β [π
Consider the following auxiliary system:
= ππ (π‘) [πππ (π‘) β πππ (π‘) ππ (π‘) β πππ (π‘) (π¦β (π ) + π1 )]
(β1)π + 1 π’ (π‘) [π3 (π‘) β π3 (π‘) π’ (π‘)] 2
π‘ β [ππ , ππ+1 ) ,
ln βππ > 0
Proof. Owing to condition (37), we have that there is a small enough constant π1 > 0, such that β«
π’Μ (π‘) = +
for (π = 1, 2, π β π) hold. Then the prey species π₯ of system (4) is permanent.
π‘+π
And then, based on assumption π»8 and Lemma 5, for arbitrary constant πΎ > 0 the following system
π¦Μ (π‘) β€ (41)
ln βππ > 0.
Then, from (41) and Lemma 5, we know system (39) has a unique globally attractive positive π-almost periodic solution β β (π‘), π2π (π‘)). ππβ1 (π‘) = (π1π 1 1
+
(β1)π + 1 π¦ (π‘) [π3 (π‘) β π3 (π‘) π¦ (π‘)] 2
(β1)π+1 + 1 π¦ (π‘) 2
β
[Μ π3 (π‘) + πΜ1 (π‘) πΎ0 + πΜ2 (π‘) πΎ0 β Μπ3 (π‘) π¦ (π‘)] , π‘ β [ππ , ππ+1 ) , β ), π¦ (ππ+1 ) = π¦ (ππ+1
π‘ = ππ+1 ,
(47)
8
Discrete Dynamics in Nature and Society
for all π‘ β₯ π‘2 . By comparison theorem of impulsive differential equations [16, 28, 29], we also have π¦(π‘) β€ π’πΎ0 (π‘) for all π‘ β₯ π‘2 , where π’πΎ0 (π‘) is the solution of system (42) satisfying πΎ = πΎ0 and with the initial value π’πΎ0 (π‘2 ) = π¦(π‘2 ). For (43), let π‘0 = π‘2 and π’(π‘0 ) = π¦(π‘2 ). From πβ1 β€ π¦(π‘) β€ π, we have for πΎ = πΎ0 1 π’πΎ0 (π‘) = π’πΎ0 (π‘, π‘2 , π¦ (π‘2 )) < π’πΎβ0 (π‘) + π1 , 2
πΎ0 exp {β [(π0 + π0 ) π‘3 + 2ππ ]} > π₯π (π‘6 ) = βππ π₯π (π‘6β ) β₯ πΎ0 exp (ln βππ ) β₯ πΎ0 exp (βππ )
(49)
which leads to a contradiction. For Case 2, we have π₯π (π‘) < πΎ0 , for all π‘ β (π‘6 , π‘7 ]. Assume π‘7 β π‘6 β€ π‘3 , then for any π‘ β [π‘6 , π‘7 ) we have π₯πΜ (π‘) β₯ π₯π (π‘) [πππ (π‘) β πππ (π‘) π1 β πππ (π‘) π]
Taking into account system (4) again, from above inequality (49) and system (4), we have
β ππππ (π‘) π₯π (π‘)
π₯πΜ (π‘) β₯ π₯π (π‘) [πππ (π‘) β πππ (π‘) π₯π (π‘) β πππ (π‘) (π¦β (π‘) + π1 )] + ππππ (π‘) [π₯π (π‘ β ππ ) β π₯π (π‘)] , π‘ β [ππ , ππ+1 ) , π₯π (ππ+1 ) =
β βππ π₯π (ππ+1 ),
β₯ π₯π (π‘) [πππ (π‘) β πππ (π‘) π1 β πππ (π‘) π β ππππ (π‘)]
β ), π₯π (ππ+1 ) = βππ π₯π (ππ+1
π‘ = ππ+1 ,
where
for all π‘ β₯ π‘2 + π‘3 . From Lemma 7, we also have π₯π (π‘) β₯ πππ1 (π‘) (π = 1, 2) for all π‘ β₯ π‘2 + π‘3 , where ππ1 (π‘) = (π1π1 (π‘), π2π1 (π‘)) is the solution of system (39) satisfying the initial value πππ1 (π‘2 + π‘3 ) = π₯π (π‘2 + π‘3 ). Because of the globally attractivity of positive π-almost periodic solution ππβ1 (π‘) of system (39), we further have that, for previous constant πΎ0 , there exists a π‘4 β₯ π‘2 + π‘3 such that
π0 (56) σ΅¨ σ΅¨ = min {σ΅¨σ΅¨σ΅¨σ΅¨πππ (π‘) β πππ (π‘) π1 β πππ (π‘) π β ππππ (π‘)σ΅¨σ΅¨σ΅¨σ΅¨} . π‘β[0,π] Integrating (55) from π‘6 to π‘, for any π‘ β [π‘6 , π‘7 ], we obtain π₯π (π‘) β₯ π₯π (π‘6 ) π‘
ππ1 (π‘) β₯ ππβ1 (π‘) β πΎ0 > πΎ0 ,
(51)
for all π‘ β₯ π‘4 . Hence, we have π₯π (π‘) > πΎ0 for all π‘ β₯ π‘4 , which leads to a contradiction. Hence, there exists a constant π‘5 β₯ π‘2 such that π₯π (π‘5 ) > πΎ0 . We now prove that βπ‘ β₯ π‘5 ,
(52)
β
exp {β« [πππ (π ) β πππ (π ) π1 β πππ (π ) π β ππππ (π )] ππ π‘6
+
β
π‘6 πΎ0 > πΎ0 exp{β[(π0 + π0 )π‘3 + 2ππ ]}, there are constants π‘6 and π‘7 satisfying π‘7 β₯ π‘6 > π‘5 , such that π₯π (π‘7 ) < πΎ0 exp{β[(π0 + π0 )π‘3 + 2ππ ]}, π₯π (π‘6β ) β₯ πΎ0 , π₯π (π‘6 ) β€ πΎ0 . There are two cases for π‘7 .
Since π₯π (π‘) β€ πΎ0 for all π‘ β (π‘6 , π‘7 ]. Then we have π¦Μ (π‘) β€ +
(β1)π + 1 π¦ (π‘) [π3 (π‘) β π3 (π‘) π¦ (π‘)] 2
(β1)π+1 + 1 π¦ (π‘) 2
β
[Μ π3 (π‘) + πΜ1 (π‘) πΎ0 + πΜ2 (π‘) πΎ0 β Μπ3 (π‘) π¦ (π‘)] , π‘ β [ππ , ππ+1 ) ,
Case 1. π‘6 = π‘7 . Case 2. π‘6 < π‘7 .
(57)
Assume π‘7 β π‘6 > π‘3 ; then when π‘ β [π‘6 , π‘6 + π‘3 ], according to above discussion we directly have
where
σ΅¨ + ππππ (π‘)σ΅¨σ΅¨σ΅¨σ΅¨} ,
(55)
β₯ π₯π (π‘) π0 , π‘ β [ππ , ππ+1 ) ,
(50)
π‘ = ππ+1 ,
π₯π (π‘) β₯ πΎ0 exp {β [(π0 + π0 ) π‘3 + 2ππ ]}
(54)
β₯ πΎ0 exp {β [(π0 + π0 ) π‘3 + 2ππ ]} ,
(48)
for all π‘ β₯ π‘2 + π‘3 . Then π¦(π‘) < π’πΎβ0 (π‘) + (1/2)π1 for all π‘ β₯ π‘2 + π‘3 . Thus, by condition (45), we further get π¦ (π‘) < π¦β (π‘) + π1 .
For Case 1, we see that π‘6 is an impulsive time. Hence, there is an integer π > 0 such that π‘6 = π‘π . We have
β ), π¦ (ππ+1 ) = π¦ (ππ+1
π‘ = ππ+1 ,
(60)
Discrete Dynamics in Nature and Society
9
for π‘ β (π‘6 , π‘7 ]. From comparison theorem of impulsive differential equations [16, 28, 29], we also have π¦(π‘) β©½ π’πΎ0 (π‘) for all π‘ β (π‘6 , π‘7 ), where π’πΎ0 (π‘) is the solution of system (42) with πΎ = πΎ0 and satisfies π’πΎ0 (π‘6 ) = π¦(π‘6 ). In (43), choose a π‘0 = π‘6 and π’0 = π¦(π‘6 ), since πβ1 β€ π¦(π‘6 ) β€ π, we have for πΎ = πΎ0 1 π’πΎ0 (π‘) = π’πΎ0 (π‘, π‘6 , π¦ (π‘6 )) β€ π’πΎβ0 (π‘) + π1 2
Hence, π1 2
βπ‘ β₯ π‘6 + π‘3 .
(62)
Thus, from (45) we further have π¦ (π‘) < π¦β (π‘) + π1
βπ‘ β₯ π‘6 + π‘3 .
(63)
From (63), we have π₯πΜ (π‘) β₯ π₯π (π‘) [πππ (π‘) β πππ (π‘) πΎ0 β πππ (π‘) (π¦β (π‘) + π1 )] β ππππ (π‘) π₯π (π‘)
π‘ β [ππ , ππ+1 ) ,
(64)
β ) , π‘ = ππ+1 , π₯π (ππ+1 ) = βππ π₯π (ππ+1
For arbitrary π‘ β [π‘6 + π‘3 , π‘7 ). We choose an integer π β₯ 0 such that π‘7 β [π‘6 + π‘3 + ππ, π‘6 + π‘3 + (π + 1)π); then integrating (64) from π‘6 + π‘3 to π‘7 , we have πΎ0 exp {β [(π0 + π0 ) π‘3 + 2ππ ]} > π₯π (π‘7 ) = π₯π (π‘6 + π‘3 ) π‘7
β
exp {β«
π‘6 +π‘3
[πππ (π ) β πππ (π ) πΎ0
β πππ (π ) (π¦β (π ) + π1 ) β ππππ (π )] ππ +
β
π‘6 +π‘3 0,
β 0 0.
β 0