Weak Singular Optimal Controls Subject to State-Control Equality Constraints

Size: px
Start display at page:

Download "Weak Singular Optimal Controls Subject to State-Control Equality Constraints"

Transcription

1 Int. Journal of Math. Analysis, Vol. 6, 212, no. 36, Weak Singular Optimal Controls Subject to State-Control Equality Constraints Javier F Rosenblueth IIMAS UNAM, Universidad Nacional Autónoma de México Apartado Postal 2-726, México DF 1 jfrl@unam.mx Gerardo Sánchez Licea Facultad de Ciencias, Universidad Nacional Autónoma de México Departamento de Matemáticas, México DF 451 gslicea@apolo.acatlan.unam.mx Abstract The aim of this paper is to derive and illustrate a new set of second order sufficient conditions for weak minima of optimal control problems involving mixed (state-control) equality constraints. The conditions include the Legendre-Clebsch necessary condition but not its strengthened form, as well as a new version of the strengthened Weierstrass condition. Being nonsingularity crucial in the literature for sufficient results, the possible singularity of the weak optimal control in the main theorem of the paper renders sufficiency theory a much wider scope. Mathematics Subject Classification: 49K15 Keywords: Optimal control, sufficient conditions for optimality, weak minima, singularity, mixed equality constraints 1 Introduction In this paper we shall be concerned with the fixed-endpoint optimal control problem, which we label (P), of minimizing the functional I(x, u) := L(t, x(t),u(t))dt over all couples x:[,t 1 ] R n absolutely continuous and u:[,t 1 ] R m essentially bounded, satisfying the constraints

2 1798 J.F. Rosenblueth and G. Sánchez Licea a. ẋ(t) =f(t, x(t),u(t)) a.e. in [,t 1 ], b. x( )=ξ, x(t 1 )=ξ 1, c. φ(t, x(t),u(t)) = (t [,t 1 ]), where the interval T := [,t 1 ] is fixed, and we are given two points ξ and ξ 1 R n, and functions L, f and φ mapping T R n R m to R, R n and R q respectively. Second order sufficient conditions for this kind of problems are well-known in the literature (see, for example, [3, 4] and references therein). The main approaches correspond to finding a bounded solution to a matrix-valued Riccati equation, or the insertion of the original problem as an optimization problem in a Banach space, and both require that the strengthened Legendre-Clebsch condition holds. In recent papers [6, 8] we derived a new set of sufficient conditions for problems involving equality constraints only in the control functions (not mixed with the state) without assuming the strengthened Legendre-Clebsch condition. Thus, the set of sufficient conditions obtained in those papers does not require the nonsingularity of the process under consideration. In [6], the main result corresponds to sufficient conditions for a weak minimum and, in [7], we provide some examples and applications. In [8], the conditions obtained are sufficient for a strong minimum. In this paper we generalize the results of [6] for weak minima involving mixed equality constraints. The technique used is similar to that of weak minima for problems involving control equality constraints, but some crucial differences arise in the insertion of the state variables. We also illustrate the main result of the paper by providing two examples that lie beyond the scope of [3, 4], due to the singularity of the weak optimal control, and of [6], due to the presence of mixed constraints. 2 Notation and the main results For simplicity of notation, let us denote by X the space of absolutely continuous functions mapping T to R n, and set U r := L (T ; R r )(r N). Elements of X U m are called processes and a process is admissible if it satisfies the constraints. An admissible process (x, u) is called a weak minimum of (P) if, for some ɛ>, I(x, u) I(y, v) for all admissible process (y, v) not equal to (x, u) and satisfying (y, v) (x, u) <ɛ. If the inequality can be replaced by a strict inequality, the minimum is said to be a strict weak minimum. We shall assume throughout the paper that the functions L, f and φ are C 2 with respect to x and u on T R n R m.

3 Weak singular optimal controls 1799 Let us consider the Hamiltonian function H given by H(t, x, u, p, μ) := p, f(t, x, u) L(t, x, u) μ, φ(t, x, u) where p R n denotes the adjoint variable and μ R q is the multiplier associated with the mixed state-control constraint φ(t, x(t),u(t)) =. Under normality assumptions (see [2]) first order necessary conditions state that, if (x, u) is a weak minimum of (P), then there exist p X and μ U q such that ṗ(t) = Hx (t, x(t),u(t),p(t),μ(t)) a.e. in T, H u (t, x(t),u(t),p(t),μ(t)) = (t T ) where denotes transpose. Usually the term extremal corresponds to admissible processes (x, u) for which the existence of p and μ satisfying those relations can be assured. We shall find convenient to express the sufficient conditions in terms of the following function associated with the Hamiltonian. Given p X and μ U q define, for all (t, x, u) T R n R m, F (t, x, u) :=L(t, x, u) p(t),f(t, x, u) + μ(t),φ(t, x, u) ṗ(t),x. Note that F (t, x, u) = H(t, x, u, p(t),μ(t)) ṗ(t),x. Following the definition of nonsingularity given by Hestenes in [2], a process (x, u) will be called nonsingular if the determinant F uu (t, x(t),u(t)) is different from zero (t T ). With respect to F (which depends on p and μ), let J(x, u) := p(t 1 ),ξ 1 p( ),ξ + F (t, x(t),u(t))dt. Consider the first variation of J along (x, u) X U m over (y, v) X L 1 (T ; R m ) given by J ((x, u); (y, v)) := {F x (t, x(t),u(t))y(t)+f u (t, x(t),u(t))v(t)}dt and the second variation of J along (x, u) X U m over (y, v) X L 2 (T ; R m ) given by J ((x, u); (y, v)) := 2Ω(t, y(t),v(t))dt, where, for all (t, y, v) T R n R m, 2Ω(t, y, v) := y, F xx (t, x(t),u(t))y +2 y, F xu (t, x(t),u(t))v + v, F uu (t, x(t),u(t))v.

4 18 J.F. Rosenblueth and G. Sánchez Licea Finally, for all u L 1 (T ; R m ) let D(u) := ϕ(u(t))dt where ϕ(c) :=(1+ c 2 ) 1/2 1. The following theorem corresponds to the main result of the paper, a sufficiency result for a strict weak minimum of problem (P). It should be noted that the conditions include the Legendre-Clebsch necessary condition, the positivity of the second variation on a certain set of admissible variations, and a condition related to the Weierstrass excess function of the auxiliary function F. Theorem 2.1 Let (x,u ) be an admissible process. p X, μ U q, and h, ɛ > such that Suppose there exist ṗ(t) = H x(t, x (t),u (t),p(t),μ(t)) a.e. in T, H u (t, x (t),u (t),p(t),μ(t)) = (t T ), and the following conditions are satisfied: i. F uu (t, x (t),u (t)) (a.e. in T ). ii. J ((x,u ); (y, v)) > for all nonnull (y, v) X L 2 (T ; R m ) satisfying a. ẏ(t) = f x (t, x (t),u (t))y(t) +f u (t, x (t),u (t))v(t) a.e. in T, and y( )=y(t 1 )=; b. φ x (t, x (t),u (t))y(t)+φ u (t, x (t),u (t))v(t) =a.e. in T. iii. For all admissible processes (x, u) satisfying (x, u) (x,u ) <ɛ, E(t, x(t),u (t),u(t))dt hd(u u ) where E denotes the Weierstrass excess function with respect to F, E(t, x, u, v) :=F (t, x, v) F (t, x, u) F u (t, x, u)(v u). Then there exist ρ, δ > such that, for all admissible processes (x, u) satisfying (x, u) (x,u ) <ρ, I(x, u) I(x,u )+δd(u u ). In particular, (x,u ) is a strict weak minimum of (P). Let us end this section by deriving a simple consequence of Theorem 2.1. It corresponds to a sufficiency result which, in certain cases, may be much easier to verify than the previous one, though both the continuity of the control and the nonsingularity of the process under consideration are assumed.

