Three Lectures on: Control of Coupled Fast and Slow Dynamics. Zvi Artstein Ravello, September 2012

Size: px
Start display at page:

Download "Three Lectures on: Control of Coupled Fast and Slow Dynamics. Zvi Artstein Ravello, September 2012"

Transcription

1 Three Lectures on: Control of Coupled Fast and Slow Dynamics Zvi Artstein Ravello, September

2 Control of Coupled Fast and Slow Dynamics Zvi Artstein Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates A Future Direction 2

3 3 Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates A Future Direction

4 4 Example from real life: Airplane

5 5 Example from real life: Helicopter

6 6 Example from real life: The hummingbird

7 The framework ODE: An ordinary differential equation 7

8 A control equation The framework CONTROL: Bolza Problem Mayer Problem Possible goals: 8

9 A reduction of Bolza to Mayer: The goal can be reduced to Adolf Mayer by adding a coordinate and an equation and seeking minimizing the additional coordinate Oskar Bolza

10 The equivalent differential inclusion: A control equation can be written as: 10

11 11 How to model coupled slow and fast motions?

12 Tikhonov s Singular Perturbations model of coupled slow and fast motions The perturbed system: The fast part can be written as: 12 Where: in the slow and in the fast, variables We are interested in the limit behavior of the system as

13 13 A mathematical example capturing reality: An elastic structure in a rapidly flowing nearly invicid fluid (with Marshall Slemrod)

14 To make the long story short: Based on a model of Iwan/Belvins and Dowel/Ilgamov, the limit (after normalization) equations: With a generator of a van der Pol oscillator 14

15 15 An example Relaxation oscillation

16 Singularly perturbed optimal control systems: Where: in the slow and in the fast, variables Of interest: The behavior of the system as 16

17 Applications A variety of natural phenomena and engineering design. The latter include: Regulation, LQ-Systems, Feedback Design, Stabilization, Robustness, Stochastics, Filters, Optimal Control, Hydropower Production, Nuclear Reactions, Aircraft Design, Flight Control, and many more! The classical approach to handle the applications is the model-reduction - after Levinson-Tikhonov in the differential equations trait and after Kokotovic in the control Setting 17

18 18 Other coupled slow and fast dynamics

19 19 Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates A Future Direction

20 The definition of a variational limit: Given a system (an equation or a control equation) with a parameter that tends to a limit. A variational limit is a system whose solutions capture the limit behavior of the solutions of the parameterized system, as the parameter tends to its limit, Capture = limit of the trajectories, limit of the optimal controls, limit of the values 20

21 21 Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates A Future Direction

22 The classical Tikhonov order reduction approach Write the perturbed system as: The limit behavior as is captured by the system: 22

23 Andrei Nikolayevich Tikhonov Norman Levinson

24 24 The geometry of the solution:

25 25 An example Relaxation oscillation

26 Recall: Singularly perturbed control systems: Of interest: The behavior of the system as What is the variational limit? 26

27 The limit of: is 27

28 28 Petar Kokotovic

29 29 The geometry of the solution:

30 30 BUT

31 The general situation: There is no reason why the optimal fast solution will converge and not, say, oscillate! 31

32 This was pointed out in the mid 1980 s, independently, By Assen Dontchev and Valadimir Veliov And by Vladimir Gaitsgory 32

33 33 A mathematical illustration: Non-stationary relaxation oscillation

34 34 Recall: A mathematical example capturing reality: An elastic structure in a rapidly flowing nearly invicid fluid

35 Numerical results: The slow dynamics The fast dynamics Computations by: 35 Zvi Artstein Jasmine Linshiz Edriss Titi:

36 36 Also in Nature: The hummingbird

37 An illustration of a control problem: 37 The questions: when should the switch be made? How should this be carried out when the speed is very fast?

38 An example after V. Veliov 1996 Applying an order reduction (i.e. plugging ) yields zero value. Clearly one can do better! 38

39 39 No equilibrium in the limit:

40 The limit solution: The limit strategy as bang-bang feedback can be expressed as a resulting in: Limit trajectories The bang-bang feedback 40

41 41 Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates A Future Direction

42 Strong limit on a space of functions: Recall the strong limit in, say L 2 : The strong limit may not work for the variational limit of singular perturbations recall the examples. 42

43 Weak limit on a space of functions: The L 2 weak-limit: The sequence converges weakly to if The weak limit may not work for the variational limit of singular perturbations recall the examples. 43

44 On a space of parameters: Consider an ordinary differential equation with parameters Strong convergence of the parameters implies continuous dependence of solution but is not compact Weak convergence of the parameters is compact but does not imply continuous dependence 44

45 Observation: Continuous dependence and compactness are opposing properties! Can we construct a convergence that will have both properties: continuous dependence of solution and compactness? 45

46 46 To the rescue

47 Laurence Chisholm Young July 14, December 24, 2000 Cambridge, England, Madison WI, USA 47

48 The price: The limit function will be out of the original space It is called: A Young measure 48

