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

Transcription

1 FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. Vol. 22, No. 1 (27), pp SUPERCRITICAL HOPF BIFURCATION IN DELAY DIFFERENTIAL EQUATIONS AN ELEMENTARY PROOF OF EXCHANGE OF STABILITY Hamad Talibi Alaoui and Radouane Yafia Abstract. The local Hopf Bifurcation theorem states that under supercritical conditions, stability passes from the trivial branch to the bifurcated one. We give here an elementary proof of this result, based on the following steps: i) reduction to a two-dimensional system via the center manifold technic; ii) estimation of the distance between solutions of the original and the reduced equations. Thus there is no reference to the Floquet theory. Incidently, we obtain an estimate of the stability region. 1. Generalities and the Reduced System We first recall the Hopf Bifurcation theorem for retarded functional differential equations. Consider the system (1.1) u(t) = Lu t + R(α, u t ), with the following hypotheses: (H ) L : C R n is a linear continuous operator, R : R C R n is C k+1, k 1, such that R(α, ) =, for all α R and D 2 R(α, ) =, for some α R; where C = C([ r, ], R n ), n N, u t (θ) = u(t + θ) for θ [ r, ]. (H 1 ) The characteristic equation (1.2) (λ, α) =, with (λ, α) = det (λid (L + D 2 R(α, ))e λ( ) Id) of the linearized equation of equation (1.1) at the equilibrium point, has only roots with a negative real parts for α α, and exactly two roots cross Received November 3, Mathematics Subject Classification. 34Kxx, 34K18. 21

2 22 H. Talibi Alaoui and R. Yafia the imaginary axis at α = α, at points iν and iν, with multiplicity 1, ν >. More precisely, if λ(α) is the characteristic root of (1.2), bifurcated from iν, we assume that: d (H 2 ) dα Re λ(α) α=α >. Denote by (T α (t)) t the semi group associated with the linearized equation of (1.1) and by A α it s infinitesimal generator. To the continuous branch of eigenvalues λ(α), we can associate a continuous branch of eigenvectors a(α) + ib(α). Let N α = span{a(α), b(α)} and Φ α = (a(α), b(α)). There exists a 2 2 matrix B(α) of the following form ( r(α) β(α) B(α) = β(α) r(α) where λ(α) = r(α) + iβ(α) such that A α Φ α = Φ α B(α). Using the formal adjoint theory for FDEs in [5], N α = {Φ α Ψ α, ϕ, ϕ C} and C can be decomposed as C = N α S α, where Ψ α = col(e(α), f(α)) is a basis of the adjoint space N α, such that Ψ α, Φ α = I (the 2 2 identity matrix) and, is the bilinear form on C C associated with the adjoint equation of the linearized equation of (1.1) around. Using the decomposition u t = Φ α x(t)+ y t, x(t) R 2, y t S α in the variation of constants formula introduced in [5], equation (1.1) reads as (1.3) x(t) = e B(α)t x + y t = T α (t)y + ), e B(α)(t τ) Ψ α ()F (α, Φ α x(τ), y τ )dτ, d[k α (t, τ) S ]F (α, Φ α x(τ), y τ ), where F (α, ϕ, ψ) = R(α, ϕ + ψ) D 2 R(α, ϕ + ψ), K α (t, τ) S = K α (t, τ) Φ α Ψ α, K α (t, τ), K α (t, τ) = τ X α(t+θ ν)dν, X α is the n n fundamental matrix solution of the linearized equation of (1.1). Looking for solutions defined on the whole real axis and uniformly bounded on their domain, the second identity in (1.3) can be written starting from any initial point σ as y t = T α (t σ)y σ + σ d[k α(t, τ) S ]F (α, Φ α x(τ), y τ ). Letting σ, the terms T α (t σ)y(σ) approach to, while the integral has a limit, which yields an expression where the projection y has

