DOUBLE HOPF BIFURCATION IN DELAYED VAN DER POL DUFFING EQUATION

Size: px
Start display at page:

Download "DOUBLE HOPF BIFURCATION IN DELAYED VAN DER POL DUFFING EQUATION"

Transcription

1 International Journal of Bifurcation and Chaos, Vol. 3, No. (3) 354 (5 pages) c World Scientific Publishing Company DOI:.4/S DOUBLE HOPF BIFURCATION IN DELAYED VAN DER POL DUFFING EQUATION YUTING DING,, WEIHUA JIANG and PEI YU, Department of Mathematics, Harbin Institute of Technology, Harbin, 5, P. R. China Department of Applied Mathematics, Western University, London, Ontario, Canada N6A 5B7 pyu@uwo.ca Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. Received October 4, In this paper, we study dynamics in delayed van der Pol Duffing equation, with particular attention focused on nonresonant double Hopf bifurcation. Both multiple time scales and center manifold reduction methods are applied to obtain the normal forms near a double Hopf critical point. A comparison between these two methods is given to show their equivalence. Bifurcations are classified in a two-dimensional parameter space near the critical point. Numerical simulations are presented to demonstrate the applicability of the theoretical results. Keywords: Van der Pol Duffing equation; double Hopf bifurcation; normal form; multiple time scales; center manifold reduction.. Introduction Recently, much attention has been focused on the study of high-codimensional bifurcations, since they may exhibit complex dynamical behaviors such as quasi-periodic solutions and chaos. Center manifold theory and normal form theory are usually applied to analyze stability and bifurcation of dynamical systems (e.g. see [Guckenheimer & Holmes, 99; Kuznetsov, 4; Wiggins, 99]). Especially, in delayed differential equations (DDEs), due to time delay, higher-codimension singularities may occur more frequently than that in ordinary differential equations (ODEs). Even for a scalar DDE, it is possible to have bifurcation of limit cycles and even higher-codimension bifurcation phenomena, (e.g. see [Choi & LeBlanc, 6; Hale, 997] and references therein). Several methods have been developed for computing the normal forms of differential equations, for example, multiple time scales (MTS) [Nayfeh, 973, 98] and center manifold reduction (CMR) [Guckenheimer & Holmes, 99; Hale, 997; Hassard et al., 98; Wiggins, 99]. For a dynamical system described by ODEs, multiple time scales method is systematic and can be directly applied to the original nonlinear system, without the application of center manifold theory [Yu, 998, ; Zheng & Wang, ]. In fact, this approach combines the two steps involved in using center manifold theory and normal form theory into one unified step to obtain the normal form and nonlinear transformation simultaneously. Based on multiple time scales, Yu [998,, ] developed Maple programs for computing the normal forms of Hopf bifurcation and other singularities. These programs can be automatically executed by using a computer algebra system for a given ODE system. The basic idea of the center manifold theory is employing successive coordinate transformations to systematically construct a simpler system which has Author for correspondence 354-

