8.1 Bifurcations of Equilibria


 Bernard O’Neal’
 1 years ago
 Views:
Transcription
1 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 discrete mappings etc Of course we are concerned with ODEs Local bifurcations refer to qualitative changes occurring in a neighborhood of an equilibrium point of a differential equation or a fixed point of an associated Poincaré map These can be studied by expanding the equations of motion in power series about the point Consider an ODE depending on a parameter : Assume this system has an equilibrium point that is a sink for where is a critical value of the parameter Graph the evls as a function of Real and imaginary axes region of evls with negative real parts for arrows showing that the real parts of some evls are increasing with A real evl or a cc pair may traverse the imaginary axis as increases through The sink becomes a source or a saddle Or the equilibrium point may simply disappear when it has a zero evl! Example Fishery model with constant harvesting Recall the logistic model of population growth with an additional constant term Interpret as a model of a fishery with the number of fish time the fish growth rate the carrying capacity and a constant rate of harvesting The model may be simplified by measuring the quantity of fish and the harvest rate in units of the carrying capacity: and Further simplify by introducing a dimensionless time variable: Then obtain [1] Equilibria correspond to Analyze graphically
2 2 sketch   The direction of the arrows on the axis now depends on the sign of When there are two equilibrium points when there is one and when there are none Alternately we can study the equilibria algebraically Let Then equilibrium points of [1] correspond to or The equilibria are These exist only for add arrows to graph We can write [1] in the form [2] where the dot now refers to the derivative with respect to dimensionless time The eigenvalues of the equilibrium points are which depends on through Since its eigenvalue is negative and is a stable equilibrium point Since its eigenvalue is positive and is an unstable equilibrium point As increases toward the equilibria converge and their eigenvalues approach 0 At the equilibria merge as a single nonhyperbolic equilibrium point and for they disappear Draw a bifurcation diagram for the logistic model with constant harvesting by plotting the equilibria as a function of sketch   the two branches of the saddlenode A saddlenode bifurcation occurs at a critical value of the parameter Simplify the logistic model with constant harvesting by centering the bifurcation at the origin Let and Then [2] becomes [3]
3 3 This is the normal form of the saddlenode bifurcation For there are two hyperbolic equilibrium points for there is a single nonhyperbolic equilibrium and for there are no equilibria The saddlenode bifurcation requires 3 conditions on the vector field: Singularity condition (the equilibrium point [3] is nonhyperbolic at ) Nondegeneracy condition (the coefficient of in [3] is nonzero) Transversality condition (that guarantees that the parameter perturbs the nonhyperbolic equilibrium point in a transverse way: in [3]) Local Bifurcations of Limit Cycles To further illustrate the meaning of local bifurcation let s briefly describe local bifurcations of limit cycles sketch limit cycle section intersection Consider a Poincaré map corresponds to a fixed point of : The eigenvalues of are the Floquet multipliers (not including the trivial unit multiplier) Real and imaginary axes region of evls within the unit circle for arrows showing that evls can traverse the circle though or a complex conjugate pair can pass out of the unit circle off the real axis Example A saddlenode bifurcation of limit cycles Suppose that for there is a semistable limit cycle and an associated Poincaré map: sketch semistable limit cycle in 2D and nearby trajectories sketch   the line the Poincaré map is concave down and tangent to the line at the origin cobweb paths inside and outside the limit cycle sketch pair of limit cycles The outer one is stable and the inner one unstable Nearby trajectories sketch   the line the Poincaré map has shifted so there are now two intersections with the line Three cobweb paths
4 4 sketch spiral flow without limit cycle no intersection between the Poincare map with the line The saddle node bifurcation corresponds to a single eigenvalue reaching the unit circle at and then disappearing Global Bifurcations of Limit Cycles Global bifurcations cannot be described by a local analysis Example Homoclinic bifurcation in sketch like Meiss Fig 818 and without the point and axis Notice the limit cycle for 82 Preservation of Equilibria Consider An equilibrium point satisfies When we linearize about the equilibrium point the Jacobian matrix is If the Jacobian matrix has a zero eigenvalue then it is singular and we say the equilibrium point is degenerate Otherwise the Jacobian matrix is nonsingular and the equilibrium point is nondegenerate Contrast this with the definition of a nonhyperbolic equilibrium point where only the real part of an eigenvalue need be zero For example if a Jacobian matrix has a pair of complex conjugate imaginary eigenvalues and all of the other eigenvalues are nonzero it is not singular Example Fishery model with constant harvesting One way that we have written the equation of motion is Let be an equilibrium point When and the equilibrium point is degenerate; corresponds to the saddlenode bifurcation For the equilibrium point disappears
5 5 If an equilibrium point is nondegenerate it cannot be removed by sufficiently small perturbations (such as a small change in the value of Theorem 81 (Implicit Function) Let be an open set in and with Suppose there is a point such that and is a nonsingular matrix Then there are open sets and and a unique function for which and The proof of this famous theorem probably appears in your favorite analysis book To gain a rough understanding of why the condition on the Jacobian is necessary expand about : Neglect the higherorder terms and solve for to obtain This can be done for arbitrary only if is nonsingular When discussing the system the application of the implicit function theorem to preservation of equilibria corresponds to and Corollary 82 (Preservation of a Nondegenerate Equilibrium) Suppose the vector field is in both and and that is a nondegenerate equilibrium point for parameter value Then there exists a unique curve of equilibria passing through at Proof The matrix governs the stability of the equilibrium and since has all of its eigenvalues nonzero is nonsingular Then theorem 81 implies that there is a neighborhood of for which there is a curve of equilibria 83 Unfolding Vector Fields Change of Variables and Topological Conjugacy Recall topological conjugacy between flows corresponds to the diagram:
6 6 where is a homeomorphism The diagram implies Bifurcation theory studies systems that depend on parameters Let be the flow of and let be the flow of If these flows are conjugate for each value of then there is the diagram: where is a homeomorphism The diagram implies Example Consider the system Change variables using to obtain Note that is a homeomorphism (in fact a diffeomorphism) The corresponding flows are topologically conjugate Meiss extends this terminology to the vector fields; he also calls and conjugate For a fixed value of gives a vector field For all possible values of gives a family of vector fields A family of vector fields is induced by a family if there is a continuous map such that If the family is induced by the family then has the same dynamics as or is simpler than Example The family is induced by the family using the map The family has the same dynamics as (In fact for corresponding values of and the vector fields are the same) Example Consider the family and the constant map The family of vector fields induced by using the map is The family is clearly simpler than the family Let be the flow of and be the flow of where the vector fields and are conjugate for corresponding values of and Then we have the diagram
7 7 where and Then the flows are related by [1] If is a diffeomorphism we can also find a relationship between the corresponding vector fields Differentiate [1] with respect to : Set to obtain the desired relationship [2] where and Compare with Eq 434 in Meiss Example Let and where and Then and [2] gives This can be verified by making the indicated substitutions for Unfolding Vector Fields and Consider a vector field that has a degenerate orbit This could be a degenerate equilibrium point or for example the homoclinic orbit in the homoclinic bifurcation Then we say fulfills a singularity condition A family of vector fields is an unfolding of if Example The singularity associated with the saddlenode bifurcation is is degenerate because it has a zero eigenvalue Two unfolding of are and
8 8 An unfolding of is versal if it contains all possible qualitative dynamics that can occur near to This means that every other unfolding in some neighborhood of will have the same dynamics as some family induced by Example Let and We will show below that is a versal unfolding of the saddle node bifurcation However is not versal For example there is no value of for which the system has two equilibrium points An unfolding is miniversal if it is a versal unfolding with the minimum number of parameters Example Let and The flows of and are diffeomorphic To see this complete the square in : The change of variables and converts into Note that is a diffeomorphism Since may be any real number the families of vector fields and have the same dynamics However has two parameters and has only one Then is a versal unfolding but not a miniversal unfolding of the saddle node bifurcation 84 SaddleNode Bifurcation in One Dimension Consider the ODE on Theorem 83 Suppose that with a nonhyperbolic equilibrium at the origin and that satisfies the nondegeneracy condition Then there is a such that when there is an open interval containing such that there is a unique extremal value There are two equilibria in when one when and zero when Proof Let The singularity and nondegeneracy conditions imply
9 9 Since the higher order term satisfies A general unfolding of will have form [3] where and Consider the function Since the implicit function theorem guarantees neighborhoods and of the origin such that when there is a unique such that and Since there must also be an open interval containing for which is monotone On is either concave up (when ) or concave down (when ) Therefore is the unique extremal value for when and determines whether the critical point is a minimum or a maximum Since when this remains true by continuity for small enough ; for Therefore when has a minimum at and is positive on This implies that if there are no zeros of if there is a single zero and if there are two zeros Similar considerations apply when Picture for :   For is a stable equilibrium and is an unstable equilibrium The bifurcation set is the set of all points in parameter space at which the bifurcation takes place A bifurcation is codimension if the bifurcation set is determined by independent conditions on parameters According to theorem 83 the saddlenode bifurcation depends on a single quantity: The bifurcation takes place when The saddlenode is a codimension 1 bifurcation Example Let Find the saddlenode bifurcation set that is the locus of saddlenode bifurcations in space is defined by From we have Then The bifurcation set corresponds to the single condition
10 Corollary 84 If f satisfies the hypotheses of theorem 83 and there is a single parameter such that the transversality condition holds then a saddlenode bifurcation occurs as crosses zero Proof Recall that was defined to satisfy We also defined Thus the chain rule implies 10 Since the sign of changes as crosses zero Remark When there is a vector of parameters application of the corollary requires that for as crosses Example Let If a saddlenode bifurcation takes place as crosses Otherwise the bifurcation takes place when Example Let an unfolding of Then and By the corollary a bifurcation occurs as crosses zero For this system and a saddlenode bifurcation also takes place when Example Let an unfolding of Note that which does not satisfy the transversality condition In fact is an equilibrium of for all values of A saddlenode bifurcation clearly does not occur in this unfolding of Theorem 85 Under the hypotheses of theorem 83 the saddlenode bifurcation has the miniversal unfolding [1] Proof Expand in Taylor series about : Identifying and gives [2] Let and Then [3]
11 11 Compare [2] and [3] with the form of section 83 Eq [2] which is reproduced below To obtain that equation we assumed that the flows of and were related by the diffeomorphism where Then we obtained the relationship between vector fields Note that the expression for given below [2] is a diffeomorphism and that Eq [3] exhibits the correct constant of proportionality [3] gives The comparison between [2] and Thus the family of vector fields is induced by a family of form [4] According to the onedimensional equivalence theorem theorem 410 there is a neighborhood of the origin for which the dynamics of [4] is topologically equivalent to those of [1] because both systems have two equilibria with the same stability types and arranged in the same order on the line Transcritical Bifurcation Consider an unfolding of This is not a versal unfolding; the equilibrium point at never disappears is a normal form of the transcritical bifurcation Sketch the bifurcation diagram   two lines of equilibria with stability exchanged at the origin   with parabolas intersecting the abscissa at and to determine stability of equilibrium points This bifurcation is sometimes referred to as an exchange of stability between the two equilibrium points and is encountered frequently in applications To see the relationship with the saddlenode bifurcation begin by completing squares: Let and to see that that flow of is homeomorphic to that of The family of vector fields is induced by using
12 12 Notice that cannot be positive As increases through zero values of take the following path through the saddlenode bifurcation diagram:   curves of equilibria corresponding to SN bifurcation path of As a bug crawls along this path it sees the transcritical bifurcation diagram shown above (modulo some distortion due to the coordinate transformations) 86 SaddleNode Bifurcation in Theorem 86 (saddle node) Let and suppose that satisfies (singularity) Choose coordinates so that is diagonal in the zero eigenvalue and set where corresponds to the zero eigenvalue and are the remaining coordinates Then where and Suppose that (nondegeneracy) Then there exists an interval containing functions and and a neighborhood of such that if there are no equilibria and if there are two Suppose that has a dimensional unstable space and an dimensional stable space Then when there are two equilibria one has a dimensional unstable and an dimensional stable manifold and the other has a dimensional unstable manifold and an dimensional stable manifold Proof The equilibria are solutions of [1] By assumption that there is a neighborhood of such that is nonsingular; thus the Implicit Function Theorem ensures where there exists a unique function
13 13 [2] and Substitute this into to obtain Consequently the problem has been reduced to the onedimensional case; we need only check that satisfies the same criteria as Theorem 83 the onedimensional case It is easy to see that Since is so is and differentiation of [2] with respect to gives Since this implies that This relationship helps compute the required derivatives of : Thus the needed hypotheses for Theorem 83 are satisfied and there exists an extremal value such that when crosses zero the number of equilibria changes from zero to two The stability of the equilibria follows by considering the stability of equilibria in the one dimensional case in conjunction with the nonhyperbolic HartmanGrobman theorem from our earlier studies of the center manifold Now we can see why this is called a saddlenode bifurcation Consider theorem 86 for sketch in 1D case; two equilibria for degenerate equilibrium for no equilibria for sketch for : for a stable node and saddle; for a degenerate node; for a converging flow without equilibrium sketch for : for a saddle and an unstable node; for a degenerate node; for a diverging flow without equilibrium Tranversality The following theorem gives a condition that guarantees that parameter is varied changes sign as a
14 14 Corollary 87 (transversality) Assume that satisfies the hypotheses of theorem 86 If is any single parameter such that (transversality) then a saddlenode bifurcation takes place when crosses zero Proof Show that Use to denote the critical point of as a function of Then and The first derivatives and both vanish by hypothesis then the transversality assumption gives Example Let s apply the systematic procedure suggested in the statement of theorem 86 to the system which has an equilibrium at the origin First rewrite this in matrix form [3] The Jacobian matrix has eigenvalues and Corresponding eigenvectors are and The matrix equation has the general form where Let the transformation matrix Then Set and We have or
15 15 This has the form given in the statement of theorem 86 (see [1]) The system satisfies the singularity conditions stated in the theorem and the nondegeneracy condition Furthermore so the transversality condition is satisfied Since has a minimum and since the minimum decreases through as increases Going back to the original system [3] we can easily solve for the equilibria to get confirming our result and 85 Normal Forms In chapter 2 we transformed linear systems in order to put them in a simple form (diagonalizing the matrices in the semisimple case) In chapter 4 we learned that these changes of coordinates were diffeomorphisms In this section we continue the program of applying diffeomorphisms to transform to simpler forms only now we seek to simplify the nonlinear terms Homological Operator Let and have an equilibrium point at the origin Further assume has as many derivatives as necessary for the manipulations below Expand in power series [1] where is a vector of homogeneous polynomials of degree Let A basis for is the set of monomials
16 16 where and This compact notation is known as multiindex notation For example is three dimensional Let ( factors) be the space of vectors of homogeneous polynomials on Example has dimension 6 and the basis Denoting the standard basis vectors of by the vector monomials provide a basis for Example let and The first degree terms in the power series [1] may be written We will construct the simplest vector field that is diffeomorphic with by a near identity transformation Let represent the new variable so that and the diffeomorphism is In a small neighborhood of the equilibrium point of at the origin is close to the identity transformation and is invertible Recall from chapter 4 that if generates the flow and if generates the flow then
17 17 If there is a diffeomorphism then between the two flows so that or Set to obtain a relationship between the vector fields [2] We will choose to eliminate as many of the nonlinear terms in as possible First consider the quadratic terms setting Write We will attempt to choose so as to eliminate so that Substituting the expansions for and into [2] we find or The linear terms on the two sides of the equation match Equating the quadratic terms gives [3] which is an equation for the unknown function is called the homological operator is a linear operator on the space of degree vector fields: (see problem 6) Eq [3] can be solved if and only if where is the range of Consequently we introduce the following direct sum decomposition of : where is a complement to and we split into two parts:
18 18 The function is the resonant part of We now reconsider the derivation of Eq [3] only aiming to eliminate the nonresonant terms Let Then insert the new expansions for and into [2] to obtain which is guaranteed to have a solution for Matrix Representation Prior to applying the near identity transformations we assume that coordinates have been transformed to bring the matrix into Jordan form Let be the standard unit basis vectors in this coordinate system Any linear operator on a finite dimensional space has a matrix representation Suppose where is a vector space of homogeneous polynomials and and let of basis vectors: represent a basis for Then and is given by a linear combination This defines the matrix as a representation of the action of on Writing and we have which implies Since the basis vectors are linearly independent this is equivalent to the matrix equation The simplest case is when the eigenvalues of are real and distinct Then in Jordan form Compute the action of on the monomial basis vectors of From [3]
19 19 Since is diagonal and is proportional to we have Only the th component of is nonzero and therefore only the th row of the Jacobian matrix is nonzero: Then is the dot product of this th row with the vector The th term is this dot product is Then Thus we have The vector monomials are eigenfunctions of on with eigenvalues If all the are nonzero we can solve by inverting to obtain Example Consider the 1D case with a hyperbolic equilibrium: : Set The homological operator is so with Since the nonlinear terms are all proportional to with there are no resonant terms and all nonlinear terms may be eliminated to obtain [4] This results is consistent with the
20 20 HartmanGrobman theorem which tells us that the dynamics of [4] in a neighborhood of the origin are topologically conjugate to those of the linearized system When one or more of the are zero is nontrivial and there may be resonances For example consider the case of for which the eigenvectors of were found to be Suppose Then the condition corresponds to For the case of and Then If then is nontrivial and may have resonant nonlinear terms that cannot be eliminated This may not be a surprise since is not hyperbolic For the case of and Then Notice that this can be zero even if is hyperbolic In other words certain nonlinear systems with hyperbolic linear parts may have resonant nonlinear terms that cannot be eliminated by normal form transformations Compare this to the statement of the HartmanGrobman theorem which assures us that the dynamics of a nonlinear system in a neighborhood of a hyperbolic equilibrium point are topologically conjugate to those of the corresponding linearized system This difference reflects the fact that normal form transformations use diffeomorphisms to eliminate the nonlinear terms This is a smaller class of transformations than the homeomorphisms used by the HartmanGrobman theorem The nonhyperbolic HartmanGrobman theorem from center manifold theory tells us that the dynamics of a nonlinear system are topologically conjugate to the linearized dynamics on the stable and unstable manifold together with the nonlinear dynamics on the center manifold Thus we are particularly interested to apply normal form transformations to systems with nonhyperbolic linear parts Example Doublezero eigenvalue Consider a 2D system of form [5] The Jacobian for this system evaluated at the origin is
21 21 is the most typical Jordan form for a system with two zero eigenvalues The alternative is identically zero We attempt to eliminate the quadratic terms on the right hand side of [5] using a near identity transformation of form where and Denote The homological operator is Recalling etc we have etc In this way we construct a matrix representation for in the basis The result is The column space of defines its range The resonant space is any subspace complementary to Since it must be an element of The second vector may be any linear combination of and that is independent of The simplest choices for are and Using the second choice leads to [6]
22 22 This second form has the advantage that it is equivalent to the nonlinear oscillator The right hand side of [6] is the singular vector field whose unfolding is known as the Takens Bogdanov bifurcation which Meiss discusses in section 810 Higher Order Normal Forms Proceed by induction Suppose that all terms in the range of order To this order have been eliminated below where contains the resonant terms through order Now let and require that the vector field for have only resonant terms through order Then Substitute these expansions into the relationship between diffeomorphic vector fields that we reproduce to obtain where the homological operator is on the left hand side of the equation Set [7] and equate terms in [7] of order to obtain This equation may be solved since the right hand side is in the range of
23 Example Continue the normal form transformations for the double zero (TakensBogdanov) bifurcation to third order to obtain See Meiss Chapter 8 problem 18 23
One 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 informationIn these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation.
1 2 Linear Systems In these chapter 2A notes write vectors in boldface to reduce the ambiguity of the notation 21 Matrix ODEs Let and is a scalar A linear function satisfies Linear superposition ) Linear
More information4 SecondOrder Systems
4 SecondOrder Systems Secondorder autonomous systems occupy an important place in the study of nonlinear systems because solution trajectories can be represented in the plane. This allows for easy visualization
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 (SpringerVerlag
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 information27. Topological classification of complex linear foliations
27. Topological classification of complex linear foliations 545 H. Find the expression of the corresponding element [Γ ε ] H 1 (L ε, Z) through [Γ 1 ε], [Γ 2 ε], [δ ε ]. Problem 26.24. Prove that for any
More informationAn Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems
An Undergraduate s Guide to the HartmanGrobman and PoincaréBendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The HartmanGrobman and PoincaréBendixon Theorems
More informationNonlinear Autonomous Dynamical systems of two dimensions. Part A
Nonlinear Autonomous Dynamical systems of two dimensions Part A Nonlinear Autonomous Dynamical systems of two dimensions x f ( x, y), x(0) x vector field y g( xy, ), y(0) y F ( f, g) 0 0 f, g are continuous
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 informationDIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS
WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Series Editor: Leon O. Chua Series A Vol. 66 DIFFERENTIAL GEOMETRY APPLIED TO DYNAMICAL SYSTEMS JeanMarc Ginoux Université du Sud, France World Scientific
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 informationAnalysis of the TakensBogdanov bifurcation on m parameterized vector fields
Analysis of the TakensBogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationANSWERS Final Exam Math 250b, Section 2 (Professor J. M. Cushing), 15 May 2008 PART 1
ANSWERS Final Exam Math 50b, Section (Professor J. M. Cushing), 5 May 008 PART. (0 points) A bacterial population x grows exponentially according to the equation x 0 = rx, where r>0is the per unit rate
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 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 informationNonlinear Control Lecture 2:Phase Plane Analysis
Nonlinear Control Lecture 2:Phase Plane Analysis Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 r. Farzaneh Abdollahi Nonlinear Control Lecture 2 1/53
More informationDynamical Systems. Pierre N.V. Tu. An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition.
Pierre N.V. Tu Dynamical Systems An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition With 105 Figures SpringerVerlag Berlin Heidelberg New York London Paris
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 informationLecture2 The implicit function theorem. Bifurcations.
Lecture2 The implicit function theorem. Bifurcations. 1 Review:Newton s method. The existence theorem  the assumptions. Basins of attraction. 2 The implicit function theorem. 3 Bifurcations of iterations.
More informationDIFFERENTIAL TOPOLOGY AND THE POINCARÉHOPF THEOREM
DIFFERENTIAL TOPOLOGY AND THE POINCARÉHOPF THEOREM ARIEL HAFFTKA 1. Introduction In this paper we approach the topology of smooth manifolds using differential tools, as opposed to algebraic ones such
More informationChapter 2 Hopf Bifurcation and Normal Form Computation
Chapter 2 Hopf Bifurcation and Normal Form Computation In this chapter, we discuss the computation of normal forms. First we present a general approach which combines center manifold theory with computation
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 informationKinematics of fluid motion
Chapter 4 Kinematics of fluid motion 4.1 Elementary flow patterns Recall the discussion of flow patterns in Chapter 1. The equations for particle paths in a threedimensional, steady fluid flow are dx
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
More informationENGI 9420 Lecture Notes 4  Stability Analysis Page Stability Analysis for Nonlinear Ordinary Differential Equations
ENGI 940 Lecture Notes 4  Stability Analysis Page 4.01 4. Stability Analysis for Nonlinear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationDynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting
International Journal of Basic & Applied Sciences IJBASIJENS Vol:13 No:05 55 Dynamics of Modified LeslieGower PredatorPrey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd
More informationUNIVERSITY OF NORTH ALABAMA MA 110 FINITE MATHEMATICS
MA 110 FINITE MATHEMATICS Course Description. This course is intended to give an overview of topics in finite mathematics together with their applications and is taken primarily by students who are not
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 informationAbstract. We consider a generalised Gause predatorprey system with a generalised Holling. ax 2. We study the cases where b is positive or negative.
BIFURCATION ANALYSIS OF A PREDATORPREY SYSTEM WIT GENERALISED OLLING TYPE III FUNCTIONAL RESPONSE YANN LAMONTAGNE, CAROLINE COUTU, AND CRISTIANE ROUSSEAU Abstract. We consider a generalised Gause predatorprey
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 informationChapter 8 Equilibria in Nonlinear Systems
Chapter 8 Equilibria in Nonlinear Sstems Recall linearization for Nonlinear dnamical sstems in R n : X 0 = F (X) : if X 0 is an equilibrium, i.e., F (X 0 ) = 0; then its linearization is U 0 = AU; A =
More informationPhysics: springmass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: springmass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationCh 10.1: Two Point Boundary Value Problems
Ch 10.1: Two Point Boundary Value Problems In many important physical problems there are two or more independent variables, so the corresponding mathematical models involve partial differential equations.
More informationAS and A level Further mathematics contents lists
AS and A level Further mathematics contents lists Contents Core Pure Mathematics Book 1/AS... 2 Core Pure Mathematics Book 2... 4 Further Pure Mathematics 1... 6 Further Pure Mathematics 2... 8 Further
More informationConics and their duals
9 Conics and their duals You always admire what you really don t understand. Blaise Pascal So far we dealt almost exclusively with situations in which only points and lines were involved. Geometry would
More informationCyclicity of common slow fast cycles
Available online at www.sciencedirect.com Indagationes Mathematicae 22 (2011) 165 206 www.elsevier.com/locate/indag Cyclicity of common slow fast cycles P. De Maesschalck a, F. Dumortier a,, R. Roussarie
More informationChapter Two Elements of Linear Algebra
Chapter Two Elements of Linear Algebra Previously, in chapter one, we have considered single first order differential equations involving a single unknown function. In the next chapter we will begin to
More informationSystems of Algebraic Equations and Systems of Differential Equations
Systems of Algebraic Equations and Systems of Differential Equations Topics: 2 by 2 systems of linear equations Matrix expression; Ax = b Solving 2 by 2 homogeneous systems Functions defined on matrices
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 informationStability of Dynamical systems
Stability of Dynamical systems Stability Isolated equilibria Classification of Isolated Equilibria Attractor and Repeller Almost linear systems Jacobian Matrix Stability Consider an autonomous system u
More informationAN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b
AN ITERATION In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b In this, A is an n n matrix and b R n.systemsof this form arise
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More information1 Review of simple harmonic oscillator
MATHEMATICS 7302 (Analytical Dynamics YEAR 2017 2018, TERM 2 HANDOUT #8: COUPLED OSCILLATIONS AND NORMAL MODES 1 Review of simple harmonic oscillator In MATH 1301/1302 you studied the simple harmonic oscillator:
More informationA matrix is a rectangular array of. objects arranged in rows and columns. The objects are called the entries. is called the size of the matrix, and
Section 5.5. Matrices and Vectors A matrix is a rectangular array of objects arranged in rows and columns. The objects are called the entries. A matrix with m rows and n columns is called an m n matrix.
More informationPerturbation Theory of Dynamical Systems
Perturbation Theory of Dynamical Systems arxiv:math.ho/0111178 v1 15 Nov 2001 Nils Berglund Department of Mathematics ETH Zürich 8092 Zürich Switzerland Lecture Notes Summer Semester 2001 Version: November
More informationIntroduction to Continuous Dynamical Systems
Lecture Notes on Introduction to Continuous Dynamical Systems Fall, 2012 Lee, Keonhee Department of Mathematics Chungnam National Univeristy  1  Chap 0. Introduction What is a dynamical system? A dynamical
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 informationAN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES
AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim
More informationMatrices and Matrix Algebra.
Matrices and Matrix Algebra 3.1. Operations on Matrices Matrix Notation and Terminology Matrix: a rectangular array of numbers, called entries. A matrix with m rows and n columns m n A n n matrix : a square
More informationMath 256: Applied Differential Equations: Final Review
Math 256: Applied Differential Equations: Final Review Chapter 1: Introduction, Sec 1.1, 1.2, 1.3 (a) Differential Equation, Mathematical Model (b) Direction (Slope) Field, Equilibrium Solution (c) Rate
More informationStratification of 3 3 Matrices
Stratification of 3 3 Matrices Meesue Yoo & Clay Shonkwiler March 2, 2006 1 Warmup with 2 2 Matrices { Best matrices of rank 2} = O(2) S 3 ( 2) { Best matrices of rank 1} S 3 (1) 1.1 Viewing O(2) S 3 (
More informationTangent spaces, normals and extrema
Chapter 3 Tangent spaces, normals and extrema If S is a surface in 3space, with a point a S where S looks smooth, i.e., without any fold or cusp or selfcrossing, we can intuitively define the tangent
More informationFundamentals of Linear Algebra. Marcel B. Finan Arkansas Tech University c All Rights Reserved
Fundamentals of Linear Algebra Marcel B. Finan Arkansas Tech University c All Rights Reserved 2 PREFACE Linear algebra has evolved as a branch of mathematics with wide range of applications to the natural
More informationStability lectures. Stability of Linear Systems. Stability of Linear Systems. Stability of Continuous Systems. EECE 571M/491M, Spring 2008 Lecture 5
EECE 571M/491M, Spring 2008 Lecture 5 Stability of Continuous Systems http://courses.ece.ubc.ca/491m moishi@ece.ubc.ca Dr. Meeko Oishi Electrical and Computer Engineering University of British Columbia,
More informationTHE GEOMETRY OF GENERIC SLIDING BIFURCATIONS
THE GEOMETRY OF GENERIC SLIDING BIFURCATIONS M. R. JEFFREY AND S. J. HOGAN Abstract. Using the singularity theory of scalar functions, we derive a classification of sliding bifurcations in piecewisesmooth
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationFIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRSTORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More informationLecture 6 Positive Definite Matrices
Linear Algebra Lecture 6 Positive Definite Matrices Prof. ChunHung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Spring 2017 2017/6/8 Lecture 6: Positive Definite Matrices
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 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 informationLesson 9 Exploring Graphs of Quadratic Functions
Exploring Graphs of Quadratic Functions Graph the following system of linear inequalities: { y > 1 2 x 5 3x + 2y 14 a What are three points that are solutions to the system of inequalities? b Is the point
More informationHighOrder Representation of Poincaré Maps
HighOrder 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 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 informationDepartment of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in. NUMERICAL ANALYSIS Spring 2015
Department of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in NUMERICAL ANALYSIS Spring 2015 Instructions: Do exactly two problems from Part A AND two
More informationCollege of Natural Science Department of Mathematics
College of Natural Science Department of Mathematics Graduate Project Report on Difference Equations and Bifurcation Analysis of Discrete Dynamical System Submitted in partial fulfilment of the requirements
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationSymmetry Properties of Confined Convective States
Symmetry Properties of Confined Convective States John Guckenheimer Cornell University 1 Introduction This paper is a commentary on the experimental observation observations of Bensimon et al. [1] of convection
More informationRelation of Pure Minimum Cost Flow Model to Linear Programming
Appendix A Page 1 Relation of Pure Minimum Cost Flow Model to Linear Programming The Network Model The network pure minimum cost flow model has m nodes. The external flows given by the vector b with m
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 informationTHE CREMONA GROUP: LECTURE 1
THE CREMONA GROUP: LECTURE 1 Birational maps of P n. A birational map from P n to P n is specified by an (n + 1)tuple (f 0,..., f n ) of homogeneous polynomials of the same degree, which can be assumed
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
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 informationIt is the purpose of this chapter to establish procedures for efficiently determining the solution of a system of linear equations.
Chapter 2. Linear Systems and their Solution It is the purpose of this chapter to establish procedures for efficiently determining the solution of a system of linear equations. We know that the data within
More information9.1  Systems of Linear Equations: Two Variables
9.1  Systems of Linear Equations: Two Variables Recall that a system of equations consists of two or more equations each with two or more variables. A solution to a system in two variables is an ordered
More informationM3A23/M4A23. Specimen Paper
UNIVERSITY OF LONDON Course: M3A23/M4A23 Setter: J. Lamb Checker: S. Luzzatto Editor: Editor External: External Date: March 26, 2009 BSc and MSci EXAMINATIONS (MATHEMATICS) MayJune 2008 M3A23/M4A23 Specimen
More information4 ORTHOGONALITY ORTHOGONALITY OF THE FOUR SUBSPACES 4.1
4 ORTHOGONALITY ORTHOGONALITY OF THE FOUR SUBSPACES 4.1 Two vectors are orthogonal when their dot product is zero: v w = orv T w =. This chapter moves up a level, from orthogonal vectors to orthogonal
More informationEigenvalues, Eigenvectors, and Diagonalization
Math 240 TA: Shuyi Weng Winter 207 February 23, 207 Eigenvalues, Eigenvectors, and Diagonalization The concepts of eigenvalues, eigenvectors, and diagonalization are best studied with examples. We will
More informationQualitative Classification of Singular Points
QUALITATIVE THEORY OF DYNAMICAL SYSTEMS 6, 87 167 (2005) ARTICLE NO. 94 Qualitative Classification of Singular Points Qibao Jiang Department of Mathematics and Mechanics, Southeast University, Nanjing
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationPerturbation Theory of Dynamical Systems
Perturbation Theory of Dynamical Systems arxiv:math/0111178v1 [math.ho] 15 Nov 2001 Nils Berglund Department of Mathematics ETH Zürich 8092 Zürich Switzerland Lecture Notes Summer Semester 2001 Version:
More informationMath 312 Lecture Notes Linearization
Math 3 Lecture Notes Linearization Warren Weckesser Department of Mathematics Colgate University 3 March 005 These notes discuss linearization, in which a linear system is used to approximate the behavior
More informationMATH 369 Linear Algebra
Assignment # Problem # A father and his two sons are together 00 years old. The father is twice as old as his older son and 30 years older than his younger son. How old is each person? Problem # 2 Determine
More informationMULTIVARIABLE CALCULUS, LINEAR ALGEBRA, AND DIFFERENTIAL EQUATIONS
T H I R D E D I T I O N MULTIVARIABLE CALCULUS, LINEAR ALGEBRA, AND DIFFERENTIAL EQUATIONS STANLEY I. GROSSMAN University of Montana and University College London SAUNDERS COLLEGE PUBLISHING HARCOURT BRACE
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 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 informationTHE POINCAREHOPF THEOREM
THE POINCAREHOPF THEOREM ALEX WRIGHT AND KAEL DIXON Abstract. Mapping degree, intersection number, and the index of a zero of a vector field are defined. The PoincareHopf theorem, which states that under
More informationSPRING 2006 PRELIMINARY EXAMINATION SOLUTIONS
SPRING 006 PRELIMINARY EXAMINATION SOLUTIONS 1A. Let G be the subgroup of the free abelian group Z 4 consisting of all integer vectors (x, y, z, w) such that x + 3y + 5z + 7w = 0. (a) Determine a linearly
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 informationHomogeneous Constant Matrix Systems, Part II
4 Homogeneous Constant Matrix Systems, Part II Let us now expand our discussions begun in the previous chapter, and consider homogeneous constant matrix systems whose matrices either have complex eigenvalues
More informationIMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET
IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET This is a (not quite comprehensive) list of definitions and theorems given in Math 1553. Pay particular attention to the ones in red. Study Tip For each
More informationChapter 8. Rigid transformations
Chapter 8. Rigid transformations We are about to start drawing figures in 3D. There are no builtin routines for this purpose in PostScript, and we shall have to start more or less from scratch in extending
More informationChapters 5 & 6: Theory Review: Solutions Math 308 F Spring 2015
Chapters 5 & 6: Theory Review: Solutions Math 308 F Spring 205. If A is a 3 3 triangular matrix, explain why det(a) is equal to the product of entries on the diagonal. If A is a lower triangular or diagonal
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 informationAbstract Vector Spaces and Concrete Examples
LECTURE 18 Abstract Vector Spaces and Concrete Examples Our discussion of linear algebra so far has been devoted to discussing the relations between systems of linear equations, matrices, and vectors.
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linearwave solutions and simplewave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationMATHEMATICS COMPREHENSIVE EXAM: INCLASS COMPONENT
MATHEMATICS COMPREHENSIVE EXAM: INCLASS COMPONENT The following is the list of questions for the oral exam. At the same time, these questions represent all topics for the written exam. The procedure for
More informationThe Jordan Canonical Form
The Jordan Canonical Form The Jordan canonical form describes the structure of an arbitrary linear transformation on a finitedimensional vector space over an algebraically closed field. Here we develop
More informationA Note on the Ramsey Growth Model with the von Bertalanffy Population Law
Applied Mathematical Sciences, Vol 4, 2010, no 65, 32333238 A Note on the Ramsey Growth Model with the von Bertalanffy Population aw uca Guerrini Department of Mathematics for Economic and Social Sciences
More informationSufficient conditions for a period incrementing big bang bifurcation in onedimensional maps.
Sufficient conditions for a period incrementing big bang bifurcation in onedimensional maps. V. Avrutin, A. Granados and M. Schanz Abstract Typically, big bang bifurcation occur for one (or higher)dimensional
More informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More information