# Generalized predator-prey model with nonlinear impulsive control strategy.

1. IntroductionDuring the last few decades, controlling insect pests of agriculture and insect vectors of important plant has been becoming an increasing important issue all over the world. The pest outbreaks often cause serious ecological and economic problems, which have been a great concern for entomologists and the society. Integrated pest management (IPM) involves choosing appropriate tactics from a range of pest control techniques including biological, cultural, and chemical methods to suit individual cropping systems, pest complexes, and local environments [1-4]. It has been proved both theoretically [2, 5] and experimentally [6, 7] that IPM has been more effective than the biological control or chemical control alone.

Mathematical models can assist in the design and understanding of IPM strategies and mathematical analyses can provide valuable information about how to control pest outbreaks. Recently, many works [7-20] have been devoted to investigate the range of possible ecological interactions between pest and natural enemy to construct and explore population models describing IPM strategies.

However, one of the major assumptions in previous works was that all control tactics such as the instant killing rate of pesticide applications with respect to the pest is a constant. In fact, every community or country has an appropriate or limited capacity for pesticides, farmers, equipment, costs, and so forth, especially for developing countries. In these above factors of limited resources, most notably are Limited Resource Farmers (LRFs); that is, the farmers used pesticides and practiced crop rotation and other practices, but a majority was not familiar or inadequately familiar with IPM. Data was collected from 90 LRFs in the Alabama Black Belt and analyzed by descriptive statistics [21-24]. It reveals that there exist many challenges including health risks to farmers as well as environmental damage.

Understanding resource limitation is critical to effective management and conservation of populations; however, resource limitation is difficult to quantify partly because it is a dynamic process [25]. Meanwhile, a saturation phenomenon of the limited resources is considered in order to better characterize the effect of resource limitation. Based on the above factors, we propose a generalized mathematical model with nonlinear pulse control tactics in order to investigate the effect of limited resources on the outbreak of pest populations. To the best of our knowledge, no work has been done for the effects of resource limitation on a generalized predator-prey model.

The organization of present paper is as follows. Section 2 introduces our model and some preliminaries. In Section 3, by using the method of the differential inequality, qualitative analysis, the discrete dynamical system determined by the stroboscopic map, and Floquet theory, some sets of sufficient conditions, which guarantee the existence and stability of the pest-free periodic solution, are obtained. We choose the Holling Type II functional response curve as an example with aims to investigate how the nonlinear pulse perturbations affect the successful pest control in Section 4. The paper ends with some interesting biological conclusions and numerical bifurcation analyses, which complement the theoretical findings.

2. The Model with Nonlinear Pulse Perturbations and Preliminaries

In this paper, we will extend the classical Lotka-Volterra model [26, 27] to a generalized form. In the following, x and y are the population abundances of the pests and natural enemies, respectively,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

with nonnegative initial conditions ([x.sub.0], [y.sub.0]), where f(x) is the per capita growth rate of the pest population in the absence of natural enemy, [phi](x) is the per capita functional response of the predator, [mu] is the conversion efficiency of the prey to predator, and [delta] is the per capital death rate of the predator population.