3 Supercritical Hopf Bifurcation in Delay Differential Equations been eliminated and the system (1.3) has the form: (1.4) x(t) = e B(α)t x + y t = e B(α)(t τ) Ψ α ()F (α, Φ α x(τ), y τ )dτ, d[k α (t, τ) S ]F (α, Φ α x(τ), y τ ). In view of the center manifold theorem (see [5]), there exists a map g(α, ) defined from a neighborhood of in N α into a neighborhood of in S α such that y t = g(α, Φ α x(t)) is the equation of a local center manifold, for α close to α, where g is C k+1, g(α, ) = and D 2 g(α, ) =. Then the system (1.4) on the center manifold can be reduced to the following problem (1.5) x(t) = e B(α)t x + e B(α)(t τ) Ψ α ()F (α, Φ α x(τ), g(α, Φ α x(τ))dτ, through this transformation: (x(t), y t ) is the solution of (1.4) on the center manifold x(t) is the solution of (1.5) and y t = g(α, Φ α x(t)). In polar coordinates x = (ρ cos φ, ρ sin φ), system (1.5) reads as { ρ = r(α)ρ + h1 (α, ρ, φ), (1.6) φ = β(α) + h 2 (α, ρ, φ), where h 1 (α, ρ, φ) and h 2 (α, ρ, φ) have the form: h 1 (α, ρ, φ) = cos φf 1 (α, ρ cos φ, ρ sin φ) + sin φf 2 (α, ρ cos φ, ρ sin φ), h 2 (α, ρ, φ) = 1 ρ [ sin φf 1(α, ρ cos φ, ρ sin φ) + cos φf 2 (α, ρ cos φ, ρ sin φ)], f 1, f 2 are the components of Ψ α ()F (α, Φ α x(t), g(α, Φ α x(t))) in R 2. For ρ near to, we have φ ν t (in view of the second equation in (1.6)). Therefore t is diffeomorphic to φ, we can eliminate t between the two equations in (1.6). We obtain (1.7) dρ dφ = r(α)ρ + h 1(α, ρ, φ) β(α) + h 2 (α, ρ, φ). Let ρ(φ, c) be the solution of the last equation with initial data c.

4 24 H. Talibi Alaoui and R. Yafia The periodic solution with initial value c correspond to finding c and α such that (1.8) G(α, c) =, where G is the Hopf bifurcation function (1.9) G(α, c) = 2π r(α)ρ(s, c) + h 1 (α, ρ(s, c), s) ds. β(α) + h 2 (α, ρ(s, c), s) From the Hopf bifurcation theorem (see [2], [3] and [4]), there exists ε >, and functions P (ε), ω(ε) and α(ε) (ε [, ε )) sufficiently regular such that P (ε) is ω(ε) periodic solution of equation (1.5) for α = α(ε). Furthermore, ω() := ω = 2π/ν, α() = α, α () = and the bifurcation is subcritical if α () <, and supercritical if α () >. Now, we introduce the assumption on supercritical Hopf bifurcation: (H 3 ) α () >. Remark 1.1. a) (H 3 ) means that the bifurcated branch lies in the parameter region where the trivial solution is unstable, so, we will find out, as is well known ([7]), that the system stabilizes around the nontrivial periodic solutions. b) (H 3 ) implies that α is increasing along the bifurcated branch, occasionally, we will use α instead of ε as an independent variable along the branch. Information about elements of bifurcation (α, ω, P ) are needed here. In fact, the following classical features will be enough for our purposes (see [8] and [9]): α = α + ε 2 α 2 + ε 4 α 4 + (1.1) ω = ω + ε 2 ω 2 + ε 4 ω 4 + P (ε)(t) = ε P (t). Note that the elements of Hopf bifurcation in (1.1) are obtained formally by inserting Taylor expansions of α, ω, ρ: α = α(ε) = α i ε i, ω = ω(ε) = ω i ε i, ρ = ε( ρ i ε i ) into the equation (1.8) and the integrated form of (1.6), with ρ() = ε, φ() =, ρ(ω) = ε, φ(ω) = 2π, and equating like powers of ε.

