Coexistence for an Almost Periodic Predator-Prey Model with ...

4 downloads 0 Views 3MB Size Report
Oct 23, 2017 - College of Mathematics and System Science, Xinjiang University, Urumqi 830046, .... Furthermore, migration movements of the population.
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