It is well known that the dynamic relationship between predator and prey has always been one of the main topics in ecology and mathematical ecology. In the past decades, many predator-prey models have been proposed and widely used to describe the food supply relationship between two species [1] [2]. At the same time, it has attracted great attention in many different fields, such as bio-economics. Recently, the interaction of predator-prey with harvesting has been studied. The effect of harvest on population is beneficial to sustainable development and renewable resource management, so many scholars take harvest into account in their models. The capture intensity depends largely on the capture strategy being implemented. Common harvest functions are: constant harvest, proportional harvest and nonlinear harvest. Gupta et al. proposed a predator-prey model with nonlinear predator in harvest [3] and discussed the dynamical properties of the following system:

{ d x = x ( t ) ( r 1 − b 1 x ( t ) ) − a x ( t ) y ( t ) d t , d y = y ( t ) ( − r 2 + η a x ( t ) y ( t ) − h y ( t ) 1 + b y ( t ) ) d t , (1.1)

On the other hand, the growth of species in nature is often limited by environmental factors. Generally speaking, there are two main types of environmental noise: white noise and colored noise. Wenjie Zuo et al. [4] considered the white noise and studied the stationary distribution and periodic solution. However, reading the literature found that studies on the non-linear harvesting of predators and prey are very few literatures [5] [6] [7]. Therefore, the following model is proposed.

{ d x ( t ) = x ( t ) [ r 1 − a 11 ( t ) x ( t ) − a 12 ( t ) y ( t ) − H ( t ) 1 + b ( t ) x ( t ) ] d t                       + σ 1 ( t ) x ( t ) d B 1 ( t ) , d y ( t ) = y ( t ) [ − r 2 + a 21 ( t ) x ( t ) − a 22 ( t ) y ( t ) − h ( t ) 1 + b ( t ) y ( t ) ] d t                       − σ 2 ( t ) y ( t ) d B 2 ( t ) − σ 3 ( t ) h ( t ) y ( t ) 1 + b ( t ) y ( t ) d B 3 ( t ) , (1.2)

In real life, however, ecosystems are often disturbed by human development or by activities related to natural factors such as drought, floods, earthquakes, and planting. In order to describe this phenomenon more accurately, impulses perturbation is added into the model. To sum up, this paper mainly studies the effects of impulse effect and nonlinear harvesting on predator and prey populations, and proposes the following interesting stochastic system.

{ d x ( t ) = x ( t ) [ r 1 ( t ) − a 11 ( t ) x ( t ) − a 12 ( t ) y ( t ) − H ( t ) 1 + b ( t ) x ( t ) ] d t                       + σ 1 ( t ) x ( t ) d B 1 ( t ) , d y ( t ) = y ( t ) [ − r 2 ( t ) + a 21 ( t ) x ( t ) − a 22 ( t ) y ( t ) − h ( t ) 1 + b ( t ) y ( t ) ] d t                       − σ 2 ( t ) y ( t ) d B 2 ( t ) − σ 3 ( t ) h ( t ) y ( t ) 1 + b ( t ) y ( t ) d B 3 ( t ) , } t ≠ t k x ( t k + ) − x ( t k ) = α k x ( t k ) , y ( t k + ) − y ( t k ) = β k y ( t k ) , t = t k , k = 1 , 2 , 3 , ⋯ (1.3)

where x ( t ) and y ( t ) represent the density of prey and predator populations respectively. The parameters r i ( t ) , a i j ( t ) , ( i , j = 1 , 2 ) are positive, and r 1 is the internal growth rate of prey, and r 2 is the mortality rate of predator. a 11 ( t ) and a 22 ( t ) represent the intra-specific competition coefficients of prey and predator populations, respectively. The coefficient a 12 ( t ) is the predator’s capture rate and a 21 ( t ) stands for the rate at which nutrients are converted to predators. In addition, H ( t ) 1 + b ( t ) x ( t ) , h ( t ) 1 + b ( t ) y ( t ) are the nonlinear harvesting.

Throughout this paper, unless otherwise specified, we suppose ( Ω , F , { F t } t ≥ 0 , ℙ ) be a complete probability space with a filtration { F t } t ≥ 0 satisfying the usual conditions and it is right continuous and increasing, while F 0 contains all ℙ -null set. All the coefficients are assumed to be T-periodic continuous functions.

The remainder of this paper is organized as follows. In Section 2, we show that the model (1.3) existence of the global positive solution. In Section 3, sufficient conditions are achieved to guarantee the existence of a positive periodic solution of the stochastic system (1.3) by using Itô’s formula. In Section 4, we discuss the globally attractive of stochastic model (1.3). In Section 5, we use numerical simulation to illustrate our results.

First, to facilitate the analysis that follows, we make the following tags. When f ( t ) is a continuous T-periodic function, we define:

f u = sup t ≥ 0 f ( t ) , f l = inf t ≥ 0 f (t)

Moreover, we assume that a product equals unity if the number of factors is zero.

Definition 2.1. [8] Consider an impulsive stochastic differential equation

{ d x ( t ) = f ( t , x ( t ) ) d t + g ( t , x ( t ) ) d B ( t ) , t ≠ t k , t > 0 , x ( t k + ) − x ( t k ) = α k x ( t k ) , t = t k , k = 1 , 2 , 3 , ⋯ . (2.1)

A stochastic process x ( t ) = ( x 1 ( t ) , x 2 ( t ) , ⋯ , x n ( t ) ) T , t ∈ [ 0 , + ∞ ) is said to be a solution of ISDE (2.1), if x ( t ) satisfies

1) x(t) is F t adapted and is continuous on ( 0 , t 1 ) and each interval ( t k , t k + 1 ) , k ∈ ℕ and f ( t , x ( t ) ) ∈ L 1 ( ℝ + , ℝ n ) , g ( t , x ( t ) ) ∈ L 2 ( ℝ + , ℝ n ) ;

2) x(t) obeys the equivalent integral equation of (2.1) for almost every t ∈ ℝ + \ t k and satisfies the impulsive conditions at each t ∈ ℝ + , k ∈ ℕ a . s . ;

3) For each t k , k ∈ ℕ , x ( t k + ) = lim t → t k + x ( t ) and x ( t k − ) ) = lim t → t k − x ( ( t ) ) exist and x ( t k − ) = x ( t k ) with probability one.

