Here f and f2 are holomorphic functions of the state variables x; _x; y; _y and parameters 2 R k. These functions are \higher order" in the sense that

Size: px
Start display at page:

Download "Here f and f2 are holomorphic functions of the state variables x; _x; y; _y and parameters 2 R k. These functions are \higher order" in the sense that"

Transcription

1 Interactions of Andronov-Hopf and Bogdanov-Takens Bifurcations W. F. Langford and K. Zhan Department of Mathematics and Statistics University of Guelph Guelph, Ontario, Canada NG 2W June 5, 998 Abstract A codimension-three bifurcation, characterized by a pair of purely imaginary eigenvalues and a nonsemisimple double zero eigenvalue, arises in the study of a pair of weakly coupled nonlinear oscillators with Z 2 Z 2 symmetry. The methodology is based on Arnold's ideas of versal deformations of matrices for the linear analysis, and Poincare normal forms for the nonlinear analysis of the system. Bifurcation submanifolds of codimension one and two are identied in the parameter space. Many dierent classes of solutions are determined in the state space, including equilibria, limit cycles, invariant tori and the possibility of homoclinic chaos. As an application, a mechanism for energy transfer between two widely-spaced oscillation modes without strong resonance is identied. Introduction This study is motivated by a pair of weakly-coupled nonlinear oscillators, modeled by second-order dierential equations of the general form x +! 2 x + f (x; _x; y; _y; ) = y+! 2 2 y+f 2(x; _x; y; _y; ) =: ()

2 Here f and f2 are holomorphic functions of the state variables x; _x; y; _y and parameters 2 R k. These functions are \higher order" in the sense that they contain only quadratic and higher order terms in these variables, and f j (; ; ; ;); j =;2. The variables and parameters are all assumed to be small. The parameters represent physically relevant coecients such as damping, detuning and the coupling between the oscillators; the role of these parameters will be claried further below. In (), the natural frequencies of the two oscillators, when linearized at (x; _x; y; _y) = (; ; ; ), are given by!;2. Systems of the form () occur frequently in the study of mechanical systems, electrical circuits, the biological sciences and elsewhere. A major concern in the study of coupled oscillators is the understanding and control of resonance phenomena. V.I. Arnold has made inuential contributions to the general theory of resonances, for example [2, 5, 6]. The ratio of the natural frequencies plays a leading role in determining resonances in the system () of two oscillators. The cases! = p, where p and q are small!2 q positive integers, p + q 5, are called strong resonances. The remaining rational cases are called weak resonances, and cases with the ratio! irrational!2 are called nonresonant. As the names suggest, coupled oscillators normally interact very strongly at or near strong resonances, and interact weakly, if at all, in the other cases. This may be explained by the fact that in the Poincare normal form corresponding to these coupled oscillator equations, terms of degree p + q in the state variables appear. These are called resonant terms. Locally, these terms have a strong inuence when p + q is small and their inuence decreases as p + q increases. This paper addresses the case of coupled oscillators () with natural frequencies which are widely separated; we assume! very small and!2 very large, relative to each other. On the basis of the above classication, this would not be considered a case of strong resonance. To x magnitudes, rescale time in () so that!2 =. Then our assumption is!!. Equivalently, we consider natural frequencies at or near the ratio! =!2. Therefore, we may refer to this as the case of : resonance. In a series of papers [28, 29, 3], Nayfeh and coworkers have investigated a system of two weakly coupled oscillators with widely separated frequencies similar to (), experimentally and analytically. They point out an interesting paradox. In both mechanical experiments and a mathematical model, when the system is excited near a high natural frequency, large low-frequency responses accompanied by slow modulations of the amplitude and phase of the high-frequency mode are observed. Thus, they found that there is a strong 2

3 interaction between these two modes, and a transfer of energy from high to low frequency modes, even though by the above classication it is not a case of strong resonance. In this paper, system () is viewed as a perturbation of a singular system with a bifurcation of codimension three. The low frequency mode, in the limit!!, is a Bogdanov-Takens bifurcation [8, 3]. This bifurcation has codimension two. The high frequency mode is close to a classical Andronov- Hopf bifurcation [, 3] which has codimension one. We verify that the combined singular system has codimension three, using the method of Arnold [3], and then investigate via Poincare normal forms the interactions of these two classical bifurcations. Cases of bifurcations involving a simple zero eigenvalue with a pair of purely imaginary eigenvalues (and no other eigenvalues with zero real part) have been studied by Arnold [5], Langford [9, 2], Guckenheimer and Holmes [4], Khazin and Shnol [8], Chow et al. [] and others cited in these references. Less is known about the case of a nonsemisimple double zero eigenvalue together with an imaginary pair. Pioneering analytical results and applications have been presented by Moson [25, 26], Yu and Huseyin [33], Geogiou et al. [], Nagata and Namachchivaya [27], and Nayfeh et al. [28, 29, 3]. To x ideas, throughout this paper it is assumed that each of the oscillators in () has an odd reectional symmetry about the rest position. We extend this symmetry to the coupled system by assuming that the system commutes with the Z 2 Z 2 symmetry group generated by two reection operators f( x; _x; y; _y) = f(x; _x; y; _y) f2( x; _x; y; _y) = f2(x; _x; y; _y) f(x; _x; y; _y) = f(x; _x; y; _y) f2(x; _x; y; _y) = f2(x; _x; y; _y) : This symmetry is very common in applications; for example, the classical pendulum equation, Dung's equation and Van der Pol's equation all have odd symmetry, as do the examples in [27, 28, 29, 3]. The symmetry (2) forbids linear coupling terms (i.e. linear x; _x terms in the y-equation or viceversa) but permits nonlinear coupling (such as x 2 y in the y-equation). Before beginning to analyze (), we rewrite it with the linear part in a (2) 3

4 more explicit form. x + _x +x+ ^f(x; _x; y; _y; ) = y + 2 _y + y + ^f2(x; _x; y; _y; ) = (3) Here, ;2 are linear damping parameters and = (!=!2) 2 is a small parameter which is now allowed to be negative or zero; therefore it is better though of as representing a restoring force rather than a frequency ratio. In the engineering literature, a change in sign of from positive to negative is associated with a transition from \utter" to \divergence". These three parameters are assumed to be small, and they represent all possible linear coecients in the system consistent with the symmetries (2), after rescaling. (Recall that t has been rescaled so that!2 = ). Since the functions ^f j commute with the two reection operators (2), they are cubic or higher order in the state variables. With the denition x _x y _y C A y y2 y3 y4 C A y system (3) can be written as the following system of rst order equations, with linear part which is a perturbation of a matrix A in real Jordan form. where A = _y =[A +C(^)] y + ^f (4) C A ; C(^)= 2 C A ;^f= ^f ^f2 C A : (5) Here C(^) is the linear perturbation matrix, and ^ = (;2;) represents the three physical parameters in the system. We may simplify notation by introducing a mixed set of real/complex state variables with (y;y2)2r 2 dened as above, and z 2 C dened by z = y3 iy4; z = y3 + iy4 : (6) 4