2 Y. Ding et al. Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. less dimension compared to the original system, and thus greatly simplifying the dynamical analysis of the system [Guckenheimer & Holmes, 99; Kuznetsov, 4; Wiggins, 99]. The multiple time scales can also be directly applied to DDEs without the application of center manifold theory [Das & Chatterjee, ; Hu & Wang, 9; Nayfeh, 8]. In contrast, by the center manifold theory, one needs to first change the retarded equations to operator differential equations, and then decompose the solution space of their linearized form into stable and center manifolds. Next, with the adjoint of the operator equations, one computes the center manifold by projecting the whole space to the center manifold, finally calculating the normal form on the center manifold [Faria & Magalhaes, 995; Jiang & Wang, ; Ma et al., 8; Revel et al., ; Wei, 7; Xu & Yu, 4; Yu et al., ; Yuan & Wei, 6]. In the research of nonlinear dynamical systems, the van der Pol Duffing equation is one of the most intensively studied equations. This celebrated equation originally was a model for an electrical circuit with a triode valve, and was extensively studied as a host of a rich class of dynamical behaviors, including relaxation oscillations, quasi-periodicity, elementary bifurcations and chaos [Kuznetsov, 4]. It is well known that the limit cycle oscillations with strong stability property are important in applications, hence, being able to modify their behaviors through feedback is a question of interest. On the other hand, most practical implementations have inherent or feedback delays, the presence of which leads to the appearance of complex phenomena in the autonomous van der Pol Duffing equation, such as Hopf pitchfork bifurcation, double Hopf bifurcation and Bogdanov Takens singularity, etc. [Jiang & Yuan, 7; Jiang & Wei, 8; Wang & Jiang, ; Xu & Chung, 3; Xu et al., 7]. In this paper, we consider the nonresonant double Hopf bifurcation in the following van der Pol equation with delayed feedback: ẍ(t)+ω x(t) [b γx (t)]ẋ(t)+βx 3 (t) = A[x(t τ) x(t)], () where ω, b, γ, β are positive real constants, A is feedback strength, and τ is time delay. Xu ω, = et al. [7] employed the perturbation-incremental scheme (PIS) to investigate the weakly resonant double Hopf bifurcation and dynamics of system (). In this paper, we will study nonresonant double Hopf bifurcation of () by using two methods to compute the normal form of the system, namely, the multiple time scales and a combination of the method of normal forms and the center manifold theorem. A comparison between the normal forms shows that the two methods are equivalent. Furthermore, we carry out the bifurcation analysis and numerical simulations. We will show that there exist a stable fixed point, stable periodic solutions and stable quasi-periodic solutions in the neighborhood of the double Hopf critical point. The rest of the paper is organized as follows. In Sec., we consider the existence of double Hopf bifurcation in the delayed van der Pol system (), and use two methods to derive the normal form associated with double Hopf bifurcation. Then, bifurcation analysis and numerical simulations are presented in Sec. 3. Finally, the conclusion isdrawninsec.4.. Analytical Study In this section, we consider the van der Pol Duffing equation with delayed feedback, described by (), and use the MTS and CMR methods to derive the normal form of the system... System formulation The characteristic equation of the linearized equation of (), evaluated at the trivial equilibrium x =ẋ =,isgivenby: λ bλ + ω + A Ae λτ =. () To find possible periodic solutions, which may bifurcate from a Hopf or double Hopf critical point, let λ =iω (i =,ω >) be a root of (). Substituting the root into () and separating the real and imaginary parts yields { ω + A ω = A cos(ωτ), (3) bω = A sin(ωτ), from which we obtain A +ω b ± (b ω A) 4(ω 4 +Aω ), 354-

3 Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. under the assumption: { A + ω >, (A b ) 4ω b >. Further, it follows from (3) that τ (j), = ω, ω, (4) [ ( arccos + ω ) ] ω, +jπ, A for A>, [ ( π arccos + ω ) ] ω, +jπ, A for A<, where j =,,,... Thus, a possible double Hopf bifurcation occurs when two such families of surfaces intersect, with τ c = τ (j) = τ (l),wherej, l =,,,... The equality τ c = τ (j) = τ (l) implies that the linearized system on the trivial equilibrium has two pairs of purely imaginary eigenvalues ±iω and ±iω. Assume ω : ω = k : k, then a possible double Hopf bifurcation with the ratio k : k appears. When k, k Z +, it is called a k : k resonant double Hopf bifurcation; otherwise, it is called a nonresonant double Hopf bifurcation. In this paper, we only consider the nonresonant double Hopf bifurcation, for which ω = ω = k.. Multiple time scales ω A + ω k, (5) k ω A + ω k. We treat the feedback strength A and the delay τ as two bifurcation parameters. From τ c = τ (j) = τ (l), l, j =,,,..., we get the critical value A c. Suppose system () undergoes a double Hopf bifurcation from the trivial equilibrium at the critical point: A = A c, τ = τ c. Further, by the MTS, the solution of () is assumed to take the form: x(t) =ɛx (T,T,T,...)+ɛ x (T,T,T,...) + ɛ 3 x 3 (T,T,T,...)+, (6) where T k = ɛ k t, k =,,,...The derivative with respect to t is now transformed into d dt = + ɛ + ɛ + T T =D + ɛd + ɛ D +, where the differential operator D i = T i, i =,,,... We take perturbations as A = A c + ɛ A ɛ and τ = τ c + ɛτ ɛ in (), where A ɛ is called a detuning parameter [Nayfeh, 8]. To deal with the delayed terms, we expand x j (T τ c ɛτ ɛ,t ɛ(τ c +ɛτ ɛ ),T ɛ (τ c +ɛτ ɛ ),...)atx j (T τ c,t,t ) for j =,, 3,... Then, substituting the solutions with the multiple scales into () and balancing the coefficients of ɛ n (n =,, 3,...) yields a set of ordered linear differential equations. First, for the ɛ -order terms, we have D x +(A c + ω )x bd x A c x τc =, (7) where x = x (T,T,T ) and x τc = x (T τ c,t,t ). Since ±iω and ±iω are the eigenvalues of the linear part of (), the solution of (7) can beexpressedintheformof x (T,T,T ) = G (T,T )sin(ω T )+G (T,T )cos(ω T ) + G 3 (T,T )sin(ω T ) + G 4 (T,T )cos(ω T ). (8) Next, for the ɛ -order terms, we obtain D x +(A c + ω )x bd x A c x τc = D D x + bd x A c (τ ɛ D x τc + τ c D x τc ). (9) Substituting solution (8) into(9) and simplifying, we obtain the following equation: D x +(A c + ω )x bd x A c x τc + P cos(ω T )+P sin(ω T ) + P 3 cos(ω T )+P 4 sin(ω T )=, () where P i (i =,, 3, 4) are given as follows: G P =ω b G + A c τ ɛ ω [G cos(ω τ c )+G sin(ω τ c )] [ G + A c τ c cos(ω τ c ) G ] sin(ω τ c ), 354-3

4 Y. Ding et al. Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. P = A c τ ɛ ω [G sin(ω τ c ) G cos(ω τ c )] ω G b G + A c τ c [ G cos(ω τ c )+ G sin(ω τ c ) ], G 3 P 3 =ω b G 4 + A c τ ɛ ω [G 3 cos(ω τ c )+G 4 sin(ω τ c )] [ G4 + A c τ c cos(ω τ c ) G ] 3 sin(ω τ c ), P 4 = A c τ ɛ ω [G 3 sin(ω τ c ) G 4 cos(ω τ c )] ω G 4 b G 3 + A c τ c [ G3 cos(ω τ c )+ G 4 sin(ω τ c ) ]. To avoid occurrence of secular terms in the solution of (), the coefficients of cos(ω i T ), sin(ω i T )(i =, ) in () must be set to zero, i.e. P i =(i =,, 3, 4). Therefore, G, G, G 3 and G 4 are solved from the four linear equations in terms of G, G, G 3 and G 4. Then, Eq. () is reduced to D x +(A c + ω )x bd x A c x τc =, () where x = x (T,T,T )andx τc = x (T τ c, T,T ), and thus, the particular solution of () is x (T,T,T )=. () Further, for the ɛ 3 -order terms, we similarly obtain D x 3 +(A c + ω )x 3 bd x 3 A c x 3τc = D D x D D x D x A ɛ x γx D x βx 3 + A ɛx τc A c (τ ɛ D x τc + τ c D x τc + τ c D x τc + τ ɛ D x τc ) + b(d x +D x )+ A c (τ ɛd x τc + τ c D x τ c +τ ɛ τ c D D x τc ), (3) where x 3 = x 3 (T,T,T ) and x 3τc = x 3 (T τ c,t,t ). Substituting the solutions (8) and() into (3) and letting the coefficients of the terms which may generate secular terms in the solution equal to zero, yields the derivatives G T, G T, G 3 T and G 4 T expressed in terms of G, G, G 3 and G 4. The above procedure can in principle continue indefinitely (to any high order). Finally, the equations for Ġ, Ġ, Ġ 3 and Ġ4 are given by {Ġ = ɛd G + ɛ D G +, (4) Ġ = ɛd G + ɛ D G +, {Ġ3 = ɛd G 3 + ɛ D G 3 +, (5) Ġ 4 = ɛd G 4 + ɛ D G 4 +. Let G = R sin(θ (t)), G = R cos(θ (t)), G 3 = R sin(θ (t)) and G 4 = R cos(θ (t)). Substituting these expressions into (4) and(5), and truncating the equations at O(ɛ 3 ) yields the following normal form in polar coordinates: R = A cω [a cos(ω τ c )+b sin(ω τ c )] (a + b )τ R τ ɛ + [b cos(ω τ c ) b a sin(ω τ c )] c a + R A ɛ b + Q τ ɛ R 3βb + ω γa (a + b ) R R 3βb + ω γa 4(a + b ) R 3, R = A cω [a cos(ω τ c )+b sin(ω τ c )] (a + b )τ R τ ɛ + [b cos(ω τ c ) b a sin(ω τ c )] c a + R A ɛ b + Q τ ɛ R 3βb + ω γa (a + b ) R R 3βb + ω γa 4(a + b ) R 3, Θ = δ + H τ ɛ T + 3a β b ω γ 4(a + b ) R + 3a β b ω γ (a + b ) R, Θ = δ + H τ ɛ T + 3a β b ω γ (a + b ) R + 3a β b ω γ 4(a + b ) R, (6) 354-4

5 where a i =ω i A c τ c sin(ω i τ c ), b i = A c τ c cos(ω i τ c ) b, Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Q i = ω i(a i A c b iτ c +a i ω i b i + A c b3 i τ c + ω i b4 i + a4 i ω i a ia c ω ib i τ c a i A c τ c ω ib i a3 i A c τ c ω i )sin(ω iτ c ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) + A cτ c (a i b i ω ib i a i ω i + a 3 i + τ cω i b 3 i + a i τ cω i b i )cos(ω i τ c ) A c ω i τ c (a i + b i ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) A cω i (a i ω i a 3 i ω ib i a ib i a i τ cω i b i τ c ω i b 3 i )sin(ω iτ c ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) + ω i A c(b 3 i a ib i ω i + a i b i a 3 i τ cω i a i b i ω iτ c )cos(ω i τ c ) τ c (a i +, b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. δ i = [a i(ω i + ω + A c A c cos(ω i τ c )) b i A c sin(ω i τ c )]τ ε (a i + b i )τ c [a i(cos(ω i τ c ) ) + b i sin(ω i τ c )]A ε a i +, b i H i = [ω3 i (a i + b i ) A cτ c ω i (a 3 i τ cω i b 3 i a i b i +a i ω i b i + a i τ c ω i b i )] cos(ω iτ c ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) + A cτ c ω i (a i + b i ) ω i A cτ c (ω i b i τ cω i b 3 i a i ω i + a 3 i + a ib i + a i τ cω i b i )sin(ω i τ c ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) A cω i (a 3 i ω i a i ω i + ω i b i + τ cω i b3 i + a iω i b i + a i τ cω i b i)cos(ω i τ c ) τ c (a i + b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )) A cω i (ω i b 3 i a iω i b i + a i ω ib i a 3 i τ cω i a iτ c ω i b i )sin(ω iτ c ) τ c (a i +, i =,. b i ) (A c τ c cos(ω i τ c )b ω i sin(ω i τ c )).3. Center manifold reduction In this section, we compute the normal form near the double Hopf bifurcation critical point (A c,τ c ) using the CMR method. First, let ẋ = y. Then, system () can be rewritten as dx dt = y(t), dy dt = ω x(t)+[b γx (t)]y(t) βx 3 (7) (t) + A[x(t τ) x(t)]. Rescale the time by t (t/τ) to normalize the delay so that system (7) becomes dx d t = τy( t), dy d t = bτy( t) ω τx( t)+aτ[x( t ) x( t)] γτx ( t)y( t) βτx 3 ( t). (8) The trivial equilibrium of (8) isx = y =.Atthe critical point (A, τ) =(A c,τ c ), we choose τ c N, θ =, η(θ) =, θ (, ), τ c N, θ =, with ( ) N = ω A, c b ( ) N =. A c Then, the linearized equation of (8) at the trivial equilibrium is dx( t) = L X t d t, where L φ = dη(θ)φ(θ),φ C=C([, ], R ), and the bilinear form on C C ( stands for adjoint) is 354-5

6 Y. Ding et al. Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. θ ψ(s),φ(θ) = ψ()φ() ψ(ξ θ)dη(θ)φ(ξ)dξ, ξ= in which φ C, ψ C. Then, the phase space C is decomposed by Λ = {±iω, ±iω } as C = P Q, where Q = {ϕ C:(ψ, ϕ) =, for all ψ P },andthebasesforp and its adjoint P are given respectively by ( ) e iω τ cθ e iω τ cθ e iω τ cθ e iω τ cθ Φ(θ) =, iω e iω τ cθ iω e iω τ cθ iω e iω τ cθ iω e iω τ cθ and h (b iω )e iω τ cs h e iω τ cs h Ψ(s) = (b +iω )e iω τ cs h e iω τ cs h (b iω )e iω τ cs h e iω τ cs, h (b +iω )e iω τ cs h e iω τ cs where h j =(b iω j A c τ c e iω jτ c ),j =,. We also use the same bifurcation parameters given by A = A c + A ε and τ = τ c + τ ε in (8), where A ε and τ ε are perturbation parameters, and denote ε =(A ε,τ ε ). Then (8) can be written as dx( t) = L(ε)X t + F (X t,ε), (9) d t where ( L(ε)X t = (τ c + τ ε )y t () ) b(τ c + τ ε )y t () ω (τ c + τ ε )x t () + (A, c + A ε )(τ c + τ ε )[x t ( ) x t ()] and ( ) F (X t,ε)= γ(τ c + τ ε )x t ()y t () β(τ. c + τ ε )x () 3 t We now consider the enlarged phase space BC of functions from [, ] to R, which are continuous on [, ) with a possible jump discontinuity at zero. This space can be identified as C R. Thus, its elements can be written in the form ϕ = ϕ + X c,whereϕ C, c R and X is a matrix-valued function, defined by X (θ) =for θ [, ) and X () = I. In the BC,(9) becomes an abstract ODE, du d t = Au + X F (u, ε), () where u C, and A is defined by A : C BC, and Au = du d t + X [ L u du() d t F (u, ε) =[L(ε) L ]u + F (u, ε). By the continuous projection π : BC P, π(φ + X c)=φ[(ψ,φ)+ψ()c], we can decompose the enlarged phase space by Λ = {±iω τ c, ±iω τ c } ], as BC = P Ker π, where Ker π = {φ + X c : π(φ + X c)=}, denoting the Kernel under the projection π. Letη =(η, η,η, η ) T, v t Q := Q C Ker π,anda Q the restriction of A as an operator from Q to the Banach space Ker π. Further, denote u t = Φη + v t. Then, Eq. () is decomposed as dη d t = Bη +Ψ() F (Φη + v t,ε), dv t d t = A Q v t +(I π)x F (Φη + v t,ε), () where B =diag{iω, iω, iω, iω }. Next, let M denote the operator defined in V 6 (C4 Ker π ), with M : V 6 (C4 ) V 6 (C4 ), (M p)(η, ε) =D ηp(η, ε)bη Bp(η, ε), where V 6 (C4 ) represents the linear space of the second-order homogeneous polynomials in six 354-6

7 Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. variables (η, η,η, η,ε) with coefficients in C 4. Then, it is easy to verify that one may choose the decomposition V 6 (C4 )=Im(M ) Im(M )c with complementary space Im(M )c spanned by the elements A ε η e, τ ε η e, A ε η e, τ ε η e, A ε η e 3, τ ε η e 3, A ε η e 4, τ ε η e 4,wheree i (i =,, 3, 4) are unit vectors. Consequently, the normal form of (9) onthe center manifold associated with the equilibrium (, ) near A ε =,τ ε = has the form dη d t = Bη + g (η,,ε)+h.o.t., where g is the function giving the quadratic terms in (η, ε) for v t =, and is determined by g (η,,ε) = Proj (Im(M f ))c (η,,ε), where f (η,,ε) is the function giving the quadratic terms in (η, ε) for v t = defined by the first equation of (). Thus, the normal form, truncated at the quadratic order terms, is given by dη d t =iω τ c η h A ε τ c (e iω τ c )η + h [ω + ω A c(e iω τ c )]τ ε η, dη d t =iω τ c η h A ε τ c (e iω τ c )η + h [ω + ω A c(e iω τ c )]τ ε η, () where h j =(b iω j A c τ c e iω jτ c ) (j =, ). To find the normal form up to third order, similarly, let M 3 denote the operator defined in V 4 3 (C4 Ker π ), with M 3 : V 4 3(C 4 ) V 4 3(C 4 ), (M 3 p)(η, ε) =D ηp(η, ε)bη Bp(η, ε), Therefore, the normal form up to third-order terms is given by dη d t = Bη +! g (η,,ε)+ 3! g 3 (η,,ε)+h.o.t., (3) where 3! g 3(η,, ) = 3! (I P I,3)f 3(η,, ), and f 3 (η,, ) is the function giving the cubic terms in (η, ε, v t )forε =,andv t = is defined by the first equation of (). Finally, the normal form on the center manifold arising from () becomes dη d t =iω τ c η h A ε τ c (e iω τ c )η + h [ω + ω A c(e iω τ c )]τ ε η + h τ c (iω γ +3β)η η +h τ c (iω γ +3β)η η η, dη d t =iω τ c η h A ε τ c (e iω τ c )η + h [ω + ω A c(e iω τ c )]τ ε η + h τ c (iω γ +3β)η η +h τ c (iω γ +3β)η η η, (4) where h j =(b iω j A c τ c e iω jτ c ), j =,. With the polar coordinates: η = R e iθ, η = R e iθ, combining with (3), we finally obtain the amplitude and phase equations of (4) on the center manifold as dr = µ R + P R 3 dt + P R R, where V 4 3 (C4 ) denotes the linear space of the thirdorder homogeneous polynomials in four variables (η, η,η, η ) with coefficients in C 4. Then, one may choose the decomposition V 4 3 (C4 )=Im(M 3 ) Im(M 3 )c with complementary space Im(M 3 )c spanned by the elements η η e, η η η e, η η e, η η η e, η η e 3, η η η e 3, η η e 4, η η η e 4,where e i (i =,, 3, 4) are unit vectors. where dr dt dθ dt dθ dt = µ R + P R R + P R 3, = ω + δ + Q R +Q R, = ω + δ +Q R + Q R, (5) µ i = [a i cos(ω i τ c )+b i sin(ω i τ c )]A c ω i τ ε (a i + b i )τ c + [b i cos(ω i τ c ) b i a i sin(ω i τ c )]A ε a i +, b i 354-7

8 Y. Ding et al. δ i = [a i(ω i + ω + A c A c cos(ω i τ c )) b i A c sin(ω i τ c )]τ ε (a i + b i )τ c P ii = 3βb i + ω i γa i 4(a i + b i ), P =P, P =P, Q i = 3a iβ b i ω i γ 4(a i + b i ), [a i(cos(ω i τ c ) ) + b i sin(ω i τ c )]A ε a i +, b i Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. with a i =ω i A c τ c sin(ω i τ c ), b i = A c τ c cos(ω i τ c ) b, i =,. Remark. In the MTS method, we assume that the solution is in a real form, while in the CMR method, we consider the solution to be in a complex form. In fact, the solution to Eq. (7) can be expressed in either a complex form or a real form, that is, x = K e iω T + K e iω T + K e iω T + K e iω T = G (T,T )sin(ω T ) + G (T,T )cos(ω T ) + G 3 (T,T )sin(ω T ) + G 4 (T,T )cos(ω T ). These two different form solutions become identical under the polar coordinate transformations: K = R e iθ, K = R e iθ and G = R sin(θ ), G = R cos(θ ), G 3 = R sin(θ ), G 4 = R cos(θ ). Therefore, in the CMR method, we use the polar coordinates: η = R e iθ, η = R e iθ in order to keep consistence..4. Comparison of the MTS and CMR methods Equation (6) is the normal form derived by using the MTS method, and Eq. (5) is the normal form derived by using the CMR method. Comparing the two normal forms, we have found that Eq. (6) has four more terms Q i τ ɛr i, H i τ ɛ (i =, ) than Eq. (5) does, which come from the Taylor expansions of x j (T τ c ɛτ ɛ,t ɛ(τ c + ɛτ ɛ ),T ɛ (τ c + ɛτ ɛ )) at x j (T τ c,t,t )forj =,, 3,...In fact, in the CMR method, we treat x(t τ) as a linear term of t, and thus no higher-order terms appear, while in the MTS method, we expand x(t τ) in T,T,... as x j (T τ c ɛτ ɛ,t ɛ(τ c + ɛτ ɛ ),T ɛ (τ c + ɛτ ɛ )) at x j (T τ c,t,t ) up to infinite orders, which results in (ɛτ ɛ ) (ɛτ ɛ) 4 + (ɛτ ɛ) 6 (ɛτ ɛ)! 4! 6! 8! =cos(ɛτ ɛ ), + ɛτ ɛ + (ɛτ ɛ) 3 (ɛτ ɛ) 5 3! 5! = sin(ɛτ ɛ ), + (ɛτ ɛ) 7 7! + but we only take terms up to ɛ -order in the MTS method. That is why Eq. (6) has additional Q i τ ɛ R i and H i τ ɛ (i =, ) terms. Or simply, in the MTS method, if we do not expand x(t τ) int,t,... (which would be the same as the CMR method), then these two normal forms are identical. Nevertheless, the small difference has very little effect on the dynamical analysis. 3. Bifurcation Analysis and Numerical Simulation In this section, we first give a bifurcation analysis based on the normal form (5),andthenpresent some numerical simulation results. 3.. Bifurcation analysis For the normal form (5),accordingtothesigns of P P, there exist two different cases, i.e. simple case (with no periodic solutions) and difficult case (with periodic solutions) [Guckenheimer & Holmes, 99]. Here, we are interested in the difficult case, i.e. when P P <. Without loss of generality, we assume P > andp <. Let r j = P jj R j, (j =, ). Then, we have the following planar system in terms of r and r : {ṙ = r (µ + r κr ), (6) ṙ = r (µ + χr r ), where κ = P P, χ = P Note that M =(, ) is always an equilibrium of (6). The two semi-trivial equilibria are given in terms of perturbation parameters as M =( µ, ) and M =(,µ ), which bifurcate from the origin on the bifurcation lines L : µ =andl : µ =, respectively. There may also exist a nontrivial equilibrium M 3 = ( µ κµ κχ, µ χµ κχ ). For this equilibrium to exist, it needs κχ. The nontrivial equilibrium M 3 collides with a semi-trivial one on P

9 Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. the bifurcation line T : µ κµ = or T : µ χµ =.If( χ)µ +( κ)µ <, the fixed point M 3 is a sink, otherwise M 3 is a source. Therefore, we further need to consider the bifurcation line T 3 :( χ)µ +( κ)µ =. In order to give a more clear bifurcation picture, choose A c =.359, τ c =6.7583, b =.774, β =.3, γ =and ω ω = 5. Then, the characteristic equation () has two pairs of purely imaginary eigenvalues Λ = {±iω, ±iω }, and all the other eigenvalues have negative real part. Assume that system () undergoes a double Hopf bifurcation from the equilibrium (, ). By a simple calculation, we obtain ω =.8, τ () =.557, τ () =6.7583, Re (λ ) >, (7) ω =.85, τ () =.9433, τ () =6.7583, Re (λ ) <, and µ =.437τ ε +.34A ε, µ =.356τ ε.63a ε, P =.54, P =.449, P =.4338, P =.69, κ =.8587 and χ = Therefore, the critical bifurcation lines become: L : A ε =.3653τ ε, L : A ε =.59τ ε, T : A ε =.756τ ε, T : A ε =.336τ ε, T 3 : A ε =.376τ ε, as shown in the bifurcation diagram (see Fig. ). Since there does not exist unstable manifold containing the equilibrium, according to the center manifold theory, the solutions on the center manifold determine the asymptotic behavior of the solutions of the original system (). Therefore, if Eq. (6) has one or two asymptotically stable (unstable) semi-trivial equilibria M and M,then() has one or two asymptotically stable (unstable) periodic solutions in the neighborhood of the trivial equilibrium. If Eq. (6) has an asymptotically stable (unstable) equilibrium M 3,then()has an asymptotically stable (unstable) quasi-periodic solution in the neighborhood of (, ). So, we shall call the periodic solution the source (respectively, saddle, sink) periodic solution of () ifthesemitrivial equilibrium of (6) is a source (respectively, saddle, sink), and call the quasi-periodic solution the source (respectively, saddle, sink) quasi-periodic solution of () when the nontrivial equilibrium of (6) is a source (respectively, saddle, sink). For the bifurcation behaviors of the original system () in the neighborhood of the trivial equilibrium, the above critical bifurcation boundaries divide the (A ε,τ ε ) parameter plane into seven regions (see Fig. ). We explain the bifurcations in the clockwise direction, starting from B and ending at B.First,inregionB, there is only one trivial equilibrium which is a saddle. When the parameters are varied across the line L from region B (a) (b) Fig.. (a) Critical bifurcation lines in the (τ ε,a ε) parameter space near (τ c,a c); (b) the corresponding phase portraits in the (r,r )plane

10 Y. Ding et al Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. Fig Simulated solution of system () for(a ε,τ ε)=(.5,.5), showing a stable fixed point. to B, the trivial equilibrium becomes a sink, and an unstable periodic solution O (saddle) appears from the trivial solution due to a Hopf bifurcation. Similarly, when the parameters are changed from region B to B 3, another periodic solution O (sink) occurs from the trivial solution due to a Hopf bifurcation while the trivial equilibrium becomes a saddle. In region B 4, a quasi-periodic solution (sink) occurs from the periodic solution O due to a Neimark Sacker bifurcation, and O is changed to a saddle from a sink. Further, the quasi-periodic solution (sink) becomes a source when the parameters are varied across line T 3 from region B 4 to B 5, and when the parameters are further changed from region B 5 to B 6 crossing the line T,thequasiperiodic solution collides with the periodic solution O and then disappears, and O becomes a source. When the parameters are further varied across line L from region B 6 to B 7, the periodic solution O collides with the trivial solution and then disappears, and the trivial solution becomes a source from a saddle. Finally, when the parameters are varied across line L from region B 7 to B,theperiodic solution O (saddle) collides with the trivial solution and then disappears, and the trivial solution becomes a saddle from a source. 3.. Numerical simulation To demonstrate the analytic results obtained in Sec. 3., here we present some numerical simulation results. We choose three groups of perturbation parameter values: (A ε,τ ε ) = (.5,.5), (.,.) and (.5,.5), belonging to the regions B, B 3 and B 4, corresponding to a stablefixedpointshowninfig., a stable periodic solution as depicted in Fig. 3, and a stable quasiperiodic solution, see Fig. 4 respectively. It is clear that the numerical simulations agree with the analytical predictions. 4. Conclusion and Discussion In this paper, we have discussed double Hopf bifurcation in delayed van der Pol Duffing equation. We derived the normal form of double Hopf bifurcation by using multiple time scales and center manifold reduction methods. A comparison between the two methods shows that the two normal forms are identical if the higher-order terms obtained in the MTS method, due to the expansion in the delayed variable with respect to time scales, are ignored. Moreover, bifurcation analysis near the double Hopf critical point is given, showing that the system may exhibit a stable fixed point, periodic solutions, and quasi-periodic solutions in the neighborhood of the critical point. Numerical simulations are given to verify the analytical predictions. The normal form method and bifurcation analysis presented in this paper are for local dynamical behaviors. But, surprisingly, we have found that even for parameter values not chosen in the 354-

11 Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only (a) Fig. 3. Simulated solution of system () for(a ε,τ ε)=(.,.): (a) the time history; (b) the phase portrait, showing a stable periodic solution (b) (a) (b) Fig. 4. Simulated solution of system () for(a ε,τ ε)=(.5,.5): (a) the time history; (b) the phase portrait, showing a stable quasi-periodic solution. 354-

12 Y. Ding et al Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only (a) Fig. 5. Simulated solution of system () for(a ε,τ ε)=( 3, 6): (a) the time history; (b) the phase portrait, showing a stable periodic solution (b) (a) 3 4 (b) Fig. 6. Simulated solution of system () for(a ε,τ ε)=( 3, 6): (a) the time history; (b) the phase portrait, showing a stable periodic solution. 354-

13 Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only (a) 4 6 Fig. 7. Simulated solution of system () for(a ε,τ ε)=(3, 6): (a) the time history; (b) the phase portrait, showing a stable quasi-periodic solution (b) (a) 4 (b) Fig. 8. Simulated solution of system () with a nonlinear delayed feedback: (a) the time history; (b) the phase portrait, showing a chaotic attractor

14 Y. Ding et al. Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only Fig. 9. The simulated Poincaré map of the chaotic attractor shown in Fig. 8(b) with Poincaré sectionx(t τ) =. neighborhood of the critical point, most simulation results are periodic or quasi-periodic solutions (see Figs. 5 7). This seemingly noncomplex dynamics is perhaps due to system () only containing time delay in a (feedback) linear term. If we introduce time delay in nonlinear terms, we may obtain more complicated dynamical behavior, such as chaos. Indeed, for example, when we change the linear delay feedback A[x(t τ) x(t)] to a nonlinear delay feedback Ax(t τ)[x(t τ) x(t)], and select A =5,τ =4,ω =,b =.774,β =.3 and γ =, we obtain chaotic motion, as shown in Figs. 8 and 9. Acknowledgments The work was supported by the National Natural Science Foundation of China (NSFC) and the National Science and Engineering Research Council of Canada (NSERC), the Heilongjiang Provincial Natural Science Foundation (No. A86), and the Program of Excellent Team in HIT. The first author also acknowledges the financial support received from the China Scholarship Council for visiting Western University. References Choi, Y. & LeBlanc, V. G. [6] Toroidal normal forms for bifurcations in retarded functional differential equations I: Multiple Hopf and transcritical/multiple Hopf interaction, J. Diff. Eqns. 7, Das, S. L. & Chatterjee, A. [] Multiple scales without center manifold reductions for delay differential equations near Hopf bifurcations, Nonlin. Dyn. 3, Faria, T. & Magalhães, L. T. [995] Normal form for retarded functional differential equations and applications to Bogdanov Takens singularity, J. Diff. Eqns., 4. Guckenheimer, J. & Holmes, P. [99] Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, 3rd edition (Springer, NY). Hale, J. [997] Theory of Functional Differential Equations (Springer, NY). Hassard, B. D., Kazarinoff, N. D. & Wan, Y. H. [98] Theory and Applications of Hopf Bifurcation (Cambridge University Press, Cambridge). Hu, H. & Wang, Z. [9] Singular perturbation methods for nonlinear dynamic systems with time delays, Chaos Solit. Fract. 4, 3 7. Jiang, W. & Yuan, Y. [7] Bogdanov Takens singularity in van der Pol s oscillator with delayed feedback, Physica D 7, Jiang, W. & Wei, J. [8] Bifurcation analysis in van der Pol s oscillator with delayed feedback, J. Comput. Appl. Math. 3, Jiang, W. & Wang, H. [] Hopf-transcritical bifurcation in retarded functional differential equations, Nonlin. Anal. 73, Kuznetsov, Y. A. [4] Elements of Applied Bifurcation Theory, 3rd edition (Springer, NY). Ma, S., Lu, Q. & Feng, Z. [8] Double Hopf bifurcation for van der Pol Duffing oscillator with parametric delay feedback control, J. Math. Anal. Appl. 338, Nayfeh, A. H. [973] Perturbation Methods (Wiley- Interscience, NY). Nayfeh, A. H. [98] Introduction to Perturbation Techniques (Wiley-Interscience, NY). Nayfeh, A. H. [8] Order reduction of retarded nonlinear systems The method of multiple scales versus center-manifold reduction, Nonlin. Dyn. 5, Revel, G., Alonso, D. M. & Moiola, J. L. [] Interactions between oscillatory modes near :3 resonant Hopf Hopf bifurcation, Chaos, 436. Wang, H. & Jiang, W. [] Hopf-pitchfork bifurcation in van der Pol s oscillator with nonlinear delayed feedback, J. Math. Anal. Appl. 368, 9 8. Wei, J. [7] Bifurcation analysis in a scalar delay differential equation, Nonlinearity, Wiggins, S. [99] Introduction to Applied Nonlinear Dynamical Systems and Chaos (Springer, NY). Xu, J. & Chung, K. W. [3] Effects of time delayed position feedback on a van der Pol Duffing oscillator, Physica D 8,

15 Int. J. Bifurcation Chaos 3.3. Downloaded from by CITY UNIVERSITY OF HONG KONG on 3/5/3. For personal use only. Double Hopf Bifurcation in Delayed van der Pol Duffing Equation Xu, J. & Yu, P. [4] Delay-induced bifurcations in nonautonomous system with delayed velocity feedbacks, Int. J. Bifurcation and Chaos 4, Xu, J., Chung, K. W. & Chan, C. L. [7] An efficient method for studying weak resonant double Hopf bifurcation in nonlinear systems with delayed feedback, SIAM J. Appl. Dyn. Syst. 6, 9 6. Yu, P. [998] Computation of normal forms via a perturbation technique, J. Sound Vibr., Yu, P. [] Symbolic computation of normal forms for resonant double Hopf bifurcations using a perturbation technique, J. Sound Vibr. 47, Yu, P. [] Analysis on double Hopf bifurcation using computer algebra with the aid of multiple scales, Nonlin. Dyn. 7, Yu, P., Yuan, Y. & Xu, J. [] Study of double Hopf bifurcation and chaos for an oscillator with time delayed feedback, Commun. Nonlin. Sci. Numer. Simul. 7, Yuan, Y. & Wei, J. [6] Multiple bifurcation analysis in a neural network model with delays, Int. J. Bifurcation and Chaos 6, Zheng, Y. & Wang, Z. [] Stability and Hopf bifurcation of a class of TCP/AQM networks, Nonlin. Anal. Real World Appl.,

An Efficient Method for Studying Weak Resonant Double Hopf Bifurcation in Nonlinear Systems with Delayed Feedbacks

An Efficient Method for Studying Weak Resonant Double Hopf Bifurcation in Nonlinear Systems with Delayed Feedbacks An Efficient Method for Studying Weak Resonant Double Hopf Bifurcation in Nonlinear Systems with Delayed Feedbacks J. Xu 1,, K. W. Chung 2 C. L. Chan 2 1 School of Aerospace Engineering and Applied Mechanics,

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

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

ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM. China

ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM. China Mathematical and Computational Applications, Vol. 9, No., pp. 84-9, 4 ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM Ping Cai,, Jia-Shi Tang, Zhen-Bo Li College of

More information

Hopf-pitchfork bifurcation of coupled van der Pol oscillator with delay

Hopf-pitchfork bifurcation of coupled van der Pol oscillator with delay ISSN 1392-5113 Nonlinear Analysis: Modelling and Control, 2017, Vol. 22, No. 5, 598 613 https://doi.org/10.15388/na.2017.5.2 Hopf-pitchfork bifurcation of coupled van der Pol oscillator with delay Yuting

More information

HOPF-ZERO BIFURCATION OF A DELAYED PREDATOR-PREY MODEL WITH DORMANCY OF PREDATORS

HOPF-ZERO BIFURCATION OF A DELAYED PREDATOR-PREY MODEL WITH DORMANCY OF PREDATORS Journal of Applied Analysis and Computation Volume 7 Number August 7 5 69 Website:http://jaac-online.com/ DOI:.948/766 HOPF-ZERO BIFURCATION OF A DELAYED PREDATOR-PREY MODEL WITH DORMANCY OF PREDATORS

More information

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

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

More information

SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM

SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM International Journal of Bifurcation and Chaos, Vol. 23, No. 11 (2013) 1350188 (7 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127413501885 SIMPLE CHAOTIC FLOWS WITH ONE STABLE EQUILIBRIUM

More information

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)

Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x

More information

Hopf Bifurcation Analysis of a Dynamical Heart Model with Time Delay

Hopf Bifurcation Analysis of a Dynamical Heart Model with Time Delay Applied Mathematical Sciences, Vol 11, 2017, no 22, 1089-1095 HIKARI Ltd, wwwm-hikaricom https://doiorg/1012988/ams20177271 Hopf Bifurcation Analysis of a Dynamical Heart Model with Time Delay Luca Guerrini

More information

Analytical estimations of limit cycle amplitude for delay-differential equations

Analytical estimations of limit cycle amplitude for delay-differential equations Electronic Journal of Qualitative Theory of Differential Equations 2016, No. 77, 1 10; doi: 10.14232/ejqtde.2016.1.77 http://www.math.u-szeged.hu/ejqtde/ Analytical estimations of limit cycle amplitude

More information

Existence of permanent oscillations for a ring of coupled van der Pol oscillators with time delays

Existence of permanent oscillations for a ring of coupled van der Pol oscillators with time delays Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 14, Number 1 (2018), pp. 139 152 Research India Publications http://www.ripublication.com/gjpam.htm Existence of permanent oscillations

More information

Chapter 14 Three Ways of Treating a Linear Delay Differential Equation

Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Chapter 14 Three Ways of Treating a Linear Delay Differential Equation Si Mohamed Sah and Richard H. Rand Abstract This work concerns the occurrence of Hopf bifurcations in delay differential equations

More information

Solutions of Nonlinear Oscillators by Iteration Perturbation Method

Solutions of Nonlinear Oscillators by Iteration Perturbation Method Inf. Sci. Lett. 3, No. 3, 91-95 2014 91 Information Sciences Letters An International Journal http://dx.doi.org/10.12785/isl/030301 Solutions of Nonlinear Oscillators by Iteration Perturbation Method A.

More information

2:2:1 Resonance in the Quasiperiodic Mathieu Equation

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,

More information

Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited

Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited Dynamics of Two Coupled van der Pol Oscillators with Delay Coupling Revisited arxiv:1705.03100v1 [math.ds] 8 May 017 Mark Gluzman Center for Applied Mathematics Cornell University and Richard Rand Dept.

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

(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

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

Time-delay feedback control in a delayed dynamical chaos system and its applications

Time-delay feedback control in a delayed dynamical chaos system and its applications Time-delay feedback control in a delayed dynamical chaos system and its applications Ye Zhi-Yong( ), Yang Guang( ), and Deng Cun-Bing( ) School of Mathematics and Physics, Chongqing University of Technology,

More information

DELAY-INDUCED BIFURCATIONS IN A NONAUTONOMOUS SYSTEM WITH DELAYED VELOCITY FEEDBACKS

DELAY-INDUCED BIFURCATIONS IN A NONAUTONOMOUS SYSTEM WITH DELAYED VELOCITY FEEDBACKS International Journal of Bifurcation and Chaos, Vol. 4, No. 8 (2004) 2777 2798 c World Scientific Publishing Company DELAY-INDUCED BIFURCATIONS IN A NONAUTONOMOUS SYSTEM WITH DELAYED VELOCITY FEEDBACKS

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

Hopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators

Hopf Bifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators The Open Acoustics Journal 8 9-3 9 Open Access Hopf ifurcation Analysis and Approximation of Limit Cycle in Coupled Van Der Pol and Duffing Oscillators Jianping Cai *a and Jianhe Shen b a Department of

More information

HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY

HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY Vol. 37 17 No. J. of Math. PRC HOPF BIFURCATION ANALYSIS OF A PREDATOR-PREY SYSTEM WITH NON-SELECTIVE HARVESTING AND TIME DELAY LI Zhen-wei, LI Bi-wen, LIU Wei, WANG Gan School of Mathematics and Statistics,

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

Direction and Stability of Hopf Bifurcation in a Delayed Model with Heterogeneous Fundamentalists

Direction and Stability of Hopf Bifurcation in a Delayed Model with Heterogeneous Fundamentalists International Journal of Mathematical Analysis Vol 9, 2015, no 38, 1869-1875 HIKARI Ltd, wwwm-hikaricom http://dxdoiorg/1012988/ijma201554135 Direction and Stability of Hopf Bifurcation in a Delayed Model

More information

Three ways of treating a linear delay differential equation

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 *,

More information

In-Plane and Out-of-Plane Dynamic Responses of Elastic Cables under External and Parametric Excitations

In-Plane and Out-of-Plane Dynamic Responses of Elastic Cables under External and Parametric Excitations Applied Mathematics 5, 5(6): -4 DOI:.59/j.am.556. In-Plane and Out-of-Plane Dynamic Responses of Elastic Cables under External and Parametric Excitations Usama H. Hegazy Department of Mathematics, Faculty

More information

Citation Acta Mechanica Sinica/Lixue Xuebao, 2009, v. 25 n. 5, p The original publication is available at

Citation Acta Mechanica Sinica/Lixue Xuebao, 2009, v. 25 n. 5, p The original publication is available at Title A hyperbolic Lindstedt-poincaré method for homoclinic motion of a kind of strongly nonlinear autonomous oscillators Author(s) Chen, YY; Chen, SH; Sze, KY Citation Acta Mechanica Sinica/Lixue Xuebao,

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

The Invariant Curve in a Planar System of Difference Equations

The Invariant Curve in a Planar System of Difference Equations Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 13, Number 1, pp. 59 71 2018 http://campus.mst.edu/adsa The Invariant Curve in a Planar System of Difference Equations Senada Kalabušić,

More information

CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION

CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION WARREN WECKESSER Department of Mathematics Colgate University Hamilton, NY 3346 E-mail: wweckesser@mail.colgate.edu Cartwright and Littlewood discovered

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

BIFURCATIONS OF PERIODIC ORBITS IN THREE-WELL DUFFING SYSTEM WITH A PHASE SHIFT

BIFURCATIONS OF PERIODIC ORBITS IN THREE-WELL DUFFING SYSTEM WITH A PHASE SHIFT J Syst Sci Complex (11 4: 519 531 BIFURCATIONS OF PERIODIC ORBITS IN THREE-WELL DUFFING SYSTEM WITH A PHASE SHIFT Jicai HUANG Han ZHANG DOI: 1.17/s1144-1-89-3 Received: 9 May 8 / Revised: 5 December 9

More information

CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS

CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems

More information

2:1:1 Resonance in the Quasi-Periodic Mathieu Equation

2:1:1 Resonance in the Quasi-Periodic Mathieu Equation Nonlinear Dynamics 005) 40: 95 03 c pringer 005 :: Resonance in the Quasi-Periodic Mathieu Equation RICHARD RAND and TINA MORRION Department of Theoretical & Applied Mechanics, Cornell University, Ithaca,

More information

Frequency locking in a forced Mathieu van der Pol Duffing system

Frequency locking in a forced Mathieu van der Pol Duffing system Nonlinear Dyn (008) 54:3 1 DOI 10.1007/s11071-007-938-x ORIGINAL ARTICLE requency locking in a forced Mathieu van der Pol Duffing system Manoj Pandey Richard H. Rand Alan T. Zehnder Received: 4 August

More information

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS

ZERO-HOPF BIFURCATION FOR A CLASS OF LORENZ-TYPE SYSTEMS This is a preprint of: Zero-Hopf bifurcation for a class of Lorenz-type systems, Jaume Llibre, Ernesto Pérez-Chavela, Discrete Contin. Dyn. Syst. Ser. B, vol. 19(6), 1731 1736, 214. DOI: [doi:1.3934/dcdsb.214.19.1731]

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

MIT Weakly Nonlinear Things: Oscillators.

MIT Weakly Nonlinear Things: Oscillators. 18.385 MIT Weakly Nonlinear Things: Oscillators. Department of Mathematics Massachusetts Institute of Technology Cambridge, Massachusetts MA 02139 Abstract When nonlinearities are small there are various

More information

Journal of Applied Nonlinear Dynamics

Journal of Applied Nonlinear Dynamics Journal of Applied Nonlinear Dynamics 1(2) (2012) 159 171 Journal of Applied Nonlinear Dynamics https://lhscientificpublishing.com/journals/jand-default.aspx Asymptotic Analysis of the Hopf-Hopf Bifurcation

More information

Problem Set Number 5, j/2.036j MIT (Fall 2014)

Problem Set Number 5, j/2.036j MIT (Fall 2014) Problem Set Number 5, 18.385j/2.036j MIT (Fall 2014) Rodolfo R. Rosales (MIT, Math. Dept.,Cambridge, MA 02139) Due Fri., October 24, 2014. October 17, 2014 1 Large µ limit for Liénard system #03 Statement:

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

NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011

NBA Lecture 1. Simplest bifurcations in n-dimensional ODEs. Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 NBA Lecture 1 Simplest bifurcations in n-dimensional ODEs Yu.A. Kuznetsov (Utrecht University, NL) March 14, 2011 Contents 1. Solutions and orbits: equilibria cycles connecting orbits other invariant sets

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

α Cubic nonlinearity coefficient. ISSN: x DOI: : /JOEMS

α Cubic nonlinearity coefficient. ISSN: x DOI: : /JOEMS Journal of the Egyptian Mathematical Society Volume (6) - Issue (1) - 018 ISSN: 1110-65x DOI: : 10.1608/JOEMS.018.9468 ENHANCING PD-CONTROLLER EFFICIENCY VIA TIME- DELAYS TO SUPPRESS NONLINEAR SYSTEM OSCILLATIONS

More information

Restrictions and unfolding of double Hopf bifurcation in functional differential equations

Restrictions and unfolding of double Hopf bifurcation in functional differential equations Restrictions and unfolding of double Hopf bifurcation in functional differential equations Pietro-Luciano Buono Jacques Bélair CRM-837 March Centre de recherches mathématiques, Université de Montréal,

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

WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY

WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY WIDELY SEPARATED FREQUENCIES IN COUPLED OSCILLATORS WITH ENERGY-PRESERVING QUADRATIC NONLINEARITY J.M. TUWANKOTTA Abstract. In this paper we present an analysis of a system of coupled oscillators suggested

More information

2:1 Resonance in the delayed nonlinear Mathieu equation

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

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

The diffusive Lotka-Volterra predator-prey system with delay

The diffusive Lotka-Volterra predator-prey system with delay University of Wollongong Research Online Faculty of Engineering and Information Sciences - Papers: Part A Faculty of Engineering and Information Sciences 5 The diffusive Lotka-Volterra predator-prey system

More information

Invariant manifolds of the Bonhoeffer-van der Pol oscillator

Invariant manifolds of the Bonhoeffer-van der Pol oscillator Invariant manifolds of the Bonhoeffer-van der Pol oscillator R. Benítez 1, V. J. Bolós 2 1 Dpto. Matemáticas, Centro Universitario de Plasencia, Universidad de Extremadura. Avda. Virgen del Puerto 2. 10600,

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

HOPF BIFURCATION CONTROL WITH PD CONTROLLER

HOPF BIFURCATION CONTROL WITH PD CONTROLLER HOPF BIFURCATION CONTROL WITH PD CONTROLLER M. DARVISHI AND H.M. MOHAMMADINEJAD DEPARTMENT OF MATHEMATICS, FACULTY OF MATHEMATICS AND STATISTICS, UNIVERSITY OF BIRJAND, BIRJAND, IRAN E-MAILS: M.DARVISHI@BIRJAND.AC.IR,

More information

HOPF BIFURCATION FOR TUMOR-IMMUNE COMPETITION SYSTEMS WITH DELAY

HOPF BIFURCATION FOR TUMOR-IMMUNE COMPETITION SYSTEMS WITH DELAY Electronic Journal of Differential Equations Vol. 4 (4) No. 7 pp. 3. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu HOPF BIFURCATION FOR TUMOR-IMMUNE

More information

2.034: Nonlinear Dynamics and Waves. Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen

2.034: Nonlinear Dynamics and Waves. Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen 2.034: Nonlinear Dynamics and Waves Term Project: Nonlinear dynamics of piece-wise linear oscillators Mostafa Momen May 2015 Massachusetts Institute of Technology 1 Nonlinear dynamics of piece-wise linear

More information

Hamad Talibi Alaoui and Radouane Yafia. 1. Generalities and the Reduced System

Hamad Talibi Alaoui and Radouane Yafia. 1. Generalities and the Reduced System FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. Vol. 22, No. 1 (27), pp. 21 32 SUPERCRITICAL HOPF BIFURCATION IN DELAY DIFFERENTIAL EQUATIONS AN ELEMENTARY PROOF OF EXCHANGE OF STABILITY Hamad Talibi Alaoui

More information

STABILITY AND HOPF BIFURCATION ON A TWO-NEURON SYSTEM WITH TIME DELAY IN THE FREQUENCY DOMAIN *

STABILITY AND HOPF BIFURCATION ON A TWO-NEURON SYSTEM WITH TIME DELAY IN THE FREQUENCY DOMAIN * International Journal of Bifurcation and Chaos, Vol. 7, No. 4 (27) 355 366 c World Scientific Publishing Company STABILITY AND HOPF BIFURCATION ON A TWO-NEURON SYSTEM WITH TIME DELAY IN THE FREQUENCY DOMAIN

More information

DynamicsofTwoCoupledVanderPolOscillatorswithDelayCouplingRevisited

DynamicsofTwoCoupledVanderPolOscillatorswithDelayCouplingRevisited Global Journal of Science Frontier Research: F Mathematics and Decision Sciences Volume 7 Issue 5 Version.0 Year 07 Type : Double Blind Peer Reviewed International Research Journal Publisher: Global Journals

More information

{ =y* (t)=&y(t)+ f(x(t&1), *)

{ =y* (t)=&y(t)+ f(x(t&1), *) Journal of Differential Equations 69, 6463 (00) doi:0.006jdeq.000.390, available online at http:www.idealibrary.com on The Asymptotic Shapes of Periodic Solutions of a Singular Delay Differential System

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

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

DYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS

DYNAMICS OF THREE COUPLED VAN DER POL OSCILLATORS WITH APPLICATION TO CIRCADIAN RHYTHMS Proceedings of IDETC/CIE 2005 ASME 2005 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference September 24-28, 2005, Long Beach, California USA DETC2005-84017

More information

Simplest Chaotic Flows with Involutional Symmetries

Simplest Chaotic Flows with Involutional Symmetries International Journal of Bifurcation and Chaos, Vol. 24, No. 1 (2014) 1450009 (9 pages) c World Scientific Publishing Company DOI: 10.1142/S0218127414500096 Simplest Chaotic Flows with Involutional Symmetries

More information

Models Involving Interactions between Predator and Prey Populations

Models Involving Interactions between Predator and Prey Populations Models Involving Interactions between Predator and Prey Populations Matthew Mitchell Georgia College and State University December 30, 2015 Abstract Predator-prey models are used to show the intricate

More information

On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction

On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction rspa.royalsocietypublishing.org Research Article submitted to journal On the analysis of the double Hopf bifurcation in machining processes via center manifold reduction T. G. Molnar, Z. Dombovari, T.

More information

Thermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space

Thermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space 1/29 Thermosolutal Convection at Infinite Prandtl Number with or without rotation: Bifurcation and Stability in Physical Space Jungho Park Department of Mathematics New York Institute of Technology SIAM

More information

Oscillation Control in Delayed Feedback Systems

Oscillation Control in Delayed Feedback Systems Oscillation Control in Delayed Feedback Systems Fatihcan M. Atay (atay@member.ams.org) Preprint. Final version in Lecture Notes in Control and Information Sciences Vol. 273, pp. 3 6 (22) Abstract. The

More information

DYNAMICS OF ASYMMETRIC NONLINEAR VIBRATION ABSORBER

DYNAMICS OF ASYMMETRIC NONLINEAR VIBRATION ABSORBER 8 Journal of Marine Science and Technology, Vol. 11, No. 1, pp. 8-19 (2003) DYNAMICS OF ASYMMETRIC NONLINEAR VIBRATION ABSORBER Chiou-Fong Chung and Chiang-Nan Chang* Key words: bifurcation, absorber,

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

Math Ordinary Differential Equations

Math Ordinary Differential Equations Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x

More information

DETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING

DETC DELAY DIFFERENTIAL EQUATIONS IN THE DYNAMICS OF GENE COPYING Proceedings of the ASME 2007 International Design Engineering Technical Conferences & Computers and Information in Engineering Conference IDETC/CIE 2007 September 4-7, 2007, Las Vegas, Nevada, USA DETC2007-3424

More information

BIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1

BIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1 KYBERNETIKA VOLUME 3 0 (1994), NUMBER 2, PAG ES 121-128 BIFURCATIONS AND CHAOS IN A PERIODICALLY FORCED PROTOTYPE ADAPTIVE CONTROL SYSTEM 1 YURI A. KUZNETSOV AND CARLO PICCARDI An adaptive control system

More information

A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS

A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS International Journal of Bifurcation and Chaos, Vol. 18, No. 5 (2008) 1567 1577 c World Scientific Publishing Company A MINIMAL 2-D QUADRATIC MAP WITH QUASI-PERIODIC ROUTE TO CHAOS ZERAOULIA ELHADJ Department

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

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

arxiv: v1 [math.ds] 20 Sep 2016

arxiv: v1 [math.ds] 20 Sep 2016 THE PITCHFORK BIFURCATION arxiv:1609.05996v1 [math.ds] 20 Sep 2016 Contents INDIKA RAJAPAKSE AND STEVE SMALE Abstract. We give development of a new theory of the Pitchfork bifurcation, which removes the

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

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

Math 266: Phase Plane Portrait

Math 266: Phase Plane Portrait Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions

More information

New approach to study the van der Pol equation for large damping

New approach to study the van der Pol equation for large damping Electronic Journal of Qualitative Theor of Differential Equations 2018, No. 8, 1 10; https://doi.org/10.1422/ejqtde.2018.1.8 www.math.u-szeged.hu/ejqtde/ New approach to stud the van der Pol equation for

More information

Stability and Hopf-bifurcation analysis of an unidirectionally coupled system

Stability and Hopf-bifurcation analysis of an unidirectionally coupled system Electronic Journal of Qualitative Theory of Differential Equations 211, No. 85, 1-17; http://www.math.u-szeged.hu/ejqtde/ Staility and Hopf-ifurcation analysis of an unidirectionally coupled system Gang

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

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y. FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the

More information

The complex non resonant double Hopf degeneracy: An alternative approach

The complex non resonant double Hopf degeneracy: An alternative approach The complex non resonant double Hopf degeneracy: An alternative approach Griselda R. Itovich Dept. of Mathematics, FaEA Universidad Nacional del Comahue Buenos Aires 1400, (8300) Neuquén, Argentina Email:

More information

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD

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

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

Higher Order Averaging : periodic solutions, linear systems and an application

Higher Order Averaging : periodic solutions, linear systems and an application Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,

More information

Shilnikov bifurcations in the Hopf-zero singularity

Shilnikov bifurcations in the Hopf-zero singularity Shilnikov bifurcations in the Hopf-zero singularity Geometry and Dynamics in interaction Inma Baldomá, Oriol Castejón, Santiago Ibáñez, Tere M-Seara Observatoire de Paris, 15-17 December 2017, Paris Tere

More information

Period-doubling cascades of a Silnikov equation

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,

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

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e. Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may

More information

A matching pursuit technique for computing the simplest normal forms of vector fields

A matching pursuit technique for computing the simplest normal forms of vector fields Journal of Symbolic Computation 35 (2003) 59 65 www.elsevier.com/locate/jsc A matching pursuit technique for computing the simplest normal forms of vector fields Pei Yu, Yuan Yuan Department of Applied

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

THE CENTER-FOCUS PROBLEM AND BIFURCATION OF LIMIT CYCLES IN A CLASS OF 7TH-DEGREE POLYNOMIAL SYSTEMS

THE CENTER-FOCUS PROBLEM AND BIFURCATION OF LIMIT CYCLES IN A CLASS OF 7TH-DEGREE POLYNOMIAL SYSTEMS Journal of Applied Analysis and Computation Volume 6, Number 3, August 016, 17 6 Website:http://jaac-online.com/ DOI:10.1194/01605 THE CENTER-FOCUS PROBLEM AND BIFURCATION OF LIMIT CYCLES IN A CLASS OF

More information

Two models for the parametric forcing of a nonlinear oscillator

Two models for the parametric forcing of a nonlinear oscillator Nonlinear Dyn (007) 50:147 160 DOI 10.1007/s11071-006-9148-3 ORIGINAL ARTICLE Two models for the parametric forcing of a nonlinear oscillator Nazha Abouhazim Mohamed Belhaq Richard H. Rand Received: 3

More information

arxiv: v1 [nlin.cd] 20 Jul 2010

arxiv: v1 [nlin.cd] 20 Jul 2010 Invariant manifolds of the Bonhoeffer-van der Pol oscillator arxiv:1007.3375v1 [nlin.cd] 20 Jul 2010 R. Benítez 1, V. J. Bolós 2 1 Departamento de Matemáticas, Centro Universitario de Plasencia, Universidad

More information