On predator-prey models
|
|
- Annice Hart
- 5 years ago
- Views:
Transcription
1 Predator-prey On models ddd Department of Mathematics College of William and Mary Math 41/CSUMS Talk February 3, 1
2 Collaborators Sze-Bi Hsu (Tsinghua University, Hsinchu, Taiwan) Junjie Wei (Harbin Institute of Technology, Harbin/Weihai, China) Rui Peng (University of New England, Armidale, Australia) Fengqi Yi (Harbin Engineering University, Harbin, China) Ph.D. 8 Jinfeng Wang (Harbin Institute of Technology/Harbin Normal University, Harbin, China) Ph.D. expected 1 Yuanyuan Liu (William & Mary, Williamsburg, USA) B.S. expected (math and econ) 1
3 Fur trade of Hudson Bay Company ( )
4 Hudson Bay Company lynx-hare data Charles Elton (194), Periodic fluctuations In the numbers of animals: their causes and effects, British Journal of Experimental Biology, was first (of MANY) publications to analyze this data set Rebecca Tyson, et.al. (9), Modelling the Canada lynx and snowshoe hare population cycle: the role of specialist predators, Theoretical Ecology, one of the latest.
5 Lotka-Volterra Predator-Prey Model Predator-prey Alfred Lotka ( ) Vito Volterra ( ) du dt = au buv, dv dt = cv + duv. see Math 3 Chapter
6 Predator-prey system with functional response du = u(a bu) cφ(u)v, dt dv = dv + f φ(u)v. dt φ(u): predator functional response φ(u) = u (Lotka-Volterra) u φ(u) = (Holling type II, m: prey handling ) 1 + mu [Holling, 1959](Michaelis-Menton biochemical kinetics) Biological work: [Rosenzweig-MacArthur, American Naturalist 1963] [Rosenzweig, Science, 1971] (Paradox of enrichment) [May, Science, 197] (Existence/uniqueness of limit cycle)
7 Rosenzweig-MacArthur model du (1 dt = u u ) muv K Nullcline: u =, v = 1 + u, dv muv = dv + dt 1 + u (K u)(1 + u) ; v =, d = mu m 1 + u. Solving d = mu d, one have u = λ 1 + u m d. Equilibria: (,), (1,), (λ,v λ ) where v λ = We take λ as a bifurcation parameter (K λ)(1 + λ) m Case 1: λ K: (K, ) is globally asymptotically stable Case : (K 1)/ < λ < K: (K,) is a saddle, and (λ,v λ ) is a locally stable equilibrium Case 3: < λ < (K 1)/: (K,) is a saddle, and (λ,v λ ) is a locally unstable equilibrium (λ = (K 1)/ is a Hopf bifurcation point)
8 Phase Portraits Left: (K 1)/ < λ < K: (K,) is a saddle, and (λ,v λ ) is a locally stable equilibrium Right: < λ < (K 1)/: (K,) is a saddle, and (λ,v λ ) is a locally unstable equilibrium; there exists a limit cycle A supercritical Hopf bifurcation occurs.
9 Global stability [Hsu, Hubble, Waltman, SIAM J. Appl. Math., 1978] [Hsu, Math. Biosci., 1978] (λ,v λ ) is globally asymptotically stable if K 1, or K > 1 and (K 1)/ < λ < K. [Cheng, SIAM J. Math. Anal., 1981] If < λ < (K 1)/, then (λ,v λ ) is unstable, and there is a unique periodic orbit which is globally asymptotically stable. More on uniqueness of limit cycle: [Zhang, 1986], [Kuang-Freedman, 1988] [Hsu-Hwang, 1995,1998], [Xiao-Zhang, 3, 8]
10 Summary of ODE du (1 dt = u u ) muv K Nullcline: u =, v = 1 + u, dv muv = dv + dt 1 + u (K u)(1 + u) ; v =, d = mu m 1 + u. Solving d = mu d, one have u = λ 1 + u m d. Equilibria: (,), (K,), (λ,v λ ) where v λ = We take λ as a bifurcation parameter (K λ)(1 + λ) m Case 1: λ K: (K, ) is globally asymptotically stable Case : (K 1)/ < λ < K: (λ,v λ ) is globally asymptotically stable Case 3: < λ < (K 1)/: unique limit cycle is globally asymptotically stable (λ = (K 1)/: Hopf bifurcation point)
11 Phase portrait (1) du dt = u (1 u K ) muv 1 + u, dv dt = dv + muv 1 + u Case 1: λ K: (K, ) is globally asymptotically stable v=f(u) d=g(u) v u (,) (k,)
12 Phase portrait () du dt = u (1 u K ) muv 1 + u, dv dt = dv + muv 1 + u Case : (K 1)/ < λ < K: (K,) is a saddle, and (λ,v λ ) is a locally stable equilibrium v=f(u) d=g(u) v (λ,f(λ)) u (,) (k,)
13 Phase portrait (3) du dt = u (1 u K ) muv 1 + u, dv dt = dv + muv 1 + u Case 3: < λ < (K 1)/: (K,) is a saddle, and (λ,v λ ) is a locally unstable equilibrium d=g(u) (λ,f(λ)) v v=f(u) u (,) (k,)
14 v A new result of this ODE [Hsu-Shi, Disc. Cont. Dyna. Syst.-B, 9] Relaxation oscillator profile of limit cycle in system. (Motivated by numerical observation) du dt = u (1 u) muv a + u, dv dt = dv + muv a + u u
15 Graph of limit cycle Parameters: a =.5, m = 1, d =.1, λ = 1/18.56, period T u and v t
16 v Small d u
17 Graph of limit cycle Parameters: a =.5, m = 1, d =.1, λ = 1/198.5, period T u and v t
18 Illustration of limit cycle v O 1 O 4 v 5 (u) v 4 (u) a/m O 5 v 6 (u) O O 3
19 Relaxation oscillation Theorem[Hsu-Shi, 9] If < a < 1 and m > are fixed, and as d (thus λ ), then C 1 λ 1 T(O 1 O ) C λ 1, T(O O 3 ) = O( lnλ ), T(O 4 O 1 ) = O( lnλ ), and T(O 3 O 4 ) = O(1). In particular, the period T as d. The shape of the graph of the limit cycle is a relaxation oscillator. Other known relaxation oscillators: Van der Pol oscillator in electrical circuits FitzHugh-Nagumo oscillator in action potentials of neurons Many other physiology models: heart beat, calcium signaling [Liu-Xiao-Yi, JDE, 3] oscillations in singularly perturbed 3-D system [Zhang-Wang-Wang-Shi, preprint] extend it to more general system
20 Relaxation oscillation (left) Van der Pol oscillator; (right) FitzHugh-Nagumo oscillator.
21 More general system du dt = g(u)(f (u) v), dv dt = v (g(u) d). (a1) f C 3 (R + ), f () >, there exists K >, such that for any u >, u K, f (u)(u K) < and f (K) = ; there exists λ (,K) such that f (u) > on [, λ), f (u) < on ( λ,k]; (a) g C (R + ), g() = ; g(u) > for u > and g (u) > for u ; there exists a unique λ (,K) such that g(λ) = d. g(u): functional response, f (u): prey isocline Rosenzweig-MacArthur model: g(u) = mu f (u) = (K u)(1 + u) m 1 + u and
22 Phase portraits du dt = g(u)(f (u) v), dv dt = v (g(u) d). v=f(u) d=g(u) d=g(u) v (λ,f(λ)) v (λ,f(λ)) v=f(u) (,) (k,) u (,) (k,) u Left: (λ,f (λ)) stable; Right: (λ,f (λ)) unstable, exist limit cycle
23 History du dt = g(u)(f (u) v), dv dt = v (g(u) d). [Hsu, 1979] If λ > λ and f (u) is concave in [,K], then (λ,f (λ)) is globally stable. ( λ is the top of the hump) [Kuang-Freedman, ( 1988] If < λ < λ, and d f ) (u)g(u) for all x [,K], then the limit cycle du g(u) d is unique and globally stable. [Hofbauer-So, 1989] counterexample: f (u) is concave, but Hopf bifurcation is subcritical, so there are two periodic orbits for λ ( λ, λ + ǫ), one of them is locally stable. More work: [Ruan-Xiao, 1] [Xiao-Zhang, 3] [Liu, 5]
24 Our result du dt = g(u)(f (u) v), dv dt = v (g(u) d). Theorem[Wang-Shi-Wei, submitted] If λ > λ, and (a8) (uf (u)), (u/g(u)) for u [,K], and (uf (u)) for u ( λ,k); or (a9) f (u), (1/g(u)) for u [,K], and f (u) for u ( λ,k), then (λ,f (λ)) is globally stable. Sharpness: If (uf (u)) u= λ > and g(u) = u/(a + u) or g(u) = u, then Hopf bifurcation is subcritical and global stability does not hold. Realistic example: du (1 dt = ru u ) ( 1 A + C ) dv Buv, = dv + Buv. K u + C dt Prey growth: weak Allee effect when A < and C > A. C large: global stability; C small: subcritical Hopf bifurcation.
25 Global stability vs. Multiple periodic orbits du (1 dt = ru u ) ( 1 A + C ) Buv, K u + C dv dt = dv + Buv. v 1. v u u Both: r = B = 1, A =.8, K = 1; (left) d =.1199, C =.5; (subcritical Hopf) (right) d =.6, C = (supercritical Hopf).
26 Habitat for predator and prey Habitat fragmentation
27 Habitat Restoration: connecting patches Crossing bridges built in Banff, Canada; (for elks) under-cross built in Florida, USA (for deers)
28 Habitat network Habitat with a network structure
29 Habitats as graph theory Network connectivity
30 Effect of network structure on dynamics Matthew Holland, Alan Hastings, (8), Strong effect of dispersal network structure on ecological dynamics, Nature. David Vasseur, Jeremy Fox, (9), Phase-locking and environmental fluctuations generate synchrony in a predator prey community, Nature. dn i dt = N i ( 1 N i K dp i dt = mn ip i 1 + N i ep i + d p ) mn ip i 1 + N i + d n n d ij (N j N i ), j=1 n d ij (P j P i ). j=1 Homogeneous patches, dispersal matrix (d ij ) (d ij = d ji ) single patch: Rosenzweig MacArthur model [Holland-Hastings]: n = 1, random network; cluster solutions, long transient dynamics.
31 Stability in network model Local dynamics: du (1 dt = u u ) muv K 1 + u, dv dt = ev + muv 1 + u Case 1: λ K: (K,) is globally stable Case : (K 1)/ < λ < K: (λ,v λ ) is globally stable Case 3: < λ < (K 1)/: unique limit cycle is globally stable (λ = (K 1)/: Hopf bifurcation point) Theorem[Wang-Shi-Liu, in preparation] If λ K, then (K,) n is globally stable; and if K 1 < λ < K, then (λ,v λ ) n is globally stable. [Li-Shuai (1), JDE] Lyapunov funct. ODE on network Note: same result holds for reaction-diffusion model.
32 Two-patch model where F(u,v) = u u = F(u,v) + a(w u), v = G(u,v) + c(x v), w = F(w,x) a(w u), x = G(w,x) c(x v), ( 1 u ) muv muv, G(u,v) = K 1 + u 1 + u ev. where a,c > are the diffusion rates. Can we completely classify the dynamics?
33 Some Partial Results for Two-patch model [Liu-Shi, in preparation] Let U(t) = u(t) + w(t), V (t) = v(t) + x(t) (sum) and W (t) = u(t) w(t), X(t) = v(t) x(t) (difference) 1. When a, c are large, then W(t), X(t) as t so the system is synchronized.. For < a < a, there exists c = c (a) such that when c < c, there are two additional Hopf bifurcation points λ H, λ+ H, and there exists a non-symmetrical periodic orbit for λ (λ H, λ+ H ). (numerical result shows the non-symmetrical periodic orbit is unique and unstable) 3. For a < a 1, there exists c 1 = c 1 (a) such that when c > c 1, there are two additional bifurcation points λ S, λ+ S such that there exists two non-symmetrical equilibrium points for λ (λ S, λ+ S ). ((a) equilibrium points can be algebraically solved with a complicated formula; (b) for any a and c, there are at most 9 equilibrium points.)
34 Numerical bifurcation diagrams Software: Matlab and MatCont (Govaerts, Kuznetsov) similar to Auto, but works under Matlab Bifurcation/continuation of equilibrium, limit cycles, homoclinic orbits of ODEs BP 3 Max(u) H LPC 1 H.5 H BP H e k = 4, m =.5, a =.6, c =.1, bifurcation parameter < e <.4; symmetric Hopf: e =.3, non-symmetric Hopf: e =.387 and e =.1173
35 Numerical solution Software: Matlab coupled population prey 1 predator 1 prey predator k = 4, m =.5, a =.6, c =.1, bifurcation parameter < e <.4; symmetric Hopf: e =.3, non-symmetric Hopf: e =.387 and e =.1173 Here e =.
36 Numerical solution Software: Matlab patch 1 population patch population coupled population coupled population prey 1 predator 1 prey predator prey 1 predator 1 prey predator total prey population total predator population k = 5, m = 9.96, a =, c =.6, bifurcation parameter < e < 1; symmetric Hopf: e = 6.64, no non-symmetric Hopf Here e = 1 and solution appears to be chaotic
37 Attractor? Chaos? An attractor exists as all positive solutions are bounded. When there are non-symmetrical equilibrium points (EQ) or periodic orbits (PO), there is only one stable state (symmetrical PO) in the attractor, but more than unstable states (symmetric EQ, non-symmetric EQ and/or PO). Hence there exist connecting orbits between the two unstable states, which become a heteroclinic loop. This may imply chaotic behavior on a lower dimensional invariant manifold. Conjecture: Chaos occurs on a zero measure set, but the generic dynamics is eventual synchronization (despite possible long transient dynamics) and converges to symmetric PO.
38 Non-symmetric patches ( heterogeneity) where for i = 1,, ) F i (u,v) = u (1 uki u = F 1 (u,v) + a(w u), v = G 1 (u,v) + c(x v), w = F (w,x) a(w u), x = G (w,x) c(x v), m iuv 1 + u, G i(u,v) = m iuv 1 + u e iv. where a,c > are the diffusion rates. 1. synchronization of coupled oscillators. synchronization of equilibrium and oscillator Goldwyn, E.E. and Hastings, A. (8) When can dispersal synchronize populations? Theoretical Population Biology 73:395-4 Goldwyn, E.E. and Hastings, A. (9) Small heterogeneity has large effects on synchronization of ecological oscillators. Bulletin of Mathematical Biology 71:
39 Non-symmetric patches: numerical results Y.Y. Liu (undergraduate research) x = [.95;.5;.;.5]; k = 3; m = ; e1 =.; e =.9 Left: a =.1; c =.; Right: a = 1; c =. coupled population coupled population patch population patch 1 population 4 prey 1 predator prey predator prey 1 predator 1 prey predator total prey population total predator population coupled population coupled population patch population patch 1 population 4 prey 1 predator prey predator prey 1 predator 1 prey predator total prey population total predator population Larger diffusion rates help to achieve synchronization; smaller diffusion rates may not result in synchronization
40 Non-symmetric patches: numerical results Y.Y. Liu (undergraduate research) x = [.95;.5;.;.5]; k = 3; m = ; e1 = 1.8; e =.45 Left: a =.1; c =.; Right: a = 1; c =. coupled population coupled population patch population patch 1 population 4 prey 1 predator prey predator prey 1 predator 1 prey predator total prey population total predator population coupled population coupled population patch population patch 1 population 4 prey 1 predator prey predator prey 1 predator 1 prey predator total prey population total predator population Connecting a prey-only (predator extinction) system to an oscillatory system could make a stable coexistence
41 Non-symmetric patches: numerical results x = [.95;.5;.;.5]; k = 3; m = ; e1 =.; e =.9 Left: a =.1; c =. The limit cycle is non-symmetrical, and not synchronized. patch 1 population patch population prey 1 predator 1 prey predator coupled population 3 1 prey 1 predator 1 prey predator coupled population total prey population total predator population
42 Future Work (for CSUMS) Predator-prey 1. Further work on -patch model (complete bifurcation diagram). Impact or network structure (small number of patches): 3-patch (linear or triangle) 3. If there is a cycle in the network, is there a periodic solution with population running around the cycle? 4. Computer work: writing better Matlab programs for bifurcation diagram calculation
43 References Sze-Bi Hsu; Relaxation oscillator profile of limit cycle in system. Discrete and Continuous Dynamical Systems B, 11, (9) no. 4, Fengqi Yi, Junjie Wei and, Bifurcation and spatiotemporal patterns in a homogeneous diffusive system. Journal of Differential Equations, 46, (9), no. 5, Rui Peng,, Non-existence of Non-constant Positive Steady States of Two Holling Type-II Predator-prey Systems: Strong Interaction Case. Journal of Differential Equations, 47, (9), no. 3,
44 Predator-prey
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
More informationReaction-Diffusion Models and Bifurcation Theory Lecture 1: Basics and Examples from ODEs
Reaction-Diffusion Models and Bifurcation Theory Lecture 1: Basics and Examples from ODEs Junping Shi College of William and Mary, USA Topic course in Fall 2007 My visits to Taiwan: 2005, 2007, 2009, 2010,
More informationLokta-Volterra predator-prey equation dx = ax bxy dt dy = cx + dxy dt
Periodic solutions A periodic solution is a solution (x(t), y(t)) of dx = f(x, y) dt dy = g(x, y) dt such that x(t + T ) = x(t) and y(t + T ) = y(t) for any t, where T is a fixed number which is a period
More informationPredator prey system with strong Allee effect in prey
J. Math. Biol. (2011 62:291 331 DOI 10.1007/s00285-010-0332-1 Mathematical Biology Predator prey system with strong Allee effect in prey Jinfeng Wang Junping Shi Junjie Wei Received: 12 June 2009 / Revised:
More informationDynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationTHE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL
THE ROSENZWEIG-MACARTHUR PREDATOR-PREY MODEL HAL L. SMITH* SCHOOL OF MATHEMATICAL AND STATISTICAL SCIENCES ARIZONA STATE UNIVERSITY TEMPE, AZ, USA 8587 Abstract. This is intended as lecture notes for nd
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationMath 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction
Math 345 Intro to Math Biology Lecture 20: Mathematical model of Neuron conduction Junping Shi College of William and Mary November 8, 2018 Neuron Neurons Neurons are cells in the brain and other subsystems
More informationMultiple focus and Hopf bifurcations in a predator prey system with nonmonotonic functional response
Multiple focus and Hopf bifurcations in a predator prey system with nonmonotonic functional response Dongmei Xiao Department of Mathematics Shanghai Jiao Tong University, Shanghai 200030, China Email:
More information14 Periodic phenomena in nature and limit cycles
14 Periodic phenomena in nature and limit cycles 14.1 Periodic phenomena in nature As it was discussed while I talked about the Lotka Volterra model, a great deal of natural phenomena show periodic behavior.
More informationModels Involving Interactions between Predator and Prey Populations
Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate
More informationHopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models
Vol.37 (SUComS 6), pp.57-66 http://dx.doi.org/.457/astl.6.37. Hopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models Shuang Guo, Xiuli Li, Jian Yu, Xiu Ren Department of Mathematics, Teacher
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationMath 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry
Math 345 Intro to Math Biology Lecture 19: Models of Molecular Events and Biochemistry Junping Shi College of William and Mary, USA Molecular biology and Biochemical kinetics Molecular biology is one of
More informationNUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD CHAIN WITH LOTKA-VOLTERRA LINEAR FUNCTIONAL RESPONSE
Journal of Sustainability Science and Management Volume 6 Number 1, June 2011: 44-50 ISSN: 1823-8556 Universiti Malaysia Terengganu Publisher NUMERICAL SIMULATION DYNAMICAL MODEL OF THREE-SPECIES FOOD
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
More informationPredator-Prey Model with Ratio-dependent Food
University of Minnesota Duluth Department of Mathematics and Statistics Predator-Prey Model with Ratio-dependent Food Processing Response Advisor: Harlan Stech Jana Hurkova June 2013 Table of Contents
More informationMULTIPLE FOCUS AND HOPF BIFURCATIONS IN A PREDATOR-PREY SYSTEM WITH NONMONOTONIC FUNCTIONAL RESPONSE
SIAM J. APPL. MATH. Vol. 66, No. 3, pp. 802 89 c 2006 Society for Industrial and Applied Mathematics MULTIPLE FOCUS AND HOPF BIFURCATIONS IN A PREDATOR-PREY SYSTEM WITH NONMONOTONIC FUNCTIONAL RESPONSE
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationCOMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS
International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL
More informationANALYSIS OF A MODEL OF TWO PARALLEL FOOD CHAINS. Sze-Bi Hsu. Christopher A. Klausmeier. Chiu-Ju Lin
DISCRETE AND CONTINUOUS doi:10.3934/dcdsb.2009.12.337 DYNAMICAL SYSTEMS SERIES B Volume 12, Number 2, September 2009 pp. 337 359 ANALYSIS OF A MODEL OF TWO PARALLEL FOOD CHAINS Sze-Bi Hsu Department of
More informationBifurcation in a Discrete Two Patch Logistic Metapopulation Model
Bifurcation in a Discrete Two Patch Logistic Metapopulation Model LI XU Tianjin University of Commerce Department of Statistics Jinba road, Benchen,Tianjin CHINA beifang xl@63.com ZHONGXIANG CHANG Tianjin
More informationKey words and phrases. Bifurcation, Difference Equations, Fixed Points, Predator - Prey System, Stability.
ISO 9001:008 Certified Volume, Issue, March 013 Dynamical Behavior in a Discrete Prey- Predator Interactions M.ReniSagaya Raj 1, A.George Maria Selvam, R.Janagaraj 3.and D.Pushparajan 4 1,,3 Sacred Heart
More informationAnalysis of a delayed predator-prey model with ratio-dependent functional response and quadratic harvesting. Peng Feng
Analysis of a delayed predator-prey model with ratio-dependent functional response and quadratic harvesting Peng Feng Journal of Applied Mathematics and Computing ISSN 1598-5865 J. Appl. Math. Comput.
More informationInteractions between predators and prey
Interactions between predators and prey What is a predator? Predator An organism that consumes other organisms and inevitably kills them. Predators attack and kill many different prey individuals over
More informationStability and bifurcation in a two species predator-prey model with quintic interactions
Chaotic Modeling and Simulation (CMSIM) 4: 631 635, 2013 Stability and bifurcation in a two species predator-prey model with quintic interactions I. Kusbeyzi Aybar 1 and I. acinliyan 2 1 Department of
More informationGlobal Stability Analysis on a Predator-Prey Model with Omnivores
Applied Mathematical Sciences, Vol. 9, 215, no. 36, 1771-1782 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.512 Global Stability Analysis on a Predator-Prey Model with Omnivores Puji Andayani
More informationA Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term
ETASR - Engineering, Technology & Applied Science Research Vol., o.,, 9-5 9 A Novel Three Dimension Autonomous Chaotic System with a Quadratic Exponential Nonlinear Term Fei Yu College of Information Science
More information1.Introduction: 2. The Model. Key words: Prey, Predator, Seasonality, Stability, Bifurcations, Chaos.
Dynamical behavior of a prey predator model with seasonally varying parameters Sunita Gakkhar, BrhamPal Singh, R K Naji Department of Mathematics I I T Roorkee,47667 INDIA Abstract : A dynamic model based
More informationNumerical bifurcation analysis of a tri-trophic food web with omnivory
Mathematical Biosciences 177&178 (2002) 201 228 www.elsevier.com/locate/mbs Numerical bifurcation analysis of a tri-trophic food web with omnivory B.W. Kooi *, L.D.J. Kuijper, M.P. Boer 1, S.A.L.M. Kooijman
More informationExact multiplicity of boundary blow-up solutions for a bistable problem
Computers and Mathematics with Applications 54 (2007) 1285 1292 www.elsevier.com/locate/camwa Exact multiplicity of boundary blow-up solutions for a bistable problem Junping Shi a,b,, Shin-Hwa Wang c a
More informationChaos and adaptive control in two prey, one predator system with nonlinear feedback
Chaos and adaptive control in two prey, one predator system with nonlinear feedback Awad El-Gohary, a, and A.S. Al-Ruzaiza a a Department of Statistics and O.R., College of Science, King Saud University,
More informationAsynchronous and Synchronous Dispersals in Spatially Discrete Population Models
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 7, No. 2, pp. 284 310 c 2008 Society for Industrial and Applied Mathematics Asynchronous and Synchronous Dispersals in Spatially Discrete Population Models Abdul-Aziz
More informationDIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS
WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Series Editor: Leon O. Chua Series A Vol. 66 DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS Jean-Marc Ginoux Université du Sud, France World Scientific
More informationInfluence of Noise on Stability of the Ecosystem
Commun. Theor. Phys. 60 (2013) 510 514 Vol. 60, No. 4, October 15, 2013 Influence of Noise on Stability of the Ecosystem LIU Tong-Bo ( ) and TANG Jun (» ) College of Science, China University of Mining
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationSome Recent Results on Diffusive Predator-prey Models in Spatially Heterogeneous Environment
Some Recent Results on Diffusive Predator-prey Models in Spatially Heterogeneous Environment Yihong Du and Junping Shi Abstract We present several recent results obtained in our attempts to understand
More informationdv dt Predator-Prey Models
Predator-Prey Models This is a diverse area that includes general models of consumption: Granivores eating seeds Parasitoids Parasite-host interactions Lotka-Voterra model prey and predator: V = victim
More informationA Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More informationMotivation and Goals. Modelling with ODEs. Continuous Processes. Ordinary Differential Equations. dy = dt
Motivation and Goals Modelling with ODEs 24.10.01 Motivation: Ordinary Differential Equations (ODEs) are very important in all branches of Science and Engineering ODEs form the basis for the simulation
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationGlobal dynamics of a predator-prey system with Holling type II functional response
4 Nonlinear Analysis: Modelling and Control, 011, Vol. 16, No., 4 53 Global dynamics of a predator-prey system with Holling type II functional response Xiaohong Tian, Rui Xu Institute of Applied Mathematics
More informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
More informationDynamical Systems in Neuroscience: Elementary Bifurcations
Dynamical Systems in Neuroscience: Elementary Bifurcations Foris Kuang May 2017 1 Contents 1 Introduction 3 2 Definitions 3 3 Hodgkin-Huxley Model 3 4 Morris-Lecar Model 4 5 Stability 5 5.1 Linear ODE..............................................
More informationFeedback-mediated oscillatory coexistence in the chemostat
Feedback-mediated oscillatory coexistence in the chemostat Patrick De Leenheer and Sergei S. Pilyugin Department of Mathematics, University of Florida deleenhe,pilyugin@math.ufl.edu 1 Introduction We study
More informationOn the stabilizing effect of specialist predators on founder-controlled communities
On the stabilizing effect of specialist predators on founder-controlled communities Sebastian J. Schreiber Department of Mathematics Western Washington University Bellingham, WA 98225 May 2, 2003 Appeared
More informationTURING AND HOPF PATTERNS FORMATION IN A PREDATOR-PREY MODEL WITH LESLIE-GOWER-TYPE FUNCTIONAL RESPONSE
Dynamics of Continuous, Discrete and Impulsive Systems Series B: Algorithms and Applications 16 2009) 479-488 Copyright c 2009 Watam Press http://www.watam.org TURING AND HOPF PATTERNS FORMATION IN A PREDATOR-PREY
More informationDYNAMIC STUDY OF A PREDATOR-PREY MODEL WITH WEAK ALLEE EFFECT AND HOLLING TYPE-III FUNCTIONAL RESPONSE
Dynamic Systems and Applications, 7, No 4 (018), 943-953 ISSN: 1056-176 DYNAMIC STUDY OF A PREDATOR-PREY MODEL WITH WEAK ALLEE EFFECT AND HOLLING TYPE-III FUNCTIONAL RESPONSE HUA LIU 1, YONG YE, YUMEI
More informationGLOBAL DYNAMICS OF A LOTKA VOLTERRA MODEL WITH TWO PREDATORS COMPETING FOR ONE PREY JAUME LLIBRE AND DONGMEI XIAO
This is a preprint of: Global dynamics of a Lotka-Volterra model with two predators competing for one prey, Jaume Llibre, Dongmei Xiao, SIAM J Appl Math, vol 742), 434 453, 214 DOI: [11137/1392397] GLOBAL
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationx 2 F 1 = 0 K 2 v 2 E 1 E 2 F 2 = 0 v 1 K 1 x 1
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 20, Number 4, Fall 1990 ON THE STABILITY OF ONE-PREDATOR TWO-PREY SYSTEMS M. FARKAS 1. Introduction. The MacArthur-Rosenzweig graphical criterion" of stability
More informationThe Hopf-van der Pol System: Failure of a Homotopy Method
DOI.7/s259--9-5 ORIGINAL RESEARCH The Hopf-van der Pol System: Failure of a Homotopy Method H. G. E. Meijer T. Kalmár-Nagy Foundation for Scientific Research and Technological Innovation 2 Abstract The
More informationTuring and Non-Turing patterns in diffusive plankton model
Computational Ecology and Software 05 5(): 6-7 Article Turing and Non-Turing patterns in diffusive plankton model N. K. Thakur R. Gupta R. K. Upadhyay National Institute of Technology Raipur India Indian
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationSPIRAL WAVE GENERATION IN A DIFFUSIVE PREDATOR-PREY MODEL WITH TWO TIME DELAYS
Bull. Korean Math. Soc. 52 (2015), No. 4, pp. 1113 1122 http://dx.doi.org/10.4134/bkms.2015.52.4.1113 SPIRAL WAVE GENERATION IN A DIFFUSIVE PREDATOR-PREY MODEL WITH TWO TIME DELAYS Wenzhen Gan and Peng
More informationDynamical Systems: Ecological Modeling
Dynamical Systems: Ecological Modeling G Söderbacka Abstract Ecological modeling is becoming increasingly more important for modern engineers. The mathematical language of dynamical systems has been applied
More informationReaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations
Reaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations Junping Shi College of William and Mary, USA Collaborators/Support Shanshan Chen (PhD, 2012; Harbin
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More informationLab 5: Nonlinear Systems
Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting
More informationA Comparison of Two Predator-Prey Models with Holling s Type I Functional Response
A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response ** Joint work with Mark Kot at the University of Washington ** Mathematical Biosciences 1 (8) 161-179 Presented by Gunog
More informationCANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION
CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION WARREN WECKESSER Department of Mathematics Colgate University Hamilton, NY 3346 E-mail: wweckesser@mail.colgate.edu Cartwright and Littlewood discovered
More informationBIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences
BIOS 6150: Ecology Dr. Stephen Malcolm, Department of Biological Sciences Week 7: Dynamics of Predation. Lecture summary: Categories of predation. Linked prey-predator cycles. Lotka-Volterra model. Density-dependence.
More informationStability analysis of a prey-predator model with a reserved area
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research, 4, 5(3):93-3 ISSN: 976-86 CODEN (USA): AASRFC Stability analysis of a prey-predator model with a reserved area Neelima
More informationDynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses
Journal of Physics: Conference Series PAPER OPEN ACCESS Dynamics Analysis of Anti-predator Model on Intermediate Predator With Ratio Dependent Functional Responses To cite this article: D Savitri 2018
More informationVITA EDUCATION. Arizona State University, , Ph.D. Huazhong Univ. of Sci. and Tech., , M.S. Lanzhou University, , B.S.
VITA Bingtuan Li Department of Mathematics University of Louisville Louisville, KY 40292 Office phone: (502) 852-6149 Fax: (502) 852-7132 E-mail: bing.li@louisville.edu EDUCATION Arizona State University,
More informationDynamical Systems in Biology
Dynamical Systems in Biology Hal Smith A R I Z O N A S T A T E U N I V E R S I T Y H.L. Smith (ASU) Dynamical Systems in Biology ASU, July 5, 2012 1 / 31 Outline 1 What s special about dynamical systems
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationNa#onal Center for Theore#cal Sciences Mathema#cs Division, Taiwan
Na#onal Center for Theore#cal Sciences Mathema#cs Division, Taiwan - TAMKANG JOURNAL OF MATHEMATICS Volume 47, Number 1, 127-141, March 2016 doi:10.5556/j.tkjm.47.2016.1984 This paper is available online
More informationBifurcation Analysis, Chaos and Control in the Burgers Mapping
ISSN 1749-3889 print, 1749-3897 online International Journal of Nonlinear Science Vol.4007 No.3,pp.171-185 Bifurcation Analysis, Chaos and Control in the Burgers Mapping E. M. ELabbasy, H. N. Agiza, H.
More information3.5 Competition Models: Principle of Competitive Exclusion
94 3. Models for Interacting Populations different dimensional parameter changes. For example, doubling the carrying capacity K is exactly equivalent to halving the predator response parameter D. The dimensionless
More informationarxiv: v1 [math.ds] 30 Jan 2012
Homoclinic Orbits of the FitzHugh-Nagumo Equation: Bifurcations in the Full System John Guckenheimer Christian Kuehn arxiv:1201.6352v1 [math.ds] 30 Jan 2012 Abstract This paper investigates travelling
More informationA NUMERICAL STUDY ON PREDATOR PREY MODEL
International Conference Mathematical and Computational Biology 2011 International Journal of Modern Physics: Conference Series Vol. 9 (2012) 347 353 World Scientific Publishing Company DOI: 10.1142/S2010194512005417
More informationDynamical Systems and Chaos Part II: Biology Applications. Lecture 6: Population dynamics. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part II: Biology Applications Lecture 6: Population dynamics Ilya Potapov Mathematics Department, TUT Room TD325 Living things are dynamical systems Dynamical systems theory
More informationCDS 101/110a: Lecture 2.1 Dynamic Behavior
CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium
More informationSpatiotemporal pattern formation in a prey-predator model under environmental driving forces
Home Search Collections Journals About Contact us My IOPscience Spatiotemporal pattern formation in a prey-predator model under environmental driving forces This content has been downloaded from IOPscience.
More informationPhenomenon: Canadian lynx and snowshoe hares
Outline Outline of Topics Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Phenomenon: Canadian lynx and snowshoe hares All began with
More informationEcology 302: Lecture VII. Species Interactions.
Ecology 302: Lecture VII. Species Interactions. (Gotelli, Chapters 6; Ricklefs, Chapter 14-15) MacArthur s warblers. Variation in feeding behavior allows morphologically similar species of the genus Dendroica
More informationLOTKA-VOLTERRA SYSTEMS WITH DELAY
870 1994 133-140 133 LOTKA-VOLTERRA SYSTEMS WITH DELAY Zhengyi LU and Yasuhiro TAKEUCHI Department of Applied Mathematics, Faculty of Engineering, Shizuoka University, Hamamatsu 432, JAPAN ABSTRACT Sufftcient
More informationPREDATOR-PREY SYSTEM WITH EXTINCTION DYNAMICS YONG-JUNG KIM
PREDATOR-PREY SYSTEM WITH EXTINCTION DYNAMICS YONG-JUNG KIM Abstract. Most of population models, if not all of them, do not have the extinction dynamics. The Lotka-Volterra ordinary differential equations
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationGerardo Zavala. Math 388. Predator-Prey Models
Gerardo Zavala Math 388 Predator-Prey Models Spring 2013 1 History In the 1920s A. J. Lotka developed a mathematical model for the interaction between two species. The mathematician Vito Volterra worked
More informationComputational Neuroscience. Session 4-2
Computational Neuroscience. Session 4-2 Dr. Marco A Roque Sol 06/21/2018 Two-Dimensional Two-Dimensional System In this section we will introduce methods of phase plane analysis of two-dimensional systems.
More informationarxiv: v1 [nlin.ao] 19 May 2017
Feature-rich bifurcations in a simple electronic circuit Debdipta Goswami 1, and Subhankar Ray 2, 1 Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA
More informationInternational Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: Vol.8, No.7, pp , 2015
International Journal of PharmTech Research CODEN (USA): IJPRIF, ISSN: 0974-4304 Vol.8, No.7, pp 99-, 05 Lotka-Volterra Two-Species Mutualistic Biology Models and Their Ecological Monitoring Sundarapandian
More informationPERMANENCE IN LOGISTIC AND LOTKA-VOLTERRA SYSTEMS WITH DISPERSAL AND TIME DELAYS
Electronic Journal of Differential Equations, Vol. 2005(2005), No. 60, pp. 1 11. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) PERMANENCE
More informationStability and Hopf bifurcation of a ratio-dependent predator prey model with time delay and stage structure
Electronic Journal of Qualitative Theory of Differential Equations 6, No 99, 3; doi: 43/ejqtde699 http://wwwmathu-szegedhu/ejqtde/ Stability and Hopf bifurcation of a ratio-dependent predator prey model
More informationThe functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z).
Review of Multi-variable calculus: The functions in all models depend on two variables: time t and spatial variable x, (x, y) or (x, y, z). The spatial variable represents the environment where the species
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationM.5 Modeling the Effect of Functional Responses
M.5 Modeling the Effect of Functional Responses The functional response is referred to the predation rate as a function of the number of prey per predator. It is recognized that as the number of prey increases,
More informationWe have two possible solutions (intersections of null-clines. dt = bv + muv = g(u, v). du = au nuv = f (u, v),
Let us apply the approach presented above to the analysis of population dynamics models. 9. Lotka-Volterra predator-prey model: phase plane analysis. Earlier we introduced the system of equations for prey
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationEnrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination
Enrichment in a Producer-Consumer Model with varying rates of Stoichiometric Elimination Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth
More informationDIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationPrey-Predator Model with a Nonlocal Bistable Dynamics of Prey. Received: 5 February 2018; Accepted: 5 March 2018; Published: 8 March 2018
mathematics Article Prey-Predator Model with a Nonlocal Bistable Dynamics of Prey Malay Banerjee 1, *,, Nayana Mukherjee 1, and Vitaly Volpert 2, 1 Department of Mathematics and Statistics, IIT Kanpur,
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationNumerical Solution of a Fractional-Order Predator-Prey Model with Prey Refuge and Additional Food for Predator
66 Numerical Solution of a Fractional-Order Predator-Prey Model with Prey Refuge Additional Food for Predator Rio Satriyantara, Agus Suryanto *, Noor Hidayat Department of Mathematics, Faculty of Mathematics
More informationNON-CONSTANT POSITIVE STEADY STATES FOR A STRONGLY COUPLED NONLINEAR REACTION-DIFFUSION SYSTEM ARISING IN POPULATION DYNAMICS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 45, Number 4, 2015 NON-CONSTANT POSITIVE STEADY STATES FOR A STRONGLY COUPLED NONLINEAR REACTION-DIFFUSION SYSTEM ARISING IN POPULATION DYNAMICS ZIJUAN WEN
More informationCOMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 2, Number, Spring 22 COMPETITION OF FAST AND SLOW MOVERS FOR RENEWABLE AND DIFFUSIVE RESOURCE SILOGINI THANARAJAH AND HAO WANG ABSTRACT. In many studies of
More informationNumerical Bifurcation Analysis of Dynamical Systems. Bob W. Kooi Faculty of Earth and Life Sciences,
MRI Master Class 2009/2010: Numerical Bifurcation Analysis of Dynamical Systems WEEK 15, Apr 13 & 15: Bob W. Kooi Faculty of Earth and Life Sciences, Department of Theoretical Biology, Vrije University,
More information