Analysis of a lumped model of neocortex to study epileptiform ac

 Branden Green
 11 months ago
 Views:
Transcription
1 of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012
2 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized by increased probability of recurring seizures Seizures The activity of braincells hypersynchronizes Observed as large oscillations
3 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized by increased probability of recurring seizures Seizures The activity of braincells hypersynchronizes Observed as large oscillations
4 EEG during a seizure
5 Goal of this research Goals Study several models of neural activity in neocortex Identify correspondences between these models Extrapolate results from simpler models to more complex models Results so far A large, detailed model describing activity of individual neurons in both normal and epileptiform states A simple, lumped model that has shown to have some similarity [Visser et al., 2010]
6 Motivation of model Simulation on Youtube
7 Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
8 Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
9 Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
10 Model Introduction Two excitatory populations: ẋ 1 (t) = µ 1 x 1 (t) F 1 (x 1 (t τ i )) + G 1 (x 2 (t τ e )) ẋ 2 (t) = µ 2 x 2 (t) F 2 (x 2 (t τ i )) + G 2 (x 1 (t τ e ))
11 Simplifications Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Make system symmetric: µ 1 = µ 2 := µ F 1 = F 2 := F G 1 = G 2 := G Choose F(x) and G(x) as rescalings of: S = (tanh(x a) tanh( a)) cosh 2 ( a)
12 Time rescaling Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Time is rescaled to nondimensionalize system: ẋ 1 (t) = x 1 (t) α 1 S(β 1 x 1 (t τ 1 )) + α 2 S(β 2 x 2 (t τ 2 )) ẋ 2 (t) = x 2 (t) α 1 S(β 1 x 2 (t τ 1 )) + α 2 S(β 2 x 1 (t τ 2 )) Introduce vector notation: ẋ(t) = f(x t ), with x t C([ h, 0], R 2 ) and h = max(τ 1, τ 2 ). Main parameters of interest: representing connection strengths between populations.
13 Time rescaling Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Time is rescaled to nondimensionalize system: ẋ 1 (t) = x 1 (t) α 1 S(β 1 x 1 (t τ 1 )) + α 2 S(β 2 x 2 (t τ 2 )) ẋ 2 (t) = x 2 (t) α 1 S(β 1 x 2 (t τ 1 )) + α 2 S(β 2 x 1 (t τ 2 )) Introduce vector notation: ẋ(t) = f(x t ), with x t C([ h, 0], R 2 ) and h = max(τ 1, τ 2 ). Main parameters of interest: representing connection strengths between populations.
14 Fixed points Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Corollary The origin is always a fixed point of the DDE. Theorem The DDE only has only symmetric fixed points, i.e. f(x) = 0 = x = (x, x ).
15 Linearization Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient The DDE is linearized about an equilibrium: in which u 1 (t) = u 1 (t) k 1 u 1 (t τ 1 ) + k 2 u 2 (t τ 2 ), u 2 (t) = u 2 (t) k 1 u 2 (t τ 1 ) + k 2 u 1 (t τ 2 ), k 1 := α 1 β 1 S (β 1 x ) k 2 := α 2 β 2 S (β 2 x ). Note that k 1 and k 2 depend on the equilibrium at which we linearize.
16 Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Characteristic equation and eigenvalues Consider solutions of the form u(t) = e λt c, c R 2. Nontrivial solutions of the system correspond with (λ)c = 0, for: [ λ k1 e (λ) = λτ 1 k 2 e λτ ] 2 k 2 e λτ 2 λ k 1 e λτ 1 Nontrivial solutions exist when det (λ) = 0, or: (λ k 1 e λτ 1 + k 2 e λτ 2 ) (λ k }{{} 1 e λτ 1 k 2 e λτ 2 ) = 0. }{{} := +(λ) := (λ) This characteristic function has an infinite (but countable) number of roots.
17 Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Characteristic equation and eigenvalues Consider solutions of the form u(t) = e λt c, c R 2. Nontrivial solutions of the system correspond with (λ)c = 0, for: [ λ k1 e (λ) = λτ 1 k 2 e λτ ] 2 k 2 e λτ 2 λ k 1 e λτ 1 Nontrivial solutions exist when det (λ) = 0, or: (λ k 1 e λτ 1 + k 2 e λτ 2 ) (λ k }{{} 1 e λτ 1 k 2 e λτ 2 ) = 0. }{{} := +(λ) := (λ) This characteristic function has an infinite (but countable) number of roots.
18 Symmetries of solutions Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Roots of correspond to symmetric solutions, whereas roots of + relate to asymetric solutions. Proof. Let Z 2 act on C 2 so that 1 Z 2 acts as ξ(x, y) : (x, y) (y, x), then: { (λ)v = 0 (λ) = 0 ξv = v { (λ)v = 0 + (λ) = 0 ξv = v
19 Symmetries of solutions Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Roots of correspond to symmetric solutions, whereas roots of + relate to asymetric solutions. Proof. Let Z 2 act on C 2 so that 1 Z 2 acts as ξ(x, y) : (x, y) (y, x), then: { (λ)v = 0 (λ) = 0 ξv = v { (λ)v = 0 + (λ) = 0 ξv = v
20 Roots of (λ) Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For λ C 0 e λh and λ > C, with C C 0, (λ) > 0. Hence, (λ) has a zero free right halfplane [Bellman and Cooke, 1963].
21 Minimal Stability Region Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For k 1 + k 2 < 1 the system has no eigenvalues with positive real part.
22 Bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Note: As in ODEs, an equilibrium of a DDE loses its stablity when either a single real eigenvalue passes through zero (fold/transcritical), or a pair of complex eigenvalues passes through the imaginary axis (Hopf).
23 Fold/transcritical bifurcation Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Equilibria undergo a fold or transcritical bifurcation on the lines 1 + k 1 ± k 2 = 0 in (k 1, k 2 )space.
24 Hopf bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem Two curves, h + (ω) and h (ω), exist along which the linearized system has a pair of complex eigenvalues λ = ±iω. Proof. For now, we consider only + (iω) = 0; (iω) is similar. iω k 1 e iωτ 1 + k 2 e iωτ 2 = 0 [ ] [ ] [ ] cos(ωτ1 ) cos(ωτ 2 ) k1 1 = sin(ωτ 1 ) sin(ωτ 2 ) ω k 2
25 Hopf bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Proof cntd. The matrix is invertible if det = sin(ω(τ 2 τ 1 )) 0, yielding: [ ] [ ] [ ] k1 1 sin(ωτ2 ) cos(ωτ = h + (ω) := 2 ) 1. sin(ω(τ 2 τ 1 )) sin(ωτ 1 ) cos(ωτ 1 ) ω k 2
26 Bifurcations in (k 1, k 2 )space Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Very sensitive, complex branch structure Intersections correspond with codim2 bifurcations. k 2 k 1 (τ 1 = 11.6, τ 2 = 20.3)
27 Returning branches Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Theorem For a given ω 0 > 0, (λ) has no roots λ = ±iω, ω > ω 0 inside the square k 1 + k 2 < 1 + ω0 2. Proof. First, assume iω is a root of + (a similar argument holds for ): 0 = 1 + iω + k 1 e iωτ 1 + k 2 e iωτ 2, 1 + iω k 1 k 2. Then, this root lies outside the designated square for ω > ω 0.
28 Revised stability region Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Conjecture For the considered parameters τ 1 and τ 2 the full stability region is bounded, connected and contains the origin.
29 The first Lyapunov coefficient Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Stability region partly bounded by curve of Hopfs The first Lyapunov coefficient is determined along the boundary Normal given by [Diekmann et al, 1995]: c 1 = 1 2 qt D 3 f(0)(φ, φ, φ) + q T D 2 f(0)(e 0 (0) 1 D 2 f(0)(φ, φ), φ) qt D 2 f(0)(e 2iω (2iω) 1 D 2 f(0)(φ, φ), φ).
30 The first Lyapunov coefficient Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Stability region partly bounded by curve of Hopfs The first Lyapunov coefficient is determined along the boundary Normal given by [Diekmann et al, 1995]: c 1 = 1 2 qt D 3 f(0)(φ, φ, φ) + q T D 2 f(0)(e 0 (0) 1 D 2 f(0)(φ, φ), φ) qt D 2 f(0)(e 2iω (2iω) 1 D 2 f(0)(φ, φ), φ).
31 Codim 2 bifurcations Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Generalized Hopf Lypunov coefficient passes through zero, Generalized Hopf bifurcation at boundary
32 Summary Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Equilibria Only symmetric equilibria Stability region identified Bifurcations of stability region Fold and transcritical bifurcations Hopfbifurcations, both sub and supercritical ZeroHopf and HopfHopf
33 Summary Introduction Fixed points and eigenvalues Linear stability Bifurcations First Lyapunov coefficient Equilibria Only symmetric equilibria Stability region identified Bifurcations of stability region Fold and transcritical bifurcations Hopfbifurcations, both sub and supercritical ZeroHopf and HopfHopf
34 Numerical bifurcation analysis Overview One parameter Two parameters One parameter Identify stable and unstable manifolds of nontrivial equilibrium and periodic solutions for varying α 2. Two parameters Continue boundaries of stable solutions in α 1 and α 2. Software DDEBIFTOOL [Engelborghs et al., 2002] PDDECONT/Knut [Roose and Szalai, 2007]
35 Numerical bifurcation analysis Overview One parameter Two parameters One parameter Identify stable and unstable manifolds of nontrivial equilibrium and periodic solutions for varying α 2. Two parameters Continue boundaries of stable solutions in α 1 and α 2. Software DDEBIFTOOL [Engelborghs et al., 2002] PDDECONT/Knut [Roose and Szalai, 2007]
36 Overview One parameter Two parameters Bifurcations in α A 2 B C D
37 Overview One parameter Two parameters Bifurcations in α 2 3 Bifurcation diagram II 2.5 I Maximal value 2 III 1.5 IV PD1 0 H1 B α 2 Detail I Detail II LPC4
38 Overview One parameter Two parameters Bifurcations in α 2 l k h Maximal value j g c i d f e α 2 a b
39 Fixed points in (k 1, k 2 )space Overview One parameter Two parameters F1 B1 k H2 H3 H4 H1 Nontrivial Trivial k 1
40 Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α GH ZH2 ZH1 HH
41 Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α CP1 0.3 GH1 0.2 FF ZH2 ZH1 HH
42 Overview One parameter Two parameters Bifurcations of stable solutions in α 1 and α CP1 0.3 GH1 0.2 FF1 CP3 0.1 R1_1 0 ZH2 ZH1 HH
43 Complex bifurcation structure Overview One parameter Two parameters α 1 R2_1 FF2 R2_2 FF3 CP4 CP2 α 2
44 Comparison Introduction The bifurcation analysis is compared with a survey on the detailed model:
45 Introduction Comparison The bifurcation analysis is compared with a survey on the detailed model:
46 Summary and conclusions Summary We studied a simple, lumped model for neural activity Parameter regions of steady states and periodic solutions are identified The behavior of both models varies similarly for changes of parameters Regions of multistability are of interest for studying epilepsy
47 Future work Introduction Details Study the model s codim2 bifurcations Proofs rather than conjectures Next steps Expand model (e.g. break symmetry) Parameter estimation
Description of the extension ddebiftool_nmfm
Description of the extension ddebiftool_nmfm Maikel Bosschaert, Bram Wage and Yuri Kuznetsov 16 September 2015 This addendum describes the normal form extension to DDEBifTool, a bifurcation analysis toolbox
More informationHOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS
PERIODICA POLYTECHNICA SER. MECH. ENG. VOL. 48, NO. 2, PP. 89 200 2004 HOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS Gábor OROSZ Bristol Centre for Applied Nonlinear Mathematics Department of Engineering
More informationActivity types in a neural mass model
Master thesis Activity types in a neural mass model Jurgen Hebbink September 19, 214 Exam committee: Prof. Dr. S.A. van Gils (UT) Dr. H.G.E. Meijer (UT) Dr. G.J.M. Huiskamp (UMC Utrecht) Contents 1 Introduction
More informationas Hopf Bifurcations in TimeDelay Systems with Bandlimited Feedback
as Hopf Bifurcations in TimeDelay Systems with Bandlimited Feedback Lucas Illing and Daniel J. Gauthier Department of Physics Center for Nonlinear and Complex Systems Duke University, North Carolina
More informationTime Delays in Neural Systems
Time Delays in Neural Systems Sue Ann Campbell 1 Department of Applied Mathematics, University of Waterloo, Waterloo ON N2l 3G1 Canada sacampbell@uwaterloo.ca Centre for Nonlinear Dynamics in Physiology
More informationOne Dimensional Dynamical Systems
16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with onedimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar
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 informationChapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point SaddleNode Bifurcation Diagram Codimension1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
More informationReliably computing all characteristic roots of delay differential equations in a given right half plane using a spectral method
Reliably computing all characteristic roots of delay differential equations in a given right half plane using a spectral method Zhen Wu Wim Michiels Report TW 596, May 2011 Katholieke Universiteit Leuven
More informationStability Analysis of UzawaLucas Endogenous Growth Model
Abstract: Stability Analysis of UzawaLucas Endogenous Growth Model William A. Barnett* University of Kansas, Lawrence and Center for Financial Stability, NY City and Taniya Ghosh Indira Gandhi Institute
More informationWeierstraßInstitut. für Angewandte Analysis und Stochastik. LeibnizInstitut im Forschungsverbund Berlin e. V. Preprint ISSN
WeierstraßInstitut für Angewandte Analysis und Stochastik LeibnizInstitut im Forschungsverbund Berlin e. V. Preprint ISSN 0946 8633 On the stability of periodic orbits in delay equations with large delay
More information2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ
1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.
More informationDynamic Feedback and the Design of Closedloop Drug Delivery Systems
Dynamic Feedback and the Design of Closedloop Drug Delivery Systems John Milton 1,2, Sue Ann Campbell 2,3,4 and Jacques Bélair 2,4,5 1) Department of Neurology, The University of Chicago Hospitals, MC
More informationTHREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS
THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the
More informationHomoclinic saddle to saddlefocus transitions in 4D systems
Faculty of Electrical Engineering, Mathematics & Computer Science Homoclinic saddle to saddlefocus transitions in 4D systems Manu Kalia M.Sc. Thesis July 2017 Assessment committee: Prof. Dr. S. A. van
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 informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationBalancing of the skateboard with reflex delay
ENOC 214, July 611, 214, Vienna, Austria Balancing of the skateboard with reflex delay Balazs Varszegi, Denes Takacs, Gabor Stepan and S. John Hogan Department of Applied Mechanics, Budapest University
More informationMaster Thesis. Equivariant Pyragas control. Isabelle Schneider Freie Universität Berlin Fachbereich Mathematik und Informatik Berlin
Master Thesis Equivariant Pyragas control Isabelle Schneider Freie Universität Berlin Fachbereich Mathematik und Informatik 14195 Berlin February 3, 2014 Abstract In this thesis we show how the concept
More informationDevelopment of novel deep brain stimulation techniques: Coordinated reset and Nonlinear delayed feedback
Development of novel deep brain stimulation techniques: Coordinated reset and Nonlinear delayed feedback Oleksandr Popovych 1, Christian Hauptmann 1, Peter A. Tass 1,2 1 Institute of Medicine and Virtual
More informationHopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 65 71 Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation Faridon AMDJADI Department
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, LinearQuadratic optimal control, Riccati equations (differential, algebraic, discretetime), controllability,
More informationDELAY DIFFERENTIAL EQUATIONS AND CONTINUATION
DELAY DIFFERENTIAL EQUATIONS AND CONTINUATION JEANPHILIPPE LESSARD Abstract. In these lecture notes, we demonstrate how rigorous numerics can help studying the dynamics of delay equations. We present
More informationDef. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)
Types of critical points Def. (a, b) is a critical point of the autonomous system Math 216 Differential Equations Kenneth Harris kaharri@umich.edu Department of Mathematics University of Michigan November
More information8.1 Bifurcations of Equilibria
1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integrodifferential equations
More informationHopf Bifurcation in a Scalar Reaction Diffusion Equation
journal of differential equations 140, 209222 (1997) article no. DE973307 Hopf Bifurcation in a Scalar Reaction Diffusion Equation Patrick Guidotti and Sandro Merino Mathematisches Institut, Universita
More informationUniversity of Utrecht
University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:
More informationBifurcations and multiple traffic jams in a carfollowing model with reactiontime delay
Physica D 211 (2005) 277 293 Bifurcations and multiple traffic jams in a carfollowing model with reactiontime delay Gábor Orosz, Bernd Krauskopf, R. Eddie Wilson Bristol Centre for Applied Nonlinear
More informationThree ways of treating a linear delay differential equation
Proceedings of the 5th International Conference on Nonlinear Dynamics NDKhPI2016 September 2730, 2016, Kharkov, Ukraine Three ways of treating a linear delay differential equation Si Mohamed Sah 1 *,
More informationLecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. PoincareBendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion PoincareBendixson theorem Hopf bifurcations Poincare maps Limit
More informationDamped harmonic motion
Damped harmonic motion March 3, 016 Harmonic motion is studied in the presence of a damping force proportional to the velocity. The complex method is introduced, and the different cases of underdamping,
More informationDynamics of a Networked Connectivity Model of Waterborne Disease Epidemics
Dynamics of a Networked Connectivity Model of Waterborne Disease Epidemics A. Edwards, D. Mercadante, C. Retamoza REU Final Presentation July 31, 214 Overview Background on waterborne diseases Introduction
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationASYMPTOTIC BEHAVIOUR AND HOPF BIFURCATION OF A THREEDIMENSIONAL NONLINEAR AUTONOMOUS SYSTEM
Georgian Mathematical Journal Volume 9 (), Number, 7 6 ASYMPTOTIC BEHAVIOUR AND HOPF BIFURCATION OF A THREEDIMENSIONAL NONLINEAR AUTONOMOUS SYSTEM LENKA BARÁKOVÁ Abstract. A threedimensional real nonlinear
More informationStability of Limit Cycle in a Delayed ISLM Business Cycle model
Applied Mathematical Sciences, Vol. 2, 28, no. 5, 2459247 Stability of Limit Cycle in a Delayed ISLM Business Cycle model A. Abta, A. Kaddar H. Talibi Alaoui Université Chouaib Doukkali Faculté des Sciences
More informationMath 312 Lecture Notes Linear Twodimensional Systems of Differential Equations
Math 2 Lecture Notes Linear Twodimensional Systems of Differential Equations Warren Weckesser Department of Mathematics Colgate University February 2005 In these notes, we consider the linear system of
More informationLinearized geometric dynamics of TobinBenhabibMiyao economic flow
Linearized geometric dynamics of TobinBenhabibMiyao economic flow Constantin Udrişte and Armando Ciancio Abstract The aim of this paper is to study the influence of the EuclideanLagrangian structure
More informationNonlinear Dynamics of Neural Firing
Nonlinear Dynamics of Neural Firing BENG/BGGN 260 Neurodynamics University of California, San Diego Week 3 BENG/BGGN 260 Neurodynamics (UCSD) Nonlinear Dynamics of Neural Firing Week 3 1 / 16 Reading Materials
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationIntroduction to the Phase Plane
Introduction to the Phase Plane June, 06 The Phase Line A single first order differential equation of the form = f(y) () makes no mention of t in the function f. Such a differential equation is called
More 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 informationHOPF BIFURCATION WITH S 3 SYMMETRY
PORTUGALIAE MATHEMATICA Vol. 63 Fasc. 006 Nova Série HOPF BIFURCATION WITH S 3 SYMMETRY Ana Paula S. Dias and Rui C. Paiva Recommended by José BastoGonçalves Abstract: The aim of this paper is to study
More informationMath 273 (51)  Final
Name: Id #: Math 273 (5)  Final Autumn Quarter 26 Thursday, December 8, 266: to 8: Instructions: Prob. Points Score possible 25 2 25 3 25 TOTAL 75 Read each problem carefully. Write legibly. Show all
More informationThe Hopfvan der Pol System: Failure of a Homotopy Method
DOI.7/s25995 ORIGINAL RESEARCH The Hopfvan der Pol System: Failure of a Homotopy Method H. G. E. Meijer T. KalmárNagy Foundation for Scientific Research and Technological Innovation 2 Abstract The
More informationHOPF BIFURCATIONS AND LIMIT CYCLES IN EVOLUTIONARY NETWORK DYNAMICS
HOPF BIFURCATIOS AD LIMIT CYCLES I EVOLUTIOARY ETWORK DYAMICS DARRE PAIS, CARLOS H. CAICEDOÚÑEZ, AD AOMI E. LEOARD Abstract. The replicatormutator equations from evolutionary dynamics serve as a model
More informationNonlinear Control Lecture 4: Stability Analysis I
Nonlinear Control Lecture 4: Stability Analysis I Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 4 1/70
More informationBEYOND INPHASE AND ANTIPHASE COORDINATION IN A MODEL OF JOINT ACTION
BEYOND INPHASE AND ANTIPHASE COORDINATION IN A MODEL OF JOINT ACTION DANIELE AVITABILE, PIOTR S lowiński, BENOIT BARDY, AND KRASIMIRA TSANEVAATANASOVA Abstract. In 1985 Haken, Kelso and Bunz proposed
More informationPHY411 Lecture notes Part 4
PHY411 Lecture notes Part 4 Alice Quillen February 1, 2016 Contents 0.1 Introduction.................................... 2 1 Bifurcations of onedimensional dynamical systems 2 1.1 Saddlenode bifurcation.............................
More informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 SYMMETRY P. Yu 1,2
More informationDiagonalization of Matrices
LECTURE 4 Diagonalization of Matrices Recall that a diagonal matrix is a square n n matrix with nonzero entries only along the diagonal from the upper left to the lower right (the main diagonal) Diagonal
More informationLocal stability and Hopf bifurcation analysis. for Compound TCP
Local stability and Hopf bifurcation analysis for Compound TCP Debayani Ghosh, Krishna Jagannathan and Gaurav Raina arxiv:604.05549v [math.ds] 9 Apr 06 Abstract We conduct a local stability and Hopf bifurcation
More informationPropagating terraces and the dynamics of frontlike solutions of reactiondiffusion equations on R
Propagating terraces and the dynamics of frontlike solutions of reactiondiffusion equations on R P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455 Abstract We consider semilinear
More informationDerivation of bordercollision maps from limit cycle bifurcations
Derivation of bordercollision maps from limit cycle bifurcations Alan Champneys Department of Engineering Mathematics, University of Bristol Mario di Bernardo, Chris Budd, Piotr Kowalczyk Gabor Licsko,...
More informationME 680 Spring Geometrical Analysis of 1D Dynamical Systems
ME 680 Spring 2014 Geometrical Analysis of 1D Dynamical Systems 1 Geometrical Analysis of 1D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
More informationTIMEDELAY IN SPRINGBLOCK MODEL FOR APERIODICITY IN EARTHQUAKES
TIMEDELAY IN SPRINGBLOCK MODEL FOR APERIODICITY IN EARTHQUAKES S Kostic & N Vasovic University of Belgrade Faculty of Mining and Geology, Serbia I Franovic University of Belgrade Faculty of Physics,
More informationMathematical Modeling in the Kidney. Characterizing a LongLooped Nephron
: Characterizing a LongLooped Nephron Quinton Neville Department of MSCS, St. Olaf College October 13, 2016 Northfield Undergraduate Math Symposium Kidney and Nephron Tubuloglomerular Feedback (TGF) System
More informationComplex Systems Workshop Lecture III: Behavioral Asset Pricing Model with Heterogeneous Beliefs
Complex Systems Workshop Lecture III: Behavioral Asset Pricing Model with Heterogeneous Beliefs Cars Hommes CeNDEF, UvA CEF 2013, July 9, Vancouver Cars Hommes (CeNDEF, UvA) Complex Systems CEF 2013, Vancouver
More informationSymmetry, Hopf Bifurcation and the Emergence of Cluster Solutions in Time Delayed Neural Networks
Symmetry, Hopf Bifurcation and the Emergence of Cluster Solutions in Time Delayed Neural Networks Zhen Wang 1 and Sue Ann Campbell a) Department of Applied Mathematics and Centre for Theoretical Neuroscience,
More informationDynamical system theory and bifurcation analysis for microscopic traffic models
Dynamical system theory and bifurcation analysis for microscopic traffic models Bodo Werner University of Hamburg L Aquila April 2010 2010 MathMods IP mailto:werner@math.unihamburg.de April 30, 2010 Abstract
More informationKatholieke Universiteit Leuven Department of Computer Science
Pseudospectra and stability radii for analytic matrix functions with application to timedelay systems Wim Michiels Kirk Green Thomas Wagenknecht SilviuIulian Niculescu Report TW 425, March 2005 Katholieke
More informationBifurcation, stability, and cluster formation of multistrain infection models
J. Math. Biol. (23) 67:57 532 DOI.7/s285263 Mathematical Biology Bifurcation, stability, and cluster formation of multistrain infection models Bernard S. Chan Pei Yu Received: January 22 / Revised:
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 informationLogistic Map, Euler & RungeKutta Method and LotkaVolterra Equations
Logistic Map, Euler & RungeKutta Method and LotkaVolterra 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 informationarxiv: v1 [nlin.ao] 19 May 2017
Featurerich bifurcations in a simple electronic circuit Debdipta Goswami 1, and Subhankar Ray 2, 1 Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA
More information2:1 Resonance in the delayed nonlinear Mathieu equation
Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s110710069165 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September
More informationMultiscale Analysis of Many Particle Systems with Dynamical Control
Michael Herrmann Multiscale Analysis of Many Particle Systems with Dynamical Control joint work with Barbara Niethammer and Juan J.L. Velázquez Kinetic description of multiscale phenomena Archimedes Center
More informationBIBO STABILITY AND ASYMPTOTIC STABILITY
BIBO STABILITY AND ASYMPTOTIC STABILITY FRANCESCO NORI Abstract. In this report with discuss the concepts of boundedinput boundedoutput stability (BIBO) and of Lyapunov stability. Examples are given to
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a nonlinear timeinvariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationEECS C128/ ME C134 Final Thu. May 14, pm. Closed book. One page, 2 sides of formula sheets. No calculators.
Name: SID: EECS C28/ ME C34 Final Thu. May 4, 25 58 pm Closed book. One page, 2 sides of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 4 2 4 3 6 4 8 5 3
More informationSimple Chaotic Oscillator: From Mathematical Model to Practical Experiment
6 J. PERŽELA, Z. KOLKA, S. HANUS, SIMPLE CHAOIC OSCILLAOR: FROM MAHEMAICAL MODEL Simple Chaotic Oscillator: From Mathematical Model to Practical Experiment Jiří PERŽELA, Zdeněk KOLKA, Stanislav HANUS Dept.
More informationDynamical Systems Generated by ODEs and Maps: Final Examination Project
Dynamical Systems Generated by ODEs and Maps: Final Examination Project Jasmine Nirody 26 January, 2010 Contents 1 Introduction and Background 1 2 Equilibria 1 3 Attracting Domain 2 4 Overview of Behavior
More information2 Qualitative theory of nonsmooth dynamical systems
2 Qualitative theory of nonsmooth dynamical systems In this chapter, we give an overview of the basic theory of both smooth and nonsmooth dynamical systems, to be expanded upon in later chapters. In
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J. APPL. MATH. Vol. 60, No. 2, pp. 503 535 c 2000 Society for Industrial and Applied Mathematics SUBCRITICAL ELLIPTIC BURSTING OF BAUTIN TYPE EUGENE M. IZHIKEVICH Abstract. Bursting behavior in neurons
More informationChapter 1 Bifurcations and Chaos in Dynamical Systems
Chapter 1 Bifurcations and Chaos in Dynamical Systems Complex system theory deals with dynamical systems containing often a large number of variables. It extends dynamical system theory, which deals with
More informationA ProducerConsumer Model With Stoichiometry
A ProducerConsumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationInvariant Manifolds of Dynamical Systems and an application to Space Exploration
Invariant Manifolds of Dynamical Systems and an application to Space Exploration Mateo Wirth January 13, 2014 1 Abstract In this paper we go over the basics of stable and unstable manifolds associated
More informationTwo Dimensional Linear Systems of ODEs
34 CHAPTER 3 Two Dimensional Linear Sstems of ODEs A firstder, autonomous, homogeneous linear sstem of two ODEs has the fm x t ax + b, t cx + d where a, b, c, d are real constants The matrix fm is 31
More informationInteraction of Canard and Singular Hopf Mechanisms in a Neural Model
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 0, No. 4, pp. 443 479 c 0 Society for Industrial and Applied Mathematics Interaction of Canard and Singular Hopf Mechanisms in a Neural Model R. Curtu and J. Rubin
More informationFirst Midterm Exam Name: Practice Problems September 19, x = ax + sin x.
Math 54 Treibergs First Midterm Exam Name: Practice Problems September 9, 24 Consider the family of differential equations for the parameter a: (a Sketch the phase line when a x ax + sin x (b Use the graphs
More informationJournal of Applied Nonlinear Dynamics
Journal of Applied Nonlinear Dynamics 1(2) (2012) 159 171 Journal of Applied Nonlinear Dynamics https://lhscientificpublishing.com/journals/janddefault.aspx Asymptotic Analysis of the HopfHopf Bifurcation
More informationAutonomous system = system without inputs
Autonomous system = system without inputs State space representation B(A,C) = {y there is x, such that σx = Ax, y = Cx } x is the state, n := dim(x) is the state dimension, y is the output Polynomial representation
More informationNumerical Computation of Canards
Numerical Computation of Canards John Guckenheimer Kathleen Hoffman Warren Weckesser January 9, 23 Shortened title: Computing Canards Abstract Singularly perturbed systems of ordinary differential equations
More informationELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky
ELEC 3035, Lecture 3: Autonomous systems Ivan Markovsky Equilibrium points and linearization Eigenvalue decomposition and modal form State transition matrix and matrix exponential Stability ELEC 3035 (Part
More informationInternet Congestion Control: Equilibrium and Dynamics
Internet Congestion Control: Equilibrium and Dynamics A. Kevin Tang Cornell University ISS Seminar, Princeton University, February 21, 2008 Networks and Corresponding Theories Power networks (Maxwell Theory)
More informationChaos in the HodgkinHuxley Equations: The TakensBodganov Cusp Bifurcation. Rami Mohieddine
CORNELL UNIVERSITY MATHEMATICS DEPARTMENT SENIOR THESIS Chaos in the HodgkinHuxley Equations: The TakensBodganov Cusp Bifurcation A THESIS PRESENTED IN PARTIAL FULFILLMENT OF CRITERIA FOR HONORS IN MATHEMATICS
More information5. THE CLASSES OF FOURIER TRANSFORMS
5. THE CLASSES OF FOURIER TRANSFORMS There are four classes of Fourier transform, which are represented in the following table. So far, we have concentrated on the discrete Fourier transform. Table 1.
More informationNonlinear differential equations  phase plane analysis
Nonlinear differential equations  phase plane analysis We consider the general first order differential equation for y(x Revision Q(x, y f(x, y dx P (x, y. ( Curves in the (x, yplane which satisfy this
More informationGlobal Bifurcation Diagram for the KernerKonhäuser Traffic Flow Model
arxiv:1311.4119v1 [math.ds] 17 Nov 2013 Global Bifurcation Diagram for the KernerKonhäuser Traffic Flow Model Joaquín Delgado Patricia Saavedra Keywords: Continuous traffic flow. Traveling waves. Bautin
More informationOptimal Control. Lecture 18. HamiltonJacobiBellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 HamiltonJacobiBellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline HamiltonJacobiBellman (HJB) Equation Iterative solution of HJB Equation
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 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 ndimensional C k differentiable manifold
More informationReduced order observer design for nonlinear systems
Applied Mathematics Letters 19 (2006) 936 941 www.elsevier.com/locate/aml Reduced order observer design for nonlinear systems V. Sundarapandian Department of Instrumentation and Control Engineering, SRM
More informationGLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationControlling biological networks by timedelayed signals
Submitted, Phil. Trans. of the Royal Society  A 29) http://www.cds.caltech.edu/~murray/papers/omm9ptrsa.html Controlling biological networks by timedelayed signals By Gábor Orosz 1,, Jeff Moehlis 1
More informationON THE PROBLEM OF LINEARIZATION FOR STATEDEPENDENT DELAY DIFFERENTIAL EQUATIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 5, May 1996 ON THE PROBLEM OF LINEARIZATION FOR STATEDEPENDENT DELAY DIFFERENTIAL EQUATIONS KENNETH L. COOKE AND WENZHANG HUANG (Communicated
More informationAMPLIFIED HOPF BIFURCATIONS IN FEEDFORWARD NETWORKS
AMPLIFIED HOPF BIFURCATIONS IN FEEDFORWARD NETWORKS BOB RINK AND JAN SANDERS Abstract In [18] the authors developed a method for computing normal forms of dynamical systems with a coupled cell network
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 information18.06 Problem Set 8  Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8  Solutions Due Wednesday, 4 November 2007 at 4 pm in 206 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More informationProblem Set Number 1, j/2.036j MIT (Fall 2014)
Problem Set Number 1, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept., Cambridge, MA 02139) September 9, 2014 Due Fri., September 19, 2014. 1 Get equation from phase line portrait problem
More information