Distributed Computation of Minimum Time Consensus for Multi-Agent Systems

Size: px
Start display at page:

Download "Distributed Computation of Minimum Time Consensus for Multi-Agent Systems"

Transcription

1 Distributed Computation of Minimum Time Consensus for Multi-Agent Systems Ameer Mulla, Deepak Patil, Debraj Chakraborty IIT Bombay, India Conference on Decision and Control 2014 Los Angeles, California, USA December 9, / 21

2 Minimum time consensus There are N double integrator agents with bounded inputs. Goal: Achieve consensus in minimum time. Strategy: Computation (in advance) of the consensus state. Key result: Computation can be distributed (By Helly s Theorem). Tool: Attainability set. 2 / 21

3 Problem Definition Outline 1 Minimum time consensus Problem Definition 2 Single integrators 3 / 21

4 Problem Definition Problem Statement Assume N identical agents with double integrator dynamics [ ] [ ] ẋ i (t) = x 0 0 i (t) + u 1 i (t) }{{}}{{} A b i = 1,...,N [ ri (t) with x i (t) = v i (t) Problem ] [ ri0, x i (0) = x i0 = v i0 Find x and min t such that, for all i, j x i (t) x j (t) 0 as t t x i ( t) = x j ( t) = x and x i (t) = x j (t) for all t t ] and u i (t) 1. 4 / 21

5 Problem Definition Definitions Attainable Set from point p at time t The set of all points that an agent can reach from initial condition p at time t using admissible input i.e. u i (t) 1 { t } A p (t) = x : x = e At p + e A(t τ) bu(τ)dτ, u(t) : u(t) 1 0 [ 1 t For double integrator e At = 0 1 Boundary of attainability set ts t x(t) = e At p ± 0 e A(t τ) Bdτ e A(t τ) Bdτ ; 0 t s t < t s Quadratic in t and t s ] 5 / 21

6 Problem Definition Definitions Figure: A p (t) 6 / 21

7 Outline 1 Minimum time consensus Problem Definition 2 Single integrators 7 / 21

8 Way to Consensus Possible if attainable sets of all agents at time t intersect Let t be minimum t such that A i(t) φ (A i (t):=a xi0 (t)) 1 i N Requires solution of coupled polynomial equations and inequalities Computation cannot be distributed 8 / 21

9 Main tool:helly s Theorem Helly s Theorem F a finite family of convex sets in R n Every n + 1 members of F have a point in common all the members of F have a point in common A i (t) R 2 is convex for all i and t Attainability set of every triplet intersect attainability set of all intersect Allows to distribute computation for all ( N 3) triplets evenly among N agents 9 / 21

10 Way to Consensus Lemma Let t ijk : Minimum time to consensus for agents {a i, a j, a k } Then, t max t ijk 1 i,j,k N If A i (t ) A j (t ) A k (t ) φ for some t > 0, then A i (t) A j (t) A k (t) φ for all t t Thus, t = max t ijk 1 i,j,k N Computing closed form expressions for three agents case suffices 10 / 21

11 Consensus: Three agents Theorem Let {a p, a q, a r } be the triple of agents and t pqr = The minimum time to consensus t = t pqr The corresponding consensus point x = x pqr Theorem For two different {a p, a q, a r } and {a p, a q, a r } If t pqr = t p q r = max t ijk then x pqr = x p q r 1 i,j,k N max 1 i,j,k N t ijk. 11 / 21

12 Consensus: Three agents Two ways towards consensus of three agents For a triple {a i, a j, a k }, t ijk max{ t ij, t jk, t ik } WLOG let t ij = max{ t ij, t jk, t ik }. Two cases: Case 1: t ijk = t ij = max{ t ij, t jk, t ik } i.e x ij A k ( t ij ) Case 2: t ijk > max{ t ij, t jk, t ik } i.e. x ij / A k ( t ij ) 12 / 21

13 Consensus: Three agents Case 1: t ijk = t ij = max{ t ij, t jk, t ik } i.e x ij A k ( t ij ) Figure: Case 1 13 / 21

14 Consensus: Three agents Case 2: t ijk > max{ t ij, t jk, t ik } i.e. x ij / A k ( t ij ) Figure: Case 2: x ij / A k ( t ij ) 14 / 21

15 Consensus: Three agents Case 2: t ijk > max{ t ij, t jk, t ik } i.e. x ij / A k ( t ij ) Figure: Case 2: Consensus of three agents 15 / 21

16 Role of Agents All N agents do not play a role in computation of t and x Triplet{a i, a j, a k } for which t = t ijk determine x and t For these agents, x A i ( t) Must use minimum-time control law to reach x at t In other words, x i0, x j0, x k0 R x ( t) For all other agents (say a q ), x inta q ( t), Must use a control law that allow them to reach x at t May use scaled down version of the input bounds 16 / 21

17 Six Agents Consensus Figure: Attainable Sets of Agents at t 17 / 21

18 Six Agents Consensus Definition The set of all initial conditions from which an agent can reach p at time t using admissible input i.e. u i (t) 1 Figure: Reachable Set to x at t 18 / 21

19 Single integrators Minimum time consensus: Single integrator N agents with single integrator dynamics ẋ i = u i, x i R and u i 1 Attainable set of agents are intervals in R given by x i (t) = x(0) + t 0 u i(t)dt Initial condition p then at time t attainable set is [p t,p + t] If every pair of the intervals intersect, then all the intervals intersect The minimum time is t ab = max ij t ij = max ij x i0 x j0 2 for some agents a and b (endpoints) That is t = max ij x i0 x j0 2 and x = max ij (x i0 +x j0 ) 2. For all other agents, any input that drives x i (t) to x in time t t ab For example, u i (t) = sign( x x i ) 19 / 21

20 Single integrators More can be said Single integrators interacting over a connected undirected graph G = (V,E) Prior computation of the t and corresponding x is not required Control law u i (t) = sign( (i,j) E (x j x i )) ensures proper input Agreement dynamics: ṗ(t) = sign(lp(t)), where L is Laplacian matrix Consensus occurs at x = max ij (x i0 +x j0 ) 2 20 / 21

21 Single integrators Conclusions An algorithm for computing minimum time consensus Computation is to be done for each possible triplet. Can be distributed evenly among each agent Final values of x and t are broad-casted to all Agreement becomes easier for single integrator agents 21 / 21

Min-max Time Consensus Tracking Over Directed Trees

Min-max Time Consensus Tracking Over Directed Trees Min-max Time Consensus Tracking Over Directed Trees Ameer K. Mulla, Sujay Bhatt H. R., Deepak U. Patil, and Debraj Chakraborty Abstract In this paper, decentralized feedback control strategies are derived

More information

A Guaranteed Cost LMI-Based Approach for Event-Triggered Average Consensus in Multi-Agent Networks

A Guaranteed Cost LMI-Based Approach for Event-Triggered Average Consensus in Multi-Agent Networks A Guaranteed Cost LMI-Based Approach for Event-Triggered Average Consensus in Multi-Agent Networks Amir Amini, Arash Mohammadi, Amir Asif Electrical and Computer Engineering,, Montreal, Canada. Concordia

More information

A Robust Event-Triggered Consensus Strategy for Linear Multi-Agent Systems with Uncertain Network Topology

A Robust Event-Triggered Consensus Strategy for Linear Multi-Agent Systems with Uncertain Network Topology A Robust Event-Triggered Consensus Strategy for Linear Multi-Agent Systems with Uncertain Network Topology Amir Amini, Amir Asif, Arash Mohammadi Electrical and Computer Engineering,, Montreal, Canada.

More information

9.5. Polynomial and Rational Inequalities. Objectives. Solve quadratic inequalities. Solve polynomial inequalities of degree 3 or greater.

9.5. Polynomial and Rational Inequalities. Objectives. Solve quadratic inequalities. Solve polynomial inequalities of degree 3 or greater. Chapter 9 Section 5 9.5 Polynomial and Rational Inequalities Objectives 1 3 Solve quadratic inequalities. Solve polynomial inequalities of degree 3 or greater. Solve rational inequalities. Objective 1

More information

Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control

Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control Outline Background Preliminaries Consensus Numerical simulations Conclusions Average-Consensus of Multi-Agent Systems with Direct Topology Based on Event-Triggered Control Email: lzhx@nankai.edu.cn, chenzq@nankai.edu.cn

More information

When Gradient Systems and Hamiltonian Systems Meet

When Gradient Systems and Hamiltonian Systems Meet When Gradient Systems and Hamiltonian Systems Meet Arjan van der Schaft Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, the Netherlands December 11, 2011 on the

More information

Cutting Planes for First Level RLT Relaxations of Mixed 0-1 Programs

Cutting Planes for First Level RLT Relaxations of Mixed 0-1 Programs Cutting Planes for First Level RLT Relaxations of Mixed 0-1 Programs 1 Cambridge, July 2013 1 Joint work with Franklin Djeumou Fomeni and Adam N. Letchford Outline 1. Introduction 2. Literature Review

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 COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES

ON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES ON COST MATRICES WITH TWO AND THREE DISTINCT VALUES OF HAMILTONIAN PATHS AND CYCLES SANTOSH N. KABADI AND ABRAHAM P. PUNNEN Abstract. Polynomially testable characterization of cost matrices associated

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

CSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization

CSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization CSCI 1951-G Optimization Methods in Finance Part 10: Conic Optimization April 6, 2018 1 / 34 This material is covered in the textbook, Chapters 9 and 10. Some of the materials are taken from it. Some of

More information

Min-Max Message Passing and Local Consistency in Constraint Networks

Min-Max Message Passing and Local Consistency in Constraint Networks Min-Max Message Passing and Local Consistency in Constraint Networks Hong Xu, T. K. Satish Kumar, and Sven Koenig University of Southern California, Los Angeles, CA 90089, USA hongx@usc.edu tkskwork@gmail.com

More information

OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL

OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL OUTPUT CONSENSUS OF HETEROGENEOUS LINEAR MULTI-AGENT SYSTEMS BY EVENT-TRIGGERED CONTROL Gang FENG Department of Mechanical and Biomedical Engineering City University of Hong Kong July 25, 2014 Department

More information

Consensus Protocols for Networks of Dynamic Agents

Consensus Protocols for Networks of Dynamic Agents Consensus Protocols for Networks of Dynamic Agents Reza Olfati Saber Richard M. Murray Control and Dynamical Systems California Institute of Technology Pasadena, CA 91125 e-mail: {olfati,murray}@cds.caltech.edu

More information

Optimal Network Topology Design in Multi-Agent Systems for Efficient Average Consensus

Optimal Network Topology Design in Multi-Agent Systems for Efficient Average Consensus 49th IEEE Conference on Decision and Control December 15-17, 010 Hilton Atlanta Hotel, Atlanta, GA, USA Optimal Network Topology Design in Multi-Agent Systems for Efficient Average Consensus Mohammad Rafiee

More information

Disconnecting Networks via Node Deletions

Disconnecting Networks via Node Deletions 1 / 27 Disconnecting Networks via Node Deletions Exact Interdiction Models and Algorithms Siqian Shen 1 J. Cole Smith 2 R. Goli 2 1 IOE, University of Michigan 2 ISE, University of Florida 2012 INFORMS

More information

10.7 Polynomial and Rational Inequalities

10.7 Polynomial and Rational Inequalities 10.7 Polynomial and Rational Inequalities In this section we want to turn our attention to solving polynomial and rational inequalities. That is, we want to solve inequalities like 5 4 0. In order to do

More information

- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs

- Well-characterized problems, min-max relations, approximate certificates. - LP problems in the standard form, primal and dual linear programs LP-Duality ( Approximation Algorithms by V. Vazirani, Chapter 12) - Well-characterized problems, min-max relations, approximate certificates - LP problems in the standard form, primal and dual linear programs

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

x y More precisely, this equation means that given any ε > 0, there exists some δ > 0 such that

x y More precisely, this equation means that given any ε > 0, there exists some δ > 0 such that Chapter 2 Limits and continuity 21 The definition of a it Definition 21 (ε-δ definition) Let f be a function and y R a fixed number Take x to be a point which approaches y without being equal to y If there

More information

MA102: Multivariable Calculus

MA102: Multivariable Calculus MA102: Multivariable Calculus Rupam Barman and Shreemayee Bora Department of Mathematics IIT Guwahati Differentiability of f : U R n R m Definition: Let U R n be open. Then f : U R n R m is differentiable

More information

Differential Equations on Cubic Julia Sets

Differential Equations on Cubic Julia Sets University of California Los Angeles June 16, 2014 Julia Sets Given a complex polynomial f, the associated filled-in Julia set is the set of points z such that the sequence z, f (z), f (f (z)),..., f

More information

CSE 254: MAP estimation via agreement on (hyper)trees: Message-passing and linear programming approaches

CSE 254: MAP estimation via agreement on (hyper)trees: Message-passing and linear programming approaches CSE 254: MAP estimation via agreement on (hyper)trees: Message-passing and linear programming approaches A presentation by Evan Ettinger November 11, 2005 Outline Introduction Motivation and Background

More information

The Transition Probability Function P ij (t)

The Transition Probability Function P ij (t) The Transition Probability Function P ij (t) Consider a continuous time Markov chain {X(t), t 0}. We are interested in the probability that in t time units the process will be in state j, given that it

More information

BCOL RESEARCH REPORT 17.07

BCOL RESEARCH REPORT 17.07 BCOL RESEARCH REPORT 17.07 Industrial Engineering & Operations Research University of California, Berkeley, CA 9470 1777 ON CAPACITY MODELS FOR NETWORK DESIGN ALPER ATAMTÜRK AND OKTAY GÜNLÜK Abstract.

More information

Definition: A "system" of equations is a set or collection of equations that you deal with all together at once.

Definition: A system of equations is a set or collection of equations that you deal with all together at once. System of Equations Definition: A "system" of equations is a set or collection of equations that you deal with all together at once. There is both an x and y value that needs to be solved for Systems

More information

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation:

Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: Chapter 7 Heat Equation Partial differential equation for temperature u(x, t) in a heat conducting insulated rod along the x-axis is given by the Heat equation: u t = ku x x, x, t > (7.1) Here k is a constant

More information

Structure-Based Comparison of Biomolecules

Structure-Based Comparison of Biomolecules Structure-Based Comparison of Biomolecules Benedikt Christoph Wolters Seminar Bioinformatics Algorithms RWTH AACHEN 07/17/2015 Outline 1 Introduction and Motivation Protein Structure Hierarchy Protein

More information

Four-coloring P 6 -free graphs. I. Extending an excellent precoloring

Four-coloring P 6 -free graphs. I. Extending an excellent precoloring Four-coloring P 6 -free graphs. I. Extending an excellent precoloring Maria Chudnovsky Princeton University, Princeton, NJ 08544 Sophie Spirkl Princeton University, Princeton, NJ 08544 Mingxian Zhong Columbia

More information

Chapter 2 Optimal Control Problem

Chapter 2 Optimal Control Problem Chapter 2 Optimal Control Problem Optimal control of any process can be achieved either in open or closed loop. In the following two chapters we concentrate mainly on the first class. The first chapter

More information

How Many Leaders Are Needed for Reaching a Consensus

How Many Leaders Are Needed for Reaching a Consensus How Many Leaders Are Needed for Reaching a Consensus Xiaoming Hu KTH, Sweden Joint work with Zhixin Liu, Jing Han Academy of Mathematics and Systems Science, Chinese Academy of Sciences Stockholm, September

More information

4.4 Noetherian Rings

4.4 Noetherian Rings 4.4 Noetherian Rings Recall that a ring A is Noetherian if it satisfies the following three equivalent conditions: (1) Every nonempty set of ideals of A has a maximal element (the maximal condition); (2)

More information

Structural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems

Structural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems Structural Consensus Controllability of Singular Multi-agent Linear Dynamic Systems M. ISAL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, sc. C, 1-3, 08038 arcelona

More information

Geometric control and dynamical systems

Geometric control and dynamical systems Université de Nice - Sophia Antipolis & Institut Universitaire de France 9th AIMS International Conference on Dynamical Systems, Differential Equations and Applications Control of an inverted pendulum

More information

Blackwell s Approachability Theorem: A Generalization in a Special Case. Amy Greenwald, Amir Jafari and Casey Marks

Blackwell s Approachability Theorem: A Generalization in a Special Case. Amy Greenwald, Amir Jafari and Casey Marks Blackwell s Approachability Theorem: A Generalization in a Special Case Amy Greenwald, Amir Jafari and Casey Marks Department of Computer Science Brown University Providence, Rhode Island 02912 CS-06-01

More information

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

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

More information

The Intersection Problem for Steiner Triple Systems. Bradley Fain

The Intersection Problem for Steiner Triple Systems. Bradley Fain The Intersection Problem for Steiner Triple Systems by Bradley Fain A thesis submitted to the Graduate Faculty of Auburn University in partial fulfillment of the requirements for the Degree of Master of

More information

Powers and Products of Monomial Ideals

Powers and Products of Monomial Ideals U.U.D.M. Project Report 6: Powers and Products of Monomial Ideals Melker Epstein Examensarbete i matematik, 5 hp Handledare: Veronica Crispin Quinonez Examinator: Jörgen Östensson Juni 6 Department of

More information

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection

Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection Stationary trajectories, singular Hamiltonian systems and ill-posed Interconnection S.C. Jugade, Debasattam Pal, Rachel K. Kalaimani and Madhu N. Belur Department of Electrical Engineering Indian Institute

More information

Duality of LPs and Applications

Duality of LPs and Applications Lecture 6 Duality of LPs and Applications Last lecture we introduced duality of linear programs. We saw how to form duals, and proved both the weak and strong duality theorems. In this lecture we will

More information

Consensus, Flocking and Opinion Dynamics

Consensus, Flocking and Opinion Dynamics Consensus, Flocking and Opinion Dynamics Antoine Girard Laboratoire Jean Kuntzmann, Université de Grenoble antoine.girard@imag.fr International Summer School of Automatic Control GIPSA Lab, Grenoble, France,

More information

Lecture Introduction. 2 Brief Recap of Lecture 10. CS-621 Theory Gems October 24, 2012

Lecture Introduction. 2 Brief Recap of Lecture 10. CS-621 Theory Gems October 24, 2012 CS-62 Theory Gems October 24, 202 Lecture Lecturer: Aleksander Mądry Scribes: Carsten Moldenhauer and Robin Scheibler Introduction In Lecture 0, we introduced a fundamental object of spectral graph theory:

More information

NP-Complete Problems and Approximation Algorithms

NP-Complete Problems and Approximation Algorithms NP-Complete Problems and Approximation Algorithms Efficiency of Algorithms Algorithms that have time efficiency of O(n k ), that is polynomial of the input size, are considered to be tractable or easy

More information

Control Systems. Frequency domain analysis. L. Lanari

Control Systems. Frequency domain analysis. L. Lanari Control Systems m i l e r p r a in r e v y n is o Frequency domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic

More information

Newton s method in eigenvalue optimization for incomplete pairwise comparison matrices

Newton s method in eigenvalue optimization for incomplete pairwise comparison matrices Newton s method in eigenvalue optimization for incomplete pairwise comparison matrices Kristóf Ábele-Nagy Eötvös Loránd University (ELTE); Corvinus University of Budapest (BCE) Sándor Bozóki Computer and

More information

Existence and uniqueness of solutions for a continuous-time opinion dynamics model with state-dependent connectivity

Existence and uniqueness of solutions for a continuous-time opinion dynamics model with state-dependent connectivity Existence and uniqueness of solutions for a continuous-time opinion dynamics model with state-dependent connectivity Vincent D. Blondel, Julien M. Hendricx and John N. Tsitsilis July 24, 2009 Abstract

More information

To get horizontal and slant asymptotes algebraically we need to know about end behaviour for rational functions.

To get horizontal and slant asymptotes algebraically we need to know about end behaviour for rational functions. Concepts: Horizontal Asymptotes, Vertical Asymptotes, Slant (Oblique) Asymptotes, Transforming Reciprocal Function, Sketching Rational Functions, Solving Inequalities using Sign Charts. Rational Function

More information

Clairvoyant scheduling of random walks

Clairvoyant scheduling of random walks Clairvoyant scheduling of random walks Péter Gács Boston University April 25, 2008 Péter Gács (BU) Clairvoyant demon April 25, 2008 1 / 65 Introduction The clairvoyant demon problem 2 1 Y : WAIT 0 X, Y

More information

Minimum Fuel Optimal Control Example For A Scalar System

Minimum Fuel Optimal Control Example For A Scalar System Minimum Fuel Optimal Control Example For A Scalar System A. Problem Statement This example illustrates the minimum fuel optimal control problem for a particular first-order (scalar) system. The derivation

More information

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES

AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES AN ELEMENTARY PROOF OF THE SPECTRAL RADIUS FORMULA FOR MATRICES JOEL A. TROPP Abstract. We present an elementary proof that the spectral radius of a matrix A may be obtained using the formula ρ(a) lim

More information

NOTES ON CALCULUS OF VARIATIONS. September 13, 2012

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

More information

Computational Aspects of Aggregation in Biological Systems

Computational Aspects of Aggregation in Biological Systems Computational Aspects of Aggregation in Biological Systems Vladik Kreinovich and Max Shpak University of Texas at El Paso, El Paso, TX 79968, USA vladik@utep.edu, mshpak@utep.edu Summary. Many biologically

More information

Multi-colouring the Mycielskian of Graphs

Multi-colouring the Mycielskian of Graphs Multi-colouring the Mycielskian of Graphs Wensong Lin Daphne Der-Fen Liu Xuding Zhu December 1, 016 Abstract A k-fold colouring of a graph is a function that assigns to each vertex a set of k colours,

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

Reverse engineering using computational algebra

Reverse engineering using computational algebra Reverse engineering using computational algebra Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4500, Spring 2015 M. Macauley (Clemson)

More information

Lecture 7 Monotonicity. September 21, 2008

Lecture 7 Monotonicity. September 21, 2008 Lecture 7 Monotonicity September 21, 2008 Outline Introduce several monotonicity properties of vector functions Are satisfied immediately by gradient maps of convex functions In a sense, role of monotonicity

More information

Lecture Note 7: Switching Stabilization via Control-Lyapunov Function

Lecture Note 7: Switching Stabilization via Control-Lyapunov Function ECE7850: Hybrid Systems:Theory and Applications Lecture Note 7: Switching Stabilization via Control-Lyapunov Function Wei Zhang Assistant Professor Department of Electrical and Computer Engineering Ohio

More information

1 Convexity explains SVMs

1 Convexity explains SVMs 1 Convexity explains SVMs The convex hull of a set is the collection of linear combinations of points in the set where the coefficients are nonnegative and sum to one. Two sets are linearly separable if

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechanics Rajdeep Sensarma! sensarma@theory.tifr.res.in Lecture #22 Path Integrals and QM Recap of Last Class Statistical Mechanics and path integrals in imaginary time Imaginary time

More information

No EFFICIENT LINE SEARCHING FOR CONVEX FUNCTIONS. By E. den Boef, D. den Hertog. May 2004 ISSN

No EFFICIENT LINE SEARCHING FOR CONVEX FUNCTIONS. By E. den Boef, D. den Hertog. May 2004 ISSN No. 4 5 EFFICIENT LINE SEARCHING FOR CONVEX FUNCTIONS y E. den oef, D. den Hertog May 4 ISSN 94-785 Efficient Line Searching for Convex Functions Edgar den oef Dick den Hertog 3 Philips Research Laboratories,

More information

Solutions for Math 217 Assignment #3

Solutions for Math 217 Assignment #3 Solutions for Math 217 Assignment #3 (1) Which of the following sets in R n are open? Which are closed? Which are neither open nor closed? (a) {(x, y) R 2 : x 2 y 2 = 1}. (b) {(x, y, z) R 3 : 0 < x + y

More information

Opinion Formation under Bounded Confidence via Gossip Algorithms and Applications

Opinion Formation under Bounded Confidence via Gossip Algorithms and Applications Opinion Formation under Bounded Confidence via Gossip Algorithms and Applications Nguyen Thi Hoai Linh 1, Takayuki Wada 1, Izumi Masubuchi 2, Toru Asai 3, Yasumasa Fujisaki 1 1 Osaka University, 2 Kobe

More information

Convex optimization examples

Convex optimization examples Convex optimization examples multi-period processor speed scheduling minimum time optimal control grasp force optimization optimal broadcast transmitter power allocation phased-array antenna beamforming

More information

Intersection of Ellipsoids

Intersection of Ellipsoids Intersection of Ellipsoids David Eberly, Geometric Tools, Redmond WA 9805 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International License. To view

More information

A Graph-Theoretic Characterization of Structural Controllability for Multi-Agent System with Switching Topology

A Graph-Theoretic Characterization of Structural Controllability for Multi-Agent System with Switching Topology Joint 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 29 FrAIn2.3 A Graph-Theoretic Characterization of Structural Controllability

More information

NCS Lecture 8 A Primer on Graph Theory. Cooperative Control Applications

NCS Lecture 8 A Primer on Graph Theory. Cooperative Control Applications NCS Lecture 8 A Primer on Graph Theory Richard M. Murray Control and Dynamical Systems California Institute of Technology Goals Introduce some motivating cooperative control problems Describe basic concepts

More information

Laplace Transform Theory - 1

Laplace Transform Theory - 1 Laplace Transform Theory - 1 Existence of Laplace Transforms Before continuing our use of Laplace transforms for solving DEs, it is worth digressing through a quick investigation of which functions actually

More information

Suppose that the approximate solutions of Eq. (1) satisfy the condition (3). Then (1) if η = 0 in the algorithm Trust Region, then lim inf.

Suppose that the approximate solutions of Eq. (1) satisfy the condition (3). Then (1) if η = 0 in the algorithm Trust Region, then lim inf. Maria Cameron 1. Trust Region Methods At every iteration the trust region methods generate a model m k (p), choose a trust region, and solve the constraint optimization problem of finding the minimum of

More information

A A x i x j i j (i, j) (j, i) Let. Compute the value of for and

A A x i x j i j (i, j) (j, i) Let. Compute the value of for and 7.2 - Quadratic Forms quadratic form on is a function defined on whose value at a vector in can be computed by an expression of the form, where is an symmetric matrix. The matrix R n Q R n x R n Q(x) =

More information

Distributed Receding Horizon Control of Cost Coupled Systems

Distributed Receding Horizon Control of Cost Coupled Systems Distributed Receding Horizon Control of Cost Coupled Systems William B. Dunbar Abstract This paper considers the problem of distributed control of dynamically decoupled systems that are subject to decoupled

More information

arxiv: v2 [math.oc] 14 Dec 2015

arxiv: v2 [math.oc] 14 Dec 2015 Cooperative Output Regulation of Discrete-Time Linear Time-Delay Multi-agent Systems Yamin Yan and Jie Huang arxiv:1510.05380v2 math.oc 14 Dec 2015 Abstract In this paper, we study the cooperative output

More information

1 Review: symmetric matrices, their eigenvalues and eigenvectors

1 Review: symmetric matrices, their eigenvalues and eigenvectors Cornell University, Fall 2012 Lecture notes on spectral methods in algorithm design CS 6820: Algorithms Studying the eigenvalues and eigenvectors of matrices has powerful consequences for at least three

More information

A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model

A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP

More information

Stein s Method for concentration inequalities

Stein s Method for concentration inequalities Number of Triangles in Erdős-Rényi random graph, UC Berkeley joint work with Sourav Chatterjee, UC Berkeley Cornell Probability Summer School July 7, 2009 Number of triangles in Erdős-Rényi random graph

More information

Reachability and Controllability

Reachability and Controllability Capitolo. INTRODUCTION 4. Reachability and Controllability Reachability. The reachability problem is to find the set of all the final states x(t ) reachable starting from a given initial state x(t ) :

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

Distributed motion constraints for algebraic connectivity of robotic networks

Distributed motion constraints for algebraic connectivity of robotic networks Distributed motion constraints for algebraic connectivity of robotic networks IEEE Conference on Decision and Control December, 2008 Michael Schuresko Jorge Cortés Applied Mathematics and Statistics University

More information

Lecture 9 Monotone VIs/CPs Properties of cones and some existence results. October 6, 2008

Lecture 9 Monotone VIs/CPs Properties of cones and some existence results. October 6, 2008 Lecture 9 Monotone VIs/CPs Properties of cones and some existence results October 6, 2008 Outline Properties of cones Existence results for monotone CPs/VIs Polyhedrality of solution sets Game theory:

More information

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi Real Analysis Math 3AH Rudin, Chapter # Dominique Abdi.. If r is rational (r 0) and x is irrational, prove that r + x and rx are irrational. Solution. Assume the contrary, that r+x and rx are rational.

More information

Sample Exam 1 KEY NAME: 1. CS 557 Sample Exam 1 KEY. These are some sample problems taken from exams in previous years. roughly ten questions.

Sample Exam 1 KEY NAME: 1. CS 557 Sample Exam 1 KEY. These are some sample problems taken from exams in previous years. roughly ten questions. Sample Exam 1 KEY NAME: 1 CS 557 Sample Exam 1 KEY These are some sample problems taken from exams in previous years. roughly ten questions. Your exam will have 1. (0 points) Circle T or T T Any curve

More information

ELEMENTARY LINEAR ALGEBRA

ELEMENTARY LINEAR ALGEBRA ELEMENTARY LINEAR ALGEBRA DAVID MEREDITH These notes are about the application and computation of linear algebra objects in F n, where F is a field. I think everything here works over any field including

More information

Dr. Tess J. Moon, (office), (fax),

Dr. Tess J. Moon, (office), (fax), THE UNIVERSITY OF TEXAS AT AUSTIN DEPARTMENT OF MECHANICAL ENGINEERING ME 380Q ENGINEERING ANALYSIS: ANALYTICAL METHODS Dr Tess J Moon 471-0094 (office) 471-877 (fax) tmoon@mailutexasedu Mathematical Preliminaries

More information

Markov Decision Processes

Markov Decision Processes Markov Decision Processes Lecture notes for the course Games on Graphs B. Srivathsan Chennai Mathematical Institute, India 1 Markov Chains We will define Markov chains in a manner that will be useful to

More information

Conflict-Free Colorings of Rectangles Ranges

Conflict-Free Colorings of Rectangles Ranges Conflict-Free Colorings of Rectangles Ranges Khaled Elbassioni Nabil H. Mustafa Max-Planck-Institut für Informatik, Saarbrücken, Germany felbassio, nmustafag@mpi-sb.mpg.de Abstract. Given the range space

More information

arxiv: v1 [cs.sy] 2 Apr 2019

arxiv: v1 [cs.sy] 2 Apr 2019 On the Existence of a Fixed Spectrum for a Multi-channel Linear System: A Matroid Theory Approach F Liu 1 and A S Morse 1 arxiv:190401499v1 [cssy] 2 Apr 2019 Abstract Conditions for the existence of a

More information

On Controllability of Linear Systems 1

On Controllability of Linear Systems 1 On Controllability of Linear Systems 1 M.T.Nair Department of Mathematics, IIT Madras Abstract In this article we discuss some issues related to the observability and controllability of linear systems.

More information

Architectures and Algorithms for Distributed Generation Control of Inertia-Less AC Microgrids

Architectures and Algorithms for Distributed Generation Control of Inertia-Less AC Microgrids Architectures and Algorithms for Distributed Generation Control of Inertia-Less AC Microgrids Alejandro D. Domínguez-García Coordinated Science Laboratory Department of Electrical and Computer Engineering

More information

Obtaining Consensus of Multi-agent Linear Dynamic Systems

Obtaining Consensus of Multi-agent Linear Dynamic Systems Obtaining Consensus of Multi-agent Linear Dynamic Systems M I GRCÍ-PLNS Universitat Politècnica de Catalunya Departament de Matemàtica plicada Mineria 1, 08038 arcelona SPIN mariaisabelgarcia@upcedu bstract:

More information

Distributed Game Strategy Design with Application to Multi-Agent Formation Control

Distributed Game Strategy Design with Application to Multi-Agent Formation Control 5rd IEEE Conference on Decision and Control December 5-7, 4. Los Angeles, California, USA Distributed Game Strategy Design with Application to Multi-Agent Formation Control Wei Lin, Member, IEEE, Zhihua

More information

L 2 -induced Gains of Switched Systems and Classes of Switching Signals

L 2 -induced Gains of Switched Systems and Classes of Switching Signals L 2 -induced Gains of Switched Systems and Classes of Switching Signals Kenji Hirata and João P. Hespanha Abstract This paper addresses the L 2-induced gain analysis for switched linear systems. We exploit

More information

Dual Decomposition.

Dual Decomposition. 1/34 Dual Decomposition http://bicmr.pku.edu.cn/~wenzw/opt-2017-fall.html Acknowledgement: this slides is based on Prof. Lieven Vandenberghes lecture notes Outline 2/34 1 Conjugate function 2 introduction:

More information

1 The Erdős Ko Rado Theorem

1 The Erdős Ko Rado Theorem 1 The Erdős Ko Rado Theorem A family of subsets of a set is intersecting if any two elements of the family have at least one element in common It is easy to find small intersecting families; the basic

More information

Autonomous navigation of unicycle robots using MPC

Autonomous navigation of unicycle robots using MPC Autonomous navigation of unicycle robots using MPC M. Farina marcello.farina@polimi.it Dipartimento di Elettronica e Informazione Politecnico di Milano 7 June 26 Outline Model and feedback linearization

More information

CURVATURE ESTIMATES FOR WEINGARTEN HYPERSURFACES IN RIEMANNIAN MANIFOLDS

CURVATURE ESTIMATES FOR WEINGARTEN HYPERSURFACES IN RIEMANNIAN MANIFOLDS CURVATURE ESTIMATES FOR WEINGARTEN HYPERSURFACES IN RIEMANNIAN MANIFOLDS CLAUS GERHARDT Abstract. We prove curvature estimates for general curvature functions. As an application we show the existence of

More information

On Backward Product of Stochastic Matrices

On Backward Product of Stochastic Matrices On Backward Product of Stochastic Matrices Behrouz Touri and Angelia Nedić 1 Abstract We study the ergodicity of backward product of stochastic and doubly stochastic matrices by introducing the concept

More information

Partial characterizations of clique-perfect graphs II: diamond-free and Helly circular-arc graphs

Partial characterizations of clique-perfect graphs II: diamond-free and Helly circular-arc graphs Partial characterizations of clique-perfect graphs II: diamond-free and Helly circular-arc graphs Flavia Bonomo a,1, Maria Chudnovsky b,2 and Guillermo Durán c,3 a Departamento de Matemática, Facultad

More information

Nowhere zero flow. Definition: A flow on a graph G = (V, E) is a pair (D, f) such that. 1. D is an orientation of G. 2. f is a function on E.

Nowhere zero flow. Definition: A flow on a graph G = (V, E) is a pair (D, f) such that. 1. D is an orientation of G. 2. f is a function on E. Nowhere zero flow Definition: A flow on a graph G = (V, E) is a pair (D, f) such that 1. D is an orientation of G. 2. f is a function on E. 3. u N + D (v) f(uv) = w ND f(vw) for every (v) v V. Example:

More information

LECTURE 4 LECTURE OUTLINE

LECTURE 4 LECTURE OUTLINE LECTURE 4 LECTURE OUTLINE Relative interior and closure Algebra of relative interiors and closures Continuity of convex functions Closures of functions Reading: Section 1.3 All figures are courtesy of

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

CATEGORIZATION GENERATED BY PROTOTYPES AN AXIOMATIC APPROACH

CATEGORIZATION GENERATED BY PROTOTYPES AN AXIOMATIC APPROACH CATEGORIZATION GENERATED BY PROTOTYPES AN AXIOMATIC APPROACH Yaron Azrieli and Ehud Lehrer Tel-Aviv University CATEGORIZATION GENERATED BY PROTOTYPES AN AXIOMATIC APPROACH p. 1/3 Outline Categorization

More information