5 Weak singular optimal controls 181 Corollary 2.2 Let (x,u ) be an admissible process with u Suppose there exist p X and μ U q such that continuous. ṗ(t) = H x(t, x (t),u (t),p(t),μ(t)) a.e. in T, H u (t, x (t),u (t),p(t),μ(t)) = (t T ), and the following conditions are satisfied: i. F uu (t, x (t),u (t)) > (t T ). ii. J ((x,u ); (y, v)) > for all nonnull (y, v) X L 2 (T ; R m ) satisfying a. ẏ(t) = f x (t, x (t),u (t))y(t) +f u (t, x (t),u (t))v(t) a.e. in T, and y( )=y(t 1 )=; b. φ x (t, x (t),u (t))y(t)+φ u (t, x (t),u (t))v(t) =a.e. in T. Then there exist ρ, δ > such that, for all admissible processes (x, u) satisfying (x, u) (x,u ) <ρ, I(x, u) I(x,u )+δd(u u ). In particular, (x,u ) is a strict weak minimum of (P). Proof: Define a restricted tube of radius ɛ> centred on (x,u )as T 1 ((x,u ); ɛ) :={(t, y, v) T R n R m : x (t) y <ɛ, u (t) v <ɛ}. By 2.2(i) and since u is continuous, there exist h, ɛ > such that c, F uu (t, x, u)c h c 2 (c R m, (t, x, u) T 1 ((x,u ); ɛ)). Hence, for (t, x, u, v) with (t, x, u) and (t, x, v) int 1 ((x,u ); ɛ), E(t, x, u, v) = (1 λ) v u, F uu (t, x, u + λ[v u])(v u) dλ 1 2 h v u 2 hϕ(v u). The conclusion follows by Theorem Proof of Theorem 2.1 In this section we shall prove Theorem 2.1. We first state an auxiliary result which is an immediate consequence of Lemma 3.1 of [5] and on which the proof of the theorem is strongly based.

6 182 J.F. Rosenblueth and G. Sánchez Licea Lemma 3.1 Let {u q } be a sequence in L 1 (T ; R m ), u L 1 (T ; R m ), and suppose that lim D(u q u )= and d q := [2D(u q u )] 1/2 > (q N). For all q N and t T define w q (t) := [ 1+ 1 ] 1/2, 2 ϕ(u q(t) u (t)) vq (t) := u q(t) u (t). d q Then the following hold: a. For some v L 2 (T ; R m ) and some subsequence of {u q }, again denoted by {u q }, {v q } converges weakly to v in L 1 (T ; R m ). b. Let A q L (T ; R n n ) and B q L (T ; R n m ) be matrix functions for which there exist constants m,m 1 > such that A q m, B q m 1 (q N), and denote by y q the solution of the initial value problem ẏ(t) =A q (t)y(t)+b q (t)v q (t) (a.e. in T ), y( )=. Then there exist σ L 2 (T ; R n ) and a subsequence of {y q }, again denoted by {y q }, such that {ẏ q } converges weakly in L 1 (T ; R n ) to σ. Moreover, if y (t) := t σ (s)ds (t T ), then y q (t) y (t) uniformly on T. c. Suppose that w q (t) 1 uniformly on T, let R q ( ) be an m m real matrix-valued measurable function on T, R ( ) L (T ; R m m ), and assume also that R q (t) R (t) uniformly on T, and R (t) a.e. in T. Then, for some subsequence of {u q }, again denoted by {u q }, lim inf R q (t)v q (t),v q (t) dt R (t)v (t),v (t) dt. Proof of Theorem 2.1: Assume that, for all ρ, δ >, there exists an admissible process (x, u) with (x, u) (x,u ) <ρsuch that J(x, u) <J(x,u )+δd(u u ). (1) We shall show that this contradicts 2.1(ii) and the statement will follow, since J(x, u) = I(x, u) for all admissible processes (x, u). Let z := (x,u ). Note that, for all admissible processes z =(x, u), J(z) =J(z )+J (z ; z z )+K(z)+Ẽ(z) (2) where Ẽ(x, u) := E(t, x(t),u (t),u(t))dt,

7 Weak singular optimal controls 183 K(x, u) := {M(t, x(t)) + u(t) u (t),n(t, x(t)) }dt, and the functions M and N are given by M(t, y) :=F (t, y, u (t)) F (t, x (t),u (t)) F x (t, x (t),u (t))(y x (t)), N(t, y) :=Fu (t, y, u (t)) Fu(t, x (t),u (t)). By Taylor s theorem we have M(t, y) = 1 2 y x (t),p(t, y)(y x (t)), N(t, y) =Q(t, y)(y x (t)), where P (t, y) :=2 (1 λ)f xx (t, x (t)+λ(y x (t)),u (t))dλ, Q(t, y) := F ux (t, x (t)+λ(y x (t)),u (t))dλ. Now, by (1), for all q N there exists an admissible process z q := (x q,u q ) such that z q z < 1 q, J(z q) J(z ) < 1 q D(u q u ). (3) Since z q is admissible, observe that the last inequality implies that u q (t) u (t) on a set of positive measure and so D(u q u ) > (q N). Since z q z asq, it follows that D(u q u ), q. Define d q, w q, v q as in Lemma 3.1 and, for all t T and q N, set y q (t) := x q(t) x (t) d q. By Lemma 3.1(a) there exist v L 2 (T ; R m ) and some subsequence of {z q } (we do not relabel) such that {v q } converges weakly in L 1 (T ; R m )tov. By Taylor s theorem, for all q N we have ẏ q (t) =A q (t)y q (t)+b q (t)v q (t) (a.e. in T ) where A q (t) = B q (t) = f x (t, x (t)+λ[x q (t) x (t)],u (t))dλ, f u (t, x q (t),u q (t)+λ[u (t) u q (t)])dλ. By continuity of f x and f u there exist m,m 1 > such that A q m and B q m 1 (q N). By Lemma 3.1(b) there exist σ L 2 (T ; R n ) and some subsequence of {z q } (we do not relabel) such that, if y (t) := t σ (s)ds (t T ), then y q (t) y (t) uniformly on T.

8 184 J.F. Rosenblueth and G. Sánchez Licea The theorem will be proved if we show that J (z ;(y,v )), (y,v ) (, ), y ( )=y (t 1 ) =, and the following two relations hold: ẏ (t) =f x (t, x (t),u (t))y (t)+f u (t, x (t),u (t))v (t) a.e. in T, φ x (t, x (t),u (t))y (t)+φ u (t, x (t),u (t))v (t) = a.e. in T. To prove it observe that, since z q is admissible, the fact that y ( ) = y (t 1 ) = follows by Lemma 3.1(b). Now, by definition of the functional K, for all q N, K(z q ) d 2 q = In view of Lemma 3.1(b), { M(t, xq (t)) + d 2 q v q (t), N(t, x q(t)) d q } dt. M(t, x q (t)) lim = 1 d 2 q 2 y (t),f xx (t, x (t),u (t))y (t), N(t, x q (t)) lim = F ux (t, x (t),u (t))y (t) d q both uniformly on T and, since {v q } converges weakly to v in L 1 (T ; R m ), 1 2 J K(z q ) (z ;(y,v )) = lim + 1 v d 2 (t),f uu (t, x (t),u (t))v (t) dt. (4) q 2 Now, by Taylor s theorem, for all t T and q N, where Clearly, 1 E(t, x q (t),u (t),u q (t)) = 1 2 v q(t),r q (t)v q (t) d 2 q R q (t) :=2 (1 λ)f uu (t, x q (t),u (t)+λ[u q (t) u (t)])dλ. lim R q(t) =R (t) :=F uu (t, x (t),u (t)) uniformly on T. Since z q z asq, it follows that w q (t) 1 uniformly on T and by 2.1(i), R (t) a.e. in T. By Lemma 3.1(c), there is a subsequence of {z q } (we do not relabel) such that lim inf Ẽ(z q ) d 2 q 1 2 v (t),f uu (t, x (t),u (t))v (t) dt. (5)