5 (The reason for the nonstandard choice of signs of the imaginary parts of z; z is the following. In the phase plane of a second order harmonic oscillator equation, the coordinate axes are normally chosen to be (x; y) =(x; _x), which forces oscillations about the origin to have a clockwise orientation. In the complex plane, the convention is that a counterclockwise rotation is positive, corresponding to the dierential equation _z = iz. The transformation (6) eects this reversal of orientation, between the real and complex phase planes.) Then system (4) is transformed to the systen _y = y2 _y2 = y y2 ^f _z = (i 2 2)z + 2 2z + i ^f2 together with a fourth equation for _z obtained by conjugation of the equation for _z given here. The matrix A in (4) is recognized as being in real Jordan form, with the four eigenvalues 2 (nonsemisimple) and i (imaginary pair). The corresponding dierential equations (4) for small amplitudes represent a peturbation of a system of two independent oscillators, one with coordinates (x; _x) and natural frequency!, and the other with coordinates (y; _y) and natural frequency!2 =. However, the results of this paper apply more generally than this. Consider any parametrized two degree-of-freedom system, which may have coecient matrix far from A ; but suppose that, for some value of the parameters, the linearization at an equilibrium has eigenvalues f 2 ; ig (after rescaling time). Then the linear part can be transformed to A by a linear change of coordinates (similarity transformation), which may be far from the identity. In this case, the state variables in (4) are the normal mode coordinates of the original system, and the symmetries (2) are symmetries of the normal modes. The state variables in this paper should be understood in that light. Even more generally, the results of this paper apply to systems of ordinary dierential equations of dimension higher than four, and to classes of innite dimensional partial dierential equations or functional dierential equations, satisfying appropriate spectral conditions. If the linearization at an equilibrium solution satises standard compactness and hyperbolicity assumptions on the linear operator and its spectrum, then the Center Manifold Theorem applies. In this case, if the spectrum contains the eigenvalues f 2 ; ig, and the complementary spectrum is bounded away from the imaginary axis 5 (7)

