Discrete Dynamics in Nature and Society
Volume 2009 (2009), Article ID 256195, 15 pages
doi:10.1155/2009/256195
Research Article

A Predator-Prey Gompertz Model with Time Delay and Impulsive Perturbations on the Prey

School of Mathematics and Computer Science, Shanxi Normal University, Shanxi, Linfen 041004, China

Received 28 June 2009; Accepted 5 November 2009

Academic Editor: Leonid Berezansky

Copyright © 2009 Jianwen Jia and Chunhua Li. 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.

Abstract

We introduce and study a Gompertz model with time delay and impulsive perturbations on the prey. By using the discrete dynamical system determined by the stroboscopic map, we obtain the sufficient conditions for the existence and global attractivity of the “predator-extinction” periodic solution. With the theory on the delay functional and impulsive differential equation, we obtain the appropriate condition for the permanence of the system.

1. Introduction

It is well known that Gompertz equation [1] describes the growth law for a single species. The model reads as

𝑑 𝑥 ( 𝑡 ) 𝐾 𝑑 𝑡 = 𝑟 𝑥 ( 𝑡 ) l n , 𝑥 ( 𝑡 ) ( 1 . 1 ) where 𝑥 ( 𝑡 ) is the density of the population, 𝑟 is a positive constant called the intrinsic growth rate, the positive constant 𝐾 is usually referred to as the environment carrying capacity or saturation level, and 𝑟 l n ( 𝐾 / 𝑥 ) denotes relative growth rate.

Many systems in physics, chemistry, biology, and information science have impulsive dynamical behavior due to abrupt jumps at certain instants during the evolving processes. This dynamical behavior can be modeled by impulsive differential equations. The theory of impulsive differential systems has been developed by numerous mathematicians [24]. But more realistic models should include some of the past states of these systems, that is, ideally, a real system should be modeled by differential equations with time delays. Recently, the investigation of impulsive delay differential equations is beginning [57].

In this paper, we need to consider a function

𝑥 𝜙 ( 𝑥 ) = 2 ( 𝑡 ) 1 + 𝑚 𝑥 2 , ( 𝑡 ) ( 1 . 2 ) which is called Holling III functional response. The coefficient 𝑚 is a positive constant and its biological meaning can be found in [8].

We introduce the model as follows:

