Tibor Krisztin. Gabriella Vas
|
|
- Teresa Oliver
- 5 years ago
- Views:
Transcription
1 Electronic Journal of Qualitative Theory of Differential Equations 11, No. 36, 1-8; ON THE FUNDAMENTAL SOLUTION OF LINEAR DELAY DIFFERENTIAL EQUATIONS WITH MULTIPLE DELAYS Tibor Krisztin Bolyai Institute, University of Szeged, Hungary Gabriella Vas Analysis and Stochastic Research Group of the Hungarian Academy of Sciences Bolyai Institute, University of Szeged, Hungary Abstract. For a class of linear autonomous delay differential equations with parameter α we give upper bounds for the integral X t, α) dt of the fundamental solution X, α). The asymptotic estimations are sharp at a critical value α where x = loses stability. We use these results to study the stability properties of perturbed equations. Key words and phrases: linear delay differential equation, fundamental solution, Laplace transform, discrete Lyapunov functional Mathematics Subject Classification: 34K6, 34K 1. Introduction Let n 1 be an integer, r 1,...,r n be real numbers with r 1 > r >... > r n, a 1,...,a n be positive numbers, α >. In case n = 1, we assume r 1 >. Consider the linear delay differential equation 1.1) ẋt) = α a j xt r j ) with a real parameter α, α ]. The natural phase space for Eq. 1.1) is C = C [ r 1, ],R), the space of all real valued continuous functions defined on [ r 1, ] equipped with the supremum norm. For each ϕ C, there exists a unique solution x ϕ : [ r 1, ) R with x ϕ t) = ϕt), r 1 t. EJQTDE, 11 No. 36, p. 1
2 The characteristic function for 1.1) is z, α) = z + α a j e rjz, z C, < α α. The zeros of the characteristic function are the eigenvalues of the generator of the strongly continuous semigroup defined by the solution operators Tt) : C ϕ x ϕ t C, t, where the segment x ϕ t C, t, is defined by x ϕ t s) = x ϕ t + s), r 1 s. We assume that the following hypothesis holds throughout the paper. H1): For α = α, there exists a unique pair of purely imaginary and simple eigenvalues iν, iν with ν >. There exist α 1, α ) and γ < such that for all α [α 1, α ), there is a unique complex conjugate pair of simple eigenvalues λ = λα), λ = λα) in {z C : γ < Rez < } with lim Reλα) = and lim Imλα) = ν. α α α α For α [α 1, α ], all the other eigenvalues are found in {z C : Rez < γ}. For each α, the fundamental solution of Eq. 1.1) is the function X, α) : [ r 1, ) R with initial condition if r 1 t < 1.) X t, α) = 1 if t =, that satisfies 1.3) Xt, α) = 1 α t a jˆ X s r j, α) ds for all t. It is clear that X exists uniquely, X [, ) is continuous and X r1, ) is continuously differentiable. It is well known [9] that X t, α) dt < provided Rez < for all zeros of z, α). Our aim is to give explicit estimations for X t, α) dt. The integral X t, α) dt has a role via the variation-of-constants formula in perturbation results. For example Eq. 1.1) can appear as the linear variational equation at a stationary point of a nonlinear delay differential equation. If the solution x = of Eq. 1.1) is asymptotically stable, then the stationary point of the original nonlinear equation is locally attracting. Integral X t, α) dt plays an important role in the estimation of the attractivity region of the stationary point [7, 9]. EJQTDE, 11 No. 36, p.
3 The technique applied to estimate X t, α) dt contains a splitting of the spectrum by the vertical line Rez = γ < so that there is no eigenvalue on Rez = γ. Then the phase space can be decomposed as C = P Q, where P is the realified generalized eigenspace of the generator corresponding to the spectrum in Rez > γ, and Q is the realified generalized eigenspace corresponding to the spectrum in Rez < γ. The solution operator T t) is easily estimated on P as it is finite dimensional. On Q it is well known that T t)ϕ M γ)e γt ϕ holds for all ϕ Q and t with some constant M γ) 1. An explicit upper bound for M γ) is crucial in our estimation for X t, α) dt. Giving an optimal upper bound for M γ) is also interesting in the construction of invariant manifolds, in particular when the size of the manifolds is of key importance. E.g., in order to prove that the local attractivity of implies global attractivity for Wright s equation, the estimates for M γ) of this paper are used to find bounds for the size of a center manifold [11]. Although the estimates for X t, α) dt seem to be a fundamental technical issue, as far as we know, not much is known except for the results of Győri and Hartung in [5, 6, 7]. In the single delay case n = 1 their estimates are sharp for small values of α, but not for α close to the critical value α. This paper is organized as follows. We present the results in Section. Section 3 estimates the location of the leading pair of eigenvalues. Sections 4 contains the proofs of the first three theorems. For the single delay case Theorem.4) a different proof is given in Section 5 yielding a sharper result. An example is shown in Section 6 with two delays. Section 7 presents two applications of the results for perturbed equations. Note that. Main results z λ, α) = n z λ, α) = 1 α a j r j e rjλ, and for λ = µ + iν, z λ, α) = 1 α a j r j e rjµ cosr j ν)) + α a j r j e rjµ sin r j ν)). According to the results [5, 6] of Győri and Hartung, ˆ ) 1 X t, α) dt = O α α). EJQTDE, 11 No. 36, p. 3
4 Our first theorem gives a sharp asymptotic estimation for X t, α) dt as α α implying that X t, α) dt = O α α) 1). Theorem.1. Under hypothesis H1) lim α α) α α 4 z X t, α) dt = iν, α ) n ) πα a n ). jr j sinr j ν ) a j sin r j ν ) From the application point of view it is more important to give an explicit upper bound for X t, α) dt on interval [α 1, α ). This is contained in the next result in terms of α, λ = λα), µ = µ α) = Reλα) and ν = ν α) = Imλα) guaranteed by H1). Theorem.. Under hypothesis H1) X t, α) dt for all α 1 α < α with.1) L 1 α) = 1 πr 1 and νe µ ν π 1 + e µ π) ν z λ, α) µ + ν ) 1 e µ π) + max {L 1 α), L α)} ν γ 1 + ω ) ˆ 8 ω + ω n ω = α a j e rjγ,.) L α) = e γr1 e αr1 P n aj + ) dω γ + iω, α) <, z ) <. λ, α) We are going to check that the upper bound given by Theorem. is sharp in the sense that this upper estimate multiplied by α α) has the same limit at α as function α α) X t, α) dt. There is a need for easily computable upper bounds. Our next aim is to give an estimate that is independent of λ = λα). For simplicity, set c 1 = 1 + α n a jr j. For each δ, π/ r 1 )), set c = c δ) = a 1 r 1 sinr 1 δ) > and fix K = K δ) > so large that K > c 1 e r1γ c δ holds. We will need the following additional hypothesis besides H1). H): There exists δ, π/ r 1 )) such that for all α [α 1, α ], we have δ Imλ π/r 1 δ and c < 1e r1γ c δ + c 1 α 1 K δ). 1 Reλ EJQTDE, 11 No. 36, p. 4
5 The upper bound given by the next result is not sharp, but it is independent of λ = λα). Theorem.3. If H1) and H) hold, then α α) X t, α) dt K δ) γe γ δ π 1 + e γ δ π ) α 1 c δ 1 e γ δ π) max { L1, L } α α 1 ) + γ for all α 1 α < α, where L 1 = sup α1 α α L 1 α), L 1 α) is defined by.1) and P n L = e γr1 e αr1 aj + ). α 1 c The particular case n = 1, r 1 = 1, a 1 = 1 is of special interest as equation.3) ẋ t) = αxt 1) is the simplest delay differential equation, and it appears as a linearization of famous equations of the form ẋt) = f xt 1)). However, surprisingly, little is known about X t, α) dt even for this simple case when α α. Theorem.1 is a generalization of a result of Krisztin in [1] saying that for this equation π ) lim α π α X t, α) dt = π π 1.5. The result of Theorem.3 can be substantially improved for.3). This is essential in the estimation of the attractivity region of x = for the Wright s equation for α near the critical value π/, see []. Theorem.4. If X, α) : [ r 1, ) R is the fundamental solution of Eq..3), then π α ) π ) X t, α) dt α for α [ 3, π ). We remark that the upper bound given in [5, 7] for the integral of the fundamental solution of Eq..3) is sharp only for small α >, in particular for α, e 1]. 3. The real part of the leading eigenvalues It is of key importance to understand the behavior of Reλα) near the critical value α. EJQTDE, 11 No. 36, p. 5
6 For all z = u + iv C and α α, set g u, v, α) = Re u + iv, α) = u + α h u, v, α) = Im u + iv, α) = v α a j e rju cosr j v), a j e rju sin r j v). Then g and h are smooth functions with the following partial derivatives: 3.1) 3.) n h u, v, α) = u v u, v, α) = 1 α a j r j e rju cosr j v), n h u, v, α) = v u u, v, α) = α a j r j e rju sin r j v), 3.3) n α u, v, α) = a j e rju cosr j v), 3.4) Note that 3.5) z z, α) = 1 α h n α u, v, α) = a j e rju sin r j v). a j r j e rjz = u, v, α) i u, v, α). u v If µ = µ α) = Reλα) and ν = ν α) = Imλα), where λ is the leading eigenvalue in H1), then g µ, ν, α) = and h µ, ν, α) = for all α 1 α α. In particular, g, ν, α ) = and h, ν, α ) =. By condition H1), λα) is a simple zero of λ, α), that is ) ) 3.6) z λ, α) = µ, ν, α) + µ, ν, α) > u v for all α 1 α α with a lower bound independent of α. The smooth dependence of µ and ν on α is easily guaranteed. Proposition 3.1. Assume that condition H1) holds. Then µ and ν are C 1 -smooth functions of α on [α 1, α ] with 3.7) µ α) = v µ, ν, α) h α µ, ν, α) u µ, ν, α) α µ, ν, α) ) ), α [α 1, α ]. u µ, ν, α) + v µ, ν, α) Proof. Choose α [α 1, α ] arbitrarily. As g µ, ν, α) =, h µ, ν, α) =, and EJQTDE, 11 No. 36, p. 6
7 det u h u µ, ν, α) µ, ν, α) h v v µ, ν, α) µ, ν, α) ) = ) ) µ, ν, α) + µ, ν, α) u v by our initial assumption, the Implicit Function Theorem yields the first assertion. Differentiating the equations with respect to α, we get u µ, ν, α)µ α) + v µ, ν, α)ν α) + µ, ν, α) =, α h u µ, ν, α)µ α) + h v µ, ν, α)ν α) + h µ, ν, α) =, α from which the formula for µ α) easily follows. Corollary 3.. If H1) holds, then n ) α a n ) jr j sin r j v ) a j sin r j v ) lim α α µ α) = z iν, α ). Proof. By hypothesis H1), lim α α λα) = iν. Proposition 3.1 with 3.1)-3.5) and relation α, ν, α ) = gives the statement of the corollary. a j cosr j ν ) = g, ν, α ) α = If ν is bounded away from, and µ is sufficiently small for all α, then we give an upper bound for α α) / µ α). The following corollary is needed in the proof of Theorem.3. Corollary 3.3. Suppose that H1) and H) hold. Then α α µ α) K δ) for each α [α 1, α ). Proof. Proposition 3.1 gives ) µ v µ, ν, α) h α µ, ν, α)) u α) = µ, ν, α) α µ, ν, α) u µ, ν, α) ) + v µ, ν, α) ), α [α 1, α ]. Using this result, we give a positive lower bound for µ α), α [α 1, α ]. It clearly follows from 3.1)-3.) that for all α [α 1, α ], µ, ν, α) 1 + α a j r j )e r1µ c 1 e r1µ, u n u µ, ν, α) α a j r j e r1µ c 1 e r1µ EJQTDE, 11 No. 36, p. 7
8 and µ, ν, α) α v a j r j )e r1µ c 1 e r1µ, with constant c 1 introduced before Theorem.3. Note that if H) holds, then sinr j ν) > for all j {1,..., n} and sin r 1 ν) > sin r 1 δ). Hence v µ, ν, α) αa 1r 1 e r1µ sin r 1 ν) αc e r1µ >, where c = a 1 r 1 sin r 1 δ). Equations 3.3), 3.4) with g µ, ν, α) = and h µ, ν, α) = give that Therefore µ µ, ν, α) = α α Conditions given in H) now yield and µ α) αc e r1µ ν α c r1µ µ 1e α c 1 e r1µ h α µ, ν, α) = ν α. = c ν + c 1 α 1 µ c 1 e r1µ for α [α 1, α ]. µ α) c δ + c 1 α 1 1 µ c 1 e r1γ > 1 K δ). The Lagrange Mean Value Theorem implies that for each α [α 1, α ), µ α) = µ α ) µ α) = µ ξ)α α) with some α < ξ < α. Thus the previous result implies µ α) α α) /K δ) for all α [α 1, α ), and the proof is complete. 4. The Proofs of Theorems.1-.3 Under hypothesis H1), the phase space C = C [ r 1, ],R) can be decomposed as C = P Q into the closed subspaces P and Q, where P is the realified generalized eigenspace of the generator associated with the leading eigenvalues λ = µ + iν, λ = µ iν, and Q is the realified generalized eigenspace associated with the rest of the spectrum of the generator. Subspace P is spanned by e µt cosνt) [ r1,] and e µt sinνt) [ r1,], therefore dimp =. Both P and Q are invariant subspaces for the solution segments of Eq. 1.1) in the sense that if x : [ r 1, ) R is a solution of Eq. 1.1) and x T P Q) for some T, then x t P Q) for all t T. The decomposition C = P Q defines a projection Pr P onto P along Q and a projection Pr Q onto Q along P. EJQTDE, 11 No. 36, p. 8
9 For all α [α 1, α ], set functions p, α) : [ r 1, ) R and q, α) : [ r 1, ) R by pt, α) = z=λ,λ e zt Res z, α), t r 1, and q t, α) = Xt, α) pt, α), t r 1. For simplicity we also use notations p = p, α) and q = q, α). As z=λ,λ e zt Res z, α) = e λt z λ, α) + z e λt e λt ) = Re λ, α z λ, α), p, α) : [ r 1, ) R is a solution of 1.1). Thus it follows from the definition of X, α) that t = is the only discontinuity of q, α), it is differentiable for t > r 1 and satisfies qt, α) = q, α) α t a jˆ q s r j, α)ds for t. It is a well known result see [4] of Diekmann et al.), that p t = Pr P X t P for all t r 1, hence q t Q for all t r 1. Moreover, formula q t, α) = 1 π eγt lim T ˆ T T e iωt γ + iω, α) dω holds for all t > by the Laplace transform technique [4, 9]. In order to estimate X t, α) dt, we estimate p t, α) dt and q t, α) dt. Proposition 4.1. Under hypothesis H1) for all α 1 α < α. Furthermore, lim α α) α α Proof. For all α [α 1, α ] and t r 1, νe µ ν π 1 + e µ π) ν p t, α) dt z λ, α) µ + ν ) 1 e µ π) ν 4 z p t, α) dt = iν, α ) n ) πα a n ). jr j sin r j ν ) a j sinr j ν ) 4.1) e λt p t, α) = Re z λ, α) e µt { ) = cosνt) 1 α a j r j e rjµ cosr j ν) z λ, α) )} + sinνt) α a j r j e rjµ sin r j ν). EJQTDE, 11 No. 36, p. 9
10 Since for all A, B R with A + B > we have Acosνt) + B sin νt) = A + B sin νt + η) with η [ π, π) satisfying sin η = A B and cosη = A + B A + B, we obtain that there exists η = η α) [ π, π) such that 4.) p t, α) = e µtsin νt + η) z λ, α), α [α 1, α ], t r 1. Choose < t 1 < t <... < t n <... such that νt j + η = jπ for all j 1. Then sin νt + η) > for all t t j, t j+1 ), j 1, and sinνt + η) < for all t t j 1, t j ), j 1. With this notation, Since it follows that t 1 p t, α) dt = ˆ t 1 for all α 1 α < α. In addition, ˆ t1 p t, α) dt = = = jˆ tj+1 1) p t, α) dt z λ, α) t j e µ ν η z λ, α) ν jˆ tj+1 1) e µt sin νt + η)dt t j jˆ j+1)π 1) e µ ν s sin sds. e µ ν s sin sds = ν e µ ν s µ + ν µ ν sins coss ), p t, α) dt = νe µ ν η e µ ν π + 1 ) z λ, α) µ + ν ) = νe µ ν η e µ ν π + 1 ) e µ ν π z λ, α) µ + ν ) 1 e µ ν π z λ, α) e µ ν η z λ, α) ν ˆ t1 { jπ e µt sin νt + η) dt ˆ π e µ ν s sinsds + = νe µ ν η e µ ν π + + e µ π) ν z λ, α) µ + ν. ) e µ ν jπ ˆ π } e µ ν s sin sds EJQTDE, 11 No. 36, p. 1
11 As e µ ν π + + e µ ν π + e µ ν π e µ ν π + 1 ) = 1 + e µ ν π 1 e µ ν π 1 e, µ ν π the last two results with η [ π, π) give the first statement of the proposition. t It is clear that 1 α. Hence lim α α) α α p t, α) dt is bounded with an upper bound independent of p t, α) dt = lim α α α α) = t 1 4 z iν, α) π lim α α p t, α) dt α α µ ) µ ν π ). 1 e µ ν π As lim x x)/1 e x ) = 1 and µ = µ ξ)α α) with some ξ α, α ), we obtain that lim α α) α α p t, α) dt = 4 z iν, α) π lim α α µ α). Now Corollary 3. yields the second statement of the proposition. It is a well known result that for each α there exists a constant L α) > such that q t, α) < L α) e γt for all t >, see [4, 9]. Next we construct an upper bound for q t, α) that is independent of α. Proposition 4.. If hypothesis H1) holds, then q t, α) < max {L 1 α), L α)} e γt for all t > and α [α 1, α ], where L 1 = L 1 α) and L = L α) are given by.1) and.), respectively. Consequently, q t, α) dt < max {L 1 α),l α)} γ Proof. Recall that for α [α 1, α ] and t >, T q t, α) = eγt π lim T ˆ T T e iωt γ + iω, α) dω, for all α [α 1, α ]. where γ + iω, α) = γ + iω + α n a je rjγ+iω). Partial integration gives that ˆ T e iωt [ ] e iωt T ˆ γ + iω, α) dω = 1 T e iωt z γ + iω, α) + it γ + iω, α) t γ + iω, α) dω. Since γ + iω, α) as ω and e iωt /t = 1/t for all t >, we conclude that ˆ T q t, α) = eγt π lim T T T T e iωt z γ + iω, α) t γ + iω, α) dω for α [α 1, α ] and t >. EJQTDE, 11 No. 36, p. 11
12 Hence It is clear that q t, α) eγt 1 + α πt ) a j r j e rjγ dω γ + iω, α). n γ + iω, α) γ + iω α a j e rjγ = n γ + ω α a j e rjγ ω if ω ω = α n a je rjγ. Therefore ˆ dω γ + iω, α) = dω γ + iω, α) {ˆ ω = 8 ω + dω γ + iω, α) + ˆ ω dω γ + iω, α). ˆ } 4 ω ω dω Thus for t r 1 and α [α 1, α ], q t, α) < L 1 α) e γt, where L 1 α) is defined by.1). We also need an estimate for t, r 1 ). It is clear from 1.)-1.3) that n )ˆ t X t, α) 1 + α a j X s, α) ds, t [, r 1 ], α [α 1, α ], ) thus X t) exp α r n 1 a j Lemma. From 4.) we see that p t, α) for t [, r 1 ] and α [α 1, α ] by Gronwall s z λ, α) for t, r 1) and α [α 1, α ], and this upper bound is finite by H1). Hence q t, α) X t) + p t, α) < L α) e γt for t, r 1 ) and α [α 1, α ] with L α) defined by.). Proof of Theorem.1. Proposition 4. implies that q t, α) dt is bounded on [α 1, α ]. Therefore Theorem.1 follows directly from the facts that Xt, α) = pt, α) + qt, α), t r 1, α [α 1, α ], and that the same limit holds for see Proposition 4.1. α α) p t, α) dt, EJQTDE, 11 No. 36, p. 1
13 then Note that if n =, r 1 = 1, r =, a 1 > a and hypothesis H1) is satisfied, a 1 cosν + a = and sin ν = 1 cos ν = 1 a. Hence Theorem.1 implies lim α α) α α a 1 X t, α) dt = 4 1 α a 1 cosν + α a 1 πα a 1 sin ν ) = α a + α a 1 πα a 1 a ). Proof of Theorem.. The upper bound for X t, α) dt in the theorem is simply the sum of the upper bounds for p t, α) dt and q t, α) dt given by Proposition 4.1 and Proposition 4., respectively. As we have already mentioned, the upper bound given by the Theorem. is sharp for parameters close to the critical value α : the upper estimate multiplied by α α) has the same limit at α as function α α) X t, α) dt. Indeed, as we see that µ = µ ξ)α α) with some ξ α, α ) and lim x x 1 e x = 1, lim α α) α α { νe µ ν π 1 + e µ π) } ν z λ, α) µ + ν ) 1 e µ π) + max {L 1, L } ν γ ν e µ ν π 1 + e µ π) ) ν α α µ ν = lim π ) α α z λ, α) µ + ν )π µ 1 e µ ν π 4 = z iν, α) π lim 1 α α µ α), which limit is the same as given by Theorem.1, see Corollary 3.. Proof of Theorem.3. Assume that not only hypothesis H1) but also H) is satisfied. In this case δ ν π/r 1 δ for all α [α 1, α ], hence sin r 1 ν) > sin r 1 δ) and sin r j ν) > for all α [α 1, α ] and j {1,...,n}. In consequence, 4.3) z λ, α) α a j r j e rjµ sin r j ν) α 1c for all α [α 1, α ] EJQTDE, 11 No. 36, p. 13
14 with c = a 1 r 1 sin r 1 δ). Also, µ > γ and by Corollary 3.3. α α µ K δ) for each α [α 1, α ] With Proposition 4.1, these estimates imply that ) α α µ ν α α) p t, α) dt π ) ν e µ ν π 1 + e µ π) ν µ 1 e µ ν π π z λ, α) µ + ν ) µ ν K δ) π ) e γ δ π 1 + e γ π) δ 1 e µ ν π πα 1 c for all α 1 α < α. The function x x/ 1 e x ) is strictly decreasing on, ), hence and α α) for all α 1 α < α. In addition, L α) = e γr1 e αr1 P n aj + for all α 1 α < α. µ ν π 1 e µ ν π < γ δ π 1 e γ δ π p t, α) dt K δ) γe γ δ π 1 + e γ π) δ α 1 c δ 1 e γ π) δ z ) P n e γr1 e αr1 aj + ) = λ, α) α 1 c L With Proposition 4. now we deduce that for all α 1 α < α, α α) X t, α) dt α α) p t, α) dt + α α) q t, α) dt ) γe γ δ π 1 + e γ δ π max { L1, L } α α 1 ) K δ) α 1 c δ 1 e γ π) +, δ γ where L 1 = sup α1 α α L 1 α) and L 1 α) is defined by.1). 5. The proof of Theorem.4 For Eq..3) and α > 1/e the eigenvalues are simple and appear in complex conjugate pairs ) λ j, λ j j= with Reλ > Reλ 1 >... > Reλ j j ) and Imλ j jπ, j + 1)π) for all j. EJQTDE, 11 No. 36, p. 14
15 In coherence with the previous sections, λ = µ + iν and λ = µ iν denote the leading eigenvalues λ and λ. As they are roots of the characteristic function, 5.1) µ + αe µ cosν = and ν αe µ sin ν =. In [1] Krisztin has verified that for α [3/, π/), 1 < µ <, 1.54 < ν < π, λ > 1, lim α π µ = and lim ν = π α π, moreover the real parts of the remaining eigenvalues are smaller than 1. 5.) With equations 5.1) Proposition 3.1 implies µ α) = α e µ sin ν) 1 αe µ cosν)e µ cosν 1 αe µ cosν) + αe µ sinν) = 1 ν µ)µ α 1 + µ), α [3/, π/]. + ν For all α [3/, π/], ν > 1 and.5 < 1 + µ)µ, hence < µ α) < 1/α. As µ α) = µ ξ)π/ α) for all α [3/, π/) with some ξ α, π/), we obtain that µ α) < 1 ξ π 3 ) < π 3 3 ) <.5, α [3/, π/). We mention that numerical approximation yields µ >.33 for all α [3/, π/), thus conditions H1) and H) hold with γ =.1, δ = 1.54 and K = 15. According to Theorem.3, X t, α) dt 9. In this section we give a better estimate without using any numerical approximation. Relation 4.1) and equations 5.1) imply that 5.3) z λ, α) = 1 + µ) + ν and e µt 5.4) p t, α) = 1 + µ) [1 + µ)cosνt) + ν sinνt)] + ν for t 1 and α [3/, π/], see also [1]. Proposition 5.1. For each α [3/, π/), π ) α p t, α) dt Proof. Proposition 4.1 and relation 5.3) implies that for all α [3/, π/) we have EJQTDE, 11 No. 36, p. 15
16 π α ) with some ξ α, π/). p t, α) dt e µ ν π 1 + e µ π) ν ν π z λ, α) µ ν π 1 + e µ π) ν e π 1 + µ) + ν µ + ν µ ν π 1 e µ ν π ) µ ν π 1 e µ ν π ) 1 µ ξ) ) π α ) As µ [.5, ] and ν [1, 54, π/] for each α [3/, π/), it follows from result 5.) that < 1 µ α) α max µ,ν) [.5,] [1,54,π/] 1 + µ ) [ 1 + µ 3 ν, α µ)µ, π ). The expression on the right hand side is strictly decreasing in ν [1, 54, π/]. In addition, as ) 1 + µ µ ν = ν µ + 1) µ)µ ν µ)µ) > for µ, ν) [.5, ] [1, 54, π/], it is strictly increasing in µ [.5, ]. We deduce that π ) α p t, α) dt 1 µ α) π ) 1.54 for all α [ 3, π ). Recall that the function x x/ 1 e x ) is strictly decreasing on, ). Hence ) e π 1 + e π π π ) π ) 1 e π 1.54 for all α [3/, π/), and this upper bound is smaller than In order to get a better estimate for q t, α) dt, we apply an approach different from that of Proposition 4.. The fact that there is only one delay is crucial here. We use the discrete Lyapunov functional V of Mallet-Paret and Sell introduced in [1]. V ϕ) counts the sign changes of ϕ C \ {} if it is an odd number or infinity, otherwise V ϕ) is the number of sign changes plus one. Then V ϕ) {1, 3,...} { }. The map t V x t ) is monotone nonincreasing along solutions of Eq..3). In addition, V is upper semi-continuous: for each ϕ C \ {} and ϕ n ) C \ {} with ϕ n ϕ as n, V ϕ) liminf n V ϕ n ). Proposition 5.. For each α [3/, π/] and t 1, V q t ) 3, that is q has at least two sign changes on each subinterval [t 1, t] of [, ). Proof. Suppose for contradiction that V q s ) = 1 for some s 1. Then the monotone property of V gives that V q t ) = 1 for all t s. By a result of Cao [3] see also Arino [1]) and by the Gronwall-Bellmann inequality there exist C 1 >, C > EJQTDE, 11 No. 36, p. 16
17 and t such that 5.5) C 1 q t q t+1 C q t for all t > t. and For n 1 define y n : [ n, ) R by y n t) = q n + t)/ q n. Then ẏ n t) = αy n t 1) for t > n + 1, C 1 yt n yt+1 n C yt n if n + t > t. The Arzelà Ascoli Theorem can be applied to find a subsequence y n k ) k=1 and a C 1 -function y : R R so that y n k t) y t), ẏ n k t) y t) as k uniformly on compact subintervals of R, moreover ẏ t) = αy t 1) for t R. As Q is closed, y = lim k y n k Q with y = 1. In addition, C 1 y t y t+1 C y t for all t R. Hence for all n {, 1,,...}, y n C1 1 y n+1 C1 y n+... C1 n y 1 + C1 1 ) n y. Let t be arbitrary and choose integer n so that n + 1) < t n holds. Then y t y n+1) + y n 1 + C 1 1 ) y n 1 + C 1 ) n+1 y. As ) 1 + C 1 n 1 = e nln1+c 1 1 ) e tln1+c 1 1 ), we conclude that y t Ae Bt for all t with A = 1 + C1 1 ) y > and B = ln 1 + C1 1 ) >. Choose c > B so that Rez c for all roots of the characteristic function. The space C has the decomposition C = ˆP ˆQ, where ˆP is the realified generalized eigenspace of the generator of the semigroup T t)) t associated with the eigenvalues having real parts greater than c, and ˆQ is the realified generalized eigenspace associated with the rest of the spectrum. By [9], there is M > so that 1 T t)ϕ Me ct ϕ for all t and ϕ ˆQ. EJQTDE, 11 No. 36, p. 17
18 Let t and σ t. Then T Pr ˆQy t = t σ)pr ˆQy σ Me ct σ) Pr ˆQy σ M Pr ˆQ e ct σ) y σ M Pr ˆQ e ct σ) Ae Bσ = M Pr ˆQ e ct Ae B c)σ as σ. It follows that Pr ˆQy t = and y t ˆP for all t. If subspace ˆP is trivial, i.e. there are no eigenvalues with real parts greater than c, then the previous result implies y =, a contradiction to y = 1. Otherwise y t) = A j e Reλjt cosimλ j t + B j ), t, j= for some real numbers A, B,...,A n, B n and integer N > so that λ = λ, λ = λ, λ 1, λ 1,...,λ N, λ N are the eigenvalues with real parts greater than c. The upper semi-continuity of V and V q t ) = 1, t s, combined yield V y t ) = 1 for all t R. As Imλ j jπ, j + 1)π), j, it follows that A 1 = A =... = A N =. This means that y t P for all t. In particular y P, a contradiction to y Q \ {}. This completes the proof. We can use this result to give an explicit estimate for the growth of q on [ 1, ). In the next proposition r denotes the integer part of the positive real number r. Proposition 5.3. For each α [3/, π/] and t 1, α ) k q t) q, where q = sup q s) and k = 1 s 1 Proof. The statement is clearly true for 1 t < 1/. It is enough to show that if t + 1). 3 t 1 and sup q s) m for some m >, s [t 3,t] then q t) α [ m for all t t, t + 3 ], so we confirm this latter assertion. EJQTDE, 11 No. 36, p. 18
19 For all t 1/ and t [t, t + 3/], there exists z [t 1/, t + 1/] with q z) =, see Proposition 5.. Hence for t [t, t + 3/], ˆ t 1 q t) = q t) q z) = α q s)ds α t z sup s [t 3,t 1 ] z 1 q s) It follows that and α sup t [t + 1,t+1] q t) α sup q s). s [t 3,t 1 ] sup q t) α t [t,t + m, ] 1 sup q t) α { t [t 1,t + max m, α } m = α m 1 ] sup q t) α sup q t) α {m, t [t +1,t + ] 3 t [t 1,t+1] max α } m = α m. The previous statement shows we need an upper bound for q. Proposition 5.4. For all α [3/, π/], q = sup 1 s 1/ q s) 1. Proof. Set α [3/, π/] arbitrarily. Differentiating 5.4) we get p t) = = e µt [ µ + µ 1 + µ) + ν ) cosνt) ν sin νt) ] + ν e µt 1 + µ) + ν ν cosνt) [ µ + µ + ν ν ] tanνt). Note that as 1.54 < ν π/ and.5 < µ, we have cosνt) > for all t [ 1, π/ ν)), and in addition µ + µ + ν >. It follows that there exists t, π/ ν)) such that tan νt ) = µ + µ + ν, ν p increases on [ 1, t ] and decreases on [t, π/ ν)). Clearly, e µ p 1) = 1 + µ) [1 + µ)cosν ν sin ν] < + ν EJQTDE, 11 No. 36, p. 19
20 as 1 + µ)cosν < 1 and ν sin ν > 1.54 sin 1.54) > 1 for µ.5, ] and ν 1.54, π/]. Therefore p 1) e µ ν 1 + µ) + ν. As the right hand side is decreasing in ν for ν 1.54, π/] and it is decreasing in µ for µ.5, ], we deduce that Also we find that < p ) = p 1) e < µ) 1 + µ) + ν 1 + ν < 1, e µ < p1) = 1 + µ) [1 + µ)cosν + ν sin ν], + ν ν 1 + µ) + ν π < 1 and e µt pt ) = [1 1 + µ) + µ) + µ + µ + ν ] ν cosνt ) = e µt cosνt ). + ν ν Clearly, p t) 1, 1) for all t [ 1, ]. In case t > 1 one has p t), 1) for all t [, 1]. Otherwise p t), ) for all t [, 1]. Using that q t) = p t) for 1 t < and q t) = 1 p t) for t 1, we obtain that q 1. Now we are able to estimate q t, α) dt. The bound given by the subsequent proposition is substantially better than bound 3e presented in paper [1] and the bound given by Proposition 4.. Proposition 5.5. For each α [3/, π/], q t, α) dt + π 4 π. Proof. By Propositions 5.3 and 5.4, sup 1 s 1/ q s, α) 1 for all α [3/, π/], and q t, α) π/4) k for all α [3/, π/] and t, where k is the greatest integer with k t + 1)/3. Thus q t, α) dt { π 4 + π 4 ) + π 4 ) } = + π 4 π. EJQTDE, 11 No. 36, p.
21 Proof of Theorem.4. The statement of the theorem follows directly from Proposition 5.1 and Proposition 5.5. Consider the linear equation 6. An example with two delays ẋt) = α xt 1) + xt )) with < α In this case the characteristic function is 6.1) z, α) = z + α e z + e z). For the eigenvalues z = u + iv, 6.) u + α e u cosv + e u cosv) ) = and 6.3) v α e u sin v + e u sin v) ) =. There are purely imaginary eigenvalues ±iπ/3 for α = π/ 3 3 ). π 3 3. We examine the location of the eigenvalues on the complex plane for < α π/ 3 3 ). Suppose < α π/ 3 3 ) and z = u + iv is an eigenvalue with u and v π. Then it follows from 6.3) that that is π v π 3 3 e u sinv + e u sin v) ) π 3 3 e u, u 1 ln 3 3 <.4. It is also known that there is exactly one pair of eigenvalues λ, λ ) in the subset {z C : Imz < π} of the complex plane, see [13]. Numerical approximation gives that for each α [.58, π/ 3 3 )), there is an eigenvalue λ = λα) in the rectangle see Fig. 1. {z C :. < Rez < and 1.3 Imz 1.5}, Hence hypotheses H1) and H) are satisfied with α = π/ 3 3 ), α 1 =.58, γ =.4, δ =.5 and K =. Theorem.1 now gives that )ˆ π lim α π α X t, α) dt = π 6 3 3π ) + π EJQTDE, 11 No. 36, p. 1
22 Figure 6.1. Curve h.58, π/ 3 i 3 α λ α) C In addition, Theorem.3 can be applied, and we get )ˆ π 3 3 α X t, α) dt 55 for all.58 α < π/ 3 3 ) by Theorem Applications In this section, let n 1 be an integer, r 1,..., r n be real numbers with r 1 > r >... > r n, a 1,..., a n be positive numbers, α > α 1 >. In case n = 1, we assume r 1 >. In papers [7, 8] Győri, Hartung and Turi has studied linear equations with perturbed time lags in the form 7.1) ẋt) = a j xt r j η j t)), where η j : [, ) [, ) are piecewise continuous functions for all j {1,...,n}. Suppose that the solution of Eq.7.1) with η j, j {1,...,n}, is asymptotically stable. Then using estimates on the integral of the fundamental solution, they gave sufficient conditions on η j to guarantee the asymptotic stability of the solution of Eq.7.1). We cite the following result as an example: if 1 a j limsup η j t) < n ) t a, j X t) dt then the trivial solution of Eq. 7.1) is asymptotically stable. EJQTDE, 11 No. 36, p.
23 We present two more applications in connection with the nonlinear equation 7.) ẏ t) = α a j y t r j ) + f y t, α) with parameter α [α 1, α ) and a continuously differentiable nonlinear function f : C [α 1, α ) R. We assume the following hypothesis: H3): To each α [α 1, α ) there correspond positive constants k 1 = k 1 α) and δ 1 = δ 1 α) such that f ϕ, α) k 1 α) ϕ for all α [α 1, α ) and ϕ C with ϕ < δ 1 α). By the variation-of-constants formula, the solution y ϕ : [ 1, ) R of Eq.7.) with y ϕ = ϕ satisfies the equation 7.3) y ϕ t) = x ϕ t) + ˆ t X t s, α) f y s, α)ds for t, where x ϕ : [ 1, ) R is the solution of the linear variational equation 1.1) with initial segment ϕ [4, 9]. Using this formula and estimates for the integral of the fundamental solution, the existence of periodic solutions of small amplitude can be excluded for Eq. 7.). Theorem 7.1. Assume that H1) holds for the linear variational equation 1.1), and f : R [α 1, α ) R satisfies H3). If α [α 1, α ) and p : R R is a nonconstant periodic solution of Eq. 7.), then { max p t) min t R 1 k 1 α) X t, α) dt, δ 1 α) Proof. Fix α [α 1, α ) and let p : R R be a periodic solution of Eq. 7.) with { } 1 A = max p t) < min t R k 1 X t, α) dt, δ 1. Then clearly f p t, α) k 1 p t k 1 A for all t R. By the periodicity of p, there is a sequence t n ) in [, ) so that t n as n and p tn = A for all n. Then the variation-of-constants formula 7.3) yields ˆ tn A = p t n ) x p t n ) + k 1 A X t n s, α) ds for all n. Hypothesis H1) implies lim t x p t) = [9]. Letting n, we obtain that A k 1 A X s, α) ds. }. EJQTDE, 11 No. 36, p. 3
24 Hence 1 A k 1 X t, α) dt, a contradiction to our initial assumption. For each ϕ C, solution x ϕ : [ 1, ) R of Eq.1.1) can be expressed in form 7.4) x ϕ t) = X t, α)ϕ) α a jˆ X t r j s, α)ϕs)ds for t, r j j= see [9]. Hence using equality 4.) and Proposition 4., one can determine explicit constants M = M α) 1 and b = b α) > such that x ϕ t M ϕ e bt for each t. As a second application, we estimate the stability region of the solution of Eq.7.) in phase space C. Theorem 7.. Assume that H1) holds for the linear variational equation 1.1), and f : R [α 1, α ) R satisfies H3). In addition, suppose that δ 1 α) in H3) is chosen so small that k 1 α) δ 1 α) X t, α) dt < 1 for all α [α 1, α ). If α [α 1, α ) and ϕ C with ϕ δ 1 / M α)), then for the solution y ϕ : [ 1, ) R of Eq. 7.) with initial function ϕ, lim t y ϕ t) =. Proof. Set α [α 1, α ) and ϕ C with ϕ < δ 1 / M). First we claim that y ϕ t) < δ 1 for all t 1. Suppose for contradiction that there is a minimal t > with y ϕ t ) = δ 1. Then the variation-of-constants formula 7.3) gives that δ 1 = y ϕ t ) x ϕ t ) + ˆ t M ϕ e bt + k 1 δ 1 X t s, α)f y s, α) ds ˆ t X t s, α) ds, where x ϕ : [ 1, ) R is the solution of the linear variational equation 1.1) with x ϕ = ϕ. It follows that δ 1 1 k 1 δ 1 ) X t, α) dt M ϕ e bt M ϕ. As 1 k 1 δ 1 X t, α) dt > 1/, we obtain that ϕ > δ 1 / M), which contradicts our initial assumption. So y ϕ t) < δ 1 for all t 1. Set A = limsup t y ϕ t). The previous step implies that A [, δ 1 ]. We have to show that A =. EJQTDE, 11 No. 36, p. 4
25 Suppose for contradiction that A = δ 1. Choose a sequence t n ) in [, ) so that t n and y ϕ t n ) δ 1 as n. Then by the variation-of-constants formula 7.3), y ϕ t n ) x ϕ t n ) + k 1 δ 1 ˆ tn X t n s, α) ds for all n. Clearly, lim t x ϕ t) =. Letting n, we obtain that δ 1 k 1 δ 1 X s, α) ds, that is 1 k 1 δ 1 X s, α) ds, which contradicts the choice of δ 1. So A < δ 1. Suppose that A, δ 1 ). Set ε, 1 ) so small that inequality 1 + ε)a < δ 1 also holds. As before, there exists a sequence t n ) in [, ) so that t n and y ϕ t n ) A as n. In addition, there is a threshold number T > so that y ϕ t) < 1 + ε)a for t T 1. Then hypothesis H3) implies that f y ϕ t, α) k 1δ 1 for all t < T, and f y ϕ t, α) k ε) A for all t T. By the variationof-constants formula, y ϕ t n ) x ϕ t n ) + k 1 δ 1 ˆ T X t n s, α) ds + ˆ tn k ε) A X t n s, α) ds = x ϕ t n ) + k 1 δ 1 T ˆ tn t n T ˆ tn T X s, α) ds + k ε) A X s, α) ds for all sufficiently large n. Since X s, α) ds < for all α [α 1, α ), we see t that n t n T X s, α) ds as n. So letting n, we conclude that A k ε) A X s, α) ds, that is 1 A k ε) X t, α) dt > 1 + ε) δ 1. This result contradicts the fact that A < δ 1. It follows that A = limsup t y ϕ t) =, and the proof is complete. As an example, consider Wright s equation ) 7.5) ẏ t) = α e yt 1) 1 with parameter α [3/, π/). Now the linear variational equation along the solution is Eq..3). For the fundamental solution X of Eq..3), Theorem.4 EJQTDE, 11 No. 36, p. 5
26 gives that for all α [3/, π/). π/ α) X t, α) dt.36 Eq.7.5) can be written in form 7.) with [ 3 f : C, π ) ) ϕ, α) α e ϕ 1) 1 ϕ 1) For α [3/, π/) and ϕ C with ϕ < δ 1 < 3, ϕ k f ϕ, α) α α ϕ k! So H3) holds with k= k= δ1 3 δ 1 α), 3) and k 1 α) π 4 ) k π δ 1. R. 3 3 δ 1 ϕ. Set δ 1 α).4 and k 1 α).8 for example. Then it follows from Theorem 7.1 that whenever p : R R is a nonconstant periodic solution of Eq.7.) for some α [3/, π/), then π ) max p t).5 t R α. We note that it is verified in [] that max p t) π ) π ) t R π α.64 α. To apply Theorem 7., we decrease δ 1 α) so that k 1 α)δ 1 α) X t, α) dt.8 δ 1 α).36 π α < 1 for all α [3/, π/); we set δ 1 α) =.6 π/ α) for each α [3/, π/). At last we need constants M α) > 1 and b α) > such that for the solution x ϕ : [ 1, ) R of Eq..3), x ϕ t M α) ϕ e bα)t for all t. Recall that for α [3/, π/),.5 < µ < and 1.54 < ν < π, where µ and ν are the real and imaginary parts of the leading eigenvalue λ = λα) in the spectrum of the generator of the semigroup defined by the solution operators of Eq..3). Hence by 5.4), 1 + µ + ν) p t, α) 1 + µ) + ν eµt ) 1 + π eµt < 1.58e µt EJQTDE, 11 No. 36, p. 6
27 for all α [3/, π/) and t 1. In addition, Proposition 5.3 and Proposition 5.4 imply q t, α) π ) t+1) 3 π ) 3 t = π e ln π 3 4) t 1.9e.16t for all α [3/, π/) and t. It follows that X t, α) p t, α) + q t, α) )e µt =.67e µt for all α [3/, π/) and t. Recall that X t, α) = for all t [ 1, ), hence the previous estimate can be used for all t 1. According to formula 7.4), x ϕ t) X t, α) ϕ) + α ˆ 1 X t 1 s, α) ϕs) ds π e.5) ϕ e µt 7.8 ϕ e µt for each t. Therefore we may set M α) = 7.8 and b α) = µ α) for α [3/, π/). According to Theorem 7., if ϕ C with ϕ.18 π/ α), then for the solution y ϕ : [ 1, ) R of Eq.7.) with initial function ϕ, lim t y ϕ t) =. Acknowledgments. This research was partially supported by the Hungarian Research Fund, Grant no. K75517, and by the TÁMOP 4.../8/1/8-8 program of the Hungarian National Development Agency. References [1] Arino, O., A note on: The discrete Lyapunov function for scalar differential delay equations. J. Differential Equations ), no. 1, [] Bánhelyi B., Csendes T., Krisztin T. and Neumaier A., Global attractivity of the zero solution for Wright s equation. Manuscript. [3] Cao, Y., The discrete Lyapunov function for scalar differential delay equations. J. Differential Equations ), no., [4] Diekmann, O., van Gils, S. A., Verduyn Lunel, S. M., and Walther, H.-O. Delay equations. Functional, complex, and nonlinear analysis. Springer-Verlag, New York 1995). [5] Győri I., Global attractivity in a perturbed linear delay differential equation. Appl. Anal ), no. 3-4, [6] Győri, I., Hartung, F., Fundamental solution and asymptotic stability of linear delay differential equations. Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 13 6), no., EJQTDE, 11 No. 36, p. 7
28 [7] Győri, I., Hartung, F., Stability in delay perturbed differential and difference equations. Topics in functional differential and difference equations Lisbon, 1999), , Fields Inst. Commun., 9, Amer. Math. Soc., Providence, RI 1). [8] Győri, I., Hartung, F., Turi, J., Preservation of stability in delay equations under delay perturbations. J. Math. Anal. Appl. 1998), no. 1, [9] Hale, J. K., and Verduyn Lunel, S. M. Introduction to functional-differential equations, Springer-Verlag, New York 1993). [1] Krisztin, T., On the fundamental solution of a linear delay differential equation. International Journal of Qualitative Theory of Differential Equations and Applications 3 9), no. 1-, [11] Krisztin T., Röst G., On the size of center manifolds with an application to Wright s equation. In preparation. [1] Mallet-Paret, J., Sell, G. R., Systems of differential delay equations: Floquet multipliers and discrete Lyapunov functions. J. Differential Equations ), no., [13] Nussbaum, R. D., Potter, A. J. B., Cyclic differential equations and period three solutions of differential-delay equations. J. Differential Equations ), no. 3, Received April, 11) EJQTDE, 11 No. 36, p. 8
Periodic Orbits and the Global Attractor for Delayed Monotone Negative Feedback
Periodic Orbits and the Global Attractor for Delayed Monotone Negative Feedback Tibor Krisztin 1 Bolyai Institute, University of Szeged Aradi vértanúk tere 1, H-6720 Szeged, Hungary E-mail: krisztin@math.u-szeged.hu
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationDifferentiability of solutions with respect to the delay function in functional differential equations
Electronic Journal of Qualitative Theory of Differential Equations 216, No. 73, 1 16; doi: 1.14232/ejqtde.216.1.73 http://www.math.u-szeged.hu/ejqtde/ Differentiability of solutions with respect to the
More informationON THE GLOBAL ATTRACTOR OF DELAY DIFFERENTIAL EQUATIONS WITH UNIMODAL FEEDBACK. Eduardo Liz. Gergely Röst. (Communicated by Jianhong Wu)
DISCRETE AND CONTINUOUS doi:10.3934/dcds.2009.24.1215 DYNAMICAL SYSTEMS Volume 24, Number 4, August 2009 pp. 1215 1224 ON THE GLOBAL ATTRACTOR OF DELAY DIFFERENTIAL EQUATIONS WITH UNIMODAL FEEDBACK Eduardo
More informationDIFFERENTIABLE SEMIFLOWS FOR DIFFERENTIAL EQUATIONS WITH STATE-DEPENDENT DELAYS
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XLI 23 DIFFERENTIABLE SEMIFLOWS FOR DIFFERENTIAL EQUATIONS WITH STATE-DEPENDENT DELAYS by Hans-Otto Walther Careful modelling of systems governed
More informationHamad Talibi Alaoui and Radouane Yafia. 1. Generalities and the Reduced System
FACTA UNIVERSITATIS (NIŠ) Ser. Math. Inform. Vol. 22, No. 1 (27), pp. 21 32 SUPERCRITICAL HOPF BIFURCATION IN DELAY DIFFERENTIAL EQUATIONS AN ELEMENTARY PROOF OF EXCHANGE OF STABILITY Hamad Talibi Alaoui
More informationON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay
ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov
More informationConvergence Rate of Nonlinear Switched Systems
Convergence Rate of Nonlinear Switched Systems Philippe JOUAN and Saïd NACIRI arxiv:1511.01737v1 [math.oc] 5 Nov 2015 January 23, 2018 Abstract This paper is concerned with the convergence rate of the
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 3. Fundamental properties IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs2012/ 2012 1 Example Consider the system ẋ = f
More informationSufficient conditions for the existence of global solutions of delayed differential equations
J. Math. Anal. Appl. 318 2006 611 625 www.elsevier.com/locate/jmaa Sufficient conditions for the existence of global solutions of delayed differential equations J. Diblík a,,1,n.koksch b a Brno University
More informationPiecewise Smooth Solutions to the Burgers-Hilbert Equation
Piecewise Smooth Solutions to the Burgers-Hilbert Equation Alberto Bressan and Tianyou Zhang Department of Mathematics, Penn State University, University Park, Pa 68, USA e-mails: bressan@mathpsuedu, zhang
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationHigher Order Averaging : periodic solutions, linear systems and an application
Higher Order Averaging : periodic solutions, linear systems and an application Hartono and A.H.P. van der Burgh Faculty of Information Technology and Systems, Department of Applied Mathematical Analysis,
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More informationCLASSIFICATIONS 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 information8 Periodic Linear Di erential Equations - Floquet Theory
8 Periodic Linear Di erential Equations - Floquet Theory The general theory of time varying linear di erential equations _x(t) = A(t)x(t) is still amazingly incomplete. Only for certain classes of functions
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationDelayed loss of stability in singularly perturbed finite-dimensional gradient flows
Delayed loss of stability in singularly perturbed finite-dimensional gradient flows Giovanni Scilla Department of Mathematics and Applications R. Caccioppoli University of Naples Federico II Via Cintia,
More informationINVARIANT SUBSPACES FOR CERTAIN FINITE-RANK PERTURBATIONS OF DIAGONAL OPERATORS. Quanlei Fang and Jingbo Xia
INVARIANT SUBSPACES FOR CERTAIN FINITE-RANK PERTURBATIONS OF DIAGONAL OPERATORS Quanlei Fang and Jingbo Xia Abstract. Suppose that {e k } is an orthonormal basis for a separable, infinite-dimensional Hilbert
More informationNonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1
Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems p. 1/1 p. 2/1 Converse Lyapunov Theorem Exponential Stability Let x = 0 be an exponentially stable equilibrium
More informationEVANESCENT SOLUTIONS FOR LINEAR ORDINARY DIFFERENTIAL EQUATIONS
EVANESCENT SOLUTIONS FOR LINEAR ORDINARY DIFFERENTIAL EQUATIONS Cezar Avramescu Abstract The problem of existence of the solutions for ordinary differential equations vanishing at ± is considered. AMS
More informationExponential stability of families of linear delay systems
Exponential stability of families of linear delay systems F. Wirth Zentrum für Technomathematik Universität Bremen 28334 Bremen, Germany fabian@math.uni-bremen.de Keywords: Abstract Stability, delay systems,
More informationOn the simplest expression of the perturbed Moore Penrose metric generalized inverse
Annals of the University of Bucharest (mathematical series) 4 (LXII) (2013), 433 446 On the simplest expression of the perturbed Moore Penrose metric generalized inverse Jianbing Cao and Yifeng Xue Communicated
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 informationExamples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling
1 Introduction Many natural processes can be viewed as dynamical systems, where the system is represented by a set of state variables and its evolution governed by a set of differential equations. Examples
More informationIdentification of Parameters in Neutral Functional Differential Equations with State-Dependent Delays
To appear in the proceedings of 44th IEEE Conference on Decision and Control and European Control Conference ECC 5, Seville, Spain. -5 December 5. Identification of Parameters in Neutral Functional Differential
More informationASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS. Tian Ma. Shouhong Wang
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 11, Number 1, July 004 pp. 189 04 ASYMPTOTIC STRUCTURE FOR SOLUTIONS OF THE NAVIER STOKES EQUATIONS Tian Ma Department of
More informationA note on linear differential equations with periodic coefficients.
A note on linear differential equations with periodic coefficients. Maite Grau (1) and Daniel Peralta-Salas (2) (1) Departament de Matemàtica. Universitat de Lleida. Avda. Jaume II, 69. 251 Lleida, Spain.
More informationNonlinear Control Lecture 5: Stability Analysis II
Nonlinear Control Lecture 5: Stability Analysis II Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2010 Farzaneh Abdollahi Nonlinear Control Lecture 5 1/41
More informationRelative Controllability of Fractional Dynamical Systems with Multiple Delays in Control
Chapter 4 Relative Controllability of Fractional Dynamical Systems with Multiple Delays in Control 4.1 Introduction A mathematical model for the dynamical systems with delayed controls denoting time varying
More informationHAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
Georgian Mathematical Journal Volume 9 (2002), Number 3, 591 600 NONEXPANSIVE MAPPINGS AND ITERATIVE METHODS IN UNIFORMLY CONVEX BANACH SPACES HAIYUN ZHOU, RAVI P. AGARWAL, YEOL JE CHO, AND YONG SOO KIM
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationOn the bang-bang property of time optimal controls for infinite dimensional linear systems
On the bang-bang property of time optimal controls for infinite dimensional linear systems Marius Tucsnak Université de Lorraine Paris, 6 janvier 2012 Notation and problem statement (I) Notation: X (the
More informationSYNCHRONIZATION OF NONAUTONOMOUS DYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 003003, No. 39, pp. 1 10. ISSN: 107-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu login: ftp SYNCHRONIZATION OF
More informationNONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS. Shouchuan Hu Nikolas S. Papageorgiou. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 29, 327 338 NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS Shouchuan Hu Nikolas S. Papageorgiou
More informationKey words. delay differential equations, characteristic matrix, stability of periodic orbits
CHARACTERISTIC MATRICES FOR LINEAR PERIODIC DELAY DIFFERENTIAL EQUATIONS JAN SIEBER AND ROBERT SZALAI Abstract. Szalai et al. (SIAM J. on Sci. Comp. 28(4), 2006) gave a general construction for characteristic
More informationGENERALIZATION OF GRONWALL S INEQUALITY AND ITS APPLICATIONS IN FUNCTIONAL DIFFERENTIAL EQUATIONS
Communications in Applied Analysis 19 (215), 679 688 GENERALIZATION OF GRONWALL S INEQUALITY AND ITS APPLICATIONS IN FUNCTIONAL DIFFERENTIAL EQUATIONS TINGXIU WANG Department of Mathematics, Texas A&M
More informationAn equivalent form, the so-called delayed logistic equation or Hutchinson s equation
GLOBAL ATTRACTIVITY OF THE ZERO SOLUTION FOR WRIGHT S EQUATION BALÁZS BÁNHELYI, TIBOR CSENDES, TIBOR KRISZTIN, AND ARNOLD NEUMAIER Abstract. In a paper published in 1955, E.M. Wright proved that all solutions
More informationBanach Journal of Mathematical Analysis ISSN: (electronic)
Banach J. Math. Anal. 1 (2007), no. 1, 56 65 Banach Journal of Mathematical Analysis ISSN: 1735-8787 (electronic) http://www.math-analysis.org SOME REMARKS ON STABILITY AND SOLVABILITY OF LINEAR FUNCTIONAL
More informationIntroduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems
p. 1/5 Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems p. 2/5 Time-varying Systems ẋ = f(t, x) f(t, x) is piecewise continuous in t and locally Lipschitz in x for all t
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationSingularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations
Singularities and Laplacians in Boundary Value Problems for Nonlinear Ordinary Differential Equations Irena Rachůnková, Svatoslav Staněk, Department of Mathematics, Palacký University, 779 OLOMOUC, Tomkova
More informationXiyou Cheng Zhitao Zhang. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 2009, 267 277 EXISTENCE OF POSITIVE SOLUTIONS TO SYSTEMS OF NONLINEAR INTEGRAL OR DIFFERENTIAL EQUATIONS Xiyou
More informationOscillation Criteria for Certain nth Order Differential Equations with Deviating Arguments
Journal of Mathematical Analysis Applications 6, 601 6 001) doi:10.1006/jmaa.001.7571, available online at http://www.idealibrary.com on Oscillation Criteria for Certain nth Order Differential Equations
More informationNONLINEAR THREE POINT BOUNDARY VALUE PROBLEM
SARAJEVO JOURNAL OF MATHEMATICS Vol.8 (2) (212), 11 16 NONLINEAR THREE POINT BOUNDARY VALUE PROBLEM A. GUEZANE-LAKOUD AND A. FRIOUI Abstract. In this work, we establish sufficient conditions for the existence
More informationEcon 204 Differential Equations. 1 Existence and Uniqueness of Solutions
Econ 4 Differential Equations In this supplement, we use the methods we have developed so far to study differential equations. 1 Existence and Uniqueness of Solutions Definition 1 A differential equation
More informationODE Final exam - Solutions
ODE Final exam - Solutions May 3, 018 1 Computational questions (30 For all the following ODE s with given initial condition, find the expression of the solution as a function of the time variable t You
More informationGLOBAL STABILITY OF SIR MODELS WITH NONLINEAR INCIDENCE AND DISCONTINUOUS TREATMENT
Electronic Journal of Differential Equations, Vol. 2015 (2015), No. 304, pp. 1 8. SSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GLOBAL STABLTY
More informationAn introduction to Mathematical Theory of Control
An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018
More informationOn periodic solutions of superquadratic Hamiltonian systems
Electronic Journal of Differential Equations, Vol. 22(22), No. 8, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) On periodic solutions
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationBIBO stabilization of feedback control systems with time dependent delays
to appear in Applied Mathematics and Computation BIBO stabilization of feedback control systems with time dependent delays Essam Awwad a,b, István Győri a and Ferenc Hartung a a Department of Mathematics,
More informationAsymptotic stability of solutions of a class of neutral differential equations with multiple deviating arguments
Bull. Math. Soc. Sci. Math. Roumanie Tome 57(15) No. 1, 14, 11 13 Asymptotic stability of solutions of a class of neutral differential equations with multiple deviating arguments by Cemil Tunç Abstract
More informationNONLINEAR DIFFERENTIAL INEQUALITY. 1. Introduction. In this paper the following nonlinear differential inequality
M athematical Inequalities & Applications [2407] First Galley Proofs NONLINEAR DIFFERENTIAL INEQUALITY N. S. HOANG AND A. G. RAMM Abstract. A nonlinear differential inequality is formulated in the paper.
More informationarxiv: v1 [math.oc] 22 Sep 2016
EUIVALENCE BETWEEN MINIMAL TIME AND MINIMAL NORM CONTROL PROBLEMS FOR THE HEAT EUATION SHULIN IN AND GENGSHENG WANG arxiv:1609.06860v1 [math.oc] 22 Sep 2016 Abstract. This paper presents the equivalence
More information6 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 informationSome Properties of NSFDEs
Chenggui Yuan (Swansea University) Some Properties of NSFDEs 1 / 41 Some Properties of NSFDEs Chenggui Yuan Swansea University Chenggui Yuan (Swansea University) Some Properties of NSFDEs 2 / 41 Outline
More informationMA5510 Ordinary Differential Equation, Fall, 2014
MA551 Ordinary Differential Equation, Fall, 214 9/3/214 Linear systems of ordinary differential equations: Consider the first order linear equation for some constant a. The solution is given by such that
More informationLOCAL FIXED POINT THEORY INVOLVING THREE OPERATORS IN BANACH ALGEBRAS. B. C. Dhage. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 24, 24, 377 386 LOCAL FIXED POINT THEORY INVOLVING THREE OPERATORS IN BANACH ALGEBRAS B. C. Dhage Abstract. The present
More informationEntrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.
Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the
More informationIn this paper we study periodic solutions of a second order differential equation
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 5, 1995, 385 396 A GRANAS TYPE APPROACH TO SOME CONTINUATION THEOREMS AND PERIODIC BOUNDARY VALUE PROBLEMS WITH IMPULSES
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear
More informationScalar curvature and the Thurston norm
Scalar curvature and the Thurston norm P. B. Kronheimer 1 andt.s.mrowka 2 Harvard University, CAMBRIDGE MA 02138 Massachusetts Institute of Technology, CAMBRIDGE MA 02139 1. Introduction Let Y be a closed,
More informationAsymptotic behavior of Ginzburg-Landau equations of superfluidity
Communications to SIMAI Congress, ISSN 1827-9015, Vol. 3 (2009) 200 (12pp) DOI: 10.1685/CSC09200 Asymptotic behavior of Ginzburg-Landau equations of superfluidity Alessia Berti 1, Valeria Berti 2, Ivana
More information(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);
STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend
More informationON THE STABILITY OF THE QUASI-LINEAR IMPLICIT EQUATIONS IN HILBERT SPACES
Khayyam J. Math. 5 (219), no. 1, 15-112 DOI: 1.2234/kjm.219.81222 ON THE STABILITY OF THE QUASI-LINEAR IMPLICIT EQUATIONS IN HILBERT SPACES MEHDI BENABDALLAH 1 AND MOHAMED HARIRI 2 Communicated by J. Brzdęk
More informationLYAPUNOV-RAZUMIKHIN METHOD FOR DIFFERENTIAL EQUATIONS WITH PIECEWISE CONSTANT ARGUMENT. Marat Akhmet. Duygu Aruğaslan
DISCRETE AND CONTINUOUS doi:10.3934/dcds.2009.25.457 DYNAMICAL SYSTEMS Volume 25, Number 2, October 2009 pp. 457 466 LYAPUNOV-RAZUMIKHIN METHOD FOR DIFFERENTIAL EQUATIONS WITH PIECEWISE CONSTANT ARGUMENT
More informationA maximum principle for optimal control system with endpoint constraints
Wang and Liu Journal of Inequalities and Applications 212, 212: http://www.journalofinequalitiesandapplications.com/content/212/231/ R E S E A R C H Open Access A maimum principle for optimal control system
More informationOscillatory Mixed Di erential Systems
Oscillatory Mixed Di erential Systems by José M. Ferreira Instituto Superior Técnico Department of Mathematics Av. Rovisco Pais 49- Lisboa, Portugal e-mail: jferr@math.ist.utl.pt Sra Pinelas Universidade
More informationATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS. Emerson A. M. de Abreu Alexandre N.
ATTRACTORS FOR SEMILINEAR PARABOLIC PROBLEMS WITH DIRICHLET BOUNDARY CONDITIONS IN VARYING DOMAINS Emerson A. M. de Abreu Alexandre N. Carvalho Abstract Under fairly general conditions one can prove that
More informationObserver design for a general class of triangular systems
1st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 014. Observer design for a general class of triangular systems Dimitris Boskos 1 John Tsinias Abstract The paper deals
More informationON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM
ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM OLEG ZUBELEVICH DEPARTMENT OF MATHEMATICS THE BUDGET AND TREASURY ACADEMY OF THE MINISTRY OF FINANCE OF THE RUSSIAN FEDERATION 7, ZLATOUSTINSKY MALIY PER.,
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More informationJae Gil Choi and Young Seo Park
Kangweon-Kyungki Math. Jour. 11 (23), No. 1, pp. 17 3 TRANSLATION THEOREM ON FUNCTION SPACE Jae Gil Choi and Young Seo Park Abstract. In this paper, we use a generalized Brownian motion process to define
More informationA Note on Some Properties of Local Random Attractors
Northeast. Math. J. 24(2)(2008), 163 172 A Note on Some Properties of Local Random Attractors LIU Zhen-xin ( ) (School of Mathematics, Jilin University, Changchun, 130012) SU Meng-long ( ) (Mathematics
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationPositive Solutions of Three-Point Nonlinear Second Order Boundary Value Problem
Positive Solutions of Three-Point Nonlinear Second Order Boundary Value Problem YOUSSEF N. RAFFOUL Department of Mathematics University of Dayton, Dayton, OH 45469-236 email:youssef.raffoul@notes.udayton.edu
More informationA fixed point theorem for multivalued mappings
Electronic Journal of Qualitative Theory of Differential Equations 24, No. 17, 1-1; http://www.math.u-szeged.hu/ejqtde/ A fixed point theorem for multivalued mappings Cezar AVRAMESCU Abstract A generalization
More informationA PROPERTY OF STRICTLY SINGULAR 1-1 OPERATORS
A PROPERTY OF STRICTLY SINGULAR - OPERATORS GEORGE ANDROULAKIS, PER ENFLO Abstract We prove that if T is a strictly singular - operator defined on an infinite dimensional Banach space X, then for every
More informationNew w-convergence Conditions for the Newton-Kantorovich Method
Punjab University Journal of Mathematics (ISSN 116-2526) Vol. 46(1) (214) pp. 77-86 New w-convergence Conditions for the Newton-Kantorovich Method Ioannis K. Argyros Department of Mathematicsal Sciences,
More informationRemark on Hopf Bifurcation Theorem
Remark on Hopf Bifurcation Theorem Krasnosel skii A.M., Rachinskii D.I. Institute for Information Transmission Problems Russian Academy of Sciences 19 Bolshoi Karetny lane, 101447 Moscow, Russia E-mails:
More informationTheorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).
Linear Systems Notes Lecture Proposition. A M n (R) is positive definite iff all nested minors are greater than or equal to zero. n Proof. ( ): Positive definite iff λ i >. Let det(a) = λj and H = {x D
More informationMath 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 informationON SOME CONSTANTS FOR OSCILLATION AND STABILITY OF DELAY EQUATIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 139, Number 11, November 2011, Pages 4017 4026 S 0002-9939(2011)10820-7 Article electronically published on March 28, 2011 ON SOME CONSTANTS FOR
More informationOptimal Control. Macroeconomics II SMU. Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112
Optimal Control Ömer Özak SMU Macroeconomics II Ömer Özak (SMU) Economic Growth Macroeconomics II 1 / 112 Review of the Theory of Optimal Control Section 1 Review of the Theory of Optimal Control Ömer
More informationEXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE. Leszek Gasiński
DISCRETE AND CONTINUOUS Website: www.aimsciences.org DYNAMICAL SYSTEMS SUPPLEMENT 2007 pp. 409 418 EXISTENCE RESULTS FOR QUASILINEAR HEMIVARIATIONAL INEQUALITIES AT RESONANCE Leszek Gasiński Jagiellonian
More information2 Statement of the problem and assumptions
Mathematical Notes, 25, vol. 78, no. 4, pp. 466 48. Existence Theorem for Optimal Control Problems on an Infinite Time Interval A.V. Dmitruk and N.V. Kuz kina We consider an optimal control problem on
More informationarxiv: v1 [math.ds] 30 Mar 2015
Hamburger Beiträge zur Angewandten Mathematik Nr. 2015-17, March 2015 A NOTE ON LOCAL CENTER MANIFOLDS FOR DIFFERENTIAL EQUATIONS WITH STATE-DEPENDENT DELAY EUGEN STUMPF arxiv:1503.08811v1 [math.ds] 30
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationTechnical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance
Technical Appendix for: When Promotions Meet Operations: Cross-Selling and Its Effect on Call-Center Performance In this technical appendix we provide proofs for the various results stated in the manuscript
More informationON APPROXIMATE HERMITE HADAMARD TYPE INEQUALITIES
ON APPROXIMATE HERMITE HADAMARD TYPE INEQUALITIES JUDIT MAKÓ AND ATTILA HÁZY Abstract. The main results of this paper offer sufficient conditions in order that an approximate lower Hermite Hadamard type
More informationPeriodic solutions of the perturbed symmetric Euler top
Periodic solutions of the perturbed symmetric Euler top Universitatea Babeş-Bolyai (Cluj-Napoca, Romania) abuica@math.ubbcluj.ro and Universitat de Lleida (Lleida, Spain) Plan of the talk 1 Our problem
More informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationProblem Set. Problem Set #1. Math 5322, Fall March 4, 2002 ANSWERS
Problem Set Problem Set #1 Math 5322, Fall 2001 March 4, 2002 ANSWRS i All of the problems are from Chapter 3 of the text. Problem 1. [Problem 2, page 88] If ν is a signed measure, is ν-null iff ν () 0.
More informationOn a non-autonomous stochastic Lotka-Volterra competitive system
Available online at www.isr-publications.com/jnsa J. Nonlinear Sci. Appl., 7), 399 38 Research Article Journal Homepage: www.tjnsa.com - www.isr-publications.com/jnsa On a non-autonomous stochastic Lotka-Volterra
More informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationSPACES ENDOWED WITH A GRAPH AND APPLICATIONS. Mina Dinarvand. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69, 1 (2017), 23 38 March 2017 research paper originalni nauqni rad FIXED POINT RESULTS FOR (ϕ, ψ)-contractions IN METRIC SPACES ENDOWED WITH A GRAPH AND APPLICATIONS
More informationPERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND
PERIODIC SOLUTIONS WITH NONCONSTANT SIGN IN ABEL EQUATIONS OF THE SECOND KIND JOSEP M. OLM, XAVIER ROS-OTON, AND TERE M. SEARA Abstract. The study of periodic solutions with constant sign in the Abel equation
More information