Introduction to Geometric Control
|
|
- Bartholomew Sharp
- 5 years ago
- Views:
Transcription
1 Introduction to Geometric Control Relative Degree Consider the square (no of inputs = no of outputs = p) affine control system ẋ = f(x) + g(x)u = f(x) + [ g (x),, g p (x) ] u () h (x) y = h(x) = (2) h p (x) (This system is assumed to have an equilibrium or operating point corresponding to u e at x e = with y e = ) Recall that for the i-th output y i, the relative degree r i is the first time - derivative of y i for which u appears in the expression for the derivative SISO For example, for a SISO system, this translates to L g L f k h(x) = dk y dt = L k f kh(x), k =, 2,, r (3) L g L f r h(x e ) dr y dt r = L f rh(x) + L gl f r h(x)u (4) For nonlinear systems, the relative degree is a local concept, defined in some neighborhood of x = x e If conditions (3) and (4) hold globally, we say that the system has a global relative degree r 2 MIMO The MIMO system ()-(2) has relative degree {r, r 2,, r p } at x = x e if L gj L f k h i (x) =, i, j p, k =, 2,, r i (5) dk y i dt k = L f kh i (x), i p, k =, 2,, r
2 for all x in a neighbourhood of x = x e The p p matrix N(x) is invertible at x = x e where [ ] y (r L i) g L f r h (x) L gp L f r h (x) i N(x) = = u (6) j i,j L g L f rp h p (x) L gp L f rp h p (x) d r y dt r L f r h (x) = +N(x)u d rp y p L dt rp f rph p (x) }{{} G(x) Condition (6) is the MIMO generalisation of the condition (4) in the SISO case If r = r 2 = = r p, we say that the system has a uniform relative degree r 2 Normal Form 2 SISO For a SISO system of relative degree r, consider the state transformation ( ) ( ) ξ Φ(x) = z =, Φ(x η e ) = where ξ R r, η R n r, defined in part as follows: ξ = Φ (x) = y = h(x) ξ 2 = Φ 2 (x) = ẏ = L f h(x) ξ r = Φ r (x) = y (r ) = L f r h(x) (7) 2
3 Since condition (4) holds, this immediately gives that ξ = ξ 2 ξ 2 = ξ 3 (8) ξ r = ξ r ξ r = L f rh(x) + L }{{} g L f r h(x) u }{{} A(x) B(x) It may be shown that the following holds Lemma (Linear independence of output derivatives): The row vectors { Φ (x x e), Φ 2 (x x e),, Φr (x x e)} are linearly independent Hence, if r = n, then Φ(x) is indeed a well-defined transformation in a neighbourhood of x = x e since its Jacobian is invertible On the other hand, if r < n, it is always possible to choose Φ r+ (x), Φ r+2 (x),, Φ n (x) so that Φ(x) has a non-singular Jacobian at x e and so is indeed a well-defined transformation in a neighbourhood of x = x e Hence the dynamics of the remaining transformed state variables η may be described by η i = Φ i (x), i = r +, r + 2,, n η i = Φ i (f(x) + g(x)u), x = L f Φ }{{} i (x) + L g Φ i D i i = r +, r + 2,, n }{{} E i (x)u, i = r +, r + 2,, n (9) Moreover it is always possible to choose Φ r+ (x), Φ r+2 (x),, Φ n (x) so that Lemma 2 (Non-dependence of remnant states on the input): L g Φ i (x) =, for all x in a neighbourhood of x e And so Eq (9 ) becomes i = r +, r + 2,, n η i = D i (x), i = r +, r + 2,, n () 3
4 Using the fact x = Φ (z), the state equations (Eq (8) and Eq () ) and output equation may be written in normal form: ξ = ξ 2 ξ 2 = ξ 3 ξ r = ξ r ξ r = A(Φ (z)) + B(Φ (z))u () η r+ = D r+ (Φ (z)) η r+2 = D r+2 (Φ (z)) η n = D n (Φ (z)) y = ξ Remark: Though Lemma 2 holds in a theoretical sense, it may be difficult to find the required transformation In which case, the state and output equations (not in normal form) are ξ = ξ 2 ξ 2 = ξ 3 ξ r = ξ r ξ r = A(Φ (z)) + B(Φ (z))u (2) η r+ = D r+ (Φ (z)) + E r+ (Φ (z))u η r+2 = D r+2 (Φ (z)) + E r+2 (Φ (z))u η n = D n (Φ (z)) + E n (Φ (z))u y = ξ 4
5 3 I/O Feedback Linearisation 3 SISO By definition of relative degree B() (see Eq ()) Hence by using the cancellation feedback control A(Φ (z)) + B(Φ (z))u = v u = A(Φ (z)) + v B(Φ (z)) (3) the ξ-subsystem can be converted into a completely controllable linear subsystem consisting of r cascaded integrators driven by input v, which then can have its dynamics redesigned by eg pole placement or LQ methods etc In normal form, the η -subsystem is not observable; it represents the remnant, and so must have stable dynamics for the approach to be useful This is investigated in the section entitled Zero Dynamics 4 Zero Dynamics 4 SISO Keeping y(t) ( zeroing the output ) requires that ξ (t) (see output equation) leading to ξ(t) Thus ξ r (t) in turn requires that v and so u must be chosen as u = u zd = A(Φ (, η)) B(Φ (, η)) The remnant or zero dynamics are then η = D(Φ (, η)) + E(Φ (, η))u zd }{{} F(η) (4) The stability of the zero (equilibrium) solution of this flow can be tested by any of the usual methods 5
6 5 Examples 5 Example Consider the 2nd order affine system ẋ = ( ) ( ) x2 + 4x 2 x 2 + u x }{{}}{{} f(x) g(x) y = x }{{} h(x) Computing the relative degree y = x ẏ = ẋ = x 2 ( + 4x 2 ) ÿ = ẋ 2 ( + 4x 2 ) + x 2 (8x ẋ ) = ( x + u)( + 4x 2 ) + 8x x 2 2( + 4x 2 ) = x ( + 4x 2 )( + 8x 2 }{{ 2) + } ( + 4x 2 }{{ ) } u L f 2 h(x,x 2 )=A(x,x 2 ) L gl f h(x,x 2 )=B(x,x 2 ) Therefore r = 2 Hence r = n and so there is no remnant Describing the Normal Form ξ = Φ (x) = h(x) = x ξ 2 = Φ 2 (x) = L f h(x) = x 2 ( + 4x 2 ) Solving to get the inverse transformation x = Φ (ξ) = ξ x 2 = Φ (ξ) = ξ 2 + 4ξ 2 The normal form is ξ = ξ 2 ξ 2 = A + Bu y = ξ 6
7 I/O Linearisation Feedback With u = A + v B the linearised system is = x ( + 4x 2 )( + 8x 2 2) + v + 4x 2 = x ( + 8x 2 v 2) + + 4x 2 ξ = ξ 2 ξ 2 = v y = ξ Stabilising feedback This double integrator system can be stabilised by where k, k 2 < Hence v = k ξ + k 2 ξ 2 + w = k x + k 2 x 2 ( + 4x 2 ) + w u = x ( + 8x 2 x 2) + k + k + 4x 2 2 x 2 + w + 4x 2 52 Example 2 Consider the 2nd order affine system (same state equation as Ex different output) ẋ = ( ) ( ) x2 + 4x 2 x 2 + u x }{{}}{{} f(x) g(x) y = x }{{} 2 h(x) but Computing the relative degree y = x ẏ = ẋ 2 = x }{{} + }{{} u L f h(x,x 2 )=A(x,x 2 ) L gh(x,x 2 )=B(x,x 2 ) 7
8 Therefore r = Hence r < n and so there is a remnant Describing the Normal Form ξ = Φ (x) = h(x) = x 2 η 2 = Φ 2 (x) η 2 = Φ 2 x (ẋ) = Φ 2 (f(x) + g(x)u) x = Φ 2 x f(x) + Φ 2 }{{} x g(x) u }{{} L f Φ 2 (x) L gφ 2 (x) From Lemma 2 it is known that Φ 2 can be chosen so that L g Φ 2 (x) =, but we calculate that L g Φ 2 (x) = Φ 2 x 2 and so Φ 2 (x) = C(x ), where C is a function of its argument Since η 2 = Φ 2 (x) = C(x ), C must also be invertible Solving to get the inverse transformation The normal form is (I/O Linearisation Feedback) Zero Dynamics For this example x = Φ (ξ) = C (η 2 ) x 2 = Φ (ξ) = ξ ξ = A + Bu = x + u η 2 = dc dx ( + 4x 2 )x 2 y = ξ y x 2 Hence the zero dynamics are (see Normal form above) η 2 = What does this imply about the behaviour of this control system? 8
9 6 Exact Feedback Linearisation The previous two examples demonstrate some features of geometric control which will be expanded on in this section: In Example, the system has relative degree equal to system order, no remnant, is transformable to a linear integrator-cascade system and is I/O stabilisable Example 2 has relative degree less than the system order, a remnant whose zero dynamics are marginally stable but which can be I/O stabilised Proposition There exists an output function y = h(x) for which the relative degree r equals the system order n is necessary and sufficient for the nonlinear system ẋ = f(x) + g(x)u to have a linearising transformation z = Φ(x) and a linearising feedback u = α(x)+β(x)v which transforms the system into a completely controllable linear system The output y may be real or synthetic Example gives a constructive proof method for sufficiency which may be generalised to n-th order system The proof of necessity needs the following two results Lemma 3 (Invariance of r under state transformation): Setting f(z) = [ Φ x f(x) ], g(z) = x=φ (z) [ ] Φ x g(x), h(z) = h(φ (z)), x=φ (z) then L f h(z) = h z f(z) = = [ ] h x [ ] h x f(x) x=φ (z) [ Φ x=φ (z) z ] [ ] Φ x f(x) x=φ (z) = [L f h(x)] x=φ (z) then iterated calculations lead to L f kh(z) = [L f kh(x)] x=φ (z) and L g L f kh(z) = [L g L f kh(x)] x=φ (z) got by rewriting A + Bu = v 9
10 Lemma 4 (Invariance of r under feedback): By induction, it may be shown that from which it follows that and provided that β(x e ) L (f+gα) kh(x) = L f kh(x), k =,,, r L gβ L (f+gα) kh(x) =, k =,,, r 2 L gβ L (f+gα) r h(x e ) To show the necessity of r = n, start by considering a linear system in Brunovsky canonical form 2 : Ã =, b = If the nonlinear system is feedback linearisable then it can be transformed to an CC linear system which in turn can be transformed into Brunovsky canonical form: ẋ = f(x) + g(x)u z=φ(x) u=α+βv ż = Az + bv z= T z v=kz+w z = Ã z + bw With this Brunovsky form system, associate the output y = (,,, ) z It is straightforward to show that r = n The invariance of r under transformation and feedback means that there is an appropriate output y = h(x) associated with the nonlinear system The problem of finding h(x) such that r = n comes down to finding a solution to the set of equations/inequation: L g L f k h(x) =, k =,,, n 2 L g L f n h(xe ) 2 This is a particular instance of ccf
11 This in turn is equivalent to the set of st order equations/inequation: L adf k g h(x) =, k =,,, n 2 (5) L adf n g h(x e ) (6) Necessary and sufficient conditions for the existence of a function h(x) satisfying these relations ( (5) and (6) ) is a consequence of Frobenius s theorem and are given by (i) (ii) F n = {g(x), ad f g(x),, ad f n 2g(x), ad f n g(x)} is linearly independent in a neighbourhood of x = x e This corresponds to complete controllability in linear systems F n = {g(x), ad f g(x),, ad f n 2g(x)} is involutive in a neighbourhood of x = x e This is the condition that guarantees the solution to the set of PDEs (Eq (5)) It automatically holds for linear systems If conditions (i) and (ii) hold, solving the relations (5) and (6) gives h(x) = y = Φ (x) = z, the first component of the state transformation The other components are given by (see normal form Eq( 7)) z i+ = L f z i, i =, 2,, n Example 3 Consider the 2nd-order affine system (the same state equation as Ex & Ex2 but with no output): ( ) ( ) x2 + 4x ẋ = 2 x 2 + u x }{{}}{{} f(x) g(x) We calculate ad f g = [f, g] = g f f(x) g(x) = x x ( 8x x 2 + 4x 2 ) ( ) = ( 4x 2 )
12 Hence F 2 = {g, ad f g} = {( ), ( )} 4x 2 F 2 is linearly independent by inspection and F is involutive 3 So, Solve L g z =, L adf gz ie L adf gz = ( z x z x 2 ) ( 4x 2 L g z = ( z z ) ( ) x x 2 ) ( z = x = z x 2 = (7) ) ( 4x 2 ) (8) To satisfy InEQ (8), choose z = (9) x Then a solution of Eq (7) and Eq (9) is z = x (2) The rest of the solution then reduces to Example What happens if conditions (i) and (ii) are not satisfied? Then it is still possible to obtain by means of state feedback a partially linear system As a matter of fact since r < n using the linearising control u = α(x) + β(x)v the system becomes ξ = ξ 2 ξ 2 = ξ 3 ξ r = ξ r ξ r = v (2) η r+ = q r+ (z)) η r+2 = q r+2 (z)) η n = q n (z)) y = ξ 3 Every singleton set of vector fields is involutive 2
13 where q i () = D i (Φ ()), i = r +,, n (see Eq ()) Figure : Normal Form 3
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.
Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control
More informationOutput Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)
Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:
More informationNonlinear Control Lecture 9: Feedback Linearization
Nonlinear Control Lecture 9: Feedback Linearization Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 9 1/75
More informationNonlinear Control Systems
Nonlinear Control Systems António Pedro Aguiar pedro@isr.ist.utl.pt 7. Feedback Linearization IST-DEEC PhD Course http://users.isr.ist.utl.pt/%7epedro/ncs1/ 1 1 Feedback Linearization Given a nonlinear
More informationFeedback Linearization
Feedback Linearization Peter Al Hokayem and Eduardo Gallestey May 14, 2015 1 Introduction Consider a class o single-input-single-output (SISO) nonlinear systems o the orm ẋ = (x) + g(x)u (1) y = h(x) (2)
More informationDisturbance Decoupling Problem
DISC Systems and Control Theory of Nonlinear Systems, 21 1 Disturbance Decoupling Problem Lecture 4 Nonlinear Dynamical Control Systems, Chapter 7 The disturbance decoupling problem is a typical example
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University
More informationStabilization and Passivity-Based Control
DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive
More information1 The Observability Canonical Form
NONLINEAR OBSERVERS AND SEPARATION PRINCIPLE 1 The Observability Canonical Form In this Chapter we discuss the design of observers for nonlinear systems modelled by equations of the form ẋ = f(x, u) (1)
More informationA nemlineáris rendszer- és irányításelmélet alapjai Relative degree and Zero dynamics (Lie-derivatives)
A nemlineáris rendszer- és irányításelmélet alapjai Relative degree and Zero dynamics (Lie-derivatives) Hangos Katalin BME Analízis Tanszék Rendszer- és Irányításelméleti Kutató Laboratórium MTA Számítástechnikai
More information1 Relative degree and local normal forms
THE ZERO DYNAMICS OF A NONLINEAR SYSTEM 1 Relative degree and local normal orms The purpose o this Section is to show how single-input single-output nonlinear systems can be locally given, by means o a
More informationChap. 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 informationZeros and zero dynamics
CHAPTER 4 Zeros and zero dynamics 41 Zero dynamics for SISO systems Consider a linear system defined by a strictly proper scalar transfer function that does not have any common zero and pole: g(s) =α p(s)
More informationRobust feedback linearization
Robust eedback linearization Hervé Guillard Henri Bourlès Laboratoire d Automatique des Arts et Métiers CNAM/ENSAM 21 rue Pinel 75013 Paris France {herveguillardhenribourles}@parisensamr Keywords: Nonlinear
More informationNonlinear Control Lecture # 1 Introduction. Nonlinear Control
Nonlinear Control Lecture # 1 Introduction Nonlinear State Model ẋ 1 = f 1 (t,x 1,...,x n,u 1,...,u m ) ẋ 2 = f 2 (t,x 1,...,x n,u 1,...,u m ).. ẋ n = f n (t,x 1,...,x n,u 1,...,u m ) ẋ i denotes the derivative
More informationLecture : Feedback Linearization
ecture : Feedbac inearization Niola Misovic, dipl ing and Pro Zoran Vuic June 29 Summary: This document ollows the lectures on eedbac linearization tought at the University o Zagreb, Faculty o Electrical
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationMathematics for Control Theory
Mathematics for Control Theory Geometric Concepts in Control Involutivity and Frobenius Theorem Exact Linearization Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials
More informationExact Linearization Of Stochastic Dynamical Systems By State Space Coordinate Transformation And Feedback I g-linearization
Applied Mathematics E-Notes, 3(2003), 99-106 c ISSN 1607-2510 Available free at mirror sites of http://www.math.nthu.edu.tw/ amen/ Exact Linearization Of Stochastic Dynamical Systems By State Space Coordinate
More informationEquilibrium points: continuous-time systems
Capitolo 0 INTRODUCTION 81 Equilibrium points: continuous-time systems Let us consider the following continuous-time linear system ẋ(t) Ax(t)+Bu(t) y(t) Cx(t)+Du(t) The equilibrium points x 0 of the system
More informationFeedback linearization control of systems with singularities: a ball-beam revisit
Chapter 1 Feedback linearization control o systems with singularities: a ball-beam revisit Fu Zhang The Mathworks, Inc zhang@mathworks.com Benito Fernndez-Rodriguez Department o Mechanical Engineering,
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationMinimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality
Minimum-Phase Property of Nonlinear Systems in Terms of a Dissipation Inequality Christian Ebenbauer Institute for Systems Theory in Engineering, University of Stuttgart, 70550 Stuttgart, Germany ce@ist.uni-stuttgart.de
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationSolution of Additional Exercises for Chapter 4
1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the
More informationLECTURE 7, WEDNESDAY
LECTURE 7, WEDNESDAY 25.02.04 FRANZ LEMMERMEYER 1. Singular Weierstrass Curves Consider cubic curves in Weierstraß form (1) E : y 2 + a 1 xy + a 3 y = x 3 + a 2 x 2 + a 4 x + a 6, the coefficients a i
More informationINVERSION IN INDIRECT OPTIMAL CONTROL
INVERSION IN INDIRECT OPTIMAL CONTROL François Chaplais, Nicolas Petit Centre Automatique et Systèmes, École Nationale Supérieure des Mines de Paris, 35, rue Saint-Honoré 7735 Fontainebleau Cedex, France,
More informationCDS Solutions to the Midterm Exam
CDS 22 - Solutions to the Midterm Exam Instructor: Danielle C. Tarraf November 6, 27 Problem (a) Recall that the H norm of a transfer function is time-delay invariant. Hence: ( ) Ĝ(s) = s + a = sup /2
More informationNetwork Analysis of Biochemical Reactions in Complex Environments
1 Introduction 1 Network Analysis of Biochemical Reactions in Complex Environments Elias August 1 and Mauricio Barahona, Department of Bioengineering, Imperial College London, South Kensington Campus,
More informationLecture 2: Controllability of nonlinear systems
DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 57 Outline 1 Lecture 13: Linear System - Stability
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationKalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and
Kalman Decomposition Controllable / uncontrollable decomposition Suppose that the controllability matrix C R n n of a system has rank n 1
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationStability 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 informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georgia Tech PHYS 612 Mathematical Methods of Physics I Instructor: Predrag Cvitanović Fall semester 2012 Homework Set #5 due October 2, 2012 == show all your work for maximum credit, == put labels, title,
More informationMS 3011 Exercises. December 11, 2013
MS 3011 Exercises December 11, 2013 The exercises are divided into (A) easy (B) medium and (C) hard. If you are particularly interested I also have some projects at the end which will deepen your understanding
More informationFirst Integrals/Invariants & Symmetries for Autonomous Difference Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 3, 1253 1260 First Integrals/Invariants & Symmetries for Autonomous Difference Equations Leon ARRIOLA Department of Mathematical
More informationOptimal Control. Lecture 18. Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen. March 29, Ref: Bryson & Ho Chapter 4.
Optimal Control Lecture 18 Hamilton-Jacobi-Bellman Equation, Cont. John T. Wen Ref: Bryson & Ho Chapter 4. March 29, 2004 Outline Hamilton-Jacobi-Bellman (HJB) Equation Iterative solution of HJB Equation
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation
High-Gain Observers in Nonlinear Feedback Control Lecture # 3 Regulation High-Gain ObserversinNonlinear Feedback ControlLecture # 3Regulation p. 1/5 Internal Model Principle d r Servo- Stabilizing u y
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
More informationBalanced realization and model order reduction for nonlinear systems based on singular value analysis
Balanced realization and model order reduction for nonlinear systems based on singular value analysis Kenji Fujimoto a, and Jacquelien M. A. Scherpen b a Department of Mechanical Science and Engineering
More informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationEE363 homework 8 solutions
EE363 Prof. S. Boyd EE363 homework 8 solutions 1. Lyapunov condition for passivity. The system described by ẋ = f(x, u), y = g(x), x() =, with u(t), y(t) R m, is said to be passive if t u(τ) T y(τ) dτ
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationMATH 220: MIDTERM OCTOBER 29, 2015
MATH 22: MIDTERM OCTOBER 29, 25 This is a closed book, closed notes, no electronic devices exam. There are 5 problems. Solve Problems -3 and one of Problems 4 and 5. Write your solutions to problems and
More informationPart II. Dynamical Systems. Year
Part II Year 2017 2016 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2017 34 Paper 1, Section II 30A Consider the dynamical system where β > 1 is a constant. ẋ = x + x 3 + βxy 2, ẏ = y + βx 2
More informationAn homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted pendulum
9 American Control Conference Hyatt Regency Riverfront, St. Louis, MO, USA June -, 9 FrA.5 An homotopy method for exact tracking of nonlinear nonminimum phase systems: the example of the spherical inverted
More informationMATH 819 FALL We considered solutions of this equation on the domain Ū, where
MATH 89 FALL. The D linear wave equation weak solutions We have considered the initial value problem for the wave equation in one space dimension: (a) (b) (c) u tt u xx = f(x, t) u(x, ) = g(x), u t (x,
More informationGreen s Functions and Distributions
CHAPTER 9 Green s Functions and Distributions 9.1. Boundary Value Problems We would like to study, and solve if possible, boundary value problems such as the following: (1.1) u = f in U u = g on U, where
More informationIntro. Computer Control Systems: F8
Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2
More informationLecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.
Lecture 8 Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture
More informationYURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL
Journal of Comput. & Applied Mathematics 139(2001), 197 213 DIRECT APPROACH TO CALCULUS OF VARIATIONS VIA NEWTON-RAPHSON METHOD YURI LEVIN, MIKHAIL NEDIAK, AND ADI BEN-ISRAEL Abstract. Consider m functions
More informationHigh-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle
High-Gain Observers in Nonlinear Feedback Control Lecture # 2 Separation Principle High-Gain ObserversinNonlinear Feedback ControlLecture # 2Separation Principle p. 1/4 The Class of Systems ẋ = Ax + Bφ(x,
More informationLyapunov Stability Theory
Lyapunov Stability Theory Peter Al Hokayem and Eduardo Gallestey March 16, 2015 1 Introduction In this lecture we consider the stability of equilibrium points of autonomous nonlinear systems, both in continuous
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationOn the Stability of the Best Reply Map for Noncooperative Differential Games
On the Stability of the Best Reply Map for Noncooperative Differential Games Alberto Bressan and Zipeng Wang Department of Mathematics, Penn State University, University Park, PA, 68, USA DPMMS, University
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationLecture 4 and 5 Controllability and Observability: Kalman decompositions
1 Lecture 4 and 5 Controllability and Observability: Kalman decompositions Spring 2013 - EE 194, Advanced Control (Prof. Khan) January 30 (Wed.) and Feb. 04 (Mon.), 2013 I. OBSERVABILITY OF DT LTI SYSTEMS
More information14 Divergence, flux, Laplacian
Tel Aviv University, 204/5 Analysis-III,IV 240 4 Divergence, flux, Laplacian 4a What is the problem................ 240 4b Integral of derivative (again)........... 242 4c Divergence and flux.................
More informationPoincaré Map, Floquet Theory, and Stability of Periodic Orbits
Poincaré Map, Floquet Theory, and Stability of Periodic Orbits CDS140A Lecturer: W.S. Koon Fall, 2006 1 Poincaré Maps Definition (Poincaré Map): Consider ẋ = f(x) with periodic solution x(t). Construct
More informationTHE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS
THE HARTMAN-GROBMAN THEOREM AND THE EQUIVALENCE OF LINEAR SYSTEMS GUILLAUME LAJOIE Contents 1. Introduction 2 2. The Hartman-Grobman Theorem 2 2.1. Preliminaries 2 2.2. The discrete-time Case 4 2.3. The
More informationImplicit Functions, Curves and Surfaces
Chapter 11 Implicit Functions, Curves and Surfaces 11.1 Implicit Function Theorem Motivation. In many problems, objects or quantities of interest can only be described indirectly or implicitly. It is then
More informationSTATE OBSERVERS FOR NONLINEAR SYSTEMS WITH SMOOTH/BOUNDED INPUT 1
K Y B E R N E T I K A V O L U M E 3 5 ( 1 9 9 9 ), N U M B E R 4, P A G E S 3 9 3 4 1 3 STATE OBSERVERS FOR NONLINEAR SYSTEMS WITH SMOOTH/BOUNDED INPUT 1 Alfredo Germani and Costanzo Manes It is known
More informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More informationPSEUDO-HAMILTONIAN REALIZATION AND ITS APPLICATION
COMMUNICATIONS IN INFORMATION AND SYSTEMS c 2002 International Press Vol. 2, No. 2, pp. 91-120, December 2002 001 PSEUDO-HAMILTONIAN REALIZATION AND ITS APPLICATION DAIZHAN CHENG, TIELONG SHEN, AND T.
More informationB5.6 Nonlinear Systems
B5.6 Nonlinear Systems 4. Bifurcations Alain Goriely 2018 Mathematical Institute, University of Oxford Table of contents 1. Local bifurcations for vector fields 1.1 The problem 1.2 The extended centre
More informationFIXED POINT ITERATIONS
FIXED POINT ITERATIONS MARKUS GRASMAIR 1. Fixed Point Iteration for Non-linear Equations Our goal is the solution of an equation (1) F (x) = 0, where F : R n R n is a continuous vector valued mapping in
More informationMathematical Systems Theory: Advanced Course Exercise Session 5. 1 Accessibility of a nonlinear system
Mathematical Systems Theory: dvanced Course Exercise Session 5 1 ccessibility of a nonlinear system Consider an affine nonlinear control system: [ ẋ = f(x)+g(x)u, x() = x, G(x) = g 1 (x) g m (x) ], where
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationSOME COMMENTS ON OBSERVABILITY NORMAL FORM AND STEP-BY-STEP SLIDING MODE OBSERVER
SOME COMMENTS ON OBSERVABILITY NORMAL FORM AND STEP-BY-STEP SLIDING MODE OBSERVER G ZHENG 1, L BOUTAT-BADDAS 2, JP BARBOT 1 AND D BOUTAT 3 1 EQUIPE COMMANDE DES SYSTÈMES (ECS), ENSEA, 6 AV DU PONCEAU,
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationThe disturbance decoupling problem (DDP)
CHAPTER 3 The disturbance decoupling problem (DDP) Consider the system 3.1. Geometric formulation { ẋ = Ax + Bu + Ew y = Cx. Problem 3.1 (Disturbance decoupling). Find a state feedback u = Fx+ v such that
More informationCOMMUTATIVE PRESEMIFIELDS AND SEMIFIELDS
COMMUTATIVE PRESEMIFIELDS AND SEMIFIELDS ROBERT S. COULTER AND MARIE HENDERSON Abstract. Strong conditions are derived for when two commutative presemifields are isotopic. It is then shown that any commutative
More informationUC Berkeley Department of Electrical Engineering and Computer Science. EECS 227A Nonlinear and Convex Optimization. Solutions 5 Fall 2009
UC Berkeley Department of Electrical Engineering and Computer Science EECS 227A Nonlinear and Convex Optimization Solutions 5 Fall 2009 Reading: Boyd and Vandenberghe, Chapter 5 Solution 5.1 Note that
More informationSome solutions of the written exam of January 27th, 2014
TEORIA DEI SISTEMI Systems Theory) Prof. C. Manes, Prof. A. Germani Some solutions of the written exam of January 7th, 0 Problem. Consider a feedback control system with unit feedback gain, with the following
More information8. Geometric problems
8. Geometric problems Convex Optimization Boyd & Vandenberghe extremal volume ellipsoids centering classification placement and facility location 8 Minimum volume ellipsoid around a set Löwner-John ellipsoid
More informationEach is equal to CP 1 minus one point, which is the origin of the other: (C =) U 1 = CP 1 the line λ (1, 0) U 0
Algebraic Curves/Fall 2015 Aaron Bertram 1. Introduction. What is a complex curve? (Geometry) It s a Riemann surface, that is, a compact oriented twodimensional real manifold Σ with a complex structure.
More informationSolutions. Problems of Chapter 1
Solutions Problems of Chapter 1 1.1 A real square matrix A IR n n is invertible if and only if its determinant is non zero. The determinant of A is a polynomial in the entries a ij of A, whose set of zeros
More informationSome recent results on controllability of coupled parabolic systems: Towards a Kalman condition
Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition F. Ammar Khodja Clermont-Ferrand, June 2011 GOAL: 1 Show the important differences between scalar and non
More informationNonlinear Control. Nonlinear Control Lecture # 25 State Feedback Stabilization
Nonlinear Control Lecture # 25 State Feedback Stabilization Backstepping η = f a (η)+g a (η)ξ ξ = f b (η,ξ)+g b (η,ξ)u, g b 0, η R n, ξ, u R Stabilize the origin using state feedback View ξ as virtual
More informationA nonlinear equation is any equation of the form. f(x) = 0. A nonlinear equation can have any number of solutions (finite, countable, uncountable)
Nonlinear equations Definition A nonlinear equation is any equation of the form where f is a nonlinear function. Nonlinear equations x 2 + x + 1 = 0 (f : R R) f(x) = 0 (x cos y, 2y sin x) = (0, 0) (f :
More informationEquivalence of dynamical systems by bisimulation
Equivalence of dynamical systems by bisimulation Arjan van der Schaft Department of Applied Mathematics, University of Twente P.O. Box 217, 75 AE Enschede, The Netherlands Phone +31-53-4893449, Fax +31-53-48938
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationElementary ODE Review
Elementary ODE Review First Order ODEs First Order Equations Ordinary differential equations of the fm y F(x, y) () are called first der dinary differential equations. There are a variety of techniques
More informationFIXED POINT ITERATION
FIXED POINT ITERATION The idea of the fixed point iteration methods is to first reformulate a equation to an equivalent fixed point problem: f (x) = 0 x = g(x) and then to use the iteration: with an initial
More informationMulti-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 informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationNonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control
Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded
More informationHence, (f(x) f(x 0 )) 2 + (g(x) g(x 0 )) 2 < ɛ
Matthew Straughn Math 402 Homework 5 Homework 5 (p. 429) 13.3.5, 13.3.6 (p. 432) 13.4.1, 13.4.2, 13.4.7*, 13.4.9 (p. 448-449) 14.2.1, 14.2.2 Exercise 13.3.5. Let (X, d X ) be a metric space, and let f
More informationClassification of Second order PDEs
Chapter 3 Classification of Second order PDEs This chapter deals with a finer classification of second order quasilinear PDEs, as compared to Chapter 1. In Section 3.1, we present a formal procedure to
More informationTHE METHOD OF CHARACTERISTICS FOR SYSTEMS AND APPLICATION TO FULLY NONLINEAR EQUATIONS
THE METHOD OF CHARACTERISTICS FOR SYSTEMS AND APPLICATION TO FULLY NONLINEAR EQUATIONS ARICK SHAO While the method of characteristics is used for solving general first-order partial differential equations
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationNonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability
Nonlinear Control Lecture # 6 Passivity and Input-Output Stability Passivity: Memoryless Functions y y y u u u (a) (b) (c) Passive Passive Not passive y = h(t,u), h [0, ] Vector case: y = h(t,u), h T =
More informationVolume 30, Issue 3. A note on Kalman filter approach to solution of rational expectations models
Volume 30, Issue 3 A note on Kalman filter approach to solution of rational expectations models Marco Maria Sorge BGSE, University of Bonn Abstract In this note, a class of nonlinear dynamic models under
More informationControllability, Observability & Local Decompositions
ontrollability, Observability & Local Decompositions Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Lie Bracket Distributions ontrollability ontrollability Distributions
More information