9 Weak singular optimal controls 185 Since ṗ(t) = Hx (t, x (t),u (t),p(t),μ(t)) a.e. in T, H u (t, x (t),u (t),p(t),μ(t)) = (t T ), it follows that J (z ;(y, v)) = for all (y, v) X L 1 (T ; R m ). With this in mind, (2), (3), (4) and (5), 1 2 J K(z q ) Ẽ(z q ) (z ;(y,v )) lim + lim inf d 2 q d 2 q = lim inf J(z q ) J(z ) d 2 q. In addition, if (y,v )=(, ), then lim K(z q )/d 2 q = and, by 2.1(iii), h 2 lim inf Ẽ(z q ) d 2 q contradicting the positivity of h. Now, observe that y q (t) y (t), A q (t) f x (t, x (t),u (t)), B q (t) f u (t, x (t),u (t)) all uniformly on T, and {v q } converges weakly to v in L 1 (T ; R m ). Therefore {ẏ q } converges weakly in L 1 (T ; R n )to f x (t, x (t),u (t))y (t)+f u (t, x (t),u (t))v (t). By Lemma 3.1(b), {ẏ q } converges weakly in L 1 (T ; R n )toσ =ẏ. Hence, ẏ (t) =f x (t, x (t),u (t))y (t)+f u (t, x (t),u (t))v (t) a.e. in T. Finally to show that φ x (t, x (t),u (t))y (t)+φ u (t, x (t),u (t))v (t) = a.e. in T, for all q N, t T and λ [, 1], set G q (t; λ) :=φ(t, x (t)+λ[x q (t) x (t)],u (t)+λ[u q (t) u (t)]). Clearly, for all q N and t T, By Taylor s theorem, if we set and G q (t;)=g q (t;1)=. φ x (t; λ; q) :=φ x (t, x (t)+λ[x q (t) x (t)],u (t)+λ[u q (t) u (t)]) φ u (t; λ; q) :=φ u (t, x (t)+λ[x q (t) x (t)],u (t)+λ[u q (t) u (t)]),

10 186 J.F. Rosenblueth and G. Sánchez Licea it follows that for all q N and t T, ( ) ( ) φ x (t; λ; q)dλ y q (t)+ φ u (t; λ; q)dλ v q (t) =. Since z q z, y q (t) y (t) uniformly on T and {v q } converges weakly to v in L 1 (T ; R m ), t {φ x (s, x (s),u (s))y (s)+φ u (s, x (s),u (s))v (s)}ds = for all t T. Consequently, φ x (t, x (t),u (t))y (t)+φ u (t, x (t),u (t))v (t) = a.e. in T and this completes the proof. 4 Examples In this section we provide two examples of a fixed-endpoint optimal control problem with equality state-control constraints for which an application of Theorem 2.1 shows that the singular extremal under consideration is in fact a strict weak minimum. Example 4.1 Consider the problem of minimizing I(x, u) = ( t/2)x(t)[u(t) 1]dt subject to a. ẋ(t) =( t/2)x(t)u(t) tx(t) a.e. in [, 1]. b. x() = 1, x(1) = e 2/3. c. x(t)u(t)e u(t) ax 2 (t)u(t) =(t [, 1]). For this case, n = m = q =1,T =[, 1], ξ =1,ξ 1 = e 2/3,<a<1, L(t, x, u) =( t/2)x[u 1], f(t, x, u) =( t/2)xu tx, φ(t, x, u) =xue u ax 2 u. Let z := (x,u ) (e ( 2/3)t3/2, ) which is admissible. We have H(t, x, u, p, μ) =(p t/2)xu p tx +( t/2)x[1 u] μ[xue u ax 2 u], H x (t, x, u, p, μ) =(p t/2)u p t +( t/2)[1 u] μ[ue u 2axu], H u (t, x, u, p, μ) =(p t/2)x ( t/2)x μ[xe u xue u ax 2 ].

11 Weak singular optimal controls 187 Let p 1/2 and t μ(t) := 4[1 ae ( 2/3)t3/2 ] (t T ). As one readily verifies (x,u,p,μ) satisfies the first order conditions of Theorem 2.1. Moreover, the function F is given by and, hence, F (t, x, u) = txu 4 t[xue u ax 2 u] 4[1 ae ( 2/3)t3/2 ] F uu (t, x (t),u (t)) = te ( 2/3)t 3/2 2[1 ae ( 2/3)t3/2 ] (t T ) implying that z is singular and the first condition of Theorem 2.1 holds. Also, observe that f x (t, x (t),u (t)) = te ( 2/3)t 3/2 t and f u (t, x (t),u (t)) = (t T ) 2 and hence (y, v) satisfies (a) and (b) of the second condition if and only if ẏ(t) = te ( 2/3)t 3/2 ty(t)+ v(t) a.e. in T, 2 y() = y(1) =, and e ( 2/3)t3/2 [1 e ( 2/3)t3/2 ]v(t) = a.e. in T. Consequently, there are no nonnull admissible variations and the second condition is fulfilled. Now observe that, by Taylor s theorem, ( ) E(t, x,,u)= (1 λ)f uu (t, x, λu)dλ u 2, and txg(u) F uu (t, x, u) = 4[1 ae ( 2/3)t3/2 ], where g(u) :=2e u ue u. Hence, since g() = 2 > and x (t) =e ( 2/3)t3/2 > (t T ), there are two positive numbers δ, ɛ such that g(w) δ for all w ( ɛ, ɛ), and T (x ; ɛ) :={(t, x) T R : x x (t) <ɛ}

12 188 J.F. Rosenblueth and G. Sánchez Licea is contained in Q δ := {(t, x) R 2 t, x δ}. Thusif(t, x, u) belongs to then T 1 ((x,u ); ɛ) ={(t, x, u) T (x ; ɛ) R : u u (t) <ɛ}, ( E(t, x,,u)= (1 λ) ) txg(λu) 4[1 ae ( 2/3)t3/2 ] dλ u 2 δ2 t 8[1 a] u2. Let us now suppose that the third condition does not hold. Then, for all q N, there exists an admissible process z q := (x q,u q ) such that z q z < 1/q, z q z <ɛ, and E(t, x q (t),u (t),u q (t))dt < 1 q D(u q). Since for all q N and almost all t T, E(t, x q (t),u (t),u q (t)) δ2 t 8[1 a] u2 q (t), for all q N we have δ 2 v t q 2(t) 8[1 a] wq 2(t)dt δ 2 tv 2 8[1 a] q (t)dt < 1 2q where w q (t) :=[ ϕ(u q(t))] 1/2, v q (t) := u q(t) d q, and d q := [2D(u q )] 1/2. From the above we conclude that lim Now note that, for all q N, v 2 t q (t) dt =. wq(t) 2 By (6) and integrating by parts, for all q N, v t q 2(t) ( t wq 2(t)dt + vq 2(s) ) dt wq 2(s)ds 2 t vq(t) 2 =1. (6) wq 2 (t)dt =1. (7) Since h q (t) := t v2 q (s)/w2 q (s)ds 1(t T, q N), by (7), for all q N, v 2 q(t) 1 t wq 2(t)dt + h 2 q(t) dt 2 t 1. (8)

