Many-on-One Pursuit. M. Pachter. AFIT Wright-Patterson A.F.B., OH 45433

Size: px
Start display at page:

Download "Many-on-One Pursuit. M. Pachter. AFIT Wright-Patterson A.F.B., OH 45433"

Transcription

1 Many-on-One Pursuit M. Pachter AFIT Wright-Patterson A.F.B., OH The views expressed on these slides are those of the author and do not reflect the official policy or position of the United States Air Force, Department of Defense, or the US Gov. October 13, 2018 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

2 AI and Control The views expressed on this slide are those of the author and do not reflect the official policy or position of the US Air Force, DoD, or the US Gov. Current state of AI Enumeration, ML/overfitting, Heuristics/hacking, ES/no AI -It s like alchemy So what is real AI? AI is decision making powered by mathematics and machine computing An attempt at a definition in the Control context AI: Optimal control in the presence of uncertainty and adversarial action. Example: Solution of a DG The adversary is an algorithm and the human player is a loser. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

3 Pursuit and Evasion M. Pachter (2018) Many-on-One Pursuit October 13, / 64

4 N-on-One Pursuit and Evasion Figure: Multi Pursuers M. Pachter (2018) Many-on-One Pursuit October 13, / 64

5 Three Agents M. Pachter (2018) Many-on-One Pursuit October 13, / 64

6 Apollonius Circle Figure: Apollonius Circle M. Pachter (2018) Many-on-One Pursuit October 13, / 64

7 Pincer Movement Pursuit Figure: Solution of Two Cutters and Fugitive Ship Game M. Pachter (2018) Many-on-One Pursuit October 13, / 64

8 No Way Out M. Pachter (2018) Many-on-One Pursuit October 13, / 64

9 Breakout M. Pachter (2018) Many-on-One Pursuit October 13, / 64

10 P 2 Redundant Figure: One Cutter Action M. Pachter (2018) Many-on-One Pursuit October 13, / 64

11 Reduced State Space Figure: Rotating Reference Frame M. Pachter (2018) Many-on-One Pursuit October 13, / 64

12 Reduced State Space The state space of the Two Cutters and Fugitive Ship differential game is reduced to the first quadrant of R 3 : R 3 1 {(x P, x E, y E ) x P 0, y E 0} Symmetry allows us to confine our attention to the case where x E 0, that is, the state will evolve in the positive orthant of R 3 : R 3 + = {(x P, x E, y E ) x P 0, x E 0, y E 0}, The three-state nonlinear dynamics of the differential game are ẋ P = 1 2 (cos χ cos ψ), x P(0) = x P0 (1) ẋ E = µ cos φ 1 2 (cos χ + cos ψ) + 1 y E 2 x P (sin χ sin ψ), x E (0) = x E0 (2) ẏ E = µ sin φ 1 2 (sin χ + sin ψ) 1 x E (sin χ sin ψ), y E (0) = y E0 (3) 2 x P M. Pachter (2018) Many-on-One Pursuit October 13, / 64

13 R 1 /R 1,2 Boundary Figure: Critical Configuration M. Pachter (2018) Many-on-One Pursuit October 13, / 64

14 Soln. of Game of Kind Theorem 1 During optimal play 1. The Evader is singlehandedly captured in PP by P 1 if the state is in the set R 1 ; the set R 1 is the wedge formed by the quarter planes {(x P, x E, y E ) x P = 0, x E 0, y E 0} and {(x P, x E, y E ) x E = µx P, x P 0, y E 0}. 2. The Evader is singlehandedly captured in PP by P 2 if the state is in the set R 2 ; the set R 2 is the mirror image of R 1 about the plane x E = The Evader is cooperatively and isochronously captured by P 1 and P 2 if the state is in the wedge shaped set R 1,2 = {(x P, x E, y E ) µx P x E µx P, x P 0, y E 0} M. Pachter (2018) Many-on-One Pursuit October 13, / 64

15 Soln. of Game of Kind M. Pachter (2018) Many-on-One Pursuit October 13, / 64

16 Soln. of Game of Kind, l > 0 Figure: Positive Orthant, l > 0 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

17 Soln. of Game of Degree in R 1,2 Theorem 2 E is isochronously captured by the P 1 &P 2 team in R 1,2 = {(x P, x E, y E ) µx P x E µx P, x P 0, y E 0} The optimal state feedback strategies are y I sin ψ =, cos ψ = xp 2 + y I 2 y I x P x 2 P + y 2 I sin χ =, cos χ = (5) xp 2 + y I 2 xp 2 + y I 2 sin φ y I y E =, cos φ x = E (6) (y I y E ) 2 + xe 2 (y I y E ) 2 + xe 2 Time-to-capture/Value function: V (x P, x E, y E ) = xp 2 + y I 2 ; the function y I (x P, x E, y E ) = 1 1 µ 2 [y E + sign(y E ) µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2 )] (7) M. Pachter (2018) Many-on-One Pursuit October 13, / 64 x P (4)

18 Proof by Isaacs Method We begin from the end. The terminal surface in R 1,2, T, is a circular arc in the plane x E = 0. T = {(x P, x E, y E ) x P = l cos ξ, x E = 0, y E = l sin ξ, 0 ξ π 2 } (8) Interested in the inward pointing normals n to the terminal surface T : They set the terminal conditions of the costate vector But although the problem formulation is physically sound (l > 0), the terminal manifold T is rank deficient The normals to the terminal surface at a point of the surface are not unique. From eq. (8) we extract the terminal costates cos ξ λ(tf ) = a b sin ξ M. Pachter (2018) Many-on-One Pursuit October 13, / 64

19 Proof by Isaacs Method The scalars a > 0, b < 0, and 0 ξ π 2 ; in the half of R 1,2 where x E < 0, b > 0 and in the plane x E = 0, b = 0. The size l of the pursuers capture set plays no role in setting the terminal costate exclusively interested in point capture: l 0. The costate λ is related to the partials of the Value function V (x P, x E, y E ): λ xp = V xp, λ xe = V xe and λ ye = V ye. The Hamiltonian H = [(λ y E + y E λ xe x E λ ye x P ) sin ψ + (λ xe + λ xp ) cos ψ + (λ ye y E λ xe x E λ ye x P ) sin χ + (λ xe λ xp ) cos χ] + µ(λ ye sin φ + λ xe cos φ) (9) Maximizing the Hamiltonian in χ and ψ and minimizing the Hamiltonian in φ helps to characterize the optimal controls. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

