Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps

Size: px
Start display at page:

Download "Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps"

Transcription

1 Numerical Continuation and Normal Form Analysis of Limit Cycle Bifurcations without Computing Poincaré Maps Yuri A. Kuznetsov joint work with W. Govaerts, A. Dhooge(Gent), and E. Doedel (Montreal) LCBIF p.1/47

2 Contents References Limit cycles of ODEs and their local bifurcations Standard approach using Poincaré maps Continuation and normal form analysis using BVPs Open problems LCBIF p.2/47

3 1. References Kuznetsov, Yu. A., Govaerts, W., Doedel, E. J., and Dhooge, A. Numerical periodic normalization for codim 1 bifurcations of limit cycles. SIAM J. Numer. Analysis 43 (2005), Govaerts, W., Kuznetsov, Yu.A., and Dhooge, A. Numerical continuation of bifurcations of limit cycles in MATLAB. SIAM J. Sci. Comput. 27 (2005), Doedel, E.J., Govaerts, W., Kuznetsov, Yu.A., and Dhooge, A. Numerical continuation of branch points of equilibria and periodic orbits. Int. J. Bifurcation & Chaos 15 (2005), Doedel, E.J., Govaerts, W., and Kuznetsov, Yu.A. Computation of periodic solution bifurcations in ODEs using bordered systems. SIAM J. Numer. Analysis 41 (2003), LCBIF p.3/47

4 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Re µ LCBIF p.4/47

5 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Re µ LCBIF p.4/47

6 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Flip (PD): µ 1 = 1; Re µ LCBIF p.4/47

7 2. Limit cycles of ODEs and their local bifurcations dx dt = f(x, α), x Rn, α R m. A limit cycle C 0 corresponds to a periodic solution x 0 (t + T 0 ) = x 0 (t) and has Floquet multipliers {µ 1, µ 2,..., µ n 1, µ n = 1} = σ(m(t 0 )), where Ṁ(t) f x (x 0 (t), α 0 )M(t) = 0, M(0) = I n. Im µ Fold (LPC): µ 1 = 1; Flip (PD): µ 1 = 1; Torus (NS): µ 1,2 = e ±iθ 0. θ 0 Re µ LCBIF p.4/47

8 2.1. Generic LPC bifurcation: µ 1 = 1 C 1 C 2 C 0 W c β W c 0 W c β β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.5/47

9 2.2. Generic PD bifurcation: µ 1 = 1 W c β C 0 W0 c Wβ c C 1 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.6/47

10 2.3. Generic NS bifurcation: µ 1,2 = e ±iθ 0 C 0 T 2 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.7/47

11 3. Standard approach using Poincaré maps Poincaré map and its fixed points Continuation of PD and LP bifurcations of maps Continuation of NS bifurcation of maps Smooth normal forms of maps on center manifolds Difficulties using Poincaré maps Reference on finite-dimensional bordering techniques: Govaerts, W. Numerical Methods for Bifurcations of Dynamical Equilibria. SIAM, LCBIF p.8/47

12 3.1. Poincaré map and its fixed points Poincaré map P : R n 1 R m R n 1 A limit cycle C 0 corresponds to a fixed point: W c α y y 0 P(y, α) P(y 0, α) y 0 = 0. C 0 Σ Write P(y 0 + η, α 0 ) = y 0 + Aη B(η, η) C(η, η, η) + O( η 4 ), where A = P y (y 0, α 0 ) and, for i = 1, 2,..., n 1, B i (η, ζ) = C i (η, ζ, ω) = n 1 j,k=1 n 1 j,k,l=1 2 P i (y, α 0 ) y j y k 3 P i (y, α 0 ) y j y k y l η j ζ k, y=y0 η j ζ k ω l. y=y0 LCBIF p.9/47

13 3.2. Continuation of PD and LC bifurcations of maps Defining system: (y, α) R n 1 R 2 P(y, α) y = 0, g(y, α) = 0, where g is defined by solving P y(y, α) ± I n 1 w 1 v1 T 0 v g = 0 1 Vectors v 1, w 1 R n 1 are adapted to make the linear system nonsingular. (g y, g α ) can be computed efficiently using the adjoint linear system. LCBIF p.10/47

14 3.3. Continuation of NS bifurcation of maps Defining system: (y, α, κ) R n 1 R 2 R P(y, α) y = 0, g 11 (y, α, κ) = 0, g 22 (y, α, κ) = 0, where κ = cos θ 0 and g ij are defined by solving I n 1 2κP y (y, α) + [P y (y, α)] 2 w 1 w 2 r s v1 T 0 0 g 11 g 12 v2 T 0 0 g 21 g 22 = Vectors v 1,2, w 1,2 R n 1 are adapted to ensure unique solvability. Efficient computation of derivatives of g ij. LCBIF p.11/47

15 3.4. Smooth normal forms of maps on center manifolds Eigenvectors Normal form Critical coefficients LP Aq = q A T p = p p, q = 1 ξ β + ξ + bξ 2 +O(ξ 3 ), ξ R b = 1 2 p, B(q, q) PD Aq = q A T p = p p, q = 1 ξ (1 + β)ξ + cξ 3 +O(ξ 4 ), ξ R c = 1 6 p, C(q, q, q) + 3B(q, h 2) h 2 = (I n A) 1 B(q, q) NS Aq = e iθ 0 q A T p = e iθ 0 p e iνθ0 1 ν = 1, 2, 3, 4 ξ ξe iθ ( 1 + β + d ξ 2 ) + O( ξ 4 ), ξ C d = 1 2 e iθ 0 p, C(q, q, q) + 2B(q, h 11 ) + B( q, h 20 ) h 11 = (I n A) 1 B(q, q) p, q = 1 h 20 = (e 2iθ 0 I n A) 1 B(q, q) LCBIF p.12/47

16 LP-coefficient b Assume y 0 = 0 and write P(H) = AH B(H, H) + O( H 3 ), and locally represent the center manifold W0 c H : R R n 1, as the graph of a function u = H(ξ) = ξq h 2ξ 2 + O(ξ 3 ), ξ R, h 2 R n 1. The restriction of the Poincaré map to W0 c is ξ G(ξ) = ξ + bξ 2 + O(ξ 3 ). The invariance of the center manifold W0 c means P(H(ξ)) = H(G(ξ)). LCBIF p.13/47