We give the main results of system (1.3) as follows.

Theorem 2.1. For any initial value ( x 0 , y 0 ) ∈ R + 2 the system (1.3) has a unique global positive solution ( x ( t ) , y ( t ) ) for t ≥ 0 and the solution remains in ℝ + with probability one.

Proof. First, we construct the following SDE without impulses:

{ d x 1 ( t ) = x 1 ( t ) [ r 1 ( t ) + 1 T ∑ j = 1 p ln ( 1 + α j ) − a 11 ( t ) A 1 ( t ) x 1 ( t ) − a 12 ( t ) A 2 ( t ) x 2 ( t )                       − H ( t ) 1 + b ( t ) A 1 ( t ) x 1 ( t ) ] d t + σ 1 ( t ) x 1 ( t ) d B 1 ( t ) , d x 2 ( t ) = x 2 ( t ) [ − r 2 ( t ) + 1 T ∑ j = 1 p ln ( 1 + β j ) + a 21 ( t ) A 1 ( t ) x 1 ( t ) − a 22 ( t ) A 2 ( t ) x 2 ( t )                       − h ( t ) 1 + b ( t ) A 2 ( t ) x 2 ( t ) ] d t − σ 2 ( t ) x 2 ( t ) d B 2 ( t ) − σ 3 ( t ) h ( t ) x 2 ( t ) 1 + b ( t ) x 2 ( t ) d B 3 ( t ) , (2.2)

with the initial value ( x 1 , x 2 ) = ( x 0 , y 0 ) , where

A 1 ( t ) = ( ∏ j = 1 p ( 1 + α j ) ) − t T ∏ 0 ≤ t k < t ( 1 + α k ) , A 2 ( t ) = ( ∏ j = 1 p ( 1 + β j ) ) − t T ∏ 0 ≤ t k < t ( 1 + β k )

Then it is obvious that A 1 ( t ) , A 2 ( t ) are positive T-periodic functions. In fact,

A 1 ( t + T ) A 2 ( t ) = ( ∏ j = 1 p ( 1 + α j ) ) − t + T T ∏ 0 ≤ t k < t + T ( 1 + α k ) ( ∏ j = 1 p ( 1 + α j ) ) − t T ∏ 0 ≤ t k < t ( 1 + α k ) = ( ∏ j = 1 p ( 1 + α j ) ) − 1 ∏ t ≤ t k < t + T ( 1 + α k ) . (2.3)

For any t ≥ 0 , there is an integer n, such that

n T ≤ t ≤ ( n + 1 ) T .

The limited mathematical induction procedures, together with t k + p = t k + T , α k + p = α k induce that

t k + n p = t k + ( n − 1 ) p + T = ⋯ = t k + n T ,   α k + n p = α k + ( n − 1 ) p = ⋯ = α k (2.4)

According to [ 0 , T ) ∩ { t k , k ∈ ℤ } = { t 1 , t 2 , ⋯ , t p } , there exists l = { 1 , 2 , ⋯ , p } such that

t l + n p , t l + 1 + n p , ⋯ , t p + n p ∈ [ t , ( n + 1 ) T ) , t 1 + ( n + 1 ) p , t 2 + ( n + 1 ) p , ⋯ , t l − 1 + ( n + 1 ) p ∈ [ ( n + 1 ) T , t + T ) . (2.5)

Thus, combining (2.2)-(2.4), we obtain

A 1 ( t + T ) = A 1 ( t ) ( ∏ j = 1 p ( 1 + α j ) ) − 1 ∏ k = l p ( 1 + α k + n p ) ∏ k = 1 l − 1 ( 1 + α k + ( n + 1 ) p ) = A 1 ( t ) ( ∏ j = 1 p ( 1 + α j ) ) − 1 ∏ k = l p ( 1 + α k + n p ) ∏ k = 1 l − 1 ( 1 + α k ) = A 1 ( t ) ( ∏ j = 1 p ( 1 + α j ) ) − 1 ∏ k = 1 p ( 1 + α k ) = A 1 ( t ) ，

Similarly, A 2 ( t + T ) = A 2 ( t ) .

By the same method as [9] and standard proof [10], Equation (2.2) has a unique global positive Solution ( x 1 ( t ) , x 2 ( t ) ) .

Next we will show that ( x ( t ) , y ( t ) ) is the solution of system (2.2), which is continuous on each interval. For any.

Let

,.

And, for every,

,

.

Similarly, we can show that,

.

Therefore, is a solution that satisfies system (1.3) Finally, we prove the nonnegative uniqueness of the solution of system (1.3) (more details see [11]).

Then the proof is completed.

In this section, we give the existence of the positive periodic solution of the stochastic system (1.3) with impulses. For convenience of readers, we first give the definition of the periodic solution of the impulsive stochastic differential equation in the sense of distribution and the results of the existence of periodic solutions (see [10] [11]).

Definition 3.1. [12] A stochastic process is said to be periodic with period T, if for every finite sequence of numbers the joint distribution of random variables is independent of h, where.

Consider the following periodic stochastic differential equation without impulse:

where is a matrix function, and the matrix are T-periodic in t.

Lemma 3.1. [12] [13] Assume that the system (3.1) has a global solution, and there exists a T-periodic function such that the following conditions hold:

1) on the outside of some compact set;