20 Proof by Isaacs Method Euler-Lagrange Equations ẋ P = 1 2 [ λ xp λ xe (λ λ xe xp )2 + (λ y E λ xe x E λ ye ye x P ) 2 + (λ xe + λ xp ) ], x P (0) = x (λ xe + λ xp ) 2 + (λ ye + y P0 E λ xe x E λ ye x P ) 2 ẋ E = λ xe µ 1 λ 2 + xe λ2 2 [ λ xp λ xe ye (λ λ xe xp )2 + (λ y E λ xe x E λ ye ye x P ) 2 (λ xe + λ xp ) ] (λ xe + λ xp ) 2 + (λ ye + y E λ xe x E λ ye x P ) y E [ 2 x P y E λ xe x E λ ye x P λ ye (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye x P ) 2 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

21 Proof by Isaacs Method λ ye ẏ E = µ 1 λ 2 + xe λ2 2 [ ye + x E λ ye y E λ xe x P y E λ xe x E λ ye x P λ ye ] (λ xe + λ xp ) 2 + (λ ye + y E λ xe x E λ ye x P ) 2 1 x E [ 2 x P λ ye (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye x P ) 2 y E λ xe x E λ ye x P λ ye (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye x P ) 2 x E λ ye y E λ xe x P λ ye ], y E (0) = y (λ xe + λ xp ) 2 + (λ ye + y E0 E λ xe x E λ ye x P ) 2 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

22 Proof by Isaacs Method λ xp = 1 2 y E λ xe x E λ ye x 2 P [ x E λ ye y E λ xe x P y E λ xe x E λ ye x P λ ye (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye x P ) 2 λ ye ], λ (λ xe + λ xp ) 2 + (λ ye + y xp (t f ) = a cos ξ E λ xe x E λ ye x P ) 2 λ xe = 1 λ ye [ 2 x P y E λ xe x E λ ye x P λ ye (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye x P ) 2 x E λ ye y E λ xe x P λ ye, λ (λ xe + λ xp ) 2 + (λ ye + y xe (t f ) = ab E λ xe x E λ ye x P ) 2 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

23 Proof by Isaacs Method λ ye = 1 2 λ xe x P [ x E λ ye y E λ xe x P λ ye (λ xe + λ xp ) 2 + (λ ye + y E λ xe x E λ ye x P ) 2 y E λ xe x E λ ye x P λ ye, λ (λ xe λ xp ) 2 + (λ ye y xe (t f ) = a sin ξ E λ xe x E λ ye x P ) 2 Inserting the expressions which characterize the optimal controls into the Hamiltonian (9) yields the optimal smooth Hamiltonian H = [ (λ xe λ xp ) 2 + (λ ye y E λ xe x E λ ye ) x 2 P + (λ xe + λ xp ) 2 + (λ ye + y E λ xe x E λ ye x P ) 2 ] µ λ 2 x E + λ 2 y E M. Pachter (2018) Many-on-One Pursuit October 13, / 64

24 Proof by Isaacs Method The Hamiltonian vanishes and evaluating the optimal Hamiltonian at t = t f yields a relationship which encompasses the parameters a, b and ξ: 1 = cos ξ [ (b cos ξ) 2 cos 2 ξ + (b cos ξ) 2 sin 2 ξ + (b + cos ξ) 2 cos 2 ξ + (b + cos ξ) 2 sin 2 ξ]a µ b + sin ξ a The three parameters, 0 ξ π 2, a > 0 and b feature in the solution of the Euler-Lagrange equations. BUT: The vanishing optimal Hamiltonian yields cos ξ a = 1 2 ( b cos ξ + b + cos ξ ) µ b + sin ξ cos ξ So The Euler-Lagrange equations yield a family of optimal trajectories parameterized by two independent parameters b and ξ. Can fill in our three dimensional state space region R 1,2. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

25 Proof by Isaacs Method The Euler-Lagrange equations are integrated in retrograde fashion starting out from the initial condition x P = 0, x E = 0, y E = 0, λ xp = a cos ξ, λ xe = ab, λ ye = a sin ξ The solution of the Euler-Lagrange equations which yield a trajectory which terminates at the initial state (x P0, x E0, y E0 ) is x P (t) = (1 t t f )x P0 x E (t) = (1 t t f )x E0 (10) where t f = 1 1 µ 2 y E (t) = (1 t t f )y E0, 0 t t f (1 µ 2 )(x 2 P 0 x 2 E 0 ) + (1 + µ 2 )y 2 E 0 + 2y E µ 2 y 2 E 0 + (1 µ 2 )(µ 2 x 2 P 0 x 2 E 0 ) M. Pachter (2018) Many-on-One Pursuit October 13, / 64

26 Proof by Isaacs Method To complete the proof it behooves on us to show that (x P0, x E0, y E0 ) R 1,2, b, 0 < ξ < π 2 and t f > 0 s.t. the following holds. x P0 = cos ξ t f (11) x E0 = b µ b 2 + sin 2 ξ t f (12) 1 y E0 = (1 µ b 2 + sin 2 ξ ) sin ξ t f (13) We must be able to solve the three equations (11)-(13) in the three unknowns b, ξ and t f. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

27 Proof by Isaacs Method It all boils down to proving the existence of a solution 0 ξ π 2 to y E0 cos ξ = x P0 sin ξ µ 2 xp 2 0 xe 2 0 cos 2 ξ (14) It is where ξ = Arccos( x) x = (1 µ 2 )(x 2 xp 2 P 0 xe 2 0 ) + (1 + µ 2 )ye y E0 µ 2 (xp ye 2 0 ) (1 µ 2 )xe xp xe ye xE 2 0 ye xP 2 0 ye 2 0 2xP 2 0 xe 2 0 and consequently (1 x)x b = µ 2 ( x P 0 x E0 ) 2 x sign(x E ), t f = x P0 x M. Pachter (2018) Many-on-One Pursuit October 13, / 64

28 The Geometric Method is Correct Theorem 2 In the Two Cutters and Fugitive Ship differential game, in the state space region R 1,2 where both Pursuers P 1 and P 2 actively engage the Evader, the players optimal state feedback strategies are sin ψ = y E + µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2) (1 µ 2 )(xp 2 x E 2) + (1 + µ2 )ye 2 + 2y E µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2) cos ψ = (1 µ 2 )x P (1 µ 2 )(xp 2 x E 2) + (1 + µ2 )ye 2 + 2y E µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2) χ = π ψ M. Pachter (2018) Many-on-One Pursuit October 13, / 64

29 The Geometric Method is Correct sin φ = 1 µ µ 2 y E + µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2 ) (1 µ 2 )(xp 2 x E 2) + (1 + µ2 )ye 2 + 2y E µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2) cos φ = 1 µ (1 µ 2 )x E (1 µ 2 )(xp 2 x E 2) + (1 + µ2 )ye 2 + 2y E µ 2 ye 2 + (1 µ2 )(µ 2 xp 2 x E 2) M. Pachter (2018) Many-on-One Pursuit October 13, / 64

30 The Geometric Method is Correct The Value function V (x P, x E, y E ) = 1 1 µ 2 (1 µ 2 )(x 2 P x 2 E ) + (1 + µ2 )y 2 E + 2y E µ 2 y 2 E + (1 µ2 )(µ 2 x 2 P x 2 E ) The field of primary optimal trajectories/regular characteristics covers the entire state space R1 3 No singular surfaces here. The geometric method yields the correct solution of the differential game. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

31 Discussion Addressed the Two Cutters and Fugitive Ship differential game. 1. Provided the solution of the Game of Kind: Partitioned the state space into the R 1, R 2 and R 1,2 regions. 2. Proved the validity of the geometric solution using Isaacs method. -Optimal Hamiltonian is smooth. -Value function is C 1. -Optimal flow field covers state space. -Regular characteristics only no singular surfaces here. Proving the validity of the geometric solution was, surprisingly, difficult: Used the known optimal solution provided by the geometric method. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

32 N-on-One Pursuit and Evasion Figure: Multi Pursuers M. Pachter (2018) Many-on-One Pursuit October 13, / 64

33 N > 2 Pursuers Safe Region (SR) Must solve 1 2N(N 1) quadratic equations. Have M N active Apollonius circles. Must construct the intersection of M Apollonius disks. SR = N 1 i M C i SR: Convex, has M corners, the sides of the M- gon are circular arcs. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

34 Safe Region N = M = 3, µ = 2 Figure: Reuleux Triangle SR M. Pachter (2018) Many-on-One Pursuit October 13, / 64

35 Optimality Principle Solve the OR constrained Obnoxious Facility Location problem: Optimality Principle The N cities are {P 1,..., P N } The constraint set is the SR Aim point is the Obnoxious Facility (P, i 1,..., i l ) = arg max P SR min 1 i 3 dist(p i, P) I = P Active Pursuers: P I 1,...,P I l M. Pachter (2018) Many-on-One Pursuit October 13, / 64