Furthermore, we assume that the IPM control tactics are implemented every period T, at which the natural enemies are released and pesticides are applied simultaneously. Moreover, to take account of the resource limitation and saturation effects, we use the Sigmoid function [28, 29] that is to assume that at every period T a nonlinear perturbation p(x) = [p.sub.1]/(1 + [e.sup.-x) or p(y) = [p.sub.2]/(1 + [e.sup.-y]) decreases on pest or natural enemies populations, respectively, and an introduction of constant a for the natural enemies, which does not depend on their population sizes, is considered. These results yield the following discrete process for every time nT:

x([t.sup.+]) = (1 - [[[p.sub.1]]/1 + [e.sup.-x]])x, y([t.sup.+]) = (1 - [[[p.sub.2]]/1 + [e.sup.-y]])y + [sigma], (2)

where n [member of] N = {0, 1, 2, ...} and [p.sub.1] and [p.sub.2] represent the maximal fatality rate for pest and natural enemy, respectively.

Combining (1) and (2) yields the following model concerning IPM strategies:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

In order to investigate the dynamics of the above generalized prey-predator model, we first provide some basic assumptions for the functions f(x) and g(x) as follows.

Hypothesis 1. Let f(x) and [phi](x) be locally Lipschitz continuous on [R.sup.+] such that

(i) f(0) = 0,

(ii) [phi](0) = 0, [phi]'(0) > 0, and [phi](x) > 0 for x > 0,

(iii) the function f(x)/[phi](x) and [phi](x)/x is upper bounded for x > 0.

The biological explanations of those conditions can be found in [30-33].

The following definition and lemma are useful throughout the paper.

Definition 1. The Lambert W function is defined to be a multivalued inverse of the function z [??] [ze.sup.z] satisfying

Lambert W(z) x [e.sup.Lambert W(z)] = z. (4)

It follows from (4) that

Lambert W'(z) = [[Lambert W(z)]/z(1 + Lambert W(z))]. (5)

Note that the function [ze.sup.z] has the positive derivative (z + 1)[e.sup.z]. Define the inverse function of [ze.sup.z] restricted on the interval (-[infinity], -1] and [-1, +[infinity]) to be W(-1, z) and W(0, z), respectively. There are two real branches of Lambert W function shown in Figure 1. For more details about the concepts and properties of the Lambert W function, see Corless et al. [34] and Waldvogel [35].

Lemma 2. Consider the following algebraic equation:

g(z) = ([a.sub.1]z + [a.sub.2])[e.sup.z] + ([a.sub.3]z) + [a.sub.2] = 0, (6)

where [a.sub.1] < 0, [a.sub.2] > 0, [a.sub.3] < 0, and [a.sub.1] < [a.sub.3]. Equation (6) with respect to z has only one positive root provided that [a.sub.1] + [a.sub.2] + [a.sub.3] > 0 and [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII].

The proof of Lemma 2 is given in Appendix A.

3. The Existence of the Pest-Free Periodic Solution and Threshold Condition

3.1. Existence of the Pest-Free Periodic Solution. In order to illustrate the existence of a pest-free periodic solution of model (3), we consider the basic properties of following subsystem firstly:

[??](t) = -[delta]y(t), t [not equal to] nT, y([t.sup.+]) = (1 - [[p.sub.2]/[1 + [e.sup.-y(t)]]])y(t) + [sigma], t = nT, y([0.sup.+]) = [y.sub.0]. (7)

It follows from model (7) that, for any time internal {nT, (n + 1)T], we integrate the first equation of model (7) and get

y{t) = y([nT.sup.+])[e.sup.-[delta](t-nT)]. (8)

It follows from the second equation of model (7) that

y((n + 1)[T.sup.+]) = [1 - [[p.sub.2]/[1 + [e.sup.-y((n+1)T)]]]]y((n + 1)T) + [sigma]. (9)

Substitution of (8) into (9) gives

y((n + 1)[T.sup.+]) = [1 - [[p.sub.2]/[1 + exp{-y([nT.sup.+])[e.sup.-[delta]T]})] x y([nT.sup.+])[e.sup.-[delta]T] + ([sigma]). (10)

That is

y((n + 1)[T.sup.+]) = [[(1 - [p.sub.2])exp(y([nT.sup.+])[e.sup.-[delta]T]) + 1]y([nT.sup.+])[e.sup.- [delta]T]/exp(y([nT.sup.+])[e.sup.-[delta]T]) + 1] + ([sigma]). (11)

Denote y([nT.sup.+]) = [y.sub.n]; then the above equation can be rewritten as the following difference equation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

which is the so-called stroboscopic map of model (7) and describes the relations of the number of nature enemies in the population between any two successive pulse points. Consequently, the existence of the positive steady state of model (12) implies the existence of a positive periodic solution of model (7). Therefore, we first discuss the existence of a positive steady state of (12). Taking the derivative of F([y.sub.n]) with respect to [y.sub.n] yields

F'([y.sub.n]) = {(1 - [p.sub.2]) + [[exp([y.sub.n][e.sup.-[delta]T]) + 1 - [y.sub.n]([e.sup.-[delta]T]) exp([y.sub.n][e.sup.- [delta]T])]/[[exp([y.sub.n][e.sup.-[delta]T]) + 1].sup.2]][p.sub.2]} x [e.sup.-[delta]T], (13)

and it is easy to see that 0 < [absolute value of F'([y.sub.n])] < 1 holds true.

Next, we discuss the positive point of the stroboscopic map (12), denoted by [y.sup.*]; then (12) is rewritten as

([y.sup.*] - [sigma])(exp [y.sup.*]([e.sup.-[delta]T]) + 1) = ((1 - [p.sub.2])exp{[y.sup.*]([e.sup.-[delta]T])} + 1) [y.sup.*]([e.sup.-[delta]T]). (14)

Denote [A.sub.2] = ([y.sup.*])[e.sup.-[delta]T] > 0; then (14) is rewritten as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

where a = 1 - [p.sub.2] - [e.sup.[delta]T] < 0, [sigma] > 0, c = 1 - [e.sup.[delta]T] < 0, and obviously a < c.

According to Lemma 2, (14) has a unique positive root [y.sup.*]. Therefore, there exists a positive fixed point of the stroboscopic map (12), and we denote the positive fixed point by [y.sup.*], which satisfies the following equation:

[y.sup.*] = F([y.sup.*]). (16)

According to the relations between a fixed point of the stroboscopic map (12) and the periodic solution of model (7), we conclude that submodel (7) has a unique nontrivial positive periodic solution, denoted by [y.sub.p](t), and

[y.sub.p](t) = [y.sup.*]([e.sup.-[delta](t-nT)]) for nT < t [less than or equal to] (n + 1)T, (17)

and it follows from 0 < [absolute value of P'([y.sub.p](t))] < 1.

Therefore, we obtain the general expression of the pest-free periodic solution of model (3) over the interval nT < t [less than or equal to] (n + 1)T for all n [member of] N denoted by

([x.sub.p](t), [y.sub.p](t)) = (0, [y.sup.*][e.sup.-[delta](t-nT)]), (18)

where [y.sup.*] is determined by the difference equation (14).

Theorem 3. Model (7) has a positive periodic solution [y.sub.p](t) with period T and for every solution y(t) of (7) such that [absolute value of y(t) - [y.sub.p](t)] [right arrow] 0 as t [right arrow] +[infinity] where [y.sub.p](t) = [y.sup.*][e.sup.-[delta](t-nT)], and [y.sup.*] is the positive root of (14).

By using Theorem 3, we have that the unique fixed point [y.sup.*] of (12) is globally stable. According to the stability of positive periodic solution, it is determined by the stability of positive equilibria of difference equation made by the sequence of impulsive points. So, the positive periodic solution [y.sub.p](t) of model (7) is globally stable.

Therefore, model (3) exists as a unique pest-free periodic solution (0, [y.sub.p](t)).

3.2. Stability of the Pest-Free Periodic Solution. The threshold conditions which guarantee the global stability of pest-free periodic solution ([x.sub.p](t), [y.sub.p](t)) play a key role in pest control. To do this, we first prove the local stability which can be determined by considering the behavior of small amplitude perturbations of the solution. So we denote ([bar.x](t), [bar.y](t)) so that

[bar.x](t) = x(t) - [x.sub.p](t), [bar.y](t) = y(t) - [y.sub.p](t). (19)

Then model (3) becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (20)

Assuming that ([bar.x], [bar.y]) is small enough, we get the linear approximation of the deviation system around the periodic solution ([x.sub.p](t), [y.sub.p](t))

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

In the following, we will present the sufficient condition for the global attractivity of pest-free periodic solution ([x.sub.p](t), [y.sub.p](t)) of model (3).

Theorem 4. The pest-free periodic solution ([x.sub.p](t), [y.sub.p](t)) of model (3) is locally stable in the first quadrant provided that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

and is globally attractive if

[R.sub.2] = [[TM.sub.s][m.sub.s]/[[m.sub.s]([[integral].sup.T.sub.0])[y.sub.p](s)ds - ln(1 - [p.sub.1])]] < 1, (23)

where [M.sub.s] = [sup.sub.x[greater than or equal to]0](f(x)/[phi](x)), [m.sub.s] = [sup.sub.x > 0]([phi](x)/x).

The proof of Theorem 4 is given in Appendix B.

4. Application of a Special Case of Model (3)

In order to show the application of the main results obtained in Theorem 4 and discuss the biological implications of the threshold conditions, we assume that the pest population follows the logistic growth in the absence of predator; that is, f(x) = ax - [bx.sup.2], and we choose the Holling Type II function response for [phi](x); that is, [phi](x) = [alpha]x/(1 + [omega]x). Thus model (3) becomes as the following special model:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

where a, b, [alpha], and [omega] are positive constants and a and b are the intrinsic growth rate of prey and the rate of intraspecific competition or density dependence, respectively.

4.1. Threshold Conditions. It follows from Theorem 4 that model (24) has a unique pest-free periodic solution which can be formulated

([x.sub.p](t), [y.sub.p](t)) = (0, ([y.sup.*])[e.sup.-[delta](t-nT)]), t [member of] (nT, (n + 1)T], (25)

where [y.sup.*] is the positive root of (14).

The threshold conditions for the local stability and global attractivity can be directly obtained from Theorem 4. Note that condition (22) is equivalent to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

and by simple calculation, the above inequality is equivalent to

aT - ([alpha][y.sup.*])[[1 - [e.sup.-[delta]T]]/[delta]] < ln[2/[2 - [p.sub.1]]]. (27)

Therefore, the threshold condition (22) can be rewritten as

[[aT[delta] - [alpha][y.sup.*](1 - [e.sup.-[delta]T])]/[delta]ln(2/(2 - [p.sub.1]))] < 1. (28)

Similarly, the threshold condition (23) can be rewritten as

[[delta]T[M.sub.s][m.sub.s]/[[m.sub.s]([y.sup.*])(1 - [e.sup.-[delta]T]) - [delta] ln (1 - [p.sub.1])]] < 1, (29)

where [M.sub.s] = [(a[omega] + b).sup.2]/4b[omega][alpha], [m.sub.s] = [alpha].

That is, if (28) holds true, then the pest-free periodic solution (25) is locally stable, and if (29) holds true, then the pest-free periodic solution (25) is globally attractive.

4.2. Bifurcation Analysis and Coexistence. One dimensional bifurcation analysis, which is a traditional approach to gain preliminary insight into the properties of a dynamic system, provides information about the dependence of the dynamics on a certain parameter. The analysis is expected to reveal the type of attractor, to which the dynamics will ultimately settle down after passing the initial transient phase, and within which the trajectory will then remain forever.

From threshold conditions (28) and (29), we know that model (24) has a globally asymptotically stable pest eradication periodic solution if the impulsive periodic T < [T.sup.max.sub.0] [??] min{[T.sup.max.sub.1], [T.sup.max.sub.2]}. Unfortunately, we cannot get its analytical expression for [T.sup.max.sub.1], [T.sub.max.sup.2] due to the complexity of [y.sup.*] in conditions (28) and (29), but it can be easily confirmed by the numerical bifurcation analysis if the impulsive period T exceeds [T.sup.max.sub.0]; both prey and predator population can oscillate periodically (see Figures 2 and 3). Moreover, here we are interested in the biological implications of these results in insect pest control.

Figure 2(a) indicates that if the pulse period T is more than [T.sup.max.sub.0], the pest-free periodic solution becomes unstable and variable x begins to oscillate with a large amplitude that corresponds to periodic outbreaks of the pest population. If the pulse period is further increased, a sequence of period adding bifurcations interchanging with regions of chaos is observed. Figures 2 and 3 show that model (24) has more complex and interesting dynamic behaviors including periodic doubling bifurcation, chaotic solutions, multistability, chaos crisis, periodic adding, periodic windows, periodic halving bifurcation, and crises with the increase of pulse period T and constant stocking number [sigma], respectively. For example, Figures 2(a) and 2(b) show that for model (24) there exists a periodic solution and a quasiperiodic solution, when T = 5 and T = 7.8, respectively; see Figures 4 and 5 for full details.

Meanwhile, bifurcation analyses also indicate that multiple attractors can coexist for a wide range of parameters. For example, the two attractors with quite different pest amplitudes can coexist, when [sigma] = 5 in Figure 3; see Figure 6 for more details. Figures 7(a)-7(d) indicate that one solution can switch to another attractor with small amplitude at a random time when small changes are introduced in initial densities of pest and natural enemies populations. That is, the proper initial densities of pest and natural enemies populations can affect the outcome of classical biological control, and the final stable states of pest and natural enemies populations depend on their initial densities, which can help us to design control strategies and to make management decisions. These results are further confirmed by basins of attraction of initial densities, as shown in Figure 8.

5. Discussion and Biological Conclusions

In this contribution we focus on a generalized predator-prey model under limited resource; the main purpose of this paper is to understand the effect of resource limitation on outbreaks of a pest population. It is worth noting that the mathematical model described here differs from the previous ones since we have taken the nonlinear impulsive function into account. The nonlinear perturbation of the impulsive functions makes the dynamical behavior of solutions change dramatically and it gets more and more complicated, which means that it is a very difficult task to control pests under resource limitation. To show applications of the main results obtained in present work and study the effects of parameter space on the threshold conditions, we employ the Holling II functional response function as an example. By choosing impulsive perturbations [sigma] and T as bifurcation parameters, we have obtained bifurcation diagrams Figures 2 and 3 for model (24).

The results demonstrate that the dynamic behavior of model (24) may be dramatically affected by small changes in the value of initial densities of pest and natural enemy with resource limitation. Bifurcation diagrams shown in Figures 2 and 3 clarify that model (24) has several different attractors which can coexist for a wide range of parameters combined with periodic doubling bifurcation, chaotic solutions, quasiperiodic solutions (see Figure 5 for details), tangent bifurcation, multistability (see Figure 6 for details), chaos crisis, periodic windows, periodic halving bifurcation, crises, and so on, which can help us to further understand the applications of nonlinear pulses in our model (for more information about bifurcation diagrams, refer to [13, 15, 36, 37]). It is worth noting that there are several hidden factors that can adversely affect our control strategy with resource limitation; that is, the increasing number of identifying complicated, possibly chaotic, dynamics in IPM models may present a major challenge for controlling the pest populations in practice.

The results presented in previous section also indicate that the successful biological control depends on the initial densities of the pest and natural enemies populations. The proper initial densities of pest and natural enemies populations can affect the outcome of pest populations and can help us to design control strategies and to make management decisions. Figures 6 and 8 show that the pest and natural enemy populations can coexist in various pest and natural enemies initial densities.

Appendices

A. Proof of Lemma 2

Taking derivative of g(z) with respect to z, one yields

g'(z) = ([a.sub.1]z + [a.sub.1] + [a.sub.2]) [e.sup.z] + [a.sub.3] = 0. (A.1)

Rearranging the above equation we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.2)

According to Definition 1, solving (A.2) with respect to z, one yields two roots [z.sup.*.sub.1] and [z.sup.*.sub.2]; that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.3)

Further, both [z.sup.*.sub.1] and [z.sup.*.sub.2] are well defined provided

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.4)

That is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.5)