2) , as.

Then Equation (3.1) has a T-periodic solution.

According to Lemma 3.2, we can obtain the main result in this section.

Theorem 3.1. Assume

(H1):

(H2):,

(H3):,

(H4):.

Then system (1.4) has a positive T-periodic solution.

Proof. We only need to prove the existence of a periodic solution of the equivalent system (2.2) without impulses. The global existence of the solution has been ensured by Theorem 1. Then, we only have to verify the conditions of by Lemma 3.1.

Define a C²-function:

where will be determined later. Here, satisfies

Which and are defined by (H1), (H2). Obviously, are T-periodic functions. And is a bounded function. Thus there is such that;

, (3.5)

In order to confirm the condition (2) of Lemma 3.1, we only need to prove that

, as.

where, here the coefficients of the quadratic term of are all positive.

Next, we verify the condition (1) of Lemma 3.2. By Itô’s formula, we have:

As so that:

Let, , , we have

Then

Then

where:

Let, we take

To confirm the condition (1) of Lemma 3.2, we choose a sufficiently small constant ε such that:

where

Define a bounded open set

.

and denote

, ,

,.

It is obvious that. Next, on must be shown.

Case 1: If, then, we have:

Using (3.8) and (3.9), we obtain

.

Case 2: If, then, we have:

By the definition (3.8) of C and the inequalities (3.9), we have:

.

By Young inequality, we have. Then by equality (3.11), the following inequality is obvious:

Case 3: If, from (3.9) and (3.10), we obtain

.

Case 4: If, from (3.9) and (3.10), we obtain

.

Thus, we obtain on, and the condition (1) of Lemma 3.2 is satisfied. Therefore, by Lemma 3.2, system (1.3) has a positive T-periodic solution.

The proof is confirmed.

Theorem 4.1. [14] [15] If system (1.3) satisfies, , then the system (1.3) is globally attractive.

Proof. Let, be two arbitrary solutions of model (1.3) with initial values.

We defined the following Lyapunov function

Then by calculating the right differential and employing Ito’s formula.

When, we have:

when, we get

Integrating both sides and then taking the expectation yields that

That is

Then, in the view of and that. Thus, it is easy to see from Lemmas 6.1 [15]

The proof is complete.

this section, we will prove our theoretical results by some examples with the help of the Matlab software [16] and reveal the influence of impulses and the white noise.

Example 1.

Let

then

.Thus, the conditions of Theorem 3.1. hold. Then the model (1.3) has a positive -periodic solution. Figure 1 confirms the results.

Example 2. Set. Making condition of the Theorem 4.1 is satisfied. We get that system (1.3) is globally attractive (see Figure 2).

In this paper, we propose a stochastic predator-prey system with nonlinear harvesting and impulsive perturbations. Firstly, we show that there is a unique positive solution in system (1.3). Secondly, the system has a positive periodic solution under a certain condition. Result shows that when the impulses are sufficiently large such that then the predator and prey will tend to exhibit periodicity. It is verified by constructing the appropriate Lyapunov functions and using Itô’s formula. Moreover, these methods used in this study can be extended to more complex and realistic models.

This work was supported by the National Natural Science Foundation of China (11861027) and Natural Science Foundation of Guangxi (2016 GXNSFAA380194).

The authors declare no conflicts of interest regarding the publication of this paper.

Yang, Y.F., Shao, Y.F. and Li, M.W. (2019) Periodic Solution for Stochastic Predator-Prey Systems with Nonlinear Harvesting and Impulses. Advances in Linear Algebra & Matrix Theory, 9, 89-103. https://doi.org/10.4236/alamt.2019.94007