x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).
|
|
- Claud Wilson
- 5 years ago
- Views:
Transcription
1
2 Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ). Oliver Junge, Institute for Mathematics, University of Paderborn 2
3 Continuation of equilibria Structure of the Solution Let u 0 = (x 0, λ 0 ) s.t. f(u 0 ) = 0 and rank(f (u 0 )) = N. Then locally f 1 (0) is a one-dimensional manifold in R d+1. x f 1 (0) u 0 λ Oliver Junge, Institute for Mathematics, University of Paderborn 3
4 Continuation of equilibria Tangent vector Let c : ( 1, 1) R d+1 be a local parametrization of f 1 (0). Then f(c(s)) = 0, i.e. In other words f (c(s)) c (s) = 0. T c(s) f 1 (0) = ker f (c(s)). f 1 (0) u 0 ker f (u 0 ) Oliver Junge, Institute for Mathematics, University of Paderborn 4
5 Continuation of equilibria Fixing the orientation We will use ker f (c(s)) to move along c. In order to fix a consistent direction, we require that c (s) = 1 and det f (c(s)) c (s) T > 0. Exercise 1 Show that f (c(s)) c (s) T is regular. Oliver Junge, Institute for Mathematics, University of Paderborn 5
6 Continuation of equilibria The induced tangent vector For a d d + 1-matrix A with rank(a) = d let t = t(a) R d+1 be the unique vector such that (i) At = 0; (ii) t = 1; (iii) det A t T > 0. We call t(a) the tangent vector induced by A. Exercise 2 Show that the set of all d (d + 1)-matrices with full rank is open. Oliver Junge, Institute for Mathematics, University of Paderborn 6
7 Continuation of equilibria The associated differential equation Consider the ODE c = t(f (c)). (1) Since d f(c(s)) = f (c(s)) c (s) = 0, f is constant along solutions of ds (1). Idea: in order to compute the path c, solve (1). Given data: vector field t(f ( )) and condition f(c(s)) = 0. predictor-corrector methods. Oliver Junge, Institute for Mathematics, University of Paderborn 7
8 Continuation of equilibria Predictor-corrector scheme Given u 0 with f(u 0 ) < ε, find sequence u 1, u 2,... of points such that f(u i ) < ε for some prescribed accuracy ε > 0. f 1 (0) u 1 u 2 u 0 (i) Predictor: solve the initial value problem u = t(f (u)), u(0) = u i, by an explicit scheme, e.g. one Euler-step: where h is the stepsize. ũ i+1 = u i + ht(f (u i )), Oliver Junge, Institute for Mathematics, University of Paderborn 8
9 Continuation of equilibria (ii) Corrector: compute u i+1 = argmin f(u)=0 u ũ i+1 by some iteration scheme, e.g. Newton s method. u i+1 t(f (u i+1 )) u i ũ i+1 A necessary condition for u i+1 is f(u i+1 ) = 0, t(f (u i+1 ))(u i+1 ũ i+1 ) = 0. Oliver Junge, Institute for Mathematics, University of Paderborn 9
10 Computation of bifurcation points Saddle-node bifurcation x x λ λ Example 1 ẋ = λ + x 2 λ 0 λ < 0 λ 0 λ = 0 0 λ > 0 Oliver Junge, Institute for Mathematics, University of Paderborn 10
11 Computation of bifurcation points Saddle-node bifurcation Condition: dim ker f x (x, λ ) = 1 (but still rank(f (x, λ )) = d). Detection: during path-following, check sign of det f x (x(s), λ(s)). Computation: Solve (Moore, Spence, 1980) f(x, λ) = 0 f x (x, λ)φ = 0 l T φ 1 = 0, l R d. Oliver Junge, Institute for Mathematics, University of Paderborn 11
12 Computation of bifurcation points Saddle-node bifurcation Let ψ f x (x, λ ) = 0. Theorem 1 If dim ker f (x, λ ) = d and dim ker f x (x, λ ) = 1, then (x, φ, λ ) is a regular zero of the extended system, if and only if ψ f xx (x, λ )φ φ 0. The point (x, λ ) is a quadratic turning point. Oliver Junge, Institute for Mathematics, University of Paderborn 12
13 Computation of bifurcation points Hopf bifurcation Consider the system ẋ = λx y x 3 ẏ = x for various λ R. Oliver Junge, Institute for Mathematics, University of Paderborn 13
14 Computation of bifurcation points Eigenvalue movement Linearization in 0, eigenvalues: Df(0) = λ 1 1 0, µ 1,2 = λ 2 ± λ µ Im µ Re Oliver Junge, Institute for Mathematics, University of Paderborn 14
15 Computation of bifurcation points The Theorem of Hopf (x 0, λ 0 ) equilibrium, iω 0 simple eigenvalue of f x (x 0, λ 0 ), ir σ(f x (x 0, λ 0 )) = {iω, iω}, C = {(x(λ), λ)} local path of equilibria, µ(λ) + iν(λ) eigenvalue of f x (x(λ), λ) with µ(λ 0 ) + iν(λ 0 ) = iω, If dµ dλ (λ 0) 0, then there exists a unique one-parameter family of periodic solutions (with positive period) in a neighborhood of (x 0, λ 0 ). Oliver Junge, Institute for Mathematics, University of Paderborn 15
16 Computation of bifurcation points Computing Hopf points Consider an eigenvector v + iw at iω of D := d x f(x 0, λ 0 ), i.e. D(v + iw) = iω(v + iw) Dv = ωw and Dw = ωv. Exercise 3 Show that v and w are linearly independent. Thus D 2 v = ω 2 v and D 2 w = ω 2 w. Oliver Junge, Institute for Mathematics, University of Paderborn 16
17 Computation of bifurcation points Computing Hopf points Consider the extended system F (x, ϕ, λ, ν) = f(x, λ) (d x f(x, λ) 2 + ν 2 I)ϕ (ϕ, ϕ) 1 l T ϕ = 0. If l is not orthogonal to span{v, w}, then there exists a unique (up to its sign) vector ϕ 0 span{v, w}, such that (x 0, ϕ 0, λ 0, ω) is a regular a solution of this system. a Roose,Hlavacek, 85 Oliver Junge, Institute for Mathematics, University of Paderborn 17
18 Goal: compute periodic solution of x = f(x), x R d, i.e. a solution x(t) with x(t) = x(t + T ) for some unknown period T = 2π/ω. The transformation y(t) = x(t/ω) yields the system y = 1 ω f(y), for which ȳ(t) = x(t/ω) is a 2π-periodic solution. Oliver Junge, Institute for Mathematics, University of Paderborn 18
19 Boundary value problem: y = 1 ω f(y) y(0) = y(2π). Let ϕ t (, ω) denote the flow of this ode. Shooting method: solve S(ξ, ω) = ϕ 2π (ξ, ω) ξ = 0. Because of the S 1 -symmetry of each periodic solution we need an additional phase condition. Simple choice: ξ 1 = ξ 1 Oliver Junge, Institute for Mathematics, University of Paderborn 19
20 The variational equation In order to solve S(ξ, ω) = 0, ξ 1 = ξ 1 by Newton s method, we need to compute DS. D ξ S(ξ, ω) = D ξ ϕ 2π (ξ, ω) I can be computed via the variational equation t D ξϕ t (ξ, ω) = 1 ω Df(ϕt (ξ, ω))d ξ ϕ t (ξ, ω), D ξ ϕ 0 (ξ, ω) = I. For S (ξ, ω) we compute ω ω ϕt = t ω 2 ẋ(t/ω) = t ω 2 f(x(t/ω)), Oliver Junge, Institute for Mathematics, University of Paderborn 20
21 thus S (ξ, ω) = ω ω ϕ2π (ξ, ω) = 2π ω 2 f(ϕ2π (ξ, ω)). Advantage of shooting: stability of the periodic solution is directly given by the eigenvalues of D ξ ϕ 2π (ξ, ω). Oliver Junge, Institute for Mathematics, University of Paderborn 21
22 Galerkin approach Let C r 2π = {u C r (R, R d ) : u(s + 2π) = u(s)}. Consider the operator F : C 1 2π R C 0 2π, F (y, ω) = dy dt 1 ω f(y). A zero of this operator is a periodic solution of the boundary value problem. Oliver Junge, Institute for Mathematics, University of Paderborn 22
23 Line of Reasoning (i) Choose a suitable, countable basis of the underlying space countable system of equations; (ii) using finitely many modes, numerically compute an approximate solution (Galerkin); (iii) construct a restricted domain for F that isolates the numerical zero; (iv) using topological arguments, show that there actually exists a zero in the restricted domain. Oliver Junge, Institute for Mathematics, University of Paderborn 23
24 Setting F (y, y,..., y (r) ) = 0, y(t) R, F : R r+1 R smooth. We look for functions y : R R, y C r, y(t) = y(t + 2π/ω) for some ω R such that F (y, y,..., y (r) ) = 0. Rescaling time, x(t) = y(t/ω), yields the map ˆF : R C r 1 C 0 1 ˆF : (ω, x) F (x, ωx,..., ω r x (r) ) Oliver Junge, Institute for Mathematics, University of Paderborn 24
25 ˆF induces a map on the Fourier coefficients (c k ) k Z of x. Note that c k = c k, k Z, so it suffices to consider { l 2 = (c k ) k=0 C N : c k 2 < }. k=0 Correspondingly l 2 r = { (c k ) k=0 l 2 u C r 1 : c k is the k-th Fourier coefficient of u }. Induced map F : R l 2 r l 2 0. Oliver Junge, Institute for Mathematics, University of Paderborn 25
26 Phase Condition Since the original ODE is autonomous: x C1 r solution x( + t) solution t [0, 1]. Numerically more favorable: regularize by phase condition ϕ(c) = 0. Full system: Φ : R l 2 r R l 2 0 Φ(ω, c) = (ϕ(c), F (ω, c)). Oliver Junge, Institute for Mathematics, University of Paderborn 26
27 We look for (ω, c) such that Φ(ω, c) = 0. Equivalently: look for fixed points of G = id + Φ. Define and consider the sets A k (ω, c) = c k P k (ω, c) = (ω, c 0,..., c k ) Q k (ω, c) = (c k+1, c k+2,...) Z k = {(ω, c) : P k G(ω, c) = P k (ω, c)} Z = k 0 Z k Oliver Junge, Institute for Mathematics, University of Paderborn 27
28 Proposition 1 If G is continuous, Z M is compact for some M and Z k for k M, then Z is nonempty and all points in Z are fixed points of G. Proof of Existence (i) Find compact restricted domain for G; (ii) Show that Z M is nonempty for some (small) M; (iii) Show that Z k is nonempty for k > M. Oliver Junge, Institute for Mathematics, University of Paderborn 28
29 (i) Construction of the Restricted Domain (i) Use Newton s method on a Galerkin projection of G in order to estimate a fixed point; (ii) F C r x C r c k = O(k r ) as k. Restricted domain: compact set Ω D R l 2 r, where D = D 0 D M 1 {c k C : c k r k }, k=m where D 0 D M 1 C M contains the numerically computed fixed point. Oliver Junge, Institute for Mathematics, University of Paderborn 29
30 (ii) Z M is nonempty Define the multivalued map G M : P M (Ω D) P M (R l 2 ) (ω, c) {P M G(ω, c) : P M (ω, c) = (ω, c) and (ω, c) Ω D}. We will show that every continuous selector of G M has a fixed point. Note that the set Q M (Ω D) determines the size of the images of G M. Oliver Junge, Institute for Mathematics, University of Paderborn 30
31 (iii) Z k is nonempty for k M Definition 1 A map f : Ω D C is linearly dominated on the ball B 0 (r) A k (Ω D) if f(ω, c) = L(ω)c k + g(ω, c), for all (ω, c) Ω D, c k r, where L : I C, Ω I, and g are continuous functions, such that sup Ω D g(ω, c) < r inf ω I L(ω) r. (2) Oliver Junge, Institute for Mathematics, University of Paderborn 31
32 r r R R f( Ω,D ) f( Ω,D ) (a) Linearly dominated. (b) Not linearly dominated. R = r inf ω I L(ω), D = D { c k = r}. Oliver Junge, Institute for Mathematics, University of Paderborn 32
33 Existence Theorem Theorem 2 Let Ω D be compact. Suppose that there exists a starshaped enclosure ĜM for G M such that (N, L) is an index pair for ĜM and Λ((N, L), ĜM) 0, where Λ is the Lefschetz number of the index pair (N, L). Assume furthermore that the maps A k G are linearly dominated for all k > M, then there exists a fixed point for the map G. Oliver Junge, Institute for Mathematics, University of Paderborn 33
34 Example Consider y + y + σy δy + y 2 = 0, σ = 2, δ = 3. Induced map on Fourier coefficients A k F : (ω, c) ( iω 3 k 3 ω 2 k 2 + σωik δ)c k + l Z c l c k l, k 0. Computations: real version. Phase condition ϕ(c) = imag(c 1 ) = 0. G(ω, c) = (ω, c) + (ϕ(c), F (ω, c)) as above. Oliver Junge, Institute for Mathematics, University of Paderborn 34
35 Numerical fixed point ω = 1.39, c 0 = 2.46, c 1 = 0.813, c 2 = i , c 3 = i , c 4 = i , c 5 = i , c 6 = i , c 7 = i , c 8 = i , c 9 = i Oliver Junge, Institute for Mathematics, University of Paderborn 35
36 Numerical fixed point: decay of coefficients c k β α k k Oliver Junge, Institute for Mathematics, University of Paderborn 36
37 Restricted domain Ω = [1.36, 1.44], D = D 0 D 1 D 2 B 0 (βα k ), with α = 0.1, β = 11.5, and D 0 = [2.43, 2.49], D 1 = [0.8, ], k=3 D 2 = [ 0.003, 0.04] [0.005, 0.06]i, imag(d 0 ) = 0, since the Fourier series is real valued. imag(d 1 ) = 0 due to the phase condition. Oliver Junge, Institute for Mathematics, University of Paderborn 37
38 Linear domiation for k > 2 Lemma: [ ] c l c k l β 2 α k 2 1 α 2 + k 1. l Z Lemma: The map A k G is linearly dominated on the ball B 0 (βα k ) for k > 2. Proof: By the estimate we need to verify that [ ] β 2 α k 2 1 α 2 + k 1 < βα k min ω Ω ( iω 3 k 3 ω 2 k 2 + σωik δ ). Rewrite the rhs, use Ω = [ω, ω], plug in ω and ω where appropriate, function k e(k, ω, ω, δ, σ) such that A k is linearly dominated, Oliver Junge, Institute for Mathematics, University of Paderborn 38
39 if e(k) > 0. Result: e(k) is increasing and e(3) > 0. Index pair for lower modes We need to numerically construct an index pair for Ĝ2 (with, hopefully, a non-zero Lefschetz number). 7d 3d by exploiting the structure of the map; Oliver Junge, Institute for Mathematics, University of Paderborn 39
40 Explicit multivalued map a 1 ( ω 2 δ + 1)a 1 + 2a 1 a 2 + 2ā 0 a 1 + c 2 I I 1 a 2 ( 4 ω 2 δ + 1)a 2 + (8 ω 3 2σ ω)b 2 + 2ā 0 a 2 + a ( c c 1 )I I 1 b 2 ( 8 ω 3 + 2σ ω)a 2 + ( 4 ω 2 δ + 1)b 2 + 2ā 0 b 2 + c I 0 + c 1 I I 1, where I 0 = [ 1, 1], I 1 = [ 1, 1], and ā 0 = ā 0 (a 1, a 2, b 2 ) and ω = ω(a 1, a 2, b 2 ) are intervals. Oliver Junge, Institute for Mathematics, University of Paderborn 40
41 Isolating neighborhood Oliver Junge, Institute for Mathematics, University of Paderborn 41
42 Index pair (N, L) := ( G M (I), G M (I)\I ), where G M : P P (P a partition of the search box ) such that for C P we have Ĝ2(C) G M (C). Induced homology map: f k = 1 k = 3, 0 else. Λ((N, L), Ĝ2) = 1. Oliver Junge, Institute for Mathematics, University of Paderborn 42
43 Theorem 3 The differential equation with parameters σ = 2 and δ = 3 possesses a periodic orbit y(t) = k Z c k e ikωt with c k = c k and ω [1.36, 1.44] c 0 [2.43, 2.49] c 1 [0.8, ] c 2 [ 0.003, 0.04] [0.005, 0.06]i c k {z C z β α k }, β = 11.5, α = 0.1, k > 2. Oliver Junge, Institute for Mathematics, University of Paderborn 43
MATH 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 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 informationNumerical Continuation of Bifurcations - An Introduction, Part I
Numerical Continuation of Bifurcations - An Introduction, Part I given at the London Dynamical Systems Group Graduate School 2005 Thomas Wagenknecht, Jan Sieber Bristol Centre for Applied Nonlinear Mathematics
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 information2.10 Saddles, Nodes, Foci and Centers
2.10 Saddles, Nodes, Foci and Centers In Section 1.5, a linear system (1 where x R 2 was said to have a saddle, node, focus or center at the origin if its phase portrait was linearly equivalent to one
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
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 informationsatisfying the following condition: If T : V V is any linear map, then µ(x 1,,X n )= det T µ(x 1,,X n ).
ensities Although differential forms are natural objects to integrate on manifolds, and are essential for use in Stoke s theorem, they have the disadvantage of requiring oriented manifolds in order for
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
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 informationMath 868 Final Exam. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each). Y (φ t ) Y lim
SOLUTIONS Dec 13, 218 Math 868 Final Exam In this exam, all manifolds, maps, vector fields, etc. are smooth. Part 1. Complete 5 of the following 7 sentences to make a precise definition (5 points each).
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationPractice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009.
Practice Qualifying Exam Questions, Differentiable Manifolds, Fall, 2009. Solutions (1) Let Γ be a discrete group acting on a manifold M. (a) Define what it means for Γ to act freely. Solution: Γ acts
More informationChanging coordinates to adapt to a map of constant rank
Introduction to Submanifolds Most manifolds of interest appear as submanifolds of others e.g. of R n. For instance S 2 is a submanifold of R 3. It can be obtained in two ways: 1 as the image of a map into
More informationHalf of Final Exam Name: Practice Problems October 28, 2014
Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half
More informationLecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)
Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic
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 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 informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationApril 13, We now extend the structure of the horseshoe to more general kinds of invariant. x (v) λ n v.
April 3, 005 - Hyperbolic Sets We now extend the structure of the horseshoe to more general kinds of invariant sets. Let r, and let f D r (M) where M is a Riemannian manifold. A compact f invariant set
More informationIntroduction LECTURE 1
LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in
More informationNumerical techniques: Deterministic Dynamical Systems
Numerical techniques: Deterministic Dynamical Systems Henk Dijkstra Institute for Marine and Atmospheric research Utrecht, Department of Physics and Astronomy, Utrecht, The Netherlands Transition behavior
More information5.2.2 Planar Andronov-Hopf bifurcation
138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω
More informationAn Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013
An Introduction to Numerical Continuation Methods with Application to some Problems from Physics Eusebius Doedel Cuzco, Peru, May 2013 Persistence of Solutions Newton s method for solving a nonlinear equation
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 informationConvexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls
1 1 Convexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls B.T.Polyak Institute for Control Science, Moscow, Russia e-mail boris@ipu.rssi.ru Abstract Recently [1, 2] the new convexity
More informationNow I switch to nonlinear systems. In this chapter the main object of study will be
Chapter 4 Stability 4.1 Autonomous systems Now I switch to nonlinear systems. In this chapter the main object of study will be ẋ = f(x), x(t) X R k, f : X R k, (4.1) where f is supposed to be locally Lipschitz
More informationRemark on Hopf Bifurcation Theorem
Remark on Hopf Bifurcation Theorem Krasnosel skii A.M., Rachinskii D.I. Institute for Information Transmission Problems Russian Academy of Sciences 19 Bolshoi Karetny lane, 101447 Moscow, Russia E-mails:
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationLinear Ordinary Differential Equations
MTH.B402; Sect. 1 20180703) 2 Linear Ordinary Differential Equations Preliminaries: Matrix Norms. Denote by M n R) the set of n n matrix with real components, which can be identified the vector space R
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationtutorial ii: One-parameter bifurcation analysis of equilibria with matcont
tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session
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 informationWIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY
WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested
More informationNBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011
NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets
More informationA review of stability and dynamical behaviors of differential equations:
A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary
More informationDIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17
DIFFERENTIAL GEOMETRY, LECTURE 16-17, JULY 14-17 6. Geodesics A parametrized line γ : [a, b] R n in R n is straight (and the parametrization is uniform) if the vector γ (t) does not depend on t. Thus,
More informationManifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds
MA 755 Fall 05. Notes #1. I. Kogan. Manifolds and tangent bundles. Vector fields and flows. 1 Differential manifolds, smooth maps, submanifolds Definition 1 An n-dimensional C k -differentiable manifold
More informationIntroduction to Bifurcation and Normal Form theories
Introduction to Bifurcation and Normal Form theories Romain Veltz / Olivier Faugeras October 9th 2013 ENS - Master MVA / Paris 6 - Master Maths-Bio (2013-2014) Outline 1 Invariant sets Limit cycles Stable/Unstable
More informationODE, part 2. Dynamical systems, differential equations
ODE, part 2 Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Dynamical systems, differential equations Consider a system of n first order equations du dt = f(u, t),
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationSeptember 21, :43pm Holm Vol 2 WSPC/Book Trim Size for 9in by 6in
1 GALILEO Contents 1.1 Principle of Galilean relativity 2 1.2 Galilean transformations 3 1.2.1 Admissible force laws for an N-particle system 6 1.3 Subgroups of the Galilean transformations 8 1.3.1 Matrix
More informationMath 273 (51) - Final
Name: Id #: Math 273 (5) - Final Autumn Quarter 26 Thursday, December 8, 26-6: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
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 informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationLECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS. 1. Lie groups
LECTURE 16: LIE GROUPS AND THEIR LIE ALGEBRAS 1. Lie groups A Lie group is a special smooth manifold on which there is a group structure, and moreover, the two structures are compatible. Lie groups are
More informationIndex Theory and Periodic Solution of Delay Differential Systems
Index Theory and Periodic Solution of Delay Differential Systems Liu Chungen School of Mathematics, Nankai University workshop on nonlinear PDE and calculus of variation, Chern Mathematical Institute 2013.9.16
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationVariational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University)
Variational construction of periodic and connecting orbits in the planar Sitnikov problem Mitsuru Shibayama(Kyoto University) 1 The three-body problem Consider the planar three-body problem which is governed
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationWeak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations
Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak
More informationComputational Methods in Dynamical Systems and Advanced Examples
and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How
More informationFinding periodic orbits in state-dependent delay differential equations as roots of algebraic equations
Finding periodic orbits in state-dependent delay differential equations as roots of algebraic equations Jan Sieber arxiv:1010.2391v8 [math.ds] 26 Feb 2015 July 6, 2017 In this paper we prove that periodic
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
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 informationRecall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm
Chapter 13 Radon Measures Recall that if X is a compact metric space, C(X), the space of continuous (real-valued) functions on X, is a Banach space with the norm (13.1) f = sup x X f(x). We want to identify
More informationMechanical Systems II. Method of Lagrange Multipliers
Mechanical Systems II. Method of Lagrange Multipliers Rafael Wisniewski Aalborg University Abstract. So far our approach to classical mechanics was limited to finding a critical point of a certain functional.
More information1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),
1.7. Stability and attractors. Consider the autonomous differential equation (7.1) ẋ = f(x), where f C r (lr d, lr d ), r 1. For notation, for any x lr d, c lr, we let B(x, c) = { ξ lr d : ξ x < c }. Suppose
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationREPRESENTATION THEORY WEEK 7
REPRESENTATION THEORY WEEK 7 1. Characters of L k and S n A character of an irreducible representation of L k is a polynomial function constant on every conjugacy class. Since the set of diagonalizable
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 373-38, 23. Copyright 23,. ISSN 68-963. ETNA A NOTE ON THE RELATION BETWEEN THE NEWTON HOMOTOPY METHOD AND THE DAMPED NEWTON METHOD XUPING ZHANG
More informationHopf bifurcation in coupled cell networks with abelian symmetry 1
Hopf bifurcation in coupled cell networks with abelian symmetry 1 Ana Paula S. Dias e Rui C. Paiva CMUP Departamento de Matemática Universidade do Porto CMUP Instituto Politécnico de Leiria Abstract We
More informationMath 598 Feb 14, Geometry and Topology II Spring 2005, PSU
Math 598 Feb 14, 2005 1 Geometry and Topology II Spring 2005, PSU Lecture Notes 7 2.7 Smooth submanifolds Let N be a smooth manifold. We say that M N m is an n-dimensional smooth submanifold of N, provided
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More information1 )(y 0) {1}. Thus, the total count of points in (F 1 (y)) is equal to deg y0
1. Classification of 1-manifolds Theorem 1.1. Let M be a connected 1 manifold. Then M is diffeomorphic either to [0, 1], [0, 1), (0, 1), or S 1. We know that none of these four manifolds are not diffeomorphic
More informationSymplectic critical surfaces in Kähler surfaces
Symplectic critical surfaces in Kähler surfaces Jiayu Li ( Joint work with X. Han) ICTP-UNESCO and AMSS-CAS November, 2008 Symplectic surfaces Let M be a compact Kähler surface, let ω be the Kähler form.
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 informationINTEGRATION ON MANIFOLDS and GAUSS-GREEN THEOREM
INTEGRATION ON MANIFOLS and GAUSS-GREEN THEOREM 1. Schwarz s paradox. Recall that for curves one defines length via polygonal approximation by line segments: a continuous curve γ : [a, b] R n is rectifiable
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationElements of Applied Bifurcation Theory
Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2
More information1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system:
1. (i) Determine how many periodic orbits and equilibria can be born at the bifurcations of the zero equilibrium of the following system: ẋ = y x 2, ẏ = z + xy, ż = y z + x 2 xy + y 2 + z 2 x 4. (ii) Determine
More informationODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class
ODE Homework Due Wed. 9 August 2009; At the beginning of the class. (a) Solve Lẏ + Ry = E sin(ωt) with y(0) = k () L, R, E, ω are positive constants. (b) What is the limit of the solution as ω 0? (c) Is
More information225A DIFFERENTIAL TOPOLOGY FINAL
225A DIFFERENTIAL TOPOLOGY FINAL KIM, SUNGJIN Problem 1. From hitney s Embedding Theorem, we can assume that N is an embedded submanifold of R K for some K > 0. Then it is possible to define distance function.
More informationISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM
ISOLATING BLOCKS NEAR THE COLLINEAR RELATIVE EQUILIBRIA OF THE THREE-BODY PROBLEM RICHARD MOECKEL Abstract. The collinear relative equilibrium solutions are among the few explicitly known periodic solutions
More informationHomoclinic saddle to saddle-focus transitions in 4D systems
Faculty of Electrical Engineering, Mathematics & Computer Science Homoclinic saddle to saddle-focus transitions in 4D systems Manu Kalia M.Sc. Thesis July 2017 Assessment committee: Prof. Dr. S. A. van
More informationNotes for Math 535 Differential Geometry Spring Francis Bonahon. Department of Mathematics, University of Southern California
Notes for Math 535 Differential Geometry Spring 2016 Francis Bonahon Department of Mathematics, University of Southern California Date of this version: April 27, 2016 c Francis Bonahon 2016 CHAPTER 1 A
More informationTHE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS
THE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS GUILLAUME LAJOIE Contents 1. Introduction 2 2. The Hartman-Grobman Theorem 2 2.1. Preliminaries 2 2.2. The discrete-time Case 4 2.3. The
More informationWe denote the derivative at x by DF (x) = L. With respect to the standard bases of R n and R m, DF (x) is simply the matrix of partial derivatives,
The derivative Let O be an open subset of R n, and F : O R m a continuous function We say F is differentiable at a point x O, with derivative L, if L : R n R m is a linear transformation such that, for
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationQUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday August 31, 2010 (Day 1)
QUALIFYING EXAMINATION Harvard University Department of Mathematics Tuesday August 31, 21 (Day 1) 1. (CA) Evaluate sin 2 x x 2 dx Solution. Let C be the curve on the complex plane from to +, which is along
More informationRigorous Numerics for Global Dynamics: a study of the Swift-Hohenberg equation
Rigorous Numerics for Global Dynamics: a study of the Swift-Hohenberg equation Sarah Day Yasuaki Hiraoka Konstantin Mischaikow Toshiyuki Ogawa February 24, 2004 Abstract This paper presents a rigorous
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationMIDTERM REVIEW AND SAMPLE EXAM. Contents
MIDTERM REVIEW AND SAMPLE EXAM Abstract These notes outline the material for the upcoming exam Note that the review is divided into the two main topics we have covered thus far, namely, ordinary differential
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More information1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal
. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal 3 9 matrix D such that A = P DP, for A =. 3 4 3 (a) P = 4, D =. 3 (b) P = 4, D =. (c) P = 4 8 4, D =. 3 (d) P
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationDUALITY, OPTIMALITY CONDITIONS AND PERTURBATION ANALYSIS
1 DUALITY, OPTIMALITY CONDITIONS AND PERTURBATION ANALYSIS Alexander Shapiro 1 School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA, E-mail: ashapiro@isye.gatech.edu
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 0946 8633 On the stability of periodic orbits in delay equations with large delay
More informationLECTURE 15: COMPLETENESS AND CONVEXITY
LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other
More informationProblem Set Number 5, j/2.036j MIT (Fall 2014)
Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:
More informationIsothermal parameters
MTH.B402; Sect. 3 (20180626) 26 Isothermal parameters A Review of Complex Analysis. Let C be the complex plane. A C 1 -function 7 f : C D z w = f(z) C defined on a domain D is said to be holomorphic if
More information2 Lecture 2: Amplitude equations and Hopf bifurcations
Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From
More informationStability of flow past a confined cylinder
Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder
More informationMA651 Topology. Lecture 10. Metric Spaces.
MA65 Topology. Lecture 0. Metric Spaces. This text is based on the following books: Topology by James Dugundgji Fundamental concepts of topology by Peter O Neil Linear Algebra and Analysis by Marc Zamansky
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 informationBifurcation of Fixed Points
Bifurcation of Fixed Points CDS140B Lecturer: Wang Sang Koon Winter, 2005 1 Introduction ẏ = g(y, λ). where y R n, λ R p. Suppose it has a fixed point at (y 0, λ 0 ), i.e., g(y 0, λ 0 ) = 0. Two Questions:
More information