Obviously, [z.sup.*.sub.1] > 0 provided

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.6)

which is equivalent to

[a.sub.1] + ([a.sub.2]) + ([a.sub.3]) > 0, ([a.sub.1]) + ([a.sub.2]) > 0. (A.7)

That is, [a.sub.1] + [a.sub.2] + [a.sub.3] > 0.

Similarly, [z.sup.*.sub.2] > 0 provided

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (A.8)

which is equivalent to

[a.sub.1] + ([a.sub.2]) + ([a.sub.3]) < 0, ([2a.sub.1]) + ([a.sub.2]) > 0. (A.9)

Since [a.sub.1] < [a.sub.3], the above two inequalities cannot be true simultaneously. This indicates that the function g(z) exists an unique extreme value just at [z.sup.*.sub.1].

Furthermore, it follows from (6) that

g(0) = [a.sub.2] > 0, g (z) [right arrow] -[infinity] as z [right arrow] +[infinity] (A.10)

Therefore, the function g(z) = 0 always exists as a unique root if the conditions of Lemma 2 are true. The proof is completed.

B. Proof of Theorem 4

To proof the local stability of the solution ([x.sub.p](t), [y.sub.p](t)) of (3), we need to investigate the difference equation determined by the impulsive period T with respect to the linear equation (20). Let [PHI](T) be the fundamental matrix of (20); thus [PHI](T) must satisfy

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.1)