13 Weak singular optimal controls 189 By (8), there exist some subsequence (we do not relabel) and some h in L 2 (T ; R) such that the sequence α q (t) :=h q (t)/ 4 t converges weakly to h in L 2 (T ; R), that is, h q (t) 1 lim 4 g(t)dt = h (t)g(t)dt for all g L 2 (T ; R). (9) Set g(t) :=1/ 4 t if <t 1, and g() :=. Since 1 tv 2 q (t)/wq 2 (t)dt as q, by using g given above, it follows by (7) and (9) that h (t) 4 t dt =2. (1) We claim that h (t) = 1/ 4 t a.e. in T. To prove it, observe first that h (t) a.e. in T for, otherwise, there would exist S subset of T such that m(s) > and h (t) < (t S). Defining g(t) :=1ift S and g(t) :=if t T \ S and replacing g in (9) by g, we obtain h q (t) lim 4 dt = h (t)dt < S t S which is a contradiction. Now, observe that h (t) dt = h 2 t (t)dt 2 Once again by (9) and (1), h (t) 1 4 dt + h 2 (t)dt = lim h q (t) h (t) 1 h (t) 4 dt t 4 dt =2, t dt 1 = h 2 (t)dt 2. and this proves our claim. By (9), it is readily seen that lim h q (t)g(t)dt = g(t)dt for all g L 2 (T ; R), that is, h q converges weakly to 1 in L 2 (T ; R). Since lim sup h q 2 1 2, by Proposition 3.32 of [1, p. 78], we have h q 1 2 as. Then by (6), = lim t v 2 q(s) w 2 q (s)ds 1 2 dt = lim t and hence for some subsequence (we do not relabel), lim t vq(s) 2 = pointwisely a.e. in T. wq 2 (s)ds v 2 q(s) w 2 q (s)ds 2 dt,

14 181 J.F. Rosenblueth and G. Sánchez Licea Consequently, lim vq(s) 2 = wq 2 (s)ds which contradicts (6). Therefore also the third condition holds and, by Theorem 2.1, (x,u ) is a strict weak minimum of problem (P). Example 4.2 Consider the problem of minimizing I(x, u) = {u 2 1(t) (1/2)u 2 2(t) (1/8)x(t)u 1 (t)}dt subject to a. ẋ(t) = sinh u 1 (t) (1/8)x 3 (t) a.e. in [, 1]. b. x() = x(1) =. c. x(t)u 2 1(t)+u 2 2(t) =4(t [, 1]). Note that, for this case, n = q =1,m =2,T =[, 1], ξ = ξ 1 =, L(t, x, u) =u 2 1 (1/2)u2 2 (1/8)xu 1, f(t, x, u) = sinh u 1 (1/8)x 3, φ(t, x, u) =xu u Let z := (x,u ) (,, 2). Clearly z is admissible, and we have H(t, x, u, p, μ) =p sinh u 1 (1/8)px 3 u 2 1 +(1/2)u2 2 +(1/8)xu 1 μxu 2 1 μu2 2 +4μ, H x (t, x, u, p, μ) = (3/8)px 2 +(1/8)u 1 μu 2 1, H u (t, x, u, p, μ) =(pcosh u 1 2u 1 +(1/8)x 2μxu 1,u 2 2μu 2 ). Therefore (x,u,p,μ) with (p, μ) (, 1/2) satisfies the first order conditions of Theorem 2.1. Now, the function F is given by Since, for all t T, F (t, x, u) =u 2 1 (1/8)xu 1 +(1/2)xu F uu (t, x (t),u (t)) = ( ) 2, it follows that F uu (t, x (t),u (t))h, h =2h 2 1, and so (x,u ) is singular and the first condition of Theorem 2.1 holds. Also, note that f x (t, x (t),u (t)) = and f u (t, x (t),u (t)) = (1, ) (t T ) so that, for any (y, v) satisfying (a) and (b) of the second condition, we must have ẏ(t) =v 1 (t) a.e. in T, y() = y(1) = and v 2 (t) = a.e. in T.

15 Weak singular optimal controls 1811 It follows that, for the second variation, for all (y, v) satisfying (a) and (b) of the second condition, J (z ;(y, v)) = {2v 2 1 (t) (1/4)y(t)v 1(t)}dt = 2v 2 1 (t)dt and the second condition of the theorem is satisfied. Suppose now that the third condition does not hold. Then, for all q N, there exists an admissible process z q := (x q,u q ) such that z q z < 1/q and E(t, x q (t),u (t),u q (t))dt < 1 q D(u q u ). Since E(t, x(t),u (t),u(t)) = E(t, x(t),, 2,u 1 (t),u 2 (t)) = u 2 1 (t)+(1/2)x(t)u2 1 (t), for all q N, u 2 1q(t)dt x q (t)u 2 1q(t)dt < 1 q D(u q u ). (11) Since x q (t)u 2 1q(t)+u 2 2q(t) =4,u 1q (t), u 2q (t) 2 both uniformly on T, ϕ(c) c 2 /2 and (2 + ϕ(c))ϕ(c) = c 2, we have the existence of some m > such that, for all q N, D(u q u ) 1 u q (t) u (t) 2 dt = 1 u 2 1q 2 2 (t)dt + 1 (u 2q (t) 2) 2 dt 2 = 1 u 2 1q 2 (t)dt + 1 x 2 q (t)u4 1q (t) 2 (2 + u 2q (t)) dt u 2 2 1q (t)dt u 2 1q(t) 2+ϕ(u 1q ) 2+ϕ(u 1q (t)) dt = 2+ϕ(u 1q ) ϕ(u 1q (t))dt m D(u 1q ). By (11) and since x q (t)u 2 1q (t) +u2 2q (t) =4(t T, q N), it follows that u 1q implying that D(u 1q ) > (q N). Consequently, D(u q u ) D(u 1q ) m (q N). Once again by (11), v1q(t) 2 w1q(t) dt + 1 w 2 2 1q(t)x 2 q (t) v2 1q(t) w1q(t) < m 2 2q (12)

