Semiconcavity and optimal control: an intrinsic approach

Size: px
Start display at page:

Download "Semiconcavity and optimal control: an intrinsic approach"

Transcription

1 Semiconcavity and optimal control: an intrinsic approach Peter R. Wolenski joint work with Piermarco Cannarsa and Francesco Marino Louisiana State University SADCO summer school, London September 5-9, 2011 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

2 Outline 1 Optimal Control problems Optimal control Value functions and semiconcavity Differential Inclusions (DI) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

3 Outline 1 Optimal Control problems Optimal control Value functions and semiconcavity Differential Inclusions (DI) 2 Smooth parameterizations? Necessary conditions Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

4 Outline 1 Optimal Control problems Optimal control Value functions and semiconcavity Differential Inclusions (DI) 2 Smooth parameterizations? Necessary conditions 3 New (DI) assumptions Examples Consequences Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

5 Outline 1 Optimal Control problems Optimal control Value functions and semiconcavity Differential Inclusions (DI) 2 Smooth parameterizations? Necessary conditions 3 New (DI) assumptions Examples Consequences 4 New idea and results A replacement for a priori estimates Semiconcavity results with (DI) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

6 Optimal control problems Control Dynamics: (CD) ẋ(s) = f ( x(s), u(s) ) a.e. s [t, T ] u(s) U a.e. s [t, T ] x(t) = x, where f : R n R n R n is continuous in (x, u) and Lipschitz in x, the admissible control set U R m is compact, and u : [t, T ] R m is measurable. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

7 Optimal control problems Control Dynamics: (CD) ẋ(s) = f ( x(s), u(s) ) a.e. s [t, T ] u(s) U a.e. s [t, T ] x(t) = x, where f : R n R n R n is continuous in (x, u) and Lipschitz in x, the admissible control set U R m is compact, and u : [t, T ] R m is measurable. Two classic problems: 1. MinTime: Given a closed target set S R n, the problem is min (T t) over ( x( ), u( ) ) satisfying (CD) and x(t ) S. The optimal value T (x) is the minimum time function. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

8 Optimal control problems Control Dynamics: (CD) ẋ(s) = f ( x(s), u(s) ) a.e. s [t, T ] u(s) U a.e. s [t, T ] x(t) = x, where f : R n R n R n is continuous in (x, u) and Lipschitz in x, the admissible control set U R m is compact, and u : [t, T ] R m is measurable. Two classic problems: 1. MinTime: Given a closed target set S R n, the problem is min (T t) over ( x( ), u( ) ) satisfying (CD) and x(t ) S. The optimal value T (x) is the minimum time function. 2. Mayer problem: Given endpoint cost l : R n R, the problem is min l ( x(t ) ) over ( x( ), u( ) ) satisfying (CD). The optimal value V (t, x) is the value function. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

9 Regularity of value functions - SemiConCavity (SCC) A natural regularity property for T ( ) and V (, ) is the property of being semiconcave. A locally Lipschitz function g : R m R is semiconcave provided there exists k > 0 so that 1[ ] g(x + z) + g(x z) g(x) k z 2 2 x, z R n. (the three point property) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

10 Regularity of value functions - SemiConCavity (SCC) A natural regularity property for T ( ) and V (, ) is the property of being semiconcave. A locally Lipschitz function g : R m R is semiconcave provided there exists k > 0 so that 1[ ] g(x + z) + g(x z) g(x) k z 2 2 x, z R n. (the three point property) A semiconcave function

11 Regularity of value functions - SemiConCavity (SCC) A natural regularity property for T ( ) and V (, ) is the property of being semiconcave. A locally Lipschitz function g : R m R is semiconcave provided there exists k > 0 so that 1[ ] g(x + z) + g(x z) g(x) k z 2 2 x, z R n. (the three point property) A semiconcave function x z x x + z

12 Regularity of value functions - SemiConCavity (SCC) A natural regularity property for T ( ) and V (, ) is the property of being semiconcave. A locally Lipschitz function g : R m R is semiconcave provided there exists k > 0 so that 1[ ] g(x + z) + g(x z) g(x) k z 2 2 x, z R n. (the three point property) A semiconcave function LHS x z x x + z

13 Regularity of value functions - SemiConCavity (SCC) A natural regularity property for T ( ) and V (, ) is the property of being semiconcave. A locally Lipschitz function g : R m R is semiconcave provided there exists k > 0 so that 1[ ] g(x + z) + g(x z) g(x) k z 2 2 x, z R n. (the three point property) A semiconcave function LHS x z x x + z Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

14 A geometric description of SCC A Lipschitz function g : R n R is (SCC) if and only if σ > 0 with g(x) = inf { q(x) : q(x) = σx 2 + bx + c, g(x) q(x) } = inf { q(x, α) : α A }, where (x, α) q(x, α) is C 1+ in x and continuous in (x, α). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

15 A geometric description of SCC A Lipschitz function g : R n R is (SCC) if and only if σ > 0 with g(x) = inf { q(x) : q(x) = σx 2 + bx + c, g(x) q(x) } = inf { q(x, α) : α A }, where (x, α) q(x, α) is C 1+ in x and continuous in (x, α).

16 A geometric description of SCC A Lipschitz function g : R n R is (SCC) if and only if σ > 0 with g(x) = inf { q(x) : q(x) = σx 2 + bx + c, g(x) q(x) } = inf { q(x, α) : α A }, where (x, α) q(x, α) is C 1+ in x and continuous in (x, α).

17 A geometric description of SCC A Lipschitz function g : R n R is (SCC) if and only if σ > 0 with g(x) = inf { q(x) : q(x) = σx 2 + bx + c, g(x) q(x) } = inf { q(x, α) : α A }, where (x, α) q(x, α) is C 1+ in x and continuous in (x, α).

18 A geometric description of SCC A Lipschitz function g : R n R is (SCC) if and only if σ > 0 with g(x) = inf { q(x) : q(x) = σx 2 + bx + c, g(x) q(x) } = inf { q(x, α) : α A }, where (x, α) q(x, α) is C 1+ in x and continuous in (x, α). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

19 Previous results yielding (SCC) There is considerable literature on value functions being (SCC). Most relevant here: Cannarsa, Frankowska, Sinestrari, McEneaney. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

20 Previous results yielding (SCC) There is considerable literature on value functions being (SCC). Most relevant here: Cannarsa, Frankowska, Sinestrari, McEneaney. Basic idea with (CD): (Illustration with Min Time) We seek an upper bound (by k z 2 ) of T (x + z) + T (x z) 2T (x). Take an optimal solution starting from x and use it to construct feasible solutions from x ± z that will yield the appropriate estimates. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