where [PHI](0) = I is the identity matrix and the term * is not necessarily computed in detail as it is not required in the following analysis.

Resetting the impulsive conditions of model (3) from third and fourth equations becomes

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.2)

Hence, according to Floquet theory [38], if the module of both eigenvalues of the matrix

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.3)

is less than one, then the periodic solution ([x.sub.p](t), [y.sub.p](t)) is locally stable. In fact, two Floquet multiplies are thus

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.4)

The pest-free solution ([x.sub.p](t), [y.sub.p](t)) is locally stable if [absolute value of [[lambda].sub.1]] < 1, [absolute value of [[lambda].sub.2]] < 1; that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.5)

and by simple calculation, we have that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.6)

holds true. The pest-free solution ([x.sub.p](t), [y.sub.p](t)) is locally stable if [absolute value of [[lambda].sub.1]] < 1, which holds true due to [R.sup.1] < 1.

In order to prove the global attractivity, we should prove that [bar.x] tends to zero, and it follows that [bar.y] tends to zero; that is, y(t) approximates to [y.sub.p](t).

It follows from [??] [greater than or equal to] -[delta][bar.y] and impulsive conditions that for any initial condition ([[bar.x].sub.0], [[bar.y].sub.0]) we have [bar.y] [greater than or equal to] min(0, [bar.[y.sub.0]])[e.sup.-[delta]t] [??] [[bar.y].sub.m](t).