17 A(ξq h 2ξ 2 ) B(ξq, ξq) + O( ξ 3 ) = (ξ + bξ 2 )q h 2ξ 2 + O( ξ 3 ) The ξ-terms give the identity: Aq = q. The ξ 2 -terms give the equation for h 2 : (A I n 1 )h 2 = B(q, q) + 2 bq. It is singular and its Fredholm solvability implies where A T p = p, p, q = 1. 1 b = p, B(q, q), 2 LCBIF p.14/47

18 3.5. Difficulties using Poincaré maps The computation of P can be unstable numerically near saddle or repelling limit cycles. Finite differences give low accuracy for B and C even if the cycle C 0 is stable. Simultaneous solving variational equations for the components of A, B, and C is possible but costly and difficult to implement. LCBIF p.15/47

19 4. Continuation and normal form analysis using BVPs Continuation of limit cycles in one parameter Simple bifurcation points Continuation of bifurcations in two parameters Periodic normalization on center manifolds LCBIF p.16/47

20 4.1. Continuation of limit cycles in one parameter Defining system (BVP with IC) [Keller & Doedel, 1981]: u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, where ũ is a reference periodic solution. Linearization w.r.t. (u, T, α): D T f x (u, α) f(u, α) T f α (u, α) δ 0 δ Int ũ 0 0 LCBIF p.17/47

