Calculating the coevolved dispersal distances
|
|
- Lucinda Hampton
- 5 years ago
- Views:
Transcription
1 A Calculating the coevolved dispersal distances The calculations for the baseline model monoculture evolution in the presence of moderate environmental variation were presented in Snyder Here I show how to calculate the coevolved dispersal distances. To find the coevolved dispersal distances, we find the evolutionarily stable ESS annual dispersal distance as a function of the perennial dispersal distance, then find the ESS perennial dispersal distance as a function of the annual dispersal distance. At the intersection of these two curves, neither species will be driven to adjust its dispersal strategy in response to the other s strategy, for each species is on its evolutionarily stable curve: this point is an evolutionarily stable coalition. Finally, we need to verify that evolutionary dynamics cause the system to converge to this point: i.e., that the intersection is not only evolutionarily stable but convergently stable, in the sense of Geritz et al To find, e.g., the ESS annual dispersal distance as a function of the perennial dispersal distance, consider an annual invader in the presence of annual and perennial residents. The invader is a mutant with dispersal distance a little larger or smaller that of the annual resident. If the invader has a positive long-run growth rate, it will replace the former resident and become subject to invasion attempts by new mutants. The process stops when the best strategy the annual invader could adopt is that of the annual resident i.e., invader growth is maximized when the invader dispersal distance equals the resident dispersal distance. Signifying the invader long run growth rate by r i d a,d p,d i, d a is the resident annual distance, d p is the resident perennial distance, and d i is the invader distance, the ESS annual dispersal distance d a is given by the condition r i d a,d p,d i d i = 0. 1 di =d a=d a Note that d a depends on d p, the perennial distance. We see that finding the ESS dispersal distance requires us to find an expression for the long-run invader growth rate r i. Invader growth is affected by competition from the residents, which will in turn be determined by the spatiotemporal distribution of the resident populations. Let us first find an expression for r i in terms of the resident populations and then find expressions for the population distributions. Before dispersal, the local population of species j at time t + 1 is equal to the local growth rate λ j x,t times the local population at time t: n j x,t + 1 = λ j x,tn j x,t. To characterize the regional dynamics, we take the spatial average of both sides: n j x t + 1 = λ j n j x t, x denotes a spatial average. Expressing the dynamics in this form is not very useful because the righthand side is not a function of n j x. However, we can rewrite the average of λ j times n j as the product of their averages plus their covariance, allowing us to factor out n j x. Thus, n j x t + 1 = λ j x t n j x t + Covλ j,n j x t = λ j x + Covλ j,n j / n j x x n j x t. 2 The quantity in square brackets is the regional growth rate the factor by which the spatially averaged population grows or shrinks in a year. We denote the regional growth 1
2 rate by λ j and rewrite it as λ j x + Covλ j,ν j x, ν j x,t = n j x,t/ n j x is referred to as the relative population density Chesson, The regional growth rate fluctuates in time, and so the long-run growth of the population is given by the time average of the logarithm of the regional growth rate: r j = ln λ j t Lewontin and Cohen, If we assume that fluctuations in fecundity and population are small Oσ relative to their means, then Covλ j,ν j x is Oσ 2 and we can Taylor expand the logarithm to Oσ 2 to obtain r j ln λ 0 j + 1 λ 0 j λ j x λ 0 j + Covλ j,ν j x 1 2λ j λ j x λ 0, 3 λ 0 j is the growth rate of species j in the absence of spatial or temporal variation. Happily, we need only deal with a single term of this expression. Our ESS condition, eq. 1, is a function of the invader long-run growth rate, r i. The invader is presumed to be of such low density that it does not contribute to competition, and thus λ i does not depend on the invader population distribution or dispersal distance. This assumption can be violated if the invader rapidly forms clusters in areas far from other sources of competition. Here, however, the invader will cluster in the same locations as the same-species resident, whose contribution to competition will dwarf that of the invader s. The only term that depends 1 on invader dispersal is Covλ i,ν i x t. Thus, the ESS dispersal condition reduces to λ 0 i Covλ i,ν i x t d i,d a,d p d i = 0. 4 di =d a=d a We now need an approximate expression for Covλ i,ν i x t, valid when fluctuations in fecundity and population density are small. Let us write F j x,t = F j x t1 + ǫ j x,t, 5 ǫ j x,t is Oσ and has spatial mean zero. Furthermore, assume that there is no spatially synchronized component to the variation in fecundity: F j x does not vary with time and is equal to F j x,t, the spatiotemporal average. Thus, By similar reasoning, let F j x,t = F j x,t 1 + ǫ j x,t. 6 n j x,t = n j x t1 + u j x,t = n j x,t 1 + u j x,t, 7 u j x,t is Oσ and has zero spatial and temporal mean. We want an Oσ 2 expression for Covλ i,ν i recall that our original expansion of r i was Oσ 2, which requires us to find Oσ expressions for ν i x,t and λ i x,t. The technical definition of Oσ n is that if gx is Oσ n, then gx decreases with σ and gx σ n can be made less than or equal to some positive constant K for σ small enough. On a more practical note, σ is some measure of smallness that we use for bookkeeping purposes. 2 j t
3 The relative population density ν i x,t is simply 1+u i x,t, which is already Oσ. What about λ i x,t? Assuming an annual invader, the local invader growth rate is given by λ i x,t = g i F i x,t γ i1 U i1 g 1 n 1 x,t + γ i2 U i2 n 2 x,t + s i1 g i, 8 the notation f g denotes a convolution: fx ygydy. Substituting our expressions for F j and n j into the expression for λ i and Taylor expanding to Oσ, we find that λ i x,t 1 + ǫ i x,t D i1 U i1 u 1 x,t D i2 U i2 u 2 x,t 9 Thus, D i1 = γ i1 g 1 n 1 x,t, D i2 = γ i2 n 2 x,t. 10 Covλ i,ν i x t = Cov ǫ i x,t D i1 U i1 u 1 D i2 U i2 u 2,u i. x,t 11 We can simplify both the covariance and the convolutions by taking the spatiotemporal Fourier transform of this expression. The Wiener-Khinchin theorem states that the Fourier transform of Covf,g x,t is given by FCovf,g x,t = lim N f N q,ω g N q,ω N 2, 12 superscript * denotes the complex conjugate and g N q,ω equals the Fourier transform of gx,t in the limit as N approaches infinity: g N q,ω = N/2 x,t= N/2 gx,t exp iqx+ ωt. Furthermore, the Fourier transform of a convolution is the product of the transforms of the factors: F f g = f g. Thus, the Fourier transform of eq. 11 is given by F Covλ i,ν i x t = F Covλ i,ν i x,t = lim ǫ N i 1 N 2 N q,ω D i1 Ũ i1 qũ N 1 q,ω D i2 Ũ i2 qũ N 2 q, ω ũ N i q,ω. 13 Honestly, this is an improvement. When we come to evaluate our ESS condition, eq. 4, we will find Covλ i,ν i x t by taking the inverse Fourier transform of eq. 13, numerically integrating 1 2π 2 F Covλ i,ν i x,t dq dω. 14 There would normally be a factor of expiqx+ωt in the integrand, x and t are the spatial and temporal lags of the covariance, but we wish to find the covariance at zero lag, so the exponential equals 1. 3
4 We have now found an Oσ 2 condition for the ESS eq. 4, which depends on Covλ i,ν i, and we have found an Oσ 2 expression for this covariance eqs. 13 and 14. We aren t finished, because we ultimately need to express our ESS condition entirely in terms of environmental conditions, as Covλ i,ν i depends on the Fourier transforms of the population densities: ũ i q,ω, ũ 1 q,ω, and ũ 2 q,ω. Our next step is to express the population densities as functions of the environment. We will begin with ũ 1 q,ω. We can obtain an equation for the annual relative population density, ν 1 x,t, by dividing both sides of the dynamical equation for n 1 eq. 3 by n 1 x t+1 and substituting λ 1 t n 1 x t for n 1 x t + 1 on the righthand side: ν 1 x,t + 1 = 1 λ 1 t k 1 x y g 1F 1 y,tν 1 y,t C 1 y,t dy + 1 λ 1 t s 11 g 1 ν 1 x,t. 15 To Oσ, however, λ 1 t can be replaced by λ 0 1, which equals 1. The annual resident is at equilibrium in the absence of environmental variation. Substituting eqs. 6 and 7 for F and n, we Taylor expand to Oσ to find g 1 F 1 x,t u 1 x,t + 1 = k 1 ǫ 1 + u 1 D 11 U 11 u 1 D 12 U 12 u 2 x,t + s 1 1 g 1 u 1 x,t, 16 D 11 = γ 11 g 1 n 1 x,t, D 12 = Taking a spatial Fourier transform, we find γ 12 n 2 x,t. 17 g 1 F 1 x,t ũ 1 q,t + 1 = k 1 q ǫ 1 q,t + ũ 1 q,t D 11 Ũ 11 qũ 1 q,t D 12 Ũ 12 qũ 2 q,t + s 1 1 g 1 ũ 1 q,t, and A u11 q = c 1 k1 q = A u11 qũ 1 q,t + A u12 qũ 2 q,t + A ǫ11 q ǫ 1 q,t + A ǫ12 q ǫ 2 q,t, 18 1 D 11 Ũ 11 q + s 1 1 g 1 A u12 q = c 1 k1 qd 12 Ũ 12 q 19 A ǫ11 q = c 1 k1 q A ǫ12 q = 0 20 c 1 = By the same process, we find g 1 F 1 x,t = 1 s 1 1 g ũ 2 q,t + 1 = A u21 qũ 1 q,t + A u22 qũ 2 q,t + A ǫ21 q ǫ 1 q,t + A ǫ22 q ǫ 2 q,t, 22 4
5 A u21 q = 1 s 2 k 2 qd 21 Ũ 21 q A u22 q = 1 s 2 k 2 q 1 D 22 Ũ 22 q + s 2 23 A ǫ21 q = 0 A ǫ22 q = 1 s 2 k 2 q. 24 Creating a vector of population fluctuations ũq,t = ũ 1,ũ 2 T q,t and a vector of environmental fluctuations ǫq,t = ǫ 1, ǫ 2 T q,t, we can express the dynamics in compact form: ũq,t + 1 = A u qũq,t + A ǫ q ǫq,t. 25 We can write ũq,t in terms of a series of past environmental conditions: t 1 ũq,t = A t 1 j u qa ǫ q ǫq,j = j=0 Mq,n = Mq,t j ǫq,j, 26 j=0 { A n 1 u qa ǫ q n > 0 0 n The second expression for ũq,t is a time convolution of Mq,t j with ǫq,t. By taking the temporal Fourier transform, we can turn this convolution into a product: ũq, ω = Mq, ω ǫq, ω, ũq, ω and ǫq, ω are the spatiotemporal Fourier transforms of ux, t and ǫx,t and M is the temporal Fourier transform of Mq,t: Mq,ω = Mq,te iωt = t= t=1 A t 1 u qa ǫ qe iωt. 28 For any matrix A whose eigenvalues have modulus less than 1, j=0 Aj = I A 1, I is the identity matrix. We only consider fluctuations about stable population equilibria, so the eigenvalues of A u q all have moduli less than 1, and the series converges. Performing the sum and multiplying through by expiω, we get ũq,ω = e iω I A u q 1 Aǫ q ǫq,ω. 29 The transfer function e iω I A u q 1 A ǫ q is a 2 2 matrix. We can write the ijth component in polar form as R ij q,ω expiφ ij q,ω, so that, for example, ũ 1 q,ω = R 11 q,ωe iφ 11q,ω ǫ 1 q,ω + R 12 q,ωe iφ 12q,ω ǫ 2 q,ω, 30 R ij q,ω equals the squareroot of the square of the real part of e iω I A u q 1 A ǫ q ij plus the square of the imaginary part, while φ ij q,ω equals the inverse tangent of the imaginary part divided by the real part. We have our ESS condition eq. 4. We have expressed the ESS condition in terms of the resident and invader population distributions eq. 13, with eq. 14. We have expressions for the resident population distributions eq. 29. All that remains is to find an expression 5
6 for the invader population distribution. We follow the same procedure as for the residents, writing ν i x,t + 1 = 1 λ i t k i x y g if i y,tν i y,t C i y,t dy + 1 λ i t s i1 g i ν i x,t. 31 We again substitute the growth rate in the absence of variation, λ 0 i, for λ i ; however because the invader is not at equilibrium, λ 0 i does not equal 1. Except for this change, the analysis remains the same. We substitute eqs. 6 and 7 for F and n, Taylor expand to Oσ, and take the spatial Fourier transform to find ũ i q,t + 1 = b ui qũ i q,t + b u1 qũ 1 q,t + b u2 qũ 2 q,t + b ǫi ǫ i q,t, 32 b ui q =c 1 ki q c 2 b ǫi q =c 1 ki q 33 b u1 q = D i1 Ũ i1 qb ǫi q b u2 q = D i2 Ũ i2 qb ǫi q, 34 and c 1 = 1 λ 0 i λ 0 i = c 2 = 1 λ 0 i s i 1 g i 35 + s i 1 g i, 36 D i1 = γ i1 g 1 n 1 x,t D i2 = γ i2 n 2 x,t. 37 In parallel with the development for the resident dynamics, we can now write t 1 ũ i q,t = b t 1 j u i qb ǫi q ǫ i q,t D i1 Ũ i1 qũ 1 q,t D i2 Ũ i2 qũ 2 q,t. 38 j=0 Again in parallel with the resident dynamics, we rewrite each term as a convolution, take the temporal Fourier transform, and find ũ i q,ω = Gq,ωe ǫ iψq,ω i q,ω D i1 Ũ i1 qũ 1 q,ω D i2 Ũ i2 qũ 2 q,ω, 39 b ǫi q Gq,ω = b 2 u i q 2b ui q cosω ψq,ω = tan 1 sinω, 41 cosω b ui q 6
7 and tan 1 has a range of π to π. Using eq. 30 and the equivalent for species 2 to substitute for ũ 1 q,ω and ũ 2 q,ω, we can express the invader distribution solely in terms of the environments experienced by the two residents and the invader: ũ i q,ω = Gq,ωe ǫ iψq,ω i q,ω D i1 Ũ i1 qr 11 q,ωe iφ11q,ω + D i2 Ũ i2 qr 21 q,ωe iφ 21q,ω ǫ 1 q,ω D i1 Ũ i1 qr 12 q,ωe iφ12q,ω + D i2 Ũ i2 qr 22 q,ωe iφ 22q,ω ǫ 2 q,ω. 42 We now have Oσ expressions for the Fourier transforms of the resident and invader populations in terms of their environments. Substituting these into our expression for the Fourier transform of Covλ i,ν i eq. 13, we can at last express Covλ i,ν i x,t entirely in terms of the environments experienced by the residents and the invader: 1 g i F i x,t F Covλ i,ν i x,t = lim ǫ N N N 2 i q, ω D i1 Ũ i1 qr 11 q,ωe iφ11q,ω + D i2 Ũ i2 qr 21 q,ωe iφ 21q,ω ǫ N 1 q, ω D i1 Ũ i1 qr 12 q,ωe iφ12q,ω + D i2 Ũ i2 qr 22 q,ωe iφ 22q,ω ǫ N 2 q, ω Gq,ωe iψq,ω ǫ N i q,ω D i1 Ũ i1 qr 11 q,ωe iφ11q,ω + D i2 Ũ i2 qr 21 q,ωe iφ 21q,ω D i1 Ũ i1 qr 12 q,ωe iφ12q,ω + D i2 Ũ i2 qr 22 q,ωe iφ 22q,ω ǫ N 1 q, ω ǫ N 2 q, ω As it stands, we would need to know precisely how the environment varied in space and time in order to calculate the ǫ N s and evaluate this expression. With a bit of cleverness, however, we can write the above in terms of the autocovariance of the environment: we don t care precisely how the environment varies, just how it is correlated in space and time. Eq. 43 contains many terms of the form stuff lim N ǫ N j k /N 2. We obtain lim N ǫ N j ǫ N k /N 2 from the Fourier transform of the covariance of ǫ j and ǫ k, which is proportional to the covariance of fecundities F j and F k : Covǫ j,ǫ k x,t x,t = CovF j,f k x,t x,t / F j x,t F k x,t = cosθv exp x /ξ exp t /τ from eq. 1. Thus, N ǫ N lim j ǫ N k References /N 2 = FCovǫ j,ǫ k x,t = cosθv = cosθv 1 e 2/τ 1 2e 1/τ cosω + e 2/τ x= t= e x /ξ e t /τ e iqx+ωt 1 e 2/ξ 1 2e 1/ξ cosq + e 2/ξ ǫ N Chesson, P General theory of competitive coexistence in spatially-varying environments. Theoretical Population Biology 58:
8 Geritz, S., É. Kisdi, G. Meszéna, and J. Metz Evolutionarily singular strategies and the adaptive growth and branching of the evolutionary tree. Evolutionary Ecology 12: Lewontin, R., and D. Cohen On population growth in a randomly varying environment. Zoology 62: Snyder, R. E Multiple risk reduction mechanisms: Can dormancy substitute for dispersal? Ecology Letters 9:
Chapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
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 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 informationKalle Parvinen. Department of Mathematics FIN University of Turku, Finland
Adaptive dynamics: on the origin of species by sympatric speciation, and species extinction by evolutionary suicide. With an application to the evolution of public goods cooperation. Department of Mathematics
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationHow the spatial scales of dispersal, competition, and environmental heterogeneity interact to affect coexistence
How the spatial scales of dispersal, competition, and environmental heterogeneity interact to affect coexistence Robin E. Snyder Department of Biology Case Western Reserve University 10900 Euclid St. Cleveland,
More information2.3 Damping, phases and all that
2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be
More information2.2 THE DERIVATIVE 2.3 COMPUTATION OF DERIVATIVES: THE POWER RULE 2.4 THE PRODUCT AND QUOTIENT RULES 2.6 DERIVATIVES OF TRIGONOMETRIC FUNCTIONS
Differentiation CHAPTER 2 2.1 TANGENT LINES AND VELOCITY 2.2 THE DERIVATIVE 2.3 COMPUTATION OF DERIVATIVES: THE POWER RULE 2.4 THE PRODUCT AND QUOTIENT RULES 25 2.5 THE CHAIN RULE 2.6 DERIVATIVES OF TRIGONOMETRIC
More informationModule BIO- M1- S06 Evolu-onary ecology. Adaptive dynamics. David Claessen Ins-tut de Biologie de l ENS Equipe Eco- Evolu-on Mathéma-que
Module BIO- M1- S06 Evolu-onary ecology Adaptive dynamics David Claessen Ins-tut de Biologie de l ENS Equipe Eco- Evolu-on Mathéma-que OBJECTIF Une théorie quantitative pour prédire la tendence evolutive
More informationD.S.G. POLLOCK: BRIEF NOTES
BIVARIATE SPECTRAL ANALYSIS Let x(t) and y(t) be two stationary stochastic processes with E{x(t)} = E{y(t)} =. These processes have the following spectral representations: (1) x(t) = y(t) = {cos(ωt)da
More informationPhysics tricks for fun and profit: A physicist s adventures in theoretical ecology p.1/44
Physics tricks for fun and profit: A physicist s adventures in theoretical ecology Robin E. Snyder robin.snyder@cwru.edu Department of Biology, Case Western Reserve University Physics tricks for fun and
More informationAssignment 8. [η j, η k ] = J jk
Assignment 8 Goldstein 9.8 Prove directly that the transformation is canonical and find a generating function. Q 1 = q 1, P 1 = p 1 p Q = p, P = q 1 q We can establish that the transformation is canonical
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More informationPreconditioners for ill conditioned (block) Toeplitz systems: facts a
Preconditioners for ill conditioned (block) Toeplitz systems: facts and ideas Department of Informatics, Athens University of Economics and Business, Athens, Greece. Email:pvassal@aueb.gr, pvassal@uoi.gr
More informationDamped Harmonic Oscillator
Damped Harmonic Oscillator Wednesday, 23 October 213 A simple harmonic oscillator subject to linear damping may oscillate with exponential decay, or it may decay biexponentially without oscillating, or
More informationChapter 11. Taylor Series. Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27
Chapter 11 Taylor Series Josef Leydold Mathematical Methods WS 2018/19 11 Taylor Series 1 / 27 First-Order Approximation We want to approximate function f by some simple function. Best possible approximation
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare http://ocw.mit.edu 5.74 Introductory Quantum Mechanics II Spring 9 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. Andrei Tokmakoff,
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationLecture 5 - Logarithms, Slope of a Function, Derivatives
Lecture 5 - Logarithms, Slope of a Function, Derivatives 5. Logarithms Note the graph of e x This graph passes the horizontal line test, so f(x) = e x is one-to-one and therefore has an inverse function.
More informationBackground ODEs (2A) Young Won Lim 3/7/15
Background ODEs (2A) Copyright (c) 2014-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any
More informationPhase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
More informationAdaptive dynamics in an individual-based, multi-resource chemostat model
Adaptive dynamics in an individual-based, multi-resource chemostat model Nicolas Champagnat (INRIA Nancy) Pierre-Emmanuel Jabin (Univ. Maryland) Sylvie Méléard (Ecole Polytechnique) 6th European Congress
More informationAn extension of the classification of evolutionarily singular strategies in Adaptive Dynamics
An extension of the classification of evolutionarily singular strategies in Adaptive Dynamics Barbara Boldin Department of Mathematics, Natural Sciences and Information Technologies University of Primorska
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationThe formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) R(x)
More informationSuperposition of electromagnetic waves
Superposition of electromagnetic waves February 9, So far we have looked at properties of monochromatic plane waves. A more complete picture is found by looking at superpositions of many frequencies. Many
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More informationChapter 9. Barotropic Instability. 9.1 Linearized governing equations
Chapter 9 Barotropic Instability The ossby wave is the building block of low ossby number geophysical fluid dynamics. In this chapter we learn how ossby waves can interact with each other to produce a
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationTime Series Analysis. Solutions to problems in Chapter 7 IMM
Time Series Analysis Solutions to problems in Chapter 7 I Solution 7.1 Question 1. As we get by subtraction kx t = 1 + B + B +...B k 1 )ǫ t θb) = 1 + B +... + B k 1 ), and BθB) = B + B +... + B k ), 1
More informationClassical Mechanics Phys105A, Winter 2007
Classical Mechanics Phys5A, Winter 7 Wim van Dam Room 59, Harold Frank Hall vandam@cs.ucsb.edu http://www.cs.ucsb.edu/~vandam/ Phys5A, Winter 7, Wim van Dam, UCSB Midterm New homework has been announced
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 informationReview of Linear Time-Invariant Network Analysis
D1 APPENDIX D Review of Linear Time-Invariant Network Analysis Consider a network with input x(t) and output y(t) as shown in Figure D-1. If an input x 1 (t) produces an output y 1 (t), and an input x
More informationSymmetric competition as a general. model for single-species adaptive. dynamics
Symmetric competition as a general model for single-species adaptive arxiv:1204.0831v1 [q-bio.pe] 3 Apr 2012 dynamics Michael Doebeli 1 & Iaroslav Ispolatov 1,2 1 Department of Zoology and Department of
More information1 Introduction. 2 Measure theoretic definitions
1 Introduction These notes aim to recall some basic definitions needed for dealing with random variables. Sections to 5 follow mostly the presentation given in chapter two of [1]. Measure theoretic definitions
More information221B Lecture Notes on Resonances in Classical Mechanics
1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small
More informationPhysics 11b Lecture #15
Physics 11b ecture #15 and ircuits A ircuits S&J hapter 3 & 33 Administravia Midterm # is Thursday If you can t take midterm, you MUST let us (me, arol and Shaun) know in writing before Wednesday noon
More information( ) f (k) = FT (R(x)) = R(k)
Solving ODEs using Fourier Transforms The formulas for derivatives are particularly useful because they reduce ODEs to algebraic expressions. Consider the following ODE d 2 dx + p d 2 dx + q f (x) = R(x)
More informationTaylor Series and stationary points
Chapter 5 Taylor Series and stationary points 5.1 Taylor Series The surface z = f(x, y) and its derivatives can give a series approximation for f(x, y) about some point (x 0, y 0 ) as illustrated in Figure
More informationLong-term adaptive diversity in Levene-type models
Evolutionary Ecology Research, 2001, 3: 721 727 Long-term adaptive diversity in Levene-type models Éva Kisdi Department of Mathematics, University of Turku, FIN-20014 Turku, Finland and Department of Genetics,
More informationz = x + iy ; x, y R rectangular or Cartesian form z = re iθ ; r, θ R polar form. (1)
11 Complex numbers Read: Boas Ch. Represent an arb. complex number z C in one of two ways: z = x + iy ; x, y R rectangular or Cartesian form z = re iθ ; r, θ R polar form. (1) Here i is 1, engineers call
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationab = c a If the coefficients a,b and c are real then either α and β are real or α and β are complex conjugates
Further Pure Summary Notes. Roots of Quadratic Equations For a quadratic equation ax + bx + c = 0 with roots α and β Sum of the roots Product of roots a + b = b a ab = c a If the coefficients a,b and c
More information2A1H Time-Frequency Analysis II Bugs/queries to HT 2011 For hints and answers visit dwm/courses/2tf
Time-Frequency Analysis II (HT 20) 2AH 2AH Time-Frequency Analysis II Bugs/queries to david.murray@eng.ox.ac.uk HT 20 For hints and answers visit www.robots.ox.ac.uk/ dwm/courses/2tf David Murray. A periodic
More informationSyllabus: for Complex variables
EE-2020, Spring 2009 p. 1/42 Syllabus: for omplex variables 1. Midterm, (4/27). 2. Introduction to Numerical PDE (4/30): [Ref.num]. 3. omplex variables: [Textbook]h.13-h.18. omplex numbers and functions,
More informationThe Feynman-Kac formula
The Feynman-Kac formula William G. Faris February, 24 The Wiener process (Brownian motion) Consider the Hilbert space L 2 (R d ) and the self-adjoint operator H = σ2, () 2 where is the Laplace operator.
More informationProblem Set 1 October 30, 2017
1. e π can be calculated from e x = x n. But that s not a good approximation method. The n! reason is that π is not small compared to 1. If you want O(0.1) accuracy, then the last term you need to include
More information5. THE CLASSES OF FOURIER TRANSFORMS
5. THE CLASSES OF FOURIER TRANSFORMS There are four classes of Fourier transform, which are represented in the following table. So far, we have concentrated on the discrete Fourier transform. Table 1.
More informationCONTINUOUS IMAGE MATHEMATICAL CHARACTERIZATION
c01.fm Page 3 Friday, December 8, 2006 10:08 AM 1 CONTINUOUS IMAGE MATHEMATICAL CHARACTERIZATION In the design and analysis of image processing systems, it is convenient and often necessary mathematically
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of under-damping,
More informationCausality. but that does not mean it is local in time, for = 1. Let us write ɛ(ω) = ɛ 0 [1 + χ e (ω)] in terms of the electric susceptibility.
We have seen that the issue of how ɛ, µ n depend on ω raises questions about causality: Can signals travel faster than c, or even backwards in time? It is very often useful to assume that polarization
More informationEvolutionary Dynamics of Seed Size and Seedling Competitive Ability
Theoretical Population Biology 55, 324343 (1999) Article ID tpbi.1998.1409, available online at http:www.idealibrary.com on Evolutionary Dynamics of Seed Size and Seedling Competitive Ability Stefan A.
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationIntroduction to Dynamical Systems
Introduction to Dynamical Systems Autonomous Planar Systems Vector form of a Dynamical System Trajectories Trajectories Don t Cross Equilibria Population Biology Rabbit-Fox System Trout System Trout System
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationVer 3808 E1.10 Fourier Series and Transforms (2014) E1.10 Fourier Series and Transforms. Problem Sheet 1 (Lecture 1)
Ver 88 E. Fourier Series and Transforms 4 Key: [A] easy... [E]hard Questions from RBH textbook: 4., 4.8. E. Fourier Series and Transforms Problem Sheet Lecture. [B] Using the geometric progression formula,
More information13.42 READING 6: SPECTRUM OF A RANDOM PROCESS 1. STATIONARY AND ERGODIC RANDOM PROCESSES
13.42 READING 6: SPECTRUM OF A RANDOM PROCESS SPRING 24 c A. H. TECHET & M.S. TRIANTAFYLLOU 1. STATIONARY AND ERGODIC RANDOM PROCESSES Given the random process y(ζ, t) we assume that the expected value
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More informationIntroduction to Aspects of Multiscale Modeling as Applied to Porous Media
Introduction to Aspects of Multiscale Modeling as Applied to Porous Media Part III Todd Arbogast Department of Mathematics and Center for Subsurface Modeling, Institute for Computational Engineering and
More information3 + 4i 2 + 3i. 3 4i Fig 1b
The introduction of complex numbers in the 16th century was a natural step in a sequence of extensions of the positive integers, starting with the introduction of negative numbers (to solve equations of
More informationChapter 4. Section Derivatives of Exponential and Logarithmic Functions
Chapter 4 Section 4.2 - Derivatives of Exponential and Logarithmic Functions Objectives: The student will be able to calculate the derivative of e x and of lnx. The student will be able to compute the
More information1. There are 8 questions spanning 9 pages total (including this cover page). Please make sure that you have all 9 pages before starting.
Instructor: K. Rotz Name: Solution PUID: 00000-00000 Instructions and tips: 1. There are 8 questions spanning 9 pages total (including this cover page). Please make sure that you have all 9 pages before
More information6.4 Logarithmic Equations and Inequalities
6.4 Logarithmic Equations and Inequalities 459 6.4 Logarithmic Equations and Inequalities In Section 6.3 we solved equations and inequalities involving exponential functions using one of two basic strategies.
More informationPHYS 3900 Homework Set #03
PHYS 3900 Homework Set #03 Part = HWP 3.0 3.04. Due: Mon. Feb. 2, 208, 4:00pm Part 2 = HWP 3.05, 3.06. Due: Mon. Feb. 9, 208, 4:00pm All textbook problems assigned, unless otherwise stated, are from the
More informationDIFFERENCE AND DIFFERENTIAL EQUATIONS COMPARED
DIFFERENCE AND DIFFERENTIAL EQUATIONS COMPARED Frequency Response Function: Continuous Systems Given the unit impulse response h(t) of a continuous-time system, the mapping from the system s input x(t)
More informationHow many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation?
How many initial conditions are required to fully determine the general solution to a 2nd order linear differential equation? (A) 0 (B) 1 (C) 2 (D) more than 2 (E) it depends or don t know How many of
More informationAnalytical Solution for Acoustic/Poroelastic Interaction
Analytical Solution for Acoustic/Poroelastic Interaction Julien Diaz,2 Abdelâaziz Ezziani 2, INRIA Bordeaux Sud-Ouest, EPI Magique 3D 2 LMA, UMR 542, Université de Pau 7 April 2009 Acoustic/Poroelastic
More informationAdditional material: Linear Differential Equations
Chapter 5 Additional material: Linear Differential Equations 5.1 Introduction The material in this chapter is not formally part of the LTCC course. It is included for completeness as it contains proofs
More informationThe Calculus of Residues
hapter 7 The alculus of Residues If fz) has a pole of order m at z = z, it can be written as Eq. 6.7), or fz) = φz) = a z z ) + a z z ) +... + a m z z ) m, 7.) where φz) is analytic in the neighborhood
More information(x 3)(x + 5) = (x 3)(x 1) = x + 5. sin 2 x e ax bx 1 = 1 2. lim
SMT Calculus Test Solutions February, x + x 5 Compute x x x + Answer: Solution: Note that x + x 5 x x + x )x + 5) = x )x ) = x + 5 x x + 5 Then x x = + 5 = Compute all real values of b such that, for fx)
More informationExplicit expression for Lindhard dielectric function at finite temperature
arxiv:1412.5705v1 [cond-mat.mtrl-sci] 18 Dec 2014 Explicit expression for Lindhard dielectric function at finite temperature A. V. Andrade-Neto Departamento de Física, Universidade Estadual de Feira de
More informationWave Phenomena Physics 15c. Lecture 11 Dispersion
Wave Phenomena Physics 15c Lecture 11 Dispersion What We Did Last Time Defined Fourier transform f (t) = F(ω)e iωt dω F(ω) = 1 2π f(t) and F(w) represent a function in time and frequency domains Analyzed
More informationIV. Covariance Analysis
IV. Covariance Analysis Autocovariance Remember that when a stochastic process has time values that are interdependent, then we can characterize that interdependency by computing the autocovariance function.
More information2.6 Logarithmic Functions. Inverse Functions. Question: What is the relationship between f(x) = x 2 and g(x) = x?
Inverse Functions Question: What is the relationship between f(x) = x 3 and g(x) = 3 x? Question: What is the relationship between f(x) = x 2 and g(x) = x? Definition (One-to-One Function) A function f
More informationLinear second-order differential equations with constant coefficients and nonzero right-hand side
Linear second-order differential equations with constant coefficients and nonzero right-hand side We return to the damped, driven simple harmonic oscillator d 2 y dy + 2b dt2 dt + ω2 0y = F sin ωt We note
More informationEvolutionarily singular strategies and the adaptive growth and branching of the evolutionary tree
Evolutionarily singular strategies and the adaptive growth and branching of the evolutionary tree S.A.H. GERITZ 1,, É. KISDI 1,3, G. MESZÉNA 4, J.A.J. METZ,5 1 Collegium Budapest, Institute for Advanced
More informationt For l = 1 a monomer cannot be destroyed or created from nothing: = b p(2,t) a p(1,t).
IITS: Statistical Physics in Biology Assignment # 5 KU Leuven 5/31/2013 Drift, Diffusion, and Dynamic Instability 1. Treadmilling Actin: Actin filaments are long, asymmetric, polymers involved in a variety
More informationLINEAR RESPONSE THEORY
MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior
More informationThe Zel dovich Approximation. Large-scale structure goes ballistic
The Zel dovich Approximation Large-scale structure goes ballistic Martin White UCB/LBNL Built on work done with Lile Wang (undergrad), Jordan Carlson (grad) and Beth Reid (postdoc). Yakov Borisovich Zel
More informationThe Derivative of a Function Measuring Rates of Change of a function. Secant line. f(x) f(x 0 ) Average rate of change of with respect to over,
The Derivative of a Function Measuring Rates of Change of a function y f(x) f(x 0 ) P Q Secant line x 0 x x Average rate of change of with respect to over, " " " " - Slope of secant line through, and,
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationAdaptive evolution : A population point of view Benoît Perthame
Adaptive evolution : A population point of view Benoît Perthame 30 20 t 50 t 10 x 0 1 0.5 0 0.5 1 x 0 1 0.5 0 0.5 1 Population formalism n(x, t) t }{{} variation of individuals = birth with mutations {}}{
More informationarxiv:hep-th/ v3 16 May 1996
BNL-63106 An Exact Solution for Quantum Tunneling in a Dissipative System arxiv:hep-th/9605081v3 16 May 1996 Li Hua Yu National Synchrotron Light Source, Brookhaven National Laboratory, N.Y.11973 Abstract
More information2.161 Signal Processing: Continuous and Discrete Fall 2008
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationApplied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation
22.101 Applied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation References -- R. M. Eisberg, Fundamentals of Modern Physics (Wiley & Sons, New York, 1961). R. L. Liboff, Introductory
More informationMath 113 (Calculus 2) Exam 4
Math 3 (Calculus ) Exam 4 November 0 November, 009 Sections 0, 3 7 Name Student ID Section Instructor In some cases a series may be seen to converge or diverge for more than one reason. For such problems
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 informationConnections and geodesics in the space of metrics The exponential parametrization from a geometric perspective
Connections and geodesics in the space of metrics The exponential parametrization from a geometric perspective Andreas Nink Institute of Physics University of Mainz September 21, 2015 Based on: M. Demmel
More informationOn the Stability of the Best Reply Map for Noncooperative Differential Games
On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: May 24, 2011. Lecture 6p: Spectral Density, Passing Random Processes through LTI Systems, Filtering Terms
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 informationMath53: Ordinary Differential Equations Autumn 2004
Math53: Ordinary Differential Equations Autumn 2004 Unit 2 Summary Second- and Higher-Order Ordinary Differential Equations Extremely Important: Euler s formula Very Important: finding solutions to linear
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationPhasor Young Won Lim 05/19/2015
Phasor Copyright (c) 2009-2015 Young W. Lim. Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version
More informationMATH 135: COMPLEX NUMBERS
MATH 135: COMPLEX NUMBERS (WINTER, 010) The complex numbers C are important in just about every branch of mathematics. These notes 1 present some basic facts about them. 1. The Complex Plane A complex
More information5.74 Introductory Quantum Mechanics II
MIT OpenCourseWare ttp://ocw.mit.edu 5.74 Introductory Quantum Mecanics II Spring 9 For information about citing tese materials or our Terms of Use, visit: ttp://ocw.mit.edu/terms. Andrei Tokmakoff, MIT
More informationP441 Analytical Mechanics - I. RLC Circuits. c Alex R. Dzierba. In this note we discuss electrical oscillating circuits: undamped, damped and driven.
Lecture 10 Monday - September 19, 005 Written or last updated: September 19, 005 P441 Analytical Mechanics - I RLC Circuits c Alex R. Dzierba Introduction In this note we discuss electrical oscillating
More informationLimits and Continuity. 2 lim. x x x 3. lim x. lim. sinq. 5. Find the horizontal asymptote (s) of. Summer Packet AP Calculus BC Page 4
Limits and Continuity t+ 1. lim t - t + 4. lim x x x x + - 9-18 x-. lim x 0 4-x- x 4. sinq lim - q q 5. Find the horizontal asymptote (s) of 7x-18 f ( x) = x+ 8 Summer Packet AP Calculus BC Page 4 6. x
More information