Geometric Control Theory

Similar documents
The disturbance decoupling problem (DDP)

Contents. Preface for the Instructor. Preface for the Student. xvii. Acknowledgments. 1 Vector Spaces 1 1.A R n and C n 2

Contents. 1 State-Space Linear Systems 5. 2 Linearization Causality, Time Invariance, and Linearity 31

MATH 304 Linear Algebra Lecture 20: The Gram-Schmidt process (continued). Eigenvalues and eigenvectors.

Rank & nullity. Defn. Let T : V W be linear. We define the rank of T to be rank T = dim im T & the nullity of T to be nullt = dim ker T.

Lecture 4 and 5 Controllability and Observability: Kalman decompositions

DIAGONALIZATION. In order to see the implications of this definition, let us consider the following example Example 1. Consider the matrix

EL2520 Control Theory and Practice

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA

Matrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein

1 Invariant subspaces

MATH 221, Spring Homework 10 Solutions

Zeros and zero dynamics

No books, no notes, no calculators. You must show work, unless the question is a true/false, yes/no, or fill-in-the-blank question.

Math 1553, Introduction to Linear Algebra

Review Notes for Linear Algebra True or False Last Updated: February 22, 2010

MCE693/793: Analysis and Control of Nonlinear Systems

1 Linearity and Linear Systems

Controllability, Observability & Local Decompositions

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #1. July 11, 2013 Solutions

MAT 242 CHAPTER 4: SUBSPACES OF R n

2. Every linear system with the same number of equations as unknowns has a unique solution.

CSL361 Problem set 4: Basic linear algebra

Worksheet for Lecture 15 (due October 23) Section 4.3 Linearly Independent Sets; Bases

series. Utilize the methods of calculus to solve applied problems that require computational or algebraic techniques..

Applied Linear Algebra in Geoscience Using MATLAB

Section 6.2: THE KERNEL AND RANGE OF A LINEAR TRANSFORMATIONS

MATH Linear Algebra

Robust Control 2 Controllability, Observability & Transfer Functions

5.) For each of the given sets of vectors, determine whether or not the set spans R 3. Give reasons for your answers.

Linear Algebra- Final Exam Review

INTRODUCTION TO LIE ALGEBRAS. LECTURE 2.

1. Let m 1 and n 1 be two natural numbers such that m > n. Which of the following is/are true?

Chapter 3. Vector spaces

Review of Some Concepts from Linear Algebra: Part 2

Disturbance Decoupling Problem

1. What is the determinant of the following matrix? a 1 a 2 4a 3 2a 2 b 1 b 2 4b 3 2b c 1. = 4, then det

1. Linear systems of equations. Chapters 7-8: Linear Algebra. Solution(s) of a linear system of equations (continued)

Matrices and Matrix Algebra.

Notes on singular value decomposition for Math 54. Recall that if A is a symmetric n n matrix, then A has real eigenvalues A = P DP 1 A = P DP T.

ẋ 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)

Lecture: Linear algebra. 4. Solutions of linear equation systems The fundamental theorem of linear algebra

Midterm solutions. (50 points) 2 (10 points) 3 (10 points) 4 (10 points) 5 (10 points)

Econ Slides from Lecture 7

Jacobians of Matrix Transformations and Functions of Matrix Argument

Lecture 9: Vector Algebra

Differential Geometry, Lie Groups, and Symmetric Spaces

Lecture 10 - Eigenvalues problem

MATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)

NORMS ON SPACE OF MATRICES

COMP 558 lecture 18 Nov. 15, 2010

Module 09 From s-domain to time-domain From ODEs, TFs to State-Space Modern Control

3.2 Frobenius Theorem

Worksheet for Lecture 15 (due October 23) Section 4.3 Linearly Independent Sets; Bases

AUTOMATIC CONTROL. Andrea M. Zanchettin, PhD Spring Semester, Introduction to Automatic Control & Linear systems (time domain)

Numerical Linear Algebra Homework Assignment - Week 2

NOTES ON BILINEAR FORMS

What is A + B? What is A B? What is AB? What is BA? What is A 2? and B = QUESTION 2. What is the reduced row echelon matrix of A =

Linear Algebra Review. Vectors

Linear Systems. Class 27. c 2008 Ron Buckmire. TITLE Projection Matrices and Orthogonal Diagonalization CURRENT READING Poole 5.4

Lecture 13: Orthogonal projections and least squares (Section ) Thang Huynh, UC San Diego 2/9/2018

