OSCILLATORS WITH DELAY COUPLING
|
|
- Dale Ball
- 1 years ago
- Views:
Transcription
1 Proceedings of DETC' ASME Design Engineering Technical Conferences September 1-15, 1999, Las Vegas, Nevada, USA D E T C 9 9 / V IB BIFURCATIONS IN THE DYNAMICS OF TWO COUPLED VAN DER POL OSCILLATORS WITH DELAY COUPLING Stephen Wirkus Center for Applied Mathematics Cornell University Ithaca, NY Richard Rand Department of Theoretical and Applied Mechanics Cornell University Ithaca, NY ABSTRACT We investigate the dynamics of a system of two van der Pol oscillators with delayed velocity coupling. We use the method of averaging to reduce the problem to the study of a slow-ow in three dimensions. We study the steady state solutions of this slow ow, with special attention given to the bifurcations accompanying their change in number and stability. Our interest in this system is due to its relevance to coupled laser oscillators. INTRODUCTION This work is concerned with the mutual interaction of two limit cycle oscillators. It is motivated by applications to laser dynamics (York and Compton, 199), (York 1993), (Lynch and York, 1995), (Lynch, 1995). When two lasers operate physically close to one another, the light output of each may aect the behavior of the other. Since the frequencies of some lasers are in the 10 GHz range, the time for light totravel from one laser to the other, a distance of the order of cm, represents a substantial portion of the period of the uncoupled laser oscillator. This immediately leads to the inclusion of delay eects in the coupling terms. In this paper we investigate the dynamics of two weakly coupled van der Pol oscillators in which thecoupling terms have time delay. We use the method of averaging to obtain an approximate simplied system of three slow ow equations and then investigate the stability and bifurcation of their equilibria (corresponding to periodic motions in the original system). This work is an extension of a previous paper by the present authors (Wirkus and Rand, 1997) and is also related to previous studies of coupled van derpol oscillators in which the coupling terms omitted delay eects (Rand and Holmes, 1980), (Storti and Rand, 198), (Storti and Rand, 1986),(Chakraborty and Rand, 1988). DERIVATION OFSLOWFLOWEQUATIONS We investigate two van der Pol oscillators with delay coupling (York 1993),(Lynch and York, 1995): x 1 + x 1 0 (10 x 1) _x 1 = _x (t 0 ); (1) x + x 0 (10 x ) _x = _x 1 (t 0 ): () where is a coupling parameter, is the delay time, and where <<1. When = 0, the system reduces to x i + x i = 0;i = 1;, with solution: x i = R i cos(t + i ); _x i = 0 R i sin(t + i ): (3) For >0, we look for a solution in the form (3) but treat R i and i as time dependent. Variation of parameters gives the following eqs. on R i and i : _R i = 0 sin(t + i ) F i ; R i _ i = 0 cos(t + i ) F i (4) where F i = (10 x i )_x i + _x j (t 0 );i;j =1; : For small we use the method of averaging, replacing the right hand sides of (4) by their averages over one period of the =0 1 Copyright c 1999 by ASME
2 system: Setting = 1 0,we then obtain in which _R i 0 1 R i _ i 0 1 Z 0 Z 0 sin(t + i ) F i dt; (5) cos(t + i ) F i dt; (6) _R 1 = [R 1 (1 0 R 1 4 )+R cos( + )]; (1) _R = [R (1 0 R 4 )+R 1 cos( 0 )]; (13) _ = [0 R R 1 sin( + ) 0 R 1 R sin( 0 )]: (14) F i =[10 R i cos (t + i )] [0R i sin(t + i )] + [0 ~ R j sin(t 0 + ~ j )]; (7) where ~ R = R(t 0 ) and ~ = (t 0 ). Evaluating the integrals H cos(t + i ) F i dt and H sin(t + i ) F i dt gives (i 6= j) _R i = R i (1 0 R i 4 )+ ~R j cos( i 0 ~ j + ); (8) These are the slow-ow equations which we will examine in the remainder of the paper. Note that R 1 and R are nonnegative and the vector eld associated with eqs.(1)-(14) is periodic in. Thus the phase space is R + R + S 1. The slow-ow isinvariant under the three transformations: (R 1 ; R ; ) 7! (R ; R 1 ; 0) (15) 7! + ; 7! 0 (16) R i _ i = 0 ~R j sin( i 0 ~ j + ): (9) Eqs. (8)-(9), show that R _ i, _ i are O(). We nowtaylor expand R ~ i and ~ i : ~R i = R i (t 0 ) =R i (t) 0 _ R i (t)+ Ri (t)+111; (10) ~ i = i (t 0 ) = i (t) 0 _ i (t)++ i (t)+111: (11) 7! + ; 7! + ; cos 7! 0 cos : (17) Eqs.(16),(17) show that another invariance involving only parameters is: 7! 0; 7! + ; cos 7! 0 cos : (18) Eqs. (10),(11) indicate that we can replace Ri ~, ~ i by R i, i in eqs.(8),(9) since R _ i (t), _ i (t) and R i (t), i (t) in (10),(11) are O() and O( ) respectively, from eqs.(8),(9). This reduces an innite dimensional problem in functional analysis to a nite dimensional problem by assuming the product is small. (Note that we do not assume asmall delay.) This key step enables us to handle the original system, adierential delay eq., as a system of dierential eqs. (York and Compton, 199),(York 1993),(Lynch and York, 1995),(Lynch, 1995). Note that if termsofo( )were retained in eqs. (10),(11), then the resulting dierential equations would be of second order. Thus extending the expansion to higher order in has the unusual eect of profoundly changing the nature of the approximate system to be solved. Nevertheless the O() truncation studied in this work is valid for small values of, as demonstrated by comparison of the subsequent slow ow analysis with numerical integration of the original dierential delay eqs. (1),(). Eq.(18) shows that we may assume >0 without loss of generality,since the phase ow for a negative value of is identical to that of the corresponding positive value of with the sign of cos reversed. STABILITY OFTHEIN-PHASE MODE Eqs.(1)-(14) possess the following equilibrium point which corresponds to the in-phase mode x 1 x in(1),(): R 1 = R = p 1+ cos ; =0;1+ cos >0: (19) Notice that eqs.(19) indicate that the amplitudes R i approach zeroas 1 + cos approaches zero, and that these Copyright c 1999 by ASME
3 amplitudes remain zero for 1 + cos 0. This result is known as amplitude death (Reddy, Sen,and Johnston, 1998), (Strogatz, 1998). Thus the in-phase mode is predicted to come into existence as we cross the curve = 01= cos in the cos 0 parameter plane. This bifurcation is accompanied by achange in stability of the trivial solution x 1 x 0. In order to show this we must return to the original dierential delay eqs.(1),(), since the slow ow(1)-(14) is singular for R 1 = R =0. Linearizing eqs.(1),() about the trivial solution, we may obtain the condition forachange in stability of x 1 x 0by assuming a solution of the form e i!t. Equating real and imaginary parts of the corresponding characteristic equation, we obtain: cos! = ! ; sin! =! : (0) In the small limit, the second of eqs.(0) gives that! =1, and then the rst of eqs.(0) gives = 01= cos,in agreement with the birth of the in-phase mode, cf.eq.(19). In order to determine the stability of the in-phase mode, we examine the eigenvalues of the Jacobian matrix of the right-hand side of eqs. (1)-(14) evaluated at the equilibrium (19). These satisfy the equation: 3 + (+4 cos ) + (1+5 cos + +4 cos ) +( cos + cos cos + 3 cos 3 )=0: (1) When = 0, eq.(1) gives conditions necessary for a saddlenode bifurcation. The resulting equation may be factored into three conditions: =0; () = 0 cos 1+cos ; (3) = 01 cos : (4) We can similarly set = i! to obtain conditions necessary for a Hopf bifurcation. Omitting extraneous roots, we obtain: = 0 1 3cos ; cos < 1 : (5) Simulation reveals that the limit cycle created in the Hopf bifurcation (5) is unstable. Eqs.()-(5) represent stability boundaries for the inphase mode x 1 x in the cos 0 parameter plane. See Fig.1. In order to check the stability results in Fig.1, we numerically simulated eqs.(1)-(14). The delay terms in eqs.(1)-(14) require that initial conditions be given for the time interval between 0 and 0. We handled this by allowing the two oscillators to run uncoupled from t = 0 to t = 0. This was accomplished by specifying x(0) and _x(0). This allowed us to obtain aninterval of initial conditions x(t) =8(t); 0 t 0. Then at time t =0,we started numerically integrating the coupled delay equations, which required that we take into account what happened time units ago. Here we useafourth order Runge-Kutta scheme with xed step size, appropriately modied to account for delay (Haireretal., 1987). In order to avoid unnecessarily long runs, we began \close" to the equilibrium point whose stability wewere investigating. For example, in the case of the in-phase mode, we took x 1 x,_x 1 _x. Apoint incos0 parameter space is saidtostable if the trajectory in four-dimensional phase space continues to spiral into the in-phase mode. In this way wewere able to conrm that the stability of the in-phase and out-of-phase modes agreed with the predictions of the analytical method shown in Fig. 1. STABILITY OFTHEOUT-OF-PHASE MODE In addition to the in-phase mode, eqs.(1)-(14) also possess an equilibrium which corresponds to the out-of-phase mode x 1 0x : R 1 = R = p 1 0 cos ; = ; 10 cos >0: (6) In order to determine the stability of this mode we could proceed in an analogous fashion to that used for the in-phase mode. However,a more direct approach is available to us due to the symmetry discussed in eqs.(16)-(18). The outof-phase mode (6) maps to the in-phase mode (19) under the transformation (17), 7! 0; 7! +. But since the entire phase ow is invariant under this transformation, the stability of (6) is seen to be identical to that of (19) with the parameter change cos 7! 0 cos. That is, the outof-phase mode has the same stability chart as the in-phase mode, reected about the -axis. See Fig.1. 3 Copyright c 1999 by ASME
4 OTHER PERIODIC MOTIONS Having investigated the stability of the slow ow equilibria corresponding to in-phase and out-of-phase modes, we now lookfor any other slow ow equilibria, each of which corresponds to a periodic motion in the original eqs.(1),(). Our task is to solve the three eqs.(1)-(14) for equilibrium values of R 1 ;R and. Wesolve (13) for R 1, substitute the result into (14) and solve for R. Call the result \eq.a". Then we solve (13) for R 1 and substitute the result into (1), giving a polynomial on R.We substitute \eq.a" into this polynomial, giving an equation with no R 1 or R in it. Algebraic and trigonometric simplication of the resulting equation gives: sin cos 4 +( sin 0 1) cos cos + cos 4 (1 + sin )=0: (7) Eq.(7) is a quadratic oncos.wemay obtain up to 4 real values for cos, corresponding to 8 values of. Of these 8, 4must be rejected because they correspond to negative R i values. The remaining 4 equilibria come in pairs which map to each other under the transformation (15). Bifurcations of these equilibria result from setting cos = 1 in eq.(7) since cos 1. This gives eq.(3) and its reection in the symmetry (17). Bifurcations also occur in eq.(7) if the discriminant vanishes. This results in the condition: = 1 8(1 0 cos ) ; cos < 1 3 : (8) Eq.(8) is displayed on Fig. along with the previously obtained bifurcation curves (3),(4),(5). Points P and Q in Fig. are dened as the points of intersection of curves (3),(5) and of curves (3),(8), respectively: P: cos = 0 1 p ;= p p 1 3 ; Q: cos = 0p ;= : (9) The total number of slow owequilibria in eqs.(1)-(14) depends on the parameters and, see Fig.. The maximum number is 6, consisting of the in-phase mode, the out-ofphase mode, and the 4 additional equilibria associated with eq.(7). Note that none of these additional equilibria can occur if the coupling is suciently large. We checked these results by numerically simulating eqs.(1)-(14). In doing so we discovered that the equilibria discussed in this section undergo Hopf bifurcations in the regions marked `4' and `6' infig.. In order to nd the location of these Hopfs in the parameter plane of Fig., we linearized the slow ow eqs.(1)-(14) about the equilibria given by eq.(7) and then set = i!in the associated eigenequation. The resulting equations were too complicated algebraically toallow us to obtain closed form conditions for the Hopfs (even using computer algebra). Nevertheless we were able to treat the resulting equations numerically, revealing that the Hopfs occur along a curve which is displayed in Fig.3. This curve is approximately given by the empirical equation: = 00:04566 cos 4 +0:8690 cos +0: ; 0:38095 cos 0:5: (30) As shown in Fig.3, the curve (30) reaches between point P and a point H on curve (8). The slow ow limit cycles which are born in these Hopfs correspond to quasiperiodic motions in the original eqs.(1),(). UNFOLDING POINT P In order to better understand the bifurcations which occur in the slow ow eqs.(1)-(14), we magnify the region of parameter space around point P. This involves (i) expanding in power series the parameters and about point P and the phase variables R 1, R and about the in-phase equilibrium point, then (ii) transforming the phase variables to local eigencoordinates (which involve a double zero eigenvalue at point P), then (iii) obtaining a power series approximation for the center manifold at P, and then (iv) transforming the resulting two-dimensional ow to normal form. This is a very complicated calculation which involved a lot of computer algebra. We give only the nal form of the ow on the center manifold: _y = z _z = y 3 ( p p )+y z( p p ) 16 + y(0 1 p p ) 6 + z( 3 p p ): (31) 4 Copyright c 1999 by ASME
5 where and are dened by the equations (cf.eq.(9)): = arccos 0 p 1 p + ; = + : (3) 3 Analysis of eq.(31) shows that there are 6 bifurcation curves which emanate from point P. See Fig.4 where these 6 curves are labeled A,B,C,D,E,F. Their description and approximate equations are found to be: p A,C: pitchfork, eq.(3): = p 54 (33) p p B: subcritical Hopf, eq.(5): = (34) D: supercritical Hopf, eq.(30): = p 6 p (35) are stable for values of cos close to zero, i.e., for delays of about 1 of the uncoupled period of the oscillators, see Fig.1. 4 In applications where it is desirable to operate the system at the in-phase mode, the presence of a stable out-of-phase mode could bethesource of unreliable behavior. The elaborate bifurcation sequences which were found to be present in this system occur over a very small region of parameter space, and hence are not likely to be exhibited by an associated physical system. Nevertheless these bifurcations are important because they explain how dynamical transitions occur between the larger regions of parameter space. The approximation used in this work is based on the assumption that << 1, and hence is not expected to give good results when this inequality is violated. In particular, the parameter always appears in the form cos in the slow ow(1)-(14), leading to the prediction that the dynamical behavior exhibits a periodic dependence on. Recent research work which does not assume that is small has shown that this prediction is only approximately true (Moon and Johnson, 1998), (Yeung and Strogatz, 1999). E: symmetry-breaking bifurcation: = 17p p (36) F: limit cycle fold : = 0:6606: (37) In curve E, the symmetry-breaking bifurcation, a pair of stable limit cycles is replaced by a single larger stable limit cycle. In curve F, the limit cycle fold, a pair oflimit cycles, one stable and one unstable, coalesce. It will be noted that the algebraic form of curve F involves an approximate decimal, in contrast to the expressions for the other bifurcation curves. In fact a closed form expression for this constant is known in terms of elliptic integrals, but is omitted here for brevity. The details by which eqs.(33)-(37) are obtained, which involve Melnikov and elliptic integrals, are omitted here, but will be included in the doctoral thesis of the rst author. Fig.5 showsaschematic of the bifurcation sequence obtained by moving through the parameter space near point P. DISCUSSION Perhaps the most signicant result of this work is the conclusion that both the in-phase and out-of-phase modes ACKNOWLEDGMENT The authors wish to thank Professor Rick Compton for proposing this problem, and Professor Steven Strogatz for helpful suggestions. REFERENCES Chakraborty,T. and Rand,R.H. The transition from phase locking to drift in a system of two weakly coupled van der Pol oscillators, Int. J. Nonlinear Mechanics 3: , Hairer,E., Norsett,S.P., and Wanner,G. Solving Ordinary Dierential Equations I: Nonsti Problems, Springer- Verlag, Lynch,J.J. Analysis and design of systems of coupled microwave oscillators, Ph.D.thesis, Dept.Electricaland Computer Engineering, University of California at Santa Barbara, Lynch,J.J. and York,R.A. Stability of mode locked states of coupled oscillator arrays, IEEE Trans. on Circuits and Systems 4: , Moon,F.C. and Johnson,M.A. Nonlinear dynamics and chaos in manufacturing processes, in Dynamics and Chaos in Manufacturing Processes, ed.f.c.moon, Wiley, pp.3-3, Rand,R.H. and Holmes,P.J. Bifurcation of periodic mo- 5 Copyright c 1999 by ASME
6 tions in two weakly coupled van der Pol oscillators, Int. J. Nonlinear Mechanics 15: , Reddy,D.V.R., Sen,A., and Johnston,G.L. Time delay induced death in coupled limit cycle oscillators, Physical Review Letters 80: , Storti,D.W. and Rand,R.H. Dynamics of two strongly coupled van der Pol oscillators, Int. J. Nonlinear Mechanics 17:143-15, 198. Storti,D.W. and Rand,R.H. Dynamics of two strongly coupled relaxation oscillators, SIAM J. Applied Math. 46:56-67, Strogatz,S.H. Death by delay, Nature 394: , Wirkus,S. and Rand,R. Dynamics of two coupled van der Pol oscillators with delay coupling, Proceedings of DETC'97, 1997 ASME Design Engineering Technical Conferences, Sept.14-17,1997, Sacramento, CA, paper no.detc97/vib Yeung,M.K.S. and Strogatz,S. Time delay in the Kuramoto model of coupled oscillators, Physical Review Letters 8: , York,R.A. Nonlinear analysis of phase relationships in quasi-optical oscillator arrays, IEEE Trans. on Microwave Theory and Tech. 41: , York,R.A. and Compton,R.C. Experimental observation and simulation of mode-locking phenomena in coupledoscillator arrays, J.Appl.Phys. 71: , Copyright c 1999 by ASME
7
8
9
10
NONLINEAR DYNAMICS OF THE BOMBARDIER BEETLE. Center for Applied Mathematics. Ithaca,NY to seconds. this increased pressure. results.
Proceedings of DETC'99 1999 ASME Design Engineering Technical Conferences September 12-15, 1999, Las Vegas, Nevada, USA D E T C 9 9 / V IB -8 0 1 1 NONLINEAR DYNAMICS OF THE BOMBARDIER BEETLE Richard H.Rand
Three ways of treating a linear delay differential equation
Proceedings of the 5th International Conference on Nonlinear Dynamics ND-KhPI2016 September 27-30, 2016, Kharkov, Ukraine Three ways of treating a linear delay differential equation Si Mohamed Sah 1 *,
THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS
THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the
2:1 Resonance in the delayed nonlinear Mathieu equation
Nonlinear Dyn 007) 50:31 35 DOI 10.1007/s11071-006-916-5 ORIGINAL ARTICLE :1 Resonance in the delayed nonlinear Mathieu equation Tina M. Morrison Richard H. Rand Received: 9 March 006 / Accepted: 19 September
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
Dynamics of a mass-spring-pendulum system with vastly different frequencies
Dynamics of a mass-spring-pendulum system with vastly different frequencies Hiba Sheheitli, hs497@cornell.edu Richard H. Rand, rhr2@cornell.edu Cornell University, Ithaca, NY, USA Abstract. We investigate
arxiv: v1 [nlin.ao] 19 May 2017
Feature-rich bifurcations in a simple electronic circuit Debdipta Goswami 1, and Subhankar Ray 2, 1 Department of Electrical and Computer Engineering, University of Maryland, College Park, MD 20742, USA
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
USING COUPLED OSCILLATORS TO MODEL THE SINO-ATRIAL NODE IN THE HEART
USING COUPLED OSCILLATORS TO MODEL THE SINO-ATRIAL NODE IN THE HEART A Thesis Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree
2:2:1 Resonance in the Quasiperiodic Mathieu Equation
Nonlinear Dynamics 31: 367 374, 003. 003 Kluwer Academic Publishers. Printed in the Netherlands. ::1 Resonance in the Quasiperiodic Mathieu Equation RICHARD RAND Department of Theoretical and Applied Mechanics,
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
1. < 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
RICH VARIETY OF BIFURCATIONS AND CHAOS IN A VARIANT OF MURALI LAKSHMANAN CHUA CIRCUIT
International Journal of Bifurcation and Chaos, Vol. 1, No. 7 (2) 1781 1785 c World Scientific Publishing Company RICH VARIETY O BIURCATIONS AND CHAOS IN A VARIANT O MURALI LAKSHMANAN CHUA CIRCUIT K. THAMILMARAN
as Hopf Bifurcations in Time-Delay Systems with Band-limited Feedback
as Hopf Bifurcations in Time-Delay Systems with Band-limited Feedback Lucas Illing and Daniel J. Gauthier Department of Physics Center for Nonlinear and Complex Systems Duke University, North Carolina
Lecture 5. Outline: Limit Cycles. Definition and examples How to rule out limit cycles. Poincare-Bendixson theorem Hopf bifurcations Poincare maps
Lecture 5 Outline: Limit Cycles Definition and examples How to rule out limit cycles Gradient systems Liapunov functions Dulacs criterion Poincare-Bendixson theorem Hopf bifurcations Poincare maps Limit
3 Stability and Lyapunov Functions
CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such
ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS
Volume 1, No. 8, August 2013 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS Tarini
DRIVEN and COUPLED OSCILLATORS. I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela
DRIVEN and COUPLED OSCILLATORS I Parametric forcing The pendulum in 1:2 resonance Santiago de Compostela II Coupled oscillators Resonance tongues Huygens s synchronisation III Coupled cell system with
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
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
Queues with Choice via Delay Differential Equations
Queues with Choice via Delay Differential Equations Jamol Pender School of Operations Research and Information Engineering Cornell University 228 Rhodes Hall, Ithaca, NY 14853 jjp274@cornell.edu Richard
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
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
Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis
Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis Mimila Prost O. Collado J. Automatic Control Department, CINVESTAV IPN, A.P. 4 74 Mexico, D.F. ( e mail:
Duffing-Type Oscillators with Amplitude-Independent Period
Duffing-Type Oscillators with Amplitude-Independent Period Ivana N. Kovacic and Richard H. Rand Abstract Nonlinear oscillators with hardening and softening cubic Duffing nonlinearity are considered. Such
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
2D-Volterra-Lotka Modeling For 2 Species
Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose
Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 65 71 Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation Faridon AMDJADI Department
Some Dynamical Behaviors In Lorenz Model
International Journal Of Computational Engineering Research (ijceronline.com) Vol. Issue. 7 Some Dynamical Behaviors In Lorenz Model Dr. Nabajyoti Das Assistant Professor, Department of Mathematics, Jawaharlal
7 Pendulum. Part II: More complicated situations
MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics
Chapter 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 =
A Study of Nonlinear Dynamics in Mathematical Biology
UNF Digital Commons UNF Theses and Dissertations Student Scholarship 2013 A Study of Nonlinear Dynamics in Mathematical Biology Joseph Ferrara University of North Florida Suggested Citation Ferrara, Joseph,
Problem Description The problem we consider is stabilization of a single-input multiple-state system with simultaneous magnitude and rate saturations,
SEMI-GLOBAL RESULTS ON STABILIZATION OF LINEAR SYSTEMS WITH INPUT RATE AND MAGNITUDE SATURATIONS Trygve Lauvdal and Thor I. Fossen y Norwegian University of Science and Technology, N-7 Trondheim, NORWAY.
The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry Friction
Send Orders for Reprints to reprints@benthamscience.ae 308 The Open Mechanical Engineering Journal, 2014, 8, 308-313 Open Access The Stick-Slip Vibration and Bifurcation of a Vibro-Impact System with Dry
Abstract. In this paper we study the performance of a symplectic numerical integrator
Geometric Numerical Integration Applied to The Elastic Pendulum at Higher Order Resonance J. M. Tuwankotta y G.R.W. Quispel E-mail: tuwankotta@math.uu.nl Email: R.Quispel@latrobe.edu.au Mathematisch Instituut,
CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD
Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen
Ridge analysis of mixture response surfaces
Statistics & Probability Letters 48 (2000) 3 40 Ridge analysis of mixture response surfaces Norman R. Draper a;, Friedrich Pukelsheim b a Department of Statistics, University of Wisconsin, 20 West Dayton
ME 680- Spring Geometrical Analysis of 1-D Dynamical Systems
ME 680- Spring 2014 Geometrical Analysis of 1-D Dynamical Systems 1 Geometrical Analysis of 1-D Dynamical Systems Logistic equation: n = rn(1 n) velocity function Equilibria or fied points : initial conditions
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
1 Ordinary points and singular points
Math 70 honors, Fall, 008 Notes 8 More on series solutions, and an introduction to \orthogonal polynomials" Homework at end Revised, /4. Some changes and additions starting on page 7. Ordinary points and
Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical
ANZIAM J. 53 (EMAC211) pp.c511 C524, 212 C511 Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical M. R. Alharthi 1 T. R. Marchant 2 M. I.
Chapter 24 BIFURCATIONS
Chapter 24 BIFURCATIONS Abstract Keywords: Phase Portrait Fixed Point Saddle-Node Bifurcation Diagram Codimension-1 Hysteresis Hopf Bifurcation SNIC Page 1 24.1 Introduction In linear systems, responses
Period-doubling cascades of a Silnikov equation
Period-doubling cascades of a Silnikov equation Keying Guan and Beiye Feng Science College, Beijing Jiaotong University, Email: keying.guan@gmail.com Institute of Applied Mathematics, Academia Sinica,
POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS
POSITIVE REALNESS OF A TRANSFER FUNCTION NEITHER IMPLIES NOR IS IMPLIED BY THE EXTERNAL POSITIVITY OF THEIR ASSOCIATE REALIZATIONS Abstract This letter discusses the differences in-between positive realness
LIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES
LIMITATIONS OF THE METHOD OF MULTIPLE-TIME-SCALES Peter B. Kahn Department of Physics and Astronomy State University of New York Stony Brook, NY 794 and Yair Zarmi Department for Energy & Environmental
1 Introduction Travelling-wave solutions of parabolic equations on the real line arise in a variety of applications. An important issue is their stabi
Essential instability of pulses, and bifurcations to modulated travelling waves Bjorn Sandstede Department of Mathematics The Ohio State University 231 West 18th Avenue Columbus, OH 4321, USA Arnd Scheel
TECHNICAL NOTE: PREDICTION OF THE THRESHOLD OF GLOBAL SURF-RIDING BY AN EXTENDED MELNIKOV METHOD
10 th International Conference 441 TECHNICAL NOTE: PREDICTION OF THE THRESHOLD OF GLOBAL SURF-RIDING BY AN EXTENDED MELNIKOV METHOD Wan Wu, Virginia Polytechnic Institute and State University, wanwu@vt.edu
Kink, singular soliton and periodic solutions to class of nonlinear equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic
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
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
expression that describes these corrections with the accuracy of the order of 4. frame usually connected with extragalactic objects.
RUSSIAN JOURNAL OF EARTH SCIENCES, English Translation, VOL, NO, DECEMBER 998 Russian Edition: JULY 998 On the eects of the inertia ellipsoid triaxiality in the theory of nutation S. M. Molodensky Joint
Chapter 3 HIGHER ORDER SLIDING MODES. L. FRIDMAN and A. LEVANT. Chihuahua Institute of Technology, Chihuahua, Mexico.
Chapter 3 HIGHER ORDER SLIDING MODES L. FRIDMAN and A. LEVANT Chihuahua Institute of Technology, Chihuahua, Mexico. Institute for Industrial Mathematics, Beer-Sheva, Israel. 3.1 Introduction One of the
Report E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache
Potsdam, August 006 Report E-Project Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory
Chaos in the Planar Two-Body Coulomb Problem with a Uniform Magnetic Field
Annual Review of Chaos Theory, Bifurcations and Dynamical Systems Vol. 3, (2013) 23-33, www.arctbds.com. Copyright (c) 2013 (ARCTBDS). ISSN 2253 0371. All Rights Reserved. Chaos in the Planar Two-Body
Modeling Low-Dimensional Submanifolds In Dynamical Systems
1 Clarkson University Modeling Low-Dimensional Submanifolds In Dynamical Systems A Dissertation by Chen Yao Department of Mathematics and Computer Science Submitted in partial fulllment of the requirements
Math 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
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
1 Introduction Self-replicating spots and pulses have been observed in excitable reaction-diusion systems [22, 17, 24, 23, 16, 9, 2, 3, 4, 25, 21, 19,
Slowly-modulated two pulse solutions and pulse splitting bifurcations Arjen Doelman Korteweg-deVries Instituut Universiteit van Amsterdam Plantage Muidergracht 24 1018TV Amsterdam, The Netherlands Wiktor
2.3. VECTOR SPACES 25
2.3. VECTOR SPACES 25 2.3 Vector Spaces MATH 294 FALL 982 PRELIM # 3a 2.3. Let C[, ] denote the space of continuous functions defined on the interval [,] (i.e. f(x) is a member of C[, ] if f(x) is continuous
Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS
Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS WALTER GALL, YING ZHOU, AND JOSEPH SALISBURY Department of Mathematics
Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type 0
Anais do CNMAC v.2 ISSN 1984-82X Characterization of the stability boundary of nonlinear autonomous dynamical systems in the presence of a saddle-node equilibrium point of type Fabíolo M. Amaral Departamento
SYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN
ELECTRONICS 7 19 1 September, Sozopol, BULGARIA SYNTHESIS OF SINEWAVE OSCILLATOR BASED ON THE MODIFIED VAN DER POL EQUATION USING MELNIKOV THEORY Zhivko Dimitrov Georgiev 1, Todor Georgiev Todorov, Emil
Monte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time
Monte Carlo simulations of harmonic and anharmonic oscillators in discrete Euclidean time DESY Summer Student Programme, 214 Ronnie Rodgers University of Oxford, United Kingdom Laura Raes University of
M445: Heat equation with sources
M5: Heat equation with sources David Gurarie I. On Fourier and Newton s cooling laws The Newton s law claims the temperature rate to be proportional to the di erence: d dt T = (T T ) () The Fourier law
Numerical techniques. Chapter Difference equations
Chapter 6 Numerical techniques The differential equations (4.61), (4.62) and (4.64), in combination with boundary conditions such as equations (4.65) (4.68), constitute a two point boundary value problem.
Stability Analysis of Uzawa-Lucas Endogenous Growth Model
Abstract: Stability Analysis of Uzawa-Lucas Endogenous Growth Model William A. Barnett* University of Kansas, Lawrence and Center for Financial Stability, NY City and Taniya Ghosh Indira Gandhi Institute
Claudio A. Ca~nizares. University ofwaterloo. E&CE Department. presented in Section III regarding the association
Panel Session: \Optimization Techniques in Voltage Collapse Analysis," IEEE/PES Summer Meeting, San Diego, July 4, 998. Applications of Optimization to Voltage Collapse Analysis Abstract This paper describes
Abteilung Theoretische Physik Universitat Ulm Semiclassical approximations and periodic orbit expansions in chaotic and mixed systems Habilitationsschrift vorgelegt an der Fakultat fur Naturwissenschaften
Our goal is to solve a general constant coecient linear second order. this way but that will not always happen). Once we have y 1, it will always
October 5 Relevant reading: Section 2.1, 2.2, 2.3 and 2.4 Our goal is to solve a general constant coecient linear second order ODE a d2 y dt + bdy + cy = g (t) 2 dt where a, b, c are constants and a 0.
An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems
An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems
A Producer-Consumer Model With Stoichiometry
A Producer-Consumer Model With Stoichiometry Plan B project toward the completion of the Master of Science degree in Mathematics at University of Minnesota Duluth Respectfully submitted by Laura Joan Zimmermann
A Finite Element Method for an Ill-Posed Problem. Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D Halle, Abstract
A Finite Element Method for an Ill-Posed Problem W. Lucht Martin-Luther-Universitat, Fachbereich Mathematik/Informatik,Postfach 8, D-699 Halle, Germany Abstract For an ill-posed problem which has its origin
Nonlinear Control Lecture 1: Introduction
Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation
Stability 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
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.
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
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS
DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University
I/O monotone dynamical systems. Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011
I/O monotone dynamical systems Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011 BEFORE: Santa Barbara, January 2003 Having handed to me a photocopied
Analysis of a lumped model of neocortex to study epileptiform ac
of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized
Nonlinear dynamics in the Cournot duopoly game with heterogeneous players
arxiv:nlin/0210035v1 [nlin.cd] 16 Oct 2002 Nonlinear dynamics in the Cournot duopoly game with heterogeneous players H. N. Agiza and A. A. Elsadany Department of Mathematics, Faculty of Science Mansoura
Xinghuo Yu. Faculty of Informatics and Communication. Yanxing Song. Zhenwei Cao. Faculty of Engineering and Physical Systems.
A Generalized OGY Method for Controlling Higher Order Chaotic Systems Xinghuo Yu Faculty of Informatics and Communication Central Queensland University, Rockhampton, Qld 472, Australia. Guanrong Chen Department
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
Secondary Math 3 Honors - Polynomial and Polynomial Functions Test Review
Name: Class: Date: Secondary Math 3 Honors - Polynomial and Polynomial Functions Test Review 1 Write 3x 2 ( 2x 2 5x 3 ) in standard form State whether the function is even, odd, or neither Show your work
Departement Elektrotechniek ESAT-SISTA/TR Dynamical System Prediction: a Lie algebraic approach for a novel. neural architecture 1
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1995-47 Dynamical System Prediction: a Lie algebraic approach for a novel neural architecture 1 Yves Moreau and Joos Vandewalle
OUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL
International Journal in Foundations of Computer Science & Technology (IJFCST),Vol., No., March 01 OUTPUT REGULATION OF RÖSSLER PROTOTYPE-4 CHAOTIC SYSTEM BY STATE FEEDBACK CONTROL Sundarapandian Vaidyanathan
Variational Embedded Solitons, And Traveling Wavetrains Generated By Generalized Hopf Bifurcations, In Some Nlpde Systems
University of Central Florida Electronic Theses and Dissertations Doctoral Dissertation (Open Access) Variational Embedded Solitons, And Traveling Wavetrains Generated By Generalized Hopf Bifurcations,
Follow links Class Use and other Permissions. For more information, send to:
COPYRIGHT NOTICE: Stephen L. Campbell & Richard Haberman: Introduction to Differential Equations with Dynamical Systems is published by Princeton University Press and copyrighted, 2008, by Princeton University
The Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
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),
Lectures on Periodic Orbits
Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete
An FKG equality with applications to random environments
Statistics & Probability Letters 46 (2000) 203 209 An FKG equality with applications to random environments Wei-Shih Yang a, David Klein b; a Department of Mathematics, Temple University, Philadelphia,
Melnikov s Method Applied to a Multi-DOF Ship Model
Proceedings of the th International Ship Staility Workshop Melnikov s Method Applied to a Multi-DOF Ship Model Wan Wu, Leigh S. McCue, Department of Aerospace and Ocean Engineering, Virginia Polytechnic
1 Introduction In recent years many eorts have been directed towards a deeper understanding of the eects caused by discretizing a continuous dynamical
Invariant Manifolds for Non-autonomous Systems with Application to One-step Methods Y.-K. Zou y W.-J. Beyn Abstract In this paper we study the existence of invariant manifolds for a special type of nonautonomous
Numerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
Stability of the Lagrange Points, L 4 and L 5
Stability of the Lagrange Points, L 4 and L 5 Thomas Greenspan January 7, 014 Abstract A proof of the stability of the non collinear Lagrange Points, L 4 and L 5. We will start by covering the basics of
Taylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
GlobalStabilityofStificrsImpactonMemeTransmissionModel. Global Stability of Stificrs Impact on Meme Transmission Model
Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 14 Issue 6 Version 1. Year 14 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global
Simple 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.
Delay Kaldor-Kalecki Model Revisited. September 2014
Discussion Paper No.234 Delay Kaldor-Kalecki Model Revisited Akio Matsumoto Chuo University Ferenc Szidarovszky University of Pécs September 2014 INSTITUTE OF ECONOMIC RESEARCH Chuo University Tokyo, Japan
Chemical Kinetics and the Rössler System. 1 Introduction. 2 The NH 3 - HCl reaction. Dynamics at the Horsetooth Volume 2, 2010.
Dynamics at the Horsetooth Volume 2, 21. Chemical Kinetics and the Rössler System Department of Mathematics Colorado State University shinn@math.colostate.edu Report submitted to Prof. P. Shipman for Math