In order to investigate the first equation of (20), that is, [??] equation, we define function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.7)

Since [phi](x) > 0, it is easy to see that the function G([bar.x]) can be seen to be an increasing function of [bar.x] from [bar.x] = 0. Thus, if we can prove that G([bar.x]) [right arrow] -[infinity] as t [right arrow] [infinity], then the pest population will die out eventually. From the definition of G([bar.x]) function we have

[[dG([bar.x])]/dt] = ([1/[[phi](bar.x)]])[??] = [[f([bar.x])]/[phi](bar.x)] - [bar.y](t) - [y.sub.p](t) [less than or equal to] ([[f(bar.x)]/[phi]([bar.x])]) - [[bar.y].sub.m](t) - [y.sub.p](t). (B.8)

Based on the definition of G([bar.x]), we integrate G([bar.x]) between two successive release predator interval, that is, the time interval (nT, (n+1)T] for all n [member of] N; that is,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.9)

For any t there exists an integer l such that t [member of] (lT, (l + 1)T], so according to the periodicity of [y.sub.p](s) and for all t > 0 we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (B.10)

The first and second terms of the right-hand side are upper bounded due to the periodicity of [y.sub.p](t) with periodic [T.sub.p]. Note that l [right arrow] [infinity] as t [right arrow] [infinity]. Therefore, if

