A Verified ODE solver and Smale s 14th Problem
|
|
- Clyde Harvey
- 6 years ago
- Views:
Transcription
1 A Verified ODE solver and Smale s 14th Problem Fabian Immler Big Isaac Newton Institute Jul 6, 2017 = Isabelle λ β α
2 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important 1 / 30
3 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important computer-assisted proof : additional value of formal proof 1 / 30
4 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification final theorem considered important computer-assisted proof : additional value of formal proof risk of Mexican hat specific lemmas general purpose lemmas Mexican Hat [Geuvers, MAP 2004] via [Boldo, 2004] 1 / 30
5 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important computer-assisted proof : additional value of formal proof risk of Mexican hat HOL-Analysis 1 / 30
6 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important computer-assisted proof : additional value of formal proof risk of Mexican hat interesting mathematics (ODEs, dynamical systems) HOL-ODE HOL-Analysis 1 / 30
7 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important computer-assisted proof : additional value of formal proof risk of Mexican hat interesting mathematics (ODEs, dynamical systems) general purpose ODE solver Numerics HOL-ODE HOL-Analysis 1 / 30
8 Smale s 14th Problem A computer-assisted proof involving ODEs. Motivation for Formal Verification considered important computer-assisted proof : additional value of formal proof risk of Mexican hat interesting mathematics (ODEs, dynamical systems) general purpose ODE solver final theorem Numerics HOL-ODE HOL-Analysis 1 / 30
9 History 2 / 30
10 Lorenz Equations (1963, E.N. Lorenz) model atmospheric flows ẋ = 10(y x) ẏ = x(28 z) y ż = xy 8 3 z 3 / 30
11 Lorenz Equations (1963, E.N. Lorenz) model atmospheric flows ẋ = 10(y x) ẏ = x(28 z) y ż = xy 8 3 z numerical simulations finite precision, approximate computations 3 / 30
12 Lorenz Equations (1963, E.N. Lorenz) model atmospheric flows ẋ = 10(y x) ẏ = x(28 z) y ż = xy 8 3 z numerical simulations finite precision, approximate computations / 30
13 Lorenz Equations (1963, E.N. Lorenz) model atmospheric flows ẋ = 10(y x) ẏ = x(28 z) y ż = xy 8 3 z numerical simulations finite precision, approximate computations chaos, butterfly effect / 30
14 Lorenz Attractor numerical observations: chaotic, strange, hyperbolic 4 / 30
15 Lorenz Attractor numerical observations: chaotic, strange, hyperbolic the number of man, woman, and computer hours spent on [the Lorenz equations] [... ] must be truly immense [Sparrow, 1982] 4 / 30
16 Lorenz Attractor numerical observations: chaotic, strange, hyperbolic the number of man, woman, and computer hours spent on [the Lorenz equations] [... ] must be truly immense [Sparrow, 1982] no proofs 4 / 30
17 Geometric Model (1979) G Σ L G Geometric Lorenz Attractor [Williams, Guckenheimer, Yorke, 1979] 5 / 30
18 Problem 14: Lorenz Attractor Is the dynamics of the ordinary differential equations of Lorenz that of the geometric Lorenz attractor of Williams, Guckenheimer, and Yorke? From Steve Smale (1998): Mathematical Problems for the Next Century 6 / 30
19 Tucker s Proof (1998) ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * Poincaré section Σ 7 / 30
20 Tucker s Proof (1998) ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * Poincaré section Σ rigorous numerics for global behavior 7 / 30
21 Tucker s Proof (1998) ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * Poincaré section Σ rigorous numerics for global behavior propagation of interval enclosures (100 cpu hours) 7 / 30
22 Tucker s Proof (1998) ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * Poincaré section Σ rigorous numerics for global behavior propagation of interval enclosures (100 cpu hours) analysis (normal form theory) locally around the origin 7 / 30
23 8 / 30
24 Verification 9 / 30
25 Overview: Verified Rigorous Numerics for ODEs High Level: Reachability Analysis Medium Level: Enclosures/Expressions Lowest Level: Numbers 10 / 30
26 High Level: Reachability Analysis 11 / 30
27 Rigorous Enclosure for Returns ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * / 30
28 Intermediate Poincaré Sections / Splitting ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * 13 / 30
29 Hyperbolicity 1 #$ '# '$ &# &$ ( %# 1 %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * / 30
30 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) 15 / 30
31 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) RES S RES T S T 15 / 30
32 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) RES S RES T S T return x := RES {x} 15 / 30
33 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) RES S RES T S T return x := RES {x} spec P := RES {x. P x} 15 / 30
34 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) RES S RES T S T return x := RES {x} spec P := RES {x. P x} program correctness: x. P x = f x spec (λy. Q y) 15 / 30
35 Program Verification Framework (Automatic) Refinement in Isabelle/HOL [Lammich] nondeterministic result monad datatype α nres = FAIL RES (α set) RES S RES T S T return x := RES {x} spec P := RES {x. P x} program correctness: x. P x = f x spec (λy. Q y) transfer rules for refinement, e.g., (λ x. return ( f x), f ) R r S nres r 15 / 30
36 Inference Algorithm (à la Shankar) Operations: split-spec X := spec (λ(a, B). X A B) 16 / 30
37 Inference Algorithm (à la Shankar) Operations: split-spec X := spec (λ(a, B). X A B) flow-spec X t := spec (λr. φ(x, t) R) 16 / 30
38 Inference Algorithm (à la Shankar) Operations: split-spec X := spec (λ(a, B). X A B) flow-spec X t := spec (λr. φ(x, t) R) poincare-spec X Σ := spec (λi. P Σ (X ) I ) 16 / 30
39 Inference Algorithm (à la Shankar) Operations: split-spec X := spec (λ(a, B). X A B) flow-spec X t := spec (λr. φ(x, t) R) poincare-spec X Σ := spec (λi. P Σ (X ) I ) Invariant: x 0 X 0. t 0.φ(x 0, t) X P Σ (x 0 ) I 16 / 30
40 Medium Level: Enclosures and Rigorous Numerics 17 / 30
41 Enclosures not intervals: 18 / 30
42 Enclosures not intervals: but zonotopes: {l 0 + i ε i l i ε i [ 1; 1]} 18 / 30
43 Enclosures not intervals: but zonotopes: {l 0 + i ε i l i ε i [ 1; 1]} (l, {l 0 + i ε i l i ε i [ 1; 1]}) zonotope r 18 / 30
44 Enclosures not intervals: but zonotopes: {l 0 + i ε i l i ε i [ 1; 1]} (l, {l 0 + i ε i l i ε i [ 1; 1]}) zonotope r e.g., (split-zonotope, split-spec ) zonotope r r zonotope r r zonotope r nres r 18 / 30
45 Rigorous Numerics aexp = Add aexp aexp Mult aexp aexp Minus aexp Inverse aexp Num R Var N... [[Add a b]] vs = [[Mult a b]] vs = [[Minus a]] vs = [[Inverse a]] vs = [[Num r]] vs = [[Var i]] vs =... [[a]] vs + [[b]] vs [[a]] vs [[b]] vs [[a]] vs 1/[[a]] vs r vs! i Rigorous Numerics: approx-spec aexp X = spec (λr. x X. [[aexp ]] x R) Refinement: (approx, approx-spec ) Id r zonotope r r zonotope r nres r 19 / 30
46 Runge-Kutta methods x 0 ϕ O(h 2 ) h f (x 0 ) 0 h ODE: φ(t) = f (φ(t)) 20 / 30
47 Runge-Kutta methods ϕ O(h 2 ) h f (x 0 ) x 0 0 h ODE: φ(t) = f (φ(t)) Expressions: [[Euler aexp]] x0,h,ξ x 0 + h f (x 0 ) + O(h 2 ) [[RK aexp]] x0,h,ξ x 0 + h + O(h 3 ) 20 / 30
48 Runge-Kutta methods ϕ O(h 2 ) h f (x 0 ) x 0 0 h ODE: φ(t) = f (φ(t)) Expressions: [[Euler aexp]] x0,h,ξ x 0 + h f (x 0 ) + O(h 2 ) [[RK aexp]] x0,h,ξ x 0 + h + O(h 3 ) Mathematics: φ(x 0, h) x 0 + h f (x 0 ) O(h 2 ) φ(x 0, h)... O(h 3 ) 20 / 30
49 Runge-Kutta methods ϕ O(h 2 ) h f (x 0 ) x 0 0 h ODE: φ(t) = f (φ(t)) Expressions: [[Euler aexp]] x0,h,ξ x 0 + h f (x 0 ) + O(h 2 ) [[RK aexp]] x0,h,ξ x 0 + h + O(h 3 ) Mathematics: φ(x 0, h) x 0 + h f (x 0 ) O(h 2 ) φ(x 0, h)... O(h 3 ) Rigorous: approx-spec (RK aexp) (X, t) flow-spec X t 20 / 30
50 Derivatives w.r.t. initial value / 30
51 Derivatives w.r.t. initial value / 30
52 Derivatives w.r.t. initial value / 30
53 Derivatives w.r.t. initial value / 30
54 Derivatives w.r.t. initial value ODE f : R n R n : ẋ(t) = f (x) 22 / 30
55 Derivatives w.r.t. initial value ODE f : R n R n : Variational Equation: ẋ(t) = f (x) W (t) = Df (x(t)) W (t) 22 / 30
56 Derivatives w.r.t. initial value ODE f : R n R n : Variational Equation: ẋ(t) = f (x) W (t) = Df (x(t)) W (t) W (t) :: R n2 22 / 30
57 Derivatives w.r.t. initial value ODE f : R n R n : Variational Equation: ẋ(t) = f (x) W (t) :: R n2 W (t) = Df (x(t)) W (t) Algorithmic: variational equation on the level of expressions 22 / 30
58 Derivatives w.r.t. initial value ODE f : R n R n : ẋ(t) = f (x) Variational Equation: W (t) :: R n2 W (t) = Df (x(t)) W (t) Algorithmic: variational equation on the level of expressions Mathematics: W (t) = φ(x, t) x :: (R n l R n ) 22 / 30
59 Lowest Level 23 / 30
60 Numbers Software Floating Point Numbers m 2 e, m, e Z Explicit Rounding p (x) := x 2 p log 2 x 2 p log 2 x 24 / 30
61 Computations 25 / 30
62 Overview ( #$ '# '$ &# &$ %# %$ "# "$ # $!"#!"$!# $ # "$!# $ # "$ +, "#!"$ ) * 26 / 30
63 Computation 300 initial rectangles 27 / 30
64 Computation 300 initial rectangles with a 20 core machine: 27 / 30
65 Computation 300 initial rectangles with a 20 core machine: 155h cpu time 27 / 30
66 Computation 300 initial rectangles with a 20 core machine: 155h cpu time factor / 30
67 Computation 300 initial rectangles with a 20 core machine: 155h cpu time factor 17 9h elapsed time 27 / 30
68 Computation 300 initial rectangles with a 20 core machine: 155h cpu time factor 17 9h elapsed time without C 1 : speedup of factor / 30
69 / 30
70 / 30
71 / 30
72 / 30
73 / 30
74 28 / 30
75 28 / 30
76 28 / 30
77 / 30
78 / 30
79 / 30
80 by eval 29 / 30
81 Thank You. 30 / 30
82 Conclusion Many compromises: intervals/zonotopes/taylor models Euler/Runge-Kutta/Taylor series (linear) variational equation with generic solver software/hardware floating point numbers But modular and good enough. 31 / 30
The 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 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 informationThe Flow of ODEs. Fabian Immler and Christoph Traut. Institut für Informatik, Technische Universität München {immler,
The Flow of ODEs Fabian Immler and Christoph Traut Institut für Informatik, Technische Universität München {immler, trautc}@in.tum.de Abstract. Formal analysis of ordinary differential equations (ODEs)
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
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 informationRound-off Error Analysis of Explicit One-Step Numerical Integration Methods
Round-off Error Analysis of Explicit One-Step Numerical Integration Methods 24th IEEE Symposium on Computer Arithmetic Sylvie Boldo 1 Florian Faissole 1 Alexandre Chapoutot 2 1 Inria - LRI, Univ. Paris-Sud
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 informationSimple ODE Solvers - Derivation
Simple ODE Solvers - Derivation These notes provide derivations of some simple algorithms for generating, numerically, approximate solutions to the initial value problem y (t =f ( t, y(t y(t 0 =y 0 Here
More informationEuler s Method applied to the control of switched systems
Euler s Method applied to the control of switched systems FORMATS 2017 - Berlin Laurent Fribourg 1 September 6, 2017 1 LSV - CNRS & ENS Cachan L. Fribourg Euler s method and switched systems September
More informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
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 informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationComputed Chaos or Numerical Errors. Lun-Shin Yao
Computed Chaos or Numerical Errors Lun-Shin Yao Department of Mechanical and Aerospace Engineering Arizona State University Tempe, Arizona E-mail: ls_yao@asu.edu Abstract Discrete numerical methods with
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 informationEdward Lorenz: Predictability
Edward Lorenz: Predictability Master Literature Seminar, speaker: Josef Schröttle Edward Lorenz in 1994, Northern Hemisphere, Lorenz Attractor I) Lorenz, E.N.: Deterministic Nonperiodic Flow (JAS, 1963)
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
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 informationOn the numerical integration of variational equations
On the numerical integration of variational equations E. Gerlach 1, S. Eggl 2 and H. Skokos 3 1 Lohrmann Observatory, Technical University Dresden 2 Institute for Astronomy, University of Vienna 3 Max
More informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationHigh-Order Representation of Poincaré Maps
High-Order Representation of Poincaré Maps Johannes Grote, Martin Berz, and Kyoko Makino Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA. {grotejoh,berz,makino}@msu.edu
More informationLecture 1: Introduction, history, dynamics, nonlinearity, 1-D problem, phase portrait
Lecture 1: Introduction, history, dynamics, nonlinearity, 1-D problem, phase portrait Dmitri Kartofelev, PhD Tallinn University of Technology, School of Science, Department of Cybernetics, Laboratory of
More informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationThe Reflection Theorem
The Reflection Theorem Formalizing Meta-Theoretic Reasoning Lawrence C. Paulson Computer Laboratory Lecture Overview Motivation for the Reflection Theorem Proving the Theorem in Isabelle Applying the Reflection
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationODE Runge-Kutta methods
ODE Runge-Kutta methods The theory (very short excerpts from lectures) First-order initial value problem We want to approximate the solution Y(x) of a system of first-order ordinary differential equations
More informationFORECASTING ECONOMIC GROWTH USING CHAOS THEORY
Article history: Received 22 April 2016; last revision 30 June 2016; accepted 12 September 2016 FORECASTING ECONOMIC GROWTH USING CHAOS THEORY Mihaela Simionescu Institute for Economic Forecasting of the
More informationSynchronization, Chaos, and the Dynamics of Coupled Oscillators. Supplemental 1. Winter Zachary Adams Undergraduate in Mathematics and Biology
Synchronization, Chaos, and the Dynamics of Coupled Oscillators Supplemental 1 Winter 2017 Zachary Adams Undergraduate in Mathematics and Biology Outline: The shift map is discussed, and a rigorous proof
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationThe family of Runge Kutta methods with two intermediate evaluations is defined by
AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will
More informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationVerification of analog and mixed-signal circuits using hybrid systems techniques
FMCAD, November 2004, Austin Verification of analog and mixed-signal circuits using hybrid systems techniques Thao Dang, Alexandre Donze, Oded Maler VERIMAG Grenoble, France Plan 1. Introduction 2. Verification
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationA Rigorous Implementation of Taylor Model-Based Integrator in Chebyshev Basis Representation 1
A Rigorous Implementation of Taylor Model-Based Integrator in Chebyshev Basis Representation 1 Tomáš Dzetkulič Institute of Computer Science Czech Academy of Sciences December 15, 2011 1 This work was
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationTaylor Model Range Bounding Schemes
Taylor Model Range Bounding Schemes I Kyoko Makino II Martin Berz Department of Physics and Astronomy Michigan State University Introduction LDB (Linear Dominated Bounder) QDB (Quadratic Dominated Bounder)
More informationDeterministic Global Optimization for Dynamic Systems Using Interval Analysis
Deterministic Global Optimization for Dynamic Systems Using Interval Analysis Youdong Lin and Mark A. Stadtherr Department of Chemical and Biomolecular Engineering University of Notre Dame, Notre Dame,
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationIs the Hénon map chaotic
Is the Hénon map chaotic Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology, Poland, galias@agh.edu.pl International Workshop on Complex Networks and Applications
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
More informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
More informationLecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.
Lecture 7 3 Ordinary differential equations (ODEs) (continued) 6 Linear equations of second order 7 Systems of differential equations Please note Please note that there is no lecture on Tuesday, 15 November
More informationCHAPTER 10 PERMUTATION ENTROPY, MULTISCALE PERMUTATION ENTROPY AND WAVELET TRANSFORM ANALYSIS OF HEART RATE VARIABILITY
157 CHAPTER 10 PERMUTATION ENTROPY, MULTISCALE PERMUTATION ENTROPY AND WAVELET TRANSFORM ANALYSIS OF HEART RATE VARIABILITY 10.1 INTRODUCTION Complexity parameters for time series are produced based on
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 informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
More informationSolving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method
World Applied Sciences Journal 21 (6): 939-944, 2013 ISSN 11-4952 IDOSI Publications, 2013 DOI: 10.529/idosi.wasj.2013.21.6.2915 Solving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method 1 1 3
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 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 informationTuning PI controllers in non-linear uncertain closed-loop systems with interval analysis
Tuning PI controllers in non-linear uncertain closed-loop systems with interval analysis J. Alexandre dit Sandretto, A. Chapoutot and O. Mullier U2IS, ENSTA ParisTech SYNCOP April 11, 2015 Closed-loop
More informationThe Big, Big Picture (Bifurcations II)
The Big, Big Picture (Bifurcations II) Reading for this lecture: NDAC, Chapter 8 and Sec. 10.0-10.4. 1 Beyond fixed points: Bifurcation: Qualitative change in behavior as a control parameter is (slowly)
More informationTowards a Mechanised Denotational Semantics for Modelica
Towards a Mechanised Denotational Semantics for Modelica Simon Foster Bernhard Thiele Jim Woodcock Peter Fritzson Department of Computer Science, University of York PELAB, Linköping University 3rd February
More informationWHAT IS A CHAOTIC ATTRACTOR?
WHAT IS A CHAOTIC ATTRACTOR? CLARK ROBINSON Abstract. Devaney gave a mathematical definition of the term chaos, which had earlier been introduced by Yorke. We discuss issues involved in choosing the properties
More informationLectures on Periodic Orbits
Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete
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 informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationOn the stability of the three classes of Newtonian three-body planar periodic orbits
. Article. SCIENCE CHINA Physics, Mechanics & Astronomy November 2014 Vol. 57 No. 11: 2121 2126 doi: 10.1007/s11433-014-5563-5 On the stability of the three classes of Newtonian three-body planar periodic
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationTime-adaptive methods for the incompressible Navier-Stokes equations
Time-adaptive methods for the incompressible Navier-Stokes equations Joachim Rang, Thorsten Grahs, Justin Wiegmann, 29.09.2016 Contents Introduction Diagonally implicit Runge Kutta methods Results with
More informationComputational project: Modelling a simple quadrupole mass spectrometer
Computational project: Modelling a simple quadrupole mass spectrometer Martin Duy Tat a, Anders Hagen Jarmund a a Norges Teknisk-Naturvitenskapelige Universitet, Trondheim, Norway. Abstract In this project
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 informationPerturbation theory, KAM theory and Celestial Mechanics 7. KAM theory
Perturbation theory, KAM theory and Celestial Mechanics 7. KAM theory Alessandra Celletti Department of Mathematics University of Roma Tor Vergata Sevilla, 25-27 January 2016 Outline 1. Introduction 2.
More informationConstructing a chaotic system with any number of equilibria
Nonlinear Dyn (2013) 71:429 436 DOI 10.1007/s11071-012-0669-7 ORIGINAL PAPER Constructing a chaotic system with any number of equilibria Xiong Wang Guanrong Chen Received: 9 June 2012 / Accepted: 29 October
More informationPropagation of Uncertainties in Nonlinear Dynamic Models
Propagation of Uncertainties in Nonlinear Dynamic Models Youdong Lin 1, Scott Ferson 2, George F. Corliss 3 and Mark A. Stadtherr 1 1 Department of Chemical and Biomolecular Engineering, University of
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 informationThe Chemical Kinetics Time Step a detailed lecture. Andrew Conley ACOM Division
The Chemical Kinetics Time Step a detailed lecture Andrew Conley ACOM Division Simulation Time Step Deep convection Shallow convection Stratiform tend (sedimentation, detrain, cloud fraction, microphysics)
More informationRound-off Error Analysis of Explicit One-Step Numerical Integration Methods
Round-off Error Analysis of Explicit One-Step Numerical Integration Methods Sylvie Boldo, Florian Faissole, Alexandre Chapoutot To cite this version: Sylvie Boldo, Florian Faissole, Alexandre Chapoutot
More informationWhat we ll do: Lecture 21. Ordinary Differential Equations (ODEs) Differential Equations. Ordinary Differential Equations
What we ll do: Lecture Ordinary Differential Equations J. Chaudhry Department of Mathematics and Statistics University of New Mexico Review ODEs Single Step Methods Euler s method (st order accurate) Runge-Kutta
More informationOn rigorous integration of piece-wise linear systems
On rigorous integration of piece-wise linear systems Zbigniew Galias Department of Electrical Engineering AGH University of Science and Technology 06.2011, Trieste, Italy On rigorous integration of PWL
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationModeling and nonlinear parameter estimation with Kronecker product representation for coupled oscillators and spatiotemporal systems
Physica D 227 (27) 78 99 wwwelseviercom/locate/physd Modeling and nonlinear parameter estimation with Kronecker product representation for coupled oscillators and spatiotemporal systems Chen Yao, Erik
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
More informationMATH 415, WEEK 12 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior
MATH 415, WEEK 1 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior 1 Higher-Dimensional Systems Consider the following system of differential equations: dx = x y dt dy dt = xy y dz dt
More informationH. Koçak - K. Palmer - B. Coomes SHADOWING IN ORDINARY DIFFERENTIAL EQUATIONS
Rend. Sem. Mat. Univ. Pol. Torino Vol. xx, x (xxxx) H. Koçak - K. Palmer - B. Coomes SHADOWING IN ORDINARY DIFFERENTIAL EQUATIONS Abstract. Shadowing deals with the existence of true orbits of dynamical
More informationBifurcation and chaos in simple jerk dynamical systems
PRAMANA c Indian Academy of Sciences Vol. 64, No. 1 journal of January 2005 physics pp. 75 93 Bifurcation and chaos in simple jerk dynamical systems VINOD PATIDAR and K K SUD Department of Physics, College
More informationCertified Roundoff Error Bounds using Semidefinite Programming
Certified Roundoff Error Bounds using Semidefinite Programming Victor Magron, CNRS VERIMAG joint work with G. Constantinides and A. Donaldson INRIA Mescal Team Seminar 19 November 2015 Victor Magron Certified
More informationORBIT 14 The propagator. A. Milani et al. Dipmat, Pisa, Italy
ORBIT 14 The propagator A. Milani et al. Dipmat, Pisa, Italy Orbit 14 Propagator manual 1 Contents 1 Propagator - Interpolator 2 1.1 Propagation and interpolation............................ 2 1.1.1 Introduction..................................
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More informationHopf Bifurcation of a Nonlinear System Derived from Lorenz System Using Centre Manifold Approach ABSTRACT. 1. Introduction
Malaysian Journal of Mathematical Sciences 10(S) March : 1-13 (2016) Special Issue: The 10th IMT-GT International Conference on Mathematics, Statistics and its Applications 2014 (ICMSA 2014) MALAYSIAN
More informationNonlinear Prediction for Top and Bottom Values of Time Series
Vol. 2 No. 1 123 132 (Feb. 2009) 1 1 Nonlinear Prediction for Top and Bottom Values of Time Series Tomoya Suzuki 1 and Masaki Ota 1 To predict time-series data depending on temporal trends, as stock price
More informationTHEODORE VORONOV DIFFERENTIABLE MANIFOLDS. Fall Last updated: November 26, (Under construction.)
4 Vector fields Last updated: November 26, 2009. (Under construction.) 4.1 Tangent vectors as derivations After we have introduced topological notions, we can come back to analysis on manifolds. Let M
More informationRelational Interfaces and Refinement Calculus for Compositional System Reasoning
Relational Interfaces and Refinement Calculus for Compositional System Reasoning Viorel Preoteasa Joint work with Stavros Tripakis and Iulia Dragomir 1 Overview Motivation General refinement Relational
More informationChapter 11 ORDINARY DIFFERENTIAL EQUATIONS
Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown
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 informationCOMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST
CHINESE JOURNAL OF GEOPHYSICS Vol.51, No.3, 2008, pp: 718 724 COMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST DING Rui-Qiang, LI Jian-Ping
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationNumerical Analysis II. Problem Sheet 12
P. Grohs S. Hosseini Ž. Kereta Spring Term 2015 Numerical Analysis II ETH Zürich D-MATH Problem Sheet 12 Problem 12.1 The Exponential Runge-Kutta Method Consider the exponential Runge Kutta single-step
More informationAffine Arithmetic: Concepts and Applications
SCAN 2002 Affine Arithmetic: Concepts and Applications Luiz Henrique de Figueiredo (IMPA) Jorge Stolfi (UNICAMP) Outline What is affine arithmetic? The dependency problem in interval arithmetic Main concepts
More informationACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES
ACCELERATION OF RUNGE-KUTTA INTEGRATION SCHEMES PHAILAUNG PHOHOMSIRI AND FIRDAUS E. UDWADIA Received 24 November 2003 A simple accelerated third-order Runge-Kutta-type, fixed time step, integration scheme
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More information