Hybrid Control and Switched Systems. Lecture #3 What can go wrong? Trajectories of hybrid systems

Similar documents
EXERCISES FOR SECTION 1.5

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Lecture 20: Riccati Equations and Least Squares Feedback Control

Chapter 1 Fundamental Concepts

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Intermediate Differential Equations Review and Basic Ideas

Crossing the Bridge between Similar Games

Predator - Prey Model Trajectories and the nonlinear conservation law

From Particles to Rigid Bodies

15. Vector Valued Functions

(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)

Solutions to Assignment 1

Ordinary Differential Equations

t dt t SCLP Bellman (1953) CLP (Dantzig, Tyndall, Grinold, Perold, Anstreicher 60's-80's) Anderson (1978) SCLP

Four Generations of Higher Order Sliding Mode Controllers. L. Fridman Universidad Nacional Autonoma de Mexico Aussois, June, 10th, 2015

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

Hybrid Control and Switched Systems. Lecture #1 Hybrid systems are everywhere: Examples

LAPLACE TRANSFORM AND TRANSFER FUNCTION

Basic Circuit Elements Professor J R Lucas November 2001

An Introduction to Malliavin calculus and its applications

4.6 One Dimensional Kinematics and Integration

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

6. Stochastic calculus with jump processes

Lecture 4 Notes (Little s Theorem)

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

Chapter 7 Response of First-order RL and RC Circuits

2. Nonlinear Conservation Law Equations

Physics 221 Fall 2008 Homework #2 Solutions Ch. 2 Due Tues, Sept 9, 2008

Hamilton Jacobi equations

f t te e = possesses a Laplace transform. Exercises for Module-III (Transform Calculus)

Math 2214 Solution Test 1A Spring 2016

arxiv: v1 [math.ca] 15 Nov 2016

Existence of multiple positive periodic solutions for functional differential equations

On-line Adaptive Optimal Timing Control of Switched Systems

Math 333 Problem Set #2 Solution 14 February 2003

4 Sequences of measurable functions

Essential Microeconomics : OPTIMAL CONTROL 1. Consider the following class of optimization problems

STABILITY OF RESET SWITCHING SYSTEMS

Guest Lectures for Dr. MacFarlane s EE3350 Part Deux

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

E β t log (C t ) + M t M t 1. = Y t + B t 1 P t. B t 0 (3) v t = P tc t M t Question 1. Find the FOC s for an optimum in the agent s problem.

Chapter 6. Systems of First Order Linear Differential Equations

CONTRIBUTION TO IMPULSIVE EQUATIONS

Math 527 Lecture 6: Hamilton-Jacobi Equation: Explicit Formulas

4.5 Constant Acceleration

10. State Space Methods

EE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

The Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales

3.6 Derivatives as Rates of Change

Math 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.

t 2 B F x,t n dsdt t u x,t dxdt

MATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence

Variational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations

Class Meeting # 10: Introduction to the Wave Equation

Math 334 Fall 2011 Homework 11 Solutions

Basic definitions and relations

= ( ) ) or a system of differential equations with continuous parametrization (T = R

THE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI

Properties Of Solutions To A Generalized Liénard Equation With Forcing Term

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

Solution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration

Chapter 2. First Order Scalar Equations

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is

Chapter 3 Boundary Value Problem

POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION

MA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions

THE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).

2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes

Math Final Exam Solutions

CHAPTER 12 DIRECT CURRENT CIRCUITS

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

1 Subdivide the optimization horizon [t 0,t f ] into n s 1 control stages,

AP CALCULUS AB 2003 SCORING GUIDELINES (Form B)

arxiv: v1 [math.fa] 9 Dec 2018

ME 391 Mechanical Engineering Analysis

5. Stochastic processes (1)

Representing a Signal. Continuous-Time Fourier Methods. Linearity and Superposition. Real and Complex Sinusoids. Jean Baptiste Joseph Fourier

Projection-Based Optimal Mode Scheduling

Direct Current Circuits. February 19, 2014 Physics for Scientists & Engineers 2, Chapter 26 1

Question 1: Question 2: Topology Exercise Sheet 3

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