([[integral].sup.T.sub.0])[[M.sub.s] - [y.sub.p](s)]ds + [[ln (1 - [p.sub.1])]/[m.sub.s]] < 0 (B.11)

hold true, then we have that G([bar.x]) [right arrow] [infinity] as t [right arrow] [infinity]. That is,

[T[M.sub.s]([m.sub.s])/[[m.sub.s]([[integral].sup.T.sub.0])[y.sub.p](s)ds - ln(1 - [p.sub.1])]] [less than or equal to] 1. (B.12)

So [bar.x](t) tends to zero as t [right arrow] [infinity] provided (B.12).

Now we prove that [bar.y](t) [right arrow] 0 as well. Since [bar.x] [right arrow] 0, so there is a finite time [t.sub.s] such that [phi]([bar.x]) [less than or equal to] [delta]/2. Therefore we have

[??] = [mu]([phi])([bar.x])[[bar.y] + [y.sub.p](t)] - [beta]([bar.y]) [less than or equal to] [mu][phi]([bar.x]) [y.sub.p](t) - [[delta]/2] (B.13)

for t > [t.sub.s]. It follows from [bar.x] [right arrow] 0 as t [right arrow] [infinity] and the periodicity of [y.sub.p](t) that we have [phi]([bar.x])[y.sub.p](t) [right arrow] 0 as t [right arrow] [infinity]. Consequently, [bar.y] [right arrow] 0 as t [right arrow] [infinity]. Those indicate that if the inequality of (B.12) holds, then the pest-free periodic solution ([x.sub.p](t), [y.sub.p](t)) is globally attractive, and in combination with local stability we can conclude the periodic solution ([x.sub.p](t), [y.sub.p](t)) of model (3).

http://dx.doi.org/10.1155/2014/919242

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (NSFCs, 11171199, 11371030, and 11301320) and the Fundamental Research Funds for the Central Universities (GK201305010). The authors would like to thank the editor and the referees for their careful reading of the original paper and many valuable comments and suggestions that greatly improved the presentation of this paper.

References

[1] M. L. Flint, "Integrated pest management for walnuts," in University of California Statewide Integrated Pest Management Project, vol. 3270, p. 3641, Division of Agriculture and Natural Resources, University of California, Oakland, Calif, USA, 2nd edition, 1987.

[2] J. C. Van Lenteren, "Integrated pest management in protected crops," in Integrated Pest Management, pp. 311-320, Chapman & Hall, London, UK, 1995.

[3] J. C. Van Lenteren and J. Woets, "Biological and integrated pest control in greenhouses," Annual Review of Entomology, vol. 33, no. 4, pp. 239-250, 1988.

[4] J. C. Van Lenteren, "Measures of success in biological control of arthropods by augmentation of natural enemies," in Measures of Success in Biological Control, S. Wratten and G. Gurr, Eds., pp. 77-89, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2000.