36 Optimal Group Pursuit N pursuers + one evader; simple motion; capture range l 0. Speed ratios: µ i > 1, i = 1,..., N. Sate Feedback Capture Strategy Algorithm Given the pursuers and the evader s instantaneous positions 1. Calculate the centers and radii of the N Apollonius circles/disks. 2. Find BSR vertex candidates by solving 1 2N(N 1) sets of bivariate quadratic equations in the variables x and y requires the solution of 1 2N(N 1) quadratic equations. This designates the set of M N active pursuers via the positivity of the quadratic equations discriminants. 3. Calculate the points antipodal to E on the circumference of each of the M active Apollonius circles. This requires the solution of a quadratic equation in total, M quadratic equations are solved. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

37 Optimal Group Pursuit 4. Check whether the antipodal point is in all the remaining M-1 Apollonius disks requires checking M(M 1) quadratic inequalities. If such an antipodal point exists and say, it resides on the circumference of the ith Apollonius circle, designate it the aim point I of Pursuer i and of the Evader the pursuit degenerates into a tail chase/pp where only Pursuer i is active. 5. If no such antipodal point exists and because a vertex V = (x, y) of the BSR/M-gon must be included in the Apollonius disks of all the M active pursuers, if M > 2, to determine the M vertices of the BSR/M-gon, the M(M 1)(M 2) quadratic inequalities (x i x Oj ) 2 + (y i y Oj ) 2 ρ 2 j, j {1,..., M,, j i}, 1 i M(M 1) must be checked only M inequalities will be satisfied, rendering the BSR s/m-gon s vertices Have constructed the SR. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

38 Optimal Group Pursuit 6. Potential aim point : Vertex V k of BSR farthest from E Formed by the intersection of the Apollonius circles C i(k ) and C j(k ) This identifies two critical pursuers, i(k ) and j(k ), which cooperatively capture the evader by heading to the aim point I = V k BSR Also E heads to the aim point I = V k BSR. Have pincer movement pursuit with P i(k ) and P j(k ) on a CC with E. Pursuer closest to E, but not P i(k ) or P j(k ), employs PP. Remaining N 3 pursuers are redundant. The Value of the game V (P 1,..., P N, E) = dist(e, V k ) (15) M. Pachter (2018) Many-on-One Pursuit October 13, / 64

39 Optimal Group Pursuit 7. Additional candidate aim point P provided by the solution of Constrained Obnoxious Facility Location The role of the cities is assumed by the pursuers, the constraint set is SR. (P, i 1,...i l ) = arg max P SR min 1 i M dist(p i, P) (16) dist(p, P i 1 ) = dist(p, P i 2 ) =...dist(p, P i l ) R, l 1 Optimal location P of the obnoxious facility: Potential aim point I Consider the condition dist(v k, P i(k )) R (17) Let R ɛ d R + ɛ. The optimal aim point is specified as follows P if dist(v k, P i(k )) < R ɛ I = V k if dist(v k, P i(k )) > R + ɛ E if dist(v k, P i(k )) = d M. Pachter (2018) Many-on-One Pursuit October 13, / 64

40 Optimal Group Pursuit If condition (17) does not hold: I = P Int(SR). The l pursuers who are closest to P head toward I = P. The remaining N l pursuers are redundant. E runs toward the aim point I = P in the interior of the SR Upon arriving to P Int(SR), E stops. Awaits isochronous capture by designated pursuers P i 1,...P i l. The Value V (P 1, P 2, P 3, E) = dist(p, P i 1 ) µ If condition (17) holds: I = V k Pursuers P i(k ) and P j(k ) head toward the vertex V k. If Condition (17) holds with equality the pursuers head toward E. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

41 Three Pursuers N = M = 3 (a) If pursuers P 1 P 2 P 3 is obtuse, E heads to the vertex k of the BSR. The vertex V k BSR being generated by C i(k ) C j(k ) Capture of E effected by pursuers P i(k ) and P j(k ) (b) If pursuers P 1 P 2 P 3 is acute, is O SR? i) If O SR eq. (17) E heads toward the vertex V k BSR and as in case (a) capture of E is effected by P i(k ) and P j(k ). ii) If O SR and dist(v k, P i(k )) < R, eq. (17) E heads toward O in the interior of the SR, and so do P 1, P 2, P 3. Since E Interior(SR) E arrives at O before the pursuers. Upon arriving to O, E stays put, awaiting his isochronous capture by the three pursuers. The Value V = R µ. If however dist(v k, P i(k )) > R, eq. (17) E V k and so do P i(k ) and P j(k ), irrespective of whether O is in the SR. The Value V = 1 µ dist(v k, P i(k )) M. Pachter (2018) Many-on-One Pursuit October 13, / 64

42 Three Pursuers Symmetric State O SR dist(v k, P i(k )) = d State symmetric E c O V k s.t. dist(e, E c ) δ, 0 < δ << 1 The pursuers team will incur a slight loss in optimality: iii) If state is symmetric the pursuers head toward E. The Evader has his pick: Head toward O or head toward V k The Value of the game is V = 1 µ R M. Pachter (2018) Many-on-One Pursuit October 13, / 64

43 Three Pursuers N = M = 3, µ = 2 P 1 C 3 P 2 I 1, 3 I 1, 2 C 1 C 2 E O P 3 I 2, 3 I * Figure: P 2, P 3 and E head to I 2,3, Redundant Pursuer 1 Employs PP M. Pachter (2018) Many-on-One Pursuit October 13, / 64

44 Figure: Hexagon in Game with Three Pursuers M. Pachter (2018) Many-on-One Pursuit October 13, / 64 Three Pursuers N = M = 3, µ = 2 Important: P 1 E < P 2 E and P 2 O < P 2 I 2,3

45 Three Pursuers N = M = 3, µ = 2 Figure: Game with Three Pursuers M. Pachter (2018) Many-on-One Pursuit October 13, / 64

46 Three Pursuers N = M = 3, µ = 2 C 3 P 1 C 2 I 1, 3 E I 1, 2 P 2 O * I 2, 3 P 3 C 1 Figure: P 1, P 2, P 3 and E head Toward O M. Pachter (2018) Many-on-One Pursuit October 13, / 64

47 Three Pursuers Symmetry N = M = 3, µ = 5 4 Figure: E is at Circumcenter of Equilateral P 1 P 2 P 3 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

48 Three Pursuers Symmetry N = M = 3, µ = 5 4 (4).png Figure: Symmetric State M. Pachter (2018) Many-on-One Pursuit October 13, / 64

49 Three Pursuers Symmetry N = M = 3, µ = 5 4 Symmetry goes hand in hand with singularity Symmetry, and E heads South A state is reached where O SR and dist(p 2, V k ) = R Without pursuers dead zone parameterized by 0 < δ << 1, E will dither. P 2 & P 3 will zigzag Time-to-Capture > V = 1 µ R Chattering is anathema to optimality Inclusion of pursuers dead zone is crucial step: E won t dither Price paid: Slight reduction in optimality M. Pachter (2018) Many-on-One Pursuit October 13, / 64

50 Three Pursuers Summary N = M = 3 Optimal Solution 1. If E P 1 P 2 P 3, E captured by 2 active pursuers. 2. If E P 1 P 2 P 3 & P 1 P 2 P 3 obtuse, E captured by 2 active pursuers. 3. If E P 1 P 2 P 3 and P 1 P 2 P 3 is acute but vertex of BSR is outside hexagon, E captured by 2 active pursuers. Interesting scenario: 4. E P 1 P 2 P 3, P 1 P 2 P 3 is acute and SR Hexagon Pursuers team might incur a slight loss in optimality, BUT: Evader strategy with no dithering. Capture strategy of team of pursuers with no chatter. When N 3 the pursuit-evasion problem is much harder M. Pachter (2018) Many-on-One Pursuit October 13, / 64