6 in the left half-plane, then the study of long-time dynamics of the full system reduces to the four dimensional system considered here. Such higher dimensional generalizations will not be pursued further here. The paper is organized as follows. In the next section, a versal deformation of the linear part of the singular vector eld is obtained, using ideas of Arnold, and the codimension is computed to be three. Section 3 presents the Poincare normal form for this system. The original four-dimensional system is reduced to a three-dimensional one, exploiting the S normal form symmetry arising from the Andronov{Hopf bifurcation. In Section 4, the generic bifurcations of codimension one and two are explored. A rich variety of solutions, including stationary, periodic and quasiperiodic (tori) are found. We present formulas for primary and secondary bifurcations of periodic solutions and tori, with the help of normal forms. In Section 5, the relevance of these results to the energy-transfer phenomenon of Nayfeh is considered, and directions of further research are indicated. 2 Linear Analysis Denitions. Let C n2 denote the set of all nn complex matrices. A family A of matrices on the base M is a holomorphic mapping A : M!C n2, where M is a neighbourhood of the origin in the parameter space C k. For brevity, we call this the family A. The germ of a family A at the point 2 M is called a deformation of the matrix A(). The elements 2 M are called parameters. An example of a family of real matrices is given by A +C(^) in equation (4) and a complex family is given by the linear part of (7). The algebra is simpler if we work with complex matrices. The real Jordan matrix A of (4) has the complex Jordan form A c = i i C A : (8) Two main questions are addressed in this section. The rst question is how to nd a canonical matrix, which is as simple as possible, and to which the linear part of all coupled-oscillator systems as described in Section 6

7 can be reduced by a linear change of coordinates. An obvious candidate is the Jordan form; however, this is not a good choice because the reduction of a matrix to its Jordan form is not a stable operation when there are multiple eigenvalues, as is the case here. Although every member of a family containing A c may be reduced to a Jordan form, the Jordan form does not in general depend continuously on the parameters, since there are matrices arbitrarily close to A c for which the Jordan form is diagonal. The mapping which transforms matrices to Jordan form is discontinous. Therefore, we seek a dierent normal form, to which members of any family containing A c may be transformed smoothly. The second question relates to the choice of parameters 2M. What is the minimum set of parameters needed, and how can we know ifany essential parameters are missing? In particular, are there systems with linear parts close to A which are not attainable by variations of the physical parameters in the given perturbation matrix C(^) in (4)? If we have in fact overlooked some essential parameters, then any study of the system in Section is likely to miss important phenomena. We would like to have assurance that the parameter set in our family of matrices is suciently rich that all perturbations of A are included, and that there are no redundant parameters. Both of these questions were answered by Arnold [3, 4, 5], who introduced a generalization of the Jordan form of a matrix A, with the minimum number of essential parameters, and to which all members of any family of A can be transformed smoothly. This generalization of the Jordan normal form may be called the Arnold normal form of a matrix. The remainder of this section summarizes Arnold's construction of this normal form and shows how it answers the two questions above. The reduction to Jordan form of any matrix A is accomplished by a similarity transformation. Consider the Lie group GL(n; C ) of all nonsingular matrices S 2 C n2. Then the group of similarity transformations on C n2 is dened by the elements Ad S (A) SAS ; A 2 C n2 ; S 2 GL(n; C ) : (9) The group orbit of a given A is O(A) = Ad S (A) =SAS j S 2 GL(n; C ) () and this is a smooth submanifold of C n2. The orbit of A thus consists of all matrices similar to A. It is well known that similar matrices share the 7

8 same eigenvalues, with the same multiplicities and the same Jordan blocks (perhaps permuted). A convenient orbit representative for any given A is its Jordan form, with Jordan blocks ordered in decreasing size, for example A c in (8). Now, combining the concepts of deformations and similarity, we say that two deformations B() and A() of the same matrix B() = A() = A are equivalent if there exists a deformation C() of the identity matrix, for 2M, such that B() =C()A()C() : () Next, this equivalence is extended to allow two dierent parametrizations for B and A. Let M and N be neighbourhoods of the origin in parameter spaces C k and C l, respectively. Let ' : M!N be holomorphic at, and '() =. Let A() be a family of matrices on N. Dene the family (deformation) induced from A by ' to be the family (deformation) ' A (' A)() =A('()): Denition. A deformation A() of A is called versal if every deformation B() of Ais equivalent to a deformation induced from A(): B() =C()A('())C (); C() = I; '() = : (2) A versal deformation is universal if mapping ' is determined uniquely by deformation B(). It is miniversal if the dimension of the parameter space N is minimal among all versal deformations. This minimum dimension is called the codimension of the matrix A. The parameters in a miniversal deformation are also called unfolding parameters, and a (uni, mini-)versal deformation is sometimes called a (uni, mini-)versal unfolding. It is clear that such a versal deformation exists and an upper bound on the codimension is n 2, since one could take as parameters all the elements of an n n matrix. A lower bound is given by n, for the case of distinct eigenvalues and a diagonal Jordan form. In general, we need a manifold transversal to the group orbits. To this end, compute tangents to the group orbits in (), for xed A = A. Write (S) Ad S (A) =SAS, then the derivative with respect to S at the identity I, for xed A, is easily found to be the linear operator dened by commutation with A D (I) :C n2!c n2 ; D (I)S =[S; A] SA AS: (3) 8

9 Denition. The centralizer Z A of matrix A 2 C n2 which commute with A is the set of all matrices Z A = f B j [B; A] BA AB =g: Thus, the nullspace of D (I) is the centralizer of A. The rank of D (I) is the dimension of the orbit of A. Since D (I) is square, the dimension of the centralizer is the codimension of the orbit, which is the minimum number of parameters in a versal deformation. This codimension was given by Arnold [3] in the following formula d = sx i= l ix (2j )n ij ( i ) (4) j= where i ; i = ;:::;s are distinct eigenvalues of A, l i is the number of Jordan blocks of i and n i( i ) n ili ( i ) are the corresponding orders of the Jordan blocks for each i. Furthermore, with respect to the Hermitian scalar product ha; Bi = tr(ab ), the orthogonal complement of the tangent to the orbit of A in the tangent space C n2 is given by fb j B 2 Z A g. Here B denotes the conjugate transpose of B. This is the form that the Fredholm Alternative takes in this setting. From this, Arnold obtains explicitly the versal deformations, see [3, 4, 5, 32]. Applying Arnold's method to the Jordan matrix A c above gives the centralizer ZA c = C() and versal deformation A c (): A c () = C() = i+3 i+4 C A (5) C A : (6) Here j ; j = ;:::;4 are arbitrary complex parameters. This versal deformation is not unique. Rather than the above orthogonal complement to the 9

10 tangent space, we could choose any transversal subspace of minimum dimension, and obtain a miniversal deformation. One convenient choice replaces the 2 entry in the upper left corner of A c () by, and thus reduces the number of nonzero entries. We adopt this choice. Furthermore, since our original system is real, we decomplexify A c (). The result of decomplexication (for this choice) is the real family of matrices A () = C A (7) where j ; j = ;:::;4 have been redened as real parameters. We call the family in equation (7) the real Arnold normal form of the original matrix A. It has the desired versality property, that any deformation B() ofa is equivalent to a deformation induced from the family A. In other words, any such B() can be obtained by means of holomorphic mappings of parameters combined with similarity transformations, as given in (2). Among all such versal deformations, it has the minimum number of parameters. Finally, we recall that our matrix is the linear part of a dierential equation, in which a rescaling of time has the eect of multiplying the matrix by a positive real number. We choose to rescale time to keep the imaginary part of the complex eigenvalue equal to i. This has the eect of making 4 = in (7) and rescaling the other parameters. Therefore, in this context the codimension is 3, and the versal deformation is A () = C A (8) where now j ; j =;2;3 represent redened real unfolding parameters. It remains to determine the mapping in (2) between the unfolding parameters j in Arnold's versal normal form (8) and the physical parameters in the original system (4). This may be done by comparing coecients of the characteristic polynomials. The relationship is found to be = + 2 = + 3 = 2 2+ (9)

11 where represents higher order terms. Since this mapping is invertible at the origin, we have proven that the original family in (4), obtained via physical arguments, in fact yields aversal deformation. The 3-parameter family of matrices (8) is the appropriate starting point for study of the phenomena associated with coupled oscillators near : resonance. As observed by Arnold in [5, page 22]:... Then the natural object of study is not the individual object (say, a vector eld with a complicated singular point), but a family so large that the singularity of the type under consideration does not disappear under a small deformation of the family. This simple argument of Poincare shows the futility of such a large number of studies in the theory of dierential equations and in other areas of analysis that it is always somewhat dangerous to mention it. 3 Poincare Normal Form The Poincare normal form may be thought of as an extension of the ideas of the previous section from the linear case to the case of nonlinear vector elds. Arnold gives an excellent account of the theory and computation of the Poincare normal form in [5]. Other useful references are [7, 9,, 4, 5, 7, 8, 23, 32]. The basic idea of Poincare is to to perform a sequence of near-identity transformations of the form x = y + h k (y); k =2;3;::: (2) where h k is a homogeneous polynomial of degree k. For each k, h k is chosen to eliminate as many as possible of the terms of degree k in the vector eld. Terms which remain lie in a complement to the tangent space, at the identity, of the manifold of group orbits corresponding to (2); just as in the linear case of Section 2. The Poincare normal form is not unique, since it depends on this choice of complement. Furthermore, the resulting power series does not converge in general; nor does the sequence of near-identity transformations (2) converge. This lack of convergence is often not a problem, since in many cases the phenomena of interest are determined by a nite truncation of the Poincare normal form, of low degree.

12 For the system of dierential equations of interest in this paper as given by (7), the Poincare normal form can be chosen to have a normal form symmetry [5], and can be written (formally) _y _y2 _z A i y y2 z A 2y y + 2y2 [3 + i4] z A (2) to which may be added the conjugate of the last equation, for _z. Here i ; i =;:::;4, are real valued formal power series (which may not converge) in the Hilbert basis of invariant polynomials f y 2; zz = y2 3 + y2 4 g, i.e. i = i (y 2 ;zz); i=;:::;4: Notice the formal resemblance between (2) and the complex Arnold normal form in the linear case, given by (6). In fact, if we were to allow the terms involving j in the vector on the right-hand side of (2) to have linear parts, then the linear part of the Poincare normal form would correspond to the Arnold normal form of its linearization. The Poincare normal form (2) can be further simplied by the nearidentity transformation: (y;y2;z)!(y;y2+y2;z); which eliminates the nonlinear terms in the _y equation of (2). Now, the Poincare normal form matches the Arnold normal form as given in (7), in the sense that we can identify j = j (; ); j =;:::;4: Combining in one equation the versal deformation of the linearization, i.e. Arnold normal form (8), and this simplied Poincare normal form, _y _y2 _z A 2 i + 3 y y2 z A y + 2y2 [3 + i4] z A (22) where y;y2 2R; z 2C, and j = j (y 2 ;zz); (; )=; j=;:::;4, and we have rescaled time to make i4 =. 2

13 Finally, we decomplexify the z-equation to real form, but in polar rather than cartesian coordinates, to better exploit the normal form symmetry of the Hopf bifurcation. Let z = re i. Then (22) can be written _y = y2 _y2 = [ + ]y +[2+2]y2 _r = [3+3]r _ = +4 : (23) Note that does not appear in the right hand sides of equations (23), so the -equation _ is decoupled from the other three. This is due to the S symmetry of the normal form, induced by the Andronov{Hopf bifurcation H, and it eectively reduces the problem from four dimensions to three. Given a solution (y(t);y2(t);r(t)) of the rst three equations of (23), then the solution (t) is given by Z t (t) =+t+ 4(y(s);r 2 2 (s)) ds : All but exotic behaviour of this system is already determined by a cubic truncation of the Poincare normal form. This can be written _y = y2 _y2 = [ + y 2 + r 2 ]y +[2+2y 2 + 2r 2 ]y2 _r = [3 + 3y 2 + 3r 2 ]r (24) where jkl is the coecient, in the power series expansion of j (y 2;r2 ), of the term y 2k r2l. The equation has been omitted in (24), but has a similar expansion. The system (24) is explored in the next section. As has been shown in [2, 2], in special cases the higher order \tail" of the normal form can dramatically change the dynamics from the behaviour of a truncated normal form such as (24), since the tail breaks the normal form symmetry which led to the decoupling of the _ equation. For example, when (24) has homoclinic or heteroclinic orbits, then the tail eects can produce homoclinic tangles, \horseshoes" and chaotic dynamics. These considerations will not be pursued in this paper. 4 Bifurcations of Codimension One and Two The collection of all possible phase portraits of this system near the codimension three bifurcation point, or even of all the equivalence classes with respect 3

14 to topological equivalence, is too large to list here. However, the stratied subvariety [5, page 227] of the bifurcation diagram in the parameter space R 3 is relatively simple, and is presented in this section. The stability of the trivial equilibrium solution is determined by the eigenvalues of the linear part A (), which are easily seen to be ;2 = 2 2 r ; 3;4 = 3 i : (25) All four eigenvalues have negative real part if and only if j ; j = ; 2; 3 are all negative. Therefore, the trivial equilibrium solution is asymptotically stable inside the negative octant of the parameter space N 2 R 3, see Figure. As the point(;2;3) leaves the negative octant, bifurcations occur. These bifurcations from the trivial equilibrium solution are called primary bifurcations. There are seven possibilities, see Figure : a point may leave by crossing through one of the three coordinate planes (codimension one), through one of the three coordinate axes (codimension two) or through the origin (codimension three). A generic exit is through a plane, codimension one. However, interesting and important mode-interactions occur around the intersections of the planes. These are best understood by studying the codimension two bifurcations on the coordinate axes. All six of these codimension one and two bifurcations have been studied in general, and are well understood. Their manifestations in the present system are described in the following subsections. 4. Primary Bifurcations of Codimension One The three primary bifurcations of codimension one are: an Andronov{Hopf bifurcation with frequency of order O() which we label H, an Andronov{ Hopf bifurcation with small frequency O(jj), labeled H, and a pitchfork bifurcation labeled P; see Figure. 4.. Primary Hopf Bifurcation H For the system (22), the z{plane is a two-dimensional invariant subspace, on which the normal form equations reduce to y = y2 = _r = [3 + 3(;r 2 )] r _ = +4(;r 2 ) : 4 (26)

15 This has the family of periodic solutions y = y2 = r = q (t) = +[+4(;r 2 )] t (mod 2) : (27) This is a classical Andronov{Hopf bifurcation, giving birth to a limit cycle in the (y3;y4){plane with period near 2 (angular frequency ) as the parameter 3 crosses through, provided ;2 6=. Generically 3 6= and the bifurcation is nondegenerate, which is assumed to be the case in this paper; then the limit cycle is locally unique. Degenerate cases of Hopf bifurcation are explored in detail in [3]. If 3 <, then the bifurcation is supercritical (the limit cycle exits for >) and if 3 > the bifurcation is subcritical (exists for <). A supercritical limit cycle is stable if ;2 <, otherwise the cycle is unstable Pitchfork Bifurcation P Any equilibrium solution of (22) lies on the y axis, and satises + (y 2 ; ) y =: (28) This is a classic pitchfork bifurcation, yielding in addition to the trivial solution y =, a symmetric pair of nontrivial solutions y = r + branching from as crosses through zero, provided 6= and 2;3 6=, which we assume. The pitchfork is supercritical if < and subcritical if >. Supercritical solutions are stable if 2;3 <, otherwise both are unstable Primary Hopf Bifurcation H The Poincare normal form (22) leaves the (y;y2){plane invariant. On this plane, the normal form reduces to _y = y2 _y2 = [ + (y 2 ; )]y +[2+2(y 2 ;)]y 2 z = : (29) 5

16 For xed <, this system has an Andronov{Hopf bifurcation as 2 crosses through zero, for any 3 6=. The new limit cycle has a small frequency, asymptotic to p. The change of variables x = y; y = y2=! where! 2 = transforms (29) to the planar system _x! x = + _y! y! (x 2 ;)x + 2(x 2 : (3) ; )y This planar system has the Andronov{Hopf normal form, with x + iy = e i' _ = 2 _' = q b 2 + : (3) From this it follows that the limit cycle has amplitude given by = r and is super-(sub-)critical if 2 < (> ). 4.2 Primary Bifurcations of Codimension Two The Poincare normal form (22) has three bifurcations of codimension two from the trivial solution in state space. In the parameter space they occur along the three coordinate axes, see Figure Bogdanov{Takens Bifurcation BT As observed for equation (29), the (y;y2){plane is invariant for the ow of the normal form. On this plane, for 3 6= and (;2) in a neighbourhood of (; ), the system has a Z 2 {symmetric Bogdanov{Takens bifurcation [8, 3], see line BT in Figure. It yields stable solutions only if 3 <. In addition to the codimension one primary bifurcations P and H described above, there are rich possibilities including simultaneous secondary Andronov{Hopf bifurcations from the pitchfork pair of equilibria, coalescence of two limit cycles, orbits homoclinic to the origin and heteroclinic orbits joining the two pitchfork equilibria. A detailed description of the versal deformations, bifurcation diagrams and phase portraits is given by Arnold [5, page 3], and need not be repeated here. 6

17 New possibilities arise in this system via the interaction of the Bogdanov{ Takens bifurcation with the primary Andronov{Hopf bifurcation H, leading toaow in a phase space which is essentially R 2 S. The simple homoclinic orbits of the planar Bogdanov{Takens bifurcation then lead to transversally intersecting two-dimensional stable and unstable manifolds of periodic orbits and thus to chaotic dynamics, as for example in [4, 23, 32]. These phenomena merit further investigation Pitchfork{Hopf Mode Interaction PH In a neighbourhood of the 2{axis for 2 6=, the normal form (22) has a pitchfork-hopf mode interaction, see PH in Figure. On the 2{axis, the system has a simple zero eigenvalue, with a Z 2 {symmetry on the corresponding eigenspace y, and a pair of purely imaginary eigenvalues i. Near (;3)=(;) with 2 <, one nds the expected pitchfork and Hopf primary bifurcations P and H, and additional mixed-mode solutions due to the nonlinear coupling of these primary bifurcations. A complete classication of nondegenerate pitchfork{hopf mode interactions is given in [22]. To investigate these here, in (23) we set _y = y2 = ; _r =, suppress the _ equation and retain only leading order terms in and 3. Then equation (24) reduces to = [ + y 2 + r 2 ] y = [3 + 3y 2 + 3r 2 ] r : (32) For nondegeneracy, we assume 3 3 6= ; 3 6= : (33) Then (32) has four branches of solutions (y 2 ;r2 ) = (; ) (y 2 ;r2 ) = ; (y 2 ;r2 ) = (y 2 ;r2 ) = 3 ; ; 3 3 : (34) 7

18 Note that the nonzero quantities on the right must be positive for the solutions to be real. The rst three of these solutions are the previously identied trivial, primary P and primary H solutions. The fourth is a mixed-mode solution, which links the P and H solutions, at the two ends of its interval of (real) existence, in secondary bifurcations. Due to symmetry, the primary bifurcations have two symmetric branches y; r, and the secondary bifurcations would seem to have four branches; however, in the full dynamics the solutions with r give the same periodic orbit with a phase shift. The secondary bifurcation from P is an Andronov{Hopf bifurcation, while the secondary bifurcation from H is a pitchfork bifurcation of two new limit cycles. While all of the above results were deduced using only the cubic truncation (32), they extend by standard arguments to the full equations, using assumptions (33), the Implicit Function Theorem and transversality. When 3 < and higher order terms are included, it is possible to show that the phase portraits may include also a ow on a two-dimensional invariant torus with a small second frequency. The eects of the symmetry{breaking normal form tail may then lead to chaos, see [4, 6, 9, 2, 2] Hopf{Hopf Mode Interaction HH In a neighbourhood of the negative semi-axis (see HH in Figure ), nonresonant Hopf-Hopf mode interactions occur, involving the H and H bifurcations. For (2;3) =(;), the full system has two pairs of purely imaginary eigenvalues, which are nonresonant because the frequency ratio is very small. The four-dimensional system can be written as two pairs of amplitude-phase equations, in (r; ) and (; '), corresponding to the two primary Andronov{ Hopf bifurcations H and H respectively. The two phase equations decouple due to nonresonance, leaving two amplitude equations whichhave the Z 2 Z 2 symmetric normal form _r = [3 + c 2 + 3r 2 + ]r _ = dr 2 + : (35) This has four families of solutions satisfying (_r; _) =(;), given to leading order by (r; ) = (; ) 8

19 (r;) = (r;) = (r;) = r 3 3 ; ;r c ; s 3d 2 23A 2 (36) where 2 = 8 23 cd 6=. The rst three of these are the previously obtained trivial, H and H solutions, respectively. The fourth gives a new family of two-frequency solutions, on an invariant 2-torus. The local persistence of this invariant torus is guaranteed by standard arguments involving normal hyperbolicity and nonresonance. This torus appears through a secondary bifurcation from either of the primary Andronov{Hopf families. To leading order, the 2-torus bifurcates from the H periodic orbit when 3d = 2 23; =; r=r 3 3 ; and bifurcates from the H periodic orbit when 2c = 4 32; = r 42 2 ; r =: Near H, the solutions on the 2-torus are close to the H periodic orbit, with primary frequency O(), but have an amplitude modulation with frequency O( p ). A description of possible dynamics in nonresonant Hopf-Hopf mode interactions is given in [4]. In certain cases, there is a further possibility ofastable 3-torus bifurcation from the 2-torus, and even of chaotic dynamics [6]. This analysis requires consideration of higher order terms. 5 Conclusions We may draw several conclusions from this study. The rst important conclusion is that the simultaneous occurrence of Andronov{Hopf and Bogdanov{Takens bifurcations is a codimension-three phenomenon, requiring three independent parameters for its exploration and 9

20 understanding. The three parameters in Arnold's versal deformation may be chosen so that they directly determine the three primary bifurcations of codimension one: determines the primary pitchfork bifurcation, 2 determines the primary Hopf bifurcation with low frequency (of order p ) and 3 determines the primary Hopf bifurcation with high frequency (of order ). These three versal unfolding parameters are related to the physical parameters of the coupled oscillators by formulae (9). In the versal deformation of the original codimension three bifurcation, there are primary bifurcations of six dierent types; three of codimension one and three of codimension two, all of which are understood. The codimension two bifurcations lead to secondary bifurcations, and produce new families of solutions which are linkages between primary bifurcation families. The double zero eigenvalue may lead to a pitchfork bifurcation or a Hopf bifurcation of low frequency, and may lead to simultaneous existence of both nontrivial equilibria and low-frequency limit cycles, as in the standard Bogdanov{Takens analysis. When coupled with the primary Andronov{Hopf bifurcation H, it may yield homoclinic chaos [32]. The nonlinear coupling between the primary high frequency and low frequency Andronov{Hopf bifurcations can lead to a quasiperiodic ow on a 2-torus, with its two frequencies close to those of the two primary Andronov{ Hopf bifurcations. This gives a plausible explanation for the observations of Nayfeh et al. [28, 29, 3], of low frequency modulation of a high frequency mode, and of energy transfer from high to low frequency modes, near : resonance. The case studied by Nayfeh corresponds to a choice of unfolding parameters < (xed) and 2;3 both small and negative; i.e., near the negative axis marked HH in Figure. This is the location where birth of a 2-torus is expected, from the normal form analysis. A more detailed investigation of this particular application, using both analytical and numerical techniques, is in progress. The present study provides indications of other directions for future research. Systems of coupled oscillators with less symmetry than required by (2) also are important in applications, and maybeinvestigated by these same normal form methods. The eect of higher order terms, particularly in the onset of more complex dynamics, is worthy of investigation. As well, the application of these results to concrete systems in mechanical or electrical engineering, and other elds, may lead to new insights. 2

21 References [] A.A. Andronov, A.A Vitt and S.E. Khaikin. Theory of Oscillators. Dover, New York, 987. [2] V.I. Arnold. Small denominators I. Mappings of the circumference onto itself. AMS Translations (2) 46 (965) [3] V.I. Arnold. On matrices depending on parameters. Russian Mathematical Surveys 26 (97) [4] V.I. Arnold. Lectures on bifurcations in versal families. Russian Mathematical Surveys 27 (972) [5] V. I. Arnold. Geometrical Methods in the Theory of Ordinary Dierential Equations. Springer-Verlag, New York, 983. [6] V. I. Arnold. Mathematical Methods of Classical Mechanics. Second edition, Springer-Verlag, New York, 989. [7] Yuri N. Bibikov. Local Theory of Nonlinear Analytic Ordinary Dierential Equations. LNM 72, Springer-Verlag, New York, 979. [8] R. I. Bogdanov. Versal deformations of a singular point on the plane in the case of zero eigenvalues. Functional Analysis and Its Applications 9 (975) [9] Alexander D. Bruno. Local Methods in Nonlinear Dierential Equations. Springer-Verlag, Berlin 989. [] S-N. Chow, C. Li and D. Wang. Normal forms and bifurcations of planar vector elds. Cambridge University Press, 994. [] I. T. Georgiou, A. K. Bajaj and M. Corless. Slow and fast invariant manifolds, and normal modes in a two-degree-of-freedom structural dynamical system with multiple equilibrium states. Int. J. Non-Linear Mechanics 33 (998) [2] S. A. van Gils, M.P. Krupa and W. F. Langford. Hopf bifurcation with nonsemisimple : resonance. Nonlinearity 3 (99)

22 [3] M. Golubitsky and W. F. Langford. Classication and unfoldings of degenerate Hopf bifurcations. J. Di. Eq. 4 (98) [4] J. Guckenheimer and P. Holmes. Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields. Springer-Verlag, New York, 983. [5] G. Iooss and M. Adelmeyer. Topics in Bifurcation Theory and Applications. World Sci., Singapore, 992. [6] G. Iooss and W. F. Langford. Conjectures on the routes to turbulence via bifurcations. In Nonlinear Dynamics, Annals N.Y. Acad. Sci. 357 (98) [7] Peter B. Kahn and Yair Zarmi. Nonlinear Dynamics: Exploration through Normal Forms. Wiley-Interscience, New York, 998. [8] L. G. Khazin and E. E. Shnol. Stability of Critical Equilibrium States. Manchester Univ. Press, Manchester, 99. [9] W. F. Langford. Periodic and steady mode interactions lead to tori. SIAM J. Appl. Math. 37 (979) [2] W. F. Langford. A review of interactions of Hopf and steady-state bifurcations. In Dynamical Systems and Turbulence. G.I. Barenblatt, G. Iooss and D.D. Joseph editors, Pitman, London, 983, [2] W. F. Langford. Numerical studies of torus bifurcations. In Numerical Methods for Bifurcation Problems. T. Kupper, H.D. Mittelmann and H. Weber editors, ISNM 7 Birkhauser Verlag, Boston, 984, [22] W. F. Langford and G. Iooss. Interactions of Hopf and pitchfork bifurcations. In Bifurcation Problems and their Numerical Solution. H.D. Mittelmann and H. Weber editors, ISNM 54, Birkhauser Verlag, Basel, 98, [23] W. F. Langford and W. Nagata. Normal Forms and Homoclinic Chaos, Fields Institute Communications 4. American Mathematical Society, Providence,

23 [24] Luo Dingjun, Wang Xian, Zhu Deming and Han Maoan. Bifurcation Theory and Methods of Dynamical Systems. World Scientic, Singapore, 997. [25] P. Moson. Quasi-periodic solutions of a special system. In Dierential Equations: Qualitative Theory. North-Holland, 984. [26] P. Moson. Local bifurcations in the case of eigenvalues 2 ; i. Z. Angew. Math. Mech. 7 (99) T69-T7. [27] W. Nagata and N. Sri Namachchivaya. Bifurcations in gyroscopic systems with an application to rotating shafts. Proc. R. Soc. Lond. A 454 (998) [28] S. A. Nayfeh and A. H. Nayfeh. Nonlinear interactions between two widely spaced modes { external excitation. International Journal of Bifurcation and Chaos 3 (993) [29] S. A. Nayfeh and A. H. Nayfeh. Energy Transfer from high- to lowfrequency modes in a exible structure via modulation. Journal of Vibration and Accoustics 6 (994) [3] A. H. Nayfeh, S. A. Nayfeh, T. A. Anderson and B. Balachandran. Transfer of energy from high-frequency to low-frequency modes. In Nonlinearity and Chaos in Engineering Dynamics. John-Wiley & Sons, Chichester, 994. [3] F. Takens. Singularities of vector elds. Publ. Math. IHES. 43 (974) 47-. [32] S. Wiggins. Introduction to Applied Nonlinear Dynamical Systems and Chaos. Springer-Verlag, New York, 99. [33] P. Yu and K. Huseyin. Bifurcations associated with a double zero and a pair of pure imaginary eigenvalues. SIAM J. Appl. Math. 48 (988)

24 α 3 α 2 HH O α H H P PH BT Figure : 24

Interactions of Andronov{Hopf. and Bogdanov{Takens Bifurcations. Abstract. A codimension-three bifurcation, characterized by a pair of purely

Interactions of Andronov{Hopf. and Bogdanov{Takens Bifurcations. Abstract. A codimension-three bifurcation, characterized by a pair of purely Interactions of Andronov{Hopf and Bogdanov{Takens Bifurcations William F. Langford and Kaijun Zhan Abstract. A codimension-three bifurcation, characterized by a pair of purely imaginary eigenvalues and

More information

8.1 Bifurcations of Equilibria

8.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 integro-differential equations

More information

Additive resonances of a controlled van der Pol-Duffing oscillator

Additive resonances of a controlled van der Pol-Duffing oscillator Additive resonances of a controlled van der Pol-Duffing oscillator This paper has been published in Journal of Sound and Vibration vol. 5 issue - 8 pp.-. J.C. Ji N. Zhang Faculty of Engineering University

More information

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay

Difference Resonances in a controlled van der Pol-Duffing oscillator involving time. delay Difference Resonances in a controlled van der Pol-Duffing oscillator involving time delay This paper was published in the journal Chaos, Solitions & Fractals, vol.4, no., pp.975-98, Oct 9 J.C. Ji, N. Zhang,

More information

Introduction to Applied Nonlinear Dynamical Systems and Chaos

Introduction to Applied Nonlinear Dynamical Systems and Chaos Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium

More information

Elements of Applied Bifurcation Theory

Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Second Edition With 251 Illustrations Springer Preface to the Second Edition Preface to the First Edition vii ix 1 Introduction to Dynamical Systems

More information

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION

CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION CHAPTER 6 HOPF-BIFURCATION IN A TWO DIMENSIONAL NONLINEAR DIFFERENTIAL EQUATION [Discussion on this chapter is based on our paper entitled Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation,

More information

Elements of Applied Bifurcation Theory

Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Introduction to Dynamical Systems 1 1.1 Definition of a dynamical system 1 1.1.1 State space 1 1.1.2

More information

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics

Chapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the

More information

A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in

A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March 22, 2000 Abstract Existence of chaotic dynamics in A Chaotic Phenomenon in the Power Swing Equation Umesh G. Vaidya R. N. Banavar y N. M. Singh March, Abstract Existence of chaotic dynamics in the classical swing equations of a power system of three interconnected

More information

One Dimensional Dynamical Systems

One Dimensional Dynamical Systems 16 CHAPTER 2 One Dimensional Dynamical Systems We begin by analyzing some dynamical systems with one-dimensional phase spaces, and in particular their bifurcations. All equations in this Chapter are scalar

More information

Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University

Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Towards a Global Theory of Singularly Perturbed Dynamical Systems John Guckenheimer Cornell University Dynamical systems with multiple time scales arise naturally in many domains. Models of neural systems

More information

27. Topological classification of complex linear foliations

27. 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 information

Hamiltonian Systems with Widely Separated. Frequencies. Mathematisch Instituut, Utrecht University

Hamiltonian Systems with Widely Separated. Frequencies.    Mathematisch Instituut, Utrecht University Hamiltonian Systems with Widely Separated Frequencies J.M. Tuwankotta F. Verhulst e-mail: tuwankotta@math.uu.nl, verhulst@math.uu.nl. Mathematisch Instituut, Utrecht University PO Box 80.010, 3508 TA Utrecht

More information

(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1,

(8.51) ẋ = A(λ)x + F(x, λ), where λ lr, the matrix A(λ) and function F(x, λ) are C k -functions with k 1, 2.8.7. Poincaré-Andronov-Hopf Bifurcation. In the previous section, we have given a rather detailed method for determining the periodic orbits of a two dimensional system which is the perturbation of a

More information

2 Lecture 2: Amplitude equations and Hopf bifurcations

2 Lecture 2: Amplitude equations and Hopf bifurcations Lecture : Amplitude equations and Hopf bifurcations This lecture completes the brief discussion of steady-state bifurcations by discussing vector fields that describe the dynamics near a bifurcation. From

More information

NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM

NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute

More information

Torus Maps from Weak Coupling of Strong Resonances

Torus Maps from Weak Coupling of Strong Resonances Torus Maps from Weak Coupling of Strong Resonances John Guckenheimer Alexander I. Khibnik October 5, 1999 Abstract This paper investigates a family of diffeomorphisms of the two dimensional torus derived

More information

Bifurcations of phase portraits of pendulum with vibrating suspension point

Bifurcations of phase portraits of pendulum with vibrating suspension point Bifurcations of phase portraits of pendulum with vibrating suspension point arxiv:1605.09448v [math.ds] 9 Sep 016 A.I. Neishtadt 1,,, K. Sheng 1 1 Loughborough University, Loughborough, LE11 3TU, UK Space

More information

Linear Algebra. Min Yan

Linear Algebra. Min Yan Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................

More information

On low speed travelling waves of the Kuramoto-Sivashinsky equation.

On low speed travelling waves of the Kuramoto-Sivashinsky equation. On low speed travelling waves of the Kuramoto-Sivashinsky equation. Jeroen S.W. Lamb Joint with Jürgen Knobloch (Ilmenau, Germany) Marco-Antonio Teixeira (Campinas, Brazil) Kevin Webster (Imperial College

More information

Im + α α. β + I 1 I 1< 0 I 1= 0 I 1 > 0

Im + α α. β + I 1 I 1< 0 I 1= 0 I 1 > 0 ON THE HAMILTONIAN ANDRONOV-HOPF BIFURCATION M. Olle, J. Villanueva 2 and J. R. Pacha 3 2 3 Departament de Matematica Aplicada I (UPC), Barcelona, Spain In this contribution, we consider an specic type

More information

7 Planar systems of linear ODE

7 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 information

Hopf bifurcation in coupled cell networks with abelian symmetry 1

Hopf bifurcation in coupled cell networks with abelian symmetry 1 Hopf bifurcation in coupled cell networks with abelian symmetry 1 Ana Paula S. Dias e Rui C. Paiva CMUP Departamento de Matemática Universidade do Porto CMUP Instituto Politécnico de Leiria Abstract We

More information

On the Takens-Bogdanov Bifurcation in the Chua s Equation

On the Takens-Bogdanov Bifurcation in the Chua s Equation 1722 PAPER Special Section on Nonlinear Theory and Its Applications On the Takens-Bogdanov Bifurcation in the Chua s Equation Antonio ALGABA, Emilio FREIRE, Estanislao GAMERO, and Alejandro J. RODRÍGUEZ-LUIS,

More information

Symmetry Properties of Confined Convective States

Symmetry 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 information

QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those

QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE. Denitizable operators in Krein spaces have spectral properties similar to those QUASI-UNIFORMLY POSITIVE OPERATORS IN KREIN SPACE BRANKO CURGUS and BRANKO NAJMAN Denitizable operators in Krein spaces have spectral properties similar to those of selfadjoint operators in Hilbert spaces.

More information

Ordinary Differential Equations and Smooth Dynamical Systems

Ordinary Differential Equations and Smooth Dynamical Systems D.V. Anosov S.Kh. Aranson V.l. Arnold I.U. Bronshtein V.Z. Grines Yu.S. Il'yashenko Ordinary Differential Equations and Smooth Dynamical Systems With 25 Figures Springer I. Ordinary Differential Equations

More information

OSCILLATORS WITH DELAY COUPLING

OSCILLATORS WITH DELAY COUPLING Proceedings of DETC'99 1999 ASME Design Engineering Technical Conferences September 1-15, 1999, Las Vegas, Nevada, USA D E T C 9 9 / V IB -8 3 1 8 BIFURCATIONS IN THE DYNAMICS OF TWO COUPLED VAN DER POL

More information

VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE

VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE PHYSCON 2013 San Luis Potosí México August 26 August 29 2013 VERSAL DEFORMATIONS OF BILINEAR SYSTEMS UNDER OUTPUT-INJECTION EQUIVALENCE M Isabel García-Planas Departamento de Matemàtica Aplicada I Universitat

More information

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space

Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) 1.1 The Formal Denition of a Vector Space Linear Algebra (part 1) : Vector Spaces (by Evan Dummit, 2017, v. 1.07) Contents 1 Vector Spaces 1 1.1 The Formal Denition of a Vector Space.................................. 1 1.2 Subspaces...................................................

More information

COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS

COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL PERIODIC FORCING TERMS International J. of Math. Sci. & Engg. Appls. (IJMSEA) ISSN 0973-9424, Vol. 9 No. III (September, 2015), pp. 197-210 COMPLEX DYNAMICS AND CHAOS CONTROL IN DUFFING-VAN DER POL EQUATION WITH TWO EXTERNAL

More information

5.2.2 Planar Andronov-Hopf bifurcation

5.2.2 Planar Andronov-Hopf bifurcation 138 CHAPTER 5. LOCAL BIFURCATION THEORY 5.. Planar Andronov-Hopf bifurcation What happens if a planar system has an equilibrium x = x 0 at some parameter value α = α 0 with eigenvalues λ 1, = ±iω 0, ω

More information

Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields

Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),

More information

QUARTERLY OF APPLIED MATHEMATICS

QUARTERLY OF APPLIED MATHEMATICS QUARTERLY OF APPLIED MATHEMATICS Volume LIV December 1996 Number 4 DECEMBER 1996, PAGES 601-607 ON EXISTENCE OF PERIODIC ORBITS FOR THE FITZHUGH NERVE SYSTEM By S. A. TRESKOV and E. P. VOLOKITIN Institute

More information

Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling

Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling Journal of Computational and Applied Mathematics 144 (2002) 359 373 wwwelseviercom/locate/cam Computation of the simplest normal forms with perturbation parameters based on Lie transform and rescaling

More information

Example of a Blue Sky Catastrophe

Example of a Blue Sky Catastrophe PUB:[SXG.TEMP]TRANS2913EL.PS 16-OCT-2001 11:08:53.21 SXG Page: 99 (1) Amer. Math. Soc. Transl. (2) Vol. 200, 2000 Example of a Blue Sky Catastrophe Nikolaĭ Gavrilov and Andrey Shilnikov To the memory of

More information

MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione

MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione MODELLING OF FLEXIBLE MECHANICAL SYSTEMS THROUGH APPROXIMATED EIGENFUNCTIONS L. Menini A. Tornambe L. Zaccarian Dip. Informatica, Sistemi e Produzione, Univ. di Roma Tor Vergata, via di Tor Vergata 11,

More information

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting

Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:05 55 Dynamics of Modified Leslie-Gower Predator-Prey Model with Predator Harvesting K. Saleh Department of Mathematics, King Fahd

More information

TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1

TWELVE 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 information

Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber

Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber Suppression of the primary resonance vibrations of a forced nonlinear system using a dynamic vibration absorber J.C. Ji, N. Zhang Faculty of Engineering, University of Technology, Sydney PO Box, Broadway,

More information

On the convergence of normal forms for analytic control systems

On the convergence of normal forms for analytic control systems Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department

More information

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)

Vector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition) Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 5. Global Bifurcations, Homoclinic chaos, Melnikov s method Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Motivation 1.1 The problem 1.2 A

More information

Stability and bifurcation of a simple neural network with multiple time delays.

Stability and bifurcation of a simple neural network with multiple time delays. Fields Institute Communications Volume, 999 Stability and bifurcation of a simple neural network with multiple time delays. Sue Ann Campbell Department of Applied Mathematics University of Waterloo Waterloo

More information

Chapter 2 Hopf Bifurcation and Normal Form Computation

Chapter 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 information

! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»»

! 4 4! o! +! h 4 o=0! ±= ± p i And back-substituting into the linear equations gave us the ratios of the amplitudes of oscillation:.»» = A p e i! +t»» Topic 6: Coupled Oscillators and Normal Modes Reading assignment: Hand and Finch Chapter 9 We are going to be considering the general case of a system with N degrees of freedome close to one of its stable

More information

A Brief Outline of Math 355

A Brief Outline of Math 355 A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting

More information

Elements of Applied Bifurcation Theory

Elements of Applied Bifurcation Theory Yuri A. Kuznetsov Elements of Applied Bifurcation Theory Third Edition With 251 Illustrations Springer Yuri A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3584 CD Utrecht The

More information

HOPF BIFURCATION AND HETEROCLINIC CYCLES IN A CLASS OF D 2 EQUIVARIANT SYSTEMS

HOPF BIFURCATION AND HETEROCLINIC CYCLES IN A CLASS OF D 2 EQUIVARIANT SYSTEMS HOPF BIFURCATION AND HETEROCLINIC CYCLES IN A CLASS OF D 2 EQUIVARIANT SYSTEMS ADRIAN C. MURZA Communicated by Vasile Brînzănescu In this paper, we analyze a generic dynamical system with D 2 constructed

More information

Hamiltonian Dynamics

Hamiltonian Dynamics Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;

More information

4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial

4.1 Eigenvalues, Eigenvectors, and The Characteristic Polynomial Linear Algebra (part 4): Eigenvalues, Diagonalization, and the Jordan Form (by Evan Dummit, 27, v ) Contents 4 Eigenvalues, Diagonalization, and the Jordan Canonical Form 4 Eigenvalues, Eigenvectors, and

More information

Application demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems

Application demonstration. BifTools. Maple Package for Bifurcation Analysis in Dynamical Systems Application demonstration BifTools Maple Package for Bifurcation Analysis in Dynamical Systems Introduction Milen Borisov, Neli Dimitrova Department of Biomathematics Institute of Mathematics and Informatics

More information

7 Two-dimensional bifurcations

7 Two-dimensional bifurcations 7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed

More information

Linear Algebra: Matrix Eigenvalue Problems

Linear Algebra: Matrix Eigenvalue Problems CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given

More information

CENTER MANIFOLD AND NORMAL FORM THEORIES

CENTER MANIFOLD AND NORMAL FORM THEORIES 3 rd Sperlonga Summer School on Mechanics and Engineering Sciences 3-7 September 013 SPERLONGA CENTER MANIFOLD AND NORMAL FORM THEORIES ANGELO LUONGO 1 THE CENTER MANIFOLD METHOD Existence of an invariant

More information

An introduction to Birkhoff normal form

An 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 information

Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems

Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems Stability Computations for Nilpotent Hopf Bifurcations in Coupled Cell Systems M. Golubitsky M. Krupa April 18, 006 Abstract Vanderbauwhede and van Gils, Krupa, and Langford studied unfoldings of bifurcations

More information

B5.6 Nonlinear Systems

B5.6 Nonlinear Systems B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre

More information

Shayne Waldron ABSTRACT. It is shown that a linear functional on a space of functions can be described by G, a

Shayne Waldron ABSTRACT. It is shown that a linear functional on a space of functions can be described by G, a UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Symmetries of linear functionals Shayne Waldron (waldron@math.wisc.edu) CMS Technical Summary Report #95-04 October 1994 ABSTRACT It

More information

On a Codimension Three Bifurcation Arising in a Simple Dynamo Model

On a Codimension Three Bifurcation Arising in a Simple Dynamo Model On a Codimension Three Bifurcation Arising in a Simple Dynamo Model Anne C. Skeldon a,1 and Irene M. Moroz b a Department of Mathematics, City University, Northampton Square, London EC1V 0HB, England b

More information

Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l = 1, 2, 3.

Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l = 1, 2, 3. June, : WSPC - Proceedings Trim Size: in x in SPT-broer Survey of strong normal-internal k : l resonances in quasi-periodically driven oscillators for l =,,. H.W. BROER and R. VAN DIJK Institute for mathematics

More information

Defining Equations for Bifurcations and Singularities

Defining Equations for Bifurcations and Singularities Defining Equations for Bifurcations and Singularities John Guckenheimer and Yongjian Xiang Mathematics Department, Ithaca, NY 14853 For Vladimir Arnold on the occasion of his 65th birthday July 1, 2002

More information

APPPHYS217 Tuesday 25 May 2010

APPPHYS217 Tuesday 25 May 2010 APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag

More information

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip

More information

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi

Abstract. Jacobi curves are far going generalizations of the spaces of \Jacobi Principal Invariants of Jacobi Curves Andrei Agrachev 1 and Igor Zelenko 2 1 S.I.S.S.A., Via Beirut 2-4, 34013 Trieste, Italy and Steklov Mathematical Institute, ul. Gubkina 8, 117966 Moscow, Russia; email:

More information

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits: equilibria cycles connecting orbits compact invariant manifolds strange

More information

The Hopf bifurcation with bounded noise

The Hopf bifurcation with bounded noise The Hopf bifurcation with bounded noise Ryan T. Botts 1, Ale Jan Homburg 2,3, and Todd R. Young 4 1 Department of Mathematical, Information & Computer Sciences, Point Loma Nazarene University, 3900 Lomaland

More information

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics

Spurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics

More information

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F :

= F ( x; µ) (1) where x is a 2-dimensional vector, µ is a parameter, and F : 1 Bifurcations Richard Bertram Department of Mathematics and Programs in Neuroscience and Molecular Biophysics Florida State University Tallahassee, Florida 32306 A bifurcation is a qualitative change

More information

Self-Excited Vibration

Self-Excited Vibration Wenjing Ding Self-Excited Vibration Theory, Paradigms, and Research Methods With 228 figures Ö Springer Contents Chapter 1 Introduction 1 1.1 Main Features of Self-Excited Vibration 1 1.1.1 Natural Vibration

More information

Introduction to Group Theory

Introduction to Group Theory Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)

More information

A Model of Evolutionary Dynamics with Quasiperiodic Forcing

A Model of Evolutionary Dynamics with Quasiperiodic Forcing paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth

More information

Generalized Shifted Inverse Iterations on Grassmann Manifolds 1

Generalized Shifted Inverse Iterations on Grassmann Manifolds 1 Proceedings of the Sixteenth International Symposium on Mathematical Networks and Systems (MTNS 2004), Leuven, Belgium Generalized Shifted Inverse Iterations on Grassmann Manifolds 1 J. Jordan α, P.-A.

More information

Dynamical Systems. Pierre N.V. Tu. An Introduction with Applications in Economics and Biology Second Revised and Enlarged Edition.

Dynamical 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 Springer-Verlag Berlin Heidelberg New York London Paris

More information

Canards at Folded Nodes

Canards at Folded Nodes Canards at Folded Nodes John Guckenheimer and Radu Haiduc Mathematics Department, Ithaca, NY 14853 For Yulij Il yashenko with admiration and affection on the occasion of his 60th birthday March 18, 2003

More information

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Boltzmanngasse 9. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Boltzmanngasse 9 Institute for Mathematical Physics A-1090 Wien, Austria Common Homoclinic Points of Commuting Toral Automorphisms Anthony Manning Klaus Schmidt

More information

Torus Doubling Cascade in Problems with Symmetries

Torus Doubling Cascade in Problems with Symmetries Proceedings of Institute of Mathematics of NAS of Ukraine 4, Vol., Part 3, 11 17 Torus Doubling Cascade in Problems with Symmetries Faridon AMDJADI School of Computing and Mathematical Sciences, Glasgow

More information

Lecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)

Lecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition) Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries

More information

Lie Groups for 2D and 3D Transformations

Lie Groups for 2D and 3D Transformations Lie Groups for 2D and 3D Transformations Ethan Eade Updated May 20, 2017 * 1 Introduction This document derives useful formulae for working with the Lie groups that represent transformations in 2D and

More information

Singular Perturbation Theory for Homoclinic Orbits in a Class of Near-Integrable Hamiltonian Systems

Singular Perturbation Theory for Homoclinic Orbits in a Class of Near-Integrable Hamiltonian Systems Singular Perturbation Theory for Homoclinic Orbits in a Class of Near-Integrable Hamiltonian Systems Gregor Kovačič Mathematical Sciences Department Rensselaer Polytechnic Institute Troy, NY 12180 Abstract

More information

In Arnold s Mathematical Methods of Classical Mechanics (1), it

In Arnold s Mathematical Methods of Classical Mechanics (1), it Near strongly resonant periodic orbits in a Hamiltonian system Vassili Gelfreich* Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom Communicated by John N. Mather, Princeton

More information

IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 8, AUGUST Control Bifurcations

IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 8, AUGUST Control Bifurcations IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 49, NO. 8, AUGUST 2004 1231 Control Bifurcations Arthur J. Krener, Fellow, IEEE, Wei Kang, Senior Member, IEEE, Dong Eui Chang Abstract A parametrized nonlinear

More information

ax 2 + bx + c = 0 where

ax 2 + bx + c = 0 where Chapter P Prerequisites Section P.1 Real Numbers Real numbers The set of numbers formed by joining the set of rational numbers and the set of irrational numbers. Real number line A line used to graphically

More information

Congurations of periodic orbits for equations with delayed positive feedback

Congurations of periodic orbits for equations with delayed positive feedback Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics

More information

A matrix over a field F is a rectangular array of elements from F. The symbol

A matrix over a field F is a rectangular array of elements from F. The symbol Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted

More information

Lectures on Dynamical Systems. Anatoly Neishtadt

Lectures on Dynamical Systems. Anatoly Neishtadt Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances

More information

Available online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics

Available online at   ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing

More information

SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in

SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING A. DELSHAMS, V. GELFREICH, A. JORBA AND T.M. SEARA At the end of the last century, H. Poincare [7] discovered the phenomenon of separatrices splitting,

More information

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria

The Erwin Schrodinger International Pasteurgasse 6/7. Institute for Mathematical Physics A-1090 Wien, Austria ESI The Erwin Schrodinger International Pasteurgasse 6/7 Institute for Mathematical Physics A-090 Wien, Austria On the Rigid Body with Two Linear Controls Mircea Craioveanu Mircea Puta Vienna, Preprint

More information

Chapter 5 Eigenvalues and Eigenvectors

Chapter 5 Eigenvalues and Eigenvectors Chapter 5 Eigenvalues and Eigenvectors Outline 5.1 Eigenvalues and Eigenvectors 5.2 Diagonalization 5.3 Complex Vector Spaces 2 5.1 Eigenvalues and Eigenvectors Eigenvalue and Eigenvector If A is a n n

More information

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,

More information

Homoclinic bifurcations in Chua s circuit

Homoclinic bifurcations in Chua s circuit Physica A 262 (1999) 144 152 Homoclinic bifurcations in Chua s circuit Sandra Kahan, Anibal C. Sicardi-Schino Instituto de Fsica, Universidad de la Republica, C.C. 30, C.P. 11 000, Montevideo, Uruguay

More information

STABILITY OF INVARIANT SUBSPACES OF COMMUTING MATRICES We obtain some further results for pairs of commuting matrices. We show that a pair of commutin

STABILITY OF INVARIANT SUBSPACES OF COMMUTING MATRICES We obtain some further results for pairs of commuting matrices. We show that a pair of commutin On the stability of invariant subspaces of commuting matrices Tomaz Kosir and Bor Plestenjak September 18, 001 Abstract We study the stability of (joint) invariant subspaces of a nite set of commuting

More information

1 Matrices and Systems of Linear Equations

1 Matrices and Systems of Linear Equations Linear Algebra (part ) : Matrices and Systems of Linear Equations (by Evan Dummit, 207, v 260) Contents Matrices and Systems of Linear Equations Systems of Linear Equations Elimination, Matrix Formulation

More information

Im(v2) Im(v3) Re(v2)

Im(v2) Im(v3) Re(v2) . (a)..1. Im(v1) Im(v2) Im(v3)....... Re(v1) Re(v2).1.1.1 Re(v3) (b) y x Figure 24: (a) Temporal evolution of v 1, v 2 and v 3 for Fluorinert/silicone oil, Case (b) of Table, and 2 =,3:2. (b) Spatial evolution

More information

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS

PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS PROOF OF TWO MATRIX THEOREMS VIA TRIANGULAR FACTORIZATIONS ROY MATHIAS Abstract. We present elementary proofs of the Cauchy-Binet Theorem on determinants and of the fact that the eigenvalues of a matrix

More information

Lecture 1. Toric Varieties: Basics

Lecture 1. Toric Varieties: Basics Lecture 1. Toric Varieties: Basics Taras Panov Lomonosov Moscow State University Summer School Current Developments in Geometry Novosibirsk, 27 August1 September 2018 Taras Panov (Moscow University) Lecture

More information

Contents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces

Contents. 2.1 Vectors in R n. Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v. 2.50) 2 Vector Spaces Linear Algebra (part 2) : Vector Spaces (by Evan Dummit, 2017, v 250) Contents 2 Vector Spaces 1 21 Vectors in R n 1 22 The Formal Denition of a Vector Space 4 23 Subspaces 6 24 Linear Combinations and

More information