Top Banner
Josh Durham Jacob Swett Dynamic Behaviors of a Harvesting Leslie-Gower Predator-Prey Model
27

Dynamic Behaviors of a Harvesting Leslie-Gower Predator-Prey Model

Feb 23, 2016

Download

Documents

Thina

Dynamic Behaviors of a Harvesting Leslie-Gower Predator-Prey Model. Josh Durham Jacob Swett. The Article. Dynamic Behaviors of a Harvesting Leslie-Gower Predator-Prey Model Na Zhang, 1 Fengde Chen, 1 Qianqian Su, 1 and Ting Wu 2 - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Josh DurhamJacob Swett

Dynamic Behaviors of a Harvesting Leslie-Gower

Predator-Prey Model

Page 2: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Dynamic Behaviors of a Harvesting Leslie-Gower Predator-Prey Model

Na Zhang,1 Fengde Chen,1 Qianqian Su,1 and Ting Wu2

1 College of Mathematics and Computer Science, Fuzhou University, Fuzhou 350002, Fujian, China

2 Department of Mathematics and Physics, Minjiang College, Fuzhou 350108, Fujian, China

Published in Discrete Dynamics in Nature and SocietyReceived 24 October 2010; Accepted 8 February 2011Academic Editor: Prasanta K. Panigrahi

The Article

Page 3: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

IntroductionStability Property of Positive EquilibriumThe Influence of HarvestingBionomic EquilibriumOptimal Harvesting PolicyNumerical ExampleConclusionsQuestions

Outline

Page 4: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

– density of prey species – density or predator species - intrinsic growth rate for prey and predator,

respectively - catch rate - competition rate - prey conversion rate is the carrying capacity of the prey is the prey-dependent carrying capacity of the

predator

Predatory Prey Model

Page 5: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

- constant effort spent by harvesting agency on the prey

- constant effort spent by harvesting agency on the predator

Assume:

Model with Harvesting

Page 6: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

The harvesting of biological resources commonly occurs in:FisheriesForestryWildlife management

Allows for predictions given various assumptions

Important for optimization

Why Study a Harvesting Model?

Page 7: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Given that the harvesting remains strictly less than the intrinsic growth rate, the system under study has a unique positive equilibrium

Satisfies the following equalities,

Stability Property of Positive Equilibrium