21 Discretization via orthogonal collocation Mesh points: 0 = τ 0 < τ 1 < < τ N = 1. Basis points: τ i,j = τ i + j m (τ i+1 τ i ), i = 0, 1,..., N 1, j = 0, 1,..., m. Approximation: u (i) (τ) = m j=0 ui,j l i,j (τ), τ [τ i, τ i+1 ], where l i,j (τ) are the Lagrange basis polynomials and u i,m = u i+1,0. Collocation [de Boor & Swartz, 1973]: ( m j=0 ui,j l i,j i,k)) (ζ T f( m j=0 ui,j l i,j (ζ i,k ), α) = 0, u 0,0 u N 1,m = 0, m j=0 σ i,j ũ i,j, u i,j = 0, N 1 i=0 where ζ i,k, k = 1, 2,..., m, are the Gauss points. LCBIF p.18/47

22 Sparse Jacobian matrix u 0,0 u 0,1 u 1,0 u 1,1 u 2,0 u 2,1 u 3,0 T α LCBIF p.19/47

23 4.2. Simple bifurcation points Φ(τ) T f x (u(τ), α 0 )Φ(τ) = 0, Φ(0) = I n, Ψ(τ) + T fx T (u(τ), α 0 )Ψ(τ) = 0, Ψ(0) = I n. LPC: (Φ(1) I n )q 0 = 0, (Φ(1) I n )q 1 = q 0, (Ψ(1) I n )p 0 = 0, (Ψ(1) I n )p 1 = p 0. PD: (Φ(1) + I n )q 2 = 0, (Ψ(1) + I n )p 2 = 0. NS: κ = cos θ 0 (Φ(1) e iθ 0 I n )(q 3 + iq 4 ) = 0, (Ψ(1) e iθ 0 I n )(p 3 + ip 4 ) = 0. We have (I n 2κΦ(1) + Φ 2 (1))q 3,4 = 0. LCBIF p.20/47

24 4.3. Continuation of bifurcations in two parameters PD and LPC: (u, T, α) C 1 ([0, 1], R n ) R R 2 u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, G[u, T, α] = 0. NS: (u, T, α, κ) C 1 ([0, 1], R n ) R R 2 R u(τ) T f(u(τ), α) = 0, τ [0, 1], u(0) u(1) = 0, 1 0 ũ(τ), u(τ) dτ = 0, G 11 [u, T, α, κ] = 0, G 22 [u, T, α, κ] = 0. LCBIF p.21/47

25 PD-continuation There exist v 01, w 01 C 0 ([0, 1], R n ), and w 02 R n, such that N 1 : C 1 ([0, 1], R n ) R C 0 ([0, 1], R n ) R n R, D T f x (u, α) w 01 N 1 = δ 0 δ 1 w 02, Int v01 0 is one-to-one and onto near a simple PD bifurcation point. Define G by solving N 1 v 0 = 0 G. 1 LCBIF p.22/47

26 The BVP for (v, G) can be written in the classical form v(τ) T f x (u(τ), α)v(τ) + Gw 01 (τ) = 0, τ [0, 1], v(0) + v(1) + Gw 02 = 0, 1 0 v 01(τ), v(τ) dτ 1 = 0. One can take and w 02 = 0 w 01 (τ) = Ψ(τ)p 2, v 01 (τ) = Φ(τ)q 2. LCBIF p.23/47

27 LPC-continuation There exist v 01, w 01 C 0 ([0, 1], R n ), w 02 R n, and v 02, w 03 R such that N 2 : C 1 ([0, 1], R n ) R 2 C 0 ([0, 1], R n ) R n R 2, D T f x (u, α) f(u, α) w 01 δ 0 δ 1 0 w 02 N 2 =, Int f(u,α) 0 w 03 Int v01 v 02 0 is one-to-one and onto near a simple LPC bifurcation point. 0 v Define G by solving N 2 S = 0. 0 G 1 LCBIF p.24/47

28 NS-continuation There exist v 01, v 02, w 11, w 12 C 0 ([0, 2], R n ), and w 21, w 22 R n, such that N 3 : C 1 ([0, 2], R n ) R 2 C 0 ([0, 2], R n ) R n R 2, D T f x (u, α) w 11 w 12 δ 0 2κδ 1 + δ 2 w 21 w 22 N 3 =, Int v Int v is one-to-one and onto near a simple NS bifurcation point. 0 0 r s Define G jk by solving N 3 G 11 G 12 = G 21 G LCBIF p.25/47

29 Remarks on continuation of bifurcations After discretization via orthogonal collocation, all linear BVPs for G s have sparsity structure that is identical to that of the linearization of the BVP for limit cycles. For each defining system holds: Simplicity of the bifurcation + Transversality Regularity of the defining BVP. Jacobian matrix of each (discretized) defining BVP can be efficiently computed using adjoint linear BVP. Border adaptation using solutions of the adjoint linear BVPs. Actually implemented in MATCONT, also with compiled C-codes for the Jacobian matrices. LCBIF p.26/47

30 4.4. Periodic normalization on center manifolds Parameter-dependent periodic normal forms for LPC, PD, and NS Critical normal form coefficients References on periodic normal forms: Iooss, G. Global characterization of the normal form for a vector field near a closed orbit. J. Diff. Equations 76 (1988), Iooss, G. and Adelmeyer, M. Topics in Bifurcation Theory and Applications. World Scientific, LCBIF p.27/47

31 LPC bifurcation T 0 -periodic normal form on Wα: c dτ dt = 1 + ν(α) ξ + a(α)ξ 2 + O(ξ 3 ), dξ dt = β(α) + b(α)ξ 2 + O(ξ 3 ), where a, b R, sign(b) = sign( b). C 1 C 2 C 0 W c β W c 0 W c β β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.28/47

32 PD bifurcation 2T 0 -periodic normal form on Wα: c dτ dt = 1 + ν(α) + a(α)ξ 2 + O(ξ 4 ), dξ dt = β(α)ξ + c(α)ξ 3 + O(ξ 4 ), where a, c R, sign(c) = sign( c). W c β C 0 W0 c Wβ c C 1 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.29/47

33 NS bifurcation T 0 -periodic normal form on Wα: c dτ dt dξ dt = 1 + ν(α) + a(α) ξ 2 + O( ξ 4 ), ( = β(α) + iθ(α) T (α) where a R, d C, sign(re(d)) = sign(re( d)). ) ξ + d(α)ξ ξ 2 + O( ξ 4 ), C 0 T 2 β(α) < 0 β(α) = 0 β(α) > 0 LCBIF p.30/47

34 Critical normal form coefficients At a codimension-one point write f(x 0 (t)+v, α 0 ) = f(x 0 (t), α 0 )+A(t)v+ 1 2 B(t; v, v)+1 6 C(t; v, v, v)+o( v 4 ), where A(t) = f x (x 0 (t), α 0 ) and the components of the multilinear functions B and C are given by and B i (t; u, v) = C i (t; u, v, w) = n j,k=1 n j,k,l=1 2 f i (x, α 0 ) x j x k 3 f i (x, α 0 ) x j x k x l for i = 1, 2,..., n. These are T 0 -periodic in t. u j v k x=x0 (t) u j v k w l, x=x0 (t) LCBIF p.31/47

35 Fold (LPC): µ 1 = 1 Critical center manifold W c 0 : τ [0, T 0], ξ R x = x 0 (τ) + ξv(τ) + H(τ, ξ), where H(T 0, ξ) = H(0, ξ), τ C 0 H(τ, ξ) = 1 2 h 2(τ)ξ 2 + O(ξ 3 ) Critical periodic normal form on W c 0 : dτ dt dξ dt = 1 ξ + aξ 2 + O(ξ 3 ), = bξ 2 + O(ξ 3 ), where a, b R, while the O(ξ 3 )-terms are T 0 -periodic in τ. ξ W c 0 LCBIF p.32/47

36 LPC: Eigenfunctions implying where ϕ satisfies v(τ) A(τ)v(τ) f(x 0 (τ), α 0 ) = 0, τ [0, T 0 ], v(0) v(t 0 ) = 0, T0 0 v(τ), f(x 0 (τ), α 0 ) dτ = 0, T0 0 ϕ (τ), f(x 0 (τ), α 0 ) dτ = 0, ϕ (τ) + A T (τ)ϕ (τ) = 0, τ [0, T 0 ], ϕ (0) ϕ (T 0 ) = 0, T0 0 ϕ (τ), v(τ) dτ 1 = 0. LCBIF p.33/47

37 LPC: Computation of b Substitute into Collect dx dt = x ξ dξ dt + x τ dτ dt ξ 0 : ẋ 0 = f(x 0, α 0 ), ξ 1 : v A(τ)v = ẋ 0, ξ 2 : ḣ 2 A(τ)h 2 = B(τ; v, v) 2af(x 0, α 0 ) + 2 v 2bv. Fredholm solvability condition b = 1 2 T0 0 ϕ (τ), B(τ; v(τ), v(τ)) + 2A(τ)v(τ) dτ. LCBIF p.34/47

38 LPC: Example in MATCONT (ABC-reaction model) u 1 = u 1 + p 1 (1 u 1 )e u 3, u 2 = u 2 + p 1 (1 u 1 p 5 u 2 )e u 3, u 3 = u 3 p 3 u 3 + p 1 p 4 (1 u 1 + p 2 p 5 u 2 )e u 3. p 3 = 1.5 p 4 = 8.0 p 5 = CPC 1.1 p b = 0 at CPC (cusp) p 1 LCBIF p.35/47

39 Flip (PD): µ 1 = 1 Critical center manifold W c 0 : τ [0, 2T 0], ξ R x = x 0 (τ) + ξw(τ) + H(τ, ξ), where H(2T 0, ξ) = H(0, ξ), τ C 0 H(τ, ξ) = 1 2 h 2(τ)ξ h 3(τ)ξ 3 + O(ξ 4 ) ξ Critical periodic normal form on W c 0 : dτ dt dξ dt = 1 + aξ 2 + O(ξ 4 ), = cξ 3 + O(ξ 4 ), where a, c R, while the O(ξ 4 )-terms are 2T 0 -periodic in τ. W c 0 LCBIF p.36/47

40 PD: Eigenfunctions w(τ) = v(τ), τ [0, T 0 ], v(τ T 0 ), τ [T 0, 2T 0 ], w (τ) = v (τ), τ [0, T 0 ], v (τ T 0 ), τ [T 0, 2T 0 ], with v(τ) A(τ)v(τ) = 0, τ [0, T 0 ], v(0) + v(t 0 ) = 0, T0 0 v(τ), v(τ) dτ 1 = 0, v (τ) + A T (τ)v (τ) = 0, τ [0, T 0 ], v (0) + v (T 0 ) = 0, T0 0 v (τ), v(τ) dτ 1/2 = 0. LCBIF p.37/47

41 PD: Quadratic terms ξ 2 : ḣ 2 A(τ)h 2 = B(τ; w, w) 2aẋ 0, τ [0, 2T 0 ]. Since Ker ( d dτ A(τ)) = span{w, ψ = ẋ 0 }, we must have where ψ satisfies 2T0 0 w (τ), B(τ; w(τ), w(τ)) 2aẋ 0 (τ) dτ = 0, 2T0 0 ψ (τ), B(τ; w(τ), w(τ)) 2aẋ 0 (τ) dτ = 0, ψ (τ) + A T (τ)ψ (τ) = 0, τ [0, T 0 ], ψ (0) ψ (T 0 ) = 0, T0 0 ψ (τ), f(x 0 (τ), α 0 ) dτ 1/2 = 0, and is extended to [T 0, 2T 0 ] by periodicity. LCBIF p.38/47

42 PD: Computation of a and h 2 a = 1 2 The first Fredholm condition holds identically for all a, while the second gives 2T0 0 ψ (τ), B(τ; w(τ), w(τ) dτ = T0 0 ψ (τ), B(τ; v(τ), v(τ) dτ. Define h 2 on [0, T 0 ] as the unique solution to ḣ 2 (τ) A(τ)h 2 (τ) B(τ; v(τ), v(τ)) + 2af(x 0 (τ), α 0 ) = 0, h 2 (0) h 2 (T 0 ) = 0, T0 0 ψ (τ), h 2 (τ) dτ = 0, and extend it by periodicity to [T 0, 2T 0 ]. LCBIF p.39/47

43 PD: Computation of c Cubic terms: ξ 3 : ḣ 3 A(τ)h 3 = C(τ; w, w, w) + 3B(τ; w, h 2 ) 6aẇ 6cw. The Fredholm solvability condition implies 6c = or 2T0 0 2T0 0 w (τ), C(τ; w(τ), w(τ), w(τ)) + 3B(τ; w(τ), h 2 (τ)) dτ w (τ), 6aA(τ)w(τ) dτ c = 1 3 T0 0 v (τ), C(τ; v(τ), v(τ), v(τ))+3b(τ; v(τ), h 2 (τ)) 6aA(τ)v(τ) dτ LCBIF p.40/47

44 Torus (NS): µ 1,2 = e ±iθ 0 with e iνθ0 1, ν = 1, 2, 3, 4 Critical center manifold W c 0 : τ [0, T 0], ξ C x = x 0 (τ) + ξv(τ) + ξ v(τ) + H(τ, ξ, ξ), H(T 0, ξ, ξ) = H(0, ξ, ξ), H(τ, ξ, ξ) = 1 2 h 20(τ)ξ 2 + h 11 (τ)ξ ξ h 02(τ) ξ h 30(τ)ξ h 21(τ)ξ 2 ξ h 12(τ)ξ ξ h 03(τ) ξ 3 + O( ξ 4 ). Critical periodic normal form on W0 c: dτ = 1 + a ξ 2 + O( ξ 4 ), dt dξ = iθ 0 ξ + dξ ξ 2 + O( ξ 4 ), dt T 0 where a R, d C, and the O( ξ 4 )-terms are T 0 -periodic in τ. LCBIF p.41/47

45 NS: Complex eigenfunctions and v(τ) A(τ)v(τ) + iθ 0 T 0 v(τ) = 0, τ [0, T 0 ], v(0) v(t 0 ) = 0, T0 0 v(τ), v(τ) dτ 1 = 0. v (τ) + A T (τ)v (τ) iθ 0 T 0 v (τ) = 0, τ [0, T 0 ], v (0) v (T 0 ) = 0, T0 0 v (τ), v(τ) dτ 1 = 0. LCBIF p.42/47

46 NS: Quadratic terms ξ 2 ξ0 : ḣ 20 A(τ)h iθ 0 T 0 h 20 = B(τ; v, v) Since e 2iθ 0 is not a multiplier of the critical cycle, the BVP ḣ 20 (τ) A(τ)h 20 (τ) + 2iθ 0 T 0 h 20 (τ) B(τ; v(τ), v(τ)) = 0, has a unique solution on [0, T 0 ]. ξ 2 : ḣ 11 A(τ)h 11 = B(τ; v, v) aẋ 0 Here Ker ( ) d dτ A(τ) h 20 (0) h 20 (T 0 ) = 0. = span(ϕ = ẋ 0 ). LCBIF p.43/47

47 NS: Computation of a and h 11 Define ϕ as the unique solution of ϕ (τ) + A T (τ)ϕ (τ) = 0, τ [0, T 0 ], ϕ (0) ϕ (T 0 ) = 0, T0 0 ϕ (τ), f(x 0 (τ), α 0 ) dτ 1 = 0. Fredholm solvability: a = T 0 0 ϕ (τ), B(τ; v(τ), v(τ) dτ. Then find h 11 on [0, T 0 ] from the BVP ḣ 11 (τ) A(τ)h 11 (τ) B(τ; v(τ), v(τ)) + af(x 0 (τ), α 0 ) = 0, h 11 (0) h 11 (T 0 ) = 0, T0 0 ϕ (τ), h 11 (τ) dτ = 0. LCBIF p.44/47

48 NS: Computation of d Cubic terms: ξ 2 ξ : ḣ21 Ah 21 + iθ 0 T 0 h 21 = 2B(τ; h 11, v) + B(τ; h 20, v) Fredholm solvability condition: + C(τ; v, v, v) 2a v 2dv. d = 1 2 T v (τ), C(τ; v(τ), v(τ), v(τ)) dτ T0 0 a v (τ), B(τ; h 11 (τ), v(τ)) + B(τ; h 20 (τ), v(τ)) dτ T0 0 v (τ), A(τ)v(τ) dτ + iaθ 0 T 0. LCBIF p.45/47

49 Remarks on numerical periodic normalization Only the derivatives of f(x, α 0 ) are used, not those of the Poincaré map P(y, α 0 ). Detection of codim 2 points is easy. After discretization via orthogonal collocation, all linear BVPs involved have the standard sparsity structure. One can re-use solutions to linear BVPs appearing in the continuation to compute the normal form coefficients. Actually implemented in MATCONT. LCBIF p.46/47

50 5. Open problems Automatic differentiation of the Poincaré map vs. BVPs. Periodic normal forms for codim 2 bifurcations of limit cycles. Branch switching at codim 2 points to limit cycle codim 1 continuation. LCBIF p.47/47

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

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

More information

BIFURCATION PHENOMENA Lecture 4: Bifurcations in n-dimensional ODEs

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

More information

Lecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL)

Lecture 5. Numerical continuation of connecting orbits of iterated maps and ODEs. Yu.A. Kuznetsov (Utrecht University, NL) Lecture 5 Numerical continuation of connecting orbits of iterated maps and ODEs Yu.A. Kuznetsov (Utrecht University, NL) May 26, 2009 1 Contents 1. Point-to-point connections. 2. Continuation of homoclinic

More information

Continuation of cycle-to-cycle connections in 3D ODEs

Continuation of cycle-to-cycle connections in 3D ODEs HET p. 1/2 Continuation of cycle-to-cycle connections in 3D ODEs Yuri A. Kuznetsov joint work with E.J. Doedel, B.W. Kooi, and G.A.K. van Voorn HET p. 2/2 Contents Previous works Truncated BVP s with projection

More information

( Bx + f(x, y) Cy + g(x, y)

( Bx + f(x, y) Cy + g(x, y) Chapter 6 Center manifold reduction The previous chaper gave a rather detailed description of bifurcations of equilibria and fixed points in generic one-parameter families of ODEs and maps with minimal

More information

An Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013

An Introduction to Numerical Continuation Methods. with Application to some Problems from Physics. Eusebius Doedel. Cuzco, Peru, May 2013 An Introduction to Numerical Continuation Methods with Application to some Problems from Physics Eusebius Doedel Cuzco, Peru, May 2013 Persistence of Solutions Newton s method for solving a nonlinear equation

More information

Numerical Continuation of Bifurcations - An Introduction, Part I

Numerical Continuation of Bifurcations - An Introduction, Part I Numerical Continuation of Bifurcations - An Introduction, Part I given at the London Dynamical Systems Group Graduate School 2005 Thomas Wagenknecht, Jan Sieber Bristol Centre for Applied Nonlinear Mathematics

More information

IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION

IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION International Journal of Bifurcation and Chaos c World Scientific Publishing Company IMPROVED HOMOCLINIC PREDICTOR FOR BOGDANOV-TAKENS BIFURCATION YU.A. KUZNETSOV, H.G.E. MEIJER Department of Applied Mathematics,

More information

The fold-flip bifurcation

The fold-flip bifurcation The fold-flip bifurcation Yu.A. Kuznetsov, H.G.E. Meijer, and L. van Veen March 0, 003 Abstract The fold-flip bifurcation occurs if a map has a fixed point with multipliers + and simultaneously. In this

More information

An Introduction to Numerical Continuation Methods. with Applications. Eusebius Doedel IIMAS-UNAM

An Introduction to Numerical Continuation Methods. with Applications. Eusebius Doedel IIMAS-UNAM An Introduction to Numerical Continuation Methods with Applications Eusebius Doedel IIMAS-UNAM July 28 - August 1, 2014 1 Persistence of Solutions Newton s method for solving a nonlinear equation [B83]

More information

computer session iv: Two-parameter bifurcation analysis of equilibria and limit cycles with matcont

computer session iv: Two-parameter bifurcation analysis of equilibria and limit cycles with matcont computer session iv: Two-parameter bifurcation analysis of equilibria and limit cycles with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht May 16,

More information

Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks

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

BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs

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

More information

LINEAR SECOND-ORDER EQUATIONS

LINEAR SECOND-ORDER EQUATIONS LINEAR SECON-ORER EQUATIONS Classification In two independent variables x and y, the general form is Au xx + 2Bu xy + Cu yy + u x + Eu y + Fu + G = 0. The coefficients are continuous functions of (x, y)

More information

Connecting Orbits with Bifurcating End Points

Connecting Orbits with Bifurcating End Points Connecting Orbits with Bifurcating End Points Thorsten Hüls Fakultät für Mathematik Universität Bielefeld Postfach 100131, 33501 Bielefeld Germany huels@mathematik.uni-bielefeld.de June 24, 2004 Abstract

More information

Improved Homoclinic Predictor for Bogdanov Takens Bifurcation

Improved Homoclinic Predictor for Bogdanov Takens Bifurcation International Journal of Bifurcation and Chaos, Vol. 24, No. 4 (214) 14557 (12 pages) c World Scientific Publishing Company DOI: 1.1142/S2181274145576 Improved Homoclinic Predictor for Bogdanov Takens

More information

Homoclinic saddle to saddle-focus transitions in 4D systems

Homoclinic saddle to saddle-focus transitions in 4D systems Faculty of Electrical Engineering, Mathematics & Computer Science Homoclinic saddle to saddle-focus transitions in 4D systems Manu Kalia M.Sc. Thesis July 2017 Assessment committee: Prof. Dr. S. A. van

More information

Computational Methods in Dynamical Systems and Advanced Examples

Computational Methods in Dynamical Systems and Advanced Examples and Advanced Examples Obverse and reverse of the same coin [head and tails] Jorge Galán Vioque and Emilio Freire Macías Universidad de Sevilla July 2015 Outline Lecture 1. Simulation vs Continuation. How

More information

Bifurcation Analysis of Large Equilibrium Systems in Matlab

Bifurcation Analysis of Large Equilibrium Systems in Matlab Bifurcation Analysis of Large Equilibrium Systems in Matlab David S. Bindel 1, James W. Demmel 2, Mark J. Friedman 3, Willy J.F. Govaerts 4, and Yuri A. Kuznetsov 5 1 Department of Electrical Engineering

More information

Elements of Applied Bifurcation Theory

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

More information

Alberto Bressan. Department of Mathematics, Penn State University

Alberto Bressan. Department of Mathematics, Penn State University Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i

More information

Department of mathematics MA201 Mathematics III

Department of mathematics MA201 Mathematics III Department of mathematics MA201 Mathematics III Academic Year 2015-2016 Model Solutions: Quiz-II (Set - B) 1. Obtain the bilinear transformation which maps the points z 0, 1, onto the points w i, 1, i

More information

tutorial ii: One-parameter bifurcation analysis of equilibria with matcont

tutorial ii: One-parameter bifurcation analysis of equilibria with matcont tutorial ii: One-parameter bifurcation analysis of equilibria with matcont Yu.A. Kuznetsov Department of Mathematics Utrecht University Budapestlaan 6 3508 TA, Utrecht February 13, 2018 1 This session

More information

Stable cycles in a Cournot duopoly model of Kopel

Stable cycles in a Cournot duopoly model of Kopel Journal of Computational and Applied Mathematics 18 (008 47 58 www.elsevier.com/locate/cam Stable cycles in a Cournot duopoly model of Kopel W. Govaerts a,, R. Khoshsiar Ghaziani a,b a Department of Applied

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,

More information

Elements of Applied Bifurcation Theory

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

More information

Transcendental cases in stability problem. Hamiltonian systems

Transcendental cases in stability problem. Hamiltonian systems of Hamiltonian systems Boris S. Bardin Moscow Aviation Institute (Technical University) Faculty of Applied Mathematics and Physics Department of Theoretical Mechanics Hamiltonian Dynamics and Celestial

More information

BIFURCATIONS AND STRANGE ATTRACTORS IN A CLIMATE RELATED SYSTEM

BIFURCATIONS AND STRANGE ATTRACTORS IN A CLIMATE RELATED SYSTEM dx dt DIFFERENTIAL EQUATIONS AND CONTROL PROCESSES N 1, 25 Electronic Journal, reg. N P23275 at 7.3.97 http://www.neva.ru/journal e-mail: diff@osipenko.stu.neva.ru Ordinary differential equations BIFURCATIONS

More information

Reproducing formulas associated with symbols

Reproducing formulas associated with symbols Reproducing formulas associated with symbols Filippo De Mari Ernesto De Vito Università di Genova, Italy Modern Methods of Time-Frequency Analysis II Workshop on Applied Coorbit space theory September

More information

CONTINUATION OF CONNECTING ORBITS IN 3D-ODES: (II) CYCLE-TO-CYCLE CONNECTIONS arxiv: v1 [math.ds] 1 Apr 2008

CONTINUATION OF CONNECTING ORBITS IN 3D-ODES: (II) CYCLE-TO-CYCLE CONNECTIONS arxiv: v1 [math.ds] 1 Apr 2008 CONTINUATION OF CONNECTING ORBITS IN 3D-ODES: (II) CYCLE-TO-CYCLE CONNECTIONS arxiv:0804.0179v1 [math.ds] 1 Apr 2008 E.J. DOEDEL 1, B.W. KOOI 2, YU.A. KUZNETSOV 3, and G.A.K. van VOORN 2 1 Department of

More information

Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions

Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions Department of Mathematics Mathematical Sciences Homoclinic Orbits of Planar Maps: Asymptotics and Mel nikov Functions A thesis submitted for the degree of Master of Science Author: Dirk van Kekem Project

More information

Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.

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

More information

Description of the extension ddebiftool_nmfm

Description of the extension ddebiftool_nmfm Description of the extension ddebiftool_nmfm Maikel Bosschaert, Bram Wage and Yuri Kuznetsov 16 September 2015 This addendum describes the normal form extension to DDE-BifTool, a bifurcation analysis toolbox

More information

Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation

Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 1, 65 71 Hopf Bifurcations in Problems with O(2) Symmetry: Canonical Coordinates Transformation Faridon AMDJADI Department

More information

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 9, Number 2, pp. 227 237 (2014) http://campus.mst.edu/adsa Elliptic Problems for Pseudo Differential Equations in a Polyhedral Cone

More information

A review of stability and dynamical behaviors of differential equations:

A review of stability and dynamical behaviors of differential equations: A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary

More information

Summary of topics relevant for the final. p. 1

Summary of topics relevant for the final. p. 1 Summary of topics relevant for the final p. 1 Outline Scalar difference equations General theory of ODEs Linear ODEs Linear maps Analysis near fixed points (linearization) Bifurcations How to analyze a

More information

Using MatContM in the study of a nonlinear map in economics

Using MatContM in the study of a nonlinear map in economics Journal of Phsics: Conference Series PAPER OPEN ACCESS Using MatContM in the stud of a nonlinear map in economics To cite this article: Niels Neirnck et al 016 J. Phs.: Conf. Ser. 69 0101 Related content

More information

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

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

More information

Global bifurcation diagram for the Kerner-Konhäuser traffic flow model

Global bifurcation diagram for the Kerner-Konhäuser traffic flow model International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.

More information

Lectures on Dynamical Systems. Anatoly Neishtadt

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

More information

x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ).

x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x, λ). Continuation of equilibria Problem Parameter-dependent ODE ẋ = f(x, λ), x R d, λ R, f smooth enough. Goal: compute ( follow ) equilibrium solutions as λ varies, i.e. compute solutions (x, λ) to 0 = f(x,

More information

On the Stability of the Best Reply Map for Noncooperative Differential Games

On the Stability of the Best Reply Map for Noncooperative Differential Games On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University

More information

LECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS. Fourth Edition, February by Tadeusz STYŠ. University of Botswana

LECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS. Fourth Edition, February by Tadeusz STYŠ. University of Botswana i LECTURE NOTES IN PARTIAL DIFFERENTIAL EQUATIONS Fourth Edition, February 2011 by Tadeusz STYŠ University of Botswana ii Contents 1 Solution of Partial Differential Equations 1 1.1 The General Solution

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

Changing coordinates to adapt to a map of constant rank

Changing coordinates to adapt to a map of constant rank Introduction to Submanifolds Most manifolds of interest appear as submanifolds of others e.g. of R n. For instance S 2 is a submanifold of R 3. It can be obtained in two ways: 1 as the image of a map into

More information

Bifurcation Analysis of Neuronal Bursters Models

Bifurcation Analysis of Neuronal Bursters Models Bifurcation Analysis of Neuronal Bursters Models B. Knowlton, W. McClure, N. Vu REU Final Presentation August 3, 2017 Overview Introduction Stability Bifurcations The Hindmarsh-Rose Model Our Contributions

More information

Optimal Transportation. Nonlinear Partial Differential Equations

Optimal Transportation. Nonlinear Partial Differential Equations Optimal Transportation and Nonlinear Partial Differential Equations Neil S. Trudinger Centre of Mathematics and its Applications Australian National University 26th Brazilian Mathematical Colloquium 2007

More information

8 A pseudo-spectral solution to the Stokes Problem

8 A pseudo-spectral solution to the Stokes Problem 8 A pseudo-spectral solution to the Stokes Problem 8.1 The Method 8.1.1 Generalities We are interested in setting up a pseudo-spectral method for the following Stokes Problem u σu p = f in Ω u = 0 in Ω,

More information

Introduction to Bifurcation and Normal Form theories

Introduction to Bifurcation and Normal Form theories Introduction to Bifurcation and Normal Form theories Romain Veltz / Olivier Faugeras October 9th 2013 ENS - Master MVA / Paris 6 - Master Maths-Bio (2013-2014) Outline 1 Invariant sets Limit cycles Stable/Unstable

More information

Stochastic Spectral Approaches to Bayesian Inference

Stochastic Spectral Approaches to Bayesian Inference Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to

More information

Construction of Lyapunov functions by validated computation

Construction of Lyapunov functions by validated computation Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The

More information

A Model of Evolutionary Dynamics with Quasiperiodic Forcing

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

More information

8.1 Bifurcations of Equilibria

8.1 Bifurcations of Equilibria 1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

More information

Distance Between Ellipses in 2D

Distance Between Ellipses in 2D Distance Between Ellipses in 2D David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To

More information

Deterministic Dynamic Programming

Deterministic Dynamic Programming Deterministic Dynamic Programming 1 Value Function Consider the following optimal control problem in Mayer s form: V (t 0, x 0 ) = inf u U J(t 1, x(t 1 )) (1) subject to ẋ(t) = f(t, x(t), u(t)), x(t 0

More information

UNIVERSIDADE DE SÃO PAULO

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

Math 489AB Exercises for Chapter 1 Fall Section 1.0

Math 489AB Exercises for Chapter 1 Fall Section 1.0 Math 489AB Exercises for Chapter 1 Fall 2008 Section 1.0 1.0.2 We want to maximize x T Ax subject to the condition x T x = 1. We use the method of Lagrange multipliers. Let f(x) = x T Ax and g(x) = x T

More information

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

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

More information

(V.C) Complex eigenvalues and other topics

(V.C) Complex eigenvalues and other topics V.C Complex eigenvalues and other topics matrix Let s take a closer look at the characteristic polynomial for a 2 2 namely A = a c f A λ = detλi A = λ aλ d bc b d = λ 2 a + dλ + ad bc = λ 2 traλ + det

More information

Model reduction for linear systems by balancing

Model reduction for linear systems by balancing Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,

More information

Global bifurcation diagram for the Kerner-Konhäuser traffic flow model

Global bifurcation diagram for the Kerner-Konhäuser traffic flow model International Journal of Bifurcation and Chaos c World Scientific Publishing Company Global bifurcation diagram for the Kerner-Konhäuser traffic flow model Joaquín Delgado Mathematics Department. UAM Iztapalapa.

More information

A proof for the full Fourier series on [ π, π] is given here.

A proof for the full Fourier series on [ π, π] is given here. niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at

More information

Math Ordinary Differential Equations

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

More information

Chapter 3. Local behavior of nonlinear systems. 3.1 Principle of linearized stability Linearized stability of fixed points of maps

Chapter 3. Local behavior of nonlinear systems. 3.1 Principle of linearized stability Linearized stability of fixed points of maps Chapter 3 Local behavior of nonlinear systems In Chapter 2 we explained how, for linear dynamical systems, the qualitative as well as the quantitative behaviour can be described in terms of eigenvalues

More information

Analytic invariants associated with a parabolic fixed point in C 2

Analytic invariants associated with a parabolic fixed point in C 2 Ergod. Th. & Dynam. Sys. (2008), 28, 1815 1848 doi:10.1017/s0143385707001046 c 2008 Cambridge University Press Printed in the United Kingdom Analytic invariants associated with a parabolic fixed point

More information

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct

More information

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE

CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE CLASSIFICATIONS OF THE FLOWS OF LINEAR ODE PETER ROBICHEAUX Abstract. The goal of this paper is to examine characterizations of linear differential equations. We define the flow of an equation and examine

More information

Hyperbolic PDEs. Chapter 6

Hyperbolic PDEs. Chapter 6 Chapter 6 Hyperbolic PDEs In this chapter we will prove existence, uniqueness, and continuous dependence of solutions to hyperbolic PDEs in a variety of domains. To get a feel for what we might expect,

More information

6 Linear Equation. 6.1 Equation with constant coefficients

6 Linear Equation. 6.1 Equation with constant coefficients 6 Linear Equation 6.1 Equation with constant coefficients Consider the equation ẋ = Ax, x R n. This equating has n independent solutions. If the eigenvalues are distinct then the solutions are c k e λ

More information

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

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

More information

Continuation and bifurcation analysis of delay differential equations

Continuation and bifurcation analysis of delay differential equations Continuation and bifurcation analysis of delay differential equations Dirk Roose 1 and Robert Szalai 2 1 Department of Computer Science, Katholieke Universiteit Leuven, Belgium 2 Department of Engineering

More information

MATH 205C: STATIONARY PHASE LEMMA

MATH 205C: STATIONARY PHASE LEMMA MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)

More information

from delay differential equations to ordinary differential equations (through partial differential equations)

from delay differential equations to ordinary differential equations (through partial differential equations) from delay differential equations to ordinary differential equations (through partial differential equations) dynamical systems and applications @ BCAM - Bilbao Dimitri Breda Department of Mathematics

More information

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

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

More information

x x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b)

x x2 2 + x3 3 x4 3. Use the divided-difference method to find a polynomial of least degree that fits the values shown: (b) Numerical Methods - PROBLEMS. The Taylor series, about the origin, for log( + x) is x x2 2 + x3 3 x4 4 + Find an upper bound on the magnitude of the truncation error on the interval x.5 when log( + x)

More information

CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION

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

More information

Lecturer: Bengt E W Nilsson

Lecturer: Bengt E W Nilsson 2009 05 07 Lecturer: Bengt E W Nilsson From the previous lecture: Example 3 Figure 1. Some x µ s will have ND or DN boundary condition half integer mode expansions! Recall also: Half integer mode expansions

More information

Sec. 1.1: Basics of Vectors

Sec. 1.1: Basics of Vectors Sec. 1.1: Basics of Vectors Notation for Euclidean space R n : all points (x 1, x 2,..., x n ) in n-dimensional space. Examples: 1. R 1 : all points on the real number line. 2. R 2 : all points (x 1, x

More information

Lecture Notes on Numerical Analysis of Nonlinear Equations

Lecture Notes on Numerical Analysis of Nonlinear Equations Lecture Notes on Numerical Analysis of Nonlinear Equations Eusebius J Doedel Department of Computer Science, Concordia University, Montreal, Canada Numerical integrators can provide valuable insight into

More information

Comparison of the Weyl integration formula and the equivariant localization formula for a maximal torus of Sp(1)

Comparison of the Weyl integration formula and the equivariant localization formula for a maximal torus of Sp(1) Comparison of the Weyl integration formula and the equivariant localization formula for a maximal torus of Sp() Jeffrey D. Carlson January 3, 4 Introduction In this note, we aim to prove the Weyl integration

More information

Numerical Solutions to Partial Differential Equations

Numerical Solutions to Partial Differential Equations Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear

More information

Coherent structures near the boundary between excitable and oscillatory media

Coherent structures near the boundary between excitable and oscillatory media Coherent structures near the boundary between excitable and oscillatory media Jeremy Bellay University of Minnesota Department of Computer Science 200 Union St. S.E. Minneapolis, MN 55455, USA Arnd Scheel

More information

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

Defining Equations for Bifurcations and Singularities

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

More information

Course Summary Math 211

Course Summary Math 211 Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.

More information

First order Partial Differential equations

First order Partial Differential equations First order Partial Differential equations 0.1 Introduction Definition 0.1.1 A Partial Deferential equation is called linear if the dependent variable and all its derivatives have degree one and not multiple

More information

Hyperbolic Systems of Conservation Laws

Hyperbolic Systems of Conservation Laws Hyperbolic Systems of Conservation Laws III - Uniqueness and continuous dependence and viscous approximations Alberto Bressan Mathematics Department, Penn State University http://www.math.psu.edu/bressan/

More information

Ahmed Abbas Mizeal and Mudhir A. Abdul Hussain

Ahmed Abbas Mizeal and Mudhir A. Abdul Hussain ARCHIVUM MATHEMATICUM (BRNO) Tomus 8 (01) 7 37 TWO-MODE BIFURCATION IN SOLUTION OF A PERTURBED NONLINEAR FOURTH ORDER DIFFERENTIAL EQUATION Ahmed Abbas Mizeal and Mudhir A. Abdul Hussain Abstract. In this

More information

Applied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.

Applied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm. Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the

More information

Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations

Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology

More information

Matcont Tutorial: ODE GUI version

Matcont Tutorial: ODE GUI version Matcont Tutorial: ODE GUI version Hil Meijer Nottingham, Oct-Nov, 2011 If you want to get credit for solving a complicated mathematical problem, you will have to provide a full proof. But if you re trying

More information

Allen Cahn Equation in Two Spatial Dimension

Allen Cahn Equation in Two Spatial Dimension Allen Cahn Equation in Two Spatial Dimension Yoichiro Mori April 25, 216 Consider the Allen Cahn equation in two spatial dimension: ɛ u = ɛ2 u + fu) 1) where ɛ > is a small parameter and fu) is of cubic

More information

OPTIMAL CONTROL OF A CLASS OF PARABOLIC

OPTIMAL CONTROL OF A CLASS OF PARABOLIC OPTIMAL CONTROL OF A CLASS OF PARABOLIC TIME-VARYING PDES James Ng a, Ilyasse Aksikas a,b, Stevan Dubljevic a a Department of Chemical & Materials Engineering, University of Alberta, Canada b Department

More information

Direct and indirect methods for optimal control problems and applications in engineering

Direct and indirect methods for optimal control problems and applications in engineering Direct and indirect methods for optimal control problems and applications in engineering Matthias Gerdts Computational Optimisation Group School of Mathematics The University of Birmingham gerdtsm@maths.bham.ac.uk

More information

Chapter 3 Second Order Linear Equations

Chapter 3 Second Order Linear Equations Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,

More information

Global Bifurcation Diagram for the Kerner-Konhäuser Traffic Flow Model

Global Bifurcation Diagram for the Kerner-Konhäuser Traffic Flow Model arxiv:1311.4119v1 [math.ds] 17 Nov 2013 Global Bifurcation Diagram for the Kerner-Konhäuser Traffic Flow Model Joaquín Delgado Patricia Saavedra Keywords: Continuous traffic flow. Traveling waves. Bautin

More information

Sparse representations from moments

Sparse representations from moments Sparse representations from moments Bernard Mourrain Inria Méditerranée, Sophia Antipolis BernardMourrain@inriafr Sparse representation of signals Given a function or signal f (t): decompose it as f (t)

More information

Scientific Computing I

Scientific Computing I Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to

More information

Stability of flow past a confined cylinder

Stability of flow past a confined cylinder Stability of flow past a confined cylinder Andrew Cliffe and Simon Tavener University of Nottingham and Colorado State University Stability of flow past a confined cylinder p. 1/60 Flow past a cylinder

More information