51 Discussion Only when we have a complete set of optimal trajectories filling a capture region separated from the escape region by a closed barrier, or possibly, the region of capturability is the whole state space this, as specified by the solution of the Game of Kind has a differential game been solved. Only a few 3-D pursuit-evasion games solved, none in a higher dimension. Pursuit-evasion differential games in 3-D solved which are an acceptable model of conflict situations of interest so that their relevance is preserved: Games with no singular surfaces, at least for a range of parameters. Examples 2 Cutters & Fugitive Ship, Guarding a Target DG, Active Target Defense. The secret sauce is provided by the following insight. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

52 Target Defense Differential Game Three States, V D V A 1 Figure: Guarding a Target M. Pachter (2018) Many-on-One Pursuit October 13, / 64

53 Defending a Coastline: Infinite Size Target Three States, V D V A 1 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

54 Active Aircraft Defense Three States y T VT n D VD w A VA m x Figure: The Dynamic Target Defense Differential Game V D V A and the speed ratio µ = V A V T > 1. D and T form a team that plays against A. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

55 Discussion Theorem 3 The solution of zero-sum differential games by solving the open-loop minmax optimal control problem using the two sided Pontryagin Maximum Principle and synthesizing the players state feedback optimal strategies in receding horizon optimal control fashion is valid iff the application of Isaacs method results in primary optimal trajectories only, a.k.a., regular characteristics. When this is the case and the players have simple motion the geometric method applies. We have proven that the application of Isaacs method to the Two cutters and Fugitive Ship Differential Game yields an optimal flow field which covers the 3-D state space. This is the reason why in the Two Cutters and Fugitive Ship Differential Game the geometric method is applicable and the correctness of Isaacs solution has been proved. Optimal strategies provide true AI action M. Pachter (2018) Many-on-One Pursuit October 13, / 64

56 Extension: Slow Pursuer P 2 µ 2 = 1 Figure: Slow Pursuer is Redundant M. Pachter (2018) Many-on-One Pursuit October 13, / 64

57 Extension: Slow Pursuer P 2 µ 2 = 1 Figure: Capture Affected by Fast Pursuer Only M. Pachter (2018) Many-on-One Pursuit October 13, / 64

58 Extension: Active Slow Pursuer µ 2 = 1 Figure: Evader is Isochronously Captured by Fast and Slow Pursuer M. Pachter (2018) Many-on-One Pursuit October 13, / 64

59 Two Slow Pursuers Figure: Escape from Two Pursuers whose Speed Ratio µ = 1 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

60 Multiple Agents Network of Slow Pursuers All Agents Have Simple Motion & Equal Speed Point Capture (l = 0) Stroboscopic Information Pattern: Pshenichnyi E Convhull({P 1,..., P N }) E Escapes M. Pachter (2018) Many-on-One Pursuit October 13, / 64

61 Three Slow Pursuers Figure: Evader Surrounded by Pursuers whose Speed Ratio µ = 1 M. Pachter (2018) Many-on-One Pursuit October 13, / 64

62 References 1 Hugo Steinhaus: Definitions for a Theory of Games of Pursuit, Naval Research Logistics Quarterly, Vol. 7, No 2, pp , Rufus Isaacs: Differential Games: A Mathematical Theory with Applications to Warfare and Pursuit, Control and Optimization, Wiley 1965, pp E. Garcia, Z. Fuchs, D. Milutinovic, D. Casbeer, and M. Pachter: A Geometric Approach for the Cooperative Two-Pursuers One-Evader Differential Game, Proceedings of the 20th World Congress of IFAC, Toulouse, France, pp , July 9-14, M. Pachter: Isaacs Two-on-One Pursuit-Evasion Game, Submitted to the Annals of the ISDG. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

63 References 1 E. Garcia, Z. Fuchs, D. Milutinovic, D. Casbeer, and M. Pachter: A Geometric Approach for the Cooperative Two-Pursuers One-Evader Differential Game, Proceedings of the 20th World Congress of IFAC, Toulouse, France, pp , July 9-14, M. Pachter, A. VonMoll, E. Garcia, D. Milutinovic, D.Casbeer: Two-on-One Pursuit, submitted to the AIAA J. of Guidance, Control and Dynamics. 3 H. Hagedorn and J. V. Breakwell: A Differential Game with Two Pursuers and One Evader, JOTA, Vol. 18, No. 1, January 1976, pp B. N. Pshenichnyi: Simple Pursuit By Several Objects, Cybernetics, May-June 1976, pp Gotfried T. Toussaint: Computing Largest Empty Circles with Location Constraints, Int. J. Computer and Information Sciences, Vol. 12, No. 5, 1983, pp M. Pachter (2018) Many-on-One Pursuit October 13, / 64

64 References 1 A. A. Chikrii and P. V. Prokopovich: Simple Pursuit of One Evader by a Group, Cybernetics and Systems Analysis, No. 3, May-June 1992, pp (in Russian). 2 S. Jin and Z. Qu: Pursuit-Evasion with Multi-Pursuer v.s. One Fast Evader, Proceedings of the 8th World Congress on Intelligent Control and Automation, July , Jinan, China, pp G. I. Ibragimov: A Game of Optimal Pursuit of One Object by Several, J. Appl. Maths. and Mechs., Vol. 62, No. 2, pp , S. Kopparty and C. V. Ravishankar: A Framework for Pursuit Evasion Games in R n, Information Processing Letters 96 (2005), pp Sergey S. Kumkov, Stephane Le Menec and Valerii S. Patsko: Zero-Sum Pursuit-Evasion Differential Games with Many Objects: Survey of Publications, Dynamic Games And Applications, 16 November 2016, DOI /s z. M. Pachter (2018) Many-on-One Pursuit October 13, / 64

Nash Equilibria in a Group Pursuit Game

Nash Equilibria in a Group Pursuit Game Applied Mathematical Sciences, Vol. 10, 2016, no. 17, 809-821 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2016.614 Nash Equilibria in a Group Pursuit Game Yaroslavna Pankratova St. Petersburg

More information

History of the Homicidal Chauffeur Problem

History of the Homicidal Chauffeur Problem History of the Homicidal Chauffeur Problem Dedicated to Pierre Bernhard on his 60th birthday Statement of the problem by R.Isaacs in the early fifties 3 R.Isaacs Rufus Isaacs has formulated his famous

More information

Phenomenon of Narrow Throats of Level Sets of Value Function in Differential Games

Phenomenon of Narrow Throats of Level Sets of Value Function in Differential Games Phenomenon of Narrow Throats of Level Sets of Value Function in Differential Games Sergey S. Kumkov 1 and Valerii S. Patsko 1 Institute of Mathematics and Mechanics, S.Kovalevskaya str., 16, Ekaterinburg,

More information

Mathematical Optimization Models and Applications

Mathematical Optimization Models and Applications Mathematical Optimization Models and Applications Yinyu Ye Department of Management Science and Engineering Stanford University Stanford, CA 94305, U.S.A. http://www.stanford.edu/ yyye Chapters 1, 2.1-2,

More information

A SOLUTION OF MISSILE-AIRCRAFT PURSUIT-EVASION DIFFERENTIAL GAMES IN MEDIUM RANGE

A SOLUTION OF MISSILE-AIRCRAFT PURSUIT-EVASION DIFFERENTIAL GAMES IN MEDIUM RANGE Proceedings of the 1th WSEAS International Confenrence on APPLIED MAHEMAICS, Dallas, exas, USA, November 1-3, 6 376 A SOLUION OF MISSILE-AIRCRAF PURSUI-EVASION DIFFERENIAL GAMES IN MEDIUM RANGE FUMIAKI