homogeneous 71 hyperplane 10 hyperplane 34 hyperplane 69 identity map 171 identity map 186 identity map 206 identity matrix 110 identity matrix 45

The Essentials of Linear State-Space Systems

Lecture 11: Finish Gaussian elimination and applications; intro to eigenvalues and eigenvectors (1)

CAAM 336 DIFFERENTIAL EQUATIONS IN SCI AND ENG Examination 1

Linear Combination. v = a 1 v 1 + a 2 v a k v k

(a) II and III (b) I (c) I and III (d) I and II and III (e) None are true.

Eigenvalues, Eigenvectors. Eigenvalues and eigenvector will be fundamentally related to the nature of the solutions of state space systems.

Family Feud Review. Linear Algebra. October 22, 2013

Definition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition

LINES IN P 3. Date: December 12,

Lecture 1: Systems of linear equations and their solutions

Input Tracking and Output Fusion for Linear Systems

MATH 583A REVIEW SESSION #1

A proof of the Jordan normal form theorem

MATH 31 - ADDITIONAL PRACTICE PROBLEMS FOR FINAL

EC Control Engineering Quiz II IIT Madras

Linear Algebra. Grinshpan

Background Mathematics (2/2) 1. David Barber

IMPORTANT DEFINITIONS AND THEOREMS REFERENCE SHEET

EK102 Linear Algebra PRACTICE PROBLEMS for Final Exam Spring 2016

MTH 464: Computational Linear Algebra

Advanced Linear Algebra Math 4377 / 6308 (Spring 2015) March 5, 2015

LINEAR MMSE ESTIMATION

Elements of linear algebra

Vector Spaces and Linear Transformations

YORK UNIVERSITY. Faculty of Science Department of Mathematics and Statistics MATH M Test #2 Solutions

Controllability, Observability, Full State Feedback, Observer Based Control

Lecture 11 The Radical and Semisimple Lie Algebras

EE731 Lecture Notes: Matrix Computations for Signal Processing

Worksheet for Lecture 25 Section 6.4 Gram-Schmidt Process

Chapter 2 Subspaces of R n and Their Dimensions

Statistical and Adaptive Signal Processing

Design of Nonlinear Control Systems with the Highest Derivative in Feedback

Mathematics Department

VECTOR SPACES & SUBSPACES

1 Last time: least-squares problems

MATH SOLUTIONS TO PRACTICE MIDTERM LECTURE 1, SUMMER Given vector spaces V and W, V W is the vector space given by

