1. Introduction
As we know, investigating the effect of toxic substances on ecological communities is significant. More and more researchers are engaged in this aspect and have obtained lots of excellent results. Smith [1] firstly introduced the effects of toxic substances into a two-species Lotka-Volterra competitive system under an assumption that each species produces a toxic substance to the other only when the other is present. Chattopadhyay [2] obtained a set of sufficient conditions which ensure the system admits a unique globally stable positive equilibrium by constructing a suitable Lyapunov function. Recently, a developing number of scholars pay attention to the effects of toxic substances on competitive model and many distinguished results have been given(see [3–9] and the references cited therein). For example, Chen et al. [3] and Xie et al. [4] considered the global stability of the interior equilibrium point of the competition model of plankton allelopathy. The authors [5–9] discussed the extinction for the competitive model with the effect of toxic substances.
However, we notice that all the above [1–10] only pay attention to the continuous or discrete system. In fact, the actual real system is easily disturbed by human exploring activities such as planting and harvesting. Different from the continuous or discrete systems, the theory of impulsive differential equation is much richer and has a wider application when analyzing many real-world phenomena [11, 12]. Recently, various kinds of systems with impulse have attracted a lot of attention [13–19], to name just a few. Chen et al. [20] considered the extinction property of the following impulsive competitive system with toxic substance:
(1)
x
˙
1
t
=
x
1
t
b
1
t
-
a
11
t
x
1
t
-
a
12
t
x
2
t
-
d
1
t
x
1
t
x
2
t
x
˙
2
t
=
x
2
t
b
2
t
-
a
21
t
x
1
t
-
a
22
t
x
2
t
-
d
2
t
x
1
t
x
2
t
t
≠
τ
k
,
k
∈
N
,
x
1
τ
k
+
=
1
+
h
1
k
x
1
τ
k
,
x
2
τ
k
+
=
1
+
h
2
k
x
2
τ
k
,
where
x
i
(
t
)
,
i
=
1,2
, denotes the population density of the
i
-th phytoplankton species at time
t
for a common pool of resources. The terms
d
1
(
t
)
x
1
2
(
t
)
x
2
(
t
)
and
d
2
(
t
)
x
1
(
t
)
x
2
2
(
t
)
denote the effects of toxic substances which means that each species produces a toxic substance to the other. And
b
i
(
t
)
and
a
i
j
(
t
)
,
i
,
j
=
1,2
, are all continuous functions which are bounded above and below by positive constants. Here the jump conditions reflect the possibility of impulsive effects on the species
x
1
and
x
2
.
h
i
k
x
i
(
τ
k
)
<
0
may present the impulsive harvesting amount of the species at
t
=
τ
k
, while
h
i
k
x
i
(
τ
k
)
>
0
, the perturbations, may stand for the impulsive stocking amount of the species at
t
=
τ
k
.
That is, the impulsive harvesting amount (or stocking amount) is proportional to the current density of the species.
Given a function
g
(
t
)
, let
g
L
and
g
M
denote
i
n
f
-
∞
<
t
<
+
∞
g
(
t
)
and
s
u
p
-
∞
<
t
<
+
∞
g
(
t
)
, respectively.
Also, throughout this paper, the following conditions are assumed.
(
A
1
)
b
i
(
t
)
,
a
i
j
(
t
)
(
i
,
j
=
1,2
)
,
d
i
(
t
)
are continuous
T
-
periodic functions such that
b
i
L
>
0
,
a
i
j
L
>
0
,
d
i
L
>
0
and
N
=
{
1,2
,
…
}
.
A
2
{
τ
k
}
k
∈
N
satisfies
0
=
τ
0
<
τ
1
<
τ
2
<
⋯
<
τ
k
<
⋯
and
l
i
m
k
→
∞
τ
k
=
+
∞
.
(
A
3
)
h
i
k
(
i
=
1,2
,
k
∈
N
)
are constants and there exists an integer
q
>
0
such that
h
i
(
k
+
q
)
=
h
i
k
,
τ
k
+
q
=
τ
k
+
T
.
A natural constraint in this case is
1
+
h
i
k
>
0
,
i
=
1,2
,
k
∈
N
.
(
A
4
)
b
1
L
¯
=
b
1
L
+
1
/
T
∑
k
=
1
q
l
n
(
1
+
h
1
k
)
>
0
and
b
2
M
¯
=
b
2
M
+
1
/
T
∑
k
=
1
q
l
n
(
1
+
h
2
k
)
>
0
.
Chen et al.[20] proved that if the coefficients of system (1) satisfy
H
1
b
1
L
¯
a
22
L
>
b
2
M
¯
a
12
M
,
b
1
L
¯
a
21
L
≥
b
2
M
¯
a
11
M
,
b
1
L
¯
d
2
L
≥
b
2
M
¯
d
1
M
,
then the species
x
2
of system (1) will be driven to extinction while the species
x
1
will be globally attractive with any positive solution of an impulsive logistic equation.
In [21], the author considered the following Lotka-Volterra type of model for two interesting phytoplankton species:
(2)
x
˙
1
t
=
x
1
t
r
1
-
a
1
x
1
t
-
b
1
x
2
t
-
c
1
x
1
t
x
2
t
,
x
˙
2
t
=
x
2
t
r
2
-
a
2
x
1
t
-
b
2
x
2
t
.
Since many biological or environmental parameters are subject to fluctuate with time, a plausible mathematical model may take the following form:
(3)
x
˙
1
t
=
x
1
t
r
1
t
-
a
1
t
x
1
t
-
b
1
t
x
2
t
-
c
1
t
x
1
t
x
2
t
,
x
˙
2
t
=
x
2
t
r
2
t
-
a
2
t
x
1
t
-
b
2
t
x
2
t
.
To the best of the authors’ knowledge, until this day, there exist no results which focus on the dynamic behaviors of the above system with impulse. All the facts stated above motivate us to investigate the extinction problem for a Lotka-Volterra competitive system with impulse and one toxin producing phytoplankton, i.e., the following system:
(4)
x
˙
1
t
=
x
1
t
b
1
t
-
a
11
t
x
1
t
-
a
12
t
x
2
t
-
d
1
t
x
1
t
x
2
t
x
˙
2
t
=
x
2
t
b
2
t
-
a
21
t
x
1
t
-
a
22
t
x
2
t
t
≠
τ
k
,
k
∈
N
,
x
1
τ
k
+
=
1
+
h
1
k
x
1
τ
k
,
x
2
τ
k
+
=
1
+
h
2
k
x
2
τ
k
.
with initial condition
(5)
x
1
0
>
0
,
x
2
0
>
0
.
Here
x
1
(
t
)
denotes the density of nontoxic phytoplankton. In other words, the first species could not produce toxic substance since
d
2
(
t
)
=
0
.
And the second species
x
2
(
t
)
is toxic liberating phytoplankton which could produce toxic substance due to
d
1
(
t
)
≠
0
. And all the coefficients in the above system have the same restriction as that of system (1). Comparing (4) with (1), one could find that (4) is the special case of (1), i.e.,
d
2
(
t
)
=
0
.
Hence, one may conjecture that the results of [20] could be applied directly to system (4) in this paper. However, unfortunately, one can notice that the inequality
(6)
b
1
L
¯
d
2
L
≥
b
2
M
¯
d
1
M
in
H
1
no longer holds since
d
2
(
t
)
=
0
. As a result, for (4), it is very necessary to figure out a new sufficient condition to guarantee the extinction of the species
x
2
and the global attractivity of the species
x
1
.
The new sufficient condition which will be presented in the next section is different from
H
1
. And the obtained results in the next section will be a good complement of those in [20].
By the basic theories of impulsive differential equations in [22, 23], system (4) has a unique solution
x
(
t
)
=
x
(
t
,
t
0
)
∈
P
C
(
[
0
,
+
∞
)
,
R
2
)
and
P
C
(
[
0
,
+
∞
)
,
R
2
)
=
{
ϕ
:
[
0
,
+
∞
)
→
R
2
,
ϕ
i
s
c
o
n
t
i
n
u
o
u
s
f
o
r
t
≠
τ
k
.
A
l
s
o
ϕ
(
τ
k
-
)
a
n
d
ϕ
(
τ
k
+
)
e
x
i
s
t
,
a
n
d
ϕ
(
τ
k
-
)
=
ϕ
(
τ
k
+
)
,
k
∈
N
}
for each initial value
x
(
0
)
=
x
0
∈
R
2
+
.
The organization of this paper is as follows. In Section 2, necessary preliminaries are presented. In Section 3, we give the sufficient conditions ensuring the extinction and global attractivity of the system. In Section 4, numerical simulations are presented to illustrate the feasibility of our main results. In the last section, we give a brief discussion.
2. Lemmas
In this section, we state the following lemmas [20] which will be useful in the proof of our main results.
Lemma 1.
Assume that
m
∈
P
C
[
R
+
,
R
]
with points of discontinuity at
t
=
τ
k
and is left continuous at
t
=
τ
k
,
k
∈
N
, and that
(7)
D
+
m
t
≤
g
t
,
m
t
,
t
≠
τ
k
,
m
τ
k
+
≤
ϕ
k
m
τ
k
,
where
g
∈
C
[
R
+
×
R
+
,
R
]
,
ϕ
k
∈
C
[
R
,
R
]
and
ϕ
k
(
u
)
is nondecreasing in
u
for each
k
∈
N
.
Let
r
(
t
)
be the maximal solution of the scalar impulsive differential equation:
(8)
u
˙
t
=
g
t
,
u
t
,
t
≠
τ
k
,
u
τ
k
+
=
ϕ
k
u
τ
k
≥
0
,
τ
k
>
t
0
,
k
∈
N
,
u
t
0
+
=
u
0
,
existing on
t
0
,
∞
, then
m
(
t
0
+
)
≤
u
0
implies
m
(
t
)
≤
r
(
t
)
,
t
≥
t
0
.
Lemma 2.
Consider the following impulsive system:
(9)
x
˙
t
=
x
t
a
t
-
b
t
x
t
,
t
≠
τ
k
,
x
τ
k
+
=
1
+
c
k
x
τ
k
,
k
∈
N
,
where
a
(
t
)
,
b
(
t
)
are continuous
T
-
periodic functions and
b
(
t
)
>
0
, and there exists an integer
q
>
0
such that
c
k
+
q
=
c
k
,
t
k
+
q
=
t
k
+
T
and
c
k
>
-
1
for
k
∈
N
.
Then (9) has a unique positive
T
-
periodic solution
θ
[
a
,
b
]
if and only if
(10)
∏
k
=
1
q
1
+
c
k
exp
∫
0
T
a
t
d
t
>
1
,
which is equivalent to the inequality
∫
0
T
a
(
t
)
d
t
+
∑
k
=
1
q
l
n
(
1
+
c
k
)
>
0
.
Moreover, the unique positive
T
-
periodic solution
θ
[
a
,
b
]
is globally asymptotically stable.
Lemma 3.
Let
X
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
)
T
be any solution of system (4) such that
x
i
(
0
+
)
>
0
, then
x
i
(
t
)
>
0
for all
t
≥
0
.
Lemma 4.
Let
X
(
t
)
=
(
x
1
(
t
)
,
x
2
(
t
)
)
T
be any solution of system (4) such that
x
i
(
0
+
)
>
0
,
i
=
1,2
, then we have
(11)
l
i
m
s
u
p
t
→
+
∞
x
i
t
≤
θ
b
i
,
a
i
i
M
,
i
=
1,2
,
where
θ
[
b
i
,
a
i
i
]
is the unique positive solution of the following Logistic equation with impulse
(12)
x
′
t
=
x
t
b
i
t
-
a
i
i
t
x
t
,
t
≠
τ
k
,
k
∈
N
,
x
τ
k
+
=
1
+
h
i
k
x
τ
k
,
3. Main Results
In this section, we present our main results for system (4). In the following, we firstly present the extinction of the species
x
2
.
Theorem 5.
If there exists
ρ
∈
[
0,1
]
such that the following condition
H
2
b
1
L
¯
a
22
L
>
b
2
M
¯
a
12
M
+
ρ
d
1
M
θ
b
1
,
a
11
M
,
b
1
L
¯
a
21
L
>
b
2
M
¯
a
11
M
+
1
-
ρ
d
1
M
θ
b
2
,
a
22
M
holds, then any positive solution
(
x
1
(
t
)
,
x
2
(
t
)
)
T
of system (4) satisfies
(13)
l
i
m
t
→
+
∞
x
2
t
=
0
.
Proof.
By
H
2
, we can choose positive constants
γ
and
η
such that
(14)
b
1
L
¯
b
2
M
¯
>
η
γ
>
a
12
M
+
ρ
d
1
M
θ
b
1
,
a
11
M
a
22
L
,
b
1
L
¯
b
2
M
¯
>
η
γ
>
a
11
M
+
1
-
ρ
d
1
M
θ
b
2
,
a
22
M
a
21
L
.
So from Lemma 4, there exist positive constants
ε
,
T
0
such that
x
i
(
t
)
≤
θ
[
b
i
,
a
i
i
]
+
ε
and
(15)
η
a
22
t
-
γ
a
12
t
+
ρ
d
1
t
θ
b
1
,
a
11
+
ε
>
0
,
η
a
21
t
-
γ
a
11
t
+
1
-
ρ
d
1
t
θ
b
2
,
a
22
+
ε
>
0
when
t
≥
T
0
.
Also there exists
δ
>
0
such that
(16)
η
b
2
t
+
1
T
∑
k
=
1
q
ln
1
+
h
2
k
-
γ
b
1
t
+
1
T
∑
k
=
1
q
ln
1
+
h
1
k
<
-
δ
T
<
0
.
From system (4) and (15), we have
(17)
d
d
t
ln
x
2
t
η
x
1
t
γ
=
η
b
2
t
-
γ
b
1
t
-
η
a
22
t
-
γ
a
12
t
x
2
t
-
η
a
21
t
-
γ
a
11
t
x
1
t
+
γ
d
1
t
x
1
t
x
2
t
≤
η
b
2
t
-
γ
b
1
t
-
η
a
22
t
-
γ
a
12
t
-
γ
ρ
d
1
t
θ
b
1
,
a
11
+
ε
x
2
t
-
η
a
21
t
-
γ
a
11
t
-
γ
1
-
ρ
d
1
t
θ
b
2
,
a
22
+
ε
x
1
t
≤
η
b
2
t
-
γ
b
1
t
,
t
≠
τ
k
,
t
>
T
0
and
(18)
ln
x
2
τ
k
+
η
x
1
τ
k
+
γ
=
ln
1
+
h
2
k
η
1
+
h
1
k
γ
+
ln
x
2
τ
k
η
x
1
τ
k
γ
.
Integrating both sides of (17) over internals
[
T
0
,
T
0
+
τ
1
)
,
[
T
0
+
τ
1
,
T
0
+
τ
2
)
,
[
T
0
+
τ
2
,
T
0
+
τ
3
)
,
[
T
0
+
τ
σ
-
1
,
T
0
+
τ
σ
)
and
[
T
0
+
τ
σ
,
t
)
, respectively, and adding the
σ
inequalities, for
t
∈
[
T
0
+
τ
σ
,
T
0
+
τ
σ
+
1
)
and
T
0
+
τ
σ
∈
[
m
1
T
,
(
m
1
+
1
)
T
)
,
m
1
∈
N
, we obtain
(19)
ln
x
2
t
η
x
1
t
γ
-
ln
x
2
T
0
η
x
1
T
0
γ
<
∫
T
0
t
η
b
2
t
-
γ
b
1
t
d
t
+
ln
∏
T
0
<
τ
k
<
t
1
+
h
2
k
η
∏
T
0
<
τ
k
<
t
1
+
h
1
k
γ
=
m
1
∫
0
T
η
b
2
t
-
γ
b
1
t
d
t
+
m
1
η
∑
k
=
1
q
ln
1
+
h
2
k
-
m
1
γ
∑
k
=
1
q
ln
1
+
h
1
k
+
∫
m
1
T
t
η
b
2
t
-
γ
b
1
t
d
t
+
ln
∏
m
1
T
<
τ
k
<
t
1
+
h
2
k
η
∏
m
1
T
<
τ
k
<
t
1
+
h
1
k
γ
≤
m
1
η
T
b
2
M
+
∑
k
=
1
q
ln
1
+
h
2
k
-
m
1
γ
T
b
1
L
+
∑
k
=
1
q
ln
1
+
h
1
k
+
B
≤
-
m
1
δ
+
B
,
where
B
=
m
a
x
0
≤
s
≤
T
(
∫
0
s
(
η
b
2
(
t
)
-
γ
b
1
(
t
)
)
d
t
+
l
n
∏
0
≤
τ
k
<
s
(
1
+
h
2
k
)
η
)
/
∏
0
≤
τ
k
<
s
(
1
+
h
1
k
)
γ
)
.
This shows that
(20)
x
2
t
η
≤
x
1
t
γ
exp
-
m
1
δ
+
B
x
2
T
0
η
x
1
T
0
γ
.
According to Lemma 4, we notice that
x
1
(
t
)
is ultimately upper bounded; hence, we obtain
l
i
m
t
→
+
∞
x
2
(
t
)
=
0
.
This completes the proof of Theorem 5.
In the above, we have provided the sufficient condition to ensure that species
x
2
is extinct. Similar to the proof of Theorems 3.2 and 3.3 in [20], we can also obtain the following result about the global attractivity of the species
x
1
.
Theorem 6.
Assume that
H
2
holds, then the species
x
1
is globally attractive; i.e., for any positive solution
(
x
1
(
t
)
,
x
2
(
t
)
)
T
of system (4) and any positive solution
x
(
t
)
of the impulsive logistic equation,
(21)
x
˙
t
=
x
t
b
1
t
-
a
11
t
x
t
,
t
≠
τ
k
,
k
∈
N
,
x
τ
k
+
=
1
+
h
1
k
x
τ
k
,
one has
(22)
l
i
m
t
→
+
∞
x
1
t
-
x
t
=
0
.
Remark 7.
For system (4), condition
H
1
is not satisfied due to
d
2
(
t
)
=
0
. As a result, we cannot easily deduce the extinction of the species
x
2
and the globally attractivity of the species
x
1
from Theorem 3.1-3.3 in [20]. However, in Theorems 5–6 of this paper, we have obtained the sufficient condition
H
2
under which the corresponding dynamic behaviors for (4) are presented. Also we can find that conditions
H
1
and
H
2
cannot be derived from each other. In other words, the above two conditions are complementary. Thus, the main results in this paper are indeed a good complement of those in [20].
Remark 8.
According to the main results in [20], for system (4), one can easily induce the extinction of the species
x
1
and the global attractivity of the species
x
2
under the condition
b
2
L
¯
a
12
L
>
b
1
M
¯
a
22
M
,
b
2
L
¯
a
11
L
≥
b
1
M
¯
a
21
M
. This paper pays attention to figure out a new sufficient condition to guarantee the extinction of the species
x
2
and the global attractivity of the species
x
1
.
The obtained results show that the species
x
2
can still be extinct even though the species
x
1
do not produce toxic substance.