Chapter 6 Controllability and Obervability

Similar documents
Linear System Theory

Chap 4. State-Space Solutions and

Control Systems Design

Module 03 Linear Systems Theory: Necessary Background

MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011

Robust Control 2 Controllability, Observability & Transfer Functions

Module 07 Controllability and Controller Design of Dynamical LTI Systems

1 Continuous-time Systems

Chap 3. Linear Algebra

Observability and Constructability

Module 08 Observability and State Estimator Design of Dynamical LTI Systems

NOTES ON LINEAR ODES

Properties of Matrices and Operations on Matrices

16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1

Observability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)

Controllability, Observability, Full State Feedback, Observer Based Control

6.241 Dynamic Systems and Control

Jordan Canonical Form

Chapter 7. Canonical Forms. 7.1 Eigenvalues and Eigenvectors

SYSTEMTEORI - ÖVNING 1. In this exercise, we will learn how to solve the following linear differential equation:

1 Introduction. 2 Determining what the J i blocks look like. December 6, 2006

CME 345: MODEL REDUCTION

Graded Questions on Matrices. 1. Reduce the following matrix to its normal form & hence find its rank Where

Reachability, Observability and Minimality for a Class of 2D Continuous-Discrete Systems

ECEN 605 LINEAR SYSTEMS. Lecture 7 Solution of State Equations 1/77

Solving Homogeneous Systems with Sub-matrices

Chap. 3. Controlled Systems, Controllability

4F3 - Predictive Control

Lecture 4 and 5 Controllability and Observability: Kalman decompositions

State space transformations

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA

Ir O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )

Review of Controllability Results of Dynamical System

LESSON 35: EIGENVALUES AND EIGENVECTORS APRIL 21, (1) We might also write v as v. Both notations refer to a vector.

Linear dynamical systems with inputs & outputs

Eigenvalues and Eigenvectors

Chapter 1. Matrix Algebra

Lecture 9: SVD, Low Rank Approximation

Nonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México

Lecture 21: The decomposition theorem into generalized eigenspaces; multiplicity of eigenvalues and upper-triangular matrices (1)

Chapter III. Stability of Linear Systems

k is a product of elementary matrices.

ME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms

Review of Linear Algebra

Elementary linear algebra

Definition (T -invariant subspace) Example. Example

Lecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore

The matrix will only be consistent if the last entry of row three is 0, meaning 2b 3 + b 2 b 1 = 0.

Lecture 11: Diagonalization

Observability. It was the property in Lyapunov stability which allowed us to resolve that

1 Similarity transform 2. 2 Controllability The PBH test for controllability Observability The PBH test for observability...

Control Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft

Math 121 Homework 4: Notes on Selected Problems

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

CHAPTER 10 Shape Preserving Properties of B-splines

Systems and Control Theory Lecture Notes. Laura Giarré

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

Control Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft

Theorem 1. ẋ = Ax is globally exponentially stable (GES) iff A is Hurwitz (i.e., max(re(σ(a))) < 0).

6 Linear Equation. 6.1 Equation with constant coefficients

A matrix over a field F is a rectangular array of elements from F. The symbol

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

