Permanence and Globally Asymptotic Stability of Cooperative System Incorporating Harvesting ()
1. Introduction
Permanence, stability and periodic solution for LotkaVolterra models had been extensively investigated by many authors (see [1-8] and the references therein). Jorge Rebaza [1] had discussed the dynamic behaviors of predator-prey model with harvesting and refuge
(1)
he obtained that harvesting and refuge affected the stability of some coexistence equilibrium and periodic solutions of model (1), where
was a continuous threshold policy harvesting function. Motivated by Jorge’s work, we consider the following cooperative system incorporating harvesting
(2)
where
and
denote the densities of two populations at time
. The parameters
are all positive constants.
Definition 1 [2]
is called asymptotically
- periodic function, if
and it satisfies
, where
is continuous periodic function with periodic
and
.
We will discuss our problems in the region
where
.
2. Permanence of System
Definition 2 [2] If there are positive constants
such that each positive solution
of system (2) satisfies
![](https://www.scirp.org/html/5-5300516\1f37ac05-8372-4d2b-b8e2-b9d78fe67add.jpg)
![](https://www.scirp.org/html/5-5300516\1ea3837a-7c1b-4ef6-850e-73ae12b487b2.jpg)
Then system (2) is persistent. If the system is not persistent, then system (2) is called non-persistent.
Lemma 1 If
, then system (2) is persistent.
Proof. By the first equation of (2) and the comparison theorem, we get
it implies that
.
For any
there exists a
, as
, it then follows
![](https://www.scirp.org/html/5-5300516\01ff856b-876f-44b9-a27a-ee1e50e98a13.jpg)
Similarly, we have
. By the discussion above, for any ![](https://www.scirp.org/html/5-5300516\a38436fd-9dbe-4d4b-acad-e06cebc92cb0.jpg)
there exists a
, as
, it yields that ![](https://www.scirp.org/html/5-5300516\570433d3-c13d-48da-ac00-58604b25f289.jpg)
On the other hand, we have
![](https://www.scirp.org/html/5-5300516\dbd543c4-c56f-4dd4-bd6f-65df732c42d9.jpg)
.
By the comparison theorem, and letting
, one gets that
![](https://www.scirp.org/html/5-5300516\59ae97e0-5306-432d-9ad1-e9e2cc19c299.jpg)
.
By Definition 2, system (2) is persistent. □
3. Equilibrium Points and Stability
If
, then the equilibrium points of (2) are
![](https://www.scirp.org/html/5-5300516\a063580b-94a3-49d7-a2ee-74d31af722ce.jpg)
![](https://www.scirp.org/html/5-5300516\69f8128d-6f87-4661-a9f1-776b814bd9d6.jpg)
![](https://www.scirp.org/html/5-5300516\5f320c54-8726-4164-9c47-b08e00c80ecd.jpg)
where
(3)
![](https://www.scirp.org/html/5-5300516\cbe3ef8d-6124-42a6-8ff4-873161eae0a5.jpg)
![](https://www.scirp.org/html/5-5300516\a3a18200-d21b-4717-a1c5-6a11f4df9cdc.jpg)
![](https://www.scirp.org/html/5-5300516\3a6902a7-11e3-494b-9b76-26b1e8dc98c7.jpg)
.
The general Jacobian matrix of (2) is given by
.
The characteristic equation of system (2) at
is
, this immediately indicates that
is always unstable.
The characteristic equation of system (2) at
is
, by the condition
, one then gets that
is a saddle point.
The characteristic equation of system (2) at
is
, we derive that
is a saddle point.
The characteristic equation of system (2) at
takes the form
![](https://www.scirp.org/html/5-5300516\41dc6f32-c1e9-4f8d-9d03-16d4784387c8.jpg)
it is easy to check that
, then
, thus
is locally asymptotically stable.
Theorem 1 If
![](https://www.scirp.org/html/5-5300516\f74ddbe9-5155-4e9f-8948-a1050387039a.jpg)
![](https://www.scirp.org/html/5-5300516\a91a8d3b-7758-4809-8b06-7738ce1d909a.jpg)
![](https://www.scirp.org/html/5-5300516\6064a329-5dd5-48d2-b49e-c206c7df15f6.jpg)
then the positive equilibrium point
of system (2) is globally asymptotically stable, where
can be found in Lemma 1.
Proof. Define a Lyapunov function
![](https://www.scirp.org/html/5-5300516\5969ec6f-4a98-4662-93cd-3dda0139e078.jpg)
it then yields that
![](https://www.scirp.org/html/5-5300516\ec643ed9-db13-4a28-8302-4e3e3337d3b8.jpg)
by the conditions of theorem 1, thus,
. The positive equilibrium point
of system (2) is globally asymptotically stable.
4. Existence and Uniqueness of Solutions
Next, we will discuss a nonautonomous system
(4)
where
are positive continuous bounded asymptotically periodic functions with period
. The initial data of (4) is given by
. (5)
The solution of (4) with initial data (5) is denoted by
,
,
.
For a continuous function
defined on ![](https://www.scirp.org/html/5-5300516\16ff4cb1-d6af-45fe-a111-c0074fddddef.jpg)
define
.
Definition 3 [2] If there exists a
, for any
,
, there exists a
![](https://www.scirp.org/html/5-5300516\8f267257-11ad-4c98-b328-b067b5c76b7e.jpg)
such that
for
, then the solution
is called ultimately bounded.
Let us consider the following asymptotically periodic system
(6)
where
. Set
,
![](https://www.scirp.org/html/5-5300516\08108bec-001a-447e-a034-bb448cfc1091.jpg)
![](https://www.scirp.org/html/5-5300516\3690f269-3044-4331-8cbe-cab0a2131a5b.jpg)
In order to discuss the existence and uniqueness of asymptotically periodic solution of system (6), we can consider the adjoint system
(7)
Lemma 2 If
![](https://www.scirp.org/html/5-5300516\f775f22d-50c2-4649-8d89-4183a4e8a959.jpg)
and
![](https://www.scirp.org/html/5-5300516\c93af1e8-4c41-4774-a2f7-7e324241162b.jpg)
![](https://www.scirp.org/html/5-5300516\43ec9103-6a87-4df6-a927-ce8e8c690079.jpg)
![](https://www.scirp.org/html/5-5300516\9b42bf6a-1e81-481c-8e63-3bb688b93845.jpg)
then the solution of system (4) is ultimately boundedness.
Proof. By the first equation of system (4) and the comparison theorem, one gets that
![](https://www.scirp.org/html/5-5300516\ffaf9704-7824-41d5-a97a-01d95f59f242.jpg)
it then implies that
.
Similarly, we have
.
By the same discussion, one thus gets that
,
![](https://www.scirp.org/html/5-5300516\e772e072-8e86-49c5-a997-e60505636e09.jpg)
Letting
, we have
![](https://www.scirp.org/html/5-5300516\00171f05-37fa-4592-b8e5-920f7e4793a9.jpg)
.
By the Definition 3, the solution of system (4) is ultimately bounded. □
Lemma 3 [2] If
satisfies the following conditions:
1)
, where
and
are continuously positively increasing functions;
2)
where
is a constant;
3) there exists a continuous non-increasing function
, such that for s > 0,
. And as
,
it then follows that
where
is a constant; furthermore, system (6) has a solution
for
and satisfies
.
Then system (6) has a unique asymptotically periodic solution, which is uniformly asymptotically stable.
Theorem 2 If conditions
![](https://www.scirp.org/html/5-5300516\1e4bfcac-955d-4672-86d1-78849e68e75b.jpg)
and
![](https://www.scirp.org/html/5-5300516\ad194e39-6c73-44bd-a8c5-408d25e85a15.jpg)
hold, the conditions of Lemma 2 are satisfied, then system (4) has a unique asymptotically periodic solution, which is uniformly asymptotically stable.
Proof. By Lemma 2, the solutions of system (4) is ultimately bounded. We consider the adjoint system
(8)
Let
![](https://www.scirp.org/html/5-5300516\39e2b4bd-6664-445d-b82b-0317c653e0fb.jpg)
and
be the solution of (8). By the fact
![](https://www.scirp.org/html/5-5300516\fe80452c-e80f-4065-b930-cc31489f6b1f.jpg)
![](https://www.scirp.org/html/5-5300516\0c26f13a-dbd9-4e24-9ceb-e9fa7e9e2ab2.jpg)
where
lies between
and
,
lies between
and
, it then follows
(9)
Define Lyapunov function
, taking
By suing of the inequality
, it is easy to check that 1) and 2) of Lemma 3 are valid. Computing the derivative of
along the solution of system (8), by (9) and
, we get that
![](https://www.scirp.org/html/5-5300516\c103c2a4-37fa-4610-ba19-8d39bf4cc3c5.jpg)
taking
, it yields
, then, system (4) has a unique positive asymptotically periodic solution, which is uniformly asymptotically stable. □
5. Examples and Numerical Simulations
Now, let us consider a autonomous cooperative system incorporating harvesting
, (10)
it is easy to check that
,
,
,
,
the conditions of Theorem 1 are valid, then the positive equilibrium point
of system (2) is globally asymptotically stable in Figures 1 and 2.
6. Conclusions
By analyzing the characteristic roots of a kind of cooperative models (2) incorporating harvesting, the stability of positive equilibrium point
to model (2) is obtained by constructing a suitable Lyapunov function. Our results have shown that the harvesting coefficient
affects the stability and the existence of equilibrium point to model (2).
The related non-autonomous asymptotically periodic cooperative model (4) has been discussed later. Under some conditions, which also depend on model parameters (see Theorem 2), model (4) has a unique asymptotically periodic solution
, which is uniformly
![](https://www.scirp.org/html/5-5300516\953da001-40f1-4839-917a-33be8196b205.jpg)
Figure 1. Positive equilibrium point
of (2) is globally asymptotically stable.
![](https://www.scirp.org/html/5-5300516\f00324df-5f82-49bb-a2cf-95abb50dea9c.jpg)
Figure 2. Solution of (2) is uniformly asymptotically stable.
asymptotically stable. Example model (10) shows the effectiveness of our results.
7. Acknowledgements
Our work is supported by Natural Science Foundation of China (11201075), the Natural Science Foundation of Fujian Province of China (2010J01005).