Kalman Decomposition B 2. z = T 1 x, where C = ( C. z + u (7) T 1, and. where B = T, and

Math 123, Week 5: Linear Independence, Basis, and Matrix Spaces. Section 1: Linear Independence

Transcription:

1 Geometric Control Theory Lecture notes by Xiaoming Hu and Anders Lindquist in collaboration with Jorge Mari and Janne Sand 2012 Optimization and Systems Theory Royal institute of technology SE-100 44 Stockholm, Sweden 1 This work is partially based on the earlier lecture notes by Lindquist, Mari and Sand.

Contents Preface Notations and definitions vii ix Chapter 1. Introduction 1 1.1. Linear and nonlinear systems 1 1.2. The geometric approach 4 Chapter 2. Invariant and controlled invariant subspaces 9 2.1. Invariant subspaces 9 2.2. Controlled invariant subspaces 10 2.3. Reachability subspaces 12 2.4. Maximal reachability subspaces 15 2.5. Reachability under state constraints 19 Chapter 3. The disturbance decoupling problem (DDP) 21 3.1. Geometric formulation 21 3.2. Computing V 23 3.3. Disturbance decoupling, observability, and zeros 25 3.4. Disturbance decoupling with eigenvalue assignment 26 3.5. Solution of the DDPP 27 3.6. Is it necessary that Im E R? 29 Chapter 4. Zeros and zero dynamics 31 4.1. Zero dynamics for SISO systems 31 4.2. Zero dynamics of MIMO systems 33 4.3. Zeros and system inversion 36 4.4. An illustration of zero dynamics: high gain control 1 38 Chapter 5. Noninteracting control and tracking 41 5.1. Noninteracting control 41 5.2. Tracking with stability 44 Chapter 6. Input-output behavior 45 6.1. State observation 45 6.2. Output tracking input 46 6.3. The partial stochastic realization problem 51 Chapter 7. Output regulation and internal model principle 55 iii

iv CONTENTS 7.1. Full information output regulation 57 7.2. Error feedback output regulation 58 7.3. Output regulation and zero dynamics 59 Chapter 8. Nonlinear systems 61 8.1. Introduction 61 8.2. Controllability 62 8.3. Stability of nonlinear systems 66 8.4. Steady state response and center manifold 69 8.5. Center Manifold Theory 71 8.6. Zero dynamics and its applications 73 8.7. Disturbance decoupling problem (DDP) 76 8.8. Output regulation 76 8.9. Exact linearization via feedback 77 Chapter 9. Multi-Robotic Systems 79 9.1. Consensus problem 79 9.2. Formation control 82 9.3. Formation control with limited sensor information 85 Appendix A. Numerical aspects in geometric control 89 A.1. Introduction 89 A.2. Representation of linear model 90 A.3. The singular value decomposition 92 A.4. A list of useful MATLAB commands 95 Appendix. Bibliography 101 Appendix. Index 103

Preface Various versions of the notes have been used for the course Geometric Control Theory given at the Royal Institute of Technology (KTH), and its predecessor Advanced Systems Theory. The first version was written by Anders Lindquist and Janne Sand, and was later revised and extended by Jorge Mari. A major revision and addition was done by Xiaoming Hu in 2002. I would like to express my gratitude to Dr Ryozo Nagamune and Docent Ulf Jönsson for their careful reading and constructive comments on the 2002 version of the notes. Some minor changes and updates have been made every year since 2002. In 2006, some new material was added and the title was changed to Geometric Control Theory ; In 2012 Chapter 9 was rewritten. Xiaoming Hu October, 2012 vii

Notations and definitions For easy reference we recall here some notations and definitions. Assume X and Y are finite dimensional vector spaces over the real field. Let n be the dimension of X, which is then isomorphic to R n. (1) Image space. GivenamapA : X Y,theimage space of A is defined as Im A := {y Y : y = Ax, for some x X}. It is a subspace of Y. We shall often use the same symbol for the map A as for its matrix representation once bases in X and Y have been chosen. (2) Linear span. Given a vector space V over the field R,letv 1,,v m V. The span of these vectors is span{v 1,,v m } = {α 1 v 1 + + α m v m α 1,,α m R}. (3) Null space. GivenamapA : X Y,thenull space,orkernel, of A is defined as ker A := {x X : Ax =0}. It is a subspace of X. (4) Preimage. LetW be any set of Y. Thepre-image of W under the map A is A 1 W := {x X : Ax W}. Observe that A need not be invertible, so beware of the distinction between preimage and inverse. (5) A-invariant subspace. A subspace V of X is A-invariant if AV V. (6) Reachable subspace. Given the pair of conformable matrices A n n and B n k,thereachable subspace of (A, B), denoted by A Im B, is defined as A Im B := Im Γ, where Γ is the reachability matrix [BAB... A n 1 B]. This is an n nk matrix. The reachable subspace is A-invariant. (7) A Im E. The minimal A-invariant subspace that contains the subspace Im E is denoted by A Im E. ix

x NOTATIONS AND DEFINITIONS (8) Reachability subspace. Given a matrix pair (A, B), a subspace R is called a reachability subspace if there are matrices F and G such that R = A + BF Im BG. (9) Unobservable subspace. Given the matrix pair (C, A), the unobservable subspace is defined as ker Ω, where Ω is the observability matrix C CA. CA n 1 The unobservable subspace is A-invariant. (10) The map A is injective if Ax 1 = Ax 2 implies x 1 = x 2. (11) If A : X Yand V X,therestriction of A to V is denoted by A V. (12) Given the two finite-dimensional subspaces V and W, thevector sum is defined as V + W := {v + w : v V,w W}. If V and W are linearly independent we write the sum as V W. Note that in these notes does not denote orthogonal vector sum! (13) Hypersurface. Suppose N is an open set in R n.thesetm is defined as M = {x N : λ i (x) =0,i =1,...,n m} where λ i are smooth functions. λ 1 x If rank. = n m x M, thenm is a (hyper)surface λ n m x of dimension m. (14) Lie derivative. In local coordinates, Lie derivative is represented by n λ L f λ := f i, x i i=1 where f is a vector and λ is a scalar function. (15) Lie bracket. Lie bracket of the two vector fields f and g is defined according to the rule: [f,g](λ) :=L f L g λ L g L f λ.

NOTATIONS AND DEFINITIONS xi In local coordinates the expression of [f,g] is given as g x f f x g.