u(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x

i L = VT L (16.34) 918a i D v OUT i L v C V - S 1 FIGURE A switched power supply circuit with diode and a switch.

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Designing Information Devices and Systems I Spring 2019 Lecture Notes Note 17

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

and v y . The changes occur, respectively, because of the acceleration components a x and a y

( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:

CHE302 LECTURE IV MATHEMATICAL MODELING OF CHEMICAL PROCESS. Professor Dae Ryook Yang

6.003 Homework 1. Problems. Due at the beginning of recitation on Wednesday, February 10, 2010.

Continuous Time. Time-Domain System Analysis. Impulse Response. Impulse Response. Impulse Response. Impulse Response. ( t) + b 0.

Notes on Kalman Filtering

Utility maximization in incomplete markets

Ann. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:

MEI Mechanics 1 General motion. Section 1: Using calculus

Transcription:

Hybrid Conrol and Swiched Sysems Lecure #3 Wha can go wrong? Trajecories of hybrid sysems João P. Hespanha Universiy of California a Sana Barbara Summary 1. Trajecories of hybrid sysems: Soluion o a hybrid sysem Execuion of a hybrid sysem 2. Degeneracies Finie escape ime Chaering Zeno rajecories Non-coninuous dependency on iniial condiions 1

Hybrid Auomaon (deerminisic) Q se of discree saes R n coninuous sae-space f : Q R n R n vecor field ϕ : Q R n Q discree ransiion ρ : Q R n R n rese map x ú ρ(q 1,x ) mode q 2 ϕ(q 1,x ) = q 2? mode q 1 mode q 3 ϕ(q 1,x ) = q 3? x ú ρ(q 1,x ) Q R n f : Q R n R n Φ : Q R n Q R n Hybrid Auomaon se of discree saes coninuous sae-space vecor field discree ransiion (& rese map) x ú Φ 2 (q 1,x ) mode q 2 Φ 1 (q 1,x ) = q 2? mode q 1 mode q 3 Φ 1 (q 1,x ) = q 3? x ú Φ 2 (q 1,x ) Compac represenaion of a hybrid auomaon 2

Soluion o a hybrid auomaon x ú Φ 2 (q 1,x ) Φ 1 (q 1,x ) = q 2? mode q 2 mode q 1 Definiion: A soluion o he hybrid auomaon is a pair of righ-coninuous signals x : [0, ) R n q : [0, ) Q such ha 1. x is piecewise differeniable & q is piecewise consan 2. on any inerval ( 1, 2 ) on which q is consan and x coninuous coninuous evoluion 3. discree ransiions Example #4: Invered pendulum swing-up u [-1,1] θ Hybrid conroller: 1 s pump/remove energy ino/from he sysem by applying maximum force, unil E 0 2 nd wai unil pendulum is close o he uprigh posiion 3 h nex o uprigh posiion use feedback linearizaion conroller remove energy E [ ε,ε]? E < ε? wai sabilize E>ε? pump energy E [ ε,ε]? ω + θ δ? 3

Example #4: Invered pendulum swing-up pump energy wai sabilize Ε ε? ω + θ δ? Q = { pump, wai, sab } R 2 = coninuous sae-space (pump,x) q = pump,e < ² (wai,x) q = pump,e ² Φ(q,x)= (wai,x) q =wai, ω + θ > δ (sab,x) q =wai, ω + θ δ (sab,x) q =sab Example #4: Invered pendulum swing-up pump energy wai sabilize Ε ε? ω + θ δ? Q = { pump, wai, sab } R 2 = coninuous sae-space δ ω + θ ε E q = wai q = sab q = pump Wha if we already have ω + θ δ when E reaches ε? 4

Example #4: Invered pendulum swing-up pump energy wai sabilize Q = { pump, wai, sab } δ ε ω + θ E q = pump * Ε ε? q = wai ω + θ δ? R 2 = coninuous sae-space (pump,x) q =pump,e < ² (wai,x) q =pump,e ² Φ(q, x) = (wai,x) q =wai, ω + θ > ² (sab,x) q =wai, ω + θ ² (sab,x) q =sab q = sab q( * ) =? In general, we may need more han one value for he sae a each insan of ime Hybrid signals Definiion: A hybrid ime rajecory is a (finie or infinie) sequence of closed inervals τ = { [τ i,τ' i ] : τ i τ' i, τ' i = τ i+1, i =1,2, } (if τ is finie he las inerval may by open on he righ) T se of hybrid ime rajecories Definiion: For a given τ = { [τ i,τ' i ] : τ i τ' i, τ i+1 = τ' i, i =1,2, } T a hybrid signal defined on τ wih values on X is a sequence of funcions x = {x i : [τ i,τ' i ] X i =1,2, } x : τ X hybrid signal defined on τ wih values on X E.g., τ ú { [0,1], [1,1], [1,1], [1,+ ) }, x ú { /2, 1/4, 3/4, 1/ } /2 3/4 1/4 1/ [0,1] [1,1] [1,1] [1,+ ) A hybrid signal can ake muliple values for he same ime-insan 5

Execuion of a hybrid auomaon x ú Φ 2 (q 1,x ) Φ 1 (q 1,x ) = q 2? mode q 2 mode q 1 Definiion: An execuion of he hybrid auomaon is a pair of hybrid signals x : τ R n q : τ Q τ = { [τ i,τ' i ] : i =1,2, } T such ha 1. on any [τ i,τ' i ] τ, q i is consan and coninuous evoluion 2. discree ransiions Example #4: Invered pendulum swing-up pump energy wai sabilize Ε ε? ω + θ δ? Q = { pump, wai, sab } R 2 = coninuous sae-space δ ε ω + θ * q = sab τ ú { [0, * ], [ *, * ], [ *,+ ) } q ú { pump, wai, sab } x ú {,, } E q = pump q = wai 1. There are oher conceps of soluion ha also avoid his problem [Teel 2005] 2. Could one fix his hybrid sysem o sill work wih he usual noion of soluion? 6

Wha can go wrong? 1. Problems in he coninuous evoluion : exisence uniqueness finie escape 2. Problems in he hybrid execuion: Chaering Zeno 3. Non-coninuous dependency on iniial condiions Exisence f(x) disconinuous x There is no soluion o his differenial equaion ha sars wih x(0) = 0 Why? one any inerval [0,ε) x canno: remain zero, become posiive, or become negaive. (x = 0 would make some sense) Theorem [Exisence of soluion] If f : R n R n is coninuous, hen x 0 R n here exiss a leas one soluion wih x(0) = x 0, defined on some inerval [0,ε) 7

Uniqueness f(x) derivaive a 0 There are muliple soluions ha sar a x(0) = 0, e.g., x x Definiions: A funcion f : R n R n is Lipschiz coninuous if in any bounded subse of S of R n here exiss a consan c such ha (f is differeniable almos everywhere and he derivaive is bounded on any bounded se) Theorem [Uniqueness of soluion] If f : R n R n is Lipschiz coninuous, hen x 0 R n here a single soluion wih x(0) = x 0, defined on some inerval [0,ε) Finie escape ime x Any soluion ha does no sar a x(0) = 0 is of he form unbounded derivaive T -1 T finie escape soluion x ends o in finie ime Definiions: A funcion f : R n R n is globally Lipschiz coninuous if here exiss a consan c such ha f grows no faser han linearly Theorem [Uniqueness of soluion] If f : R n R n is globally Lipschiz coninuous, hen x 0 R n here a single soluion wih x(0) = x 0, defined on [0, ) 8

Degenerae execuions due o problems in he coninuous evoluion x 1? q = 1 q = 2 τ = { [0,1], [1,2] } q = {q 1 () = 1, q 2 () = 2 } x = {x 1 () =, () = 2 } x x 1 () () 1 2 τ = { [0,1], [1,2) } q = {q 1 () = 1, q 2 () = 2 } x = {x 1 () =, () = 1/(2 ) } x 1? q = 1 q = 2 x x 1 () () 1 2 Chaering There is no soluion pas x = 0 x () 1 q = 1 x 0? q = 2 bu here is an execuion x < 0? q 1 () q 2k+1 () k 1 τ = { [0,1], [1], [1], [1], } q = {q 1 () = 2, q 2 () = 1, q 1 () = 2, } x = {x 1 () = 1, () = 0, () = 0, } x 1 () q 2k () k 1 x k () k 2 1 Chaering execuion τ is infinie bu afer some ime, all inervals are singleons 9

Example #3: Semi-auomaic ransmission v() { up, down, keep } drivers inpu (discree) v = up or ω ω 2? v = up or ω ω 3? v = up or ω ω 4? g = 1 g = 2 g = 3 g = 4 v = down or ω ϖ 1? v = down or ω ϖ 2? v = down or ω ϖ 3? ω 2 g = 1 ϖ 1 g = 2 ω 3 ϖ 2 g = 3 ϖ 3 g = 4 ω 4 If he driver ses v() = up * and ω( ) ϖ 1 one ges chaering. For ever? Example #1: Bouncing ball (Zeno execuion) y g Free fall Collision c [0,1) energy absorbed a impac x 1 = 0 & <0? ú c 10

Zeno soluion x 1 = 0 & <0? 1 2 3 ú c Zeno soluion x 1 = 0 & = 0? x 1 = 0 & <0? 1 2 3 ú c 11

Zeno execuion x 1 = 0 & <0? 1 2 3 ú c Zeno execuion τ is infinie bu he execuion does no exend o = + Zeno ime τ ú sup s τ s Example #9: Waer ank sysem (regularizaion) inflow λ goal preven boh anks from becoming empying h 2? q = 1 q = 2 x 1 h 1 h 2 ouflow μ 1 ouflow μ 2 x 1 h 1? x 1 h 1 = h 2 q = 1 q = 2 q = 1 12

inflow λ Temporal regularizaion Example #9: Waer ank sysem h 2? q = 2 x 1 h 1 h 2 x 1 h 1? ouflow μ 1 ouflow μ 2 only jump afer a minimum ime has elapsed h 2, τ > ε? q = 1 τ ú 0 q = 2 τ ú 0 x 1 h 1, τ > ε? inflow λ Spaial regularizaion Example #9: Waer ank sysem h 2? q = 2 x 1 ouflow μ 1 h 1 h 2 ouflow μ 2 only jump afer a minimum change in he coninuous sae has occurred h 2, x 1 x 3 + x 4 > ε? q = 1 x 1 h 1? x 3 ú x 1, x 4 ú q = 2 x 3 ú x 1, x 4 ú x 1 h 1, x 1 x 3 + x 4 > ε? 13

Coninuiy wih respec o iniial condiions Theorem [Uniqueness & coninuiy of soluion] If f : R n R n is Lipschiz coninuous, hen x 0 R n here a single soluion wih x(0) = x 0, defined on some inerval [0,ε) Moreover, given any T <, and wo soluions x 1, ha exis on [0,T]: ε > 0 δ > 0 : x 1 (0) (0) δ x 1 () () ε [0,T] value of he soluion on he inerval [0,T] is coninuous wih respec o he iniial condiions Disconinuiy wih respec o iniial condiions x 1 0 & > 0? mode q 1 ú 1 mode q 2 disconinuiy of he rese map x 1 0 & 0? ú 2 2 ε 1 x 1 1 1 no maer how close o zero (0) = ε > 0 is, (1) = 1 x 1 if (0) = 0 hen (1) = 2 14

Disconinuiy wih respec o iniial condiions x 1 0 & > 0? mode q 1 x 1 0 & 0? ẋ 1 =1 ẋ 2 =0 mode q 2 mode q 3 ẋ 1 =0 ẋ 1 =1 ẋ 2 = 1 ẋ 2 =1 ε 1 x 1 ε 1 x 1 no maer how close o zero (0) = ε > 0 is, (2) = ε 1 if (0) = 0 hen (2) = 1 problem arises from disconinuiy of he ransiion funcion Nex class 1. Numerical simulaion of hybrid auomaa simulaions of ODEs zero-crossing deecion 2. Simulaors Simulink Saeflow SHIFT Modelica Follow-up homework Find condiions for he exisence of soluion o a hybrid sysem 15