[5] J. C. Van Lenteren, "Environmental manipulation advantageous to natural enemies of pests," in Integrated Pest Management, V. Delucchi, Ed., pp. 123-166, Parasitis, Geneva, Switzerland, 1987.

[6] Y. N. Xiao and F. Van den Bosch, "The dynamics of an eco-epidemic mode- with biological control," Ecological Modelling, vol. 168, no. 1-2, pp. 203-214, 2003.

[7] H. J. Barclay, "Models for pest control using predator release, habitat management and pesticide release in combination," Journal of Applied Ecology, vol. 19, no. 2, pp. 337-348, 1982.

[8] H. J. Barclay and P. Van den Diressche, "Predator-prey models with added mortality," Canadian Entomologist, vol. 109, no. 5, pp. 763-768, 1977.

[9] N. D. Barlow, H. Moller, and J. R. Beggs, "A model for the effect of Sphecophaga vesparum as a biological control agent of the common wasp in New Zealand," Journal of Applied Ecology, vol. 33, no. 1, pp. 31-44, 1996.

[10] S. Y. Tang and L. S. Chen, "Modelling and analysis of integrated pest management strategy," Discrete and Continuous Dynamical Systems. Series B, vol. 4, no. 3, pp. 759-768, 2004.

[11] S. Y. Tang, Y. N. Xiao, L. S. Chen, and R. A. Cheke, "Integrated pest management models and their dynamical behaviour," Bulletin of Mathematical Biology, vol. 67, no. 1, pp. 115-135, 2005.

[12] S. Y. Tang, Y. N. Xiao, and D. Clancy, "New modelling approach concerning integrated disease control and cost-effectivity," Nonlinear Analysis: Theory, Methods & Applications, vol. 63, no. 3, pp. 439-471, 2005.

[13] S. Y. Tang, Y. N. Xiao, and R. A. Cheke, "Multiple attractors of host-parasitoid models with integrated pest management strategies: eradication, persistence and outbreak," Theoretical Population Biology, vol. 73, no. 2, pp. 181-197, 2008.

[14] S. Y. Tang, G. Y. Tang, and R. A. Cheke, "Optimum timing for integrated pest management: modelling rates of pesticide application and natural enemy releases," Journal of Theoretical Biology, vol. 264, no. 2, pp. 623-638, 2010.

[15] S. Y. Tang, J. H. Liang, Y. S. Tan, and R. A. Cheke, "Threshold conditions for integrated pest management models with pesticides that have residual effects," Journal of Mathematical Biology, vol. 66, no. 1-2, pp. 1-35, 2013.

[16] J. H. Liang, S. Y. Tang, and R. A. Cheke, "An integrated pest management model with delayed responses to pesticide applications and its threshold dynamics," Nonlinear Analysis: Real World Applications, vol. 13, no. 5, pp. 2352-2374, 2012.

[17] J. H. Liang, S. Y. Tang, R. A. Cheke, and J. H. Wu, "Adaptive release of natural enemies in a pest-natural enemy system with pesticide resistance," Bulletin of Mathematical Biology, vol. 75, no. 11, pp. 2167-2195, 2013.

[18] J. H. Liang, S. Y. Tang, J. J. Nieto, and R. A. Cheke, "Analytical methods for detecting pesticide switches with evolution of pesticide resistance," Mathematical Biosciences, vol. 245, no. 2, pp. 249-257, 2013.

[19] P. Georgescu, H. Zhang, and L. S. Chen, "Bifurcation of nontrivial periodic solutions for an impulsively controlled pest management model," Applied Mathematics and Computation, vol. 202, no. 2, pp. 675-687, 2008.

[20] B. Liu, L. S. Chen, and Y. J. Zhang, "The dynamics of a prey-dependent consumption model concerning impulsive control strategy," Applied Mathematics and Computation, vol. 169, no. 1, pp. 305-320, 2005.

[21] N. Baharanyi, N. O. Tackie, A. Pierce, C. Woolery, and R. Hopkinson, "The impact of an integrated marketing and educational strategy mix on the Montgomery State Farmers Market-Phase I," Prepared for the College of Agricultural, Environmental and Natural Sciences, Tuskegee University, Tuskegee, AL and the Alabama Department of Agriculture and Industries, Montgomery, Ala, USA, 1993.

[22] N. O. Tackie, H. J. Findlay, and N. Baharanyi, "Farm products marketing practices by limited resource farmers," Journal of Agribusiness, vol. 16, no. 1, pp. 43-51, 1998.