16 1812 J.F. Rosenblueth and G. Sánchez Licea where, for all t T and q N, w 1q (t) :=[ ϕ(u 1q(t))] 1/2, v 1q (t) := u 1q(t) d 1q and d 1q := [2D(u 1q )] 1/2. But since for all q N, v1q 2 (t) dt =1, w1q(t) 2 w 1q (t) 1, x q (t) both uniformly on T, by letting q in (12), one obtains 1 which is a contradiction. Therefore also the third condition holds and, by Theorem 2.1, (x,u ) is a strict weak minimum of problem (P). References [1] Brezis H (21) Functional Analysis, Sobolev Spaces and Partial Differential Equations, Universitext, Springer [2] Hestenes MR (1966) Calculus of Variations and Optimal Control Theory, John Wiley, New York [3] Maurer H, Pickenhain S (1995) Second order sufficient conditions for control problems with mixed control-state constraints, Journal of Optimization Theory and Applications, 86: [4] Milyutin AA, Osmolovskiǐ NP (1998) Calculus of Variations and Optimal Control, Translations of Mathematical Monographs 18, Providence, Rhode Island: American Mathematical Society [5] Rosenblueth JF, Sánchez Licea G (21) A direct sufficiency proof for a weak minimum in optimal control, Applied Mathematical Sciences, 4: [6] Rosenblueth JF, Sánchez Licea G (211) Sufficiency for singular controls with equality constraints, Proceedings of the 13th IASTED International Conference on Intelligent Systems & Control, Cambridge, England, doi: /P , [7] Rosenblueth JF, Sánchez Licea G (in press) Singular weak optimal controls, Control & Intelligent Systems [8] Rosenblueth JF, Sánchez Licea G (in press) Sufficiency and singularity in optimal control, IMA Journal of Mathematical Control and Information Received: March, 212

Sufficiency for Essentially Bounded Controls Which Do Not Satisfy the Strengthened Legendre-Clebsch Condition

Sufficiency for Essentially Bounded Controls Which Do Not Satisfy the Strengthened Legendre-Clebsch Condition Applied Mathematical Sciences, Vol. 12, 218, no. 27, 1297-1315 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/ams.218.81137 Sufficiency for Essentially Bounded Controls Which Do Not Satisfy the Strengthened

More information

Sufficiency in optimal control without. the strengthened condition of Legendre

Sufficiency in optimal control without. the strengthened condition of Legendre Journal of Applied Mathematics & Bioinformatics, vol.1, no.1, 2011, 1-20 IN: 1792-6602 (print), 1792-6939 (online) c International cientific Press, 2011 ufficiency in optimal control without the strengthened

More information

Strong and Weak Augmentability in Calculus of Variations

Strong and Weak Augmentability in Calculus of Variations Strong and Weak Augmentability in Calculus of Variations JAVIER F ROSENBLUETH National Autonomous University of Mexico Applied Mathematics and Systems Research Institute Apartado Postal 20-126, Mexico

More information

Critical Cones for Regular Controls with Inequality Constraints

Critical Cones for Regular Controls with Inequality Constraints International Journal of Mathematical Analysis Vol. 12, 2018, no. 10, 439-468 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2018.8856 Critical Cones for Regular Controls with Inequality Constraints

More information

Conjugate Journey in Optimal Control

Conjugate Journey in Optimal Control International Mathematical Forum, 2, 2007, no. 14, 633-674 Conjugate Journey in Optimal Control Javier F Rosenblueth IIMAS-UNAM, Apartado Postal 20-726 México DF 01000, México jfrl@servidor.unam.mx Abstract

More information

LECTURES ON OPTIMAL CONTROL THEORY

LECTURES ON OPTIMAL CONTROL THEORY LECTURES ON OPTIMAL CONTROL THEORY Terje Sund May 24, 2016 CONTENTS 1. INTRODUCTION 2. FUNCTIONS OF SEVERAL VARIABLES 3. CALCULUS OF VARIATIONS 4. OPTIMAL CONTROL THEORY 1 INTRODUCTION In the theory of

More information

Existence Of Solution For Third-Order m-point Boundary Value Problem

Existence Of Solution For Third-Order m-point Boundary Value Problem Applied Mathematics E-Notes, 1(21), 268-274 c ISSN 167-251 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Existence Of Solution For Third-Order m-point Boundary Value Problem Jian-Ping

More information

2 Statement of the problem and assumptions

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

On perturbations in the leading coefficient matrix of time-varying index-1 DAEs

On perturbations in the leading coefficient matrix of time-varying index-1 DAEs On perturbations in the leading coefficient matrix of time-varying index-1 DAEs Institute of Mathematics, Ilmenau University of Technology Darmstadt, March 27, 2012 Institute of Mathematics, Ilmenau University

More information

Continuous Functions on Metric Spaces

Continuous Functions on Metric Spaces Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0

More information

Nonlinear Control Systems

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

Lu u. µ (, ) the equation. has the non-zero solution

Lu u. µ (, ) the equation. has the non-zero solution MODULE 18 Topics: Eigenvalues and eigenvectors for differential operators In analogy to the matrix eigenvalue problem Ax = λx we shall consider the eigenvalue problem Lu = µu where µ is a real or complex

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

Existence of homoclinic solutions for Duffing type differential equation with deviating argument

Existence of homoclinic solutions for Duffing type differential equation with deviating argument 2014 9 «28 «3 Sept. 2014 Communication on Applied Mathematics and Computation Vol.28 No.3 DOI 10.3969/j.issn.1006-6330.2014.03.007 Existence of homoclinic solutions for Duffing type differential equation

More information

BANACH SPACES IN WHICH EVERY p-weakly SUMMABLE SEQUENCE LIES IN THE RANGE OF A VECTOR MEASURE

