Set oriented numerics
|
|
- Irma Peters
- 5 years ago
- Views:
Transcription
1 Set oriented numerics Oliver Junge Center for Mathematics Technische Universität München Applied Koopmanism, MFO, February 2016
2 Motivation Complicated dynamics ( chaos ) Edward N. Lorenz,
3 Sensitive dependence on initial conditions Henri Poincaré Science et Méthodes, 1908 Il peut arriver que de petites différences dans les conditions initiales en engendrent de trés grandes dans les phénomènes finaux. Un petite erreur sur les premières produirait une erreur énorme sur les dernières. La prédiction devient impossible et nous avons le phénomène fortuit... Henri Poincaré,
4 Motivation Conformation dynamics Picture: L.E. Kavraki (alanine dipeptide) Dynamics of the distance between two CH 2 groups in an MD simulation of a polymer chain of 100 CH 2 groups (Grubmüller, Tavan, 94).
5 Set oriented numerics M. Dellnitz, A. Hohmann, A subdivision algorithm for the computation of unstable manifolds and global attractors, Num. Math., 75(3): , Related approaches: Hsu ( cell mapping ), Osipenko
6 Attractors X R d state space f : X X continuous (e.g., f is a time-sampling of the flow of an ODE). For x 0 X, consider the iteration Q X compact, the set x k+1 = f (x k ), k = 0, 1, 2,... A Q = f k (Q) k=0 is the global attractor (relative to Q).
7 Attractors Properties of the (relative) global attractor A Q contains all invariant sets in Q A Q contains (part of) the unstable set of an invariant set in Q
8 Set oriented computation of the global attractor input: Q X compact output: sequence B 0, B 1,... of finite collections of compact subsets of Q. Set B 0 = {Q}. Given B k 1, B k is constructed in two steps: Subdivide each cell B B k 1 to cells ˆB ˆB k such that diam( ˆB) θ diam(b), θ (0, 1). Define B k = { ˆB ˆB k : f ( ˆB ) ˆB for some ˆB ˆB k }.
9 y Example: Global attractor of the Hénon-map f (x, y) = (1 ax 2 + y, bx) x
10 Example: Global attractor in the Lorenz system
11 Convergence Let obviously Q k+1 Q k. Q k = B B k B, Theorem (Dellnitz, Hohmann, 97) The sets Q k converge to A Q in Hausdorff-distance as k, i.e. lim h(q k, A Q ) = 0. k Equivalently let Q = k=0,1,... Q k, then Q = A Q.
12 Speed of convergence Theorem (Dellnitz, Hohmann, 97) Let A Q be the global attractor relative to Q with respect to the diffeomorphism f q, q N. If A Q is attractive and hyperbolic and [...] then h(a Q, Q k ) diam(b k )(1 + α + + α k ), where If α < 1 then α = C λq θ min. 1 h(a Q, Q k ) diam(b k ) 1 α.
13 Set oriented computation of chain recurrent sets M. Dellnitz, O. Junge, M. Rumpf, R. Strzodka: The computation of an unstable invariant set inside a cylinder containing a knotted flow, Proc. EQUADIFF, 1999.
14 Software: GAIO Global Analysis of Invariant Objects C-library and Matlab toolbox set collections: set of boxes B(c, r) = {x : x i c i r i, i = 1,..., n} subdivison by bisection storage in a binary tree Root
15 The transfer operator on L 1 : T µ(a) = µ(f 1 (A)) T h (x) = h ( f 1 x ) Df 1 x, properties (Markov operator) (i) T is linear, (ii) T h 0 if h 0, (iii) T h 1 = h 1 for all h 0. T 1 1, hence σ(t ) B 1 (0) fixed points are invariant measures (resp. densities)
16 Theory Example: invariant measure and almost invariant sets in the Lorenz system
17 Numerics Ulam s method partition B 1,..., B n of M (B i connected, positive volume) approximation space: n = span{χ B1,..., χ Bn } projection π n : L 1 n, π n f = n i=1 1 f dm m(b i ) B i approximate Frobenius-Perron operator: T t n := π n T t matrix representation (Tn t ) ij = m(b i Φ t B j ). m(b j )
18 Convergence h n fixed point of T n, Theorem (Li, 76) Let f : [0, 1] [0, 1] be a piecewise C 2 -map, such that inf x f (x) > 2. If f has a unique absolutely continuous invariant measure µ with density h, then h n h 1 0 as n.
19 Random perturbations Replace the deterministic evolution x k+1 = f (x k ), k = 0, 1,... by the randomly perturbed evolution x k+1 = f (x k )+ξ k, k = 0, 1,... where ξ k is chosen randomly from some ball B ε (0) according to (e.g) a uniform distribution.
20 Random perturbations Stochastic transition function p : X A [0, 1], where for every x the probability measure p(x, ) determines the distribution of the perturbation ξ k (for all k). Example: p(x, A) = p ε (x, A) = m(a B ε(f (x))). m(b ε (0))
21 Small random perturbations Definition A family p ε of stochastic transition functions is a small random perturbation of (the deterministic system) f, if lim sup ε 0 g(y) p ε (x, dy) g(f (x)) = 0. x X Theorem (Kifer, 86) Let p ε be a small random perturbation of f. For every ε let µ ε be the unique invariant measure of p ε, supported on a neighbourhood of some attractive hyperbolic invariant set of f. Let µ be the corresponding SRB-measure of f. Then (under additional assumptions) µ ε µ weakly as ε 0.
22 Convergence: the hyperbolic case Fix ε, consider p ε, transfer operator P ε, compact on L 2. discretization: P ε P ε d, fixed point uε d Theorem (Dellnitz, J., 99) Let µ be the unique SRB-measure supported on a topologically transitive hyperbolic attractor. Then lim lim ε 0 d uε d µ.
23 Background Almost invariant sets [ 1 ε ε T = ε 1 ε ] the two nodes are (1 ε)-almost invariant eigenvalues : 1 and = 0.98 eigenvectors: (1, 1) at 1 (1, 1) at 0.98 Note that 1 ε = 1 + λ 2. 2
24 Theory Example: invariant measure and almost invariant sets in the Lorenz system
25 Background Spectrum of transfer operator macroscopic dynamics eigenvalue 1 invariant measure eigenvalues λ 1 almost invariant sets A identification via eigenfunction: invariance ratio: A joint nodal domain of eigenfunction m(a Φ t (A)) m(a) λ Dellnitz, Junge: On the approximation of complicated dynamical behavior. SIAM J. Numerical Analysis Froyland, Dellnitz: Detecting and Locating Near-Optimal Almost-Invariant Sets and Cycles, SIAM J. Scientific Computing, 2003.
26 The transfer operator and its generator The infinitesimal generator {T t } t 0 continuous semigroup infinitesimal generator Af := lim t 0 T t f f t for f D(A), D(A) the subspace of X where the limit exists if the vector field F is C 1 Af = div(f F )
27 Numerics Ulam s method for the generator same partition, approximation space, projection n D(A), thus A n := π n A does not work instead define matrix representation lim t 0 (A n ) ij = Tn t f f A n f := lim, f n t 0 t lim t 0 m(b i Φ t B j ), t m(b j ) i j; m(b i Φ t B i ) m(b i ), t m(b i ) otherwise.
28 Numerics Spectral collocation need explicit noise for a well posed problem A ε = ε + A eigenfunctions of A ε are C approximation spaces V n M = S 1 : trigonometric polynomials collocation nodes: uniform grid M = [ 1, 1]: Chebyshev polynomials collocation nodes: Chebyshev grid I n : C V n interpolation in collocation nodes approximate generator I n A ε f, f V n
29 Example 1: A flow on the circle invariant density vector field F (x) = sin(4πx) unique invariant density with C such that f 1 = 1. f (x) = C/F (x)
30 Example 1: A flow on the circle Efficiency L 1 error histogram Ulam s method Ulam for the generator spectral collocation # of evaluations
31 Example 2: An area-preserving cylinder flow The vector field 8 x y x x 10 5 periodic in x, 0 at the y-boundaries
32 Example 2: An area-preserving cylinder flow The fourth eigenvector
33 Example 2: An area-preserving cylinder flow The largest eigenvalues method λ 2 λ 3 λ 4 Ulam (log(λ i )/T ) Ulam (generator) spectral collocation Numerical effort method # of rhs evals time matrix time eigs Ulam sec. 1.0 sec. Ulam (generator) sec. 0.8 sec. spectral collocation sec. 3 sec.
34 Example 3: A volume-preserving 3D flow: the ABC-flow The vector field ẋ = a sin(2πz) + c cos(2πy) ẏ = b sin(2πx) + a cos(2πz) ż = c sin(2πy) + b cos(2πx), with a = 3, b = 2 and c = 1.
35 Example 3: A volume-preserving 3D flow: the ABC-flow Almost invariant sets Ulam s method for the generator Second eigenvector Fifth eigenvector
36 Example 3: A volume-preserving 3D flow: the ABC-flow Almost invariant sets spectral collocation Second eigenvector Fifth eigenvector
37 Example 3: A volume-preserving 3D flow: the ABC-flow Efficiency Ulam s method for the generator boxes Gauss quadrature with 16 nodes on each face rhs evaluations error in the invariant density spectral collocation grid 10 3 rhs evaluations error in the invariant density 10 14
38 Motivation Challenges in controller design highly nonlinear hybrid, quantized, event-based
39 Optimization based controller design time-discrete control system x k+1 = f (x k, u k ), k = 0, 1,..., cost function: c(x, u) > 0 for x T, c(t, u) = 0. accumulated cost: J(x 0, (u k ) k ) = c(x k, u k ) k=0 value function: V (x) = inf (u k ) k J(x, (u k ) k ) optimality principle/bellman operator V (x) = inf {c(x, u) + V (f (x, u))}, V (0) = 0. u U optimal feedback: F (x) = argmin u U {c(x, u) + V (f (x, u))}
40 Quantization of the state information P N P 1 P 2 x [x] Partition P = {P 1,..., P N } of state space X x [x] P quantized state approximation space for the value function: A = span{χ P1,..., χ PN }
41 Discrete problem & solution (I) idea: projection onto A discrete optimality principle on A V ([x]) = inf inf {c(ˆx, u) + V (f (ˆx, u))}, V ([0]) = 0 ˆx [x] u U = shortest path problem on a directed weighted graph solution: Dijkstra s algorithm Theorem (J., Osinga, 04) The discrete value function converges pointwise to the true one.
42 Discrete problem & solution (II) idea: dynamic game, discretization as perturbation discrete optimality principle on A V ([x]) = inf sup {c(ˆx, u) + V (f (ˆx, u))}, V ([0]) = 0 u U ˆx [x] = shortest path problem on a directed weighted hypergraph solution: minmax-dijkstra algorithm [Grüne, J., 08] Theorem (Grüne, J., 07) The discrete value function is a Lyapunov function for the closed loop system.
43 Convergence Theorem (Grüne, J., 08) P i nested sequence of partitions with diam P i 0 c continuous, c(x, u) > 0 for x T V continuous on T Then V Pi V uniformly on sets on which V is continuous. If the set of discontinuities of V has Lebesgue measure 0, then V Pi V in L 1 on compact subsets of the domain of V.
44 Example: An inverted pendulum The model m ϕ u state: x = (ϕ, ϕ) control system: f (x, u) = Φ T (x, u) cost function c(x, u) = T 0 q 1 ϕ 2 (t) + q 2 ϕ 2 (t) dt + Tq 3 u 2
45 Example: An inverted pendulum Value function and feedback trajectory
46 Hybrid systems x k+1 = f c (x k, y k, u k ) y k+1 = f d (x k, y k, u k ) k = 0, 1,... x k X R d compact, y k Y finite 1. no discrete dynamics, control discrete (i.e. U finite) x k+1 = f c (x k, u k ) 2. state dependent switched system x k+1 = f c (x k, y k, u k ) y k+1 = f d (x k ) 3. x controlled by automaton, which is controlled by u x k+1 = f c (x k, y k ) y k+1 = f d (y k, u k )
47 Example: a switched voltage converter The model [Lincoln, Rantzer, 06] ẋ 1 = 1 C (x 2 I load ) ẋ 2 = 1 L x 1 R L x L uv in control: u { 1, 1}.
48 Example: a switched voltage converter Feedback trajectory 0.6 Voltage Current Switch sign Sample
49 Set oriented numerics is good for... computing low-dimensional sets (even in high dimensional state spaces) attractors, invariant sets (dynamical systems) manifolds (mission design) Pareto sets (global multi-objective optimization) computing non-smooth functions on low-dimensional and irregular state spaces value functions (controller design) stability regions transport rates (astrodynamics) incorporating bounded uncertainties in a worst case way software: GAIO (new version 3.0 out soon) book in preparation
Zur Konstruktion robust stabilisierender Regler mit Hilfe von mengenorientierten numerischen Verfahren
Zur Konstruktion robust stabilisierender Regler mit Hilfe von mengenorientierten numerischen Verfahren Oliver Junge Zentrum Mathematik Technische Universität München mit Lars Grüne, Bayreuth Outline optimal
More informationAn adaptive subdivision technique for the approximation of attractors and invariant measures. Part II: Proof of convergence
An adaptive subdivision technique for the approximation of attractors and invariant measures. Part II: Proof of convergence Oliver Junge Department of Mathematics and Computer Science University of Paderborn
More informationSet Oriented Numerical Methods for Dynamical Systems
Set Oriented Numerical Methods for Dynamical Systems Michael Dellnitz and Oliver Junge Department of Mathematics and Computer Science University of Paderborn D-3395 Paderborn http://www.upb.de/math/~agdellnitz
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 informationTOWARDS AUTOMATED CHAOS VERIFICATION
TOWARDS AUTOMATED CHAOS VERIFICATION SARAH DAY CDSNS, Georgia Institute of Technology, Atlanta, GA 30332 OLIVER JUNGE Institute for Mathematics, University of Paderborn, 33100 Paderborn, Germany KONSTANTIN
More informationA tutorial on numerical transfer operator methods. numerical analysis of dynamical systems
A tutorial on transfer operator methods for numerical analysis of dynamical systems Gary Froyland School of Mathematics and Statistics University of New South Wales, Sydney BIRS Workshop on Uncovering
More informationDynamical system theory and numerical methods applied to Astrodynamics
Dynamical system theory and numerical methods applied to Astrodynamics Roberto Castelli Institute for Industrial Mathematics University of Paderborn BCAM, Bilbao, 20th December, 2010 Introduction Introduction
More informationDynamic programming using radial basis functions and Shepard approximations
Dynamic programming using radial basis functions and Shepard approximations Oliver Junge, Alex Schreiber Fakultät für Mathematik Technische Universität München Workshop on algorithms for dynamical systems
More informationPHY411 Lecture notes Part 5
PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
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 informationProblemas abiertos en dinámica de operadores
Problemas abiertos en dinámica de operadores XIII Encuentro de la red de Análisis Funcional y Aplicaciones Cáceres, 6-11 de Marzo de 2017 Wikipedia Old version: In mathematics and physics, chaos theory
More informationRepresentation of Markov chains. Stochastic stability AIMS 2014
Markov chain model Joint with J. Jost and M. Kell Max Planck Institute for Mathematics in the Sciences Leipzig 0th July 204 AIMS 204 We consider f : M! M to be C r for r 0 and a small perturbation parameter
More informationLecture 4. Entropy and Markov Chains
preliminary version : Not for diffusion Lecture 4. Entropy and Markov Chains The most important numerical invariant related to the orbit growth in topological dynamical systems is topological entropy.
More informationSRB measures for non-uniformly hyperbolic systems
SRB measures for non-uniformly hyperbolic systems Vaughn Climenhaga University of Maryland October 21, 2010 Joint work with Dmitry Dolgopyat and Yakov Pesin 1 and classical results Definition of SRB measure
More informationMean field coupling of expanding circle maps
Department of Stochastics, TU Budapest October 2, 205. Motivation José Koiller and Lai-Sang Young. Coupled map networks. Nonlinearity, 23(5):2, 200. Bastien Fernandez. Breaking of ergodicity in expanding
More informationOne dimensional Maps
Chapter 4 One dimensional Maps The ordinary differential equation studied in chapters 1-3 provide a close link to actual physical systems it is easy to believe these equations provide at least an approximate
More informationThe Existence of Chaos in the Lorenz System
The Existence of Chaos in the Lorenz System Sheldon E. Newhouse Mathematics Department Michigan State University E. Lansing, MI 48864 joint with M. Berz, K. Makino, A. Wittig Physics, MSU Y. Zou, Math,
More informationRigorous approximation of invariant measures for IFS Joint work April 8, with 2016 S. 1 Galat / 21
Rigorous approximation of invariant measures for IFS Joint work with S. Galatolo e I. Nisoli Maurizio Monge maurizio.monge@im.ufrj.br Universidade Federal do Rio de Janeiro April 8, 2016 Rigorous approximation
More informationPeriodic Sinks and Observable Chaos
Periodic Sinks and Observable Chaos Systems of Study: Let M = S 1 R. T a,b,l : M M is a three-parameter family of maps defined by where θ S 1, r R. θ 1 = a+θ +Lsin2πθ +r r 1 = br +blsin2πθ Outline of Contents:
More informationSet oriented construction of globally optimal controllers
at / Set oriented construction of globally optimal controllers Mengenorientierte Konstruktion global optimaler Regler Oliver Junge, Technische Universität München and Lars Grüne, Universität Bayreuth In
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 informationDuality in Stability Theory: Lyapunov function and Lyapunov measure
Forty-Fourth Annual Allerton Conference Allerton House, UIUC, Illinois, USA Sept 27-29, 2006 WID.150 Duality in Stability Theory: Lyapunov function and Lyapunov measure Umesh Vaidya Abstract In this paper,
More informationProperties for systems with weak invariant manifolds
Statistical properties for systems with weak invariant manifolds Faculdade de Ciências da Universidade do Porto Joint work with José F. Alves Workshop rare & extreme Gibbs-Markov-Young structure Let M
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 informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
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 informationDynamical Systems and Ergodic Theory PhD Exam Spring Topics: Topological Dynamics Definitions and basic results about the following for maps and
Dynamical Systems and Ergodic Theory PhD Exam Spring 2012 Introduction: This is the first year for the Dynamical Systems and Ergodic Theory PhD exam. As such the material on the exam will conform fairly
More informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationComputational Issues in Nonlinear Dynamics and Control
Computational Issues in Nonlinear Dynamics and Control Arthur J. Krener ajkrener@ucdavis.edu Supported by AFOSR and NSF Typical Problems Numerical Computation of Invariant Manifolds Typical Problems Numerical
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
More informationLozi-like maps. M. Misiurewicz and S. Štimac. May 13, 2017
Lozi-like maps M. Misiurewicz and S. Štimac May 13, 017 Abstract We define a broad class of piecewise smooth plane homeomorphisms which have properties similar to the properties of Lozi maps, including
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 informationConstrained Ulam Dynamic Mode Decomposition: Approximation of the Perron-Frobenius Operator for Deterministic and Stochastic Systems
Constrained Ulam Dynamic Mode Decomposition: Approximation of the Perron-Frobenius Operator for Deterministic and Stochastic Systems Debdipta Goswami, Emma Thackray 2, and Derek A. Paley 3 Abstract Dynamical
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 informationCalculating the domain of attraction: Zubov s method and extensions
Calculating the domain of attraction: Zubov s method and extensions Fabio Camilli 1 Lars Grüne 2 Fabian Wirth 3 1 University of L Aquila, Italy 2 University of Bayreuth, Germany 3 Hamilton Institute, NUI
More informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationSet oriented methods, invariant manifolds and transport
C C Dynamical A L T E C S H Set oriented methods, invariant manifolds and transport Shane D. Ross Control and Dynamical Systems, Caltech www.cds.caltech.edu/ shane Caltech/MIT/NASA: W.S. Koon, F. Lekien,
More informationLorenz like flows. Maria José Pacifico. IM-UFRJ Rio de Janeiro - Brasil. Lorenz like flows p. 1
Lorenz like flows Maria José Pacifico pacifico@im.ufrj.br IM-UFRJ Rio de Janeiro - Brasil Lorenz like flows p. 1 Main goals The main goal is to explain the results (Galatolo-P) Theorem A. (decay of correlation
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 1. Linear systems Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Linear systems 1.1 Differential Equations 1.2 Linear flows 1.3 Linear maps
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationScaling limits of anisotropic random growth models
Scaling limits of anisotropic random growth models Amanda Turner Department of Mathematics and Statistics Lancaster University (Joint work with Fredrik Johansson Viklund and Alan Sola) Overview 1 Generalised
More informationHybrid Control and Switched Systems. Lecture #8 Stability and convergence of hybrid systems (topological view)
Hybrid Control and Switched Systems Lecture #8 Stability and convergence of hybrid systems (topological view) João P. Hespanha University of California at Santa Barbara Summary Lyapunov stability of hybrid
More informationRigorous computation of topological entropy with respect to a finite partition
Physica D 154 (2001) 68 84 Rigorous computation of topological entropy with respect to a finite partition Gary Froyland a,1, Oliver Junge a,, Gunter Ochs b,2 a Department of Mathematics and Computer Science,
More informationTheory and applications of random Poincaré maps
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ Random Dynamical Systems Theory and applications of random Poincaré maps Nils Berglund MAPMO, Université d'orléans,
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 informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
More information1 Introduction Definitons Markov... 2
Compact course notes Dynamic systems Fall 2011 Professor: Y. Kudryashov transcribed by: J. Lazovskis Independent University of Moscow December 23, 2011 Contents 1 Introduction 2 1.1 Definitons...............................................
More informationComparison of systems with complex behavior
Physica D 197 (2004) 101 133 Comparison of systems with complex behavior Igor Mezić a,, Andrzej Banaszuk b a Department of Mechanical and Environmental Engineering, University of California, Santa Barbara,
More informationHeat kernels of some Schrödinger operators
Heat kernels of some Schrödinger operators Alexander Grigor yan Tsinghua University 28 September 2016 Consider an elliptic Schrödinger operator H = Δ + Φ, where Δ = n 2 i=1 is the Laplace operator in R
More informationIf Λ = M, then we call the system an almost Anosov diffeomorphism.
Ergodic Theory of Almost Hyperbolic Systems Huyi Hu Department of Mathematics, Pennsylvania State University, University Park, PA 16802, USA, E-mail:hu@math.psu.edu (In memory of Professor Liao Shantao)
More informationA FAMILY OF PIECEWISE EXPANDING MAPS HAVING SINGULAR MEASURE AS A LIMIT OF ACIM S
A FAMILY OF PIECEWISE EXPANDING MAPS HAVING SINGULAR MEASURE AS A LIMIT OF ACIM S ZHENYANG LI, PAWE L GÓ, ABRAHAM BOYARSKY, HARALD PROPPE, AND PEYMAN ESLAMI Abstract Keller [9] introduced families of W
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 informationThe Conley Index and Rigorous Numerics of Attracting Periodic Orbits
The Conley Index and Rigorous Numerics of Attracting Periodic Orbits Marian Mrozek Pawe l Pilarczyk Conference on Variational and Topological Methods in the Study of Nonlinear Phenomena (Pisa, 2000) 1
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
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 informationSHADOWING AND INVERSE SHADOWING IN SET-VALUED DYNAMICAL SYSTEMS. HYPERBOLIC CASE. Sergei Yu. Pilyugin Janosch Rieger. 1.
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 32, 2008, 151 164 SHADOWING AND INVERSE SHADOWING IN SET-VALUED DYNAMICAL SYSTEMS. HYPERBOLIC CASE Sergei Yu. Pilyugin
More informationInfinite-Dimensional Dynamical Systems in Mechanics and Physics
Roger Temam Infinite-Dimensional Dynamical Systems in Mechanics and Physics Second Edition With 13 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vii ix GENERAL
More informationApplied Dynamical Systems. 6 Statistical properties. 6.1 Probability measures 6 STATISTICAL PROPERTIES. This file was processed March 7, 2017.
Applied Dynamical Systems This file was processed March 7, 2017. 6 Statistical properties 6.1 Probability measures We now consider statistical properties of dynamical systems. Since for chaotic systems
More informationarxiv:chao-dyn/ v1 18 Dec 1997
1 arxiv:chao-dyn/9712016v1 18 Dec 1997 Random perturbations of chaotic dynamical systems. Stability of the spectrum Michael Blank, Gerhard Keller Mathematisches Institut, Universitat Erlangen-Nurnberg
More informationSAMPLE PATH PROPERIES OF RANDOM TRANSFORMATIONS.
SAMPLE PATH PROPERIES OF RANDOM TRANSFORMATIONS. DMITRY DOLGOPYAT 1. Models. Let M be a smooth compact manifold of dimension N and X 0, X 1... X d, d 2, be smooth vectorfields on M. (I) Let {w j } + j=
More informationarxiv: v3 [math.ds] 12 Jun 2013
Isostables, isochrons, and Koopman spectrum for the action-angle representation of stable fixed point dynamics A. Mauroy, I. Mezic, and J. Moehlis Department of Mechanical Engineering, University of California
More informationSOME EXAMPLES OF NONUNIFORMLY HYPERBOLIC COCYCLES. L.-S. Young 1
SOME EXAMPLES OF NONUNIFORMLY HYPERBOLIC COCYCLES L.-S. Young Abstract. We consider some very simple examples of SL(2, R)-cocycles and prove that they have positive Lyapunov exponents. These cocycles form
More informationUsing Ulam s method to calculate entropy and other dynamical invariants
Nonlinearity 2 (999) 79 0. Printed in the UK PII: S095-775(99)92923-7 Using Ulam s method to calculate entropy and other dynamical invariants Gary Froyland Department of Mathematical Engineering and Information
More informationMechanisms of Chaos: Stable Instability
Mechanisms of Chaos: Stable Instability Reading for this lecture: NDAC, Sec. 2.-2.3, 9.3, and.5. Unpredictability: Orbit complicated: difficult to follow Repeatedly convergent and divergent Net amplification
More informationDYNAMICAL SYSTEMS PROBLEMS. asgor/ (1) Which of the following maps are topologically transitive (minimal,
DYNAMICAL SYSTEMS PROBLEMS http://www.math.uci.edu/ asgor/ (1) Which of the following maps are topologically transitive (minimal, topologically mixing)? identity map on a circle; irrational rotation of
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 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 informationEigenfunctions for smooth expanding circle maps
December 3, 25 1 Eigenfunctions for smooth expanding circle maps Gerhard Keller and Hans-Henrik Rugh December 3, 25 Abstract We construct a real-analytic circle map for which the corresponding Perron-
More informationEigenvalues and eigenfunctions of the Laplacian. Andrew Hassell
Eigenvalues and eigenfunctions of the Laplacian Andrew Hassell 1 2 The setting In this talk I will consider the Laplace operator,, on various geometric spaces M. Here, M will be either a bounded Euclidean
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 informationdynamical zeta functions: what, why and what are the good for?
dynamical zeta functions: what, why and what are the good for? Predrag Cvitanović Georgia Institute of Technology November 2 2011 life is intractable in physics, no problem is tractable I accept chaos
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 informationDynamical Systems and Stochastic Processes.
Dynamical Systems and Stochastic Processes. P. Collet Centre de Physique Théorique CNRS UMR 7644 Ecole Polytechnique F-91128 Palaiseau Cedex France e-mail: collet@cpht.polytechnique.fr April 23, 2008 Abstract
More informationA conjecture on sustained oscillations for a closed-loop heat equation
A conjecture on sustained oscillations for a closed-loop heat equation C.I. Byrnes, D.S. Gilliam Abstract The conjecture in this paper represents an initial step aimed toward understanding and shaping
More informationPopulation Games and Evolutionary Dynamics
Population Games and Evolutionary Dynamics William H. Sandholm The MIT Press Cambridge, Massachusetts London, England in Brief Series Foreword Preface xvii xix 1 Introduction 1 1 Population Games 2 Population
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationAnalytic Continuation of Analytic (Fractal) Functions
Analytic Continuation of Analytic (Fractal) Functions Michael F. Barnsley Andrew Vince (UFL) Australian National University 10 December 2012 Analytic continuations of fractals generalises analytic continuation
More informationFeedback stabilization methods for the numerical solution of ordinary differential equations
Feedback stabilization methods for the numerical solution of ordinary differential equations Iasson Karafyllis Department of Environmental Engineering Technical University of Crete 731 Chania, Greece ikarafyl@enveng.tuc.gr
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 informationCoexistence of Zero and Nonzero Lyapunov Exponents
Coexistence of Zero and Nonzero Lyapunov Exponents Jianyu Chen Pennsylvania State University July 13, 2011 Outline Notions and Background Hyperbolicity Coexistence Construction of M 5 Construction of the
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationMARKOV PARTITIONS FOR HYPERBOLIC SETS
MARKOV PARTITIONS FOR HYPERBOLIC SETS TODD FISHER, HIMAL RATHNAKUMARA Abstract. We show that if f is a diffeomorphism of a manifold to itself, Λ is a mixing (or transitive) hyperbolic set, and V is a neighborhood
More informationPersistent Chaos in High-Dimensional Neural Networks
Persistent Chaos in High-Dimensional Neural Networks D. J. Albers with J. C. Sprott and James P. Crutchfield February 20, 2005 1 Outline: Introduction and motivation Mathematical versus computational dynamics
More informationMATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation.
MATH 614 Dynamical Systems and Chaos Lecture 24: Bifurcation theory in higher dimensions. The Hopf bifurcation. Bifurcation theory The object of bifurcation theory is to study changes that maps undergo
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationEfficient approximation methods for the global long-term behavior of dynamical systems Theory, algorithms and examples
Technische Universität München Zentrum Mathematik Efficient approximation methods for the global long-term behavior of dynamical systems Theory, algorithms and examples Péter Koltai Vollständiger Abdruck
More informationTHEOREMS, ETC., FOR MATH 516
THEOREMS, ETC., FOR MATH 516 Results labeled Theorem Ea.b.c (or Proposition Ea.b.c, etc.) refer to Theorem c from section a.b of Evans book (Partial Differential Equations). Proposition 1 (=Proposition
More informationConvergence of Feller Processes
Chapter 15 Convergence of Feller Processes This chapter looks at the convergence of sequences of Feller processes to a iting process. Section 15.1 lays some ground work concerning weak convergence of processes
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More information(Non-)Existence of periodic orbits in dynamical systems
(Non-)Existence of periodic orbits in dynamical systems Konstantin Athanassopoulos Department of Mathematics and Applied Mathematics University of Crete June 3, 2014 onstantin Athanassopoulos (Univ. of
More informationLecture 1 Monday, January 14
LECTURE NOTES FOR MATH 7394 ERGODIC THEORY AND THERMODYNAMIC FORMALISM INSTRUCTOR: VAUGHN CLIMENHAGA Lecture 1 Monday, January 14 Scribe: Vaughn Climenhaga 1.1. From deterministic physics to stochastic
More informationA short introduction with a view toward examples. Short presentation for the students of: Dynamical Systems and Complexity Summer School Volos, 2017
A short introduction with a view toward examples Center of Research and Applications of Nonlinear (CRANS) Department of Mathematics University of Patras Greece sanastassiou@gmail.com Short presentation
More informationIntroduction to fractal analysis of orbits of dynamical systems. ZAGREB DYNAMICAL SYSTEMS WORKSHOP 2018 Zagreb, October 22-26, 2018
Vesna Županović Introduction to fractal analysis of orbits of dynamical systems University of Zagreb, Croatia Faculty of Electrical Engineering and Computing Centre for Nonlinear Dynamics, Zagreb ZAGREB
More informationTHE POINT SPECTRUM OF FROBENIUS-PERRON AND KOOPMAN OPERATORS
PROCEEDINGS OF THE MERICN MTHEMTICL SOCIETY Volume 126, Number 5, May 1998, Pages 1355 1361 S 0002-9939(98)04188-4 THE POINT SPECTRUM OF FROBENIUS-PERRON ND KOOPMN OPERTORS J. DING (Communicated by Palle
More informationBasins of Attraction and Perturbed Numerical Solutions using Euler s Method
Basins of Attraction and Perturbed Numerical Solutions using Euler s Method Hendrik Orem 8 May 2008 1 Introduction 1.1 Motivation Any physical system corresponding to a two-dimensional vector field has
More information