1. The Transition Matrix (Hint: Recall that the solution to the linear equation ẋ = Ax + Bu is

Control Systems. Linear Algebra topics. L. Lanari

Linear Algebra. Solving SLEs with Matlab. Matrix Inversion. Solving SLE s by Matlab - Inverse. Solving Simultaneous Linear Equations in MATLAB

MATHEMATICS 217 NOTES

6.241 Dynamic Systems and Control

Linear Equations in Linear Algebra

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

Fundamentals of Engineering Analysis (650163)

Optimal Sensor and Actuator Location for Descriptor Systems using Generalized Gramians and Balanced Realizations

Solution of Linear State-space Systems

Criterions on periodic feedback stabilization for some evolution equations

Lecture 2 and 3: Controllability of DT-LTI systems

Math 240 Calculus III

MAT Linear Algebra Collection of sample exams

Jordan Canonical Form

Canonical lossless state-space systems: staircase forms and the Schur algorithm

Putzer s Algorithm. Norman Lebovitz. September 8, 2016

Module 9: State Feedback Control Design Lecture Note 1

(f + g)(s) = f(s) + g(s) for f, g V, s S (cf)(s) = cf(s) for c F, f V, s S

Lecture Notes for Math 524

SUPPLEMENT TO CHAPTERS VII/VIII

MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction

1 Some Facts on Symmetric Matrices

Introduction to Modern Control MT 2016

EEE582 Homework Problems

ME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ

Lecture 22: Jordan canonical form of upper-triangular matrices (1)

Full State Feedback for State Space Approach

33A Linear Algebra and Applications: Practice Final Exam - Solutions

SYSTEMTEORI - ÖVNING 5: FEEDBACK, POLE ASSIGNMENT AND OBSERVER

Stabilisation of Linear Time-Invariant Systems subject to Output Saturation

Cayley-Hamilton Theorem

Homework 6 Solutions. Solution. Note {e t, te t, t 2 e t, e 2t } is linearly independent. If β = {e t, te t, t 2 e t, e 2t }, then

ACM 104. Homework Set 5 Solutions. February 21, 2001

Linear Algebra. Linear Algebra. Chih-Wei Yi. Dept. of Computer Science National Chiao Tung University. November 12, 2008

LECTURE 6: VECTOR SPACES II (CHAPTER 3 IN THE BOOK)

ECE504: Lecture 8. D. Richard Brown III. Worcester Polytechnic Institute. 28-Oct-2008

3 (Maths) Linear Algebra

Transcription:

Chapter 6 Controllability and Obervability Controllability: whether or not the state-space equation can be controlled from input. Observability: whether or not the initial state can be observed from output.

6.2 Controllability Consider the n-dimensional p-input equation & = + Bu Definition 6.2 The pair (, B) is said to be controllable if for any initial state () = and any final state 1, there eists an input that transfer to 1 in a finite time. Eample 6.1: Figure 6.2 (a) and (b) is not controllable.

Theorem 6.1 The following statement are equivalent. 1. The n-dimensional pair (, B) is controllable 2. The n n matri W (t) = c t e τ τ t (t τ) (t τ) BB e dτ is nonsingular for any t >. 3. The n np controllability matri C =[B B 2 B n-1 B] has rank n (full row rank). = e BB e dτ

4. The n (n+p) matri [-λi B] has full row rank at every eigenvalue, λ, of. 5. If all eigenvalues of have negative real parts, then the unique solution of Wc + Wc T = -BB T is positive definite. The solution is called the controllability Gramian and can be epressed as W c = τ e BB e τ dτ

6.2.1 Controllability indices Let and B be n n and n p constant matri. If (, B) is controllable, its controllability matri C = [B B 2 B n-1 B] = [b 1 b p b 1 b p n-1 b 1 n-1 b p ] has rank n and consequently, n linearly independent column.

If C has rank n, then µ 1 + µ 2 +... +µ p = n. The set {µ 1, µ 2,,µ p } is called the controllability indices and µ = ma(µ 1, µ 2,,µ p ) is called the controllability inde of (, B). Consequently, if (, B) is controllable, the controllability inde µ is the least integer such that ρ(c µ ) =ρ([b B 2 B µ-1 B]) = n.

The controllability inde satifies n / p µ min(n, n p + 1) where r(b) = p, and n B {B,B,..., n 1 B}(Linear combination) Corollary 6.1 The dimensional pair (, B) is controllable if and only if the matri C n-p+1 := [B B n-p B] where ρ(b) = p, has rank n or the n n matri C n-p+1 C T n-p+1 is nonsingular.

Theorem 6.2 The controllability property is invariant under equivalence transformation. Theorem 6.3 The set of controllability indices (, B) is invariant under any equivalence transformation and any reordering of the columns of B.

6.3 Observability The concept of observability is dual to that of controllability. Controllability studies the possibility of steering the state from input. Observability studies the possibility of estimating the state from output.

Consider the n-dimensional p-input q- output state equation & = + Bu y = C + Du Definition 6.O1 The above state equation is said to be onservable if for any unknow initial state (), there eists a finite t 1 > such that the knowledge of the input u and the output y over [, t 1 ] suffices t determine uniquely the initial state ().

Theorem 6.4 The state equation is observable if and only if the n n matri t τ τ W = e C Ce dτ is nonsingular for any t >. Theorem 6.5 (Theorem of duality) The pair (, B) is controllable if and only if the pair ( T, B T ) is observable.

Theorem 6.O1 The following statement are equivalent. 1. The n-dimensional pair (, C) is observable. 2. The n n matri t τ τ W = τ is e C Ce d nonsingular for any t >. 3. The nq n observability has rank n (full column rank). 4. The (n+q) n matri λi has full column rank at C every eigenvalue. O = C C. C n 1

5. If all eigenvalues of have negative real parts, then the unique solution of T W + W T = -C T C is positive definite. If (, C) is observable, its observability matri O has rank n. Let v m be the number of the linearly independent rows associated with c m. If O has rank n, then v 1 + v 2 + + v q = n

The set {v 1,v 2,,v q } is called observablity indices and v = ma(v 1,v 2,,v q ) is called the observablility inde of (, C). If (, C) is observable, it is the least integer such that ρ(o C C ) : = C. v C 2 v = Dual to the controllability part, 1 n / q v min(n, n q + 1) where r(c)=q and n is the degree of the minimal polynomial of. n

Corollary 6.O1 The dimensional pair (, C) is observable if and only if the matri has full rank or the n n matri O T n-q+1o n-q+1 is nonsingular. = + q n 2 1 q n C. C C C O

Theorem 6.O2 The observability property is invariant under any equivalence transformation. Theorem 6.O3 The set of observability indices of (, C) is invariant under any equivalence transformation and any reordering of the rows of C.

6.4 Canonical decomposition Theorem 6.6 Consider the n-dimensional state equation with ρ(c) = ρ([b B n-1 B]) = n 1 < n. We form the n n matri P -1 := [q 1 q n 1 q n ] where the first n 1 columns are any linearly independent columns of C, and the remaining columns can arbitrarily be chosen as long as P is nonsingular.

Then the equivalent transformation will obtain where is n 1 n 1 and is (n-n 1 ) (n-n 1 ), and the n 1 -dimensional subequation is controllable. u B c c c c 12 c c c + = & & [ ] Du C C y c c c c + = c c

Theorem 6.O6 Consider the n-dimensional state equation with We form the nn matri where n 2 rows are any n 2 linearly independent rows of O. n n C. C C O) ( 2 1 n < = = ρ ρ = n n 1 p. p. p P 2

Then the equivalence transformation will obtain where is n 2 n 2 and is (n-n 2 ) (n-n 2 ), and the n 2 -dimensional subequation is observable. u B B o o o o o 12 o o o + = & & [ ] Du C y o o o + = o o

Theorem 6.7 Every state-space equation can be transformed into the canonical form u B B co co co co co co co 43 co 24 23 co 21 13 co co co co co + = & & & & [ ] Du C C y o c co + =

6.5 Conditions in Jordan-form equations If a state equation is transformed into Jordan form, then the controllability and observability conditions become very simple. Consider the state equation & = J + Bu y = C ssume that J has only two distinct eigenvalues and can be written as J = diag(j 1, J 2 )

gain assume that J 1 has three Jordan blocks and J 2 has two Jordan blocks or J 1 = diag(j 11, J 12, J 13 ) J 2 = diag(j 21, J 22 ) The row of B corresponding to the last row of J ij is denoted by b lij. The column of C corresponding to the first column of J ij is denoted by c fij.

Theorem 6.8 1. The state equation in Jordan form is controllable if and only if {b l11, b l12, b l13 } are linearly independent and {b l21, b l22 } are linearly independent. 2. That is observable if and only if {c f11, c f12, c f13 } are linearly independent and {c f21, c f22 } are linearly independent.

6.6 Discrete-time state equations Consider the n-dimensional p-input q-output state equation [k+1] = [k] + Bu[k] y[k] = C[k] Definition 6.D1 The discrete-time state equation or the pair (, B) is said to be controllable if for any initial state () = and any final state 1, there eists an input sequence of finite length that transfers to 1.

Theorem 6.D1 The following statements are equivalent: 1. The n-dimensional pair (, B) is controllable. 2. The n n matri W = n 1 m m= m dc [n 1] () BB ( ) is nonsingular. 3. n np controllability matri C d = [B B... has rank n (full row rank). 2 B n 1 B]

4. The n (n+p) matri [-λi B] has full row rank at every eigenvalue, λ, of. 5. If alleigenvalues of has magnitudes less than 1, then the unique solution of W dc -W dc T = BB T (Controllability Gramian) is positive definite.

Definition 6.D2 The discrete-time state equation or the pair (, C) is said to be observable if for any unknown initial state [], there eists a finite integer k 1 >, such that the knowledge of the input sequence u[k] and the output sequence y[k] from k = to k 1 suffices to determine uniquely the initial state [].

Theorem 6.DO1 The following statements are equivalent: 1. The n-dimensional pair (, C) is observable. 2. The n n matri W do n m is nonsingular (positive definite). 3. The nq n observability matri O d = has rank n (full column rank). m [n 1] = 1 = ( ) C C C C. C n 1 m

4. The (n+q) n matri λi C has full column rank at every eigenvalue, λ, of. 5. If all eigenvalues of have magnitude less than 1, then the unique solution of W do - T W do = C T C (observability Gramian) is positive definite.

There are three different controllability definitions: 1. Definition 6.D1. 2. Controllability to the origin: transfer any state to the zero state. 3. Reachability (controllability from the origin): Transfer the zero state to any state. In continuous-time case, because e t is nonsingular, the three definitions are equivalent.

In discrete-time case, if is nonsingular, the three definitions are again equivalent. But if is singular, then 1 and 3 are equivalent, but not 2 and 3.

6.7 Controllability after sampling Consider a continuous-time state equation &(t) = (t) + If the input is piecewise constant or u[k] := u[kt] = u(t) for kt t<(k+1)t then the equation can be described by [k + 1] = [k] + with = T e B = ( T Bu(t) Bu[k] a t dt)b = : MB

Theorem 6.9 suppose (6.65) is controllable. sufficient condition for its discretized equation in (6.66), with sampling period T, to be controllable is that Im[λ i -λ j ] 2πm/T for m=1, 2,, whenever Re[λ i -λ j ] =. Theorem 6.1 If a continuous-time linear timeinvariant state equation is not controllable, then its discretized state equation, with any sampling period, is not controllable.

6.8 LTV state equations Consider the n-dimensional p-input q- output state equation & = (t) + y = C(t) B(t)u In the time-varying case, the specification of t and t 1 is crucial.

Theorem 6.11 The n-dimensional pair ((t), B(t)) is controllable at time t if and only if there eists a finite t 1 > t such that n n matri W (t c, t 1 t τ τ τ Φ t 1, )B( )B ( ) (t1 ) = 1 Φ(t, τ) dτ where Φ(t, τ): the state transition matri is nonsingular.

The controllability condition without involving Φ(t, τ): Define M (t) = B(t), and then define recursively a sequence of n p matrices M m (t) as Mm + 1(t) : = (t)mm(t) + we have Using Φ(t, t) = (t, t)(t), we compute d dt M Φ( t2, t)b(t) = Φ(t2, t)m(t) t 2 Φ 2 m (t) d [ Φ(t2, t)b(t)] = [ Φ(t2, t)]b(t) + Φ(t2, t) B(t) t t dt d = Φ( t2, t)[ (t)m(t) + M(t)] = Φ(t2, t)m 1(t) dt

Thus, we have t m m Φ(t 2, t)b(t) = Φ(t 2, t)m m (t) Theorem 6.12 Let (t) and B(t) be n-1 times continuously differentiable. Then the n- dimensional pair ((t), B(t)) is controllable at t if there eists a finite t 1 >t such that [ M (t ) M (t )... M (t )] n rank 1 1 1 n 1 1 = Theorem 6.O11 The pair((t), C(t)) is observable at time t if there eists a finite t 1 >t such that the n n matri t t W (t, t ) = 1 o 1 Φ ( τ, τ )C ( τ)c( τ) Φ( τ, t) dτ

Theorem 6.O12 Let (t) and C(t) be n-1 times continuously differentiable. Then the n-dimensional pair ((t), C(t)) is observable at t if there eists a finite t 1 >t such that N N rank Nn where 1 (t1) (t ) 1 =. 1(t1) Nm + 1(t) = Nm(t)(t) + with N = C(t) n d dt N m (t)