Hopf-zero bifurcation of Oregonator oscillator with delay
|
|
- Erica Cole
- 5 years ago
- Views:
Transcription
1 Cai et al. Advances in Difference Equations 08 08:438 R E S E A R C H Open Access Hopf-zero bifurcation of Oregonator oscillator with delay Yuting Cai Liqin Liu * and Chunrui Zhang * Correspondence: fymm007@63.com Department of Mathematics Northeast Forestry University Heilongjiang China Abstract In this paper we study the Hopf-zero bifurcation of Oregonator oscillator with delay. The interaction coefficient and time delay are taken as two bifurcation parameters. Firstly we get the normal form by performing a center manifold reduction and using the normal form theory developed by Faria and Magalhães. Secondly we obtain a critical value to predict the bifurcation diagrams and phase portraits. Under some conditions saddle-node bifurcation and pitchfork bifurcation occur along M and N respectively; Hopf bifurcation and heteroclinic bifurcation occur along H and S respectively. Finally we use numerical simulations to support theoretical analysis. MSC: 34K8; 35B3 Keywords: Oregonator model; Delay; Hopf-zero bifurcation; Normal form Introduction An oscillatory chemical reaction refers to the reaction system of certain antileather concentration showing relatively stable cyclical changes. In 9 Bray achieved a liquid oscillatory reaction in the experiment. In 964 Zhabotinsky reported some other oscillatory reactions of this nature [4 5 ]. Until the 970s Field Koros and Noyes proposed the Oregonator model based on an in-depth study of the BZ reaction. According to Field and Noyes the BZ reaction is simplified. In 979 Tyson assumed that the concentration of the reactants A =[BrO 3 ]andb =[BrCHCOOH ] is independent of time. Therefore we can give the following reactant concentration equation: dp dt = k 3AQ k PQ + k 5 AP k 4 P dq = k dt 3 AQ k PQ + fk 0BW dw =k dt 5 AP k 0 BW P =[HBrO ] Q =[Br]W = [CeIV]. We will make the following changes in this system: x = αp y = βq z = γ W t = δt α k 4 k 3 A 06 mol/l β = k k 3 A 07 mol/l The Authors 08. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License which permits unrestricted use distribution and reproduction in any medium provided you give appropriate credit to the original authors and the source provide a link to the Creative Commons license and indicate if changes were made.
2 Cai et al. Advances in Difference Equations 08 08:438 Page of γ k 4k 5 B k 3 A 3 0 mol/l δ = k 5 B s. So we can obtain another form of the Oregonator oscillator the so-called Tyson-type oscillator: ε dx = qy xy + x x dt δ dy dt = qy xy + h z dz dt = x z ε = k 5B k 3 A 4 0 δ = ε α β 0 6 q = k k 4 k k h =h. Because δ is much smaller than ε the second formula in the system can be approximated as y h z q + x. This yields a simplified two-dimensional Oregonator model [7]withrespecttox and z: ε dx dt = x x h z x q x+q dz dt = x z. x =[HBrO ] z =CeIV. When electric current is applied the catalyst CeIV is perturbed and other species are not affected see [5]. Since the perturbation term is introduced in the equation dz dt = x z we rewrite this equation in the following form: dz = x z + kzt τ. dt We consider the Oregonator model with delay: ε dx dt = x x h z x q x+q dz dt = x z + kzt τ. ε =4 0 δ =4 0 4 q =8 0 4 andh 0 is an adjustable parameter. Nowadays many scholars study the Hopf bifurcation or Hopf-zero bifurcation in delay differential equations and some results have been obtained see [ ]. However to the best of our knowledge there are no studies on the Hopf-zero bifurcation of Oregonator oscillator with time delay. Therefore it is the far-reaching significance to research the Hopf-zero bifurcation of Oregonator model. The remainder of the paper is organized as follows. In Sect. we provide stability and conditions of existence of the Hopf-zero bifurcation by taking the interaction coefficient and delay as two parameters. In Sect. 3 we use the center manifold theory and normal
3 Cai et al. Advances in Difference Equations 08 08:438 Page 3 of form method [6 3] to investigate the Hopf-zero bifurcation with original parameters. In Sect. 4 we give several numerical simulations to support the analytic results. Finally we draw the conclusion in Sect. 5. Stability and existence of Hopf-zero bifurcation Let x z be an equilibrium point of system.. Obviously x x h z x q z = x k. x+q =0 Then we have x x h k xx q = 0.. x + q There are three roots x =0x = x + andx = x of Eq.. x ± = h k q ± h k q +4q + h k.. Therefore we obtain that system. has three steady-state solutions x z 0 0 and x + z +. Obviously there is a unique positive steady state. Theorem. For any ε >0q >0and h >0x + z + is the unique positive steady state of system.. Proof Let Hx=x x h k x x q x+q.from.and.wegetthatonlyx + satisfies Hx= 0 Hx >0for0<x < x + andhx <0forx > x +.FurthermorewehaveHx + =0and H x + <0. We further mainly study the dynamics of the equilibrium point x + z +. If the characteristic equation of system. has a simple pair of purely imaginary eigenvalues ±iω asimple root 0 and all other roots of the characteristic equation have negative real parts then the Hopf-zero bifurcation will occur. Let x = x x + and z = z z +.Thenwecanvary. as the following equivalent system: dx dt = ε x + x + x x + h z + z + x+x + q x+x + +q dz dt = x z + kzt τ..3 The linearization equation of system.3 at 0 0 is dx dt = a x + a z dz dt = x z + kzt τ.4
4 Cai et al. Advances in Difference Equations 08 08:438 Page 4 of a = ε qh z + q+x +.4is + x + anda = ε qh h x + q+x +. The characteristic equation of system Eλ=λ b λ b kλe λτ + ka e λτ = 0.5 b = a andb = a + a. If λ =0istherootofEq..5 then ka b =0.Ifτ =0thensystem.5becomes Eλ=λ b + kλ =0. Therefore we obtain that if τ =0andb + k < 0 then excluding a single zero eigenvalue all the roots of Eq..5havenegativerealparts. Next we consider the case of τ 0.Letiωwith ω >0bearootofλ b λ b kλe λτ + ka e λτ =0.Then ω b iω b kiωe iωτ + ka e iωτ =0. Separating the real and imaginary parts we get ω b = kω sin ωτ ka cos ωτ b ω = kω cos ωτ + ka sin ωτ..6 It follows that ω satisfies ω 4 + b + b k ω + b k a = 0..7 Suppose m = ω and denote u =b + b k r = b k a.theneq..7becomes m + u m + r = 0..8 Following [7] we consider the following cases: B r <0. Then we find that system.5 has a unique positive root m = u + u 4r. B r >0 u >0. Then system.5 has no positive root. B3 r >0 u <0. In this case if system.5 has real positive roots then k is very large and h is infinitely close to one which is a contradiction. Theorem. For the quadratic Eq..8 we have: i if r <0 then Eq..5 has a unique positive root m = u + u 4r. ii if r >0 then Eq..5 has no positive root. Suppose that Eq..8 has positive roots. Without loss of generality we assume that it has a positive root defined by m. ThenEq..7 hasapositiverootω andω must satisfy
5 Cai et al. Advances in Difference Equations 08 08:438 Page 5 of the equation ω + a b ω 3 + ωa kω + a + + a kω + a =. According to system.6 we obtain Denote cosωτ= ω + a b kω + a sinωτ= ω3 + ωa + a kω + a. ω τ j = arccos β +jπ α 0 ω π arccos β +jπ α 0 β = ω +a b α kω +a = ω3 +ωa +a j =0...Thensystem.5hasapairofpurely kω +a imaginary roots ±iω with τ = τ j andτ = τ j j =0...satisfytheequationsinωτ>0. We get k <0whenB holds and then τ j = ω arccos β +jπ j {0...}. We obtain the transversality conditions as follows. Theorem.3 If r <0then d Re{λτ j} dτ 0. Proof Substituting λτ τ = τ j intoeq..5 we get So dλ dτ = λka e λτ kλ e λτ λ b λ + τkλe λτ ke λτ τa λe. λτ dλ = τk a dτ ka λ λa λ + b e λτ ka λ. Consequently we obtain dre λτ dτ { τk a = Re τ=τ j ka λ λa λ + b e λτ } ka λ = τa k a ka + ω + ω ω 4 + a + b ω + b ω k a + 0. ω Theorem.4 If ka = b b + k <0and r <0then for τ = τ j j =0...system. undergoes a Hopf-zero bifurcation at equilibrium x + z +. τ=τ j
6 Cai et al. Advances in Difference Equations 08 08:438 Page 6 of 3 Normal form for Hopf-zero bifurcation In this section we use the center manifold theory and normal form method [6 3] to study Hopf-zero bifurcations. The normal form of a Hopf-zero bifurcation for a general delay-differential equations has been given in the following two papers: one is for a saddlenode-hopf bifurcation [] and theother is fora steady-state Hopf bifurcation [0]. After scaling t t/τ system.3becomes dx dt = τ ε x + x + x x + h z + z + x+x + q x+x + +q dz dt = τx z + kzt. 3. Let τ = τ +μ k =+ a a +μ μ and μ are bifurcation parameters. Then system 3.canbewrittenas dx dt =τ + μ a x + a z + M dz dt =τ 3. + μ x z ++ a a + μ zt M = ε [ h z + q x + + ]x + qh q+x + ε q+x + xz + qh x + x 3 + ε x + +q 4 ε qh x + +q 3 x z. Choose the phase space C = C[ 0]; R 4 with supremum norm and define X t C by X t θ=xt + θ τ θ 0 and X t = sup X t θ.thensystem3.becomes Ẋt=LμX t + FX t μ 3.3 a x + a z LμX t =τ + μ x z ++ a a + μ zt and τ + μ M FX t μ= 0 Lμϕ = 0 dηθ μϕξ dξ for ϕ [ 0] R4 0 θ =0 ηθ μ= τ + μ A θ 0 τ + μ A + B θ = with A = a a 0 0 and B = 0 + a. a Consider the linear system Ẋt=L0X t.
7 Cai et al. Advances in Difference Equations 08 08:438 Page 7 of Between C and C = C[0 τ] C n the bilinear form is defined by ψs ϕθ = ψ0ϕ0 0 = ψ0ϕ0 = ψ0ϕ0 θ 0 0 θ 0 0 ψξ θ dηθ0ϕξ dξ ψξ θ d Aϕθ+Bϕθ + ϕξ dξ ψξ +Bϕξ dξ ψ C ϕ C ϕθ=ϕ θ ϕ θ ϕ 3 θ C ψs= ψ s ψ s C. ψ 3 s We know that L0 has a pair of purely imaginary eigenvalues ±iω ω >0asimple0 and all other eigenvalues have negative real parts. Let Λ = {iω iω0}letp be the generalized eigenspace associated with ΛandletP is the space adjoint with P.ThenC can be decomposed as C = P QQ = {ϕ C :ψ ϕ = 0 for all ψ P }.Wecanchoose the bases Φ and Ψ for P and P such that Ψ s Φθ = I Φ = ΦJ and Ψ = JΨ J = diagiω iω0. We calculate ΦθandΨ s asfollows: Φθ= a iω a e iwτ θ a iω a e iwτ θ a e iwτ θ e iwτ θ a and Ψ s= D iω a e iwτ s D iω a e iwτ s D e iwτ s D e iwτ s D a D D = a iω a ++τ D = a + a + a τ b. + a a Let usenlarge the space C to the following space: { } BC = ϕ is a continuous function on [ 0 and lim ϕθexists. θ 0 Itselementscanbewrittenasψ = ϕ + X 0 α with ϕ C α C n and 0 θ [ 0 X 0 θ= I θ =0. In BC system3.3 varies an abstract ODE: d dt X t = Au + X 0 Fu μ 3.4
8 Cai et al. Advances in Difference Equations 08 08:438 Page 8 of u CandA is defined by A : C BC Au = u + X 0 [ L0u u0 ] and Fu μ= [ Lμ L 0 ] u + Fu μ. Then the enlarged phase space BC can be decomposed as BC = P Ker π.letx t = Φxt+ ỹθ xt=x x x 3 T namely xθ= a iω a e iwτθ x + a iω a e iwτθ x a x 3 + y θ zθ=e iwτθ x + e iwτθ x + a x 3 + y θ. Let ψ ψ Ψ 0 = ψ ψ = ψ 3 ψ 3 D iω a D D iω a D D a D. Equation 3.4can be expressed as ẋ = Jz + Ψ 0 F Φz + ỹθ μ ỹ = A Q ỹ +I πy 0 F Φz + ỹ0 μ 3.5 ỹθ Q := Q C Ker πanda Q is the restriction of A as an operator from Q to the Banach space Ker π. System 3.5canberewrittenas ẋ = Jx +! f x y μ+ 3! f 3 x y μ+h.o.t. ẏ = A Q y +! f x y μ+ 3! f 3 x y μ+h.o.t. f f 3 ψ F x y μ+ψ F x y μ x y μ= ψ F x y μ+ψ F x y μ ψ 3 F x y μ+ψ 3F x y μ ψ F3 x y μ+ψ F3 x y μ x y μ= ψ F3 x y μ+ψ F3 x y μ ψ 3 F3 x y μ+ψ 3F3 x y μ
9 Cai et al. Advances in Difference Equations 08 08:438 Page 9 of with [ F = μ a a a x + x a x 3 + y 0 + a x + x + a x 3 + y 0 ] iω a iω a [ [ ] hz+ q x + + a a + τ x ε q + x + + x a x 3 + y 0 iω a iω a + qh a a x x + x a x 3 + y 0 + x + a x 3 + y 0 ] ε q + x + iω a iω a F = μ a a x + x a x 3 + y 0 μ x + x + a x 3 + y 0 iω a iω a + + a μ e iwτ x + e iwτ x + a x 3 + y a + τ μ e iwτ x + e iwτ x + a x 3 + y 3! F 3 = τ qhx + ε q + x + 4 F 3 =0. + τ ε qh q + x + 3 a x + iω a a x + iω a x + x + a x 3 + y 0 According to [5] ImM c is spanned by a iω a x a x 3 + y 0 3 a iω a x a x 3 + y 0 { z e z z 3 e z μ i e μ μ e z z e z μ i e z z 3 e 3 z 3 μ i e 3 } i = with e = 0 0 T e = 0 0 T e 3 = 0 0 T. ImM 3 c is spanned by { z 3 e z z z 3 e z z e z z 3e z z 3e 3 z z3 e } 3. Then we get g x0μ=proj ImM cf x0μ=proj S f x0μ+o μ g3 x0μ=proj ImM3 c f 3 x0μ=proj S f 3 x00+o μ x + μ x S and S are spanned respectively by and z μ i e z μ i e z 3 μ i e 3 i = z 3 e z z z 3 e z z e z z 3e z z 3e 3 z z 3 e 3. System 3.5 canbetransformedonthecentermanifoldinthefollowingnormalform: ẋ = Jx + g x0μ+ 6 g 3 x0μ+h.o.t. 3.6
10 Cai et al. Advances in Difference Equations 08 08:438 Page 0 of We need to compute g x0μandg 3 x0μin3.6. We can compute g x0μ: g x0μ= Proj S f x0μ+ϑ u a μ + a μ x + a 3 x x 3 = a μ + a μ x + a 3 x x 3 + ϑ u a μ + a μ x 3 + a 3 x x + a 4 x 3 a = D τ e iωτ a a a = D iω a + a + iω a [ D τ hz + q x + + a 3 = iω a ε q + x + a = a τ D a = a b D [ τ hz + q x + + a a 3 = D ε q + x + a ω a 4 = τ D ε + a a a iω a τ ε [ hz+ q x + + q + x + a + qh q + x + a a e iωτ qhτ εq + x + qh q + x + ]. ] a a a iω a a iω a + a iω a ] Next we compute g 3 x0μ: 6 g 3 x0μ= 6 Proj KerM f 3 x0μ = 6 Proj S f 3 x00+ϑ x u + x μ = 6 Proj S f 3 x Proj S [ Dx f + ϑ x u + x μ. x00u x0+ D y f x00u x0 ] We can get Proj S f 3 x00.since 6 f 3 x00= τ ε qhx D + a x + +q 4 iω a x + a iω a x a x qh a x + +q 3 iω a x + a iω a x a x 3 qhx D + a x + +q 4 iω a x + a iω a x a x qh a x + +q 3 iω a x + a iω a x a x 3 x + x + a x 3 qhx a D + a x + +q 4 iω a x + a iω a x a x qh a iω a x + a iω a x a x 3 x + x + a x 3 x + +q 3
11 Cai et al. Advances in Difference Equations 08 08:438 Page of we have 6 Proj S f3 x00= b x x + b x x 3 b x x + b x x 3 + ϑ u b x x x 3 + b x 3 3 b = τ ε D qhx + a a 3 x + + q 4 ω + a + ω + a iω a + τ ε D qh a a + x + + q 3 iω a ω + a b = τ ε D qhx + a 3 τ x + + q 4 iω a ε D qh a a x + + q 3 iω a b = τ ε a qhx + 6a 3 D x + + q 4 ω + a + τ ε a qh a a D x + + q 3 ω + a b = τ ε a qhx + D x + + q 4 a3 + τ ε a qh D x + + q a a a iω a Next we compte Proj S [D x f x00u x0]. From [5] weknowthatbecausej is a diagonal matrix the operators Mj j are defined in Vj 5C3 so we have a diagonal representation relative to the canonical basis {μ p x q e k : k = 3p N 0 q N3 5 0 p + q = j} of Vj C3 e = 0 0 T e = 00 T e 3 = 0 0 T. Clearly we get Mj μ p x q e k = iω q q + k μ p x q e k k = μ p x q e 3 = iωq q μ p x q e 3 p + q = j. M j So Ker Mj { = span μ p x q e k :q λ=λ k k =3p N 0 q N3 0 p + q = j} with λ =λ λ λ 3 =iω iω 0. The elements of the canonical basis of V 5 C3 are μ μ e μ i e μ i x e μ i x e μ i x 3 e x x e x x 3 e x x 3 e x e x e x 3 e μ μ e μ i e μ i x e μ i x e μ i x 3 e x x e x x 3 e x x 3 e x e x e x 3 e μ μ e 3 μ i e 3 μ i x e 3 μ i x e 3 μ i x 3 e 3 x x e 3 x x 3 e 3 x x 3 e 3 x e 3 x e 3 x 3 e 3 i = the images of which under iω M are μ μ e μ i e 0 μ i x e μ i x 3 e x x e 0 x x 3 e x e 3x e x 3 e μ μ e μ i e μ i x e 0μ i x 3 e x x e x x 3 e 03x e x e x 3 e 0 0 μ i x e 3 μ i x e 3 00x x 3 e 3 x x 3 e 3 x e 3 x e 30 i =.
12 Cai et al. Advances in Difference Equations 08 08:438 Page of Hence U x0=u x μ μ=0 = M ProjImM f x00 D h[a hx + a hx a x 3 m + n a hx + a hx a x 3 x + x + a x 3 ] D h[a hx + a hx a x 3 m + n a hx + a hx a x 3 x + x + a x 3 ] D [a hx + a hx a x 3 m + n a hx + a hx a x 3 x + x + a x 3 ] D h[ma h x a hhx x + a hx x 3 3 a h a x 3 + na hx a h + hx x 3 a hx a a h +x x 3 + a a x 3 ] D h[m 3 a h x a h x + a x 3 +a hhx x a hx x 3 + n 3 a hx + a a h x x 3 a hx + a a hx x 3 a a x 3 + a h + hx x ] D [m a h x a h x a hx x 3 +a hx x 3 + n a hx + a a hx x 3 a hx a a hx x 3 a x x 3 + a x x 3 = M ProjImM τ ε = τ iωε m = hz +q x + + n = qh q+x + q+x + andh = iω a. Therefore we obtain 4 Proj [ S Dx f x00 U x0] = c x x + c x x 3 c x x + c x x 3 c x x x 3 + c x 3 3 c = τ [ D iωε h m a 4 h3 h + D h mn a 3 h3 + D h mn 3a 3 h h 6 + D h n a h h h+d D hhm 3 a4 h h 4 + D D hhmn 3 a3 h h a3 hh + a h + a hh + D D hhn a hh + D D hm a 4 h h + D D hmn a 3 hh D D hmn a a 3 h h a3 h D D hn a a hh a h ]
13 Cai et al. Advances in Difference Equations 08 08:438 Page 3 of c = τ [ D iωε h m a 4 h +D h mna 3 a h + h + D h n a a h + D D hhm 4a 4 hh + D D hhmn 4a a 3 hh +a3 h +a3 h + D D hhn 3 a a h + h+ a a hh + + D D hm 4a 4 h D D hmn 6a a 3 h a3 +D D hn a a h a a ] c = τ iωε [ D D hm 6a 4 h h D D hmn 3a a 3 h h +6a 3 hh +3a3 h D D hn a a hh +a h a a h + a h + D D hm 6a 4 hh D D hmn 3a a 3 hh 3a 3 h 6a 3 hh D D hn a a hh a h + a a h a h] c = τ [ D D iωε hm a 4 h D D hmn a 3 3a a 3 h D D hn a a h a a D D hm a 4 h D D hmn 3a a 3 h a3 + D D hn a a h a a ]. Finally we can compute Proj S [D y f x00u x0].defineh = hxθ=u x0and write h θ hθ= = h h 00 x θ + h 00x + h 00x 3 + h 0x x + h 0 x x 3 + h 0 x x 3 h 00 h 00 h 00 h 0 h 0 h 0 Q.M hx=f x00 decides the coefficients of hwhichisequivalentto D x hjx A Q h=i πx 0 F Φx0. We use the definitions of A Q and π to obtain ḣ D x hjx = ΦθΨ 0F Φx0 h0 Lh = F Φx0 ḣ is the derivative of hθwithrespecttoθ.let F Φx0=A 00 x + A 00x + A 00x 3 + A 0x x + A 0 x x 3 + A 0 x x 3 A ijk C0 i j k i+j+k =. We can compare the coefficients of x x x 3 x x x x 3 x x 3 andwegetthath 00 = h 00 h 0 = h 0 and the following differential equations are satisfied by h 00 h 0 h 0 h 00 respectively: ḣ 00 iωτ h 00 = ΦθΨ 0A 00 ḣ 00 0 Lh 00 =A
14 Cai et al. Advances in Difference Equations 08 08:438 Page 4 of ḣ 0 iωτ h 0 = ΦθΨ 0A 0 ḣ 0 0 Lh 0 =A 0 ḣ 0 = ΦθΨ 0A 0 ḣ 0 0 Lh 0 =A 0 ḣ 00 = ΦθΨ 0A 00 ḣ 00 0 Lh 00 =A τ a hma h + n τ a hma h + n A 00 = A 00 = 0 0 τ a ma na τ a hmh h + nh + hn A 00 = A 0 = 0 0 τ a ma h + nha n τ a ma h + a n h n A 0 = A 0 =. 0 0 Since we have τε mx + nxz F u t 0= 0 f x y0=ψ 0F Φx + y0 which gives = τ ε = τ ε D h D D h D D a D ma hx + a hx a x 3 + y 0 + na hx + a hx a x 3 + y 0 x + x + a x 3 + y 0 0 D h[ma hx + a hx a x 3 + y 0 + na hx + a hx a x 3 + y 0x + x + a x 3 + y 0] D h[ma hx + a hx a x 3 + y 0 + na hx + a hx a x 3 + y 0x + x + a x 3 + y 0] D [ma hx + a hx a x 3 + y 0 + na hx + a hx a x 3 + y 0x + x + a x 3 + y 0] D h[ma hx + a hx a x 3 + na hx + a hx a x 3 + x + x + a x 3 ]h 0 4 D yf y=0μ=0h= τ D h[ma hx + a hx a x 3 ε + na hx + a hx a x 3 + x + x + a x 3 ]h 0. D [ma hx + a hx a x 3 + na hx + a hx a x 3 + x + x + a x 3 ]h 0
15 Cai et al. Advances in Difference Equations 08 08:438 Page 5 of Thus d x 4 Proj S D y f y=0μ=0u = x + d x x 3 d x x + d x x 3 d x x x 3 + d x x 3 3 d = τ ε D hm + na h h h h τ ε D hn h h 00 0 d = τ ε D hm + na h h 00 0 h τ ε D hn h a h 0 0 d = τ ε D hm + na h h h h 0 0 h τ ε D hn a h h h 0 0 d = D hm + na h D hna h So we can obtain b + c + d x 6 g 3 x0μ= x +b + c + d x x 3 b + c + d x x +b + c + d x x 3 + ϑ x μ + x μ. b + c + d x x x 3 +b + c + d x 3 3 Therefore on the center manifold the system ẋ = Jx + g x0μ + 6 g 3 x0μ +h.o.t. becomes ẋ =a μ + a μ x + a 3 x x 3 +b + c + d x x +b + c + d x x 3 + h.o.t. ẋ =a μ + a μ x + a 3 x x 3 +b + c + d x x +b + c + d x x 3 + h.o.t. ẋ 3 =a μ + a μ x 3 + a 3 x x + a 4 x 3 +b + c + d x x x 3 +b + c + d x h.o.t. 3. By changing variables x = ρ iρ x = ρ + iρ x 3 = ρ 3 and introducing the cylindrical coordinates ρ = r cos θ ρ = r sin θ ρ 3 = γ r >0system3.becomes ṙ = α μr + β rγ + β 30 r 3 + β rγ + h.o.t. γ = α μγ + m 0 r + m 0 γ + m r γ + m 03 γ 3 + h.o.t. θ = ω +Im[a ]μ + Im[a ]μ α μ=re[a ]μ + Re[a ]μ β = Re[a 3 ] β 30 = Re[b + c + d ] β = Re[b + c + d ] α μ=a μ m 0 = a 3 m 0 = a 4 m = b + c + d m 03 = b + c + d.
16 Cai et al. Advances in Difference Equations 08 08:438 Page 6 of Therefore we can get the system in the plane r γ : ṙ = α μr + β rγ + β 30 r 3 + β rγ + h.o.t. ζ = α μγ + m 0 r + m 0 γ + m r γ + m 03 γ 3 + h.o.t. 3. From [5]we know that Eq.3.becomes ṙ =α μ+β δ + β δ r +β +β δrγ + β 30 r 3 + β rγ γ =α μδ + m 0 δ + m 03 δ 3 +α μ+m 0 δ +3m 03 δ γ +m 0 + m δr +m 0 +3m 03 δγ + m r γ + m 03 γ Choose δ = δμsuchthat α μ+m 0 δ +3m 03 δ =0. To simplify the above system we only discuss the case of m 0 0m Clearly for small α μ the equation has two real roots. We take 3m δ = 03 [ m 0 + m 0 3m 03α μ] if m 0 >0 3m 03 [ m 0 m 0 3m 03α μ] if m 0 <0. Then δ = δμ is differentiable at μ = 0 andδ0 = 0. Define k = α μ+β δ + β δ k = α μδ + m 0 δ + m 03 δ 3 a = β +β δ b = m 0 + m δ c = m 0 +3m 03 δ and choose x = r y = γ.theneq.3.3becomes ẋ = k x + axy + β 30 x 3 + β xy ẏ = k + bx + cy + γ x y + γ 03 y Let x c x y b y t c b t and η = k c b η = k c b. Then system 3.4becomes ẋ = η x + Bxy + d x 3 + d xy ẏ = η + ηx y y + d 3 x y + d 4 y B = a c 0 η = sgnbc
17 Cai et al. Advances in Difference Equations 08 08:438 Page 7 of and d = β 30 c c b d b β = d 3 = m c c c b d b m03 4 =. c According to [5] we assume that K 3 = η B + d + B d + ηd 3 +3d 4 0. For small η and η the qualitative behavior of 3.5 near 0 0 is the same as that of the following system see [9]: ẋ = η x + Bxy + xy ẏ = η x y. 3.6 In Eq. 3.6 there are two trivial equilibrium points E =0± η η > 0 and two nontrivial equilibrium points E 34 = B B ± B 4η +η + η B ± B 4η. In [3]and[6] wecan find thecompletebifurcationdiagrams ofsystem 3.3. Herewe list some of them. Theorem 3. a If B <0 then the bifurcation diagram of system 3.3 consists of the origin and the following curves: M = { η η :η =0η 0 } { N = η η :η = } B η + ϑ η 3 η 0. Along M and N a saddle-node bifurcation and pitchfork bifurcation occur respectively. System 3.3 has no periodic orbits. Moreover if η η is in the region between M and N then the solution of system 3.3 goes asymptotically to one of the equilibrium points E E and E 3. b If B >0 then the bifurcation diagram of system 3.3 consists of the origin the curves M and N and the following curves: H = { η η :η =0η >0 } { S = η η :η = B 3B + η + ϑ η 3/ } η >0. Along M and N we have exactly the same bifurcation as in a. Along H and S a Hopf bifurcation and a heteroclinic bifurcation occur respectively. If η η lies between the curves H and S then system 3.3 has a unique limit cycle which is unstable and becomes a heteroclinic orbit when η η S. Figures and showaandboftheorem3.respectively.
18 Cai et al. Advances in Difference Equations 08 08:438 Page 8 of Figure When B < 0 the bifurcation diagrams and phase portraits of system 3.3see[8] Figure When B > 0 the bifurcation diagrams and phase portraits of system 3.3see[8] 4 Numerical simulations In this section we give some examples to explain the theoretical results. Set q =8 0 4 h = 3 k =.5andε =4 0 and consider the following system: dx dt = ε dz dt x u x x hz x+u = x z + kzt τ. By calculations we obtain a = i a = a 3 = i a =.9430 a =.453 a 3 = a 4 =.654 b = i b = i b = i b =.8473 c = i c = i c = 0.07 c = d = i d =
19 Cai et al. Advances in Difference Equations 08 08:438 Page 9 of Figure 3 The system is asymptotically stable around the equilibrium point for μ μ = The green line represents x the red line represents z. Waveform diagram for variable of x z left. Phase diagram for variable x zright Figure 4 The system has an asymptotically stable periodic orbit near τ for μ μ = The green line represents x the red line represents z. Waveform diagram for variable of x z left. Phase diagram for variable x zright i d = d = ; the equilibrium point is and τ = For small μweobtaink Conclusions In this article we have discussed the Hopf-zero bifurcation of Oregonator oscillator with delay. We thoroughly analyze the distribution of the eigenvalues of the corresponding characteristic equation and find some specific conditions ensuring that all the eigenvalues have negative real parts. We also can discover the factors that make system. undergo a Hopf-zero bifurcation at equilibrium x + z +. Meanwhile by using the normal form method and the center manifold theorem we have derived the normal form of the reduced system on the center manifold and discussed the Hopf-zero bifurcation with parameters in system.. Besides we have obtained bifurcation diagrams and phase portraits of system 3.3whenB >0andB < 0 respectively. We also note that a saddle-node bifurcation and pitchfork bifurcation occur along M and Nrespectivelyand a Hopf bifurcation and a heteroclinic bifurcation occur along H and S respectively. Finally numerical stimulations see Figure 3 4 and 5 have been given to illustrate the theoretical results. Our work is a further study of the Oregonator oscillator which will be useful in the research of the complex phenomenon caused by high codimensional bifurcation of a delaydifferential equation.
20 Cai et al. Advances in Difference Equations 08 08:438 Page 0 of Figure 5 Waveform diagram for variable of x z. The first two figures show that when τ =.50<τ 0 the system is stable around the equilibrium point. When τ =.80>τ 0 the next two figures show that the system is unstable around the equilibrium point Acknowledgements The authors wish to express their gratitude to the editors and the reviewers for the helpful comments. Funding This research is supported by the Heilongjiang Provincial Natural Science Foundation No. A0506. Competing interests The authors declare that they have no competing interests. Authors contributions The idea of this research was introduced by YC LL and CZ. All authors contributed to the main results and numerical simulations and approved the final manuscript. Publisher s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Received: 9 September 08 Accepted: 9 November 08 References. Cao X. Song Y. Zhang T.: Hopf bifurcation and delay-induced Turing instability in a diffusive lac operon model. Int. J. Bifurc. Chaos Chang X. Wei J.: Hopf bifurcation and optimal control in a diffusive predator prey system with time delay and prey harvesting. Nonlinear Anal. Model. Control Chow S.N. Li C. Wang D.: Normal forms and bifurcation of planar vector fields. Cambridge Ding W. Liao X.F. Dong T.: Hopf bifurcation in a love-triangle model with time delays. Neurocomputing Ding Y. Jiang W. Yu P.: Hopf-zero bifurcation in a generalized Gopalsamy neural network model. Nonlinear Dyn Faria T. Magalhaes L.T.: Normal forms for retarded functional differential equation with parameters and applications to Hopf bifurcation. J. Differ. Equ Gazor M. Sadri N.: Bifurcation control and universal unfolding for Hopf-zero singularities with leading solenoidal terms. SIAM J. Appl. Dyn. Syst Guo Y. Jiang W.: Hopf bifurcation in two groups of delay-coupled Kuramoto oscillators. Int. J. Bifurc. Chaos
21 Cai et al. Advances in Difference Equations 08 08:438 Page of 9. Isaac A.G. Jaume L. Susanna M.: On the periodic orbit bifurcating from a zero Hopf bifurcation in systems with two slow and one fast variables. Appl. Math. Comput Jaume L. Zhang X.: On the Hopf-zero bifurcation of the Michelson system. Nonlinear Anal. Real World Appl Jiang H. Jiang J. Song Y.: Normal form of saddle-node-hopf bifurcation in retarded functional differential equations and applications. Int. J. Bifurc. Chaos Jiang W. Wang H.: Hopf-transcritical bifurcation in retarded functional differential equations. Nonlinear Anal Jiang W. Wang J.: Hopf-zero bifurcation of a delayed predator prey model with dormancy of predators. J. Appl. Anal. Comput Jimenez-Prieto R. Silva M. Perez-Bendito D.: Analytical assessment of the oscillating chemical reactions by use chemiluminescence detection. Talanta Kaplan D.T. Glass L.: In Understanding Nonlinear Dynamics. Springer Berlin Kuznetsov Yu.: Elements of Applied Bifurcation Theory 3rd edn. Springer Berlin Liu Z. Yuan R.: Zero-Hopf bifurcation for an infection-age structured epidemic model with a nonlinear incidence rate. Sci. China Math Marsden J.E. Sirovich L. John F.: Nonlinear Oscillations Dynamical System and Bifurcations of Vector Fields. Nonlinear Mathematical Sciences vol Rodrigo D. Jaume L.: Zero-Hopf bifurcation in a Chua system. Nonlinear Anal. Real World Appl Song Y. Jiang J.: Steady-state Hopf and steady-state-hopf bifurcations in delay differential equations with applications to a damped harmonic oscillator with delay feedback. Int. J. Bifurc. Chaos Takens F.: Lecture Notes in Math vol. 3 pp Springer Berlin 98. Tian X. Xu R.: The Kaldor Kalecki stochastic model of business cycle. Nonlinear Anal. Model. Control Wang H. Wang J.: Hopf-pitchfork bifurcation in a two-neuron system with discrete and distributed delays. Math. Methods Appl. Sci Wang Y. Wang H. Jiang W.: Hopf-transcritical bifurcation in toxic phytoplankton zooplankton model with delay. J. Math. Anal. Appl Wu X. Wang L.: Zero-Hopf bifurcation for van der Pol s oscillator with delayed feedback. J. Comput. Appl. Math Wu X. Wang L.: Zero-Hopf bifurcation analysis in delayed differential equations with two delays. J. Franklin Inst Wu X. Zhang C.R.: Dynamic properties of the coupled Oregonator model with delay. Nonlinear Anal. Model. Control Yang J. Zhao L.: Bifurcation analysis and chaos control of the modified Chua s circuit system. Chaos Solitons Fractals Zhao H. Lin Y. Dai Y.: Hopf bifurcation and hidden attractors of a delay-coupled Duffing oscillator. Int. J. Bifurc. Chaos
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 informationHOPF-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 informationDynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response
Liu Advances in Difference Equations 014, 014:314 R E S E A R C H Open Access Dynamical analysis of a delayed predator-prey system with modified Leslie-Gower and Beddington-DeAngelis functional response
More informationSTABILITY AND HOPF BIFURCATION IN A SYMMETRIC LOTKA-VOLTERRA PREDATOR-PREY SYSTEM WITH DELAYS
Electronic Journal of Differential Equations Vol. 013 (013) No. 09 pp. 1 16. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu STABILITY AND HOPF BIFURCATION
More informationHOPF 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 informationContinuous Threshold Policy Harvesting in Predator-Prey Models
Continuous Threshold Policy Harvesting in Predator-Prey Models Jon Bohn and Kaitlin Speer Department of Mathematics, University of Wisconsin - Madison Department of Mathematics, Baylor University July
More informationDOUBLE HOPF BIFURCATION IN DELAYED VAN DER POL DUFFING EQUATION
International Journal of Bifurcation and Chaos, Vol. 3, No. (3) 354 (5 pages) c World Scientific Publishing Company DOI:.4/S8743544 DOUBLE HOPF BIFURCATION IN DELAYED VAN DER POL DUFFING EQUATION YUTING
More informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationHopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models
Vol.37 (SUComS 6), pp.57-66 http://dx.doi.org/.457/astl.6.37. Hopf-Fold Bifurcation Analysis in a Delayed Predator-prey Models Shuang Guo, Xiuli Li, Jian Yu, Xiu Ren Department of Mathematics, Teacher
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,
More informationSTABILITY. 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 informationResearch Article Hopf Bifurcation Analysis and Anticontrol of Hopf Circles of the Rössler-Like System
Abstract and Applied Analysis Volume, Article ID 3487, 6 pages doi:.55//3487 Research Article Hopf Bifurcation Analysis and Anticontrol of Hopf Circles of the Rössler-Like System Ranchao Wu and Xiang Li
More informationStability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks
International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationResearch Article On the Stability Property of the Infection-Free Equilibrium of a Viral Infection Model
Hindawi Publishing Corporation Discrete Dynamics in Nature and Society Volume, Article ID 644, 9 pages doi:.55//644 Research Article On the Stability Property of the Infection-Free Equilibrium of a Viral
More informationPD control at Neimark Sacker bifurcations in a Mackey Glass system
Wang Advances in Difference Equations (2018) 2018:411 https://doi.org/10.1186/s13662-018-1864-8 R E S E A R C H Open Access PD control at Neimark Sacker bifurcations in a Mackey Glass system Yuanyuan Wang
More informationHOPF 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 informationDirection 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 informationSTABILITY 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 informationExistence 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{ =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 informationTWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY. P. Yu 1,2 and M. Han 1
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 3, Number 3, September 2004 pp. 515 526 TWELVE LIMIT CYCLES IN A CUBIC ORDER PLANAR SYSTEM WITH Z 2 -SYMMETRY P. Yu 1,2
More informationB5.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 informationHopf bifurcation analysis for a model of plant virus propagation with two delays
Li et al. Advances in Difference Equations 08 08:59 https://doi.org/0.86/s366-08-74-8 R E S E A R C H Open Access Hopf bifurcation analysis for a model of plant virus propagation with two delays Qinglian
More informationHopf Bifurcation of a Nonlinear System Derived from Lorenz System Using Centre Manifold Approach ABSTRACT. 1. Introduction
Malaysian Journal of Mathematical Sciences 10(S) March : 1-13 (2016) Special Issue: The 10th IMT-GT International Conference on Mathematics, Statistics and its Applications 2014 (ICMSA 2014) MALAYSIAN
More informationResearch Article Stability and Hopf Bifurcation in a Computer Virus Model with Multistate Antivirus
Abstract and Applied Analysis Volume, Article ID 84987, 6 pages doi:.55//84987 Research Article Stability and Hopf Bifurcation in a Computer Virus Model with Multistate Antivirus Tao Dong,, Xiaofeng Liao,
More informationHopf 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 informationCOMPLEX 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 informationEE222 - Spring 16 - Lecture 2 Notes 1
EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued
More informationBIFURCATION 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 informationSTABILITY AND HOPF BIFURCATION OF A MODIFIED DELAY PREDATOR-PREY MODEL WITH STAGE STRUCTURE
Journal of Applied Analysis and Computation Volume 8, Number, April 018, 573 597 Website:http://jaac-online.com/ DOI:10.11948/018.573 STABILITY AND HOPF BIFURCATION OF A MODIFIED DELAY PREDATOR-PREY MODEL
More informationANALYSIS 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 informationResearch Article Adaptive Control of Chaos in Chua s Circuit
Mathematical Problems in Engineering Volume 2011, Article ID 620946, 14 pages doi:10.1155/2011/620946 Research Article Adaptive Control of Chaos in Chua s Circuit Weiping Guo and Diantong Liu Institute
More informationHopf Bifurcation and Limit Cycle Analysis of the Rikitake System
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(0) No.,pp.-5 Hopf Bifurcation and Limit Cycle Analysis of the Rikitake System Xuedi Wang, Tianyu Yang, Wei Xu Nonlinear
More information7 Two-dimensional bifurcations
7 Two-dimensional bifurcations As in one-dimensional systems: fixed points may be created, destroyed, or change stability as parameters are varied (change of topological equivalence ). In addition closed
More informationNew results on the existences of solutions of the Dirichlet problem with respect to the Schrödinger-prey operator and their applications
Chen Zhang Journal of Inequalities Applications 2017 2017:143 DOI 10.1186/s13660-017-1417-9 R E S E A R C H Open Access New results on the existences of solutions of the Dirichlet problem with respect
More informationDynamical analysis and circuit simulation of a new three-dimensional chaotic system
Dynamical analysis and circuit simulation of a new three-dimensional chaotic system Wang Ai-Yuan( 王爱元 ) a)b) and Ling Zhi-Hao( 凌志浩 ) a) a) Department of Automation, East China University of Science and
More informationTracking Control of a Class of Differential Inclusion Systems via Sliding Mode Technique
International Journal of Automation and Computing (3), June 24, 38-32 DOI: 7/s633-4-793-6 Tracking Control of a Class of Differential Inclusion Systems via Sliding Mode Technique Lei-Po Liu Zhu-Mu Fu Xiao-Na
More informationBIFURCATIONS 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 informationResearch Article Stability Switches and Hopf Bifurcations in a Second-Order Complex Delay Equation
Hindawi Mathematical Problems in Engineering Volume 017, Article ID 679879, 4 pages https://doi.org/10.1155/017/679879 Research Article Stability Switches and Hopf Bifurcations in a Second-Order Complex
More informationSIMPLE 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 informationHopf 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 informationBifurcation of Fixed Points
Bifurcation of Fixed Points CDS140B Lecturer: Wang Sang Koon Winter, 2005 1 Introduction ẏ = g(y, λ). where y R n, λ R p. Suppose it has a fixed point at (y 0, λ 0 ), i.e., g(y 0, λ 0 ) = 0. Two Questions:
More informationTWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations
TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1
More informationTime-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 informationGlobal dynamics of two systems of exponential difference equations by Lyapunov function
Khan Advances in Difference Equations 2014 2014:297 R E S E A R C H Open Access Global dynamics of two systems of exponential difference equations by Lyapunov function Abdul Qadeer Khan * * Correspondence:
More informationHOPF 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 informationA Comparison of Two Predator-Prey Models with Holling s Type I Functional Response
A Comparison of Two Predator-Prey Models with Holling s Type I Functional Response ** Joint work with Mark Kot at the University of Washington ** Mathematical Biosciences 1 (8) 161-179 Presented by Gunog
More informationCONTROLLING 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 informationA Unified Lorenz-Like System and Its Tracking Control
Commun. Theor. Phys. 63 (2015) 317 324 Vol. 63, No. 3, March 1, 2015 A Unified Lorenz-Like System and Its Tracking Control LI Chun-Lai ( ) 1, and ZHAO Yi-Bo ( ) 2,3 1 College of Physics and Electronics,
More informationLYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR
LYAPUNOV EXPONENTS AND STABILITY FOR THE STOCHASTIC DUFFING-VAN DER POL OSCILLATOR Peter H. Baxendale Department of Mathematics University of Southern California Los Angeles, CA 90089-3 USA baxendal@math.usc.edu
More informationBoundedness of solutions to a retarded Liénard equation
Electronic Journal of Qualitative Theory of Differential Equations 21, No. 24, 1-9; http://www.math.u-szeged.hu/ejqtde/ Boundedness of solutions to a retarded Liénard equation Wei Long, Hong-Xia Zhang
More informationAnalysis of the Takens-Bogdanov bifurcation on m parameterized vector fields
Analysis of the Takens-Bogdanov bifurcation on m parameterized vector fields Francisco Armando Carrillo Navarro, Fernando Verduzco G., Joaquín Delgado F. Programa de Doctorado en Ciencias (Matemáticas),
More informationNew Homoclinic and Heteroclinic Solutions for Zakharov System
Commun. Theor. Phys. 58 (2012) 749 753 Vol. 58, No. 5, November 15, 2012 New Homoclinic and Heteroclinic Solutions for Zakharov System WANG Chuan-Jian ( ), 1 DAI Zheng-De (à ), 2, and MU Gui (½ ) 3 1 Department
More informationTHE 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 informationDYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS
International Journal of Bifurcation and Chaos c World Scientific Publishing Company DYNAMICAL SYSTEMS WITH A CODIMENSION-ONE INVARIANT MANIFOLD: THE UNFOLDINGS AND ITS BIFURCATIONS KIE VAN IVANKY SAPUTRA
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationBOUNDARY VALUE PROBLEMS OF A HIGHER ORDER NONLINEAR DIFFERENCE EQUATION
U.P.B. Sci. Bull., Series A, Vol. 79, Iss. 4, 2017 ISSN 1223-7027 BOUNDARY VALUE PROBLEMS OF A HIGHER ORDER NONLINEAR DIFFERENCE EQUATION Lianwu Yang 1 We study a higher order nonlinear difference equation.
More informationA Novel Hyperchaotic System and Its Control
1371371371371378 Journal of Uncertain Systems Vol.3, No., pp.137-144, 009 Online at: www.jus.org.uk A Novel Hyperchaotic System and Its Control Jiang Xu, Gouliang Cai, Song Zheng School of Mathematics
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationAnti-synchronization of a new hyperchaotic system via small-gain theorem
Anti-synchronization of a new hyperchaotic system via small-gain theorem Xiao Jian( ) College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China (Received 8 February 2010; revised
More informationA DELAY-DEPENDENT APPROACH TO DESIGN STATE ESTIMATOR FOR DISCRETE STOCHASTIC RECURRENT NEURAL NETWORK WITH INTERVAL TIME-VARYING DELAYS
ICIC Express Letters ICIC International c 2009 ISSN 1881-80X Volume, Number (A), September 2009 pp. 5 70 A DELAY-DEPENDENT APPROACH TO DESIGN STATE ESTIMATOR FOR DISCRETE STOCHASTIC RECURRENT NEURAL NETWORK
More informationBifurcation in a Discrete Two Patch Logistic Metapopulation Model
Bifurcation in a Discrete Two Patch Logistic Metapopulation Model LI XU Tianjin University of Commerce Department of Statistics Jinba road, Benchen,Tianjin CHINA beifang xl@63.com ZHONGXIANG CHANG Tianjin
More informationADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM
ADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM WITH UNKNOWN PARAMETERS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University
More informationControlling a Novel Chaotic Attractor using Linear Feedback
ISSN 746-7659, England, UK Journal of Information and Computing Science Vol 5, No,, pp 7-4 Controlling a Novel Chaotic Attractor using Linear Feedback Lin Pan,, Daoyun Xu 3, and Wuneng Zhou College of
More informationEntrance 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 informationAdditive 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 informationExistence of Positive Periodic Solutions of Mutualism Systems with Several Delays 1
Advances in Dynamical Systems and Applications. ISSN 973-5321 Volume 1 Number 2 (26), pp. 29 217 c Research India Publications http://www.ripublication.com/adsa.htm Existence of Positive Periodic Solutions
More informationBifurcation and Stability Analysis of a Prey-predator System with a Reserved Area
ISSN 746-733, England, UK World Journal of Modelling and Simulation Vol. 8 ( No. 4, pp. 85-9 Bifurcation and Stability Analysis of a Prey-predator System with a Reserved Area Debasis Mukherjee Department
More informationLecture 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
More informationAPPPHYS217 Tuesday 25 May 2010
APPPHYS7 Tuesday 5 May Our aim today is to take a brief tour of some topics in nonlinear dynamics. Some good references include: [Perko] Lawrence Perko Differential Equations and Dynamical Systems (Springer-Verlag
More informationResearch Article Existence and Uniqueness of Homoclinic Solution for a Class of Nonlinear Second-Order Differential Equations
Applied Mathematics Volume 2012, Article ID 615303, 13 pages doi:10.1155/2012/615303 Research Article Existence and Uniqueness of Homoclinic Solution for a Class of Nonlinear Second-Order Differential
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationHopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd
Hopf bifurcations, and Some variations of diffusive logistic equation JUNPING SHIddd College of William and Mary Williamsburg, Virginia 23187 Mathematical Applications in Ecology and Evolution Workshop
More informationPermanence and global stability of a May cooperative system with strong and weak cooperative partners
Zhao et al. Advances in Difference Equations 08 08:7 https://doi.org/0.86/s366-08-68-5 R E S E A R C H Open Access ermanence and global stability of a May cooperative system with strong and weak cooperative
More informationLYAPUNOV-RAZUMIKHIN METHOD FOR DIFFERENTIAL EQUATIONS WITH PIECEWISE CONSTANT ARGUMENT. Marat Akhmet. Duygu Aruğaslan
DISCRETE AND CONTINUOUS doi:10.3934/dcds.2009.25.457 DYNAMICAL SYSTEMS Volume 25, Number 2, October 2009 pp. 457 466 LYAPUNOV-RAZUMIKHIN METHOD FOR DIFFERENTIAL EQUATIONS WITH PIECEWISE CONSTANT ARGUMENT
More informationOn Nonlinear Dynamics of Predator-Prey Models with Discrete Delay
Math. Model. Nat. Phenom. Vol. 4, No. 2, 2009, pp. 40-88 DOI: 0.05/mmnp/20094207 On Nonlinear Dynamics of Predator-Prey Models with Discrete Delay S. Ruan Department of Mathematics, University of Miami,
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp DOI: /m2an:
Mathematical Modelling and Numerical Analysis ESAIM: M2AN Modélisation Mathématique et Analyse Numérique M2AN, Vol. 37, N o 2, 2003, pp. 339 344 DOI: 10.1051/m2an:2003029 PERSISTENCE AND BIFURCATION ANALYSIS
More informationSTABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 20, Number 4, Winter 2012 STABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS XIHUI LIN AND HAO WANG ABSTRACT. We use an algebraic
More informationA New Fractional-Order Chaotic System and Its Synchronization with Circuit Simulation
Circuits Syst Signal Process (2012) 31:1599 1613 DOI 10.1007/s00034-012-9408-z A New Fractional-Order Chaotic System and Its Synchronization with Circuit Simulation Diyi Chen Chengfu Liu Cong Wu Yongjian
More informationBifurcations of Traveling Wave Solutions for a Generalized Camassa-Holm Equation
Computational and Applied Mathematics Journal 2017; 3(6): 52-59 http://www.aascit.org/journal/camj ISSN: 2381-1218 (Print); ISSN: 2381-1226 (Online) Bifurcations of Traveling Wave Solutions for a Generalized
More informationDynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting
Dynamics and Bifurcations in Predator-Prey Models with Refuge, Dispersal and Threshold Harvesting August 2012 Overview Last Model ẋ = αx(1 x) a(1 m)xy 1+c(1 m)x H(x) ẏ = dy + b(1 m)xy 1+c(1 m)x (1) where
More informationA Study of the Van der Pol Equation
A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple
More informationHopf bifurcation analysis of Chen circuit with direct time delay feedback
Chin. Phys. B Vol. 19, No. 3 21) 3511 Hopf bifurcation analysis of Chen circuit with direct time delay feedback Ren Hai-Peng 任海鹏 ), Li Wen-Chao 李文超 ), and Liu Ding 刘丁 ) School of Automation and Information
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: January 14, 2019, at 08 30 12 30 Johanneberg Kristian
More informationMidterm 2 Nov 25, Work alone. Maple is allowed but not on a problem or a part of a problem that says no computer.
Math 416 Name Midterm 2 Nov 25, 2008 Work alone. Maple is allowed but not on a problem or a part of a problem that says no computer. There are 5 problems, do 4 of them including problem 1. Each problem
More informationComplex Behavior in Coupled Nonlinear Waveguides. Roy Goodman, New Jersey Institute of Technology
Complex Behavior in Coupled Nonlinear Waveguides Roy Goodman, New Jersey Institute of Technology Nonlinear Schrödinger/Gross-Pitaevskii Equation i t = r + V (r) ± Two contexts for today: Propagation of
More informationConstruction of a New Fractional Chaotic System and Generalized Synchronization
Commun. Theor. Phys. (Beijing, China) 5 (2010) pp. 1105 1110 c Chinese Physical Society and IOP Publishing Ltd Vol. 5, No. 6, June 15, 2010 Construction of a New Fractional Chaotic System and Generalized
More informationAn 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 informationResearch Article Mean Square Stability of Impulsive Stochastic Differential Systems
International Differential Equations Volume 011, Article ID 613695, 13 pages doi:10.1155/011/613695 Research Article Mean Square Stability of Impulsive Stochastic Differential Systems Shujie Yang, Bao
More informationComputers and Mathematics with Applications. Adaptive anti-synchronization of chaotic systems with fully unknown parameters
Computers and Mathematics with Applications 59 (21) 3234 3244 Contents lists available at ScienceDirect Computers and Mathematics with Applications journal homepage: www.elsevier.com/locate/camwa Adaptive
More informationChaos suppression of uncertain gyros in a given finite time
Chin. Phys. B Vol. 1, No. 11 1 1155 Chaos suppression of uncertain gyros in a given finite time Mohammad Pourmahmood Aghababa a and Hasan Pourmahmood Aghababa bc a Electrical Engineering Department, Urmia
More informationA Model of Evolutionary Dynamics with Quasiperiodic Forcing
paper presented at Society for Experimental Mechanics (SEM) IMAC XXXIII Conference on Structural Dynamics February 2-5, 205, Orlando FL A Model of Evolutionary Dynamics with Quasiperiodic Forcing Elizabeth
More informationAndronov Hopf and Bautin bifurcation in a tritrophic food chain model with Holling functional response types IV and II
Electronic Journal of Qualitative Theory of Differential Equations 018 No 78 1 7; https://doiorg/10143/ejqtde018178 wwwmathu-szegedhu/ejqtde/ Andronov Hopf and Bautin bifurcation in a tritrophic food chain
More informationResearch Article New Exact Solutions for the 2 1 -Dimensional Broer-Kaup-Kupershmidt Equations
Hindawi Publishing Corporation Abstract and Applied Analysis Volume 00, Article ID 549, 9 pages doi:0.55/00/549 Research Article New Exact Solutions for the -Dimensional Broer-Kaup-Kupershmidt Equations
More informationGenerating a Complex Form of Chaotic Pan System and its Behavior
Appl. Math. Inf. Sci. 9, No. 5, 2553-2557 (2015) 2553 Applied Mathematics & Information Sciences An International Journal http://dx.doi.org/10.12785/amis/090540 Generating a Complex Form of Chaotic Pan
More informationUNIVERSIDADE DE SÃO PAULO
UNIVERSIDADE DE SÃO PAULO Instituto de Ciências Matemáticas e de Computação ISSN 010-577 GLOBAL DYNAMICAL ASPECTS OF A GENERALIZED SPROTT E DIFFERENTIAL SYSTEM REGILENE OLIVEIRA CLAUDIA VALLS N o 41 NOTAS
More informationCONVERGENCE BEHAVIOUR OF SOLUTIONS TO DELAY CELLULAR NEURAL NETWORKS WITH NON-PERIODIC COEFFICIENTS
Electronic Journal of Differential Equations, Vol. 2007(2007), No. 46, pp. 1 7. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) CONVERGENCE
More informationControlling chaos in Colpitts oscillator
Chaos, Solitons and Fractals 33 (2007) 582 587 www.elsevier.com/locate/chaos Controlling chaos in Colpitts oscillator Guo Hui Li a, *, Shi Ping Zhou b, Kui Yang b a Department of Communication Engineering,
More information