More information

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

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

More information

Dynamic Games. Meir Pachter. March 11, Department of Electrical and Computer Engineering. Air Force Institute of Technology/AFIT

Dynamic Games. Meir Pachter. March 11, Department of Electrical and Computer Engineering. Air Force Institute of Technology/AFIT Dynamic Games Meir achter March, 2009 Department of Electrical and Computer Engineering Air Force Institute of Technology/AFIT Wright atterson AFB, OH 45433 Introduction Objective: Include adversarial

More information

Pursuit, herding and evasion: A three-agent model of caribou predation

Pursuit, herding and evasion: A three-agent model of caribou predation Proceedings of the American Control Conference, Washington, DC, 2013 Pursuit, herding and evasion: A three-agent model of caribou predation William Scott and Naomi Ehrich Leonard Abstract Motivated by

More information

Suboptimality of minmax MPC. Seungho Lee. ẋ(t) = f(x(t), u(t)), x(0) = x 0, t 0 (1)

Suboptimality of minmax MPC. Seungho Lee. ẋ(t) = f(x(t), u(t)), x(0) = x 0, t 0 (1) Suboptimality of minmax MPC Seungho Lee In this paper, we consider particular case of Model Predictive Control (MPC) when the problem that needs to be solved in each sample time is the form of min max

More information

Fixed Duration Pursuit-Evasion Differential Game with Integral Constraints

Fixed Duration Pursuit-Evasion Differential Game with Integral Constraints Journal of Physics: Conference Series OPEN ACCESS Fixed Duration Pursuit-Evasion Differential Game with Integral Constraints To cite this article: Ibragimov G I and Kuchkarov A Sh 213 J. Phys.: Conf. Ser.

More information

arxiv: v1 [math.oc] 30 Apr 2015

arxiv: v1 [math.oc] 30 Apr 2015 arxiv:155.54v1 [math.oc] 3 Apr 215 DIFFERENTIAL GAME OF MANY PURSUERS WITH INTEGRAL CONSTRAINTS ON A CONVEX SET IN THE PLANE IDHAM ARIF ALIAS, GAFURJAN IBRAGIMOV, MASSIMILIANO FERRARA, MEHDI SALIMI, AND

More information

pacman A classical arcade videogame powered by Hamilton-Jacobi equations Simone Cacace Universita degli Studi Roma Tre

pacman A classical arcade videogame powered by Hamilton-Jacobi equations Simone Cacace Universita degli Studi Roma Tre HJ pacman A classical arcade videogame powered by Hamilton-Jacobi equations Simone Cacace Universita degli Studi Roma Tre Control of State Constrained Dynamical Systems 25-29 September 2017, Padova Main

More information

An infinite dimensional pursuit evasion differential game with finite number of players

An infinite dimensional pursuit evasion differential game with finite number of players An infinite dimensional pursuit evasion differential game with finite number of players arxiv:156.919v1 [math.oc] 2 Jun 215 Mehdi Salimi a Massimiliano Ferrara b a Center for Dynamics, Department of Mathematics,

More information

Reach Sets and the Hamilton-Jacobi Equation

Reach Sets and the Hamilton-Jacobi Equation Reach Sets and the Hamilton-Jacobi Equation Ian Mitchell Department of Computer Science The University of British Columbia Joint work with Alex Bayen, Meeko Oishi & Claire Tomlin (Stanford) research supported

More information

I ve Got a Three-Body Problem

I ve Got a Three-Body Problem I ve Got a Three-Body Problem Gareth E. Roberts Department of Mathematics and Computer Science College of the Holy Cross Mathematics Colloquium Fitchburg State College November 13, 2008 Roberts (Holy Cross)

More information

Dynamic and Adversarial Reachavoid Symbolic Planning

Dynamic and Adversarial Reachavoid Symbolic Planning Dynamic and Adversarial Reachavoid Symbolic Planning Laya Shamgah Advisor: Dr. Karimoddini July 21 st 2017 Thrust 1: Modeling, Analysis and Control of Large-scale Autonomous Vehicles (MACLAV) Sub-trust

More information

Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems

Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems 1 Decentralized Stabilization of Heterogeneous Linear Multi-Agent Systems Mauro Franceschelli, Andrea Gasparri, Alessandro Giua, and Giovanni Ulivi Abstract In this paper the formation stabilization problem

More information

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Anton Shiriaev 1,2, Leonid Freidovich 1, Sergey Gusev 3 1 Department

More information

Online Appendix for Coordination of Outsourced Operations at a Third-Party Facility Subject to Booking, Overtime, and Tardiness Costs

Online Appendix for Coordination of Outsourced Operations at a Third-Party Facility Subject to Booking, Overtime, and Tardiness Costs Submitted to Operations Research manuscript OPRE-2009-04-180 Online Appendix for Coordination of Outsourced Operations at a Third-Party Facility Subject to Booking, Overtime, and Tardiness Costs Xiaoqiang

More information

Review Notes for IB Standard Level Math

Review Notes for IB Standard Level Math Review Notes for IB Standard Level Math 1 Contents 1 Algebra 8 1.1 Rules of Basic Operations............................... 8 1.2 Rules of Roots..................................... 8 1.3 Rules of Exponents...................................

More information

Alberto Bressan. Department of Mathematics, Penn State University

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

More information

Smooth Path Generation Based on Bézier Curves for Autonomous Vehicles

Smooth Path Generation Based on Bézier Curves for Autonomous Vehicles Smooth Path Generation Based on Bézier Curves for Autonomous Vehicles Ji-wung Choi, Renwick E. Curry, Gabriel Hugh Elkaim Abstract In this paper we present two path planning algorithms based on Bézier

More information

United Arab Emirates University

United Arab Emirates University United Arab Emirates University University Foundation Program - Math Program ALGEBRA - COLLEGE ALGEBRA - TRIGONOMETRY Practice Questions 1. What is 2x 1 if 4x + 8 = 6 + x? A. 2 B. C. D. 4 E. 2. What is

More information

arxiv:math/ v1 [math.ds] 8 Oct 2006

arxiv:math/ v1 [math.ds] 8 Oct 2006 arxiv:math/0610257v1 [math.ds] 8 Oct 2006 On Estimates of the Number of Collisions for Billiards in Polyhedral Angles Lizhou Chen Abstract We obtain an upper bound of the number of collisions of any billiard

More information

minimize x subject to (x 2)(x 4) u,

minimize x subject to (x 2)(x 4) u, Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for

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

An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game

An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game An Application of Pontryagin s Maximum Principle in a Linear Quadratic Differential Game Marzieh Khakestari (Corresponding author) Institute For Mathematical Research, Universiti Putra Malaysia, 43400

More information

CSE 554 Lecture 7: Alignment

CSE 554 Lecture 7: Alignment CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex

More information

Multi-Pursuer Single-Evader Differential Games with Limited Observations*

Multi-Pursuer Single-Evader Differential Games with Limited Observations* American Control Conference (ACC) Washington DC USA June 7-9 Multi-Pursuer Single- Differential Games with Limited Observations* Wei Lin Student Member IEEE Zhihua Qu Fellow IEEE and Marwan A. Simaan Life

More information

arxiv: v1 [cs.sy] 29 Sep 2017

arxiv: v1 [cs.sy] 29 Sep 2017 Beacon-referenced Mutual Pursuit in Three Dimensions Kevin S. Galloway and Biswadip Dey arxiv:709.037v cs.sy 9 Sep 07 Abstract Motivated by station-keeping applications in various unmanned settings, this

More information