21 Previous results yielding (SCC) There is considerable literature on value functions being (SCC). Most relevant here: Cannarsa, Frankowska, Sinestrari, McEneaney. Basic idea with (CD): (Illustration with Min Time) We seek an upper bound (by k z 2 ) of T (x + z) + T (x z) 2T (x). Take an optimal solution starting from x and use it to construct feasible solutions from x ± z that will yield the appropriate estimates. Assume that x f (x, u) is C 1+ and take ( x( ), ū( ) ) optimal. Use a priori estimates on the ODEs { ẋ ± (s) = f ( x ± (s), ū(s) ) a.e. s [t, T ] (ODE ± ) x ± (t) = x ± z. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

22 x S

23 S { x(t) = f ( x(t), ū(t) ) x(0) = x Optimal x

24 S x + z { x(t) = f ( x(t), ū(t) ) x(0) = x Optimal x x z

25 S { ẋ+(t) = f ( x +(t), ū(t) ) x +(0) = x + z Same control ū( ) x + z { x(t) = f ( x(t), ū(t) ) x(0) = x Optimal x { ẋ (t) = f ( x (t), ū(t) ) x (0) = x z Same control ū( ) x z Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

26 This proof depends on the parameterization A priori estimates rely on the specific parameterization (x, u) f (x, u) of the dynamics, but the value functions T ( ) and V (, ) do not. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

27 This proof depends on the parameterization A priori estimates rely on the specific parameterization (x, u) f (x, u) of the dynamics, but the value functions T ( ) and V (, ) do not. Important: To obtain second order (SCC) estimates, one definitely requires more than mere Lipschitz of the map x f (x, u). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

28 This proof depends on the parameterization A priori estimates rely on the specific parameterization (x, u) f (x, u) of the dynamics, but the value functions T ( ) and V (, ) do not. Important: To obtain second order (SCC) estimates, one definitely requires more than mere Lipschitz of the map x f (x, u). The question is What? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

29 This proof depends on the parameterization A priori estimates rely on the specific parameterization (x, u) f (x, u) of the dynamics, but the value functions T ( ) and V (, ) do not. Important: To obtain second order (SCC) estimates, one definitely requires more than mere Lipschitz of the map x f (x, u). The question is What? Recall the previous work assumed x f (x, u) is C 1+. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

30 (Very) simple example Note that T ( ) and V (, ) depend only on the trajectories x( ) and NOT in the parameterization of the admissible velocity set: F (x) := { f (x, u) : u U }. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

31 (Very) simple example Note that T ( ) and V (, ) depend only on the trajectories x( ) and NOT in the parameterization of the admissible velocity set: F (x) := { f (x, u) : u U }. Note the[ admissible ] velocity multifunction F : R R given by F (x) = x, x can be parameterized two ways: F (x) = { {x u : u 1} { x u : u 1} The former is smoothly parameterized whereas the latter is not. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

32 (Very) simple example Note that T ( ) and V (, ) depend only on the trajectories x( ) and NOT in the parameterization of the admissible velocity set: F (x) := { f (x, u) : u U }. Note the[ admissible ] velocity multifunction F : R R given by F (x) = x, x can be parameterized two ways: F (x) = { {x u : u 1} { x u : u 1} The former is smoothly parameterized whereas the latter is not. Trajectories coincide, but theorems only apply to the former. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

33 Differential Inclusions and Filippov s Lemma The set of solutions to the Differential Inclusion { ẋ(s) F ( x(s) ) a.e. s [t, T ] (DI) x(t) = x does not depend on the particular parameterization. This is (essentially) the content of the well-known Filippov s Lemma. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

34 Differential Inclusions and Filippov s Lemma The set of solutions to the Differential Inclusion { ẋ(s) F ( x(s) ) a.e. s [t, T ] (DI) x(t) = x does not depend on the particular parameterization. This is (essentially) the content of the well-known Filippov s Lemma. Natural question: For which F is the value function semiconcave? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

35 Differential Inclusions and Filippov s Lemma The set of solutions to the Differential Inclusion { ẋ(s) F ( x(s) ) a.e. s [t, T ] (DI) x(t) = x does not depend on the particular parameterization. This is (essentially) the content of the well-known Filippov s Lemma. Natural question: For which F is the value function semiconcave? A satisfactory answer should be given in terms of F, or equivalently, by the Hamiltonian H : R n R n R defined by: H(x, p) = v, p. sup v F (x) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

36 Equivalence of F and H We assume throughout that F : R n R n satisfies the following Standard Hypotheses: 1) F (x) is nonempty, convex, and compact x, (SH) + 2) F is Lipschitz on bounded sets w.r.t. Hausdorff metric; 3) r > 0 so that max{ v : v F (x)} r(1 + x ). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

37 Equivalence of F and H We assume throughout that F : R n R n satisfies the following Standard Hypotheses: 1) F (x) is nonempty, convex, and compact x, (SH) + 2) F is Lipschitz on bounded sets w.r.t. Hausdorff metric; 3) r > 0 so that max{ v : v F (x)} r(1 + x ). Such assumptions on F give way to equivalent conditions on H because v F (x) v, p H(x, p) p R n Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

38 Equivalence of F and H We assume throughout that F : R n R n satisfies the following Standard Hypotheses: 1) F (x) is nonempty, convex, and compact x, (SH) + 2) F is Lipschitz on bounded sets w.r.t. Hausdorff metric; 3) r > 0 so that max{ v : v F (x)} r(1 + x ). Such assumptions on F give way to equivalent conditions on H because v F (x) v, p H(x, p) p R n 1) x R n, H(x, p) is finite and convex, positively homogeneous in p; (SH) + 2) M > 0, k > 0 so that x, y M, p R n, H(x, p) H(y, p) k p x y ; 3) r > 0 so that H(x, p) r p (1 + x ) x, p R n. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

39 Smooth parameterizations? Perhaps one can characterize those multifunctions that have a smooth parameterization: Question: Given F : R n R n, when does there exist f : R n U R n that is C 1 in the first coordinate and satisfies F (x) := { f (x, u) : u U }? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

40 Smooth parameterizations? Perhaps one can characterize those multifunctions that have a smooth parameterization: Question: Given F : R n R n, when does there exist f : R n U R n that is C 1 in the first coordinate and satisfies F (x) := { f (x, u) : u U }? This seems virtually impossible to answer. Worse: Even sufficient conditions for smooth selections seems intractable: A simpler (?) question: Under what conditions on F does there exist a C 1 function f : R n R n so that f (x) F (x) x? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

41 Necessary conditions for smooth parameterizations We can say when a smooth parameterization does NOT exist! Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

42 Necessary conditions for smooth parameterizations We can say when a smooth parameterization does NOT exist! Suppose F : R n R n is smoothly parameterized and H has the form { } f H(x, p) = sup (x, u), p : u U, where f : R n U R n has f and f x both continuous in (x, u). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

43 Necessary conditions for smooth parameterizations We can say when a smooth parameterization does NOT exist! Suppose F : R n R n is smoothly parameterized and H has the form { } f H(x, p) = sup (x, u), p : u U, where f : R n U R n has f and f x Then for 0 p R n, we have (H1) The map x H(x, p) is semiconvex; both continuous in (x, u). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