𝐻1∗ = (𝑟1−𝑐1)𝑎2𝑎1(𝑟2−𝑟1+𝑎2𝑏1,

𝑃1∗= (𝑟1−𝑐1)(𝑟2−𝑐2)𝑎1ሺ𝑟2−𝑐2ሻ+𝑎2𝑏1

Page 8: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

It means that the positive equilibrium of the system under study is locally asymptotically stableStability is the same as what we have

discussed in classThe proof of this is very similar to examples

that we have done in class.First, find the Jacobian

What does this mean?

Page 9: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Then we find the characteristic equation

Given the above information, we can see that the unique positive equilibrium of the system is stable

Proof Continued

𝜆2 + 𝑎𝜆+ 𝑏= 0, 𝑤ℎ𝑒𝑟𝑒 𝑎 = 𝑎2𝑏1ሺ𝑟1 − 𝑐1ሻ𝑎1ሺ𝑟2 − 𝑐2 + 𝑎2𝑏1ሻ+ 𝑟2 − 𝑐2ሺ> 0ሻ

𝑎𝑛𝑑 𝑏= ሺ𝑟1 − 𝑐1ሻሺ𝑟2 − 𝑐2ሻ(> 0)

Page 10: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

The positive equilibrium is globally stableThe proof of this fact is beyond the scope of

this courseInstead, there will be a brief summary of major

pointsEquilibrium dependent only on coefficients of

systemLyapunov FunctionLyapunov’s asymptotic stability theorem

Global Stability

Page 11: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Case 1: Harvesting only prey species

Case 2: Harvesting only predator species

Influence of Harvesting𝑑𝐻1∗𝑑𝑐1 = −𝑎2𝑎1𝑟2 + 𝑎2𝑏1 < 0 𝑎𝑛𝑑

𝑑𝑃1∗𝑑𝑐1 = −𝑟2𝑎1𝑟2 + 𝑎2𝑏1 < 0

𝑑𝐻1∗𝑑𝑐2 = 𝑟1𝑎2𝑎1(𝑎1൫𝑟2 − 𝑐2) + 𝑎2𝑏1൯2 > 0 𝑎𝑛𝑑 𝑑𝑃1∗𝑑𝑐2 = −𝑟1𝑎2𝑏1(𝑎1(𝑟2 − 𝑐2) + 𝑎2𝑏1)2 < 0

Page 12: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Case 3: Harvesting predator and prey togetherDifficult to give a detailed analysis of all

possible cases so the focus will be on answeringWhether or not it is possible to choose harvesting

parameters such that the harvesting of predator and prey will not cause a change in density of the prey species over time

If it is possible, what will the dynamic behaviors of the predator species be?

We find that:

allows the first question to be answered with, yes

We also find that by substituting the above equality into the predator equation, it will lead to a decrease in the predator species

𝑐2 = ൬𝑎1𝑟2 + 𝑎2𝑏1𝑟1𝑎1 ൰𝑐1

Page 13: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Biological equilibrium + Economic equilibrium =Bionomic equilibrium

Biological equilibrium:

Economic equilibrium:TR = TC

TR is the total revenue obtained by selling the harvested predators and prey

TC is the total cost for the effort of harvesting both predators and prey

Bionomic Equilibrium

Page 14: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

We will define four new variables:p1 is the price per unit biomass of the prey Hp2 is the price per unit biomass of the predator

Pq1 is the fishing cost per unit effort of the prey

Hq2 is the fishing cost per unit effort of the

predator P

The revenue from harvesting can written as:

Where: and And: and

Bionomic Equilibrium (Cont.)

Page 15: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

The revenue from harvesting equation and the predator and prey equations must all be considered together:

Price per unit biomass (p1,2) and the fishing cost per unit effort (q1,2) are assumed to be constant

Since the total revenue (TR) and total cost (TC) are not determined, four cases will be considered to determine bionomic equilibrium

Bionomic Equilibrium (Cont.)

Page 16: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

That is: In other words the revenue () is less than the cost () for

harvesting prey and it will be stopped (i.e. c1 = 0) Thus, Predator harvesting will continue as long as

To determine equilibrium: Solve for P Substitute it into the predator equation Substitute the predator equation and P into the prey equation Simplify

Thus, if r1 > a2(q2/p2) and r2 > (a2b1q2/r1q2) both hold, then bionomic equilibrium is obtained.

Case I

Page 17: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

That is: In other words the revenue () is less than the cost () for harvesting

predators and it will be stopped (i.e. c2 = 0) Thus, Prey harvesting will continue as long as

Example:

Sub H into predator equation , Yields: Substituting H and P into the prey equation, Yields:

Thus, if r1 > ((a1r2 – a2b1)q1/a2p1) holds then bionomic equilibrium is obtained for this case.

Case II

Page 18: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

and That is to say, the total costs exceeds the

revenue for both predators and preyNo profit

Clearly c1 = c2 = 0No bionomic equilibrium

Case III

Page 19: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

and That is: and Solve as before

Thus if and hold then bionomic equilibrium is obtained.

It becomes obvious that bionomic equilibrium may occur if the intrinsic growth rates of the predators and prey exceed the values calculated

Case IV

Page 20: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

To determine an optimal harvesting policy, a continuous time-stream of revenue function, J, is maximized:

-δ is the instantaneous annual rate of discount

c1(t) and c2(t) are the control variablesThe assumption: ; still holds

Pontryagin's Maximum Principle is invoked to maximize the equation

Optimal Harvesting Policy

Page 21: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

A method for the computation of optimal controls

The Maximum Principle can be thought of as a far reaching generalization of the classical subject of the calculus of variations

Pontryagin's Maximum Principle

Page 22: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Applying Pontryagin's Maximum Principle to the revenue function J, shows that optimal equilibrium effort levels (c1 and c2) are obtained when:

Recall:- intrinsic growth rate for prey and predator,

respectively - catch rate - competition rate - prey conversion rate

Results and Implications

Page 23: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Using the following values as inputs into the optimized equation:

Gives:

Solving with Maple, the authors obtained:

Numerical Example

Page 24: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

From the previous results only one results meets the following conditions:

Namely:

These values can then be entered into the predator-prey harvesting equation

Numerical Example (Cont.)

Page 25: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Substituting the values for H and P into the below equations:

And rearranging to solve for c1 and c2, gives:

Numerical Example (Cont.)

Page 26: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Introduced a harvesting Leslie-Gower predator-prey model

The system discussed was globally stableProvided an analysis of some effects of

different harvesting policiesConsidered economic profit of harvestingResults show that optimal harvesting policies

may existDemonstrate that the optimal harvesting

policy is attainable

Conclusions

Page 27: Dynamic Behaviors of a Harvesting  Leslie-Gower Predator-Prey  Model

Questions?