BANACH SPACES IN WHICH EVERY p-weakly SUMMABLE SEQUENCE LIES IN THE RANGE OF A VECTOR MEASURE PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 124, Number 7, July 1996 BANACH SPACES IN WHICH EVERY p-weakly SUMMABLE SEQUENCE LIES IN THE RANGE OF A VECTOR MEASURE C. PIÑEIRO (Communicated by

More information

Relative Controllability of Fractional Dynamical Systems with Multiple Delays in Control

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

2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.

2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers. Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following

More information

Half of Final Exam Name: Practice Problems October 28, 2014

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

Partial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations

Partial Differential Equations, 2nd Edition, L.C.Evans The Calculus of Variations Partial Differential Equations, 2nd Edition, L.C.Evans Chapter 8 The Calculus of Variations Yung-Hsiang Huang 2018.03.25 Notation: denotes a bounded smooth, open subset of R n. All given functions are

More information

MATH 220: MIDTERM OCTOBER 29, 2015

MATH 220: MIDTERM OCTOBER 29, 2015 MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and

More information

NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS

NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS Fixed Point Theory, Volume 9, No. 1, 28, 3-16 http://www.math.ubbcluj.ro/ nodeacj/sfptcj.html NONTRIVIAL SOLUTIONS TO INTEGRAL AND DIFFERENTIAL EQUATIONS GIOVANNI ANELLO Department of Mathematics University

More information

Formula Sheet for Optimal Control

Formula Sheet for Optimal Control Formula Sheet for Optimal Control Division of Optimization and Systems Theory Royal Institute of Technology 144 Stockholm, Sweden 23 December 1, 29 1 Dynamic Programming 11 Discrete Dynamic Programming

More information

EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS

EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC SCHRÖDINGER EQUATIONS Electronic Journal of Differential Equations, Vol. 017 (017), No. 15, pp. 1 7. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO ASYMPTOTICALLY PERIODIC

More information

A note on linear differential equations with periodic coefficients.

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

GENERATORS WITH INTERIOR DEGENERACY ON SPACES OF L 2 TYPE

GENERATORS WITH INTERIOR DEGENERACY ON SPACES OF L 2 TYPE Electronic Journal of Differential Equations, Vol. 22 (22), No. 89, pp. 3. ISSN: 72-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu GENERATORS WITH INTERIOR

More information

Ordinary Differential Equation Theory

Ordinary Differential Equation Theory Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit

More information

INITIAL AND BOUNDARY VALUE PROBLEMS FOR NONCONVEX VALUED MULTIVALUED FUNCTIONAL DIFFERENTIAL EQUATIONS

INITIAL AND BOUNDARY VALUE PROBLEMS FOR NONCONVEX VALUED MULTIVALUED FUNCTIONAL DIFFERENTIAL EQUATIONS Applied Mathematics and Stochastic Analysis, 6:2 23, 9-2. Printed in the USA c 23 by North Atlantic Science Publishing Company INITIAL AND BOUNDARY VALUE PROBLEMS FOR NONCONVEX VALUED MULTIVALUED FUNCTIONAL

More information

Abstract In this paper, we consider bang-bang property for a kind of timevarying. time optimal control problem of null controllable heat equation.

Abstract In this paper, we consider bang-bang property for a kind of timevarying. time optimal control problem of null controllable heat equation. JOTA manuscript No. (will be inserted by the editor) The Bang-Bang Property of Time-Varying Optimal Time Control for Null Controllable Heat Equation Dong-Hui Yang Bao-Zhu Guo Weihua Gui Chunhua Yang Received:

More information

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation Advances in Dynamical Systems and Applications ISSN 973-532, Volume 6, Number 2, pp. 24 254 (2 http://campus.mst.edu/adsa Nontrivial Solutions for Boundary Value Problems of Nonlinear Differential Equation

More information

NOTES ON CALCULUS OF VARIATIONS. September 13, 2012

NOTES ON CALCULUS OF VARIATIONS. September 13, 2012 NOTES ON CALCULUS OF VARIATIONS JON JOHNSEN September 13, 212 1. The basic problem In Calculus of Variations one is given a fixed C 2 -function F (t, x, u), where F is defined for t [, t 1 ] and x, u R,

More information

Second Order Sufficient Conditions for Optimal Control Problems with Non-unique Minimizers

Second Order Sufficient Conditions for Optimal Control Problems with Non-unique Minimizers 2 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 2 WeA22. Second Order Sufficient Conditions for Optimal Control Problems with Non-unique Minimizers Christos Gavriel

More information

ATTRACTORS 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. 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 information

arxiv: v1 [math.oc] 22 Sep 2016

arxiv: 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 information

The local equicontinuity of a maximal monotone operator

The local equicontinuity of a maximal monotone operator arxiv:1410.3328v2 [math.fa] 3 Nov 2014 The local equicontinuity of a maximal monotone operator M.D. Voisei Abstract The local equicontinuity of an operator T : X X with proper Fitzpatrick function ϕ T

More information

Notes on uniform convergence

Notes on uniform convergence Notes on uniform convergence Erik Wahlén erik.wahlen@math.lu.se January 17, 2012 1 Numerical sequences We begin by recalling some properties of numerical sequences. By a numerical sequence we simply mean

More information

Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.

Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological

More information

BIHARMONIC WAVE MAPS INTO SPHERES

BIHARMONIC WAVE MAPS INTO SPHERES BIHARMONIC WAVE MAPS INTO SPHERES SEBASTIAN HERR, TOBIAS LAMM, AND ROLAND SCHNAUBELT Abstract. A global weak solution of the biharmonic wave map equation in the energy space for spherical targets is constructed.

More information

Advanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x

Advanced Calculus Math 127B, Winter 2005 Solutions: Final. nx2 1 + n 2 x, g n(x) = n2 x . Define f n, g n : [, ] R by f n (x) = Advanced Calculus Math 27B, Winter 25 Solutions: Final nx2 + n 2 x, g n(x) = n2 x 2 + n 2 x. 2 Show that the sequences (f n ), (g n ) converge pointwise on [, ],

More information

Maximal monotone operators are selfdual vector fields and vice-versa

Maximal monotone operators are selfdual vector fields and vice-versa Maximal monotone operators are selfdual vector fields and vice-versa Nassif Ghoussoub Department of Mathematics, University of British Columbia, Vancouver BC Canada V6T 1Z2 nassif@math.ubc.ca February

More information

WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)

WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2) WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H

More information

Some asymptotic properties of solutions for Burgers equation in L p (R)

Some asymptotic properties of solutions for Burgers equation in L p (R) ARMA manuscript No. (will be inserted by the editor) Some asymptotic properties of solutions for Burgers equation in L p (R) PAULO R. ZINGANO Abstract We discuss time asymptotic properties of solutions

More information

Nonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:

Nonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage: Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm

More information

Initial value problems for singular and nonsmooth second order differential inclusions

Initial value problems for singular and nonsmooth second order differential inclusions Initial value problems for singular and nonsmooth second order differential inclusions Daniel C. Biles, J. Ángel Cid, and Rodrigo López Pouso Department of Mathematics, Western Kentucky University, Bowling

More information

ALMOST PERIODIC SOLUTIONS OF HIGHER ORDER DIFFERENTIAL EQUATIONS ON HILBERT SPACES