[23] N. O. Tackie, J. E. N. Jackai, R. Ankumah, R. Noble, and M. Collins, "Perceptions of small farmers on marketing farm products," in Applied Economics and Related Sciences, Number 0604-01, George Washington Carver Agricultural Experiment Station, Tuskegee University, Tuskegee, Ala, USA, 2004.

[24] N. O. Tackie, L. E. N. Jackai, R. Ankumah, B. N. Dingha, A. W. Salifu, and O. Ojumu, "Integrated pest management and protection practices by limited resource farmers," Journal of Extension, vol. 47, no. 1, Article ID 1RIB6, 2009.

[25] M. J. Reynolds-Hogland, L. B. Pacifici, and M. S. Mitchell, "Linking resources with demography to understand resource limitation for bears," Journal of Applied Ecology, vol. 44, no. 6, pp. 1166-1175, 2007.

[26] A. J. Lotka, "Undamped oscillations derived from the law of mass action," Journal of the American Chemical Society, vol. 42, no. 8, pp. 1595-1599, 1920.

[27] V. Volterra, "Variations and fluctuations of a number of individuals in animal species living together," in Animal Ecology, R. N. Chapman, Ed., pp. 440-448, McGraw-Hill, New York, NY, USA, 1931.

[28] Y. Ito, "Representation of functions by superpositions of a step or sigmoid function and their applications to neural network theory," Neural Networks, vol. 4, no. 3, pp. 385-394, 1991.

[29] L. E. Yelle, "The learning curve: historical review and comprehensive survey," Decision Sciences, vol. 10, no. 2, pp. 302-328, 1979.

[30] L. Mailleret and F. Grognard, "Global stability and optimisation of a general impulsive biological control model," Mathematical Biosciences, vol. 221, no. 2, pp. 91-100, 2009.

[31] S. Nundloll, L. Mailleret, and F. Grognard, "The effect of partial crop harvest on biological pest control," Rocky Mountain Journal of Mathematics, vol. 38, no. 5, pp. 1633-1661, 2008.

[32] S. Nundloll, L. Mailleret, and F. Grognard, "Two models of interfering predators in impulsive biological control," Journal of Biological Dynamics, vol. 4, no. 1, pp. 102-114, 2010.

[33] S. Nundloll, L. Mailleret, and F. Grognard, "Influence of Intrapredatory Interferences on Impulsive Biological Control Efficiency," Bulletin of Mathematical Biology, vol. 72, no. 8, pp. 2113-2138, 2010.

[34] R. M. Corless, G. H. Gonnet, D. E. G. Hare, and D. J. Jeffrey, "On the Lambert W function," Advances in Computational Mathematics, vol. 5, no. 1, pp. 329-359, 1996.

[35] J. Waldvogel, "The period in the Volterra-Lotka predator-prey modle," SIAM Journal on Numerical Analysis, vol. 20, no. 6, pp. 1264-1272, 1983.

[36] S. Y. Tang and L. S. Chen, "Multiple attractors in stage-structured population models with birth pulses," Bulletin of Mathematical Biology, vol. 65, no. 3, pp. 479-495, 2003.

[37] Y. N. Xiao and S. Y. Tang, "The effect of initial density and parasitoid intergenerational survival rate on classical biological control," Chaos, Solitons & Fractals, vol. 37, no. 4, pp. 1048-1058, 2008.

[38] A. Lakmeche and O. Arino, "Bifurcation of non trivial periodic solutions of impulsive differential equations arising chemotherapeutic treatment," Dynamics of Continuous, Discrete and Impulsive Systems, vol. 7, no. 2, pp. 265-287, 2000.

Wenjie Qin, (1) Guangyao Tang, (1,2) and Sanyi Tang (1)

(1) College of Mathematics and Information Science, Shaanxi Normal University, Xi'an 710062, China

(2) School of Science, Hubei Minzu University, Enshi 445000, China

Correspondence should be addressed to Wenjie Qin; wenjieqin@hotmail.com

Received 8 October 2013; Accepted 16 December 2013; Published 9 February 2014

Academic Editor: XianHua Tang

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Qin, Wenjie; Tang, Guangyao; Tang, Sanyi |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5909 |

Previous Article: | A novel method for target navigation and mapping based on laser ranging and MEMS/GPS navigation. |

Next Article: | Modified nonradial supper efficiency models. |

Topics: |