9 Frequency-dependent selection and game theory
|
|
- Camron Marshall
- 6 years ago
- Views:
Transcription
1 9 Frequency-dependent selection and game theory 9.1 Frequency-dependent selection in a haploid population In many biological situations an individual s fitness is not constant but depends on the genetic structure of its own population or populations of other species the first population interacts with. In such situations, we say that fitnesses are frequency-dependent or that the population is subject to frequency-dependent selection. For example, in the case of warning coloration (e.g. in butterflies) individuals that have a rare coloration pattern may suffer from predation to a larger degree than common individuals. Another example is mating behavior of Drosophila where males with rare phenotypes often have mating advantage. Yet another example would be a population with within-population competition for two renewable resources where rare individuals specializing to one resource would have advantage. Let us consider a single haploid locus with two alleles: A 1 and A 2. Let p and q be the frequencies of alleles A 1 and A 2, respectively (p + q = 1), and w 1 and w 2 be the corresponding fitnesses. The dynamics of the frequency of allele A 1 are described by a difference equation p = pq(w 1 w 2 ). (81) w Note that a polymorphic equilibrium w 1 = w 2. Let us assume that fitnesses depend linearly (the simplest case!) on frequencies w 1 = ap + bq, w 2 = cp + dp, where a, b, c and d are coefficients (chosen in a way that guarantees that fitnesses are positive). One should be able to see that the resulting equation is mahmetically identical to that describing constant viability selection in a diploid one-locus two-allele population. [Notice that in the diploid case, the induced fitnesses of alleles are w 1 = w 11 p + w 12 q, w 2 = w 21 p + w 22 q. While in the diploid model we always assume w 12 = w 21, there are no a priori reason to set b = c! Homework Problem: rare allele disadvantage. 3. Symmetric frequency-dependent selection against rare alleles can be modeled by assuming that w A = 1 + sp, w a = 1 + tq, where s and t are positive coefficients. Show that in this case depending on the initial allele frequency p the system evolves towards fixation of allele A or allele a. Find a threshold value p separating the domains of attraction of the two fixation states Major points frequency-dependent selection can easily maintain genetic variation; at a polymorphic equilibrium the fitnesses of all genotypes present are identical. 71
2 9.2 Coevolution of two species with applications to mimicry Mimicry is a number of phenomena by which a species imitates another one. There exist numerous examples of mimicry in butterflies, birds, fish, reptiles, amphibians, mammals, moths, beetles, bugs, flies, snails and fungi. There are two general types of mimicry. In Batesian mimicry a palatable species mimics a species (the model) unpalatable to predators. In Müllerian mimicry two (or more) unpalatable species resemble each other. The standard explanation for the resemblance both between Batesian model and mimic and between Müllerian mimics is natural selection by predation. The Batesian mimic increases its fitness by evolving to a phenotypic pattern that the predator tends to avoid. In the Müllerian situation both species improve their protection by evolving to a similar phenotypic pattern that the predator learns to avoid more effectively. We consider two haploid populations of different species with non-overlapping generations. We will concentrate on a single locus with alleles A and a in species 1 and a single locus with alleles B and b in species 2. Let w A, w a, w B and w b be fitnesses (viabilities) of the corresponding morphs, and p 1 and p 2 be the frequencies of morph A within species 1 and morph B within species 2, respectively. Let q i = 1 p i, i = 1, 2. The changes in allele frequencies are described by the standard equations p 1 = p 1 q 1 (w A w a )/w 1, p 2 = p 2 q 2 (w B w b )/w 2, (82a) (82b) where w i is the mean fitness of the i-th species (for example, w 1 = w A p 1 + w a q 1 ). We will assume that the fitness of an individual depends on the genetic composition of its own species as well on the genetic composition of the other species, i.e., w A = w A (p 1, p 2 ), w a = w a (p 1, p 2 ) etc. We will make the symmetry assumption that w A (p 1, p 2 ) = w a (q 1, q 2 ), w B (p 1, p 2 ) = w b (q 1, q 2 ) and use linear (i.e., the simplest) fitness functions. Our symmetry assumption implies that there are no differences in fitness between two morphs of a species not related to their frequencies. These assumptions result in the fitnesses w A = C 1 + ap 1 + bp 2, w a = C 1 + aq 1 + bq 2, w B = C 2 + cp 1 + dp 2, w b = C 2 + cq 1 + dq 2. (83a) (83b) Here C 1 and C 2 are positive constants equal to the fitnesses of morphs A and B (or morphs a and b) when fixed in the populations. Parameters a and d characterize the strength of (indirect) within-species interactions. Positive values of a (d) imply that individuals of a given morph benefit when they are common. Negative values of a (d) imply that individuals of a given morph benefit when they are rare. Parameters b and c characterize the strength of (indirect) between-species interactions. Positive values of b (c) imply that a given morph benefits when its counterpart in another species is common. Negative values of b (c) imply that a given morph benefits when its counterpart in another species is rare. [Add a game-theoretical interpretation of the fitness scheme.] We will assume that selection is weak (in the sense that a, b << C 1, c, d << C 2 ) which allows us to describe the dynamics using differential equations: ṗ 1 = p 1 q 1 [a(p 1 q 1 ) + b(p 2 q 2 )], ṗ 2 = p 2 q 2 [c(p 1 q 1 ) + d(p 2 q 2 )], (84a) (84b) where ṗ i dp i /dt. Note that parameters a, b and c, d in (84) are different from those in (83) by the factors C 1 and C 2, respectively. The main reason for all the simplifying assumptions leading to (84) is that they allow us to study the system analytically. Before analyzing the whole system, it is illuminating to consider several special cases. 72
3 a. Only a single species is allowed to evolve. If allele frequency in a species, say species 1, is somehow fixed (p 1 = const), then ṗ 2 = p 2 q 2 [c(p 1 q 1 ) + d(p 2 q 2 )] This equation is similar to equation (??) analyzed above. b. No between species interactions. If b = c = 0, the resulting equation is similar to equation (??) analyzed above. c. Fitness depends only on the genetic composition of another species. This may be the case in many coevolving systems including host-parasite systems, competitive and cooperative interactions. In this case, dp 1 dt = p 1q 1 b(p 2 q 2 ), dp 2 dt = p 2q 2 c(p 1 q 1 ). One can get some qualitative ideas about the dynamics of this system by inspecting the right-hand sides of the equations. For example, if both b and c are positive, p 1 increases (dp 1 /dt > 0) or decreases (dp 1 /dt < 0) depending on whether p 2 is large or smaller than 0.5. Similar conclusions can be made about p 2 and for other choices of the signs of parameters b and c. This knowledge will be used below after we find the exact solutions. We will use the separation of variables method. Dividing (85a) by (85b) one finds that Separating the variables dp 1 = bp 1q 1 (p 2 q 2 ) dp 2 cp 2 q 2 (p 1 q 1 ). c p 1 q 1 p 1 q 1 dp 1 = b p 2 q 2 p 2 q 2 dp 2, and integrating, we finds that all solutions of (85) satisfy to (85a) (85b) p1 q 1 c p 1 q 1 dp 1 = b p2 q 2 p 2 q 2 dp 2. (86) Because x (1 x) x(1 x) dx = equation (86) can be integrated to 2x 1 d[x(1 x)] du x(1 x) dx = = = ln u, where u x(1 x), x(1 x) u c ln(p 1 q 1 ) = b ln(p 2 q 2 ) + const. Thus, all solutions of (85) satisfy to the following equality (p 1 q 1 ) c = const. (87) (p 2 q 2 ) b Now, we can combine our knowledge of the exact trajectories (given by equation (87)) with the qualitative conclusions about the areas in the phase-plane (p 1, p 2 ) where allele frequencies increase or decrease (see above). There are three cases to consider. (i) Parameters b and c are positive. In this case the fitness of each morph increases with an increase in the frequency of its counterpart in the other species. This choice of parameters can be interpreted as describing cooperative interactions between the species. (ii) Parameters b and c are negative. In this case the fitness of each morph decreases with an increase in the frequency of its counterpart in the other species. This 73
4 choice of parameters can be interpreted as describing competitive interactions between the species. (iii) Parameters b and c have different signs. In this case the fitness of a morph in one species increases with an increase in the frequency of its counterpart in the other species, but the fitness of the counterpart decreases with an increase in the frequency of the first morph. This choice of parameters can be interpreted as describing victim-exploiter type interactions between the species. In case (i), the system evolves to a state with both species monomorphic (there are two such states: p 1 = p 2 = 0 and p 1 = p 2 = 1). In case (ii), the system evolves to a state with both species monomorphic (there are two such states: p 1 = 0, p 2 = 1 and p 1 = 1, p 2 = 0). In case (iii), the solutions are neutrally stable periodic orbits encircling the polymorphic equilibrium (1/2, 1/2); any perturbation moves the system to a different periodic orbit. The dynamic behavior of (84) is similar to the behavior of the classical Lotka-Volterra predator-prey system (e.g., Hofbauer & Sigmund, 1988) Major points With linear frequency dependent fitnesses and no within-species interactions cooperative and competitive systems do not maintain genetic variation; victim-exploiter type systems maintain genetic variation; these systems exhibit cycling. i. ii. iii. p2 Figure 8: Dynamic regimes with no p1 within-species interactions. (i) Cooperative p1 between-species interactions. p1 (ii) Competitive between-species interactions. (iii). Victim-exploiter type between-species interactions. The two species whose coevolutionary dynamics are described by (84) can be considered as Batesian mimic and model or as two species belonging to the same Müllerian ring. Morph A is similar to morph B and morph a is similar to morph b. Biologically the symmetry assumption means that we include only the effects of mimicry in our model and no costs or other aspects of fitness. Hence morphs within species are alike in the sense that there is no difference between A and a (B and b) except that they look similar to B and b (A and a), respectively. Positive values of a and/or d imply that individuals of a given morph benefit when they are common. This is the case when species 1 and/or 2 is unpalatable. Negative values of a and/or d imply that individuals of a given morph benefit when they are rare. This is the case when species 1 and/or 2 is palatable. Parameters b and c characterize the strength of indirect between-species interactions. These interactions are mediated through the predator and thereby indirect. Positive values of b and/or c imply that a given morph benefits when its counterpart in another species is common. This is the case when the other species is unpalatable. Negative values of b and/or c imply that a given morph benefits when its counterpart in another species is rare. This is the case when the other species is palatable or it is unpalatable but has a smaller degree of unpalatability (a weaker Müllerian mimic). 74
5 Our results on dynamic behavior of (84) can be applied to understand the evolutionary dynamics of different mimicry systems. We will consider three different parameter configurations: Classical Müllerian mimicry Classical Batesian mimicry Two unpalatable species have different abundances and/or different degrees of unpalatability Parameter Interpretation configuration a, b, c, d > 0 Within each species a morph benefits when its own frequency or the frequency of its counterpart in the other species increases a, c < 0; b, d > 0 Both species 1 (palatable mimic) and species 2 (unpalatable model) benefit from increasing the model frequency and suffer from increasing the mimic frequency a, b, d > 0, c < 0 Within each species a morph benefits when its own frequency increases. Species 1 (the weaker mimic) also benefits when the frequency of its stronger counterpart (from species 2) increases, but the stronger mimic suffers when its weaker counterpart increases in frequency - Homework Problem: stability of equilibria of the coevolutionary model (84). Find all equilibria of the coevolutionary model (84). Summarize the conditions for existence and stability of these equilibria in a table: Equilibrium (p 1, p 2 ) Conditions for existence Conditions for stability Use Maple (or your favorite software) to plot phase portraits corresponding to the situations when the equilibria are stable and unstable. Can genetic variation be maintained in the system? Are there any situations when none of the equilibria are stable? - 75
6 9.3 Linear frequency-dependent selection in a diploid population We consider a deterministic model of a large randomly mating diploid population with discrete generations concentrating on a single diallelic locus with alleles A and a. Let w AA, w Aa and w aa be the fitnesses (viabilities) of genotypes AA, Aa and aa, respectively, and p be the frequency of allele A, with q = 1 p. With random mating the population is in Hardy-Weinberg proportions so that the frequencies of the three genotypes are p 2, 2pq and q 2. The change in p in one generation is described by the standard equation p = p(w A w). (88) w Here w A = pw AA + qw Aa is the average fitness of allele A and w = p 2 w AA + 2pqw Aa + q 2 w aa is the mean fitness of the population. If the fitnesses are constant, the population gradually evolves to a polymorphic equilibrium (with 0 < p < 1) provided there is overdominance (i.e., if w Aa > w AA, w aa ) or to a fixation state (with p = 0 or p = 1), otherwise. We assume that selection is frequency-dependent and will use linear (that is the simplest) functions to model the frequency-dependence: w AA = p 2 W pqW 12 + q 2 W 13, w Aa = p 2 W pqW 22 + q 2 W 23, w aa = p 2 W pqW 32 + q 2 W 33, where W ij are parameters (i, j = 1, 2, 3) characterizing the extent to which changes in the frequencies of three genotypes influence their fitnesses. Asmussen and Basnayake (1990): a b b d c b, c b a where all coefficients are positive. Dynamics: multiple stable equilibria; plenty of opportunities for the maintenance of variation. Altenber (1991): a (a + c)/2 c b b b, c (a + c)/2 a where some coefficients can be negative. Gavrilets and Hastings (1995): δ β α γ η γ α β δ Under this symmetric model, the dependence of fitness of a heterozygote w Aa on the allele frequency p is described by a (quadratic) function symmetric about 1/2, while w AA and w aa considered as (quadratic) functions of p are reflections of each other about 1/2. For this symmetric model to produce feasible (i.e., non-negative) fitnesses, one has to assume that α, γ, δ > 0, β > αδ, η > γ. [Note that β and η are allowed to be negative.] The dynamics are not changed if fitnesses are multiplied or divided by a constant. This allows one to assume without loss of generality that δ = 1. For the symmetric model, the dynamic equation (106) can be represented as p =. pq(p q)(1 γ Ωpq), (89) w 76
7 where the mean fitness w can be written as w = 1 2(2 β γ)pq + 2Ωp 2 q 2 with Ω = 1 + α 2β 2γ + 2η. In the Altenberg model, Ω = 0. The allele frequency p does not change if p = 0, q = 0, p = q or pq = (1 γ)/ω. Thus, equation (107) always has two monomorphic equilibria at p = 0 and p = 1 and a polymorphic equilibrium at p = 1/2. If 0 < (1 γ)/ω < 1/4, it has two additional polymorphic equilibria with allele frequencies satisfying p(1 p) = (1 γ)/ω. We shall denote these equilibria p and p +. An equilibrium of (107) is stable if the corresponding eigenvalue lies between -2 and 0. These eigenvalues can be found in a straightforward manner. 8 7 γ>1 6 5 USU C UUU USUSU γ<1 5 C 4 SUSUS 3 2 SUUUS SUS beta Figure 9: Areas in parameter space corresponding to different patterns of existence and stability of equilibria in a model of linear frequency-dependent selection. Only areas with C 0 are shown. Figure 14 summarizes conditions for existence and stability of different equilibria in terms of γ, β and a parameter C that combines several parameters. The parameter C is defined as C = c 1 + c 2 with c 1 α β 2 and c 2 2(γ + η) or, alternatively, as C = Ω + 4γ (1 β) 2. Parameter c 2 is the minimal possible value of the fitness of heterozygote. If β < 1, c 1 determines the minimal possible value of the fitness of homozygotes, c 1 /(c 1 + (1 β 2 )), while if β > 1, c 1 determines the maximum possible value of the fitness of homozygotes, c 1 + β 2. Thus, parameter C characterizes the overall strength of selection. For fitnesses to be feasible c 2 must be non-negative and c 1 must be non-negative if β < 0. These conditions imply that C must be larger or equal to 0 if β is negative, and must be larger or equal to β 2 if β is positive. Figure 1 shows areas in parameter space corresponding to different patters of existence and stability of equilibria in a model of linear frequency-dependent selection. Each pattern is described by a string of S s (for stable) and U s (for unstable). The leftmost, the middle and the rightmost entries indicate the stability of the monomorphic equilibrium at p = 0, the polymorphic equilibrium at p = 1/2 and the monomorphic equilibrium at p = 1, while the remaining entries (if any) indicate the stability of the polymorphic equilibria p and p +. The left parabola is described by equation C = (β + 5)(β + 1). The right parabola is described by equation C = (β + 1)(β 3). - Homework: Verify the conditions for existence and stability of different equilibria presented in Figure 2. - Figure 14 shows that the system can have up to three different stable equilibria simultaneously, that a polymorphic 77
8 equilibrium can be stable simultaneously with two monomorphic equilibria, and that two different polymorphic equilibria can be stable simultaneously. Simultaneous stability of different equilibria implies that the outcome of evolution strongly depends on the initial conditions and population history. If parameters change in such a way that the system moves from one area to another, the dynamic system undergoes a bifurcation. For example, a change from USU to USUSU corresponds to a pitchfork bifurcation. Figure 14 also shows that there are two areas with non-standard patterns of stability of equilibria. In the first area (marked UUU), none of the three equilibria (two monomorphic and one polymorphic at p = 1/2) are stable. If parameters change in such a way that the system moves from from USU to UUU, the dynamic system undergoes a period-doubling bifurcation. In the second area (marked SUUUS), the two monomorphic equilibria are stable, while none of the three polymorphic equilibria are stable. Numerical iterations of (107) with parameter values corresponding to these areas reveal a variety of complex dynamic behaviors (e.g., cycles and chaos that arises via period-doubling route) similar to those observed in classical ecological models (e.g., May, 1974, 1976; May and Oster, 1976; Hastings et al., 1993). Figure 14 shows that sufficient conditions for the complex dynamic behavior to occur is sufficiently strong overall selection (i.e., small C) and sufficiently strong deleterious effect of heterozygote on homozygotes (i.e., β < 1). 1.0 allele frequency a allele frequency b generation number Figure 10: Dynamics of allele frequency with parameter values from the area marked SUUUS in Figure 1. Both monomorphic equilibria (p = 0 and p = 1) are stable to small perturbations. a) β = 1.002, γ =.999, C = 0. Three trajectories are shown with initial allele frequencies p(0) =.09, p(0) =.55 and p(0) =.92. b) β = 1.001, γ =.9985, C = 0, p(0) =.55. There are also two unusual types of behavior described in Figure 15. In Figure 15a depending on the initial conditions the population evolves to a fixation state or remains polymorphic indefinitely. In the latter case, the gradual changes in the allele frequency towards p = 1/2 are interrupted by apparently chaotic fluctuations that move p away from 1/2. Such alterations between an apparently deterministic behavior and apparently chaotic fluctuations, repeated at apparently random intervals, are called intermittency (Pomeau and Manneville, 1979; Olsen and Degn, 1985). In Figure 15b, these alterations end at some time point with the population settling down to a monomorphic state. In this case the system exhibits transient chaos (Grebogi et al., 1983; Tél, 1990). In the examples presented in Figure 15, the deterministic phase of the dynamics can last for hundreds of generations, the chaotic phase is extremely short and transient chaos (in Figure 15b) persists for a very long time. In general, the durations of all 78
9 these stages depend both on the parameter values and initial allele frequency (see below the discussion of Figure 16). Both intermittency and transient chaos are known to occur in various dynamical models including the simplest nonlinear model, the single logistic map. Our model allows a very simple explanation of the mechanisms underlying intermittency and transient chaos using the graphical cobwebbing method (May and Oster, 1976). 1.0 a p + max b p + max p min p p min p p Figure 11: Graphs of the allele frequency in the next generation, p = p + p, as function of the allele frequency at this generation, p. Also shown are the diagonal and the lines corresponding to the unstable polymorphic equilibria p and p +. Parameter values are γ =.9, C = 0, β = 1.15 (in Figure 15a) and β = 1.05 (in Figure 15b). Figure 16 shows graphs of the allele frequency in the next generation, p = p+ p, as function of the allele frequency at this generation, p. Also shown are the diagonal and the lines corresponding to the unstable polymorphic equilibria p and p +. First note that for p in the neighborhood of p or p +, the graph of p lies very close to the diagonal (at which p = p). That means that in these areas the changes in the allele frequency p are very small. For p < p, p < p, and p slowly moves towards fixation of allele a. For p > p +, p > p, and the allele frequency slowly moves towards fixation of allele A. For p values slightly larger than p, p is slightly larger than p, and, thus, p slowly moves towards 1/2. In the neighborhood of p = 1/2, however, the dynamics are chaotic as suggested by the fact that the slope of the graph of p at p = 1/2 is smaller than -1. It takes many generations to leave the neighborhood of p or p + and once the system has left this neighborhood, the dynamics can abruptly become fully chaotic in the neighborhood of p = 1/2, only to get caught in the neighborhood of p or p + again, sooner or later. Between 0 and 1, the graph of p has a minimum, marked min, and a maximum, marked max. During the chaotic phase the allele frequency remains between these points that represent the boundaries of the chaotic attractor. In Figure 15a these boundaries lie closer to 1/2 than the unstable equilibria and the allele frequency cannot cross the values p + and p. In this case the intermittent chaos in the system (in the form similar to that one in Figure 15a) is present forever. A different situation is described in Figure 15b, where the boundaries of the chaotic attractor lie further from 1/2 than the unstable equilibria p and p +. Now the allele frequency can cross the values p + and p during the chaotic phase. Once this has happened, p slowly evolves to a fixation state (as in Figure 15b). The situation when the boundaries of the chaotic attractor coincide exactly with the unstable equilibria is called a crisis (Grebogi et al., 1983). In the model considered here, if C = 0, the crisis occurs when γ β + 2. Note that in general, the dependence of the length of chaotic transients on the system parameter is proportional to (a a c ) b, where a is the parameter value, a c is the parameter value at which the crisis occurs, and b is the critical exponent, which is equal to 0.5 for a broad class of one-dimensional systems (Grebogi et al., 1987). The parameter values for Figure 15 and Figure 17 below were chosen to result in long transients. The parameters for Figure 16 were chosen slightly different from those for Figure 15 in order to produce a smoother graph of p as function of p. Numerical analysis of (107) has also shown, perhaps surprisingly, that complex dynamic behavior occurs even 79
10 allele frequency allele frequency generation number Figure 12: Dynamics of allele frequency with parameter values from the area marked SUSUS in Figure 14. Both monomorphic equilibria (p = 0 and p = 1) and the polymorphic equilibrium at p = 1/2 are stable to small perturbations. a) β = 5, γ =.9, C =.1, p(0) =.55. b) β = 5, γ =.5, C =.1, p(0) =.55. outside areas marked UUU and SUUUS in Figure 14. In the areas marked USU and SUSUS, cycles and chaos can not only exist simultaneously with stable equilibria, but the former can have much larger domains of attraction than the latter (see Figure 17). For parameter values used in computing the dynamics in Figure 17, both monomorphic equilibria and the polymorphic equilibrium at p = 1/2 are stable to small perturbations. In Figure 17a, the growing regular oscillations in allele frequency are interrupted by apparently chaotic fluctuations that move p back to the neighborhood of 1/2, i.e., one observes intermittency. In Figure 17b, the alterations between apparently deterministic behavior and apparent chaos end at some time point with the population settling down to a polymorphic state at p = 1/2, i.e., one observes transient chaos. The graphical cobwebbing method can be used to understand these kinds of behavior as well. A necessary condition for dramatic changes in allele frequency described in Figures 15 and 17 is strong (at least occasionally) selection. If selection is very weak (i.e., if the differences among coefficients α, β, γ, δ and η are very small) then the difference equation (107) can be approximated by the corresponding differential equation and the only possible outcome of the dynamics is gradual evolution towards an equilibrium. 80
11 9.4 Cycling in systems of ordinary differential equations Let us consider a system of two ordinary differential equations ẋ = f(x, y), ẏ = g(x, y), (90a) (90b) defined on an open set G R Dulac s criterion Theorem. If there exists a continuously differentiable function h(x, y) such that (hf) x + (hg) y is continuous and non-zero on some simply connected domain D G, then no periodic orbit can lie entirely in D. Example. Consider the Lotka-Volterra model with h = 1/(xy) Poincaré-Bendixon theorem ẋ = x(a a 1 x + b 1 y), ẏ = y(b a 2 y + b 2 x), Theorem. Assume that a trajectory γ(x 0, y 0 ) enters and does not leave some closed domain D and that there are no equilibrium points in D. Then there is at least one periodic orbit in D, and this orbit is in the ω-limit set of (x 0, y 0 ). To apply the Poincaré-Bendixon theorem we need to find a region D which contains no equilibrium points and which trajectories enter but do not leave. Both the Poincaré-Bendixon theorem and Dulac s criterion do not work in higher dimensions (n > 2) Non-existence of periodic orbits Theorem (Hofbauer&Sigmund). If the partial derivatives f/ x 0 and g/ y 0 for all (x, y) G, cycling is impossible. The same result holds if f/ x 0 and g/ y 0 for all (x, y) G. No cycling means the trajectories converge either to an equilibrium or to infinity. Example. For the coevolutionary model (85), ṗ 1 / p 2 = 2bp 1 q 1, ṗ 2 / p 1 = 2cp 2 q 2. Thus, if b and c have the same sign so do ṗ 1 / p 2 and ṗ 2 / p 1, and stable cycling is impossible Poincaré-Andronov-Hopf bifurcation (birth of cycles from an equilibrium) Example. Let us consider a system of non-linear ODE: ẋ = αx βy (x 2 + y 2 )x, ẏ = βx + αy (x 2 + y 2 )y. (91a) (91b) 81
12 The origin (0, 0) is an equilibrium; the corresponding stability matrix is ( ) α β S = β α with eigenvalues λ = α ± iβ. Thus, (0, 0) is a stable focus for α < 0 and an unstable focus for α > 0. Using polar coordinates x = r cos θ, y = r sin θ the system can be rewritten as ṙ = αr r 3, θ = β. Thus, if β 0, θ changes with constant speed β. ṙ = 0 if r = 0 or r = α. If α > 0, ṙ > 0 for 0 < r < α and ṙ < 0 for r > α. In this case, there is a periodic solution which is stable. If α < 0, ṙ < 0 always and the only (stable) equilibrium is r = 0. Feb.2. General case. Theorem. Let us consider a system of n ordinary differential equations ż = F µ (z), (92) depending on some parameter µ and defined on an open subset of R n. Let z be an equilibrium of (92). Assume that all eigenvalues of the Jacobian have negative real parts, with the exception of one pair of complex conjugate eigenvalues α(µ) ± iβ(µ). Let for some µ 0, these two eigenvalues be pure imaginary: α(µ 0 ) = 0, β(µ 0 ) 0, and let d dµ α(µ 0) 0. Then for µ sufficiently close to µ 0, a periodic solution bifurcates from equilibrium z. If the first Lyapunov value L 1 (defined below) is negative, the bifurcation is supercritical (i.e. the periodic orbit exists for α(µ) > 0 and is stable). If L 1 is positive, the bifurcation is subcritical (i.e. the periodic orbit exists for α(µ) < 0 and is unstable). For small values of α(µ), the period and the amplitude of the periodic orbit are approximately 2π/ β(µ) and α(µ), respectively. First Lyapunov value L 1. Let us consider a system of two ordinary differential equations dx =ax + by + P (x, y), dt (93a) dy =cx + dy + Q(x, y). dt (93b) Let TrS = a + d and det S = ad bc be the trace and the determinant of the stability matrix for the equilibrium (0,0). Let P (x) and Q(x) be represented as power series P (x, y) =P 2 (x, y) + P 3 (x, y) +..., Q(x, y) =Q 2 (x, y) + Q 3 (x, y) +..., (94a) (94b) where P 2 (x, y) =a 20 x 2 + a 11 xy + a 02 y 2, P 3 (x, y) =a 30 x 3 + a 21 x 2 y + a 12 xy 2 + a 03 y 3, Q 2 (x, y) =b 20 x 2 + b 11 xy + b 02 y 2, Q 3 (x, y) =b 30 x 3 + b 21 x 2 y + b 12 xy 2 + b 03 y 3. (95a) (95b) (95c) (95d) 82
13 The first Lyapunov value L 1 is defined as (Bautin 1938). L 1 = π 1 4(det S) 3/2 b [ ac(a a 11 b 02 + a 02 b 11 ) + ab(b a 20b 11 + a 11 b 20 ) + c 2 (a 11 a a 02 b 02 ) 2ac(b 2 02 a 20 a 02 ) 2ab(a 2 20 b 20 b 02 ) b 2 (2a 20 b 20 + b 11 b 20 ) + (bc 2a 2 )(b 11 b 02 a 11 a 20 ) (a 2 + bc) (3(cb 03 ba 30 ) + 2a(a 21 + b 12 ) + ca 12 bb 21 ) ] The expression for L 1 looks horrible and almost impossible to use practically. However, using modern symbolic manipulation packages such as Maple and Mathematica makes it a trivial exercise. Example 1. The first Lyapunov value corresponding to the equilibrium (0, 0) of system (91) is L 1 = 8. Thus, the bifurcation is supercritical (i.e. the periodic orbit exists for α(µ) > 0 and is stable) which is what we already found using exact analysis. Example 2. Let us consider a system of two ordinary differential equations x ẋ = rx(1 x/k) β x + C y, x ẏ = (kβ x + C m)y, where all coefficients are positive. One can interpret x as the density of a prey population which grows logistically in the absence of predator (if y = 0, then ẋ = rx(1 x/k)), and y as the density of predator which experiences densityindependent mortality (with rate m). Terms β x x+c y and kβ x x+c y stand for the rates at which predation decreases the prey population and increases the predator population, respectively (a Holling type II functional response). The trivial equilibrium (0, 0) is always unstable. [It is a saddle point with eigenvalues r and m.] The equilibrium with only the prey population present (K, 0) is stable if m > m where m = Kkβ K+C, and is unstable otherwise. [The eigenvalues are r and m m.] The model also has a single coexistence state with population densities x = mc kβ m, y = rkc(k + C)(m m) K(kβ m) 2. This equilibrium is feasible if m < m (predator s mortality is not too large). The trace and determinant of the corresponding stability matrix are mr(kkβ Km Ckβ mc) T rs =, kβ(kβ m)k det S = (K + C)(m m)mr, Kkβ that is det S > 0 always, but the TrS can change its sign. The coexistence state is stable if m > m kβ(k C)/(K + C) and is unstable otherwise. [Note that m < m.] If m = m, the system undergoes a Poincaré- Andronov-Hopf bifurcation. For this bifurcation, the first Lyapunov value is 4Cβk(K C)r2 L 1 = K 2 (K + C) 2, 83
14 and is positive (because at m = m, K/C = (kβ + m)/(kβ m) > 1). Thus, the bifurcation is supercritical (the periodic orbit exists for m < m and is stable). Summarizing, if m > m, the predator population cannot maintain itself whereas if m < m both species coexist. If m < m < m, the coexistence is non-oscillatory whereas if m < m, there are stable oscillations of the densities of both populations. Example 3. For the coevolutionary model (84), the first Lyapunov value corresponding to Poincaré-Andronov- Hopf bifurcation from the doubly polymorphic equilibrium (1/2, 1/2) is L 1 = 2a(b + c), b which can be positive or negative depending on parameter values. Let us consider the case corresponding to classical Batesian mimicry (a < 0, b > 0, c < 0, d > 0). Assume that a < b, d < c (i.e. between species interactions are stronger than within species interactions). These conditions guarantee that (i) monomorphic equilibria are unstable, (ii) singly polymorphic equilibria do not exist, and (iii) the determinant of the stability matrix evaluated at the doubly polymorphic equilibrium (1/2, 1/2) is positive. Thus, if a + d < 0, the equilibrium (1/2, 1/2) is a stable focus. If a + d > 0, this equilibrium is an unstable focus. The sign of L 1 coincides with that of b+c. Thus, for small a+d, the system should have a stable limit cycle if a+d > 0, b+c < 0. If a + d < 0, b + c > 0, then there is an unstable cycle encircling the equilibrium (1/2, 1/2). Numerical examples. Stable cycle: a = 1, b = 1.2, c = 1.3, d = Unstable cycle: a = 1.01, b = 1.3, c = 1.2, d = 1. Illuminating initial conditions are p 1 = 5, p 2 = 0.75 and p 1 = 0.5, p 2 =. 84
A simple genetic model with non-equilibrium dynamics
J. Math. Biol. (1998) 36: 550 556 A simple genetic model with non-equilibrium dynamics Michael Doebeli, Gerdien de Jong Zoology Institute, University of Basel, Rheinsprung 9, CH-4051 Basel, Switzerland
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
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 informationProblems for 3505 (2011)
Problems for 505 (2011) 1. In the simplex of genotype distributions x + y + z = 1, for two alleles, the Hardy- Weinberg distributions x = p 2, y = 2pq, z = q 2 (p + q = 1) are characterized by y 2 = 4xz.
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 informationCoevolutionary Chase in Two-species Systems with Applications to Mimicry
J. theor. Biol. (1998) 191, 415 427 Coevolutionary Chase in Two-species Systems with Applications to Mimicry SERGEY GAVRILETS* AND ALAN HASTINGS *Departments of Ecology and Evolutionary Biology and Mathematics,
More informationA SIMPLE MATHEMATICAL MODEL FOR BATESIAN MIMICRY
A SIMPLE MATHEMATICAL MODEL FOR BATESIAN MIMICRY TERENCE R. BLOWS AND BARRY J. WIMMER Received 24 September 2003 A simple mathematical model is presented for Batesian mimicry, which occurs when a harmless
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
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: January 14, 2019, at 08 30 12 30 Johanneberg Kristian
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationMath 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations
Math 345 Intro to Math Biology Lecture 7: Models of System of Nonlinear Difference Equations Junping Shi College of William and Mary, USA Equilibrium Model: x n+1 = f (x n ), here f is a nonlinear function
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 informationSelection Page 1 sur 11. Atlas of Genetics and Cytogenetics in Oncology and Haematology SELECTION
Selection Page 1 sur 11 Atlas of Genetics and Cytogenetics in Oncology and Haematology SELECTION * I- Introduction II- Modeling and selective values III- Basic model IV- Equation of the recurrence of allele
More informationHomework 2 Modeling complex systems, Stability analysis, Discrete-time dynamical systems, Deterministic chaos
Homework 2 Modeling complex systems, Stability analysis, Discrete-time dynamical systems, Deterministic chaos (Max useful score: 100 - Available points: 125) 15-382: Collective Intelligence (Spring 2018)
More information6 Price equation and Selection in quantitative characters
6 Price equation and Selection in quantitative characters There are several levels of population description. At the most fundamental level, e describe all genotypes represented in the population. With
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationAndronov Hopf and Bautin bifurcation in a tritrophic food chain model with Holling functional response types IV and II
Electronic Journal of Qualitative Theory of Differential Equations 018 No 78 1 7; https://doiorg/10143/ejqtde018178 wwwmathu-szegedhu/ejqtde/ Andronov Hopf and Bautin bifurcation in a tritrophic food chain
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 informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
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 informationMULTISTABILITY IN A BUTTERFLY FLOW
International Journal of Bifurcation and Chaos, Vol. 23, No. 12 (2013) 1350199 (10 pages) c World Scientific Publishing Company DOI: 10.1142/S021812741350199X MULTISTABILITY IN A BUTTERFLY FLOW CHUNBIAO
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
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 informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
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 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 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 information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
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 informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More information2 One-dimensional models in discrete time
2 One-dimensional models in discrete time So far, we have assumed that demographic events happen continuously over time and can thus be written as rates. For many biological species with overlapping generations
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 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 informationChaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB
Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the
More informationEvolution Module. 6.2 Selection (Revised) Bob Gardner and Lev Yampolski
Evolution Module 6.2 Selection (Revised) Bob Gardner and Lev Yampolski Integrative Biology and Statistics (BIOL 1810) Fall 2007 1 FITNESS VALUES Note. We start our quantitative exploration of selection
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationOutline of lectures 3-6
GENOME 453 J. Felsenstein Evolutionary Genetics Autumn, 007 Population genetics Outline of lectures 3-6 1. We want to know what theory says about the reproduction of genotypes in a population. This results
More informationCoevolutionary Chase in Exploiter Victim Systems with Polygenic Characters
J. theor. Biol. (1997) 186, 57 534 Coevolutionary Chase in Exploiter Victim Systems with Polygenic Characters SERGEY GAVRILETS* Departments of Ecology and Evolutionary Biology and Mathematics, University
More informationCHAOS FROM LINEAR FREQUENCY-DEPENDENT SELECTION
1 CHAOS FROM LINEAR FREQUENCY-DEPENDENT SELECTION LEE ALTENBERG Department of Zoology, Duke University, Durham, NC 27706 Submitted February 13, 1989; Revised March 9, 1990; Accepted May 20, 1990 1 Abstract.
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 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 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 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: January 08, 2018, at 08 30 12 30 Johanneberg Kristian
More informationChapter 4. Transition towards chaos. 4.1 One-dimensional maps
Chapter 4 Transition towards chaos In this chapter we will study how successive bifurcations can lead to chaos when a parameter is tuned. It is not an extensive review : there exists a lot of different
More information2 Discrete growth models, logistic map (Murray, Chapter 2)
2 Discrete growth models, logistic map (Murray, Chapter 2) As argued in Lecture 1 the population of non-overlapping generations can be modelled as a discrete dynamical system. This is an example of an
More informationChapitre 4. Transition to chaos. 4.1 One-dimensional maps
Chapitre 4 Transition to chaos In this chapter we will study how successive bifurcations can lead to chaos when a parameter is tuned. It is not an extensive review : there exists a lot of different manners
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 informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
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 information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
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 informationPHY411 Lecture notes Part 4
PHY411 Lecture notes Part 4 Alice Quillen February 1, 2016 Contents 0.1 Introduction.................................... 2 1 Bifurcations of one-dimensional dynamical systems 2 1.1 Saddle-node bifurcation.............................
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 informationIntroduction to Dynamical Systems Basic Concepts of Dynamics
Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
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 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 informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A
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 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 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 informationBoulder School for Condensed Matter and Materials Physics. Laurette Tuckerman PMMH-ESPCI-CNRS
Boulder School for Condensed Matter and Materials Physics Laurette Tuckerman PMMH-ESPCI-CNRS laurette@pmmh.espci.fr Dynamical Systems: A Basic Primer 1 1 Basic bifurcations 1.1 Fixed points and linear
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 informationOutline of lectures 3-6
GENOME 453 J. Felsenstein Evolutionary Genetics Autumn, 009 Population genetics Outline of lectures 3-6 1. We want to know what theory says about the reproduction of genotypes in a population. This results
More informationSummary of topics relevant for the final. p. 1
Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
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: April 06, 2018, at 14 00 18 00 Johanneberg Kristian
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationBIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs
BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange
More information(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,
2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations
More informationSimplest Chaotic Flows with Involutional Symmetries
International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries
More informationMath 1060 Linear Algebra Homework Exercises 1 1. Find the complete solutions (if any!) to each of the following systems of simultaneous equations:
Homework Exercises 1 1 Find the complete solutions (if any!) to each of the following systems of simultaneous equations: (i) x 4y + 3z = 2 3x 11y + 13z = 3 2x 9y + 2z = 7 x 2y + 6z = 2 (ii) x 4y + 3z =
More information... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré
Chapter 2 Dynamical Systems... it may happen that small differences in the initial conditions produce very great ones in the final phenomena. Henri Poincaré One of the exciting new fields to arise out
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
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 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 informationLotka-Volterra Models Nizar Ezroura M53
Lotka-Volterra Models Nizar Ezroura M53 The Lotka-Volterra equations are a pair of coupled first-order ODEs that are used to describe the evolution of two elements under some mutual interaction pattern.
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationSpring /30/2013
MA 138 - Calculus 2 for the Life Sciences FINAL EXAM Spring 2013 4/30/2013 Name: Sect. #: Answer all of the following questions. Use the backs of the question papers for scratch paper. No books or notes
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 informationDerivation of Itô SDE and Relationship to ODE and CTMC Models
Derivation of Itô SDE and Relationship to ODE and CTMC Models Biomathematics II April 23, 2015 Linda J. S. Allen Texas Tech University TTU 1 Euler-Maruyama Method for Numerical Solution of an Itô SDE dx(t)
More informationZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS
This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]
More informationUnit Ten Summary Introduction to Dynamical Systems and Chaos
Unit Ten Summary Introduction to Dynamical Systems Dynamical Systems A dynamical system is a system that evolves in time according to a well-defined, unchanging rule. The study of dynamical systems is
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 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 informationComputational Aspects of Aggregation in Biological Systems
Computational Aspects of Aggregation in Biological Systems Vladik Kreinovich and Max Shpak University of Texas at El Paso, El Paso, TX 79968, USA vladik@utep.edu, mshpak@utep.edu Summary. Many biologically
More informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
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 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 informationSolution to Homework #5 Roy Malka 1. Questions 2,3,4 of Homework #5 of M. Cross class. dv (x) dx
Solution to Homework #5 Roy Malka. Questions 2,,4 of Homework #5 of M. Cross class. Bifurcations: (M. Cross) Consider the bifurcations of the stationary solutions of a particle undergoing damped one dimensional
More informationNonlinear Dynamics and Chaos
Ian Eisenman eisenman@fas.harvard.edu Geological Museum 101, 6-6352 Nonlinear Dynamics and Chaos Review of some of the topics covered in homework problems, based on section notes. December, 2005 Contents
More informationIntroduction to the Phase Plane
Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called
More informationFunctional divergence 1: FFTNS and Shifting balance theory
Functional divergence 1: FFTNS and Shifting balance theory There is no conflict between neutralists and selectionists on the role of natural selection: Natural selection is the only explanation for adaptation
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 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 information