Optimization: Insights and Applications. Jan Brinkhuis Vladimir Tikhomirov PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD

Optimization: Insights and Applications. Jan Brinkhuis Vladimir Tikhomirov PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Optimization: Insights and Applications Jan Brinkhuis Vladimir Tikhomirov PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Preface 0.1 Optimization: insights and applications xiii 0.2 Lunch, dinner,

More information

Multi-Robotic Systems

Multi-Robotic Systems CHAPTER 9 Multi-Robotic Systems The topic of multi-robotic systems is quite popular now. It is believed that such systems can have the following benefits: Improved performance ( winning by numbers ) Distributed

More information

DRAFT Formulation and Analysis of Linear Programs

DRAFT Formulation and Analysis of Linear Programs DRAFT Formulation and Analysis of Linear Programs Benjamin Van Roy and Kahn Mason c Benjamin Van Roy and Kahn Mason September 26, 2005 1 2 Contents 1 Introduction 7 1.1 Linear Algebra..........................

More information

Wheelchair Collision Avoidance: An Application for Differential Games

Wheelchair Collision Avoidance: An Application for Differential Games Wheelchair Collision Avoidance: An Application for Differential Games Pooja Viswanathan December 27, 2006 Abstract This paper discusses the design of intelligent wheelchairs that avoid obstacles using

More information

Intruder Isolation on a General Road Network under Partial Information

Intruder Isolation on a General Road Network under Partial Information Intruder Isolation on a General Road Network under Partial Information Hua Chen, Krishna Kalyanam, Wei Zhang and David Casbeer 1 Abstract The intruder isolation problem (IIP) involves using an unmanned

More information

Georgia Tech High School Math Competition

Georgia Tech High School Math Competition Georgia Tech High School Math Competition Multiple Choice Test February 28, 2015 Each correct answer is worth one point; there is no deduction for incorrect answers. Make sure to enter your ID number on

More information