𝑑 𝑥 ( 𝑡 ) 𝐾 𝑑 𝑡 = 𝑟 𝑥 ( 𝑡 ) l n 𝑥 𝑥 ( 𝑡 ) 2 ( 𝑡 ) 1 + 𝑚 𝑥 2 ( 𝑡 ) 𝑦 ( 𝑡 ) , 𝑑 𝑦 ( 𝑡 ) 𝑥 𝑑 𝑡 = 𝛽 2 ( 𝑡 𝜏 ) 1 + 𝑚 𝑥 2 ( 𝑡 𝜏 ) 𝑦 ( 𝑡 𝜏 ) 𝑑 𝑦 ( 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , Δ 𝑥 ( 𝑡 ) = 𝑝 𝑥 ( 𝑡 ) , Δ 𝑦 ( 𝑡 ) = 0 , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 , ( 1 . 3 ) where 𝑥 ( 𝑡 ) and 𝑦 ( 𝑡 ) are densities of the prey and the predator, respectively, 𝑟 is the Gomportz intrinsic growth rate of the prey in the absence of the predator, 𝛽 is the conversion rate, 𝑑 is the death rate of the predator, 𝜏 is mean length of the digest period, and 𝑝 ( 0 < 𝑝 < 1 ) represents impulsive harvest to preys by catching or pesticides at 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 .

The initial conditions for system (1.3) are

𝜙 1 ( 𝑠 ) , 𝜙 2 ( 𝑠 ) 𝐶 + [ ] = 𝐶 𝜏 , 0 , 𝑅 2 + , 𝜙 𝑖 ( 0 ) > 0 ( 𝑖 = 1 , 2 ) . ( 1 . 4 )

From the biological point of view, we only consider system (1.3) in the biological meaning region 𝐷 = { ( 𝑥 , 𝑦 ) 𝑥 , 𝑦 0 } .

All outline of this paper is as follows. We give some basic knowledge in Section 2. In Section 3, using discrete dynamical system determined by the stroboscopic map, we obtain the existence and global attractivity of the “predator-extinction” periodic solution. In Section 4, with the theory of delay and impulsive different equations, we obtain the sufficient condition for the permanence of the system. In the last section, we give the numerical simulation and discussion.

2. Basic Knowledge

Let 𝑅 + = [ 0 , ) , 𝑅 2 + = { 𝑥 𝑅 2 | 𝑥 0 } , Ω = i n t 𝑅 2 + , and 𝑁 be the set of all nonnegative integers. Denote that 𝑓 = ( 𝑓 1 , 𝑓 2 ) 𝑇 is the map defined by the right-hand side of the first two equations of system (1.3). Let 𝑉 0 = { 𝑉 𝑅 + × 𝑅 2 + 𝑅 + , continuous on ( 𝑛 𝑇 , ( 𝑛 + 1 ) 𝑇 ] × 𝑅 2 + , and l i m ( 𝑡 , 𝑧 ) ( 𝑛 𝑇 + , 𝑥 ) 𝑉 ( 𝑡 , 𝑧 ) = 𝑉 ( 𝑛 𝑇 + , 𝑥 ) exist } .

Definition 2.1. Let 𝑉 𝑉 0 , then for ( 𝑡 , 𝑋 ) ( 𝑛 𝑇 , ( 𝑛 + 1 ) 𝑇 ] × 𝑅 2 + , the upper right derivative of 𝑉 ( 𝑡 , 𝑥 ) with respect to the impulsive differential system (1.3) is defined as 𝐷 + 𝑉 ( 𝑡 , 𝑋 ) = l i m s u p 0 + 1 [ ] . 𝑉 ( 𝑡 + , 𝑋 + 𝑓 ( 𝑡 , 𝑥 ) ) 𝑉 ( 𝑡 , 𝑋 ) ( 2 . 1 ) The solution of system (1.3) is a piecewise continuous function. The smoothness properties of 𝑓 guarantee the global existence and uniqueness of the solution of system (1.3) [9].

Lemma 2.2 (see [9]). Considering the following impulsive differential inequalities: 𝑑 𝑚 ( 𝑡 ) 𝑑 𝑡 𝑝 ( 𝑡 ) 𝑚 ( 𝑡 ) + 𝑞 ( 𝑡 ) , 𝑡 𝑡 𝑘 𝑚 𝑤 , 𝑘 = 1 , 2 , + 𝑘 𝑑 𝑘 𝑚 𝑡 𝑘 + 𝑏 𝑘 , 𝑡 = 𝑡 𝑘 , 𝑡 𝑡 0 , ( 2 . 2 ) where 𝑝 , 𝑞 𝑃 𝐶 [ 𝑅 + , 𝑅 ] and 𝑑 𝑘 0 , 𝑏 𝑘 are constants, then 𝑚 𝑡 ( 𝑡 ) 𝑚 0 𝑡 0 < 𝑡 𝑘 < 𝑡 𝑑 𝑘 e x p 𝑡 𝑡 0 𝑝 + ( 𝑠 ) 𝑑 𝑠 𝑡 0 < 𝑡 𝑘 < 𝑡 𝑡 𝑘 < 𝑡 𝑗 < 𝑡 𝑑 𝑗 e x p 𝑡 𝑡 0 𝑝 𝑏 ( 𝑠 ) 𝑑 𝑠 𝑘 + 𝑡 𝑡 0 𝑠 < 𝑡 𝑘 < 𝑡 𝑑 𝑘 e x p 𝑡 s 𝑝 ( 𝜎 ) 𝑑 𝜎 𝑞 ( 𝑠 ) 𝑑 𝑠 , 𝑡 𝑡 0 . ( 2 . 3 )

Lemma 2.3. There exists a constant 𝐿 = 𝛽 𝑟 𝐾 ( 𝐾 𝑒 𝑑 / 𝑟 1 ) / 𝑑 > 0 such that 𝑥 ( 𝑡 ) 𝐾 and 𝑦 ( 𝑡 ) 𝐿 for each solution of system (1.3) with 𝑡 large enough.

Proof. Define 𝑉 ( 𝑡 ) = 𝛽 𝑥 ( 𝑡 ) + 𝑦 ( 𝑡 + 𝜏 ) , ( 2 . 4 ) then 𝑉 𝑉 0 . Since 𝑑 𝑥 / 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n ( 𝐾 / 𝑥 ( 𝑡 ) ) , 𝑑 𝑥 / 𝑑 𝑡 𝑥 = 𝐾 = 0 ; in addition, 𝑥 ( 𝑛 𝑇 + ) 𝑥 ( 𝑛 𝑇 ) . Thus 𝑥 ( 𝑡 ) 𝐾 for 𝑡 large enough.
We calculate the upper right derivative of 𝑉 ( 𝑡 ) along a solution of system (1.3) as 𝐷 + 𝑉 ( 𝑡 ) + 𝑑 𝑉 ( 𝑡 ) 𝛽 𝑥 ( 𝑡 ) 𝑟 l n 𝐾 𝑒 𝑑 / 𝑟 𝑥 ( 𝑡 ) = 𝛽 𝑥 ( 𝑡 ) 𝑟 l n 𝐾 𝑒 𝑑 / 𝑟 𝛽 𝑟 𝐾 𝐾 𝑒 𝑑 / 𝑟 . 1 ( 2 . 5 )
So we have that 𝐷 + 𝑉 𝑡 𝑉 ( 𝑡 ) 𝑑 𝑉 ( 𝑡 ) + 𝑀 , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , + 𝑉 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 . ( 2 . 6 ) By Lemma 2.2, for 𝑡 ( 𝑛 𝑇 , ( 𝑛 + 1 ) 𝑇 ] , we have that 𝑉 0 𝑉 ( 𝑡 ) = + 𝛽 𝑟 𝐾 𝐾 𝑒 𝑑 / 𝑟 1 𝑑 e x p ( 𝑑 𝑡 ) + 𝛽 𝑟 𝐾 𝐾 𝑒 𝑑 / 𝑟 1 𝑑 , l i m 𝑡 𝑉 ( 𝑡 ) 𝛽 𝑟 𝐾 𝐾 𝑒 𝑑 / 𝑟 1 𝑑 . ( 2 . 7 ) So, 𝑉 ( 𝑡 ) is ultimately bounded. Therefore by the definition of 𝑉 ( 𝑡 ) , there exists a constant 𝐿 = 𝛽 𝑟 𝐾 ( 𝐾 𝑒 𝑑 / 𝑟 1 ) / 𝑑 > 0 such that 𝑥 ( 𝑡 ) 𝐾 and 𝑦 ( 𝑡 ) 𝐿 for each solution of system (1.3) with 𝑡 large enough. This completes the proof.

Lemma 2.4 (see [10]). Consider a delay equation ̇ 𝑥 ( 𝑡 ) = 𝑎 𝑥 ( 𝑡 𝜏 ) 𝑏 𝑥 ( 𝑡 ) , ( 2 . 8 ) where 𝑎 , 𝑏 , 𝜏 are all positive constants and 𝑥 ( 𝑡 ) > 0 for 𝜏 𝑡 0 . Thus one has the following. (1)If 𝑎 > 𝑏 , then l i m 𝑡 𝑥 ( 𝑡 ) = + .(2)If 𝑎 < 𝑏 , then l i m 𝑡 𝑥 ( 𝑡 ) = 0 .

3. “Predator-Extinction” Periodic Solution

3.1. Existence of a “Predator-Extinction” Periodic Solution

First, we begin analyzing the existence of a “predator-extinction” solution, in which predator is absent from the system, that is,

𝑦 ( 𝑡 ) = 0 , 𝑡 0 . ( 3 . 1 )

In this condition, we know the growth of the prey in the time-interval 𝑛 𝑇 𝑡 ( 𝑛 + 1 ) 𝑇 and give some basic properties of the following subsystem of (1.3):

𝑑 𝑥 ( 𝑡 ) 𝐾 𝑑 𝑡 = 𝑟 𝑥 ( 𝑡 ) l n 𝑥 ( 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , Δ 𝑥 ( 𝑡 ) = 𝑝 𝑥 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 . ( 3 . 2 )

For system (3.2), we let 𝜛 ( 𝑡 ) = l n 𝑥 ( 𝑡 ) and obtain a linear nonhomogeneous impulsive equation

𝑑 𝜛 ( 𝑡 ) 𝑑 𝑡 = 𝑟 𝜛 ( 𝑡 ) + 𝑟 l n 𝐾 , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , Δ 𝜛 ( 𝑡 ) = l n ( 1 𝑝 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 . ( 3 . 3 )

Solving the first equation of system (3.3) between pulses yields

𝜛 ( 𝑡 ) = l n 𝐾 + ( 𝜛 ( ( 𝑛 1 ) 𝑇 ) l n 𝐾 ) 𝑒 𝑟 ( 𝑡 ( 𝑛 1 ) 𝑇 ) , ( 𝑛 1 ) 𝑇 < 𝑡 𝑛 𝑇 . ( 3 . 4 )

For the second equation of system (3.3), using discrete dynamical system determined by the stroboscopic map yields

𝜛 ( 𝑛 𝑇 ) = l n 𝐾 + ( 𝜛 ( ( 𝑛 1 ) 𝑇 ) l n 𝐾 ) 𝑒 𝑟 𝑇 + l n ( 1 𝑝 ) = 𝑓 ( 𝜛 ( 𝑛 1 ) 𝑇 ) , ( 3 . 5 ) where 𝑓 ( 𝜛 ) = l n 𝐾 + ( 𝜛 l n 𝐾 ) 𝑒 𝑟 𝑇 + l n ( 1 𝑝 ) . From (3.5), we can see that this difference system has an equilibrium 𝜛 0 = l n 𝐾 + l n ( 1 𝑝 ) / ( 1 𝑒 𝑟 𝑇 ) , which implies that system (3.3) has a unique T-periodic solution

𝜛 𝜛 ( 𝑡 ) = l n 𝐾 + 0 𝑒 l n 𝐾 𝑟 ( 𝑡 ( 𝑛 1 ) 𝑇 ) , ( 𝑛 1 ) 𝑇 < 𝑡 𝑛 𝑇 . ( 3 . 6 )

Since 𝑥 0 = 𝑒 𝜛 0 , then 𝑥 ( 𝑡 ) = 𝑒 𝜛 ( 𝑡 ) is the unique positive T-periodic solution of (3.2).

In the following, we will prove that 𝑥 ( 𝑡 ) is globally asymptotically stable. By Lemma 2.3, we find that any solution of (3.2) is ultimately upper bounded, so we need only to prove that

l i m 𝑡 | | 𝜛 ( 𝑡 ) 𝜛 | | ( 𝑡 ) = 0 , ( 3 . 7 ) where 𝜛 ( 𝑡 ) is the periodic solution of system (3.3) Suppose that with l n 𝑥 ( 0 ) = l n 𝑥 0 = 𝜛 ( 0 ) = 𝜛 0 .

| | 𝜛 ( 𝑡 ) 𝜛 | | = | | ( 𝑡 ) 𝜛 ( 0 ) 𝜛 0 | | 𝑒 𝑟 𝑡 𝑒 ( 𝑛 1 ) 𝑟 𝑇 , ( 3 . 8 ) since l i m 𝑡 𝑒 ( 𝑛 1 ) 𝑟 𝑇 = 0 ; thus, l i m 𝑡 | 𝜛 ( 𝑡 ) 𝜛 ( 𝑡 ) | = 0 .

Theorem 3.1. System (3.2) has unique a positive periodic solution 𝑥 ( 𝑡 ) which is globally asymptotically stable. That is, system (1.3) has a “predator-extinction” periodic solution ( 𝑥 ( 𝑡 ) , 0 ) for 𝑡 ( 𝑛 𝑇 , ( 𝑛 + 1 ) 𝑇 ] , 𝑛 𝑁 .

3.2. Global Attractivity of the “Predator-Extinction” Periodic Solution

Denote that

1 = 𝛽 𝐾 2 𝑑 𝑚 𝐾 2 . + e x p ( 2 l n ( 1 / 1 𝑝 ) ) / ( e x p ( 𝑟 𝑇 ) 1 ) ( 3 . 9 )

Theorem 3.2. If 1 < 1 , then the “predator-extinction” periodic solution ( 𝑥 ( 𝑡 ) , 0 ) of system (1.3) is globally attractive.

Proof. We denote that ( 𝑥 ( 𝑡 ) , 𝑦 ( 𝑡 ) ) be any solution of system (1.3) with initial condition (1.4). From the second equation of system (1.3), we have that 𝑑 𝑦 ( 𝑡 ) 𝛽 𝑑 𝑡 𝑚 𝑦 ( 𝑡 𝜏 ) 𝑑 𝑦 ( 𝑡 ) . ( 3 . 1 0 )
We consider the following comparison equation: 𝑑 𝑧 ( 𝑡 ) = 𝛽 𝑑 𝑡 𝑚 𝑧 ( 𝑡 𝜏 ) 𝑑 𝑧 ( 𝑡 ) . ( 3 . 1 1 ) If 𝛽 / 𝑚 𝑑 < 0 , then 1 < 1 ; according to Lemma 2.4, we obtain l i m 𝑡 𝑧 ( 𝑡 ) = 0 .
Since 𝑦 ( 𝑠 ) = 𝑧 ( 𝑠 ) = 𝜙 2 ( 𝑠 ) > 0 for all 𝑠 [ 𝜏 , 0 ] , by the comparison theorem in differential equation and the positivity of the solution, we have that 𝑦 ( 𝑡 ) 0 as 𝑡 .
In the following, we suppose that 𝛽 / 𝑚 𝑑 . Since 1 < 1 , we have that 𝛽 𝐾 2 e x p ( 2 l n ( 1 𝑝 ) / ( e x p ( 𝑟 𝑇 ) 1 ) ) 1 + 𝑚 𝐾 2 e x p ( 2 l n ( 1 𝑝 ) / ( e x p ( 𝑟 𝑇 ) 1 ) ) 𝑑 < 0 . ( 3 . 1 2 )
Because the function 𝑥 2 / ( 1 + 𝑚 𝑥 2 ) is monotonically increasing with respect to 𝑥 , we choose a sufficiently small positive constant 𝜖 such that 𝛽 𝐾 e x p ( l n ( 1 𝑝 ) / ( e x p ( 𝑟 𝑇 ) 1 ) ) + 𝜖 2 1 + 𝑚 𝐾 e x p ( l n ( 1 𝑝 ) / ( e x p ( 𝑟 𝑇 ) 1 ) ) + 𝜖 2 𝑑 < 0 . ( 3 . 1 3 )
Noting that 𝑑 𝑥 ( 𝑡 ) / 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n ( 𝐾 / 𝑥 ( 𝑡 ) ) , Δ 𝑥 ( 𝑡 ) = 𝑝 𝑥 ( 𝑡 ) for 𝑛 𝑇 < 𝑡 ( 𝑛 + 1 ) 𝑇 , then we consider the following comparison system: 𝑑 𝑧 ( 𝑡 ) 𝐾 𝑑 𝑡 = 𝑟 𝑧 ( 𝑡 ) l n 𝑧 𝑡 𝑧 ( 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , + 𝑧 0 = ( 1 𝑝 ) 𝑧 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 , + 0 = 𝑥 + . ( 3 . 1 4 )
From Section 3.1, we see that 𝑧 ( 𝑡 ) = e x p l n 𝐾 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 ( 𝑡 ( 𝑛 1 ) 𝑇 ) ( 𝑛 1 ) 𝑇 < 𝑡 𝑛 𝑇 , ( 3 . 1 5 ) which is unique a globally asymptotically stable positive T-periodic solution of system (3.14).
There exist a positive integer 𝑛 1 and an arbitrarily small positive constant 𝜖 for all 𝑡 𝑛 1 𝑇 such that 𝑥 ( 𝑡 ) 𝑧 ( 𝑡 ) + 𝜖 e x p l n 𝐾 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 𝑇 + 𝜖 = 𝜎 . ( 3 . 1 6 )
From (3.16) and the second equation of system (1.3), for 𝑡 > 𝑛 1 𝑇 + 𝜏 we have 𝑑 𝑦 ( 𝑡 ) 𝑑 𝑡 𝛽 𝜎 2 1 + 𝑚 𝜎 2 𝑦 ( 𝑡 𝜏 ) 𝑑 𝑦 ( 𝑡 ) . ( 3 . 1 7 )
Considering the following comparison equation: 𝑑 𝑧 ( 𝑡 ) = 𝑑 𝑡 𝛽 𝜎 2 1 + 𝑚 𝜎 2 𝑧 ( 𝑡 𝜏 ) 𝑑 𝑧 ( 𝑡 ) , ( 3 . 1 8 ) we can see (3.13), that is, 𝛽 𝜎 2 / ( 1 + 𝑚 𝜎 2 ) < 𝑑 ; according to Lemma 2.4, we obtain l i m 𝑡 𝑧 ( 𝑡 ) = 0 .
Since 𝑦 ( 𝑠 ) = 𝑧 ( 𝑠 ) = 𝜙 2 ( 𝑠 ) > 0 for all 𝑠 [ 𝜏 , 0 ] , by the comparison theorem in differential equation and the positivity of the solution, we have that 𝑦 ( 𝑡 ) 0 as 𝑡 .
Next, we will prove that 𝑥 ( 𝑡 ) 𝑥 ( 𝑡 ) as 𝑡 . Without loss of generality, we may assume that 0 < 𝑦 ( 𝑡 ) 𝜖 for all 𝑡 0 . By the first equation of system (1.3), we have that 𝑑 𝑥 ( 𝑡 ) 𝐾 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n 𝑥 ( 𝑡 ) 𝑥 2 𝐾 ( 𝑡 ) 𝑦 ( 𝑡 ) 𝑟 𝑥 ( 𝑡 ) l n 𝑥 ( 𝑡 ) 𝐾 𝑥 ( 𝑡 ) 𝜖 𝑟 𝑥 ( 𝑡 ) l n 𝐾 𝑒 𝐾 𝜖 / 𝑟 . 𝑥 ( 𝑡 ) ( 3 . 1 9 )
Then we have that ̃ 𝑧 1 ( 𝑡 ) 𝑥 ( 𝑡 ) as 𝜖 0 as 𝑡 , where ̃ 𝑧 1 ( 𝑡 ) is unique a positive solution of 𝑑 𝑧 1 ( 𝑡 ) 𝑑 𝑡 = 𝑟 𝑧 1 ( 𝑡 ) l n 𝐾 𝑒 𝐾 𝜖 / 𝑟 𝑧 1 𝑧 ( 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , 1 𝑡 + = ( 1 𝑝 ) 𝑧 1 𝑧 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 , 1 0 + 0 = 𝑥 + . ( 3 . 2 0 ) From Section 3.1, we obtain ̃ 𝑧 1 ( 𝑡 ) = e x p l n 𝐾 𝑒 𝐾 𝜖 / 𝑟 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 ( 𝑡 𝑛 𝑇 ) , 𝑛 𝑇 < 𝑡 ( 𝑛 + 1 ) 𝑇 . ( 3 . 2 1 )
By using comparison theorem of impulsive differential equation [11], for any 𝜖 1 > 0 , there exists 𝑇 1 > 0 and such that 𝑡 > 𝑇 1 𝑥 ( 𝑡 ) > ̃ 𝑧 1 ( 𝑡 ) 𝜖 1 . ( 3 . 2 2 ) On the other hand, from the first equation of (1.3), we have that 𝑑 𝑥 ( 𝑡 ) / 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n ( 𝐾 / 𝑥 ( 𝑡 ) ) .
Considering the following comparison system: 𝑑 𝑧 2 ( 𝑡 ) 𝑑 𝑡 = 𝑟 𝑧 2 𝐾 ( 𝑡 ) l n 𝑧 2 ( 𝑧 𝑡 ) , 𝑡 𝑛 𝑇 , 2 𝑡 + = ( 1 𝑝 ) 𝑧 2 𝑧 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 2 0 + 0 = 𝑥 + , ( 3 . 2 3 ) we have that ̃ 𝑧 2 ( 𝑡 ) = e x p l n 𝐾 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 ( 𝑡 𝑛 𝑇 ) 𝑥 , 𝑛 𝑇 < 𝑡 ( 𝑛 + 1 ) 𝑇 , ( 3 . 2 4 ) ( 𝑡 ) < ̃ 𝑧 2 ( 𝑡 ) + 𝜖 1 , ( 3 . 2 5 ) as 𝑡 and ̃ 𝑧 2 ( 𝑡 ) = 𝑥 ( 𝑡 ) . Let 𝜖 1 0 , then it follows from (3.22) and (3.25), that we can see 𝑥 ( 𝑡 ) 𝜖 1 < 𝑥 ( 𝑡 ) < 𝑥 ( 𝑡 ) + 𝜖 1 for 𝑡 large enough, which implies 𝑥 ( 𝑡 ) 𝑥 ( 𝑡 ) as 𝑡 . The proof is completed.

Now, we give the following theorem.

Theorem 3.3. ( 1 ) If 𝛽 / 𝑚 < 𝑑 , then the “predator-extinction” periodic solution ( 𝑥 ( 𝑡 ) , 0 ) is globally attractive.
( 2 ) If 𝛽 / 𝑚 𝑑 and 𝑝 > 𝑝 , 𝑇 < 𝑇 , 𝑚 > 𝑚 , then the “predator-extinction” periodic solution ( 𝑥 ( 𝑡 ) , 0 ) is globally attractive, where the critical values 𝑝 , 𝑇 , 𝑚 are listed in Table 1.

tab1
Table 1: Critical values of some parameters of system (1.3) ( 𝑅 𝟏 < 𝟏 must be satisfied).

4. Permanence

In the above section, we have proved that, when 𝑝 > 𝑝 , 𝑇 < 𝑇 , or 𝑚 > 𝑚 , the “predator-extinction” periodic solution ( 𝑥 ( 𝑡 ) , 0 ) is globally attractive. But in natural world, the predator cannot be eradicated totally. In order to save resources, we need to keep the prey and predator coexisting when the prey does not bring about immense economic losses. Next, we will discuss the permanence of system (1.3).

Definition 4.1. System (1.3) is said to be uniformly persistent if there are positive constants 𝑚 𝑖 ( 𝑖 = 1 , 2 ) and a finite time 𝑇 0 such that for all solutions ( 𝑥 ( 𝑡 ) , 𝑦 ( 𝑡 ) ) with initial values 𝑥 ( 0 + ) > 0 , 𝑦 ( 0 + ) > 0 , 𝑥 ( 𝑡 ) 𝑚 1 , and 𝑦 ( 𝑡 ) 𝑚 2 hold for all 𝑡 > 𝑇 0 .

Definition 4.2. System (1.3) is said to be permanent if there exists a compact region 𝐷 i n t Ω , such that every solution of system (1.3) with initial condition (1.4) will eventually enter and remain in region 𝐷 .
Denote that 2 = 𝛽 𝐾 2 𝑑 𝑚 𝐾 2 . + e x p ( 2 l n ( 1 / 1 𝑝 ) / ( 1 e x p ( 𝑟 𝑇 ) ) ) ( 4 . 1 )

Theorem 4.3. If 2 > 1 , then there exist two positive constants 𝑚 1 and 𝑚 2 such that 𝑥 ( 𝑡 ) 𝑚 1 , 𝑦 ( 𝑡 ) 𝑚 2 for 𝑡 large enough, that is, system (1.3) is uniformly persistent.

Proof. Suppose that ( 𝑥 ( 𝑡 ) , 𝑦 ( 𝑡 ) ) is any positive solution of system (1.3) with initial condition (1.4). From Lemma 2.3, we know that 𝑥 ( 𝑡 ) 𝐾 , 𝑦 ( 𝑡 ) 𝐿 .
Firstly, from the first equation of system (1.3), we know that 𝑑 𝑥 ( 𝑡 ) / 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n ( 𝐾 𝑒 𝐾 𝐿 / 𝑟 / 𝑥 ( 𝑡 ) ) . Now, we consider the following equation: 𝑑 𝑧 3 ( 𝑡 ) 𝑑 𝑡 = 𝑟 𝑧 3 ( 𝑡 ) l n 𝐾 𝑒 𝐾 𝐿 / 𝑟 𝑧 3 𝑧 ( 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , 3 𝑡 + = ( 1 𝑝 ) 𝑧 3 𝑧 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 , 3 0 + 0 = 𝑥 + . ( 4 . 2 )
therefore we have that ̃ 𝑧 3 ( 𝑡 ) = e x p l n 𝐾 𝑒 𝐾 𝐿 / 𝑟 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 ( 𝑡 𝑛 𝑇 ) 𝑛 𝑇 < 𝑡 ( 𝑛 + 1 ) 𝑇 . ( 4 . 3 ) There exists a 𝜖 > 0 small enough such that, for sufficiently large 𝑡 , 𝑥 ( 𝑡 ) > ̃ 𝑧 3 ( 𝑡 ) 𝜖 e x p l n 𝐾 𝑒 𝐾 𝐿 / 𝑟 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝜖 = 𝑚 1 . ( 4 . 4 )
Secondly, we will find 𝑚 2 such that 𝑦 ( 𝑡 ) 𝑚 2 .
The second equation of system (1.3) may be rewritten as follows: 𝑑 𝑦 ( 𝑡 ) = 𝑑 𝑡 𝛽 𝑥 2 ( 𝑡 ) 1 + 𝑚 𝑥 2 𝑑 ( 𝑡 ) 𝑑 𝑦 ( 𝑡 ) 𝑑 𝑡 𝑡 𝑡 𝜏 𝛽 𝑥 2 ( 𝜃 ) 1 + 𝑚 𝑥 2 ( 𝜃 ) 𝑦 ( 𝜃 ) . ( 4 . 5 )
Define 𝑉 ( 𝑡 ) = 𝑦 ( 𝑡 ) + 𝑡 𝑡 𝜏 𝛽 𝑥 2 ( 𝜃 ) 1 + 𝑚 𝑥 2 ( 𝜃 ) 𝑦 ( 𝜃 ) . ( 4 . 6 ) Calculating the derivative of 𝑉 ( 𝑡 ) along the solution of (1.3), it follows from (4.5) that 𝑑 𝑉 ( 𝑡 ) 𝑥 𝑑 𝑡 = 𝛽 2 ( 𝑡 𝜏 ) 1 + 𝑚 𝑥 2 ( 𝑡 𝜏 ) 𝑦 ( 𝑡 𝜏 ) 𝑑 𝑦 ( 𝑡 ) + 𝛽 𝑥 2 ( 𝑡 ) 1 + 𝑚 𝑥 2 ( 𝑡 ) 𝑦 ( 𝑡 ) 𝛽 𝑥 2 ( 𝑡 𝜏 ) 1 + 𝑚 𝑥 2 𝛽 ( 𝑡 𝜏 ) 𝑦 ( 𝑡 𝜏 ) = 𝑑 𝑑 𝑥 2 ( 𝑡 ) 1 + 𝑚 𝑥 2 ( 𝑡 ) 1 𝑦 ( 𝑡 ) . ( 4 . 7 )
Since 2 > 1 , then we can deduce that 𝑟 2 𝐾 2 l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝐾 + l n 2 𝛽 𝑑 𝑚 > 0 . ( 4 . 8 ) Hence, we choose two positive constants 𝑚 2 and 𝜖 1 small enough such that 𝛽 𝑑 𝜁 2 1 + 𝑚 𝜁 2 > 1 , ( 4 . 9 ) where 𝜁 = 𝐾 e x p 𝐾 𝑚 2 𝑟 e x p l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝜖 1 > 0 , 0 < 𝑚 2 < 𝑟 2 𝐾 2 l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝐾 + l n 2 𝛽 𝑑 . 𝑚 ( 4 . 1 0 ) In the following, we will find 𝑚 2 such that 𝑦 ( 𝑡 ) 𝑚 2 . There are two cases.
(1) We claim that the inequality 𝑦 ( 𝑡 ) < 𝑚 2 cannot hold for all 𝑡 > 𝑡 1 ( 𝑡 1 > 0 ) ; otherwise, there is a positive constant 𝑡 1 such that 𝑦 ( 𝑡 ) < 𝑚 2 for all 𝑡 > 𝑡 1 .
From the first equation of system (1.3), we have that 𝑑 𝑥 ( 𝑡 ) 𝑑 𝑡 𝑟 𝑥 ( 𝑡 ) l n 𝐾 𝑒 𝐾 𝑚 2 / 𝑟 . 𝑥 ( 𝑡 ) ( 4 . 1 1 ) We have that 𝑥 ( 𝑡 ) ̃ 𝑧 4 ( 𝑡 ) , where ̃ 𝑧 4 ( 𝑡 ) is unique a positive solution of 𝑑 𝑧 4 ( 𝑡 ) 𝑑 𝑡 = 𝑟 𝑧 4 ( 𝑡 ) l n 𝐾 𝑒 𝐾 𝑚 2 / 𝑟 𝑧 4 ( 𝑧 𝑡 ) , 𝑡 𝑛 𝑇 , 𝑛 𝑁 , 4 𝑡 + = ( 1 𝑝 ) 𝑧 4 𝑧 ( 𝑡 ) , 𝑡 = 𝑛 𝑇 , 𝑛 𝑁 , 4 0 + 0 = 𝑥 + . ( 4 . 1 2 ) From Section 3.1, we have that ̃ 𝑧 4 ( 𝑡 ) = e x p l n 𝐾 𝑒 𝐾 𝑚 2 / 𝑟 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝑒 𝑟 ( 𝑡 𝑛 𝑇 ) 𝑛 𝑇 < 𝑡 ( 𝑛 + 1 ) 𝑇 . ( 4 . 1 3 )
By comparison theory [11], for any 𝜖 1 > 0 , there exists a 𝑇 1 > 0 , for 𝑡 > 𝑇 1 , such that 𝑥 ( 𝑡 ) > ̃ 𝑧 4 ( 𝑡 ) 𝜖 1 e x p l n 𝐾 𝑒 𝐾 𝑚 2 / 𝑟 + l n ( 1 𝑝 ) 1 𝑒 𝑟 𝑇 𝜖 1 = 𝜁 . ( 4 . 1 4 ) From (4.7) and (4.14), we have that 𝑑 𝑉 ( 𝑡 ) 𝛽 𝑑 𝑡 𝑑 𝑑 𝜁 2 ( 𝑡 ) 1 + 𝑚 𝜁 2 ( 𝑡 ) 1 𝑦 ( 𝑡 ) , 𝑡 𝑇 1 . ( 4 . 1 5 )
Denote that 𝑦 𝑙 = m i n 𝑇 𝑡 1 , 𝑇 1 + 𝜏 𝑦 ( 𝑡 ) . ( 4 . 1 6 )
We show that 𝑦 ( 𝑡 ) 𝑦 𝑙 for all 𝑡 𝑇 1 . Otherwise, there exists a nonnegative constant 𝑇 2 such that 𝑦 ( 𝑡 ) 𝑦 𝑙 , for 𝑡 [ 𝑇 1 , 𝑇 1 + 𝜏 + 𝑇 2 ) , 𝑦 ( 𝑇 1 + 𝜏 + 𝑇 2 ) = 𝑦 𝑙 and ̇ 𝑦 ( 𝑇 1 + 𝜏 + 𝑇 2 ) 0 .
Thus, from the second equation of (1.3) and (4.12), (4.14), we can see that 𝑇 𝑑 𝑦 1 + 𝜏 + 𝑇 2 > 𝛽 𝜁 𝑑 𝑡 2 1 + 𝑚 𝜁 2 𝑦 𝑑 𝑙 𝛽 = 𝑑 𝑑 𝜁 2 1 + 𝑚 𝜁 2 𝑦 1 𝑙 > 0 , ( 4 . 1 7 ) which is a contradiction to 𝑦 ( 𝑡 ) 𝐿 . Hence we get that 𝑦 ( 𝑡 ) 𝑦 𝑙 > 0 for all 𝑡 𝑇 1 . Meanwhile, we can see that 𝑑 𝑉 ( 𝑡 ) / 𝑑 𝑡 > 0 , which implies that 𝑉 ( 𝑡 ) , 𝑡 .
This is a contradiction to 𝑉 ( 𝑡 ) ( 1 + 𝛽 𝜏 𝐾 ) 𝐿 . Thus, for any positive constant 𝑡 1 , the inequality 𝑦 ( 𝑡 ) 𝑚 2 cannot hold for all 𝑡 𝑡 1 .
(2) If 𝑦 ( 𝑡 ) 𝑚 2 holds for all 𝑡 large enough, then our aim is obtained. Otherwise, 𝑦 ( 𝑡 ) is oscillatory about 𝑚 2 . Let 𝑚 2 = m i n { 𝑚 2 / 2 , 𝑚 2 𝑒 𝑑 𝜏 } .
We will prove that 𝑦 ( 𝑡 ) 𝑚 2 . There exist two positive constants 𝑡 , 𝜒 such that 𝑦 ( 𝑡 ) = 𝑦 ( 𝑡 + 𝜒 ) = 𝑚 2 and 𝑦 ( 𝑡 ) < 𝑚 2 for 𝑡 [ 𝑡 , 𝑡 + 𝜒 ] .
When 𝑡 is large enough, the inequality 𝑥 ( 𝑡 ) > 𝜎 holds true for 𝑡 [ 𝑡 , 𝑡 + 𝜒 ] .
Since 𝑦 ( 𝑡 ) is continuous and bounded and does not have impulsive effort, we conclude that 𝑦 ( 𝑡 ) is uniformly continuous. There exist a constant 𝑇 3 ( 0 < 𝑇 3 < 𝜏 and 𝑇 3 is independent to 𝑡 ) such that 𝑦 ( 𝑡 ) > 𝑚 2 / 2 for all 𝑡 [ 𝑡 , 𝑡 + 𝑇 3 ] .
If 𝜒 < 𝑇 3 , then our aim is obtained.
If 𝑇 3 < 𝜒 𝜏 , from the second equation of (1.3), we have that 𝑑 𝑦 ( 𝑡 ) / 𝑑 𝑡 𝑑 𝑦 ( 𝑡 ) for 𝑡 [ 𝑡 , 𝑡 + 𝜒 ] , and then we have that 𝑦 ( 𝑡 ) 𝑚 2 𝑒 𝑑 𝜏 for 𝑡