ALMOST PERIODIC SOLUTIONS OF HIGHER ORDER DIFFERENTIAL EQUATIONS ON HILBERT SPACES Electronic Journal of Differential Equations, Vol. 21(21, No. 72, pp. 1 12. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu ALMOST PERIODIC SOLUTIONS

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

A Concise Course on Stochastic Partial Differential Equations

A Concise Course on Stochastic Partial Differential Equations A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original

More information

An Integral-type Constraint Qualification for Optimal Control Problems with State Constraints

An Integral-type Constraint Qualification for Optimal Control Problems with State Constraints An Integral-type Constraint Qualification for Optimal Control Problems with State Constraints S. Lopes, F. A. C. C. Fontes and M. d. R. de Pinho Officina Mathematica report, April 4, 27 Abstract Standard

More information

EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction

EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.

More information

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics

Lecture Notes in Mathematics. Arkansas Tech University Department of Mathematics Lecture Notes in Mathematics Arkansas Tech University Department of Mathematics Introductory Notes in Ordinary Differential Equations for Physical Sciences and Engineering Marcel B. Finan c All Rights

More information

u(0) = u 0, u(1) = u 1. To prove what we want we introduce a new function, where c = sup x [0,1] a(x) and ɛ 0:

u(0) = u 0, u(1) = u 1. To prove what we want we introduce a new function, where c = sup x [0,1] a(x) and ɛ 0: 6. Maximum Principles Goal: gives properties of a solution of a PDE without solving it. For the two-point boundary problem we shall show that the extreme values of the solution are attained on the boundary.

More information

Lecture 9. Systems of Two First Order Linear ODEs

Lecture 9. Systems of Two First Order Linear ODEs Math 245 - Mathematics of Physics and Engineering I Lecture 9. Systems of Two First Order Linear ODEs January 30, 2012 Konstantin Zuev (USC) Math 245, Lecture 9 January 30, 2012 1 / 15 Agenda General Form

More information

Math 273, Final Exam Solutions

Math 273, Final Exam Solutions Math 273, Final Exam Solutions 1. Find the solution of the differential equation y = y +e x that satisfies the condition y(x) 0 as x +. SOLUTION: y = y H + y P where y H = ce x is a solution of the homogeneous

More information

Chapter 2 Convex Analysis

Chapter 2 Convex Analysis Chapter 2 Convex Analysis The theory of nonsmooth analysis is based on convex analysis. Thus, we start this chapter by giving basic concepts and results of convexity (for further readings see also [202,

More information

A TWO PARAMETERS AMBROSETTI PRODI PROBLEM*

A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* PORTUGALIAE MATHEMATICA Vol. 53 Fasc. 3 1996 A TWO PARAMETERS AMBROSETTI PRODI PROBLEM* C. De Coster** and P. Habets 1 Introduction The study of the Ambrosetti Prodi problem has started with the paper

More information

OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS

OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

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

PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS

PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS Annales Academiæ Scientiarum Fennicæ Mathematica Volumen 28, 2003, 207 222 PERTURBATION THEORY FOR NONLINEAR DIRICHLET PROBLEMS Fumi-Yuki Maeda and Takayori Ono Hiroshima Institute of Technology, Miyake,

More information

HOMOCLINIC SOLUTIONS FOR SECOND-ORDER NON-AUTONOMOUS HAMILTONIAN SYSTEMS WITHOUT GLOBAL AMBROSETTI-RABINOWITZ CONDITIONS

HOMOCLINIC SOLUTIONS FOR SECOND-ORDER NON-AUTONOMOUS HAMILTONIAN SYSTEMS WITHOUT GLOBAL AMBROSETTI-RABINOWITZ CONDITIONS Electronic Journal of Differential Equations, Vol. 010010, No. 9, pp. 1 10. ISSN: 107-6691. UL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu HOMOCLINIC SOLUTIONS FO

More information

Fundamental Solutions of Stokes and Oseen Problem in Two Spatial Dimensions

Fundamental Solutions of Stokes and Oseen Problem in Two Spatial Dimensions J. math. fluid mech. c 6 Birkhäuser Verlag, Basel DOI.7/s-5-9-z Journal of Mathematical Fluid Mechanics Fundamental Solutions of Stokes and Oseen Problem in Two Spatial Dimensions Ronald B. Guenther and

More information

POSITIVE SOLUTIONS FOR NONLOCAL SEMIPOSITONE FIRST ORDER BOUNDARY VALUE PROBLEMS

POSITIVE SOLUTIONS FOR NONLOCAL SEMIPOSITONE FIRST ORDER BOUNDARY VALUE PROBLEMS Dynamic Systems and Applications 5 6 439-45 POSITIVE SOLUTIONS FOR NONLOCAL SEMIPOSITONE FIRST ORDER BOUNDARY VALUE PROBLEMS ERBIL ÇETIN AND FATMA SERAP TOPAL Department of Mathematics, Ege University,

More information

MATH 5640: Fourier Series

MATH 5640: Fourier Series MATH 564: Fourier Series Hung Phan, UMass Lowell September, 8 Power Series A power series in the variable x is a series of the form a + a x + a x + = where the coefficients a, a,... are real or complex

More information

Layer structures for the solutions to the perturbed simple pendulum problems

Layer structures for the solutions to the perturbed simple pendulum problems Layer structures for the solutions to the perturbed simple pendulum problems Tetsutaro Shibata Applied Mathematics Research Group, Graduate School of Engineering, Hiroshima University, Higashi-Hiroshima,

More information

On the minimum of certain functional related to the Schrödinger equation

On the minimum of certain functional related to the Schrödinger equation Electronic Journal of Qualitative Theory of Differential Equations 2013, No. 8, 1-21; http://www.math.u-szeged.hu/ejqtde/ On the minimum of certain functional related to the Schrödinger equation Artūras

More information

Chapter 4. Inverse Function Theorem. 4.1 The Inverse Function Theorem

Chapter 4. Inverse Function Theorem. 4.1 The Inverse Function Theorem Chapter 4 Inverse Function Theorem d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d d dd d d d d This chapter

More information

Nonlinear equations. Norms for R n. Convergence orders for iterative methods

Nonlinear equations. Norms for R n. Convergence orders for iterative methods Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector

More information

Oscillation by Impulses for a Second-Order Delay Differential Equation

Oscillation by Impulses for a Second-Order Delay Differential Equation PERGAMON Computers and Mathematics with Applications 0 (2006 0 www.elsevier.com/locate/camwa Oscillation by Impulses for a Second-Order Delay Differential Equation L. P. Gimenes and M. Federson Departamento

More information

NONHOMOGENEOUS ELLIPTIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT AND WEIGHT

NONHOMOGENEOUS ELLIPTIC EQUATIONS INVOLVING CRITICAL SOBOLEV EXPONENT AND WEIGHT Electronic Journal of Differential Equations, Vol. 016 (016), No. 08, pp. 1 1. ISSN: 107-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu NONHOMOGENEOUS ELLIPTIC

More information

Joint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018

Joint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018 EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,

More information

PERIODIC SOLUTIONS OF THE FORCED PENDULUM : CLASSICAL VS RELATIVISTIC

PERIODIC SOLUTIONS OF THE FORCED PENDULUM : CLASSICAL VS RELATIVISTIC LE MATEMATICHE Vol. LXV 21) Fasc. II, pp. 97 17 doi: 1.4418/21.65.2.11 PERIODIC SOLUTIONS OF THE FORCED PENDULUM : CLASSICAL VS RELATIVISTIC JEAN MAWHIN The paper surveys and compares some results on the

More information

B. Appendix B. Topological vector spaces

B. Appendix B. Topological vector spaces B.1 B. Appendix B. Topological vector spaces B.1. Fréchet spaces. In this appendix we go through the definition of Fréchet spaces and their inductive limits, such as they are used for definitions of function

More information

Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions

Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Existence and Uniqueness Results for Nonlinear Implicit Fractional Differential Equations with Boundary Conditions Mouffak Benchohra a,b 1 and Jamal E. Lazreg a, a Laboratory of Mathematics, University

More information

I. The space C(K) Let K be a compact metric space, with metric d K. Let B(K) be the space of real valued bounded functions on K with the sup-norm

I. The space C(K) Let K be a compact metric space, with metric d K. Let B(K) be the space of real valued bounded functions on K with the sup-norm I. The space C(K) Let K be a compact metric space, with metric d K. Let B(K) be the space of real valued bounded functions on K with the sup-norm Proposition : B(K) is complete. f = sup f(x) x K Proof.

More information

CMAT Centro de Matemática da Universidade do Minho

CMAT Centro de Matemática da Universidade do Minho Universidade do Minho CMAT Centro de Matemática da Universidade do Minho Campus de Gualtar 471-57 Braga Portugal wwwcmatuminhopt Universidade do Minho Escola de Ciências Centro de Matemática Blow-up and

More information

Nonlinear Control Lecture 5: Stability Analysis II

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

PUTNAM PROBLEMS DIFFERENTIAL EQUATIONS. First Order Equations. p(x)dx)) = q(x) exp(

PUTNAM PROBLEMS DIFFERENTIAL EQUATIONS. First Order Equations. p(x)dx)) = q(x) exp( PUTNAM PROBLEMS DIFFERENTIAL EQUATIONS First Order Equations 1. Linear y + p(x)y = q(x) Muliply through by the integrating factor exp( p(x)) to obtain (y exp( p(x))) = q(x) exp( p(x)). 2. Separation of

More information

Applied Differential Equation. November 30, 2012

Applied Differential Equation. November 30, 2012 Applied Differential Equation November 3, Contents 5 System of First Order Linear Equations 5 Introduction and Review of matrices 5 Systems of Linear Algebraic Equations, Linear Independence, Eigenvalues,

More information

Exercises: Brunn, Minkowski and convex pie

Exercises: Brunn, Minkowski and convex pie Lecture 1 Exercises: Brunn, Minkowski and convex pie Consider the following problem: 1.1 Playing a convex pie Consider the following game with two players - you and me. I am cooking a pie, which should

More information

************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan

************************************* Applied Analysis I - (Advanced PDE I) (Math 940, Fall 2014) Baisheng Yan ************************************* Applied Analysis I - (Advanced PDE I) (Math 94, Fall 214) by Baisheng Yan Department of Mathematics Michigan State University yan@math.msu.edu Contents Chapter 1.

More information

Stability and Instability for Dynamic Equations on Time Scales

Stability and Instability for Dynamic Equations on Time Scales PERGAMON Computers and Mathematics with Applications 0 (2005) 1 0 www.elsevier.com/locate/camwa Stability and Instability for Dynamic Equations on Time Scales J. Hoffacker Department of Mathematical Sciences,

More information

Theory of Ordinary Differential Equations

Theory of Ordinary Differential Equations Theory of Ordinary Differential Equations Existence, Uniqueness and Stability Jishan Hu and Wei-Ping Li Department of Mathematics The Hong Kong University of Science and Technology ii Copyright c 24 by

More information

Math 699 Reading Course, Spring 2007 Rouben Rostamian Homogenization of Differential Equations May 11, 2007 by Alen Agheksanterian

Math 699 Reading Course, Spring 2007 Rouben Rostamian Homogenization of Differential Equations May 11, 2007 by Alen Agheksanterian . Introduction Math 699 Reading Course, Spring 007 Rouben Rostamian Homogenization of ifferential Equations May, 007 by Alen Agheksanterian In this brief note, we will use several results from functional

More information

NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS. Shouchuan Hu Nikolas S. Papageorgiou. 1. Introduction

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

Chapter 4 Optimal Control Problems in Infinite Dimensional Function Space

Chapter 4 Optimal Control Problems in Infinite Dimensional Function Space Chapter 4 Optimal Control Problems in Infinite Dimensional Function Space 4.1 Introduction In this chapter, we will consider optimal control problems in function space where we will restrict ourselves

More information

The complex periodic problem for a Riccati equation

The complex periodic problem for a Riccati equation The complex periodic problem for a Riccati equation Rafael Ortega Departamento de Matemática Aplicada Facultad de Ciencias Universidad de Granada, 1871 Granada, Spain rortega@ugr.es Dedicated to Jean Mawhin,

More information

MATH107 Vectors and Matrices

MATH107 Vectors and Matrices School of Mathematics, KSU 20/11/16 Vector valued functions Let D be a set of real numbers D R. A vector-valued functions r with domain D is a correspondence that assigns to each number t in D exactly

More information

An invariance result for Hammersley s process with sources and sinks

An invariance result for Hammersley s process with sources and sinks An invariance result for Hammersley s process with sources and sinks Piet Groeneboom Delft University of Technology, Vrije Universiteit, Amsterdam, and University of Washington, Seattle March 31, 26 Abstract

More information

Some approximation theorems in Math 522. I. Approximations of continuous functions by smooth functions

Some approximation theorems in Math 522. I. Approximations of continuous functions by smooth functions Some approximation theorems in Math 522 I. Approximations of continuous functions by smooth functions The goal is to approximate continuous functions vanishing at ± by smooth functions. In a previous homewor

More information

Convex Optimization Theory. Chapter 5 Exercises and Solutions: Extended Version

Convex Optimization Theory. Chapter 5 Exercises and Solutions: Extended Version Convex Optimization Theory Chapter 5 Exercises and Solutions: Extended Version Dimitri P. Bertsekas Massachusetts Institute of Technology Athena Scientific, Belmont, Massachusetts http://www.athenasc.com

More information

Examples include: (a) the Lorenz system for climate and weather modeling (b) the Hodgkin-Huxley system for neuron modeling

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

PERIODIC SOLUTIONS FOR A KIND OF LIÉNARD-TYPE p(t)-laplacian EQUATION. R. Ayazoglu (Mashiyev), I. Ekincioglu, G. Alisoy

PERIODIC SOLUTIONS FOR A KIND OF LIÉNARD-TYPE p(t)-laplacian EQUATION. R. Ayazoglu (Mashiyev), I. Ekincioglu, G. Alisoy Acta Universitatis Apulensis ISSN: 1582-5329 http://www.uab.ro/auajournal/ No. 47/216 pp. 61-72 doi: 1.17114/j.aua.216.47.5 PERIODIC SOLUTIONS FOR A KIND OF LIÉNARD-TYPE pt)-laplacian EQUATION R. Ayazoglu

More information

HeadMedia Interaction in Magnetic Recording

HeadMedia Interaction in Magnetic Recording Journal of Differential Equations 171, 443461 (2001) doi:10.1006jdeq.2000.3844, available online at http:www.idealibrary.com on HeadMedia Interaction in Magnetic Recording Avner Friedman Departament of

More information

Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation

Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Symmetry, Integrability and Geometry: Methods and Applications Vol. (5), Paper 3, 9 pages Transient Phenomena in Quantum Bound States Subjected to a Sudden Perturbation Marcos MOSHINSKY and Emerson SADURNÍ

More information

THE NEARLY ADDITIVE MAPS

THE NEARLY ADDITIVE MAPS Bull. Korean Math. Soc. 46 (009), No., pp. 199 07 DOI 10.4134/BKMS.009.46..199 THE NEARLY ADDITIVE MAPS Esmaeeil Ansari-Piri and Nasrin Eghbali Abstract. This note is a verification on the relations between

More information

GEOMETRY OF THE PROJECTIVE UNITARY GROUP OF A C -ALGEBRA

GEOMETRY OF THE PROJECTIVE UNITARY GROUP OF A C -ALGEBRA REVISTA DE LA UNIÓN MATEMÁTICA ARGENTINA Vol. 58, No. 2, 2017, Pages 319 329 Published online: June 5, 2017 GEOMETRY OF THE PROJECTIVE UNITARY GROUP OF A C -ALGEBRA ESTEBAN ANDRUCHOW Abstract. Let A be

More information

SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS

SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS J. Appl. Math. & Computing Vol. 23(2007), No. 1-2, pp. 243-256 Website: http://jamc.net SWEEP METHOD IN ANALYSIS OPTIMAL CONTROL FOR RENDEZ-VOUS PROBLEMS MIHAI POPESCU Abstract. This paper deals with determining

More information

OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL

OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL OPTIMALITY CONDITIONS AND ERROR ANALYSIS OF SEMILINEAR ELLIPTIC CONTROL PROBLEMS WITH L 1 COST FUNCTIONAL EDUARDO CASAS, ROLAND HERZOG, AND GERD WACHSMUTH Abstract. Semilinear elliptic optimal control

More information

ON WEAKLY NONLINEAR BACKWARD PARABOLIC PROBLEM

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

Periodic solutions of weakly coupled superlinear systems

Periodic solutions of weakly coupled superlinear systems Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize

More information

EXISTENCE RESULTS FOR OPERATOR EQUATIONS INVOLVING DUALITY MAPPINGS VIA THE MOUNTAIN PASS THEOREM

EXISTENCE RESULTS FOR OPERATOR EQUATIONS INVOLVING DUALITY MAPPINGS VIA THE MOUNTAIN PASS THEOREM EXISTENCE RESULTS FOR OPERATOR EQUATIONS INVOLVING DUALITY MAPPINGS VIA THE MOUNTAIN PASS THEOREM JENICĂ CRÎNGANU We derive existence results for operator equations having the form J ϕu = N f u, by using

More information

van Rooij, Schikhof: A Second Course on Real Functions

van Rooij, Schikhof: A Second Course on Real Functions vanrooijschikhof.tex April 25, 2018 van Rooij, Schikhof: A Second Course on Real Functions Notes from [vrs]. Introduction A monotone function is Riemann integrable. A continuous function is Riemann integrable.

More information