5 Supercritical Hopf Bifurcation in Delay Differential Equations Stability, Asymptotic Stability Along the Branch of Bifurcation From now on, we assume all of the hypotheses introduced before: (H ) through (H 3 ). Our first step is to give (1.4) in another form, by centering it around y t = g(α, Φ α x(t)). We set y t = z t + g(α, Φ α x(t)), that is, we consider (x, z) instead of (x, y). Inserting this into system (1.4), we obtain: x(t) = e B(α)t x + e B(α)(t τ) R 1 (α, x(τ), z τ )dτ, (2.1) z t = d[k α (t, τ) S ]R 2 (α, x(τ), z τ ), where (2.2) and R 1 (α, x(t), z t ) = Ψ α ()F (α, Φ α x(t), z t + g(α, Φ α x(t))), R 2 (α, x(t), z t ) = F (α, Φ α x(t), z t +g(α, Φ α x(t))) F (α, Φ α x(t), g(α, Φ α x(t))). Note that the bifurcating periodic solutions of (2.1) lie on z =, that is, are in fact the same as bifurcating solutions of (1.5). Extending the notations of Section 1, we denote such solutions by p = (P, ). We will normalize these solutions, considering that P () = P = (ε, ), ε >. Let (x, z ) be a data near (P, ), with x = (c, ), c >. We will denote by x the solution of (1.5) such that x () = x, and by (x #, z) the solution of (2.1) such that x # () = x, z t= = z, and let ω = ω (α, x ) denotes the first return-time for x, that is x (ω ) = (ζ, ), ζ > and ω > is the solution of the equation (2.3) φ(ω) = 2π; ω # = ω # (α, x, z ) denotes the first return-time for x #, that is the solution of the analogous equation of (2.3) resulting from polar coordinates in the first equation in (2.1). The main result of the paper is the following one. Theorem 2.1. Under hypotheses (H ) through (H 3 ), there exists ε 1 >, η >, K > such that for each ε, < ε < ε 1 ; θ, θ 1, the relations P = (ε, ), ε c θε 2 and z ηθε 3 imply i) x # (ω # ) P θ(ε 2 Kε 4 ); ii) z ω # ηθε 3.

6 26 H. Talibi Alaoui and R. Yafia Corollary 2.1. Assume (H ) through (H 3 ). Then the bifurcated orbits of equation (1.1) for α close enough to α are asymptotically stable. Proof. In view of the transformations considered above, we only have to prove the same fact for the equation (2.1). Starting with (x, z ) and θ such that ε c θε 2, z ηθε 3, denoting by (c 1, ) = x # (ω # ), z 1 = z ω # and cj, z j the corresponding terms after j rotations (c j, ) = x # (ω # j ) and φ# (ω # j ) = 2jπ, we obtain: ε c j θ(1 Kε 2 ) j ε 2, ηθε 3. z j Thus (x(t), z t ) approaches exponentially the orbit of p. Theorem 2.1 is a consequence of two intermediate results. We will first state these as lemmas, and derive the Theorem from them. Lemma 2.1. Assume the hypotheses of Theorem 2.1. Then there exists K 1 > such that x (ω ) P θ(ε 2 K 1 ε 4 ). Remark 2.1. Lemma 2.1 states the stability result for a two-dimensional system. There no doubt however that the estimation in Lemma 2.1 is the control estimate for Theorem 2.1. Precisely, our purpose here is to connect in an elementary way the result for such systems to the corresponding stability result for higher dimensional systems. Lemma 2.2. Assume the hypotheses of Theorem 2.1 and fix Then there exists K 2 >, such that: ω 1 > max(ω, ω, ω # ). i) x # (t) P (t) θk2 ε 2, t ω 1 ; ii) z t θε 3 K 2 η, t ω 1 ; iii) x # (t) x (t) θk2 ε 4 η, t ω 1. Proof of Lemma 2.1. In polar coordinates, let (ρ, φ) corresponding to the ω periodic solution P (of (1.6)) guaranteed by Hopf Bifurcation theorem and (ρ, φ ) corresponding to x with ρ() = ε and ρ () = c. From (1.9) we obtain x (ω ) P = ρ (2π) ε = c ε + G(α, c).

7 Supercritical Hopf Bifurcation in Delay Differential Equations On one hand, G(α, c) = G(α, ε) + (c ε)d c G(α, ε) (c ε)2 D 2 ccg(α, ε + ν(c ε)), for some ν [, 1]. In view of (1.8), G(α, ε) =. By differentiation with respect to ε (with α = α(ε)), we obtain D c G(α, ε) = α (ε)d α G(α, ε) = εα (ε)d α G (α, ε), where G (α, ε) = 1 G(α, ε), ε for ε, D c G(α, ), for ε =. Setting K(c, ε) = α (ε) D α G (α, ε) 1 ε 2ε 2 (c ε)d2 ccg(α, ε + ν(c ε)), we have G(α, c) = (c ε)ε 2 K(c, ε). Our next task is to look at limk(c, ε). c ε Claim 2.1. lim c ε K(c, ε) = α ()2πr (α ) ω >. Proof. In view of (1.1) and the hypotheses of Theorem 2.1 we have α (ε) = O(ε) and c ε = O(ε 2 ). Moreover, by (1.9), we have: α (ε) D α G (α, ε) α () 2πr (α ) as ε ε ω (because for c = the corresponding solution ρ = and h i = o(ρ), i = 1, 2). Using the Taylor expansion of h 1 and h 2 in terms of ρ, we have h 1 (α, ρ, s) = ρ 2 C 3 (α, s) + ρ 3 C 4 (α, s) + h 2 (α, ρ, s) = ρd 3 (α, s) + ρ 2 D 4 (α, s) + where C j and D j are polynomials of degree j with respect to cos(s), sin(s).

8 28 H. Talibi Alaoui and R. Yafia Substituting these expansions for h 1 and h 2 in the integral expression of (1.7) we obtain (2.4) ρ (φ) = c + 1 β(α) φ where for each j, D j = D j/β(α). In view of (2.4), we have D 2 ccg(α, c) = 1 β(α) [ r(α)ρ + ( C 3 D 3r(α) ) (ρ ) 2 + ( C 4 C 3 D 3 + r(α) ( ] (D 3) 2 D 4)) (ρ ) 3 + ds, 2π [ r(α)d 2 ccρ + ( C 3 D 3r(α) ) D 2 cc(ρ ) 2 ( + C 4 C 3 D 3 + r(α) ( (D 3) 2 D 4) ) ] Dcc(ρ 2 ) 3 + ds. From (2.4) once again, we can deduce that limd c ρ = 1; c ε limdcc(ρ 2 ) 2 = 2, c ε limdcc(ρ 2 ) j = for j 3 c ε (since r(α ) = and limρ = ). c ε Therefore, we may conclude that limdccg(α, 2 c) = c ε = 1 2π β(α ) 2 β(α ) 2π C 3 (α, s)dcc(ρ 2 ) 2 ε= ds c= C 3 (α, s)ds = (since C 3 is a polynomial in (cos(s), sin(s)) of degree 3). Coming back to K, this yields the result. By continuity property of K, we may deduce that ( ) K(c, ε) K 1 >, for (c, ε) in a neighborhood of (, ) and some K 1 >. We may conclude the proof of Lemma 2.1: taking ε 1 small enough, we can assert that for each pair (c, ε) satisfying the requirements stated in Theorem 2.1, condition ( ) holds.

9 Supercritical Hopf Bifurcation in Delay Differential Equations Therefore, we have x (ω ) P c ε (1 ε 2 K 1 ). This complete the proof of Lemma 2.1. Proof of Lemma 2.2. According to (1.5) and the first equation of (2.1), we have (2.5) x # (t) P (t) = e B(α)t (x P ) [ ] + e B(α)(t τ) R 1 (α, x # (τ), z τ ) R 1 (α, P (τ), ) dτ, (2.6) R 1 (α, x # (τ), z τ ) R 1 (α, P (τ), ) R 1 (α, x # (τ), z τ ) R 1 (α, x # (τ), ) + R 1 (α, x # (τ), ) R 1 (α, P (τ), ), R 1 (α, x # (τ), z τ ) R 1 (α, x # (τ), ) D 3 R 1 (α, x # (2.7) (τ), µz τ ) z τ, and (2.8) R 1 (α, x # (τ), ) R 1 (α, P (τ), ) D 2 R 1 (α, P (τ) + ν(x # (τ) P (τ)), ) x # (τ) P (τ) for some µ, ν [, 1]. From the hypotheses of Theorem 2.1, if ε then (x, z ) (, ) and lim ε x# (t) =, limz t =. Therefore, x # (t) = o(1), z t = o(1) as ε and ε (2.9) { D3 R 1 (α, x # (t), µz t ) = o(1), D 2 R 1 (α, P (t) + ν(x # (t) P (t)), ) = o(1). According to (2.5) through (2.9), with t ω 1, we obtain } x # (t) P (t) M 1 { x P +o(1) max τ t x# (τ) P (τ) +o(1) max z τ τ t with M 1 = max t ω 1 e B(α)t. We conclude that (2.1) x # (t) P (t) M { } x P + o(1) max z τ, τ t

10 3 H. Talibi Alaoui and R. Yafia where M > is a constant. The second equation in (2.1) gives z t T α (t)z + R 2 (α, x # (τ), z τ ) = d[k α (t, τ) S ]R 2 (α, x # (τ), z τ ). R 2 (α, x # (τ), z τ ) R 2 (α, x # (τ), ) D 3 R 2 (α, x # (τ), µz τ ) z τ. } Then z t N 1 { z + o(1) max z τ, where N 1 = max T α (t), and τ t t ω 1 we deduce (2.11) z t N z. With the estimate given in Theorem 2.1, (2.1) and (2.11) lead to i) and ii) of Lemma 2.2. Now, from (1.5) and the first equation in (2.1), we obtain: x # (t) x (2.12) (t) e B(α)(t τ) (R 1 (α, x # (τ), z τ ) R 1 (α, x (τ), )) dτ. R 1 (α, x # (τ), z τ ) R 1 (α, x (τ), ) D 3 R 1 (α, x # (τ), µz τ ) z τ + D 2 R 1 (α, x (τ) + ν(x # (τ) x (τ)), ) x # (τ) x (τ) for some µ, ν [, 1]. Finally, from (2.9), there exists K > such that R 1 (α, x # (τ), z τ ) R 1 (α, x (τ), ) [ K ε max x # (τ) x (τ) τ t ] z τ τ t + max with a constant K >, and from (2.12), we obtain x # (t) x (t) εk max z τ (with K > ), τ t

11 Supercritical Hopf Bifurcation in Delay Differential Equations which, together with ii), gives iii). Proof of Theorem 2.1. i) In inequality (2.13) x # (ω # ) P x # (ω # ) x (ω # ) + x (ω # ) x (ω ) + x (ω ) P the first and the third addend on the right hand side can be estimated using Lemma 2.1 and Lemma 2.2, respectively. In order to estimate the second addend, we need to evaluate ω # ω. Claim 2.2. ω ω # ηθk ε 3. Proof. Near the bifurcation point, equation (1.5) varies more slowly in magnitude than in direction. This is a classical observation, best perceived d in polar coordinates representation: dt φ β(α d ); dtρ = o(ρ). So, for some ε >, C, C, C > we have: C x x (t) C x, d dt φ (2.14) C, for t ω 1, ε ε. Define ψ = angle(x # (ω # ), x (ω # )) = angle(x (ω ), x (ω # )) = φ (ω # ). Using Lemma 2.2 iii) and (2.14), we have tan ψ = tan φ (ω # ) x (ω # ) x # (ω # ) x (ω ) θk 2ε 3 η C. On the other hand tan ψ φ (ω # ) C ω ω #. This yields the desired estimate: ω ω # ηθk ε 3. The second addend of (2.13) can now be estimated using x (ω # ) x (ω ) sup dx (t) dt ω ω #. <t<ω 1 From equation (2.2) we have R 1 (α, x, ) = o(x), and from the proof of Claim 2.2 we have, dx (t) dt C C ε, for some constant C. Thus, in view of the claim 2.2, we obtain x (ω # ) x (ω ) ηθε 4 K 4.