49 An implicit Definition of a Young Measure: Let be a sequence of parameter functions (say bounded from an interval to. There exist a subsequence (say the sequence itself) and a family of probability measures on parameterized by such that for every right hand side converges weakly to 49

50 Proof Based on (simple) functional analysis arguments (incorporating weak* convergence and Alaoglu compactness Theorem). 50

51 A consequence: Solutions of the ordinary differential equation with parameters converge to solution of the ordinary differential equation with the Young measure 51

52 A constructive Definition of a Young Measure: Let be a metric space Denote by the family of probability measures on Let be another metric space endowed with a measure (say Lebesgue measure on an interval) Definition: A mapping from to is a Young Measure 52

53 53 A Pictorial Definition:

54 The structure of the space : The elements: -additive set-functions from the Borel subsets onto the unit interval. Convergence of if For every continuous and bounded. Prohorov metric between measures and is the smallest such that for every Borel set 54

55 Consequences concerning : If If is complete and separable so is is compact so is 55

56 The structure of Young measures : Can be viewed as a probability measure on Consequences: If is complete and separable so is the space of Young measures If is compact so is the space of Young measures 56

57 A major property of Young measures: An ordinary function can be viewed as a Diracvalued Young measure. 57

58 58 The nature of the convergence

59 For instance: The sequence converges to a Young measure with a constant value, namely the measure on given by 59

60 Functions are dense in the space of Young Measures! When the underlying space is without atoms then any Young measure can be approximated by a function 60

61 When the limit is a function: The space of Young Measures completes the space of functions. What convergence does it reflect if the limit Young measure happens to be a function? 61 The limit is then strong (,, but not )

62 Key properties: Existence of the limit Keeping information about the location of the values Possibility to approximate by an ordinary function 62

63 Recall the case of a space of parameters: Consider an ordinary differential equation with parameters Strong convergence implies continuous dependence of solution but is not compact Weak convergence is compact but does not imply continuous dependence What happens if measure? converges to a Young 63

64 Definition: If and is a probability measure then Likewise, for an ordinary differential equation we mean 64

65 Main application: Consider an ordinary differential equation with parameters And converges to a Young measure Then the solutions of the odes with parameters converge to the solution of the ode with the Young measure Thus, the convergence to the Young measure is both compact and implies continuous dependence 65

66 The key tool: Consider the right hand side of the ordinary differential equation with parameters and Then converges to a Young measure converges weakly to 66

67 67 Plan: Modeling Variational Limits Classical Approach to slow-fast dynamics What limits are appropriate? Young Measures Modern Approach to slow-fast dynamics Other chattering limits and averaging techniques Control Invariant Measures Stabilization Optimal Control Some special cases Computations, error estimates Future Directions

68 68 Applications to Control and the Calculus of Variations

69 69 An application (after L.C. Young):

70 70 A problem without a solution:

71 71 Approximate solutions:

72 Better approximations: 72 No ordinary limit which is useful!

73 73 A generalized limit of a Young measure type: called Relaxed Control

74 The effect of Relaxed Control: Convexification of the vector field Is there a control that makes a solution? Yes, the control that averages +1 and

75 Jack Warga

76 A prior appearance in the calculus of variations: The general problem: The particular problem without a solution: 76

77 Generalized curves in the calculus of variations: For the problem: A generalized curve is a pair: satisfying: The goal: 77

78 Recall the illustration of a control problem: The questions: when should the switch be made? How should this be carried out when the speed is very fast? 78

79 79 In the limit: A distribution arises

80 The limit solution: The limit strategy as bang-bang feedback can be expressed as a resulting in: Limit measure The bang-bang feedback 80

81 81 Recall: A mathematical example capturing reality: An elastic structure in a rapidly flowing nearly invicid fluid

82 Numerical results: The slow dynamics The fast dynamics 82

83 Recall: Singular perturbations as a model of coupled slow and fast motions: The perturbed system: Equivalently: We are interested in the limit behavior of the system as 83

84 The classical Tikhonov approach Write the perturbed system as: The limit behavior as is captured by the system: This type of variational limit does not capture the general situation 84

85 A program to exploit Young Measures this started by Zvi Artstein and Alexander Vigodner,

86 The general situation: The solution oscillates. The limit (as ) may be described as a Young measure 86

87 The general situation: The Young measure is defined on the x-space with values being probability measures on the y-space 87

88 The variational limit solution in the new formulation : where and is a Young measure solves the averaging equation 88

89 89 Now to control systems

90 Recall: Singularly perturbed control systems: Where: in the slow and in the fast, variables Of interest: The behavior of the system as 90

91 The order reduction method (Petar Kokotovic et al.) The limit as is depicted by namely, by: 91

92 The general situation: There is no reason why the optimal fast solution will converge and not, say, oscillate! 92

93 The general situation: The values of the Young measure are: measures of the (fast state, control) dynamics! 93

94 The general variational limit solution is of the form: Where: solves the averaging equation is a Young measure (parameterized by ) and the limit cost is based on averaging: Notice, the values of the Young measure are the control variables, (replacing the equilibrium points in the classical case 94

95 The equivalent differential inclusion: solves the differential inclusion where is a Young measure (parameterized by ) Notice, the values of the Young measure are the control variables, here they determine the velocity of the slow variable 95

96 A question: Could any probability measure be a value for the Young Measure of the variational limit? If not, how can the possible values be classified and identified? A promise: We shall soon give a characterization of the probability measures that may appear as values in the variational limit. We denote this family by 96

97 Recall: A variational limit What do we want from a variational limit? 1. Convergence of the values 2. Convergence of trajectories 3. Convergence of optimal controls and 4. Possibility to construct near optimal solution for the perturbed system given an optimal solution to the variational limit. 97

98 A question: Under what conditions is a variational limit of 98

99 A Theorem: The conditions are: Regularity (modest) of and Uniform boundedness and controllability of solutions of The set-valued map is Lipschitz 99

100 The Lipschitz condition cannot be dropped: Example (Olivier Alvarez, Martino Bardi): is a polar coordinate 100

101 An issue: How to relate trajectories (say optimal solutions) of the limit problem to the perturbed problem? The answer: If is designed such that approximates the limit Young Measure (in the space of Young Measures), the outcome of the perturbed equation will be a good approximation of the limit (hence of the optimal solution to the perturbed equation. Under the conditions of the theorem his can be done! 101

102 The End of lecure 1 Thanks for the attention See you tomorrow 102

Control of Coupled Slow and Fast Dynamics Zvi Artstein

Control of Coupled Slow and Fast Dynamics Zvi Artstein Control of Coupled Slow and Fast Dynamics Zvi Artstein Presentations in: DIMACS Workshop on Perspectives and Future Directions in Systems and Control Theory = Sontagfest May 23, Piscataway 1 Happy Birthday

More information

Thuong Nguyen. SADCO Internal Review Metting

Thuong Nguyen. SADCO Internal Review Metting Asymptotic behavior of singularly perturbed control system: non-periodic setting Thuong Nguyen (Joint work with A. Siconolfi) SADCO Internal Review Metting Rome, Nov 10-12, 2014 Thuong Nguyen (Roma Sapienza)

More information

An introduction to Mathematical Theory of Control

An introduction to Mathematical Theory of Control An introduction to Mathematical Theory of Control Vasile Staicu University of Aveiro UNICA, May 2018 Vasile Staicu (University of Aveiro) An introduction to Mathematical Theory of Control UNICA, May 2018

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear

More information

A unified framework for input-to-state stability in systems with two time scales

A unified framework for input-to-state stability in systems with two time scales 1 A unified framework for input-to-state stability in systems with two time scales Andrew R. Teel and Luc Moreau and Dragan Nešić Abstract This paper develops a unified framework for studying robustness

More information

The Hopf equation. The Hopf equation A toy model of fluid mechanics

The Hopf equation. The Hopf equation A toy model of fluid mechanics The Hopf equation A toy model of fluid mechanics 1. Main physical features Mathematical description of a continuous medium At the microscopic level, a fluid is a collection of interacting particles (Van

More information

Computing Singularly Perturbed Differential Equations

Computing Singularly Perturbed Differential Equations Computing Singularly Perturbed Differential Equations Sabyasachi Chatterjee Amit Acharya Zvi Artstein June 30, 2017 Abstract A computational tool for coarse-graining nonlinear systems of ordinary differential

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

STOCHASTIC DIFFERENTIAL SYSTEMS WITH MEMORY. Salah-Eldin A. Mohammed

STOCHASTIC DIFFERENTIAL SYSTEMS WITH MEMORY. Salah-Eldin A. Mohammed STOCHASTIC DIFFERENTIAL SYSTEMS WITH MEMORY Salah-Eldin A. Mohammed Research monograph. Preliminary version. Introduction and List of Contents. Tex file sfdebkintrocont.tex. Research supported in part

More information

DIFFERENTIAL GEOMETRY APPLICATIONS TO NONLINEAR OSCILLATORS ANALYSIS

DIFFERENTIAL GEOMETRY APPLICATIONS TO NONLINEAR OSCILLATORS ANALYSIS DIFFERENTIAL GEOMETR APPLICATIONS TO NONLINEAR OSCILLATORS ANALSIS NDES6, J.M. Ginoux and B. Rossetto Laboratoire P.R.O.T.E.E., I.U.T. de Toulon, Université du Sud, B.P., 8957, La Garde Cedex, France e-mail:

More information

Selecting Efficient Correlated Equilibria Through Distributed Learning. Jason R. Marden

Selecting Efficient Correlated Equilibria Through Distributed Learning. Jason R. Marden 1 Selecting Efficient Correlated Equilibria Through Distributed Learning Jason R. Marden Abstract A learning rule is completely uncoupled if each player s behavior is conditioned only on his own realized

More information

A relaxation of the strangeness index

A relaxation of the strangeness index echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se

More information

Control Theory and its Applications

Control Theory and its Applications Control Theory and its Applications E.O. Roxin The University of Rhode Island USA GORDON AND BREACH SCIENCE PUBLISHERS Australia Canada. China France Germany India. Japan Luxembourg Malaysia The Netherlands.

More information

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij

Weak convergence and Brownian Motion. (telegram style notes) P.J.C. Spreij Weak convergence and Brownian Motion (telegram style notes) P.J.C. Spreij this version: December 8, 2006 1 The space C[0, ) In this section we summarize some facts concerning the space C[0, ) of real

More information

Topic 6: Projected Dynamical Systems

Topic 6: Projected Dynamical Systems Topic 6: Projected Dynamical Systems John F. Smith Memorial Professor and Director Virtual Center for Supernetworks Isenberg School of Management University of Massachusetts Amherst, Massachusetts 01003

More information

3. (a) What is a simple function? What is an integrable function? How is f dµ defined? Define it first

3. (a) What is a simple function? What is an integrable function? How is f dµ defined? Define it first Math 632/6321: Theory of Functions of a Real Variable Sample Preinary Exam Questions 1. Let (, M, µ) be a measure space. (a) Prove that if µ() < and if 1 p < q

More information

(Refer Slide Time: 00:32)

(Refer Slide Time: 00:32) Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 12 Scilab simulation of Lotka Volterra

More information

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

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

More information

Finite propagation operators which are Fredholm

Finite propagation operators which are Fredholm Finite propagation operators which are Fredholm Vladimir Rabinovich IPN, Mexico Steffen Roch Darmstadt John Roe Penn State September 20, 2003 The translation algebra Let X be a metric space. We will assume

More information

Relaxation of Euler-Type Discrete-Time Control System 1

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

More information

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II

Ch 5.7: Series Solutions Near a Regular Singular Point, Part II Ch 5.7: Series Solutions Near a Regular Singular Point, Part II! Recall from Section 5.6 (Part I): The point x 0 = 0 is a regular singular point of with and corresponding Euler Equation! We assume solutions

More information

Lecture 4 An Introduction to Stochastic Processes

Lecture 4 An Introduction to Stochastic Processes Lecture 4 An Introduction to Stochastic Processes Prof. Massimo Guidolin Prep Course in Quantitative Methods for Finance August-September 2017 Plan of the lecture Motivation and definitions Filtrations

More information

Basic Measure and Integration Theory. Michael L. Carroll

Basic Measure and Integration Theory. Michael L. Carroll Basic Measure and Integration Theory Michael L. Carroll Sep 22, 2002 Measure Theory: Introduction What is measure theory? Why bother to learn measure theory? 1 What is measure theory? Measure theory is

More information

Mathematical Theory of Control Systems Design

Mathematical Theory of Control Systems Design Mathematical Theory of Control Systems Design by V. N. Afarias'ev, V. B. Kolmanovskii and V. R. Nosov Moscow University of Electronics and Mathematics, Moscow, Russia KLUWER ACADEMIC PUBLISHERS DORDRECHT

More information

Introduction to Empirical Processes and Semiparametric Inference Lecture 08: Stochastic Convergence

Introduction to Empirical Processes and Semiparametric Inference Lecture 08: Stochastic Convergence Introduction to Empirical Processes and Semiparametric Inference Lecture 08: Stochastic Convergence Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations

More information

Some Background Material

Some Background Material Chapter 1 Some Background Material In the first chapter, we present a quick review of elementary - but important - material as a way of dipping our toes in the water. This chapter also introduces important

More information

FROM EQUILIBRIUM TO CHAOS

FROM EQUILIBRIUM TO CHAOS FROM EQUILIBRIUM TO CHAOS Practica! Bifurcation and Stability Analysis RÜDIGER SEYDEL Institut für Angewandte Mathematik und Statistik University of Würzburg Würzburg, Federal Republic of Germany ELSEVIER

More information

Continuous dependence estimates for the ergodic problem with an application to homogenization

Continuous dependence estimates for the ergodic problem with an application to homogenization Continuous dependence estimates for the ergodic problem with an application to homogenization Claudio Marchi Bayreuth, September 12 th, 2013 C. Marchi (Università di Padova) Continuous dependence Bayreuth,

More information

Method of Averaging for Differential Equations on an Infinite Interval

Method of Averaging for Differential Equations on an Infinite Interval Method of Averaging for Differential Equations on an Infinite Interval Theory and Applications SUB Gottingen 7 222 045 71X ;, ' Vladimir Burd Yaroslavl State University Yaroslavl, Russia 2 ' 08A14338 Contents

More information

Solutions Definition 2: a solution

Solutions Definition 2: a solution Solutions As was stated before, one of the goals in this course is to solve, or find solutions of differential equations. In the next definition we consider the concept of a solution of an ordinary differential

More information

Objective. 1 Specification/modeling of the controlled system. 2 Specification of a performance criterion

Objective. 1 Specification/modeling of the controlled system. 2 Specification of a performance criterion Optimal Control Problem Formulation Optimal Control Lectures 17-18: Problem Formulation Benoît Chachuat Department of Chemical Engineering Spring 2009 Objective Determine the control

More information

Lecture Notes Introduction to Ergodic Theory

Lecture Notes Introduction to Ergodic Theory Lecture Notes Introduction to Ergodic Theory Tiago Pereira Department of Mathematics Imperial College London Our course consists of five introductory lectures on probabilistic aspects of dynamical systems,

More information

18.175: Lecture 2 Extension theorems, random variables, distributions

18.175: Lecture 2 Extension theorems, random variables, distributions 18.175: Lecture 2 Extension theorems, random variables, distributions Scott Sheffield MIT Outline Extension theorems Characterizing measures on R d Random variables Outline Extension theorems Characterizing

More information

MAT 578 FUNCTIONAL ANALYSIS EXERCISES

MAT 578 FUNCTIONAL ANALYSIS EXERCISES MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.

More information

Five Mini-Courses on Analysis

Five Mini-Courses on Analysis Christopher Heil Five Mini-Courses on Analysis Metrics, Norms, Inner Products, and Topology Lebesgue Measure and Integral Operator Theory and Functional Analysis Borel and Radon Measures Topological Vector

More information

Morse Theory and Applications to Equivariant Topology

Morse Theory and Applications to Equivariant Topology Morse Theory and Applications to Equivariant Topology Morse Theory: the classical approach Briefly, Morse theory is ubiquitous and indomitable (Bott). It embodies a far reaching idea: the geometry and

More information

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)

8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) 8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)

More information

Monotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005

Monotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005 Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 005 Foreword The aim of this presentation is to give a (primitive) overview of monotone systems and monotone

More information

Population models from PEPA descriptions

Population models from PEPA descriptions Population models from PEPA descriptions Jane Hillston LFCS, The University of Edinburgh, Edinburgh EH9 3JZ, Scotland. Email: jeh@inf.ed.ac.uk 1 Introduction Stochastic process algebras (e.g. PEPA [10],

More information

Auxiliary signal design for failure detection in uncertain systems

Auxiliary signal design for failure detection in uncertain systems Auxiliary signal design for failure detection in uncertain systems R. Nikoukhah, S. L. Campbell and F. Delebecque Abstract An auxiliary signal is an input signal that enhances the identifiability of a

More information

Chapter One. Introduction

Chapter One. Introduction Chapter One Introduction A system is a combination of components or parts that is perceived as a single entity. The parts making up the system may be clearly or vaguely defined. These parts are related

More information

Weak Mean Stability in Random Holomorphic Dynamical Systems Hiroki Sumi Graduate School of Human and Environmental Studies Kyoto University Japan

Weak Mean Stability in Random Holomorphic Dynamical Systems Hiroki Sumi Graduate School of Human and Environmental Studies Kyoto University Japan Weak Mean Stability in Random Holomorphic Dynamical Systems Hiroki Sumi Graduate School of Human and Environmental Studies Kyoto University Japan E-mail: sumi@math.h.kyoto-u.ac.jp http://www.math.h.kyoto-u.ac.jp/

More information

Variational Methods and Optimal Control

Variational Methods and Optimal Control Variational Methods and Optimal Control A/Prof. Matthew Roughan July 26, 2010 lecture01 Introduction What is the point of this course? Revision Example 1: The money pit. Example 2: Catenary: shape of a

More information

On Kusuoka Representation of Law Invariant Risk Measures

On Kusuoka Representation of Law Invariant Risk Measures MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 1, February 213, pp. 142 152 ISSN 364-765X (print) ISSN 1526-5471 (online) http://dx.doi.org/1.1287/moor.112.563 213 INFORMS On Kusuoka Representation of

More information

Measure and integration

Measure and integration Chapter 5 Measure and integration In calculus you have learned how to calculate the size of different kinds of sets: the length of a curve, the area of a region or a surface, the volume or mass of a solid.

More information

Multidimensional Divide and Conquer 1 Skylines

Multidimensional Divide and Conquer 1 Skylines Yufei Tao ITEE University of Queensland The next few lectures will be dedicated to an important technique: divide and conquer. You may have encountered the technique in an earlier algorithm course, but

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

The Euler Method for Linear Control Systems Revisited

The Euler Method for Linear Control Systems Revisited The Euler Method for Linear Control Systems Revisited Josef L. Haunschmied, Alain Pietrus, and Vladimir M. Veliov Research Report 2013-02 March 2013 Operations Research and Control Systems Institute of

More information

Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms

Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms Almost Sure Convergence of Two Time-Scale Stochastic Approximation Algorithms Vladislav B. Tadić Abstract The almost sure convergence of two time-scale stochastic approximation algorithms is analyzed under

More information

Lecture 10: Singular Perturbations and Averaging 1

Lecture 10: Singular Perturbations and Averaging 1 Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.243j (Fall 2003): DYNAMICS OF NONLINEAR SYSTEMS by A. Megretski Lecture 10: Singular Perturbations and

More information

Hausdorff Measure. Jimmy Briggs and Tim Tyree. December 3, 2016

Hausdorff Measure. Jimmy Briggs and Tim Tyree. December 3, 2016 Hausdorff Measure Jimmy Briggs and Tim Tyree December 3, 2016 1 1 Introduction In this report, we explore the the measurement of arbitrary subsets of the metric space (X, ρ), a topological space X along

More information

A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS

A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS International Journal of Scientific and Research Publications Volume 5 Issue 9 September 2015 1 A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS Mansour Hassan Mansour * M A Bashir ** * Department

More information

Near-Potential Games: Geometry and Dynamics

Near-Potential Games: Geometry and Dynamics Near-Potential Games: Geometry and Dynamics Ozan Candogan, Asuman Ozdaglar and Pablo A. Parrilo January 29, 2012 Abstract Potential games are a special class of games for which many adaptive user dynamics

More information

Consistency and Asymptotic Normality for Equilibrium Models with Partially Observed Outcome Variables

Consistency and Asymptotic Normality for Equilibrium Models with Partially Observed Outcome Variables Consistency and Asymptotic Normality for Equilibrium Models with Partially Observed Outcome Variables Nathan H. Miller Georgetown University Matthew Osborne University of Toronto November 25, 2013 Abstract

More information

1.4 The Jacobian of a map

1.4 The Jacobian of a map 1.4 The Jacobian of a map Derivative of a differentiable map Let F : M n N m be a differentiable map between two C 1 manifolds. Given a point p M we define the derivative of F at p by df p df (p) : T p

More information

Procedia Computer Science 00 (2011) 000 6

Procedia Computer Science 00 (2011) 000 6 Procedia Computer Science (211) 6 Procedia Computer Science Complex Adaptive Systems, Volume 1 Cihan H. Dagli, Editor in Chief Conference Organized by Missouri University of Science and Technology 211-

More information

Overview of the Seminar Topic

Overview of the Seminar Topic Overview of the Seminar Topic Simo Särkkä Laboratory of Computational Engineering Helsinki University of Technology September 17, 2007 Contents 1 What is Control Theory? 2 History

More information

16 1 Basic Facts from Functional Analysis and Banach Lattices

16 1 Basic Facts from Functional Analysis and Banach Lattices 16 1 Basic Facts from Functional Analysis and Banach Lattices 1.2.3 Banach Steinhaus Theorem Another fundamental theorem of functional analysis is the Banach Steinhaus theorem, or the Uniform Boundedness

More information

Overview of normed linear spaces

Overview of normed linear spaces 20 Chapter 2 Overview of normed linear spaces Starting from this chapter, we begin examining linear spaces with at least one extra structure (topology or geometry). We assume linearity; this is a natural

More information

Models for Control and Verification

Models for Control and Verification Outline Models for Control and Verification Ian Mitchell Department of Computer Science The University of British Columbia Classes of models Well-posed models Difference Equations Nonlinear Ordinary Differential

More information

MA651 Topology. Lecture 9. Compactness 2.

MA651 Topology. Lecture 9. Compactness 2. MA651 Topology. Lecture 9. Compactness 2. This text is based on the following books: Topology by James Dugundgji Fundamental concepts of topology by Peter O Neil Elements of Mathematics: General Topology

More information

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0

More information

Optimal Control of Weakly Coupled Systems and Applications

Optimal Control of Weakly Coupled Systems and Applications Optimal Control of Weakly Coupled Systems and Applications HIGH ACCURACY TECHNIQUES Z. Gajić, M-T. Lim, D. Skatarić W-C. Su, and V. Kecman Taylor & Francis (CRC Press, Dekker) 2008 Preface This book is

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

1.3.1 Definition and Basic Properties of Convolution

1.3.1 Definition and Basic Properties of Convolution 1.3 Convolution 15 1.3 Convolution Since L 1 (R) is a Banach space, we know that it has many useful properties. In particular the operations of addition and scalar multiplication are continuous. However,

More information

ME751 Advanced Computational Multibody Dynamics. September 14, 2016

ME751 Advanced Computational Multibody Dynamics. September 14, 2016 ME751 Advanced Computational Multibody Dynamics September 14, 2016 Quote of the Day My own business always bores me to death; I prefer other people's. -- Oscar Wilde 2 Looking Ahead, Friday Need to wrap

More information

Convergence results for solutions of a first-order differential equation

Convergence results for solutions of a first-order differential equation vailable online at www.tjnsa.com J. Nonlinear ci. ppl. 6 (213), 18 28 Research rticle Convergence results for solutions of a first-order differential equation Liviu C. Florescu Faculty of Mathematics,

More information

Dynamical Systems and Ergodic Theory PhD Exam Spring Topics: Topological Dynamics Definitions and basic results about the following for maps and

Dynamical Systems and Ergodic Theory PhD Exam Spring Topics: Topological Dynamics Definitions and basic results about the following for maps and Dynamical Systems and Ergodic Theory PhD Exam Spring 2012 Introduction: This is the first year for the Dynamical Systems and Ergodic Theory PhD exam. As such the material on the exam will conform fairly

More information

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3.

In terms of measures: Exercise 1. Existence of a Gaussian process: Theorem 2. Remark 3. 1. GAUSSIAN PROCESSES A Gaussian process on a set T is a collection of random variables X =(X t ) t T on a common probability space such that for any n 1 and any t 1,...,t n T, the vector (X(t 1 ),...,X(t

More information

arxiv: v1 [math.ds] 23 Jan 2009

arxiv: v1 [math.ds] 23 Jan 2009 Discontinuous Dynamical Systems A tutorial on solutions, nonsmooth analysis, and stability arxiv:0901.3583v1 [math.ds] 23 Jan 2009 Jorge Cortés January 23, 2009 Discontinuous dynamical systems arise in

More information

Mathematical Biology - Lecture 1 - general formulation

Mathematical Biology - Lecture 1 - general formulation Mathematical Biology - Lecture 1 - general formulation course description Learning Outcomes This course is aimed to be accessible both to masters students of biology who have a good understanding of the

More information

Introduction to Atmosphere, Ocean, and Climate Dynamics

Introduction to Atmosphere, Ocean, and Climate Dynamics 12.003 Introduction to Atmosphere, Ocean, and Climate Dynamics Instructor Paul O Gorman, 54-1712 pog@mit.edu Teaching Assistant Elezhan Zhakiya (e.zhakiya@gmail.com) Webpage http://www.mit.edu/~pog/12.003

More information

Fragmentability and σ-fragmentability

Fragmentability and σ-fragmentability F U N D A M E N T A MATHEMATICAE 143 (1993) Fragmentability and σ-fragmentability by J. E. J a y n e (London), I. N a m i o k a (Seattle) and C. A. R o g e r s (London) Abstract. Recent work has studied

More information

CHAPTER 1. Introduction

CHAPTER 1. Introduction CHAPTER 1 Introduction Linear geometric control theory was initiated in the beginning of the 1970 s, see for example, [1, 7]. A good summary of the subject is the book by Wonham [17]. The term geometric

More information

I R TECHNICAL RESEARCH REPORT. Risk-Sensitive Probability for Markov Chains. by Vahid R. Ramezani and Steven I. Marcus TR

I R TECHNICAL RESEARCH REPORT. Risk-Sensitive Probability for Markov Chains. by Vahid R. Ramezani and Steven I. Marcus TR TECHNICAL RESEARCH REPORT Risk-Sensitive Probability for Markov Chains by Vahid R. Ramezani and Steven I. Marcus TR 2002-46 I R INSTITUTE FOR SYSTEMS RESEARCH ISR develops, applies and teaches advanced

More information

Chapter 18. Remarks on partial differential equations

Chapter 18. Remarks on partial differential equations Chapter 8. Remarks on partial differential equations If we try to analyze heat flow or vibration in a continuous system such as a building or an airplane, we arrive at a kind of infinite system of ordinary

More information

CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION

CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION CANARDS AND HORSESHOES IN THE FORCED VAN DER POL EQUATION WARREN WECKESSER Department of Mathematics Colgate University Hamilton, NY 3346 E-mail: wweckesser@mail.colgate.edu Cartwright and Littlewood discovered

More information

An introduction to some aspects of functional analysis

An introduction to some aspects of functional analysis An introduction to some aspects of functional analysis Stephen Semmes Rice University Abstract These informal notes deal with some very basic objects in functional analysis, including norms and seminorms

More information

FEEDBACK DIFFERENTIAL SYSTEMS: APPROXIMATE AND LIMITING TRAJECTORIES

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

More information

Lecture 1: Pragmatic Introduction to Stochastic Differential Equations

Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic

More information

The Codimension of the Zeros of a Stable Process in Random Scenery

The Codimension of the Zeros of a Stable Process in Random Scenery The Codimension of the Zeros of a Stable Process in Random Scenery Davar Khoshnevisan The University of Utah, Department of Mathematics Salt Lake City, UT 84105 0090, U.S.A. davar@math.utah.edu http://www.math.utah.edu/~davar

More information

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively);

(x k ) sequence in F, lim x k = x x F. If F : R n R is a function, level sets and sublevel sets of F are any sets of the form (respectively); STABILITY OF EQUILIBRIA AND LIAPUNOV FUNCTIONS. By topological properties in general we mean qualitative geometric properties (of subsets of R n or of functions in R n ), that is, those that don t depend

More information

Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis

Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis Phase Model for the relaxed van der Pol oscillator and its application to synchronization analysis Mimila Prost O. Collado J. Automatic Control Department, CINVESTAV IPN, A.P. 4 74 Mexico, D.F. ( e mail:

More information

Stochastic-Process Limits

Stochastic-Process Limits Ward Whitt Stochastic-Process Limits An Introduction to Stochastic-Process Limits and Their Application to Queues With 68 Illustrations Springer Contents Preface vii 1 Experiencing Statistical Regularity

More information

SINGULAR PERTURBATIONS IN OPTIMAL CONTROL*

SINGULAR PERTURBATIONS IN OPTIMAL CONTROL* ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 6, Number 4, Fall 1976 SINGULAR PERTURBATIONS IN OPTIMAL CONTROL* PETAR V. KOKOTOVIC ABSTRACT. Optimal control is a new area for applications of singular perturbation

More information

Introduction. Christophe Prange. February 9, This set of lectures is motivated by the following kind of phenomena:

Introduction. Christophe Prange. February 9, This set of lectures is motivated by the following kind of phenomena: Christophe Prange February 9, 206 This set of lectures is motivated by the following kind of phenomena: sin(x/ε) 0, while sin 2 (x/ε) /2. Therefore the weak limit of the product is in general different

More information

Mathematical Economics (ECON 471) Lecture 3 Calculus of Several Variables & Implicit Functions

Mathematical Economics (ECON 471) Lecture 3 Calculus of Several Variables & Implicit Functions Mathematical Economics (ECON 471) Lecture 3 Calculus of Several Variables & Implicit Functions Teng Wah Leo 1 Calculus of Several Variables 11 Functions Mapping between Euclidean Spaces Where as in univariate

More information

OPTIMAL CONTROL AND RELATED TOPICS Czesław Olech in memoriam Warsaw, October Abstracts

OPTIMAL CONTROL AND RELATED TOPICS Czesław Olech in memoriam Warsaw, October Abstracts OPTIMAL CONTROL AND RELATED TOPICS Czesław Olech in memoriam Warsaw, 20-21 October 2016 Abstracts Zvi Artstein (Weizmann Institute) Olech s extremals in the singular perturbations limit We shall examine

More information

ON STOCHASTIC STRUCTURAL TOPOLOGY OPTIMIZATION

ON STOCHASTIC STRUCTURAL TOPOLOGY OPTIMIZATION ON STOCHASTIC STRUCTURAL TOPOLOGY OPTIMIZATION A. EVGRAFOV and M. PATRIKSSON Department of Mathematics, Chalmers University of Technology, SE-4 96 Göteborg, Sweden J. PETERSSON Department of Mechanical

More information

HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES ABSTRACT

HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES ABSTRACT HYPERBOLIC DYNAMICAL SYSTEMS AND THE NONCOMMUTATIVE INTEGRATION THEORY OF CONNES Jan Segert ABSTRACT We examine hyperbolic differentiable dynamical systems in the context of Connes noncommutative integration

More information

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems. Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics

More information

SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES

SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES SUMMARY OF RESULTS ON PATH SPACES AND CONVERGENCE IN DISTRIBUTION FOR STOCHASTIC PROCESSES RUTH J. WILLIAMS October 2, 2017 Department of Mathematics, University of California, San Diego, 9500 Gilman Drive,

More information

Notes for Functional Analysis

Notes for Functional Analysis Notes for Functional Analysis Wang Zuoqin (typed by Xiyu Zhai) November 6, 2015 1 Lecture 18 1.1 The convex hull Let X be any vector space, and E X a subset. Definition 1.1. The convex hull of E is the

More information

Effective Field Theory for Nuclear Physics! Akshay Vaghani! Mississippi State University!

Effective Field Theory for Nuclear Physics! Akshay Vaghani! Mississippi State University! Effective Field Theory for Nuclear Physics! Akshay Vaghani! Mississippi State University! Overview! Introduction! Basic ideas of EFT! Basic Examples of EFT! Algorithm of EFT! Review NN scattering! NN scattering

More information

A Study of the Van der Pol Equation

A Study of the Van der Pol Equation A Study of the Van der Pol Equation Kai Zhe Tan, s1465711 September 16, 2016 Abstract The Van der Pol equation is famous for modelling biological systems as well as being a good model to study its multiple

More information

1 )(y 0) {1}. Thus, the total count of points in (F 1 (y)) is equal to deg y0

1 )(y 0) {1}. Thus, the total count of points in (F 1 (y)) is equal to deg y0 1. Classification of 1-manifolds Theorem 1.1. Let M be a connected 1 manifold. Then M is diffeomorphic either to [0, 1], [0, 1), (0, 1), or S 1. We know that none of these four manifolds are not diffeomorphic

More information

ROBUSTNESS COMPARISON OF CONTROL SYSTEMS FOR A NUCLEAR POWER PLANT

ROBUSTNESS COMPARISON OF CONTROL SYSTEMS FOR A NUCLEAR POWER PLANT Control 004, University of Bath, UK, September 004 ROBUSTNESS COMPARISON OF CONTROL SYSTEMS FOR A NUCLEAR POWER PLANT L Ding*, A Bradshaw, C J Taylor Lancaster University, UK. * l.ding@email.com Fax: 0604

More information

2 A Model, Harmonic Map, Problem

2 A Model, Harmonic Map, Problem ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or

More information

MTNS Conference, Polynomial Inequalities and Applications. Kyoto, July

MTNS Conference, Polynomial Inequalities and Applications. Kyoto, July MTNS Conference, Polynomial Inequalities and Applications. Kyoto, July 24-28 2006. July 20, 2006 Salma Kuhlmann 1, Research Center Algebra, Logic and Computation University of Saskatchewan, McLean Hall,

More information

MAT 544 Problem Set 2 Solutions

MAT 544 Problem Set 2 Solutions MAT 544 Problem Set 2 Solutions Problems. Problem 1 A metric space is separable if it contains a dense subset which is finite or countably infinite. Prove that every totally bounded metric space X is separable.

More information