(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.

(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2. 1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of

More information

State Math Contest Senior Exam SOLUTIONS

State Math Contest Senior Exam SOLUTIONS State Math Contest Senior Exam SOLUTIONS 1. The following pictures show two views of a non standard die (however the numbers 1-6 are represented on the die). How many dots are on the bottom face of figure?

More information

How curvature shapes space

How curvature shapes space How curvature shapes space Richard Schoen University of California, Irvine - Hopf Lecture, ETH, Zürich - October 30, 2017 The lecture will have three parts: Part 1: Heinz Hopf and Riemannian geometry Part

More information

For math conventions used on the GRE, refer to this link:

For math conventions used on the GRE, refer to this link: GRE Review ISU Student Success Center Quantitative Workshop One Quantitative Section: Overview Your test will include either two or three 35-minute quantitative sections. There will be 20 questions in

More information

High School Math Contest

High School Math Contest High School Math Contest University of South Carolina February th, 017 Problem 1. If (x y) = 11 and (x + y) = 169, what is xy? (a) 11 (b) 1 (c) 1 (d) (e) 8 Solution: Note that xy = (x + y) (x y) = 169

More information

The Locations of the Brocard Points

The Locations of the Brocard Points Forum Geometricorum Volume 6 (2006) 71 77. FORUM GEOM ISSN 1534-1178 The Locations of the Brocard Points Christopher J. Bradley and Geoff C. Smith Abstract. We fix and scale to constant size the orthocentroidal

More information

Control of Mobile Robots

Control of Mobile Robots Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and

More information

Stability of the two queue system

Stability of the two queue system Stability of the two queue system Iain M. MacPhee and Lisa J. Müller University of Durham Department of Mathematical Science Durham, DH1 3LE, UK (e-mail: i.m.macphee@durham.ac.uk, l.j.muller@durham.ac.uk)

More information

Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations

Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations Martino Bardi Italo Capuzzo-Dolcetta Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations Birkhauser Boston Basel Berlin Contents Preface Basic notations xi xv Chapter I. Outline

More information

A Residual Gradient Fuzzy Reinforcement Learning Algorithm for Differential Games

A Residual Gradient Fuzzy Reinforcement Learning Algorithm for Differential Games International Journal of Fuzzy Systems manuscript (will be inserted by the editor) A Residual Gradient Fuzzy Reinforcement Learning Algorithm for Differential Games Mostafa D Awheda Howard M Schwartz Received:

More information

Differential Game Missile Guidance with Impact Angle and Time Constraints

Differential Game Missile Guidance with Impact Angle and Time Constraints Differential Game Missile Guidance with Impact Angle and Time Constraints Seonhyeok Kang H. Jin Kim School of Mechanical and Aerospace Engineering, Seoul National University, Seoul, Korea, (e-mail: hkkang@snu.ac.kr)

More information

Energy-based Swing-up of the Acrobot and Time-optimal Motion

Energy-based Swing-up of the Acrobot and Time-optimal Motion Energy-based Swing-up of the Acrobot and Time-optimal Motion Ravi N. Banavar Systems and Control Engineering Indian Institute of Technology, Bombay Mumbai-476, India Email: banavar@ee.iitb.ac.in Telephone:(91)-(22)

More information

An Optimal Evader Strategy in a Two-Pursuer One-Evader Problem

An Optimal Evader Strategy in a Two-Pursuer One-Evader Problem An Optimal Evader Strateg in a Two-Pursuer One-Evader Problem Wei Sun and Panagiotis Tsiotras 2 Abstract We consider a rela pursuit-evasion problem with two pursuers and one evader. We reduce the problem

More information

Submitted to IEEE TR, Nov 10, 2007 On Discrete-time Pursuit-evasion Games with Sensing Limitations

Submitted to IEEE TR, Nov 10, 2007 On Discrete-time Pursuit-evasion Games with Sensing Limitations Submitted to IEEE TR Nov 10 007 On Discrete-time Pursuit-evasion Games with Sensing Limitations Shaunak D. Bopardikar (Corresponding author Department of Mechanical Engineering University of California

More information

Algebra One Dictionary

Algebra One Dictionary Algebra One Dictionary Page 1 of 17 A Absolute Value - the distance between the number and 0 on a number line Algebraic Expression - An expression that contains numbers, operations and at least one variable.

More information

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013

Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard June 15, 2013 Abstract As in optimal control theory, linear quadratic (LQ) differential games (DG) can be solved, even in high dimension,

More information

Integrals. D. DeTurck. January 1, University of Pennsylvania. D. DeTurck Math A: Integrals 1 / 61

Integrals. D. DeTurck. January 1, University of Pennsylvania. D. DeTurck Math A: Integrals 1 / 61 Integrals D. DeTurck University of Pennsylvania January 1, 2018 D. DeTurck Math 104 002 2018A: Integrals 1 / 61 Integrals Start with dx this means a little bit of x or a little change in x If we add up

More information

14 th Annual Harvard-MIT Mathematics Tournament Saturday 12 February 2011

14 th Annual Harvard-MIT Mathematics Tournament Saturday 12 February 2011 14 th Annual Harvard-MIT Mathematics Tournament Saturday 1 February 011 1. Let a, b, and c be positive real numbers. etermine the largest total number of real roots that the following three polynomials

More information

Finding Optimal Strategies for Influencing Social Networks in Two Player Games. MAJ Nick Howard, USMA Dr. Steve Kolitz, Draper Labs Itai Ashlagi, MIT

Finding Optimal Strategies for Influencing Social Networks in Two Player Games. MAJ Nick Howard, USMA Dr. Steve Kolitz, Draper Labs Itai Ashlagi, MIT Finding Optimal Strategies for Influencing Social Networks in Two Player Games MAJ Nick Howard, USMA Dr. Steve Kolitz, Draper Labs Itai Ashlagi, MIT Problem Statement Given constrained resources for influencing

More information

A Receding Horizon Control Law for Harbor Defense

A Receding Horizon Control Law for Harbor Defense A Receding Horizon Control Law for Harbor Defense Seungho Lee 1, Elijah Polak 2, Jean Walrand 2 1 Department of Mechanical Science and Engineering University of Illinois, Urbana, Illinois - 61801 2 Department

More information

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 31 OCTOBER 2018

SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY 31 OCTOBER 2018 SOLUTIONS FOR ADMISSIONS TEST IN MATHEMATICS, COMPUTER SCIENCE AND JOINT SCHOOLS WEDNESDAY OCTOBER 8 Mark Scheme: Each part of Question is worth marks which are awarded solely for the correct answer. Each

More information

Individual Round CHMMC November 20, 2016

Individual Round CHMMC November 20, 2016 Individual Round CHMMC 20 November 20, 20 Problem. We say that d k d k d d 0 represents the number n in base 2 if each d i is either 0 or, and n d k ( 2) k + d k ( 2) k + + d ( 2) + d 0. For example, 0

More information

Pursuit Differential Game Described by Infinite First Order 2-Systems of Differential Equations

Pursuit Differential Game Described by Infinite First Order 2-Systems of Differential Equations Malaysian Journal of Mathematical Sciences 11(2): 181 19 (217) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Journal homepage: http://einspem.upm.edu.my/journal Pursuit Differential Game Described by Infinite

More information

Pontryagin s Minimum Principle 1

Pontryagin s Minimum Principle 1 ECE 680 Fall 2013 Pontryagin s Minimum Principle 1 In this handout, we provide a derivation of the minimum principle of Pontryagin, which is a generalization of the Euler-Lagrange equations that also includes

More information

Hyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems. Joanna Zwierzchowska

Hyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems. Joanna Zwierzchowska Decision Making in Manufacturing and Services Vol. 9 05 No. pp. 89 00 Hyperbolicity of Systems Describing Value Functions in Differential Games which Model Duopoly Problems Joanna Zwierzchowska Abstract.

More information

ON CALCULATING THE VALUE OF A DIFFERENTIAL GAME IN THE CLASS OF COUNTER STRATEGIES 1,2

ON CALCULATING THE VALUE OF A DIFFERENTIAL GAME IN THE CLASS OF COUNTER STRATEGIES 1,2 URAL MATHEMATICAL JOURNAL, Vol. 2, No. 1, 2016 ON CALCULATING THE VALUE OF A DIFFERENTIAL GAME IN THE CLASS OF COUNTER STRATEGIES 1,2 Mikhail I. Gomoyunov Krasovskii Institute of Mathematics and Mechanics,

More information

DSCC TIME-OPTIMAL TRAJECTORIES FOR STEERED AGENT WITH CONSTRAINTS ON SPEED AND TURNING RATE

DSCC TIME-OPTIMAL TRAJECTORIES FOR STEERED AGENT WITH CONSTRAINTS ON SPEED AND TURNING RATE Proceedings of the ASME 216 Dynamic Systems and Control Conference DSCC216 October 12-14, 216, Minneapolis, Minnesota, USA DSCC216-9892 TIME-OPTIMAL TRAJECTORIES FOR STEERED AGENT WITH CONSTRAINTS ON SPEED

More information

2018 Best Student Exam Solutions Texas A&M High School Students Contest October 20, 2018

2018 Best Student Exam Solutions Texas A&M High School Students Contest October 20, 2018 08 Best Student Exam Solutions Texas A&M High School Students Contest October 0, 08. You purchase a stock and later sell it for $44 per share. When you do, you notice that the percent increase was the

More information

Some Properties of the Augmented Lagrangian in Cone Constrained Optimization

Some Properties of the Augmented Lagrangian in Cone Constrained Optimization MATHEMATICS OF OPERATIONS RESEARCH Vol. 29, No. 3, August 2004, pp. 479 491 issn 0364-765X eissn 1526-5471 04 2903 0479 informs doi 10.1287/moor.1040.0103 2004 INFORMS Some Properties of the Augmented

More information

Integrated Math II. IM2.1.2 Interpret given situations as functions in graphs, formulas, and words.

Integrated Math II. IM2.1.2 Interpret given situations as functions in graphs, formulas, and words. Standard 1: Algebra and Functions Students graph linear inequalities in two variables and quadratics. They model data with linear equations. IM2.1.1 Graph a linear inequality in two variables. IM2.1.2

More information

Lecture 2 The Core. 2.1 Definition and graphical representation for games with up to three players

Lecture 2 The Core. 2.1 Definition and graphical representation for games with up to three players Lecture 2 The Core Let us assume that we have a TU game (N, v) and that we want to form the grand coalition. The Core, which was first introduced by Gillies [1], is the most attractive and natural way

More information

1 Implicit Differentiation

1 Implicit Differentiation MATH 1010E University Mathematics Lecture Notes (week 5) Martin Li 1 Implicit Differentiation Sometimes a function is defined implicitly by an equation of the form f(x, y) = 0, which we think of as a relationship

More information

Application of dynamic programming approach to aircraft take-off in a windshear

Application of dynamic programming approach to aircraft take-off in a windshear Application of dynamic programming approach to aircraft take-off in a windshear N. D. Botkin, V. L. Turova Technische Universität München, Mathematics Centre Chair of Mathematical Modelling 10th International

More information

Min-Max Certainty Equivalence Principle and Differential Games

Min-Max Certainty Equivalence Principle and Differential Games Min-Max Certainty Equivalence Principle and Differential Games Pierre Bernhard and Alain Rapaport INRIA Sophia-Antipolis August 1994 Abstract This paper presents a version of the Certainty Equivalence

More information

MATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH,

MATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH, MATH 189, MATHEMATICAL METHODS IN CLASSICAL AND QUANTUM MECHANICS. HOMEWORK 2. DUE WEDNESDAY OCTOBER 8TH, 2014 HTTP://MATH.BERKELEY.EDU/~LIBLAND/MATH-189/HOMEWORK.HTML 1. Lagrange Multipliers (Easy) Exercise

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

Quadratic and Polynomial Inequalities in one variable have look like the example below.

Quadratic and Polynomial Inequalities in one variable have look like the example below. Section 8 4: Polynomial Inequalities in One Variable Quadratic and Polynomial Inequalities in one variable have look like the example below. x 2 5x 6 0 (x 2) (x + 4) > 0 x 2 (x 3) > 0 (x 2) 2 (x + 4) 0

More information

Asteroid Rendezvous Missions

Asteroid Rendezvous Missions Asteroid Rendezvous Missions An Application of Optimal Control G. Patterson, G. Picot, and S. Brelsford University of Hawai`i at Manoa 25 June, 215 1/5 G. Patterson, G. Picot, and S. Brelsford Asteroid

More information

EN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018

EN Applied Optimal Control Lecture 8: Dynamic Programming October 10, 2018 EN530.603 Applied Optimal Control Lecture 8: Dynamic Programming October 0, 08 Lecturer: Marin Kobilarov Dynamic Programming (DP) is conerned with the computation of an optimal policy, i.e. an optimal

More information

Classical Mechanics Review (Louisiana State University Qualifier Exam)

Classical Mechanics Review (Louisiana State University Qualifier Exam) Review Louisiana State University Qualifier Exam Jeff Kissel October 22, 2006 A particle of mass m. at rest initially, slides without friction on a wedge of angle θ and and mass M that can move without

More information

Noncooperative Games, Couplings Constraints, and Partial Effi ciency

Noncooperative Games, Couplings Constraints, and Partial Effi ciency Noncooperative Games, Couplings Constraints, and Partial Effi ciency Sjur Didrik Flåm University of Bergen, Norway Background Customary Nash equilibrium has no coupling constraints. Here: coupling constraints

More information

Cooperative Pursuit with Voronoi Partitions

Cooperative Pursuit with Voronoi Partitions Cooperative Pursuit with Voronoi Partitions Zhengyuan Zhou a, Wei Zhang b, Jerry Ding c, Haomiao Huang d, Dušan M. Stipanović e, Claire Tomlin f a Department of Electrical Engineering, Stanford University,

More information

Stability Analysis of Optimal Adaptive Control under Value Iteration using a Stabilizing Initial Policy

Stability Analysis of Optimal Adaptive Control under Value Iteration using a Stabilizing Initial Policy Stability Analysis of Optimal Adaptive Control under Value Iteration using a Stabilizing Initial Policy Ali Heydari, Member, IEEE Abstract Adaptive optimal control using value iteration initiated from

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

Linear Quadratic Zero-Sum Two-Person Differential Games

Linear Quadratic Zero-Sum Two-Person Differential Games Linear Quadratic Zero-Sum Two-Person Differential Games Pierre Bernhard To cite this version: Pierre Bernhard. Linear Quadratic Zero-Sum Two-Person Differential Games. Encyclopaedia of Systems and Control,

More information

Navigation and Obstacle Avoidance via Backstepping for Mechanical Systems with Drift in the Closed Loop

Navigation and Obstacle Avoidance via Backstepping for Mechanical Systems with Drift in the Closed Loop Navigation and Obstacle Avoidance via Backstepping for Mechanical Systems with Drift in the Closed Loop Jan Maximilian Montenbruck, Mathias Bürger, Frank Allgöwer Abstract We study backstepping controllers

More information

Geometry Final Exam Review

Geometry Final Exam Review Name: Date: Period: Geometry Final Exam Review 1. Fill in the flow chart below with the properties that belong to each polygon. 2. Find the measure of each numbered angle: 3. Find the value of x 4. Calculate

More information

Formally Analyzing Adaptive Flight Control

Formally Analyzing Adaptive Flight Control Formally Analyzing Adaptive Flight Control Ashish Tiwari SRI International 333 Ravenswood Ave Menlo Park, CA 94025 Supported in part by NASA IRAC NRA grant number: NNX08AB95A Ashish Tiwari Symbolic Verification

More information

A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations

A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations 2005 American Control Conference June 8-10, 2005. Portland, OR, USA WeB10.3 A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations Arik Melikyan, Andrei Akhmetzhanov and Naira

More information

3. Which of these numbers does not belong to the set of solutions of the inequality 4

3. Which of these numbers does not belong to the set of solutions of the inequality 4 Math Field Day Exam 08 Page. The number is equal to b) c) d) e). Consider the equation 0. The slope of this line is / b) / c) / d) / e) None listed.. Which of these numbers does not belong to the set of

More information

Official Solutions 2014 Sun Life Financial CMO Qualifying Rêpechage 1

Official Solutions 2014 Sun Life Financial CMO Qualifying Rêpechage 1 Official s 2014 Sun Life Financial CMO Qualifying Rêpechage 1 1. Let f : Z Z + be a function, and define h : Z Z Z + by h(x, y) = gcd(f(x), f(y)). If h(x, y) is a two-variable polynomial in x and y, prove

More information

MATH 32 FALL 2012 FINAL EXAM - PRACTICE EXAM SOLUTIONS

MATH 32 FALL 2012 FINAL EXAM - PRACTICE EXAM SOLUTIONS MATH 3 FALL 0 FINAL EXAM - PRACTICE EXAM SOLUTIONS () You cut a slice from a circular pizza (centered at the origin) with radius 6 along radii at angles 4 and 3 with the positive horizontal axis. (a) (3

More information

Optimal Control. McGill COMP 765 Oct 3 rd, 2017

Optimal Control. McGill COMP 765 Oct 3 rd, 2017 Optimal Control McGill COMP 765 Oct 3 rd, 2017 Classical Control Quiz Question 1: Can a PID controller be used to balance an inverted pendulum: A) That starts upright? B) That must be swung-up (perhaps

More information

Pontryagin s maximum principle

Pontryagin s maximum principle Pontryagin s maximum principle Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2012 Emo Todorov (UW) AMATH/CSE 579, Winter 2012 Lecture 5 1 / 9 Pontryagin

More information

ON LEVEL SETS WITH NARROW THROATS IN LINEAR DIFFERENTIAL GAMES

ON LEVEL SETS WITH NARROW THROATS IN LINEAR DIFFERENTIAL GAMES Fourth International ISDG Workshop, May 19 to May 21, 2003, Goslar Proceedings edited by Michael H. Breitner, pp. 89 114 ON LEVEL SETS WITH NARROW THROATS IN LINEAR DIFFERENTIAL GAMES SERGEY S. KUMKOV,

More information

8th Grade Math Definitions

8th Grade Math Definitions 8th Grade Math Definitions Absolute Value: 1. A number s distance from zero. 2. For any x, is defined as follows: x = x, if x < 0; x, if x 0. Acute Angle: An angle whose measure is greater than 0 and less

More information

Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization

Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012

More information

CS188: Artificial Intelligence, Fall 2009 Written 2: MDPs, RL, and Probability

CS188: Artificial Intelligence, Fall 2009 Written 2: MDPs, RL, and Probability CS188: Artificial Intelligence, Fall 2009 Written 2: MDPs, RL, and Probability Due: Thursday 10/15 in 283 Soda Drop Box by 11:59pm (no slip days) Policy: Can be solved in groups (acknowledge collaborators)

More information

Robust Predictions in Games with Incomplete Information

Robust Predictions in Games with Incomplete Information Robust Predictions in Games with Incomplete Information joint with Stephen Morris (Princeton University) November 2010 Payoff Environment in games with incomplete information, the agents are uncertain

More information

INVERSION IN INDIRECT OPTIMAL CONTROL

INVERSION IN INDIRECT OPTIMAL CONTROL INVERSION IN INDIRECT OPTIMAL CONTROL François Chaplais, Nicolas Petit Centre Automatique et Systèmes, École Nationale Supérieure des Mines de Paris, 35, rue Saint-Honoré 7735 Fontainebleau Cedex, France,

More information

Uniqueness of Optimal Cube Wrapping

Uniqueness of Optimal Cube Wrapping CCCG 014, Halifax, Nova Scotia, August 11 13, 014 Uniqueness of Optimal Cube Wrapping Qinxuan Pan EECS, MIT qinxuan@mit.edu Abstract Consider wrapping the unit cube with a square without stretching or

More information

Outline. Roadmap for the NPP segment: 1 Preliminaries: role of convexity. 2 Existence of a solution

Outline. Roadmap for the NPP segment: 1 Preliminaries: role of convexity. 2 Existence of a solution Outline Roadmap for the NPP segment: 1 Preliminaries: role of convexity 2 Existence of a solution 3 Necessary conditions for a solution: inequality constraints 4 The constraint qualification 5 The Lagrangian

More information

MEAM 520. More Velocity Kinematics

MEAM 520. More Velocity Kinematics MEAM 520 More Velocity Kinematics Katherine J. Kuchenbecker, Ph.D. General Robotics, Automation, Sensing, and Perception Lab (GRASP) MEAM Department, SEAS, University of Pennsylvania Lecture 12: October

More information

Glossary. Glossary 981. Hawkes Learning Systems. All rights reserved.

Glossary. Glossary 981. Hawkes Learning Systems. All rights reserved. A Glossary Absolute value The distance a number is from 0 on a number line Acute angle An angle whose measure is between 0 and 90 Addends The numbers being added in an addition problem Addition principle

More information