12 32 H. Talibi Alaoui and R. Yafia To reach the conclusion of i) in Theorem 2.1, choose η small enough for K 2 η < K 1 /4 (Lemma 2.1 and Lemma 2.2) and K 4 η < K 1 /4; define K = K 1 /2, and replace addend of right hand side of (2.13) by their corresponding estimates. From Lemma 2.2 we conclude ii). R E F E R E N C E S 1. A. Casal and M. Freedman: Linsted-Poincaré approach to bifurcation problems for differential-delay equations. IEEE Trans. Automatic Control 25 (198). 2. S. N. Chow and J. K. Hale: Methods of Bifurcation Theory. Springer- Verlag New-York, Heidelberg, Berlin, S. N. Chow and J. Mallet-Paret: Integral averaging and bifurcation. J. Differential Equations 26 (1977), T. Faria and L. T. Magalhas: Normal form for retarded functional differential equations with parameters and applications to hopf singularity. J. Differential Equations 122 (1995), J. K. Hale and S. M. Verduyn Lunel: Introduction to Functional Differential Equations. Springer-Verlag, New York, B. D. Hassard, N. D. Kazarinoff and Y-H. Wan: Theory and Applications of Hopf Bifurcation. Cambridge University Press, J. E. Marsden and M. McCracken: The Hopf Bifurcation and its Applications. Applied Math. Series, Vol. 19, Springer-Verlag, New York, H. Talibi and Casal: A perturbation approach to bifurcation branches of periodic solutions for some functional differential equations. Facta Univ. Ser. Math. Inform. 4 (1989), H. Talibi: Contribution à l Etude du Problème de Bifurcation de Hopf dans le Cadre des Equations Differentielles à Retard. Thèse d état, Universit Mohammed V-Rabat, Université Chouaib Doukkali Faculté des Sciences Département de Mathématiques et Informatique B.P: 2, El Jadida, Morocco talibi@math.net Université Ibn Zohr Faculté Polydisciplinaire B.P: 638, Ouarzazate, Morocco yafia1@yahoo.fr