44 Necessary conditions for smooth parameterizations We can say when a smooth parameterization does NOT exist! Suppose F : R n R n is smoothly parameterized and H has the form { } f H(x, p) = sup (x, u), p : u U, where f : R n U R n has f and f x Then for 0 p R n, we have (H1) The map x H(x, p) is semiconvex; and (NC) If H(x, p) = H(x, p), then both continuous in (x, u). x H(x, p) = x H(x, p). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

45 Necessary conditions for smooth parameterizations We can say when a smooth parameterization does NOT exist! Suppose F : R n R n is smoothly parameterized and H has the form { } f H(x, p) = sup (x, u), p : u U, where f : R n U R n has f and f x Then for 0 p R n, we have (H1) The map x H(x, p) is semiconvex; and (NC) If H(x, p) = H(x, p), then both continuous in (x, u). x H(x, p) = x H(x, p). If the assumption of a smooth parameterization is replaced by the existence of a smooth selection, then the conclusion of (NC) is x H(x, p) x H(x, p). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

46 Illustration of (NC) with n = 1 F (x) H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

47 Illustration of (NC) with n = 1 F (x) F (x) = [ h(x), H(x) ] H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

48 Illustration of (NC) with n = 1 F (x) F (x) = [ h(x), H(x) ] H(x) F (x 3) F (x) x 1 x x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

49 Illustration of (NC) with n = 1 F (x) no smooth parameterization H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 x 1 : No smooth parameterization since x H(x) not semiconvex. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

50 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 x 2, x 3 : Smooth parameterizations are possible between x 1 and x 4. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

51 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

52 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

53 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

54 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 x 2, x 3 : Smooth parameterizations are possible between x 1 and x 4. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

55 Illustration of (NC) with n = 1 F (x) smooth parameterizations H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 x 2, x 3 : Smooth parameterizations are possible between x 1 and x 4. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

56 Illustration of (NC) with n = 1 F (x) H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

57 Illustration of (NC) with n = 1 F (x) no smooth selection H(x) x 1 x 2 h(x) Graph of F (x) x x 3 x 4 x 4 : No smooth selection. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

58 Proof of (NC) Proof. Note that the assumption means that H(x, p) = H(x, p) sup f (x, u), p = sup v, p = inf v, p = inf f (x, u), p u U v F (x) v F (x) u U That is, the assumption is that every u U both minimizes and maximizes the quantity f (x, u), p. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

59 Proof of (NC) Proof. Note that the assumption means that H(x, p) = H(x, p) sup f (x, u), p = sup v, p = inf v, p = inf f (x, u), p u U v F (x) v F (x) u U That is, the assumption is that every u U both minimizes and maximizes the quantity f (x, u), p. By a theorem on nonsmooth differentiation of max functions, one has { } x H(x, p) = co x f (x, u) p : u U and { } x H(x, p) = co x f (x, u) p : u U, from which (NC) follows: Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

60 Proof of (NC) Proof. Note that the assumption means that H(x, p) = H(x, p) sup f (x, u), p = sup v, p = inf v, p = inf f (x, u), p u U v F (x) v F (x) u U That is, the assumption is that every u U both minimizes and maximizes the quantity f (x, u), p. By a theorem on nonsmooth differentiation of max functions, one has { } x H(x, p) = co x f (x, u) p : u U and { } x H(x, p) = co x f (x, u) p : u U, from which (NC) follows: x H(x, p) = x H(x, p). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

61 New DI assumptions We abandon looking for smooth parameterizations, and introduce: { 1) x H(x, p) is semiconvex, and (H) 2) The gradient p H(x, p) exists and is locally Lipschitz in x. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

62 New DI assumptions We abandon looking for smooth parameterizations, and introduce: { 1) x H(x, p) is semiconvex, and (H) 2) The gradient p H(x, p) exists and is locally Lipschitz in x. Class of examples: One can generate a class of examples that satisfy (H) but do not satisfy (NC), and therefore could not have a C 1 parameterization: Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

63 New DI assumptions We abandon looking for smooth parameterizations, and introduce: { 1) x H(x, p) is semiconvex, and (H) 2) The gradient p H(x, p) exists and is locally Lipschitz in x. Class of examples: One can generate a class of examples that satisfy (H) but do not satisfy (NC), and therefore could not have a C 1 parameterization: n = 1: Let F (x) := [ h(x), H(x) ] where h( ) and H( ) are semiconvex. These always satisfy (H), and could satisfy (NC) only if h(x) = H(x) h(x) = H(x). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

64 New DI assumptions We abandon looking for smooth parameterizations, and introduce: { 1) x H(x, p) is semiconvex, and (H) 2) The gradient p H(x, p) exists and is locally Lipschitz in x. Class of examples: One can generate a class of examples that satisfy (H) but do not satisfy (NC), and therefore could not have a C 1 parameterization: n = 1: Let F (x) := [ h(x), H(x) ] where h( ) and H( ) are semiconvex. These always satisfy (H), and could satisfy (NC) only if h(x) = H(x) h(x) = H(x). n > 1: Let F (x) := f (x) + r(x)b where f ( ) is C 2 and r : R n [0, ) is semiconvex. Then H(x, p) = f (x), p + r(x) p, and so (H) is satisfied. Then (NC) is satisfied only if r(x) = 0 implies r(x) = r(x). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

65 Consequences, part I Consequences of (H1): The semiconvexity of x H(x, p) implies x,p H(x, p) x H ( x, p ) p H ( x, p ). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

66 Consequences, part I Consequences of (H1): The semiconvexity of x H(x, p) implies x,p H(x, p) x H ( x, p ) p H ( x, p ). The significance of this result is in utilizing a nonsmooth maximum principle (Clarke 1975): Suppose x( ) is optimal in one of the classical problems with (DI) dynamics. Then there exists an adjoint arc p( ) for which ( p(s), x(s) ) x,p H ( x(s), p(s) ) (plus transversality conditions). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

67 Consequences, part I Consequences of (H1): The semiconvexity of x H(x, p) implies x,p H(x, p) x H ( x, p ) p H ( x, p ). The significance of this result is in utilizing a nonsmooth maximum principle (Clarke 1975): Suppose x( ) is optimal in one of the classical problems with (DI) dynamics. Then there exists an adjoint arc p( ) for which ( p(s), x(s) ) x,p H ( x(s), p(s) ) (plus transversality conditions). Thus the dynamics of a Hamiltonian arc ( x( ), p( ) ) splits into a much more usable form: p(s) x H ( x(s), p(s) ) and x(s) p H ( x(s), p(s) ) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

68 Consequences, part II Consequences of (H2): That the gradient p H(x, p) exists means that the argmax of v, p is unique - we denote it by f p (x) F (x). sup v F (x) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

