Continuous time population models
|
|
- Samuel McLaughlin
- 5 years ago
- Views:
Transcription
1 Continuous time population models Jaap van der Meer Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations by means of a coupled set of non-linear differential equations to the population sizes. Various classical models are discussed, including a DEB model for a substrate eating V1-morph in a chemostat.
2 April 11, /39 Content One population Linear model: exponential growth Non-linear model: logistic growth Two populations Uncoupled linear model Coupled linear model Non-linear model: logistic predator-prey model DEB V1-morphs Chemostat models
3 April 11, /39 Exponential growth Rate of change is a linear function of population size itself, dx/dt = rx. Equilibria can be obtained by setting dx/dt = 0, which gives rise to one equilibrium x = 0. A general solution of this linear differential equation can be obtained by separation of variables and subsequent integration, 1 x dx = rdt, and is given by lnx = rt + c or x = x 0 e rt. If time t goes to infinity, population size x will either be infinitely large if r > 0 or will approach zero if r < 0, whatever the initial value x 0 is.
4 April 11, /39 Exponential growth Population size x Time t Figure 1: Exponential growth either means unlimited growth in case the instantaneous growth rate r is positive, or extinction when the rate is negative.
5 April 11, /39 Exponential growth Rate of change dx/dt Population size x Figure 2: The equilibrium of the exponential growth equation, which is given by x = 0 is either unstable in case the slope r of the linear relationship between population rate of change dx/dt and population size x is positive, or stable when the slope r is negative. The slope might be called the eigenvalue of the equation.
6 April 11, /39 Logistic growth The logistic equation, given by dx/dt = rx(1 x/k), has no unbounded growth. The two equilibria are x = 0 and x = k. Generally, a tangent in x = a is given by y = f(a) + f (a)(x a), where f (a) is the first derivative of the function f(x) in x = a. Here, the tangent in x = k is given by dx/dt = r(x k), and the linear approximation in the equilibrium x = k can be re-written as dy/dt = ry, where y = x k. Nothing else than the exponential growth equation.
7 April 11, /39 Logistic growth Rate of change dx/dt Population size x Figure 3: The logistic growth equation has two equilibria. The red and blue lines show linear approximations of the growth equation in the neighbourhood of the equilibria. One equilibrium is locally unstable, the other is locally stable.
8 April 11, /39 Logistic growth Population size x Time t Figure 4: Four particular solutions for the logistic growth equation are given by the black lines. The red and blue lines show various other specific solutions in the neighbourhood of the equilibria. They resemble exponential growth.
9 April 11, /39 Allee effect Rate of change dx/dt Population size x Figure 5: A cubic relationship between dx/dt and x may give rise to multiple equilibria, which are either stable or unstable. As a result of a slight change in the parameter values the equilibria can easily disappear, as depicted by the green lines.
10 April 11, /39 Two populations If two populations grow exponentially, we may write the differential equations in matrix form [ dx0 dt dx 1 dt = [ λ0 0 0 λ 1 [ x0 x 1 (in short notation x = Ax) and plot the trajectories in a phase plane diagram, in which x 1 is plotted versus x 0. If the two growth coefficients (or eigenvalues) are negative the equilibrium is stable. The two populations are decoupled and the analysis is exactly the same as discussed for a single population.
11 April 11, /39 Figure 6: Node, where λ 2 < λ 1 < 0.
12 April 11, /39 Coupled populations But what do we do when the growth of the two populations is coupled, such that the differential equations in matrix form look like [ dx0 dt dx 1 dt = [ a b c d [ x0 x 1 where the elements b and c are not necessarily zero? We need some matrix algebra.
13 April 11, /39 A new basis (3, 2) (1,1) Figure 7: A new basis.
14 April 11, /39 Consider now the multiplication of a vector (expressed with respect to the standard basis) with a matrix. For example [ [ [ Ax = = This matrix A does with respect to the standard basis β the same thing as the matrix B does with respect to the basis η: [ [ [ By = = Note that [ [ = [ 7 13
15 April 11, /39 Eigenvectors as the new basis Figure 8: A new basis is shown with blue vectors. The red vectors are the standard basis. The purple arrow refers to the example. Points on a blue vector (or its elongation), such as the purple point in the upper right, will stay on that vector.
16 April 11, /39 Eigenvectors and eigenvalues The fundamental equation needed is Ax = λx, where x, which is not the zero vector x = 0, is an eigenvector and λ an eigenvalue of A. Hence, if one multiplies the matrix A with its eigenvector, the result is simply a scalar multiplication of the eigenvector. A two by two matrix has at most two eigenvectors. The equation Ax = λx is re-written as (A λi)x = 0, where is called the identity matrix. I = [
17 April 11, /39 The determinant The equation Ax = [ a b c d [ x1 x 2 is true if ax 1 +bx 2 = 0 and cx 1 +dx 2 = 0, hence if x 2 = a b x 1 and if x 2 = c d x 1. A solution, other than the zero vector, only exists if a/b = c/d, which results in ad bc = 0. The term ad bc is called the determinant of the matrix A. Hence, the equation (A λi)x = 0 has only a solution, other than x = 0, if the determinant of the matrix A λi is zero, hence if (a λ)(d λ) bc = 0. = [ 0 0
18 April 11, /39 The characteristic equation This leads to the so-called characteristic equation: The eigenvalues are thus λ 2 (a + d)λ + (ad bc) = λ 2 T λ + D = 0 λ 1,2 = (a + d) ± (a + d) 2 4(ad bc) 2 = T ± T 2 4D 2 where T = a + d is the trace and D = ad bc the determinant of the matrix A. For the case that the part under the square root sign, which is called the discriminant, is positive, it can easily be seen that both eigenvalues are negative if and only if T < 0 and D > 0.
19 April 11, /39 Returning to the example A = [ gives the characteristic equation λ 2 + 3λ 4 = (λ + 4)(λ 1) = 0, resulting in the eigenvalues λ 1 = 4 and λ 1 = 1. The trace T = 3, the determinant D = 4 and the discriminant equals 25. The eigenvalue λ = 4 reveals (A + 4I)e 1 = [ e 1 = 0 Hence e 1 = [2 3 t is an accompanying eigenvector. Similarly, for the eigenvalue λ = 1.
20 April 11, /39 Coupled populations Hence, when the growth of the two populations is coupled, such that [ dx0 dt dx 1 dt = [ a b c d [ x0 x 1 we can now easily re-write it in the form of [ dy0 dt dy 1 dt = [ λ0 0 0 λ 1 [ y0 y 1 for which we know the solution (or don t we?).
21 April 11, /39 Figure 9: Saddle, where λ 1 < 0 and λ 2 > 0. The orange lines are the eigenvectors. The purple lines the isoclines, for which dx i /dt = 0
22 April 11, /39 Complex numbers But what if the discriminant not positive? λ 1,2 = (a + d) ± (a + d) 2 4(ad bc) 2 = T ± T 2 4D 2 We enter the realm of imaginary numbers i = 1, but thanks to Euler can we re-write them in terms of sines and cosines. Noting that i 2 = 1, i 3 = i, i 4 = 1, i 5 = i, i 6 = 1, etc. gives you some intuitive understanding of this link. The practical implication is that the stability properties are determined by the real part, and that the equilibrium is either approached in spirals (if the real part is negative) or that the system moves away from the equilibrium in ever extending spirals (if the real part is positive).
23 April 11, /39 A logistic predator-prey model A predator-prey model, in which the prey population follows logistic growth in the absence of predation, and where the predation term is based on the laws of mass-action, can be written (in a dimensionless form) as dx 0 dτ = x 0(1 x 0 ) qx 0 x 1 dx 1 dτ = pqx 0x 1 mx 1 where x 0 is the prey and x 1 the predator population size. ( ) The biologically interesting equilibrium is given by x 0 = m pq and x 1 = 1 q 1 m pq. Again we can use a linear approximation around the equilibrium, using the proposition that a non-linear function f(x 0, x 1 ) around the point x 0 = a and
24 April 11, /39 x 1 = b can be approximated by the tangent plane, which is given by y = f(a, b) + f x 0 (x 0 a) + f x 1 (x 1 b) where the partial derivatives are taken at x 0 = a and x 1 = b. For our purpose this leads to the matrix equation: [ d(x0 x 0 ) dτ d(x 1 x 1 ) dτ = f (x,x ) [ f0 (x 0,x 1 ) x 0 f 1 (x 0,x 1 ) x 1 f 1 (x 0,x 1 ) x 0 x 1 x=x [ x0 x 0 x 1 x 1 where f i (x 0, x 1 ) is the differential equation for x i.
25 April 11, /39 The Jacobian equals [ x 0 qx 0 p(1 x 0) 0 = [ m pq m p p(1 m pq ) 0 which leads to a trace T = x 0 = m pq that is always negative (given that the parameters m, p and q are positive) and a determinant D = mqx 1 = m(1 m pq ) that is positive as long as x 1 > 0. Hence the equilibrium is stable. The discriminant T 2 4D = x 2 0 4m( x 0) can be either positive or negative, implying that the equilibrium is or is not approached in spirals. For example, if q = 10, p = 0.5 and m = 0.1, then T = 0.02, D = and the discriminant equals , implying complex eigenvalues. The stable equilibrium x 0 = 0.02 and x 1 = is approached in spirals.
26 April 11, /39 x e 05 1e 04 1e 03 1e 02 1e 01 x0 Figure 10: A specific solution of the logistic predator-prey model. lines give the isoclines, which are defined by dx i /dτ = 0. The purple
27 April 11, /39 V1-morphs Reserve density follows first order dynamics d[e dt = ṗa V c[e = [ṗ Amf k E [E (1) where k E has been given the name specific energy conductance. It has the physical dimension per time. The mobilization rate equals ṗ C = k E [EV [E dv dt (2)
28 April 11, /39 For organisms that simply divide into two daughter cells, and which are classified as juveniles in DEB terminology, κ can be set equal to one. Hence the allocation is given by ṗ C = [E G dv dt + [ṗ MV (3) Substituting equation 2 in equation 3 gives the growth equation dv dt = k E [E [ṗ M V (4) [E + [E G
29 April 11, /39 Under constant food conditions, when the reserve density is in equilibrium and proportional to the scaled functional response, the growth equation simplifies to dv dt = [ṗ Amf [ṗ M 1 [ṗ Am f k E + [E V = ṙv (5) G where ṙ is the specific growth rate. Hence, the growth rate of the structural volume is proportional to the structural volume itself. V1-morphs show exponential growth at constant food density.
30 April 11, /39 A mass-mass framework One might choose to write the reserve dynamics and the growth equation in a mass-mass framework. Equation 1 becomes dm E dt = d[e dt 1 µ E [M V = j EAmf k E m E (6) and Equation 4 becomes dm V dt = dv dt [M V = k E m E j EM m E + y EV M V (7)
31 April 11, /39 Table 1: Primary parameters of the DEB model for V1-morphs, specific for a mass-mass framework. The last column indicates the relationship with the parameters from the energy-length framework that has been replaced. µ E is the potential energy of the reserves expressed in energy per C-mole, and [M V is the specific density of the structural body expressed in C-moles per volume. Symbol Dimension Interpretation Relationship j EAm ## 1 t 1 Mass-specific [ṗ Am = µ E [M V j EAm maximum assimilation rate y EX ## 1 Yield of reserve µ AX = µ E y EX on food j EM ## 1 t 1 Mass-specific [ṗ M = µ E [M V j EM maintenance rate y EV ## 1 Mass-specific costs of growth [E G = µ E [M V y EV
32 April 11, /39 Yield The ratio of yield on substrate use is given by the rate of increase of structural mass and reserves divided by the ingestion rate: Y = dm V /dt + d(m E M V )/dt (8) J X Assuming constant substrate availability, and hence that the reserve density is in equilibrium, it follows that Y = (1 + m E )ṙm V J X (9)
33 April 11, /39 The yield on substrate can be expressed as a function of the food availability only Y = y V X y EV f + g f + g f l d f (10) where y V X is known as the true yield on substrate in microbiology. The compound parameter g is the energy investment ratio. It is defined as the ratio of the energetic costs of structure [E G relative to the maximum available energy for growth and maintenance κ[e m. For dividing V1-morphs, with κ = 1, g can be written as k E [E G /[ṗ Am within an energy-length framework. The compound parameter l d is defined as the ratio between the volume-specific maintenance rate [ṗ M and the volume-specific maximum assimilation rate [ṗ Am.
34 April 11, /39 The chemostat Total food intake rate of the V1-morphs equals J X = [ J Xm fv, where [ J Xm is the maximum volume-specific intake rate, f the scaled functional response, and V the overall volume of the population of V1-morphs. The rates of change of the food density and of the biomass density of V1- morphs X 1 = M V /V C look like: dx 0 dt = ḣ (X J r X 0 ) X dx 1 V C dt = Y J X V ḣx 1 C where ḣ is the throughput rate, and V C the constant volume of the chemostat.
35 April 11, /39 A dimensionless form The model can be expressed in a dimensionless form resulting in dx 0 dτ = x r x 0 αfx 1 dx 1 dτ = Y αfx 1 x 1
36 April 11, /39 Alternative models Y = y V X g f l d f + g f (11) DEB V1-morphs Droop model has no maintenance: l d = 0 Monod model has no reserves: g Lotka-Volterra kind of model has no handling time: f x 0
37 April 11, /39 The Lotka-Volterra kind of model Equilibria x 0 = (αy ) 1 x 1 = 1 α (x rαy 1) The Jacobian of this equilibrium equals: f (x,x ) [ f0 (x 0,x 1 ) x 0 f 1 (x 0,x 1 ) x 1 f 1 (x 0,x 1 ) x 0 x 1 x=x = [ αx 1 1 αx 0 Y αx 1 Y αx 0 1 = [ x r αy Y (x r αy 1) 0 1 Y where f i (x 0, x 1 ) is the differential equation for x i. Eigenvalues equal λ 1 = (x r αy 1) and λ 2 = 1. Both eigenvalues are real and negative, pointing to a so-called node.
38 April 11, /39 Figure 11: Simulated trajectories of the Lotka-Volterra kind of model. Isoclines in purple.
39 April 11, /39 Table 2: A numerical example of the chemostat models. The first part gives the parameter values, the second the equilibrium, the third the elements of the Jacobian, and the fourth the eigenvalues. Lotka Monod Droop DEB x r y V X α g l d x x a b c d λ λ
40 April 11, /39 References [1 Gurney, M. and Nisbet, R. (1998) Ecological dynamics. Publisher, Place [2 Kooijman, S.A.L.M. (2010) Dynamic Energy Budget theory for metabolic organization. Cambridge University Press, Cambridge [3 Strang, G. (1980) Linear algebra and its applications. Second edition. Academic Press, New York
Physics: 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 informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationNon-Dimensionalization of Differential Equations. September 9, 2014
Non-Dimensionalization of Differential Equations September 9, 2014 1 Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side
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 informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
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 informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More information2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
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 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 informationFundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More informationWorkshop on Theoretical Ecology and Global Change March 2009
2022-3 Workshop on Theoretical Ecology and Global Change 2-18 March 2009 Stability Analysis of Food Webs: An Introduction to Local Stability of Dynamical Systems S. Allesina National Center for Ecological
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationMath 232, Final Test, 20 March 2007
Math 232, Final Test, 20 March 2007 Name: Instructions. Do any five of the first six questions, and any five of the last six questions. Please do your best, and show all appropriate details in your solutions.
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationModeling Prey and Predator Populations
Modeling Prey and Predator Populations Alison Pool and Lydia Silva December 15, 2006 Abstract In this document, we will explore the modeling of two populations based on their relation to one another. Specifically
More informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
More informationMath Lecture 46
Math 2280 - Lecture 46 Dylan Zwick Fall 2013 Today we re going to use the tools we ve developed in the last two lectures to analyze some systems of nonlinear differential equations that arise in simple
More information6.3. Nonlinear Systems of Equations
G. NAGY ODE November,.. Nonlinear Systems of Equations Section Objective(s): Part One: Two-Dimensional Nonlinear Systems. ritical Points and Linearization. The Hartman-Grobman Theorem. Part Two: ompeting
More informationMA 138: Calculus II for the Life Sciences
MA 138: Calculus II for the Life Sciences David Murrugarra Department of Mathematics, University of Kentucky. Spring 2016 David Murrugarra (University of Kentucky) MA 138: Section 11.4.2 Spring 2016 1
More informationCalculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1
Calculus I Homework: The Derivatives of Polynomials and Exponential Functions Page 1 Questions Example Differentiate the function y = ae v + b v + c v 2. Example Differentiate the function y = A + B x
More informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
More informationMath 312 Lecture Notes Linearization
Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationMatrices, Linearization, and the Jacobi matrix. y f x g y g J = dy/dt = g(x, y) Theoretical Biology, Utrecht University
Matrices, Linearization, and the Jacobi matrix { dx/dt f(x, y dy/dt g(x, y ( x f J y f x g y g λ 1, tr ± tr 4 det Theoretical Biology, Utrecht University i c Utrecht University, 018 Ebook publically available
More informationMA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, Exam Scores. Question Score Total
MA 138 Calculus 2 for the Life Sciences Spring 2016 Final Exam May 4, 2016 Exam Scores Question Score Total 1 10 Name: Section: Last 4 digits of student ID #: No books or notes may be used. Turn off all
More information+ i. cos(t) + 2 sin(t) + c 2.
MATH HOMEWORK #7 PART A SOLUTIONS Problem 7.6.. Consider the system x = 5 x. a Express the general solution of the given system of equations in terms of realvalued functions. b Draw a direction field,
More informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More informationSection 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d
Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d July 6, 2009 Today s Session Today s Session A Summary of This Session: Today s Session A Summary of This
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
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 informationInteractions. Yuan Gao. Spring Applied Mathematics University of Washington
Interactions Yuan Gao Applied Mathematics University of Washington yuangao@uw.edu Spring 2015 1 / 27 Nonlinear System Consider the following coupled ODEs: dx = f (x, y). dt dy = g(x, y). dt In general,
More informationPhase Plane Analysis
Phase Plane Analysis Phase plane analysis is one of the most important techniques for studying the behavior of nonlinear systems, since there is usually no analytical solution for a nonlinear system. Background
More informationCalculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y
Calculus for the Life Sciences II Assignment 6 solutions Find the tangent plane to the graph of the function at the point (0, π f(x, y = 3π 3 cos 2x + 2 sin 3y Solution: The tangent plane of f at a point
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 informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations Scott A. McKinley October 22, 2013 In these notes, which replace the material in your textbook, we will learn a modern view of analyzing systems of differential
More informationMath 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Two-dimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationMatrices, Linearization, and the Jacobi matrix
Matrices, Linearization, and the Jacobi matrix { dx/dt f(x, y dy/dt g(x, y ( x f J y f x g y g λ 1, tr ± tr 4 det Alexander V. Panfilov, Kirsten H.W.J. ten Tusscher & Rob J. de Boer Theoretical Biology,
More informationLecture 3. Dynamical Systems in Continuous Time
Lecture 3. Dynamical Systems in Continuous Time University of British Columbia, Vancouver Yue-Xian Li November 2, 2017 1 3.1 Exponential growth and decay A Population With Generation Overlap Consider a
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
More informationMAT 22B - Lecture Notes
MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationE209A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions
E09A: Analysis and Control of Nonlinear Systems Problem Set 3 Solutions Michael Vitus Stanford University Winter 007 Problem : Planar phase portraits. Part a Figure : Problem a This phase portrait is correct.
More informationNon-Linear Models. Non-Linear Models Cont d
Focus on more sophistiated interaction models between systems. These lead to non-linear, rather than linear, DEs; often not soluble exactly in analytical form so use Phase-Plane Analysis. This is a method
More information1 (t + 4)(t 1) dt. Solution: The denominator of the integrand is already factored with the factors being distinct, so 1 (t + 4)(t 1) = A
Calculus Topic: Integration of Rational Functions Section 8. # 0: Evaluate the integral (t + )(t ) Solution: The denominator of the integrand is already factored with the factors being distinct, so (t
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 informationMath 1280 Notes 4 Last section revised, 1/31, 9:30 pm.
1 competing species Math 1280 Notes 4 Last section revised, 1/31, 9:30 pm. This section and the next deal with the subject of population biology. You will already have seen examples of this. Most calculus
More informationOutline. Learning Objectives. References. Lecture 2: Second-order Systems
Outline Lecture 2: Second-order Systems! Techniques based on linear systems analysis! Phase-plane analysis! Example: Neanderthal / Early man competition! Hartman-Grobman theorem -- validity of linearizations!
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
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 informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)
ENGI 940 4.06 - Linear Approximation () Page 4. 4.06 Linear Approximation to a System of Non-Linear ODEs () From sections 4.0 and 4.0, the non-linear system dx dy = x = P( x, y), = y = Q( x, y) () with
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 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 informationHOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION)
HOMEWORK PROBLEMS FROM STRANG S LINEAR ALGEBRA AND ITS APPLICATIONS (4TH EDITION) PROFESSOR STEVEN MILLER: BROWN UNIVERSITY: SPRING 2007 1. CHAPTER 1: MATRICES AND GAUSSIAN ELIMINATION Page 9, # 3: Describe
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
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 informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More informationModeling Microbial Populations in the Chemostat
Modeling Microbial Populations in the Chemostat 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) Modeling Microbial Populations in the Chemostat MBI, June 3, 204 / 34 Outline Why
More information4 Second-Order Systems
4 Second-Order Systems Second-order autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
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 informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationA Discrete Model of Three Species Prey- Predator System
ISSN(Online): 39-8753 ISSN (Print): 347-670 (An ISO 397: 007 Certified Organization) Vol. 4, Issue, January 05 A Discrete Model of Three Species Prey- Predator System A.George Maria Selvam, R.Janagaraj
More informationDEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS. MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 2018
DEPARTMENT OF MATHEMATICS AND STATISTICS UNIVERSITY OF MASSACHUSETTS MATH 233 SOME SOLUTIONS TO EXAM 2 Fall 208 Version A refers to the regular exam and Version B to the make-up. Version A. A particle
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 information1 Functions and Graphs
1 Functions and Graphs 1.1 Functions Cartesian Coordinate System A Cartesian or rectangular coordinate system is formed by the intersection of a horizontal real number line, usually called the x axis,
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
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 informationApril 9, 2009 Name The problems count as marked. The total number of points available is 160. Throughout this test, show your work.
April 9, 009 Name The problems count as marked The total number of points available is 160 Throughout this test, show your work 1 (15 points) Consider the cubic curve f(x) = x 3 + 3x 36x + 17 (a) Build
More informationMath 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv
Math 1270 Honors ODE I Fall, 2008 Class notes # 1 We have learned how to study nonlinear systems x 0 = F (x; y) y 0 = G (x; y) (1) by linearizing around equilibrium points. If (x 0 ; y 0 ) is an equilibrium
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More informationIntegration - Past Edexcel Exam Questions
Integration - Past Edexcel Exam Questions 1. (a) Given that y = 5x 2 + 7x + 3, find i. - ii. - (b) ( 1 + 3 ) x 1 x dx. [4] 2. Question 2b - January 2005 2. The gradient of the curve C is given by The point
More information24. Partial Differentiation
24. Partial Differentiation The derivative of a single variable function, d f(x), always assumes that the independent variable dx is increasing in the usual manner. Visually, the derivative s value at
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 informationSome Notes on Linear Algebra
Some Notes on Linear Algebra prepared for a first course in differential equations Thomas L Scofield Department of Mathematics and Statistics Calvin College 1998 1 The purpose of these notes is to present
More informationM469, Fall 2010, Practice Problems for the Final
M469 Fall 00 Practice Problems for the Final The final exam for M469 will be Friday December 0 3:00-5:00 pm in the usual classroom Blocker 60 The final will cover the following topics from nonlinear systems
More informationSolutions of Spring 2008 Final Exam
Solutions of Spring 008 Final Exam 1. (a) The isocline for slope 0 is the pair of straight lines y = ±x. The direction field along these lines is flat. The isocline for slope is the hyperbola on the left
More informationChapter 7. Nonlinear Systems. 7.1 Introduction
Nonlinear Systems Chapter 7 The scientist does not study nature because it is useful; he studies it because he delights in it, and he delights in it because it is beautiful. - Jules Henri Poincaré (1854-1912)
More informationCopyright (c) 2006 Warren Weckesser
2.2. PLANAR LINEAR SYSTEMS 3 2.2. Planar Linear Systems We consider the linear system of two first order differential equations or equivalently, = ax + by (2.7) dy = cx + dy [ d x x = A x, where x =, and
More information1 Simple one-dimensional dynamical systems birth/death and migration processes, logistic
NOTES ON A Short Course and Introduction to Dynamical Systems in Biomathematics Urszula Foryś Institute of Applied Math. & Mech. Dept. of Math., Inf. & Mech. Warsaw University 1 Simple one-dimensional
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationEigenvalues and Eigenvectors, More Direction Fields and Systems of ODEs
Eigenvalues and Eigenvectors, More Direction Fields and Systems of ODEs First let us speak a bit about eigenvalues. Defn. An eigenvalue λ of an nxn matrix A means a scalar (perhaps a complex number) such
More informationNonlinear differential equations - phase plane analysis
Nonlinear differential equations - phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, y-plane which satisfy this
More informationFinal Exam Review. Review of Systems of ODE. Differential Equations Lia Vas. 1. Find all the equilibrium points of the following systems.
Differential Equations Lia Vas Review of Systems of ODE Final Exam Review 1. Find all the equilibrium points of the following systems. (a) dx = x x xy (b) dx = x x xy = 0.5y y 0.5xy = 0.5y 0.5y 0.5xy.
More informationOrdinary Differential Equations (ODEs)
c01.tex 8/10/2010 22: 55 Page 1 PART A Ordinary Differential Equations (ODEs) Chap. 1 First-Order ODEs Sec. 1.1 Basic Concepts. Modeling To get a good start into this chapter and this section, quickly
More information7.1. Calculus of inverse functions. Text Section 7.1 Exercise:
Contents 7. Inverse functions 1 7.1. Calculus of inverse functions 2 7.2. Derivatives of exponential function 4 7.3. Logarithmic function 6 7.4. Derivatives of logarithmic functions 7 7.5. Exponential
More informationThe Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor
, pp. 35-46 http://dx.doi.org/10.14257/ijbsbt.2017.9.3.04 The Dynamic Behaviour of the Competing Species with Linear and Holling Type II Functional Responses by the Second Competitor Alemu Geleta Wedajo
More informationModeling Prey-Predator Populations
Modeling Prey-Predator Populations Alison Pool and Lydia Silva December 13, 2006 Alison Pool and Lydia Silva () Modeling Prey-Predator Populations December 13, 2006 1 / 25 1 Introduction 1 Our Populations
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 informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More informationPhase portraits in two dimensions
Phase portraits in two dimensions 8.3, Spring, 999 It [ is convenient to represent the solutions to an autonomous system x = f( x) (where x x = ) by means of a phase portrait. The x, y plane is called
More information6 EIGENVALUES AND EIGENVECTORS
6 EIGENVALUES AND EIGENVECTORS INTRODUCTION TO EIGENVALUES 61 Linear equations Ax = b come from steady state problems Eigenvalues have their greatest importance in dynamic problems The solution of du/dt
More information