Double Hopf Bifurcation of a Hematopoietic System with State Feedback Control
|
|
- Elfrieda Townsend
- 5 years ago
- Views:
Transcription
1 Journal of Mathematics Research; Vol. 1, No. 4; August 218 ISSN E-ISSN Published by Canadian Center of Science and Education Double Hopf Bifurcation of a Hematopoietic System with State Feedback Control 1 College of Science, China Agricultural University, China Suqi Ma 1 Correspondence: Suqi Ma, College of Science, China Agricultural University, Beijing 183, China Received: March 13, 218 Accepted: March 28, 218 Online Published: June 28, 218 doi:1.5539/jmr.v1n4p116 Abstract URL: The dynamics of a system composed of hematopoietic stem cells and its relationship with neutrophils is ubiquitous due to periodic oscillating behavior induce cyclical neutropenia. Underlying the methodology of state feedback control with two time delays, double Hopf bifurcation occurs as varying time delay to reach its threshold value. By applying center manifold theory, the analytical analysis of system exposed the different dynamical feature in the classified regimes near double Hopf point. The novel dynamics as periodical solution and quasi-periodical attractor coexistence phenomena are explored and verified by numerical simulation. Keywords: Double Hopf bifurcation, cyclical neutropenia, time delay 1. Introduction Mathematically and Biologically, the discussion of the complex dynamics of hematopoietic cell model is ubiquitous and interesting due to its inherent highly nonlinear attribute. The hematopoietic stem cells(hscs originates from the bone marrow to proliferate and differentiate into all other types of blood cells. Underlying biomedical discpline, people introduced population dynamics to describe the relationship between cells compartments and the application of delay differential equations is the necessary useful tool on hand. The definition of a HSC has two properties: self-renewal and multipotent differentiation to cells with different function capabilities. People set forth all kinds of mathematcical models with the good understanding of how introduce compartment cells of proliferating phase into comparment cells of nonproliferating phase. Therefore, HSCs system is build on the basis that both nonproliferating phase and proliferating phase can produce, by successive divisions, all types of blood cells(white cells, red blood cells and platelets. Mackeys model proposed at the end of the seventies describe a decoupled quiescent phase of HSCS to study the dynamics of HSCs and its related diseases (Mackey, Since then it has been improved and analyzed by many authors, including Mackey and coauthors (Pujo Menjouet & Mackey, 24; Pujo Menjouet et al., 25; Fowler & Mackey, 22. In their models, delay is an salient factor in introducing HSCs population models. Proliferating cells can die with apoptosis and divide after a given time to produce two daughter cells, then immedieatelly enter a new nonproliferating phase. Due to the dysfunction in regulatory control process of blood cell production, some haematological disease emerge and cyclical neutropenia(cn is the most extensively studied with low level oscillating phenomena happening. Medical observation have exhibited that the neutrophils fall down into abnormally level in a period about 19 to 21 days, and longer period happen in some patients with 4 to 8 days in human being. With the exception of period oscillating from 11 to 15 days, the similar oscillation character of cyclical neutropenia have been observed in the contrast experiment result of grey collie. With the assumption that neutrophils experience maturation time to release into circulation through the body, people have put forth the cyclical neutropenia dynamical system which is based on the negative feedack regulation mechanism (Fowler & Mackey, 22; Bernard et al., 23; Santillan et al., 2; Ma et al.,21. Haurie suggests that oscillation mechanism in cyclical neutropenia is due to destabilization of HSCs regulatory would explain the fact that other cell lineages oscillate with the same period as cyclical neutropenia (Haurie et al., 1998; Haurie et al., 2; Hearn et al.,1998. Hopf bifurcation phenomenon has been observed in cyclical neutropenia model which lead to periodical oscillation of neutrophil numbers and the stable periodic solution branches out from equilibrium (Haurie et al., Biological meaningfully, Lei have developed CN model by introducing granulocyte-colony stimulating factors (G-CSF into CN administration with the assumption that G-CSF may decrease the amplification coefficient in neutrophil line for the treatment of cyclical neutropenia (Lei & Mackey, 214; Zhuge et al., 212. Recently, the dependence of neutrophil response on the period of simulated chemotherapy and the secondary response to G-CSF adminstaration is also reported physiologically and mathematically (Haurie et al.,
2 Lei developed cyclical neutropenia model including G-CSF administration which can be described as follows: dq = (β(q + k N (N + k δ Q + 2e r sτ s β(q τs Q τs, dt dn (1.1 = γ N N + A N k N (N τn Q τn dt Wherein, HSCs exist in its proliferating phase with nonlinear rate β, and enter into CN compartment with nonlinear rate k N. Commonly, both β and k N are given by the Holling function which listed f θ1 m k N (N = θ1 m + β(q = k cθ2 n Nm θ2 n + Qn f θ1 m k c θ2 n k N (N τn = θ1 m + N(t τ N m β(q τs = θ2 n + Q(t τ s n In model (1.1, HSCs can be self-renewal and differentiate to produce all kinds of blood cells including white blood cells, platelets and erythrocytes with rate k δ. γ s denotes apoptosis rate and τ s denotes proliferation time. The delay in circulating neutrophil is τ N which can be decomposed as different time period experienced, that is, the proliferation time τ NP (days of its progenitor cells and the maturation time τ NM (days of cyclical neutropenia. η NP expresses the proliferation rate and γ is the death rate in its maturation phase. Therefore, the amplification coefficient in circulating neutropenia is A N = e η NPτ NP γ τ NM In this paper, we introduce state feedback control mechanism to describe system (1.1 as dq = (β(q + k N (N + k δ Q + 2e r sτ s β(q τs Q τs + K(Q τn + Q τs 2Q, dt dn (1.2 = γ N N + A N k N (N τn Q τn dt Hopf bifurcation and the instability of system (1.2 may bring out the complex dynamical behavior as varying parameter continuously. However, it is tremendous rough work to work out the threshold of Hopf bifurcation as varying time delay. Recently, DDE-Biftool is developed to effectively solve the root attribution of the related characteristic equation and compute the stability of the equlibrium solution of system (1.2. Hopf bifurcation occurs as a pair of imaginary roots of the characteristic equation cross the imaginary axis from left half plane to right half plane transversally and gives rise to oscillating periodic solutions with small amplitudes. In system (1.1, Hopf bifurcation point is tracked by regarding amplification coefficient as bifurcation parameters and bistable regime of steady states and periodic solution is found by using DDE-biftool software (Engelborghs et al., 21; Green et al., 22. In this paper, we explore double Hopf bifurcation point of system (1.2 by tracking the secondary Hopf bifurcation lines which form the resonant region near death island (Zhang & Xu, 213; Zhang & Xu, 211; Ma et al., 29. Different features as asymptotically stability in death island and periodic or chaotic attractor after double Hopf bifurcation exhibit the richness of system dynamics inherently. In section 2, we discuss the stability of the related linearized system and double Hopf bifurcation point which is highlighted by death island is explored by using DDE-Biftool. In section 3, by using center manifold method and normal form theory, dynamics of the linearized system near the double Hopf bifurcation point is analyzed. In section 4, dynamics of nonlinear system near the double Hopf point is analyzed. Numerical simulation further verifies the classified dynamical results, and coexistence of periodical solutions and torus are detected near double Hopf bifurcation point. 2. Double Hopf Bifurcation Suppose E(Q, N is a positive equilibrium solution of (1.1, by using DDE-Biftool, stability of equilibrium solution E(Q, N is computed with accuracy and Hopf bifurcation line is tracked by varying parameter continuously. We choose γ N = 2.5, η NP = 2.542, γ =.25, f =.8, k δ =.12, θ 1 = 36, θ 2 =.16, r s =.631, τ NP = 5, ϵ =.583, τ s = fixed and set parameter k, τ N free, instability of the equilibium solution is detected after Hopf bifurcation happens. Therefore, Hopf bifurcation lines is drawn on (k, τ N plane as shown in Figure 1. It is seen that the secondary Hopf bifurcation happens at double Hopf bifurcation points which brings out the so-called death island. The linearized equation of system (1.2 is written as x = a 11 x + a 12 y + b 1 x(t τ s + K(x(t τ N + x(t τ s 2x y = γ N y + b 2 x(t τ N + b 3 y(t τ N (
3 where a 11 = β(q k N (N k δ β (Q Q, a 12 = k N (N Q, b 1 = 2e r sτ s (β(q + β (Q Q, b 2 = e τ NPη NP γ τ NM k(n, b 3 = e τ NPη NP γ τ NM k (N Q, The related characteristic equation of linearized system (2.1 is (λ, γ s, τ s = (a 11 2K λ + (b 1 + Ke λτ s + Ke λτ N a 12 b 2 e λτ N γ N λ + b 3 e λτ N = ( N 2 N 15 DH 15 DH k (a k (b Figure 1. Double Hopf bifurcation occurs in system(1.2. (ahopf bifurcation lines drawn on (k, τ N parameter plane (b The so-called death island is found near double Hopf bifurcation point and denoted by the shaded regime Roots of the characteristic equation (2.2 with zero real parts appears along Hopf bifurcation line, and one denotes λ = ±iω 1 which satisfies Eq(2.2 with ω 1 >. The transversallity of Hopf bifurcation can be verified which expresses the necessary condition of stability switching of the equilibrium solution E(Q, N, that is, equilibrium solution E(Q, N becomes unstable when entering into the regime enclosed and bounded by Hopf bifurcation line. The secondary Hopf bifurcation line is also detected with roots λ = ±iω 2 (ω 2 > and the double Hopf bifurcation point is explored at critical parameters k = , τ N = At the double Hopf point DH, both transversallity condition as dr(λ λ=iω1 < and dr(λ λ=iω2 > are satisfied. As shown in Figure 1(a, ahead of the double Hopf bifurcation point dτ N dτ N DH, equilibrium solution E(Q, N changes from local asymptotically stable to be unstable as increasing time delay to cross over Hopf bifurcation line, then retrieve its stability again. The stability switching of equilibrium solution E(Q, N brings out the death island which is shown by shaded regime in Figure 1(b which is the amplification of Figure 1(a. 3. Double Hopf Bifurcation As analyzed in section 2, Double Hopf bifurcation happened at pairs of critical values (k c, τ N. For convenience, set β = k c β 1 hereinafter. By setting x ϵx, y ϵy and k = k c + ϵk ϵ, τ N = τ N + ϵτ ϵ, system (1.2 is equivalent to x = a 11 x + a 12 y + b 1 x(t τ s + K(x(t τ N + x(t τ s + 2x + ϵl 11 x + ϵl 12 x(t τ s +K(x(t τ N x(t τ N + f (x, x(t τ s, y, y = γ N y + b 2 x(t τ N + b 3 y(t τ N + b 2 (x(t τ N x(t τ N +b 3 (y(t τ N y(t τ N + ϵl 21 x(t τ N +ϵl 22 y(t τ N + g(x(t τ N, y, y(t τ N, (3.1 where a 11 = k c β 1 (Q k N (N k δ k c β 1 (Q Q, a 12 = k N (N Q, b 1 = 2e r sτ s (k c β 1 (Q + k c β (Q Q, b 2 = e τ NPη NP γ τ NM k(n, b 3 = e τ NPη NP γ τ NM k (N Q, l 11 = k ϵ β 1 (Q k ϵ β 1 (Q Q, l 12 = 2e r sτ s (k ϵ β 1 (Q + k ϵ β (Q Q, l 21 = γ τ ϵ e τ NPη NP γ τ NM k N (N, l 22 = γ τ ϵ e τ NPη NP γ τ NM k N (N Q 118
4 and f (x, x(t τ s, y = 2 ϵk c e r sτ s β 1 (Q x 2 (t τ s + ϵk c e r sτ s β 1 (Q Q x 2 (t τ s +2ϵ ϵk ϵ e r sτ s β 1 (Q x 2 (t τ s 1 2 ϵkc β 1 (Q Q x 2 ϵ 1 6 k cβ 1 (Q Q x 3 ϵk c β 1 (Q x ϵk cβ 1 (Q x ϵ ϵk ϵ β 1 (Q Q x 2 ϵ ϵk ϵ β 1 (Q x 2 ϵk N (N xy 1 2 ϵk N (N xy ϵk N (N Q y ϵk N (N Q y 2 +ϵ ϵk ϵ e r sτ s β 1 (Q Q x 2 (t τ s + ϵk c e r sτ s β 1 (Q x 3 (t τ s +ϵ 1 3 k ce r sτ s β 1 (Q Q x 3 (t τ s, g(x(t τ N, y, y(t τ N = ϵe τ NPη NP γ τ NM k N (N x(t τ N y(t τ N ϵeτ NPη NP γ τ NM k N (N x(t τ N y 2 (t τ N ϵ ϵe τ NPη NP γ τ NM k N (N x(t τ N y(t τ N ϵe τ NP η NP γ τ NM k N (N Q y 2 (t τ N ϵeτ NPη NP γ τ NM k N (N Q y 3 (t τ N Choose phase space as C = C([ τ N, ], R 2, which is a Banach space composed as all continuous and differential function by mapping from [ τ N, ] R 2. For every ϕ C, one can define the linear operator as L(ϕ = (3.2 τ N [dη(θ]ϕ(θ, (3.3 where η : [ τ N, ] R 2 R 2 is a bounded variation function in [ τ N, ] and [ (a11 2Kδ(θ + (b dη(θ = 1 + Kδ(θ + τ s + Kδ(θ + τ N a 12 δ(θ b 2 δ(θ + τ N γ N δ(θ + b 3 δ(θ + τ N Furthermore, we define with and [ dη 2 (θ, τ ϵ = L(k ϵ, τ ϵ = dη 1 (θ, k ϵ = ] (3.4 τ N [dη 1 (θ, k ϵ + dη 2 (θ, τ ϵ ]ϕ(θ, (3.5 [ l11 δ(θ + l 12 δ(θ + τ s Kδ(θ + τ N Kδ(θ + τ N b 2 δ(θ + τ N (b 2 + ϵl 21 δ(θ + τ N b 3 δ(θ + τ N (b 3 + ϵl 22 δ(θ + τ N ] ] (3.6 (3.7 Suppose the linear operator at double Hopf bifurcation point (k c, τ N is dϕ H(ϕ =, θ [ τ N,, L(ϕ, θ =. (3.8 and { H(k ϵ, τ ϵ =, θ [ τ N,, L(k ϵ, τ ϵ ϕ, θ = (3.9 The nonlinear operator R is defined as ( and F(ϕ = f (ϕ 1 (, ϕ 1 (τ s, ϕ 2 ( g(ϕ 1 ( τ N, ϕ 2 (, ϕ 2 ( τ N R(ϕ = {, θ [ τn,, F(ϕ, θ =, System (3.1 can be rewritten as its operator form as (Hale, 23. (3.1 u (t = H(u t + H(k ϵ, τ ϵ u t + Ru t (3.11 with u = (x, y T, and u t = u(t + θ, for τ N θ <. For Ψ C = C([, τ N ], R 2 define the adjoint operator H ( of operator H( and adjoint operator H (k ϵ, τ ϵ of H(k ϵ, τ ϵ,respectively as H (ψ = dψ ds, < s τ N, τ N dη T (sψ( s, s =, (
5 and H (k ϵ, τ ϵ ψ = For ϕ C, ψ C, define the belinear form (, : C C R as dψ ds, < s τ N, τ N [dη T 1 (s, k ϵ + dη T 2 (s, τ ϵ]ψ( s, s =, (3.13 < ψ, ϕ >= ψ T (ϕ( + ψ T (ξ + τ s R 1 ϕ(ξdξ + ψ T (ξ + τ N R 2 ϕ(ξdξ (3.14 τ s τ N By the above discussion, there are two pairs of imaginary roots at the double Hopf bifurcation point (k ϵ, τ ϵ, and other characteristic roots with negative real parts. We denote the collection set Λ = {iω 1, iω 1, iω 2, iω 2 } and decompose the phase space C to be the direct sum of its subspaces as C = P Λ Q Λ, with P Λ being the characteristic space of Λ and Q Λ being the corresponding complement subspace. System (3.1 is an ordinary differential system on the Banach space C of functions from [ τ N, ] to R 2, which is bounded and continuous on ( τ N, with a possible jump discontinuity at (Buono et al., 23. We choose Φ(θ = (q 1 (θ, q 1 (θ, q 2 (θ, q 2 (θ, Ψ(θ = (q 1 (θ, q 1 (θ, q 2 (θ, q 2 (θ, (3.15 as the bases of P Λ and P Λ, where ( γn + iω q1(θ = 1 b 3 e iω 1τ N b 2 e iω 1τ N e iω1θ, ( γn + iω q 2 (θ = 2 b 3 e iω 2τ N b 2 e iω 2τ N e iω2θ, for τ N θ <, are respectively the eigenvectors of operator H( with respect to the characteristic root iω 1, iω 2, and the eigenvectors corresponding to adjoint operator A ( with respect to iω 1 and iω 2 are respectively as ( q 1 (s = l b 2 e iω 1τ N 1 a 11 + iω 1 + (b 1 + Ke iω 1τ s + Ke iω 1τ N e iω1s, ( q 2 (s = l b 2 e iω 2τ N ( a 11 + iω 2 + (b 1 + Ke iω 2τ s + Ke iω 2τ N e iω 2s for s τ N,where with It can be verified that l 1 = is satisfied and further we can prove that e iω τ N s 1 τ N + s 2 τ s + s 3, l 2 = 1 v 1 τ N + v 2 τ s + v 3 s 1 = Kω ( (2iKγ N + (2iKb 3 e iω 1τ N ia 12 b 2 ω 1 + Kγ 2 N +a 12 b 2 Kb 3 e iω 1τ N γ N + Kb 2 3 e2iω 1τ N, s 2 = (ib 2 Ke iω 1τ s + ib 2 b 1 e iω 1τ s ω 1 + ( b 2 Ke iω 1τ s b 2 b 1 e iω 1τ s γ N +b 2 Kb 3 e iω 1(τ N +τ s + b 2 b 1 b 3 e iω 1(τ N +τ s, s 3 = e iω 1τ N ω (( 2ie iω 1τ N γ N + 2ib 3 ω 1 +e iω 1τ N γ 2 N 2γ Nb 3 + b 2 3 eiω 1τ N + a 12 b 2, v 1 = Kω ( (2iKγ N + (2iKb 3 e iω 2τ N ia 12 b 2 ω 2 + Kγ 2 N +a 12 b 2 Kb 3 e iω 2τ N γ N + Kb 2 3 e2iω 2τ N, v 2 = (ib 2 Ke iω 2τ s + ib 2 b 1 e iω 2τ s ω 2 + ( b 2 Ke iω 2τ s b 2 b 1 e iω 2τ s γ N +b 2 Kb 3 e iω 2(τ N +τ s + b 2 b 1 b 3 e iω 2(τ N +τ s, v 3 = e iω 2τ N ω (( 2ie iω 2τ N γ N + 2ib 3 ω 2 +e iω 2τ N γ 2 N 2γ Nb 3 + b 2 3 eiω 2τ N + a 12 b 2, < q 1, q 1 >=< q 2, q 2 >= 1 < q 1, q 1 >=< q 1, q 2 >=< q 1, q 2 >=, < q 2, q 2 >=< q 2, q 1 >=< q 2, q 1 >=. (3.16 It is noticed that Aq 1 (θ = iω 1 q 1 (θ, Aq 2 (θ = iω 2 q 2 (θ, τ N θ, A q 1 (s = iω 1q 1 (s, A q 2 (s = iω 2q 2 (s, s τ N. 12
6 Let z = (z 1, z 1, z 2, z 2 T, for every u t C, and v t Q Λ C ([ τ N, ], R 2, we have the expression u t = Φz + v t and the linearized system corresponds to (3.11 is u t = H(u t + H(k ϵ, τ ϵ u t, (3.18 with the reduction form on its center manifold as z = Bz+ < Ψ, H(k ϵ, τ ϵ u t > = Bz + Ψ T (A(k ϵ, τ ϵ (Φz + v t, (3.19 where B = iω 1 iω 1 iω 2 iω 2 The matrix B generates the torus group T 2 = S 1 S 1 whose action on C 2 is given by (θ 1, θ 2 (z 1, z 2 = (e iθ 1 z 1, e iθ 2 z 2 Then the T 2 -equivalent normal form, which is truncated to the quadratic order, becomes ( ( ( ( ż1 iω1 z1 ϵkϵ c = z 1 + ϵτ ϵ c 1 z 1 ż 2 iω 2 z 2 ϵk ϵ c 2 1 z 1 + ϵτ ϵ c 2 2 z 1 (3.19 with c 1 1 = (( (2iω 1 2γ N exp(( r s iω 1 τ s iω 1 τ N + 2exp(( r s iω 1 τ s (2iω 1 τ N b 3 +(γ N + iω 1 exp( iω 1 τ N exp( (2iω 1 τ N b 3 l 1 ϵb 2 (β 1 (Q Q + β 1 (Q, c 1 2 = ϵl 1(γ A N b 3 (Q b 2 k N (N b 3 k N (N (b 1 + Kexp( iω 1 (2τ N + τ s +(k N A N γ (γ N + iω 1 b 3 + ω 1 (iγ N ω 1 b 2 (b 1 + Kexp( iω 1 (τ N + τ s +( k N (N A N ( ω (iγ N ia 11 ω 1 a 11 γ N γ b 3 +ω 1 (iω ω 1(γ N a 11 + ia 11 γ N b 2 exp( iω 1 τ N + (A N (k N (N (iω 1 a 11 b 3 +( ik N (N Q b 2 + ikk N (N ω 1 + k N (N Q b 2 a 11 + Kγ N k N (N γ exp( (2iω 1 τ N +Kexp( (3iω 1 τ N ( γ A N b 3 k N (N + (k N (N γ Q A N + iω 1 b 2 b 3, c 2 1 = (( (2iω 1 2γ N exp(( rs iω 1 τ s iω 1 τ N + 2exp(( r s iω 1 τ s (2iω 1 τ N b 3 +(γ N + iω 1 exp( iω 1 τ N exp( (2iω 1 τ N b 3 l 1 ϵb 2 (β 1 (Q Q + β 1, c 2 2 = l 2ϵ(k N (N A N b 3 (b 1 + Kγ (γ N + iω 1 exp( iω 1 τ N iω 2 τ s +γ A N b 3 (Q b 2 k N (N b 3 k N (N (b 1 + Kexp( (2iω 1 τ N iω 2 τ s +Kγ A N b 3 (Q b 2 k N (N b 3 k N (N exp( iτ N (2ω 1 + ω 2 +ω 2 (iγ N ω 2 (b 1 + Kb 2 exp( iω 2 (τ N + τ s +k N (N A N Kb 3 γ (γ N + iω 1 exp( iτ N (ω 1 + ω 2 A N b 3 (iω 2 a 11 (Q b 2 k N (N b 3 k N (N γ exp( (2iω 1 τ N k N (N A N ((iγ N ω 1 ω 2 a 11 γ N + iω 1 b 3 γ exp( iω 1 τ N +ω 2 ((iω ω 2(γ N a 11 + ia 11 γ N exp( iω 2 τ N + ikexp( (3iω 2 τ N b 3 b 2, In polar coordinate (ρ 1, θ 1, ρ 2, θ 2, we express the associated amplitude equation generated from system (3.19 as ρ 1 = ϵr(k ϵc τ ϵc 1 2 ρ 1, ρ 2 = ϵr(k ϵc τ ϵc 2 2 ρ 2 (
7 4. Center Manifold Reduction Furthermore, consider the nonlinear system u t = H(u t + Ru t (4.1 It is directly computed that the reduction form derived from Eq.(4.1 is z 1 = iω 1z + q T 1 (F(Φz + v t = iω 1 z + q T 1 ( f ˆ(z 1, z 1, z 2, z 2 z 2 = iω 2z + q T 2 (F(Φz + v t = iω 2 z + q T 2 ( f ˆ(z 1, z 1, z 2, z 2 v t = Av t + (I ΠX (F(Φz + v t (4.2 with I identity and Π is the projection mapping from C([ τ, ], R 2 to the subspace P Λ, and X denotes {, θ [ τn,, X = I, θ = where ˆ f (z 1, z 1, z 2, z 2, v t = F(2R(z 1 q 1 (θ + 2R(z 2 q 2 (θ + v t (θ (4.3 We write v t into its Taylor series v t = w(z 1, z 1, z 2, z 2 = w i jkl z i 1 z j 1 zk 2 zl 2 (4.4 i+ j+k+l 2 with coefficient vector w i jkl only depend on θ. We now expand ˆ f into the Taylor series and g 1 (z 1, z 1, z 2, z 2 = q T 1 ( f ˆ(z 1, z 1, z 2, z 2, v t = 3 i+ j+k+l 2 g1 i jkl zi 1 z j 1 zk 2 zl 2, g 2 (z 1, z 1, z 2, z 2 = q T 2 ( f ˆ(z 1, z 1, z 2, z 2, v t = 3 i+ j+k+l 2 g2 i jkl zi 1 z j 1 zk 2 zl 2 (4.5 From system (4.1, we obtain v t = Av t + 2R{p T 1 ( ˆ f (z 1, z 1, z 2, z 2 q 1 (θ} 2R{p T 2 ( f ˆ(z 1, z 1, z 2, z 2 q 2 (θ}, τ N θ <, 2R{p T 2 ( f ˆ(z 1, z 1, z 2, z 2 q 2 (θ} 2R{p T 1 ( f ˆ(z 1, z 1, z 2, z 2 q 1 (θ} + f ˆ(z 1, z 1, z 2, z 2, θ = (4.6 The first equation in Eq.(4.6 can be rewritten as w (z 1, z 1, z 2, z 2 = Aw(z 1, z 1, z 2, z 2 + B(z 1, z 1, z 2, z 2, τ N θ < (4.7 where B(z 1, z 1, z 2, z 2 = 2R{p T 1 ( f ˆ(z 1, z 1, z 2, z 2 q 1 (θ} 2R{p T 2 ( f ˆ(z 1, z 1, z 2, z 2 q 2 (θ} = i+ j+k+l 2 B i jkl z i 1 z j 1 zk 2 zl 2, with the coefficients B i jkl = (B 1 i jkl (θ, B2 i jkl (θt. It can be derived from Eq.(4.6 that (A 2iω 1 w 2 = B 2 (θ, Aw 11 = B 11 (θ, (A + 2iω 1 w 2 = B 2 (θ, (A iω 1 iω 2 w 11 = B 11 (θ, (A iω 1 + iω 2 w 11 = B 11 (θ, (A + iω 1 iω 2 w 11 = B 11 (θ, (A + iω 1 + iω 2 w 11 = B 11 (θ, (A 2iω 2 w 2 = B 2 (θ, Aw 11 = B 11 (θ, (A + 2iω 2 w 2 = B 2 (θ, (
8 Taking account of formula of operator A, we have dw 2 = 2iω 1 w 2 (θ B 2 (θ, dw 11 = B 11 (θ, dw 2 = 2iω 1 w 2 (θ B 2 (θ, dw 11 = i(ω 1 + ω 2 w 11 (θ B 11 (θ, dw 11 = i(ω 1 ω 2 w 11 (θ B 11 (θ, dw 11 = i( ω 1 + ω 2 w 11 (θ B 11 (θ, dw 11 = i(ω 1 + ω 2 w 11 (θ B 11 (θ, dw 2 = 2iω 2 w 2 (θ B 2 (θ, dw 11 = B 11 (θ, dw 2 = 2iω 2 w 2 (θ B 2 (θ, Integrating Eq.(4.9 under the initial condition given be the second equation in Eq.(4.6, we can compute the coefficient vectors w i jkl (θ with i + j + k + l = 2. Using the standard technique as invertible parameter-dependent complex coordinates transformation which is introduced in (Buono et al., 23, the normal form at (k c, τ N of Eq.(4.1 is described as (4.9 z 1 = iω 1z 1 + B 11 z 2 1 z 1 + B 12 z 1 z 2 z 2, z 2 = iω 2z 2 + B 21 z 1 z 1 z 2 + B 22 z 2 2 z 2 (4.1 In polar coordination (ρ 1, θ 1, ρ 2, θ 2, the amplitude equation generated from Eq.(4.1 is { ρ 1 = ρ 1 (R(B 11 ρ R(B 12ρ 2 2, ρ 2 = ρ 2(R(B 21 ρ R(B 22ρ 2 2, (4.11 Combining the results of Eqs(3.19 and Eqs(4.1, the formal normal form arising from system (4.11 is z 1 = iω 1z 1 + ϵc 1 1 k ϵz 1 + ϵc 2 1 τ ϵz 1 + B 11 z 2 1 z 1 + B 12 z 1 z 2 z 2, z 2 = iω 2z 2 + ϵc 1 2 k ϵz 2 + ϵc 2 2 τ ϵz 2 + B 21 z 1 z 1 z 2 + B 22 z 2 2 z 2, (4.12 Using polar coordinate transformation z 1 = ρ 1 e iθ 1, z 2 = ρ 2 e iθ 2, system (4.12 is transformed into ρ 1 = ρ 1(µ 1 + R(B 11 ρ R(B 12ρ 2, ρ 2 = ρ 2(µ 2 + R(B 21 ρ R(B 22ρ 2 2, θ 1 = ω 1 + v 1 + I(B 11 ρ I(B 12ρ 2, (4.13 θ 2 = ω 2 + v 2 + I(B 21 ρ I(B 22ρ 2 2, with µ i = ϵ(r(c 1 i k ϵ + R(c 2 i τ ϵ, v i = ϵ(i(c 1 i k ϵ + I(c 2 i τ ϵ for i = 1, 2. Using Maple software, we compute the coefficients g i jlk with i + j + k + l = 2, 3,and further obtain the following formulas (Zhang, S. & Xu, J., 211. B 11 = g i ω 1 g 1 11 g1 2 + i ω 2 (g 1 11 g2 11 g1 11ḡ2 11 i 2ω 1 +ω 2 g 1 11ḡ2 2 i 2ω 1 ω 2 g 1 11 g2 2 i ω 1 g i 3ω 1 g 1 2 2, B 12 = g i ω 2 (g 1 11 g2 11 g1 11ḡ i ω 1 (2g 1 2 g1 11 g1 11ḡ1 11 g1 11ḡ2 11 g 1 11 g2 11 2i ω 1 +2ω 2 g 1 2ḡ2 11 2i ω 1 2ω 2 g 1 2 g2 11 i 2ω 1 ω 2 g i 2ω 1 +ω 2 g , B 21 = g i ω 1 (g 1 11 g2 11 ḡ1 11 g i ω 2 (2g 2 2 g2 11 g2 11ḡ2 11 g1 11 g2 11 ḡ 2 11 g i 2ω 1 ω 2 g 1 11 g2 2 2i 2ω 1 +ω 2 (ḡ 1 11 g2 2 i 2ω 2 ω 1 g i ω 1 +2ω 2 g , B 22 = g i ω 1 (g 1 11 g2 11 ḡ1 11 g i ω 2 g 2 11 g2 2 i 2ω 2 ω 1 g 1 2 g2 11 i 2ω 2 +ω 1 ḡ 1 2 g2 11 i ω 2 g i 3ω 2 g At double Hopf point DH with Q = , N = , ω 1 = ,ω 2 = , we compute the coefficients g 1,2 i jlk, (i + j + k + l = 2, 3. Therefore, by formula(4.14,choose ϵ =.1, (
9 then we have We also have B 11 = e e 4i e 2k ϵ +( e 1ik ϵ + ( e 2ik 2 ϵ e 1k 2 ϵ, B 12 = e i e 1k ϵ +( e 1ik ϵ e 2k 2 ϵ +( e 2ik 2 ϵ, B 21 = e i e 1k ϵ ( e 1ik ϵ ( e 2ik 2 ϵ e 2k 2 ϵ, B 22 = e e 1i e 1k ϵ ( e 2ik ϵ ( e 3ik 2 ϵ e 3k 2 ϵ µ 1 = e 2k ϵ + ( e 2ik ϵ ( e 2iτ ϵ e 2τ ϵ ( e 4iτ 2 ϵ e 5τ 2 ϵ, µ 2 = e 3k ϵ + ( e 3ik ϵ ( e 2iτ ϵ e 2τ ϵ ( e 4iτ 2 ϵ e 4τ 2 ϵ By time continuous transformation with direction preserved, system (4.13 is equivalent to (4.15 (4.16 with a = R(µ 1 B 11, b = R(µ 2. B 11 ρ 1 = ρ 1(a + R(B 11 B 11 ρ 2 = ρ 2(b + R(B 21 B 11 ρ R(B 12 ρ 2 2 B 11, ρ R(B ( ρ 2 2 B 11. System (4.17 may have equilibrium solutions (ρ 1,, (, ρ 2 and (ρ 1, ρ 2 and stability of equilibrium solution can be derived via computation of its corresponding Jacobin determinant. Therefore, based on system(4.16, we deduce the following results (INeutral saddle bifurcation of equilibrium (ρ 1, happens if µ 2 + (2B 11 + B 21 ρ 2 1 = ; (IINeutral saddle bifurcation of equilibrium solution (, ρ 2 happens if µ 1 + (B B 22 ρ 2 2 = ; (III The equlibrium solution (ρ 1, ρ 2 loses its stability through line HL3 : B 11ρ B 22ρ 2 2 = and correspondingly, the associated torus of system (4.17 getting into torus instability state. As shown in Figure2, system (4.17 has the equilibrium solution (, ρ 2 with parameter pairs(k c, τ N chosen in the regime above boundary line L1 in green color, and the regime bounded by line in red color and line in black color exhibits torus solution. We discuss system dynamics in detail with the partition of plane into different regimes. The stable periodic solution observed in regime I, V I is shown in Figure3(a-(f. The periodic solution exhibited in regime I further bifurcates into quasi-periodical solution after crossing neutral saddle line HL1 shown in cyan color. The scenarios of the bifurcating quasi-periodical solutions observed in regimes from (II to V is plotted as shown in Figure4. Due to neutral line HL2 in purple color, system may have multi attractors coexistence in regime IV. We also obtain the observed coexisted attractors in regime II and V alike, as shown in Figure5. 5. Discussion A blood cell model which composed of two compartment cells known as hematopoietic stem cells and its relationship with neutrophil was investigated. The dysfunction in regulatory control process of blood cell production induce the observed hematological disease and cyclical neutropenia model including G-CSF adminstration was developed. As system state dependent delay reflects different time stage during feedback regulation process, we introduced state feedback control with time delay originally arising in system which is subtle disturbance around stationary equilibrium solution. We use DDE-Biftool software to explore double Hopf bifurcation of system via tracking second Hopf bifurcaton lines as varying parameters continuously and dynamics near double Hopf point was discussed. A classified method of dynamics near double Hopf point was developed by using center manifold theory. The analytical norm form near double Hopf point is computed based on fundamental theory of functional differential equation and the observed numerical simulation results were in coincidence with dynamics of the classified method which further verified the correctness of theoretical analytical results. 124
10 15.3 N VI V IV II III 14.6 I k c Figure 2. The dynamical classify results near the double Hopf point DH in parameter space (k, τ N plane Q(t N(t- s t Q(t (a t (c Q(t 7 6 N(t- 5 s (b Q(t (d Figure 3. Periodical solution of system (1.2.(a Phase portraits of solution with parameter k c = 1.1, τ N = 14.6; (b Time series solution with parameter k c = 1.1, τ N = 14.6; (c Phase portraits of solution with parameter k c =.9, τ N = 15.1; (d Time series solution with parameter k c =.9, τ N =
11 N(t- s 8 1 N(t- s Q(t N(t- 1 s (a Q(t (c N(t- s Q(t (b Q(t (d Figure 4. The numerical simulation results of system (1.2 near the double Hopf point DH. (atorus at parameter pairs (k c, τ N lying in regime II; (b Torus at parameter pairs (k c, τ N lying in regime III;(cTorus at parameter pairs (k c, τ N lying in regime IV;(dTorus at parameter pairs (k c, τ N lying in regime V 15 2 N(t- s 1 15 N(t- s Q(t (a Q(t (b Figure 5. The coexistence phenomena of torus solutions. (achosen parameter pairs (k c, τ N lying in regime III; (b Chosen parameter pairs (k c, τ N lying regime V 126
12 References Bernard, S., Bélair, J., & Mackey, M. C. (23. Oscillation in cyclical neutropenia: new evidence based on mathematical modelling. Journal of Theory Biology, 223, Buono, P. L., & Bélair, J. (23. Restrictions and unfolding of double Hopf bifurcation in functional differentail equations. J. Differentail Equations, 189, Engelborghs, K., Luzyanina, T., & Samacy, G. (21. DDE-BIFTOOL v2.: a Matlab package for bifurcation analysis of delay differential equations. Technical Report TW-33, Department of Compute Science, K.U. Leuven, Belgium. koen/delay/ddfbiftool.shtml. Fowler, A. C., & Mackey, M. C. (22. Relaxation oscillations in a class of delay differential equations. SIAM J. Appl. Math, 63(1, Green, K., Krauskopf, B., & Engelborghs, K. (22. Bistability and torus break-up in a semiconductor lase with phaseconjugate feedback. Phys.D, 173, Hale, J. (23. Theory of Functional Differential Equations, World Publishing Corporation, Beijing, China. Fowler, A. C. (1997. Mathematical Models in the Applied Sciences. Cambridge Texts in Applied Mathematics, Cambridge University Press, Cambridge. Haurie, C., Dale, D. C., & Mackey, M. C. (1998. Cyclical neutropenia and other periodic hematological diseases: a review of mechanisms and mathematical models. Blood, 92, Haurie, C., Dale, D. C., Rudnicki, R., & Mackey, M. C. (2. Modelling complex neutrophil dynamics in the grey collie. Journal of Thoery Biology 24, Hearn, T., Haurie, C., & Mackey, M. C.(1998. Cyclical neutropenia and the peripherial control of white blood cell production. Journal of Theory Biology, 192, Lei, J. Z., & Mackey, M. C. (214. Understanding and treating cytopenia through mathematical modelling, 844, , Springer New York. Mackey, M. C. (1978. Unified hypothesis for the origin of aplastic anemia and periodic hematopoiesis. Blood 51(5, 941C956. Ma, S. Q., Wang, X. H., Lei, J. H., & Feng, Z. S. (21. Dynamics of the delay hematological cell model. International Journal of Biomathematics, 3, Ma, S. Q., Feng, Z. S., & Lu, Q. S. (29. Dynamics and double Hopf bifurtions of the Rose-Hindmarsh model with time delay. I nternational Journal of Bifurcation and Chaos, 19, Pujo-Menjouet, L., & Mackey, M. C. (24. Contribution to the study of periodic chromic myelogenous leukemia. Comptes Rendus Biol., 327, Pujo-Menjouet, L., Bernard, S., & Mackey, M. C. (25. Long period oscillation in a G model of hematopoietic stem cells. SIAM J. Appl. Dyn.Syst., 4(2, Santillan, M., Mahaffy, J. M., Bélair, J., & Mackey, M. C. (2. Regulation of platelet production: The normal response to perturbation and cyclical platelet disease. Journal of Theory Biology, 26, Zhang, S., & Xu, J. (213. Quasiperiodic motion induced by heterogeneous delays in a simplified internet congestion control model. Nonlineat Analysis: Real World Applications, 14, Zhang, S., & Xu, J. (211. Time-varying delayed feedback control for an internet congestion control model. Discrete and Continuous Dynamical Systems Series B., 211, Zhuge, C. J., Lei, J. Z., & Mackey, M. C. (212. Neutrophil dynamics in response to chemotherapy and G-SCF. Journal of Theory Biology, 293(1, Copyrights Copyright for this article is retained by the author(s, with first publication rights granted to the journal. This is an open-access article distributed under the terms and conditions of the Creative Commons Attribution license ( 127
Stability and Hopf bifurcation in a mathematical model of pluripotent stem cell dynamics
Nonlinear Analysis: Real World Applications 6 (2005) 651 670 www.elsevier.com/locate/na Stability and Hopf bifurcation in a mathematical model of pluripotent stem cell dynamics Mostafa Adimy a, Fabien
More informationDelay Differential Equations and Autonomous Oscillations in Hematopoietic Stem Cell Dynamics Modeling
Math. Model. Nat. Phenom. Vol. 7, No. 6, 212, pp. 1 22 DOI: 1.151/mmnp/212761 Delay Differential Equations and Autonomous Oscillations in Hematopoietic Stem Cell Dynamics Modeling M. Adimy 1,2, F. Crauste
More informationPeriodic Solutions in a Mathematical Model for the Treatment of Chronic Myelogenous Leukemia
Math. Model. Nat. Phenom. Vol. 7, No.,, pp. 35-44 DOI:.5/mmnp/7 Periodic Solutions in a Mathematical Model for the Treatment of Chronic Myelogenous Leukemia A. Halanay Department of Mathematics I, Politehnica
More informationChapter 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 informationAdding Self-Renewal in Committed Erythroid Progenitors Improves the Biological Relevance of a Mathematical Model of Erythropoiesis
Adding Self-Renewal in Committed Erythroid Progenitors Improves the Biological Relevance of a Mathematical Model of Erythropoiesis Fabien Crauste a, Laurent Pujo-Menjouet a, Stéphane Génieys a, Clément
More informationarxiv: v1 [math.ap] 16 Apr 2009
A Mathematical Study of the Hematopoiesis Process with Applications to Chronic Myelogenous Leukemia arxiv:94.2493v1 [math.ap] 16 Apr 29 Mostafa Adimy, Fabien Crauste Year 24 and Shigui Ruan Laboratoire
More informationKrauskopf, B., Erzgraber, H., & Lenstra, D. (2006). Dynamics of semiconductor lasers with filtered optical feedback.
Krauskopf, B, Erzgraber, H, & Lenstra, D (26) Dynamics of semiconductor lasers with filtered optical feedback Early version, also known as pre-print Link to publication record in Explore Bristol Research
More informationThree 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 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 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 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 informationHOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS
PERIODICA POLYTECHNICA SER. MECH. ENG. VOL. 48, NO. 2, PP. 89 200 2004 HOPF BIFURCATION CALCULATIONS IN DELAYED SYSTEMS Gábor OROSZ Bristol Centre for Applied Nonlinear Mathematics Department of Engineering
More informationUniversity of Bristol - Explore Bristol Research. Early version, also known as pre-print
Erzgraber, H, Krauskopf, B, & Lenstra, D (2004) Compound laser modes of mutually delay-coupled lasers : bifurcation analysis of the locking region Early version, also known as pre-print Link to publication
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 informationDifference 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 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 informationModels 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 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 informationAdding self-renewal in committed erythroid progenitors improves the biological relevance of a mathematical model of erythropoiesis
Journal of Theoretical Biology 25 (28) 322 338 www.elsevier.com/locate/yjtbi Adding self-renewal in committed erythroid progenitors improves the biological relevance of a mathematical model of erythropoiesis
More informationPREDATOR-PREY MODELS WITH DISTRIBUTED TIME DELAY
PREDATOR-PREY MODELS WITH DISTRIBUTED TIME DELAY PREDATOR-PREY MODELS WITH TIME DISTRIBUTED DELAY By ALEXANDRA TESLYA, M.SC. A Thesis Submitted to the School oif Graduate Studies in Partial Fulfillment
More informationNBA 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 informationStability 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 informationAsynchronous oscillations due to antigenic variation in Malaria Pf
Asynchronous oscillations due to antigenic variation in Malaria Pf Jonathan L. Mitchell and Thomas W. Carr Department of Mathematics Southern Methodist University Dallas, TX SIAM LS, Pittsburgh, 21 Outline
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 informationNonlinear Stability and Bifurcation of Multi-D.O.F. Chatter System in Grinding Process
Key Engineering Materials Vols. -5 (6) pp. -5 online at http://www.scientific.net (6) Trans Tech Publications Switzerland Online available since 6//5 Nonlinear Stability and Bifurcation of Multi-D.O.F.
More informationAnalytical 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 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 informationDynamic Feedback and the Design of Closed-loop Drug Delivery Systems
Dynamic Feedback and the Design of Closed-loop Drug Delivery Systems John Milton 1,2, Sue Ann Campbell 2,3,4 and Jacques Bélair 2,4,5 1) Department of Neurology, The University of Chicago Hospitals, MC
More informationBifurcations in Delayed Lotka-Volterra Intraguild Predation Model
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 5-696 Vol. 37 Nos. - (4) pp. - Bifurcations in Delayed Lotka-Volterra Intraguild Predation Model Juancho A. Collera Department
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 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 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 informationDELAY-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 informationStability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations
Stability and Hopf bifurcation analysis of the Mackey-Glass and Lasota equations Sreelakshmi Manjunath Department of Electrical Engineering Indian Institute of Technology Madras (IITM), India JTG Summer
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 informationBIFURCATION TO TRAVELING WAVES IN THE CUBIC-QUINTIC COMPLEX GINZBURG LANDAU EQUATION
BIFURCATION TO TRAVELING WAVES IN THE CUBIC-QUINTIC COMPLEX GINZBURG LANDAU EQUATION JUNGHO PARK AND PHILIP STRZELECKI Abstract. We consider the 1-dimensional complex Ginzburg Landau equation(cgle) which
More informationMathematical Models for Erythropoiesis
Mathematical Models for Erythropoiesis Joseph M. Mahaffy Department of Mathematical and Statistical Sciences San Diego State University January 2003 p.1/37 Collaborators Jacques Bélair - Université de
More information8. MATHEMATICAL MODELS OF HEMATOPOIETIC CELL REPLICATION AND CONTROL
8. MATHEMATICAL MODELS OF HEMATOPOIETIC CELL REPLICATION AND CONTROL Michael C. Mackey 8.1 INTRODUCTION The biological sciences offer an abundance of examples in which the dynamics are dependent on the
More informationPhase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
Commun. Theor. Phys. (Beijing, China) 47 (2007) pp. 265 269 c International Academic Publishers Vol. 47, No. 2, February 15, 2007 Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect
More informationAnalysis of a lumped model of neocortex to study epileptiform ac
of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized
More informationas Hopf Bifurcations in Time-Delay Systems with Band-limited Feedback
as Hopf Bifurcations in Time-Delay Systems with Band-limited Feedback Lucas Illing and Daniel J. Gauthier Department of Physics Center for Nonlinear and Complex Systems Duke University, North Carolina
More informationDynamics of Erythroid Progenitors and Erythroleukemia
Math. Model. Nat. Phenom. Vol. 4, No. 3, 2009, pp. 210-232 DOI: 10.1051/mmnp/20094309 Dynamics of Erythroid Progenitors and Erythroleukemia N. Bessonov a,b, F. Crauste b1, I. Demin b, V. Volpert b a Institute
More informationNON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 016, 91 95 NON-STATIONARY RESONANCE DYNAMICS OF THE HARMONICALLY FORCED PENDULUM Leonid I. Manevitch Polymer and Composite Materials Department N. N. Semenov Institute
More informationUniversity of Bristol - Explore Bristol Research. Early version, also known as pre-print
Krauskopf, B., Lenstra, D., Fischer, A. P. A., Erzgraber, H., & Vemuri, G. (6). Feedback phase sensitivity of a semiconductor laser subject to filtered optical feedback. Early version, also known as pre-print
More informationPhase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
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 informationarxiv: v1 [q-bio.cb] 21 Dec 2015
arxiv:1512.6896v1 [q-bio.cb] 21 Dec 215 A mathematical model of granulopoiesis incorporating the negative feedback dynamics and kinetics of G-CSF/neutrophil binding and internalisation M. Craig A.R. Humphries
More informationNonchaotic random behaviour in the second order autonomous system
Vol 16 No 8, August 2007 c 2007 Chin. Phys. Soc. 1009-1963/2007/1608)/2285-06 Chinese Physics and IOP Publishing Ltd Nonchaotic random behaviour in the second order autonomous system Xu Yun ) a), Zhang
More informationNumerical bifurcation analysis of delay differential equations
Numerical bifurcation analysis of delay differential equations Dirk Roose Dept. of Computer Science K.U.Leuven Dirk.Roose@cs.kuleuven.be Acknowledgements Many thanks to Koen Engelborghs Tatyana Luzyanina
More informationStability and Hopf Bifurcation Analysis of the Delay Logistic Equation
1 Stability and Hopf Bifurcation Analysis of the Delay Logistic Equation Milind M Rao # and Preetish K L * # Department of Electrical Engineering, IIT Madras, Chennai-600036, India * Department of Mechanical
More informationStability and hybrid synchronization of a time-delay financial hyperchaotic system
ISSN 76-7659 England UK Journal of Information and Computing Science Vol. No. 5 pp. 89-98 Stability and hybrid synchronization of a time-delay financial hyperchaotic system Lingling Zhang Guoliang Cai
More informationarxiv: v2 [physics.optics] 15 Dec 2016
Bifurcation analysis of the Yamada model for a pulsing semiconductor laser with saturable absorber and delayed optical feedback Abstract Soizic Terrien 1, Bernd Krauskopf 1, Neil G.R. Broderick 2 arxiv:1610.06794v2
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 informationModelling of cell choice between differentiation and apoptosis on the basis of intracellular and extracellular regulations and stochasticity
Modelling of cell choice between differentiation and apoptosis on the basis of intracellular and extracellular regulations and stochasticity M. Banerjee 1, V. Volpert 2 arxiv:163.8712v1 [q-bio.cb] 29 Mar
More informationOn 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 informationAnalysis of Duopoly Output Game With Different Decision-Making Rules
Management Science and Engineering Vol. 9, No. 1, 015, pp. 19-4 DOI: 10.3968/6117 ISSN 1913-0341 [Print] ISSN 1913-035X [Online] www.cscanada.net www.cscanada.org Analysis of Duopoly Output Game With Different
More informationEffects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of Non-degenerate Cascade Two-Photon Lasers
Commun. Theor. Phys. Beijing China) 48 2007) pp. 288 294 c International Academic Publishers Vol. 48 No. 2 August 15 2007 Effects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of
More informationα 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 informationGeneralized projective synchronization between two chaotic gyros with nonlinear damping
Generalized projective synchronization between two chaotic gyros with nonlinear damping Min Fu-Hong( ) Department of Electrical and Automation Engineering, Nanjing Normal University, Nanjing 210042, China
More informationAvailable online at ScienceDirect. Procedia IUTAM 19 (2016 ) IUTAM Symposium Analytical Methods in Nonlinear Dynamics
Available online at www.sciencedirect.com ScienceDirect Procedia IUTAM 19 (2016 ) 11 18 IUTAM Symposium Analytical Methods in Nonlinear Dynamics A model of evolutionary dynamics with quasiperiodic forcing
More informationDYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD
Journal of Sustainability Science Management Volume 10 Number 2, December 2015: 16-23 ISSN: 1823-8556 Penerbit UMT DYNAMICS OF A PREDATOR-PREY INTERACTION IN CHEMOSTAT WITH VARIABLE YIELD SARKER MD SOHEL
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 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 informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationSuppression 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 informationExact Solutions for the Nonlinear (2+1)-Dimensional Davey-Stewartson Equation Using the Generalized ( G. )-Expansion Method
Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Published by Canadian Center of Science and Education Exact Solutions for the Nonlinear +-Dimensional Davey-Stewartson Equation
More informationDynamics 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 informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationSome Dynamical Behaviors In Lorenz Model
International Journal Of Computational Engineering Research (ijceronline.com) Vol. Issue. 7 Some Dynamical Behaviors In Lorenz Model Dr. Nabajyoti Das Assistant Professor, Department of Mathematics, Jawaharlal
More informationDELAY DIFFERENTIAL EQUATION MODELING FOR THE ANALYSIS OF HEMATOLOGICAL DISEASES BY ALTERNATIVE G-CSF TREATMENT
DELAY DIFFERENTIAL EQUATION MODELING FOR THE ANALYSIS OF HEMATOLOGICAL DISEASES BY ALTERNATIVE G-CSF TREATMENT A THESIS REPORT Submitted by BALAMURALITHARAN S. Under the guidance of Dr. S.RAJASEKARAN in
More informationHopf Bifurcation Analysis of Pathogen-Immune Interaction Dynamics With Delay Kernel
Math. Model. Nat. Phenom. Vol. 2, No. 1, 27, pp. 44-61 Hopf Bifurcation Analysis of Pathogen-Immune Interaction Dynamics With Delay Kernel M. Neamţu a1, L. Buliga b, F.R. Horhat c, D. Opriş b a Department
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 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 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 informationOSCILLATIONS IN A MATURATION MODEL OF BLOOD CELL PRODUCTION
SIAM J. APPL. MATH. Vol. 66, No. 6, pp. 227 248 c 26 Society for Industrial and Applied Mathematics OSCILLATIONS IN A MATURATION MODEL OF BLOOD CELL PRODUCTION IVANA DROBNJAK, A. C. FOWLER, AND MICHAEL
More informationSTABILITY AND BIFURCATION ANALYSIS ON DELAY DIFFERENTIAL EQUATIONS. Xihui Lin. Master of Science. Applied Mathematics
University of Alberta STABILITY AND BIFURCATION ANALYSIS ON DELAY DIFFERENTIAL EQUATIONS by Xihui Lin A thesis submitted to the Faculty of Graduate Studies and Research in partial fulfillment of the requirements
More informationApproximating the Stability Region for a Differential Equation with a Distributed Delay
Approximating the Stability Region for a Differential Equation with a Distributed Delay S.A. Campbell a and R. Jessop a a Department of Applied Mathematics, University of Waterloo, Waterloo, NL 3G, Canada
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature14242 Supplementary Methods Mathematical modelling Hematopoietic fluxes. Consider two successive cell compartments in hematopoiesis, an upstream and a reference
More informationThe exponential asymptotic stability of singularly perturbed delay differential equations with a bounded lag
J. Math. Anal. Appl. 270 (2002) 143 149 www.academicpress.com The exponential asymptotic stability of singularly perturbed delay differential equations with a bounded lag Hongjiong Tian Department of Mathematics,
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 informationLocal stability and Hopf bifurcation analysis for Compound TCP
This article has been accepted for publication in a future issue of this journal, but has not been fully edited Content may change prior to final publication Citation information: DOI 1119/TCNS17747839,
More informationA Mathematical Model of Chronic Myelogenous Leukemia
A Mathematical Model of Chronic Myelogenous Leukemia Brent Neiman University College Oxford University Candidate for Master of Science in Mathematical Modelling and Scientific Computing MSc Dissertation
More informationRate Equation Model for Semiconductor Lasers
Rate Equation Model for Semiconductor Lasers Prof. Sebastian Wieczorek, Applied Mathematics, University College Cork October 21, 2015 1 Introduction Let s consider a laser which consist of an optical resonator
More informationThe 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 informationReaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations
Reaction-Diffusion Models and Bifurcation Theory Lecture 11: Bifurcation in delay differential equations Junping Shi College of William and Mary, USA Collaborators/Support Shanshan Chen (PhD, 2012; Harbin
More informationTHREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS
THREE PROBLEMS OF RESONANCE IN COUPLED OR DRIVEN OSCILLATOR SYSTEMS A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the
More informationDynamical modelling of systems of coupled oscillators
Dynamical modelling of systems of coupled oscillators Mathematical Neuroscience Network Training Workshop Edinburgh Peter Ashwin University of Exeter 22nd March 2009 Peter Ashwin (University of Exeter)
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
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 informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
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 informationApplication 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 informationQualitative Analysis of a Discrete SIR Epidemic Model
ISSN (e): 2250 3005 Volume, 05 Issue, 03 March 2015 International Journal of Computational Engineering Research (IJCER) Qualitative Analysis of a Discrete SIR Epidemic Model A. George Maria Selvam 1, D.
More informationNumerical Solution of Second Order Singularly Perturbed Differential Difference Equation with Negative Shift. 1 Introduction
ISSN 749-3889 print), 749-3897 online) International Journal of Nonlinear Science Vol.8204) No.3,pp.200-209 Numerical Solution of Second Order Singularly Perturbed Differential Difference Equation with
More informationControlling the Period-Doubling Bifurcation of Logistic Model
ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.20(2015) No.3,pp.174-178 Controlling the Period-Doubling Bifurcation of Logistic Model Zhiqian Wang 1, Jiashi Tang
More informationAnalysis of Bifurcations in a Power System Model with Excitation Limits
Analysis of Bifurcations in a Power System Model with Excitation Limits Rajesh G. Kavasseri and K. R. Padiyar Department of Electrical Engineering Indian Institute of Science, Bangalore, India Abstract
More informationA Structured Population Model for Erythropoiesis
A Structured Population Model for Erythropoiesis Doris H. Fuertinger 1, Franz Kappel 1, Peter Kotanko 2, Nathan W. Levine 2, Stephan Thijssen 2 1 Institute for Mathematics and Scientific Computation, University
More informationLong Period Oscillations in a G 0 Model of Hematopoietic Stem Cells
SIAM J. APPLIED DYNAMICAL SYSTEMS Vol. 0, No. 0, pp. 000 000 c 2004 Society for Industrial and Applied Mathematics Long Period Oscillations in a G 0 Model of Hematopoietic Stem Cells Laurent Pujo-Menjouet,
More informationFormation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas )
Formation and Long Term Evolution of an Externally Driven Magnetic Island in Rotating Plasmas ) Yasutomo ISHII and Andrei SMOLYAKOV 1) Japan Atomic Energy Agency, Ibaraki 311-0102, Japan 1) University
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 information