69 Consequences, part II Consequences of (H2): That the gradient p H(x, p) exists means that the argmax of v, p is unique - we denote it by f p (x) F (x). sup v F (x) We have also assumed x f p (x) is Lipschitz. If p( ) is any continuous function, then the ODE (ODE) x { ẋ(t) = f p(t) ( x(t) ) x(0) = x, satisfies standard Carathéodory-type assumptions a.e. t [0, T ] Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

70 Consequences, part II Consequences of (H2): That the gradient p H(x, p) exists means that the argmax of v, p is unique - we denote it by f p (x) F (x). sup v F (x) We have also assumed x f p (x) is Lipschitz. If p( ) is any continuous function, then the ODE (ODE) x { ẋ(t) = f p(t) ( x(t) ) x(0) = x, a.e. t [0, T ] satisfies standard Carathéodory-type assumptions, and has the properties: A unique solution x( ; x) of (ODE) x exists on [0, ); Each solution x( ; x) is a solution of (DI); The function x x(t; x) is locally Lipschitz; If ( x( ), p( ) ) is a Hamiltonian arc, then x( ) satisfies (ODE) with p( ) = p( ). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

71 A replacement for a priori estimates S x

72 A replacement for a priori estimates S { x(t) = f p(t) ( x(t) ) x(0) = x Optimal with adjoint p( ) x

73 A replacement for a priori estimates S x + z { x(t) = f p(t) ( x(t) ) x(0) = x Optimal with adjoint p( ) x x z

74 A replacement for a priori estimates S { ẋ+(t) = f p(t) ( x+(t) ) x +(0) = x + z Same adjoint arc p( ) x + z { x(t) = f p(t) ( x(t) ) x(0) = x Optimal with adjoint p( ) x { ẋ (t) = f p(t) ( x (t) ) x (0) = x z Same adjoint arc p( ) x z Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

75 Summary of new results We assume F satisfies (SH) + and (H). Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

76 Summary of new results We assume F satisfies (SH) + and (H). Theorem Consider the minimum time problem. Suppose the target S is compact and satisfies the Petrov condition and the Interior Sphere Property. Then there exists ρ > 0 so that T ( ) is semiconcave on each convex subset of S + ρb \ S. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

77 Summary of new results We assume F satisfies (SH) + and (H). Theorem Consider the minimum time problem. Suppose the target S is compact and satisfies the Petrov condition and the Interior Sphere Property. Then there exists ρ > 0 so that T ( ) is semiconcave on each convex subset of S + ρb \ S. Theorem Consider the Mayer problem. Assume the endpoint cost function l( ) is semiconcave. Then the value function V (, ) is locally semiconcave on (, T ] R n. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

78 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

79 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. We offered an alternative; It is not yet clear how generic these assumptions are. For example, could they be invoked in an auxiliary problem that approximates an original one? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

80 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. We offered an alternative; It is not yet clear how generic these assumptions are. For example, could they be invoked in an auxiliary problem that approximates an original one? Local versions? Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

81 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. We offered an alternative; It is not yet clear how generic these assumptions are. For example, could they be invoked in an auxiliary problem that approximates an original one? Local versions? Perhaps a large segment of optimal control can be cast in this way.(?) Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

82 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. We offered an alternative; It is not yet clear how generic these assumptions are. For example, could they be invoked in an auxiliary problem that approximates an original one? Local versions? Perhaps a large segment of optimal control can be cast in this way.(?) Students: you re urged and invited to join the fun. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

83 Conclusions and future work General Philosophy: Assumptions should not be required only in the proof. Here, smooth parameterizations of control systems should not be assumed if the conclusion does not reflect them. We offered an alternative; It is not yet clear how generic these assumptions are. For example, could they be invoked in an auxiliary problem that approximates an original one? Local versions? Perhaps a large segment of optimal control can be cast in this way.(?) Students: you re urged and invited to join the fun. Thank you for your attention and for sticking around. Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

84 Thank you for your attention! Grazie per l attenzione! Dziękuję za uwagę! Merci de l attention Obrigado pela atençāo Vielen Dank für Ihre Aufmerksamkeit Gracias por su atención Gràcies per la vostra atenció Cám on Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

85 And finally... Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

86 And finally... Many thanks to Richard, Rosario, Hasnaa, Estelle, and all the SADCO people! Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

87 And finally... Many thanks to Richard, Rosario, Hasnaa, Estelle, and all the SADCO people! It s been a great week! Peter R. Wolenski (LSU) Semiconcavity and control SADCO summer / 22

Value Function in deterministic optimal control : sensitivity relations of first and second order

Value Function in deterministic optimal control : sensitivity relations of first and second order Séminaire de Calcul Scientifique du CERMICS Value Function in deterministic optimal control : sensitivity relations of first and second order Hélène Frankowska (IMJ-PRG) 8mars2018 Value Function in deterministic

More information

arxiv: v1 [math.oc] 6 Oct 2011

arxiv: v1 [math.oc] 6 Oct 2011 EXTERIOR SPHERE CONDITION AND TIME OPTIMAL CONTROL FOR DIFFERENTIAL INCLUSIONS PIERMARCO CANNARSA AND KHAI T. NGUYEN arxiv:111.1387v1 [math.oc] 6 Oct 211 Abstract. The minimum time function T( ) of smooth

More information

Second-order sensitivity relations and regularity of the value function for Mayer s problem in optimal control

Second-order sensitivity relations and regularity of the value function for Mayer s problem in optimal control Second-order sensitivity relations and regularity of the value function for Mayer s problem in optimal control Piermarco Cannarsa, Hélène Frankowska, Teresa Scarinci To cite this version: Piermarco Cannarsa,

More information

Value Function and Optimal Trajectories for some State Constrained Control Problems

Value Function and Optimal Trajectories for some State Constrained Control Problems Value Function and Optimal Trajectories for some State Constrained Control Problems Hasnaa Zidani ENSTA ParisTech, Univ. of Paris-saclay "Control of State Constrained Dynamical Systems" Università di Padova,

More information

Neighboring feasible trajectories in infinite dimension

Neighboring feasible trajectories in infinite dimension Neighboring feasible trajectories in infinite dimension Marco Mazzola Université Pierre et Marie Curie (Paris 6) H. Frankowska and E. M. Marchini Control of State Constrained Dynamical Systems Padova,

More information

Hamilton-Jacobi theory for optimal control problems on stratified domains

Hamilton-Jacobi theory for optimal control problems on stratified domains Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 2010 Hamilton-Jacobi theory for optimal control problems on stratified domains Richard Charles Barnard Louisiana

More information

Recent Trends in Differential Inclusions

Recent Trends in Differential Inclusions Recent Trends in Alberto Bressan Department of Mathematics, Penn State University (Aveiro, June 2016) (Aveiro, June 2016) 1 / Two main topics ẋ F (x) differential inclusions with upper semicontinuous,

More information

Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control

Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control Duality and dynamics in Hamilton-Jacobi theory for fully convex problems of control RTyrrell Rockafellar and Peter R Wolenski Abstract This paper describes some recent results in Hamilton- Jacobi theory

More information

The Minimum Speed for a Blocking Problem on the Half Plane

The Minimum Speed for a Blocking Problem on the Half Plane The Minimum Speed for a Blocking Problem on the Half Plane Alberto Bressan and Tao Wang Department of Mathematics, Penn State University University Park, Pa 16802, USA e-mails: bressan@mathpsuedu, wang

More information

Weak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations

Weak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak

More information

An introduction to Mathematical Theory of Control

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

Necessary optimality conditions for optimal control problems with nonsmooth mixed state and control constraints

Necessary optimality conditions for optimal control problems with nonsmooth mixed state and control constraints Necessary optimality conditions for optimal control problems with nonsmooth mixed state and control constraints An Li and Jane J. Ye Abstract. In this paper we study an optimal control problem with nonsmooth

More information

Recent Advances in State Constrained Optimal Control

Recent Advances in State Constrained Optimal Control Recent Advances in State Constrained Optimal Control Richard B. Vinter Imperial College London Control for Energy and Sustainability Seminar Cambridge University, 6th November 2009 Outline of the talk

More information

Existence of viscosity solutions for a nonlocal equation modelling polymer

Existence of viscosity solutions for a nonlocal equation modelling polymer Existence of viscosity solutions for a nonlocal modelling Institut de Recherche Mathématique de Rennes CNRS UMR 6625 INSA de Rennes, France Joint work with P. Cardaliaguet (Univ. Paris-Dauphine) and A.

More information

Minimal time mean field games

Minimal time mean field games based on joint works with Samer Dweik and Filippo Santambrogio PGMO Days 2017 Session Mean Field Games and applications EDF Lab Paris-Saclay November 14th, 2017 LMO, Université Paris-Sud Université Paris-Saclay

More information

Isodiametric problem in Carnot groups

Isodiametric problem in Carnot groups Conference Geometric Measure Theory Université Paris Diderot, 12th-14th September 2012 Isodiametric inequality in R n Isodiametric inequality: where ω n = L n (B(0, 1)). L n (A) 2 n ω n (diam A) n Isodiametric

More information

Maximum Principle and Sensitivity Relations in the Infinite Horizon Optimal Control

Maximum Principle and Sensitivity Relations in the Infinite Horizon Optimal Control Relations in the Infinite Horizon Optimal Control Hélène Frankowska CNRS and UNIVERSITÉ PIERRE et MARIE CURIE In collaboration with Piermarco Cannarsa INdAM Workshop - New Trends in Control Theory and

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

Lecture 3. Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation

Lecture 3. Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Lecture 3. Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Albert Fathi Padova, 14 February 2018 We will present results which are a joint work with Piermarco Cannarsa

More information

Dynamic Blocking Problems for a Model of Fire Propagation

Dynamic Blocking Problems for a Model of Fire Propagation Dynamic Blocking Problems for a Model of Fire Propagation Alberto Bressan Department of Mathematics, Penn State University, University Park, Pa 16802, USA bressan@mathpsuedu February 6, 2012 Abstract This

More information

Introduction. 1 Minimize g(x(0), x(1)) +

Introduction. 1 Minimize g(x(0), x(1)) + ESAIM: Control, Optimisation and Calculus of Variations URL: http://www.emath.fr/cocv/ Will be set by the publisher, UNMAXIMIZED INCLUSION TYPE CONDITIONS FOR NONCONVEX CONSTRAINED OPTIMAL CONTROL PROBLEMS

More information

Local semiconvexity of Kantorovich potentials on non-compact manifolds

Local semiconvexity of Kantorovich potentials on non-compact manifolds Local semiconvexity of Kantorovich potentials on non-compact manifolds Alessio Figalli, Nicola Gigli Abstract We prove that any Kantorovich potential for the cost function c = d / on a Riemannian manifold

More information

LECTURE 15: COMPLETENESS AND CONVEXITY

LECTURE 15: COMPLETENESS AND CONVEXITY LECTURE 15: COMPLETENESS AND CONVEXITY 1. The Hopf-Rinow Theorem Recall that a Riemannian manifold (M, g) is called geodesically complete if the maximal defining interval of any geodesic is R. On the other

More information

Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation

Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Albert Fathi IAS Princeton March 15, 2016 In this lecture, a singularity for a locally Lipschitz real valued function

More information

Mañé s Conjecture from the control viewpoint

Mañé s Conjecture from the control viewpoint Mañé s Conjecture from the control viewpoint Université de Nice - Sophia Antipolis Setting Let M be a smooth compact manifold of dimension n 2 be fixed. Let H : T M R be a Hamiltonian of class C k, with

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 AND UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH LINEAR PROGRAMS EMBEDDED

EXISTENCE AND UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH LINEAR PROGRAMS EMBEDDED EXISTENCE AND UNIQUENESS THEOREMS FOR ORDINARY DIFFERENTIAL EQUATIONS WITH LINEAR PROGRAMS EMBEDDED STUART M. HARWOOD AND PAUL I. BARTON Key words. linear programs, ordinary differential equations, embedded

More information

Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities

Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Pseudo-Poincaré Inequalities and Applications to Sobolev Inequalities Laurent Saloff-Coste Abstract Most smoothing procedures are via averaging. Pseudo-Poincaré inequalities give a basic L p -norm control

More information

Nonpathological Lyapunov functions and discontinuous Carathéodory systems

Nonpathological Lyapunov functions and discontinuous Carathéodory systems Nonpathological Lyapunov functions and discontinuous Carathéodory systems Andrea Bacciotti and Francesca Ceragioli a a Dipartimento di Matematica del Politecnico di Torino, C.so Duca degli Abruzzi, 4-9

More information

The Pontryagin Maximum Principle and a Unified Theory of Dynamic Optimization

The Pontryagin Maximum Principle and a Unified Theory of Dynamic Optimization The Pontryagin Maximum Principle and a Unified Theory of Dynamic Optimization Francis Clarke February 10, 2009 Abstract The Pontryagin maximum principle is the central result of optimal control theory.

More information

BOLZA PROBLEMS WITH GENERAL TIME CONSTRAINTS P. D. LOEWEN AND R. T. ROCKAFELLAR

BOLZA PROBLEMS WITH GENERAL TIME CONSTRAINTS P. D. LOEWEN AND R. T. ROCKAFELLAR PREPRINT, Prepared on 5 February 1996 c 1996 P. D. LOEWEN AND R. T. ROCKAFELLAR BOLZA PROBLEMS WITH GENERAL TIME CONSTRAINTS P. D. LOEWEN AND R. T. ROCKAFELLAR Abstract. This work provides necessary conditions

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

Reformulations for control systems and optimization problems with impulses

Reformulations for control systems and optimization problems with impulses Louisiana State University LSU Digital Commons LSU Doctoral Dissertations Graduate School 214 Reformulations for control systems and optimization problems with impulses Jacob Blanton Louisiana State University

More information

HJ equations. Reachability analysis. Optimal control problems

HJ equations. Reachability analysis. Optimal control problems HJ equations. Reachability analysis. Optimal control problems Hasnaa Zidani 1 1 ENSTA Paris-Tech & INRIA-Saclay Graz, 8-11 September 2014 H. Zidani (ENSTA & Inria) HJ equations. Reachability analysis -

More information

Nonsmooth Analysis in Systems and Control Theory

Nonsmooth Analysis in Systems and Control Theory Nonsmooth Analysis in Systems and Control Theory Francis Clarke Institut universitaire de France et Université de Lyon [January 2008. To appear in the Encyclopedia of Complexity and System Science, Springer.]

More information

On Discontinuous Differential Equations

On Discontinuous Differential Equations On Discontinuous Differential Equations Alberto Bressan and Wen Shen S.I.S.S.A., Via Beirut 4, Trieste 3414 Italy. Department of Informatics, University of Oslo, P.O. Box 18 Blindern, N-316 Oslo, Norway.

More information

Dynamic Blocking Problems for Models of Fire Propagation

Dynamic Blocking Problems for Models of Fire Propagation Dynamic Blocking Problems for Models of Fire Propagation Alberto Bressan Department of Mathematics, Penn State University bressan@math.psu.edu Alberto Bressan (Penn State) Dynamic Blocking Problems 1 /

More information

DERIVED CONES TO REACHABLE SETS OF A CLASS OF SECOND-ORDER DIFFERENTIAL INCLUSIONS

DERIVED CONES TO REACHABLE SETS OF A CLASS OF SECOND-ORDER DIFFERENTIAL INCLUSIONS DERIVED CONES TO REACHABLE SETS OF A CLASS OF SECOND-ORDER DIFFERENTIAL INCLUSIONS AURELIAN CERNEA We consider a second-order differential inclusion and we prove that the reachable set of a certain second-order

More information

Hamilton-Jacobi-Bellman equations for optimal control processes with convex state constraints

Hamilton-Jacobi-Bellman equations for optimal control processes with convex state constraints Hamilton-Jacobi-Bellman equations for optimal control processes with convex state constraints Cristopher Hermosilla a, Richard Vinter b, Hasnaa Zidani c a Departamento de Matemática, Universidad Técnica

More information

Relaxation of Euler-Type Discrete-Time Control System 1

Relaxation of Euler-Type Discrete-Time Control System 1 Relaxation of Euler-Type Discrete-Time Control System 1 Vladimir Veliov 2 Abstract O(h) estimation is obtained for the relaxation of the Euler discretization with step h of some classes of differential

More information

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points Nonlinear Control Lecture # 3 Stability of Equilibrium Points The Invariance Principle Definitions Let x(t) be a solution of ẋ = f(x) A point p is a positive limit point of x(t) if there is a sequence

More information

History (Minimal Dynamics Theorem, Meeks-Perez-Ros 2005, CMC Dynamics Theorem, Meeks-Tinaglia 2008) Giuseppe Tinaglia

History (Minimal Dynamics Theorem, Meeks-Perez-Ros 2005, CMC Dynamics Theorem, Meeks-Tinaglia 2008) Giuseppe Tinaglia History (Minimal Dynamics Theorem, Meeks-Perez-Ros 2005, CMC Dynamics Theorem, Meeks-Tinaglia 2008) Briefly stated, these Dynamics Theorems deal with describing all of the periodic or repeated geometric

More information

Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach

Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Extremal Solutions of Differential Inclusions via Baire Category: a Dual Approach Alberto Bressan Department of Mathematics, Penn State University University Park, Pa 1682, USA e-mail: bressan@mathpsuedu

More information

Stochastic optimal control with rough paths

Stochastic optimal control with rough paths Stochastic optimal control with rough paths Paul Gassiat TU Berlin Stochastic processes and their statistics in Finance, Okinawa, October 28, 2013 Joint work with Joscha Diehl and Peter Friz Introduction

More information

Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary

Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Deforming conformal metrics with negative Bakry-Émery Ricci Tensor on manifolds with boundary Weimin Sheng (Joint with Li-Xia Yuan) Zhejiang University IMS, NUS, 8-12 Dec 2014 1 / 50 Outline 1 Prescribing

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

Necessary conditions in optimal control and in the calculus of variations

Necessary conditions in optimal control and in the calculus of variations Necessary conditions in optimal control and in the calculus of variations Francis Clarke Institut universitaire de France et Institut Camille Jordan Université Claude Bernard Lyon 1 69622 Villeurbanne,

More information

Theory and Applications of Constrained Optimal Control Proble

Theory and Applications of Constrained Optimal Control Proble Theory and Applications of Constrained Optimal Control Problems with Delays PART 1 : Mixed Control State Constraints Helmut Maurer 1, Laurenz Göllmann 2 1 Institut für Numerische und Angewandte Mathematik,

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

SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS

SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS SUPERCONVERGENCE PROPERTIES FOR OPTIMAL CONTROL PROBLEMS DISCRETIZED BY PIECEWISE LINEAR AND DISCONTINUOUS FUNCTIONS A. RÖSCH AND R. SIMON Abstract. An optimal control problem for an elliptic equation

More information

FEEDBACK DIFFERENTIAL SYSTEMS: APPROXIMATE AND LIMITING TRAJECTORIES

FEEDBACK DIFFERENTIAL SYSTEMS: APPROXIMATE AND LIMITING TRAJECTORIES STUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLIX, Number 3, September 2004 FEEDBACK DIFFERENTIAL SYSTEMS: APPROXIMATE AND LIMITING TRAJECTORIES Abstract. A feedback differential system is defined as

More information

Introduction to Real Analysis Alternative Chapter 1

Introduction to Real Analysis Alternative Chapter 1 Christopher Heil Introduction to Real Analysis Alternative Chapter 1 A Primer on Norms and Banach Spaces Last Updated: March 10, 2018 c 2018 by Christopher Heil Chapter 1 A Primer on Norms and Banach Spaces

More information

Functional Analysis I

Functional Analysis I Functional Analysis I Course Notes by Stefan Richter Transcribed and Annotated by Gregory Zitelli Polar Decomposition Definition. An operator W B(H) is called a partial isometry if W x = X for all x (ker

More information

Differential Games II. Marc Quincampoix Université de Bretagne Occidentale ( Brest-France) SADCO, London, September 2011

Differential Games II. Marc Quincampoix Université de Bretagne Occidentale ( Brest-France) SADCO, London, September 2011 Differential Games II Marc Quincampoix Université de Bretagne Occidentale ( Brest-France) SADCO, London, September 2011 Contents 1. I Introduction: A Pursuit Game and Isaacs Theory 2. II Strategies 3.

More information

Introduction to Empirical Processes and Semiparametric Inference Lecture 13: Entropy Calculations

Introduction to Empirical Processes and Semiparametric Inference Lecture 13: Entropy Calculations Introduction to Empirical Processes and Semiparametric Inference Lecture 13: Entropy Calculations Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations Research

More information

asymptotic behaviour of singularly perturbed control systems in the non-periodic setting

asymptotic behaviour of singularly perturbed control systems in the non-periodic setting UNIVERSITÀ DI ROMA LA SAPIENZA Doctoral Thesis asymptotic behaviour of singularly perturbed control systems in the non-periodic setting Author Nguyen Ngoc Quoc Thuong Supervisor Prof. Antonio Siconolfi

More information

Optimization and Optimal Control in Banach Spaces

Optimization and Optimal Control in Banach Spaces Optimization and Optimal Control in Banach Spaces Bernhard Schmitzer October 19, 2017 1 Convex non-smooth optimization with proximal operators Remark 1.1 (Motivation). Convex optimization: easier to solve,

More information

Lecture 5: The Bellman Equation

Lecture 5: The Bellman Equation Lecture 5: The Bellman Equation Florian Scheuer 1 Plan Prove properties of the Bellman equation (In particular, existence and uniqueness of solution) Use this to prove properties of the solution Think

More information

PROPERTIES OF C 1 -SMOOTH FUNCTIONS WITH CONSTRAINTS ON THE GRADIENT RANGE Mikhail V. Korobkov Russia, Novosibirsk, Sobolev Institute of Mathematics

PROPERTIES OF C 1 -SMOOTH FUNCTIONS WITH CONSTRAINTS ON THE GRADIENT RANGE Mikhail V. Korobkov Russia, Novosibirsk, Sobolev Institute of Mathematics PROPERTIES OF C 1 -SMOOTH FUNCTIONS WITH CONSTRAINTS ON THE GRADIENT RANGE Mikhail V. Korobkov Russia, Novosibirsk, Sobolev Institute of Mathematics e-mail: korob@math.nsc.ru Using Gromov method of convex

More information

LINEAR-CONVEX CONTROL AND DUALITY

LINEAR-CONVEX CONTROL AND DUALITY 1 LINEAR-CONVEX CONTROL AND DUALITY R.T. Rockafellar Department of Mathematics, University of Washington Seattle, WA 98195-4350, USA Email: rtr@math.washington.edu R. Goebel 3518 NE 42 St., Seattle, WA

More information

The optimal partial transport problem

The optimal partial transport problem The optimal partial transport problem Alessio Figalli Abstract Given two densities f and g, we consider the problem of transporting a fraction m [0, min{ f L 1, g L 1}] of the mass of f onto g minimizing

More information

1 Cheeger differentiation

1 Cheeger differentiation 1 Cheeger differentiation after J. Cheeger [1] and S. Keith [3] A summary written by John Mackay Abstract We construct a measurable differentiable structure on any metric measure space that is doubling

More information

On reduction of differential inclusions and Lyapunov stability

On reduction of differential inclusions and Lyapunov stability 1 On reduction of differential inclusions and Lyapunov stability Rushikesh Kamalapurkar, Warren E. Dixon, and Andrew R. Teel arxiv:1703.07071v5 [cs.sy] 25 Oct 2018 Abstract In this paper, locally Lipschitz

More information

1 Directional Derivatives and Differentiability

1 Directional Derivatives and Differentiability Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=

More information

Journal of Differential Equations

Journal of Differential Equations J. Differential Equations 252 2012 1912 1933 Contents lists available at SciVerse ScienceDirect Journal of Differential Equations www.elsevier.com/locate/jde L estimates on trajectories confined to a closed

More information

DUALIZATION OF SUBGRADIENT CONDITIONS FOR OPTIMALITY

DUALIZATION OF SUBGRADIENT CONDITIONS FOR OPTIMALITY DUALIZATION OF SUBGRADIENT CONDITIONS FOR OPTIMALITY R. T. Rockafellar* Abstract. A basic relationship is derived between generalized subgradients of a given function, possibly nonsmooth and nonconvex,

More information

ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS

ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS MATHEMATICS OF OPERATIONS RESEARCH Vol. 28, No. 4, November 2003, pp. 677 692 Printed in U.S.A. ON A CLASS OF NONSMOOTH COMPOSITE FUNCTIONS ALEXANDER SHAPIRO We discuss in this paper a class of nonsmooth

More information

Locally Lipschitzian Guiding Function Method for ODEs.

Locally Lipschitzian Guiding Function Method for ODEs. Locally Lipschitzian Guiding Function Method for ODEs. Marta Lewicka International School for Advanced Studies, SISSA, via Beirut 2-4, 3414 Trieste, Italy. E-mail: lewicka@sissa.it 1 Introduction Let f

More information

Math 4317 : Real Analysis I Mid-Term Exam 1 25 September 2012

Math 4317 : Real Analysis I Mid-Term Exam 1 25 September 2012 Instructions: Answer all of the problems. Math 4317 : Real Analysis I Mid-Term Exam 1 25 September 2012 Definitions (2 points each) 1. State the definition of a metric space. A metric space (X, d) is set

More information

MATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES

MATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES MATH 4200 HW: PROBLEM SET FOUR: METRIC SPACES PETE L. CLARK 4. Metric Spaces (no more lulz) Directions: This week, please solve any seven problems. Next week, please solve seven more. Starred parts of

More information

Generic Aubry sets on surfaces

Generic Aubry sets on surfaces Université de Nice - Sophia Antipolis & Institut Universitaire de France Nanjing University (June 10th, 2013) Setting Let M be a smooth compact manifold of dimension n 2 be fixed. Let H : T M R be a Hamiltonian

More information

Subdifferential representation of convex functions: refinements and applications

Subdifferential representation of convex functions: refinements and applications Subdifferential representation of convex functions: refinements and applications Joël Benoist & Aris Daniilidis Abstract Every lower semicontinuous convex function can be represented through its subdifferential

More information

McGill University Math 354: Honors Analysis 3

McGill University Math 354: Honors Analysis 3 Practice problems McGill University Math 354: Honors Analysis 3 not for credit Problem 1. Determine whether the family of F = {f n } functions f n (x) = x n is uniformly equicontinuous. 1st Solution: The

More information

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION

PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION PARTIAL REGULARITY OF BRENIER SOLUTIONS OF THE MONGE-AMPÈRE EQUATION ALESSIO FIGALLI AND YOUNG-HEON KIM Abstract. Given Ω, Λ R n two bounded open sets, and f and g two probability densities concentrated

More information

Geometry and topology of continuous best and near best approximations

Geometry and topology of continuous best and near best approximations Journal of Approximation Theory 105: 252 262, Geometry and topology of continuous best and near best approximations Paul C. Kainen Dept. of Mathematics Georgetown University Washington, D.C. 20057 Věra

More information

Tutorial on Control and State Constrained Optimal Control Pro. Control Problems and Applications Part 3 : Pure State Constraints

Tutorial on Control and State Constrained Optimal Control Pro. Control Problems and Applications Part 3 : Pure State Constraints Tutorial on Control and State Constrained Optimal Control Problems and Applications Part 3 : Pure State Constraints University of Münster Institute of Computational and Applied Mathematics SADCO Summer

More information

The Rocket Car. UBC Math 403 Lecture Notes by Philip D. Loewen

The Rocket Car. UBC Math 403 Lecture Notes by Philip D. Loewen The Rocket Car UBC Math 403 Lecture Notes by Philip D. Loewen We consider this simple system with control constraints: ẋ 1 = x, ẋ = u, u [ 1, 1], Dynamics. Consider the system evolution on a time interval

More information

On fine differentiability properties of horizons and applications to Riemannian geometry

On fine differentiability properties of horizons and applications to Riemannian geometry On fine differentiability properties of horizons and applications to Riemannian geometry Piotr T. Chruściel Département de Mathématiques Faculté des Sciences Parc de Grandmont F37200 Tours, France Joseph

More information

Regularity of a kind of marginal functions in Hilbert spaces

Regularity of a kind of marginal functions in Hilbert spaces Regularity of a kind of marginal functions in Hilbert spaces Fátima F. Pereira and Vladimir V. Goncharov IMA-UE, Rua Romão Ramalho 59, 7000-671, Évora, Portugal e-mail: fmfp@uevora.pt and goncha@uevora.pt

More information

Iowa State University. Instructor: Alex Roitershtein Summer Homework #5. Solutions

Iowa State University. Instructor: Alex Roitershtein Summer Homework #5. Solutions Math 50 Iowa State University Introduction to Real Analysis Department of Mathematics Instructor: Alex Roitershtein Summer 205 Homework #5 Solutions. Let α and c be real numbers, c > 0, and f is defined

More information

LECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE

LECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE LECTURE 25: REVIEW/EPILOGUE LECTURE OUTLINE CONVEX ANALYSIS AND DUALITY Basic concepts of convex analysis Basic concepts of convex optimization Geometric duality framework - MC/MC Constrained optimization

More information

A convergence result for an Outer Approximation Scheme

A convergence result for an Outer Approximation Scheme A convergence result for an Outer Approximation Scheme R. S. Burachik Engenharia de Sistemas e Computação, COPPE-UFRJ, CP 68511, Rio de Janeiro, RJ, CEP 21941-972, Brazil regi@cos.ufrj.br J. O. Lopes Departamento

More information

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation Instructor: Moritz Hardt Email: hardt+ee227c@berkeley.edu Graduate Instructor: Max Simchowitz Email: msimchow+ee227c@berkeley.edu

More information

The optimization problem of interest in this paper, which will be referred to as the

The optimization problem of interest in this paper, which will be referred to as the SIAM J. CONTROL OPTIM. Vol. 44, No. 3, pp. 939 968 c 2005 Society for Industrial Applied Mathematics MINIMAX OPTIMAL CONTROL R. B. VINTER Abstract. This paper provides a framework for deriving necessary

More information

Derivatives. Differentiability problems in Banach spaces. Existence of derivatives. Sharpness of Lebesgue s result

Derivatives. Differentiability problems in Banach spaces. Existence of derivatives. Sharpness of Lebesgue s result Differentiability problems in Banach spaces David Preiss 1 Expanded notes of a talk based on a nearly finished research monograph Fréchet differentiability of Lipschitz functions and porous sets in Banach

More information

Convex Geometry. Carsten Schütt

Convex Geometry. Carsten Schütt Convex Geometry Carsten Schütt November 25, 2006 2 Contents 0.1 Convex sets... 4 0.2 Separation.... 9 0.3 Extreme points..... 15 0.4 Blaschke selection principle... 18 0.5 Polytopes and polyhedra.... 23

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

CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.

CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008. 1 ECONOMICS 594: LECTURE NOTES CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS W. Erwin Diewert January 31, 2008. 1. Introduction Many economic problems have the following structure: (i) a linear function

More information

Normality and Nondegeneracy for Optimal Control Problems with State Constraints

Normality and Nondegeneracy for Optimal Control Problems with State Constraints Normality and Nondegeneracy for Optimal Control Problems with State Constraints Fernando Fontes, Hélène Frankowska To cite this version: Fernando Fontes, Hélène Frankowska. Normality and Nondegeneracy

More information

Calculating the domain of attraction: Zubov s method and extensions

Calculating the domain of attraction: Zubov s method and extensions Calculating the domain of attraction: Zubov s method and extensions Fabio Camilli 1 Lars Grüne 2 Fabian Wirth 3 1 University of L Aquila, Italy 2 University of Bayreuth, Germany 3 Hamilton Institute, NUI

More information

Second order forward-backward dynamical systems for monotone inclusion problems

Second order forward-backward dynamical systems for monotone inclusion problems Second order forward-backward dynamical systems for monotone inclusion problems Radu Ioan Boţ Ernö Robert Csetnek March 6, 25 Abstract. We begin by considering second order dynamical systems of the from

More information

Sequential Unconstrained Minimization: A Survey

Sequential Unconstrained Minimization: A Survey Sequential Unconstrained Minimization: A Survey Charles L. Byrne February 21, 2013 Abstract The problem is to minimize a function f : X (, ], over a non-empty subset C of X, where X is an arbitrary set.

More information

MATH 819 FALL We considered solutions of this equation on the domain Ū, where

MATH 819 FALL We considered solutions of this equation on the domain Ū, where MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,

More information

Peak Point Theorems for Uniform Algebras on Smooth Manifolds

Peak Point Theorems for Uniform Algebras on Smooth Manifolds Peak Point Theorems for Uniform Algebras on Smooth Manifolds John T. Anderson and Alexander J. Izzo Abstract: It was once conjectured that if A is a uniform algebra on its maximal ideal space X, and if

More information

Graph Completions for Impulsive Feedback Controls

Graph Completions for Impulsive Feedback Controls Graph Completions for Impulsive Feedback Controls Alberto Bressan and Marco Mazzola * Department of Mathematics, Penn State University. ** Université Pierre et Marie Curie, Paris VI. e-mails: bressan@math.psu.edu,

More information

Computation of Local ISS Lyapunov Function Via Linear Programming

Computation of Local ISS Lyapunov Function Via Linear Programming Computation of Local ISS Lyapunov Function Via Linear Programming Huijuan Li Joint work with Robert Baier, Lars Grüne, Sigurdur F. Hafstein and Fabian Wirth Institut für Mathematik, Universität Würzburg

More information

MAC College Algebra

MAC College Algebra MAC 05 - College Algebra Name Review for Test 2 - Chapter 2 Date MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Find the exact distance between the

More information

Implications of the Constant Rank Constraint Qualification

Implications of the Constant Rank Constraint Qualification Mathematical Programming manuscript No. (will be inserted by the editor) Implications of the Constant Rank Constraint Qualification Shu Lu Received: date / Accepted: date Abstract This paper investigates

More information

Weak convergence and large deviation theory

Weak convergence and large deviation theory First Prev Next Go To Go Back Full Screen Close Quit 1 Weak convergence and large deviation theory Large deviation principle Convergence in distribution The Bryc-Varadhan theorem Tightness and Prohorov

More information

The Caratheodory Construction of Measures

The Caratheodory Construction of Measures Chapter 5 The Caratheodory Construction of Measures Recall how our construction of Lebesgue measure in Chapter 2 proceeded from an initial notion of the size of a very restricted class of subsets of R,

More information