MEM633 Lectures 7&8. Chapter 4. Descriptions of MIMO Systems 4-1 Direct Realizations. (i) x u. y x

Similar documents
APPENDIX A Some Linear Algebra

5 The Rational Canonical Form

2.3 Nilpotent endomorphisms

Perron Vectors of an Irreducible Nonnegative Interval Matrix

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

DIFFERENTIAL FORMS BRIAN OSSERMAN

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

Google PageRank with Stochastic Matrix

D.K.M COLLEGE FOR WOMEN (AUTONOMOUS), VELLORE DEPARTMENT OF MATHEMATICS

Quantum Mechanics I - Session 4

18.781: Solution to Practice Questions for Final Exam

NOTES ON SIMPLIFICATION OF MATRICES

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

DISCRIMINANTS AND RAMIFIED PRIMES. 1. Introduction A prime number p is said to be ramified in a number field K if the prime ideal factorization

First day August 1, Problems and Solutions

FACTORIZATION IN KRULL MONOIDS WITH INFINITE CLASS GROUP

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Dynamic Systems on Graphs

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

STAT 309: MATHEMATICAL COMPUTATIONS I FALL 2018 LECTURE 16

MAXIMAL INVARIANT SUBSPACES AND OBSERVABILITY OF MULTIDIMENSIONAL SYSTEMS. PART 2: THE ALGORITHM

9 Characteristic classes

Math 217 Fall 2013 Homework 2 Solutions

Homework Notes Week 7

MTH 819 Algebra I S13. Homework 1/ Solutions. 1 if p n b and p n+1 b 0 otherwise ) = 0 if p q or n m. W i = rw i

GELFAND-TSETLIN BASIS FOR THE REPRESENTATIONS OF gl n

ALGEBRA HW 7 CLAY SHONKWILER

Affine transformations and convexity

CSCE 790S Background Results

U.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Zeros and Zero Dynamics for Linear, Time-delay System

6. Stochastic processes (2)

6. Stochastic processes (2)

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Distribution det x s on p-adic matrices

SL n (F ) Equals its Own Derived Group

Continuous Time Markov Chain

Vector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.

7. Products and matrix elements

Matrix Approximation via Sampling, Subspace Embedding. 1 Solving Linear Systems Using SVD

Norms, Condition Numbers, Eigenvalues and Eigenvectors

MATH Sensitivity of Eigenvalue Problems

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

An explicit solution to polynomial matrix right coprime factorization with application in eigenstructure assignment

PRIMES 2015 reading project: Problem set #3

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Math 101 Fall 2013 Homework #7 Due Friday, November 15, 2013

FINITELY-GENERATED MODULES OVER A PRINCIPAL IDEAL DOMAIN

NOTES FOR QUANTUM GROUPS, CRYSTAL BASES AND REALIZATION OF ŝl(n)-modules

Refined Coding Bounds for Network Error Correction

1 (1 + ( )) = 1 8 ( ) = (c) Carrying out the Taylor expansion, in this case, the series truncates at second order:

Mixed-integer vertex covers on bipartite graphs

), it produces a response (output function g (x)

332600_08_1.qxp 4/17/08 11:29 AM Page 481

MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS

MEM Chapter 4b. LMI Lab commands

Random Walks on Digraphs

Problem Set 9 Solutions

Lorentz Group. Ling Fong Li. 1 Lorentz group Generators Simple representations... 3

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

Linear Algebra and its Applications

Polynomials. 1 What is a polynomial? John Stalker

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

Math1110 (Spring 2009) Prelim 3 - Solutions

Representation theory and quantum mechanics tutorial Representation theory and quantum conservation laws

Errata to Invariant Theory with Applications January 28, 2017

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

REAL ANALYSIS I HOMEWORK 1

Lecture 21: Numerical methods for pricing American type derivatives

Solutions to the 71st William Lowell Putnam Mathematical Competition Saturday, December 4, 2010

On Finite Rank Perturbation of Diagonalizable Operators

From the Euclidean Algorithm for Solving a Key Equation for Dual Reed Solomon Codes to the Berlekamp-Massey Algorithm

A summation on Bernoulli numbers

w ). Then use the Cauchy-Schwartz inequality ( v w v w ).] = in R 4. Can you find a vector u 4 in R 4 such that the

hapter 6 System Norms 6. Introducton s n the matrx case, the measure on a system should be nduced from the space of sgnals t operates on. hus, the sze

28 Finitely Generated Abelian Groups

Solutions to exam in SF1811 Optimization, Jan 14, 2015

NP-Completeness : Proofs

topology If p = q and det(s + ), the pencl s + s called a regular pencl Pencls not satsfyng these condtons are referred to as sngular pencls A complet

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

Min Cut, Fast Cut, Polynomial Identities

Errors for Linear Systems

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Math 594. Solutions 1

6) Derivatives, gradients and Hessian matrices

Outline and Reading. Dynamic Programming. Dynamic Programming revealed. Computing Fibonacci. The General Dynamic Programming Technique

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

Foundations of Arithmetic

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Lecture 4: Constant Time SVD Approximation

Computing π with Bouncing Balls

2 More examples with details

COMPUTING THE NORM OF A MATRIX

Modelli Clamfim Equazioni differenziali 7 ottobre 2013

A graph is a pair G= (V,E) with V a nonempty finite set of nodes or vertices V { v1 a set of edges or arcs E V V. We assume ( vi, vi) E,

Transcription:

MEM633 Lectures 7&8 Chapter 4 Descrptons of MIMO Systems 4- Drect ealzatons y() s s su() s y () s u () s ( s)( s) s y() s u (), s y() s u() s s s y() s u(), s y() s u() s ( s)( s) s () ( s ) y ( s) u ( s) x x u y x () ( s ) y( s) u( s) x x u y x () ( s 3s) y ( s) u ( s) x 3 3 x3 u x 4 0 x 4 0 3 y 0 x x 4

(v) ( s ) y( s) u( s) x x u y 5 5 x 5 Combne (), (), (), and (v), we have x x 0 x x 0 u x 3 3 x3 0 u x 4 0 x 4 0 0 x x 0 5 5 x x 0 0 0 y x 3 0 0 0 x4 x 5

Block Controller Form 3 us ( ) Hs ( ) y( s) Ns Hs ( ) ( ) ds ( ) Ns ( ) s a polynomal matrx n s and ds ( ) s the least common multple of the denomnator polynomals n Hs. ( ) ewrte Hs ( ) as Hs ( ) Ns ( ) D( s ) where D ( s) d( s ). m Now, Defne then ys Ns D s us ( ) ( ) ( ) ( ) () s D () s u() s y() s N() s () s D () s () s u() s

Example: 4 u _ s ds ( ) s asa Ns ( ) Ns N a a u y N N -a m -a m. _ s N N y x am amx m u x m 0 x 0 y N N x x The dmenson of the state space s r m where r s the degree of ds ( ).

Block Observer Form 5 us ( ) Hs ( ) y( s) H() s Ns () ds () Ns ( ) s a polynomal matrx n s and ds ( ) s the least common multple of the denomnator polynomals n Hs. ( ) ewrte Hs ( ) as H s D s N s () L () () where D () s d() s. Now, Defne then L p ys D s Nsus () L() ()() () s D () s y() s L ys D s s () L() () () s Nsus ()()

Example: 6 Defne then ds ( ) s asa Ns ( ) NsN () s DL () sys () ( s as a) ys () y s as y a s y ( sn N ) u u N -a p _ s N -a p _ s y x a p px N u x a p 0 x N x y p 0 x The dmenson of the state space s r p where r s the degree of ds ( ).

Glbert Dagonal ealzatons 7 H() s Ns () ds () Suppose ds ( ) has dstnct roots,.e., r, j ds () ( s ) Then Hs ( ) can be wrtten as where Let H() s r s lm( s ) H( s) s rank and decompose as CB, C : p, B : m Then t s easy to verfy that A = block dag {,,,..., r } T T T T B B B.. B r, C C C C.. r s a realzaton of H () s. The dmenson of the state r space s.

Example: 8 s s H() s s s ( s)( s) s = = A = B = C = Dmenson of A = + = 3.

4- Observablty and Controllablty 9 Def: observablty The lnear tme-nvarant system x () t Ax() t Bu() t yt () Cxt () A: nxn, B: nxm, C: pxn s sad to be observable f the ntal state x(0 - ) can be unquely determned from measurements of the output y(t) over a fnte nterval of tme, 0 - t t f. Theorem: The n-dmensonal lnear tme-nvarant system x () t Ax() t Bu() t yt () Cxt () s observable f and only f the row vectors of the npxn observablty matrx O C CA : n CA span the n-dmensonal space,.e., rank O = n.

Def: controllablty 0 The lnear tme-nvarant system x () t Ax() t Bu() t s sad to be controllable f the state of the system can be transferred from the ntal zero state x(0 - ) = 0 to any fnal state x(t f ) n a fnte tme t f 0, by some control u(t). Theorem: The lnear tme-nvarant system x () t Ax() t Bu() t s controllable f and only f t can be transferred from any ntal state x(0 - ) to any fnal state x(t f ) n a fnte tme t f 0, by some control u(t). Theorem: The n-dmensonal lnear tme-nvarant system x () t Ax() t Bu() t s controllable f and only f the column vectors of the nxnm controllablty matrx C n B AB.. A B span the n-dmensonal space,.e., rank C = n.

The Controllablty Index and the Observablty Index Suppose that the controllablty matrx has rank n. Let C n B AB.. A B q C q B AB.. A B, q n Then the smallest q, say, such that C q has rank n s called the controllablty ndex of {A,B}. Suppose that the observablty matrx O C CA : n CA has rank n. Let O q C CA : q CA Then the smallest q, say, such that O q has rank n s called the observablty ndex of {C,A}.

Smlarty Transformaton Same as scalar case. Theorem: The transfer functon, characterstc polynomal, and ranks of the observablty and controllablty matrces are preserved under smlarty transformaton.

Noncontrollable Canoncal Form 3 Theorem: Let { ABC,, } be such that rank C ( A, B) = r < n Then there exsts a nonsngular T such that the realzaton has the form { A T AT, B T B, C CT } A A 0 A A, B B 0 C C C and the rxr subsystem A, B, C s controllable. Proof: Hnt: Choose T as T T T wth T f f f... r where f, f,..., f r are lnearly ndependent and span the nm columns of C ( AB, ). The columns of T together wth the columns of T span the n- dmensonal space.

Nonobservable Canoncal Form 4 Theorem: Let { ABC,, } be such that ranko ( C, A) = r < n Then there exsts a nonsngular T such that the realzaton { A T AT, B T B, C CT } has the form A 0 A A A, B B B C C 0 and the rxr subsystem A, B, C s observable. Proof: Hnt: Choose T such that U T U where the rows U of are lnearly ndependent and span the np rows of O ( CA, ). The rows of U together wth those of U span the n-dmensonal space.

Mnmal ealzatons 5 Def: mnmal realzaton A mnmal realzaton of a gven transfer functon Hs ( ) s one that has the smallest-sze A matrx for all trples {A,B,C} satsfyng Theorem: Cs ( A) BHs. ( ) A realzaton {A,B,C} s mnmal f and only f t s controllable and observable. Proof: "Necessty": If {A,B,C} s ether noncontrollable or nonobservable, then by the theorems n p.3 and 4, we could obtan another realzaton A, B, C wth the same transfer functon but a smaller number of states. Therefore {A,B,C} s not mnmal f t s noncontrollable or nonobservable. "Suffcency": To prove that a controllable and observable realzaton { ABC,, } s mnmal. Suppose t s not mnmal, and let { ABC,, } be another controllable and observable realzaton wth a lower state dmenson, n<n. Snce CA B CA B all we have

OC OC n n 6 From Sylvester s nequalty, we have ( O ) ( C )- n ( OC ) mn{ ( O), ( C )} Snce O and C have rank n, hence ( OC ) n But { ABC,, } s also controllable and observable, t must follow smlarly that ( OC ) n n n Thus, n n. It s a contradcton. Therefore a controllable and observable realzaton {A,B,C} must be mnmal. Theorem: Suppose that A, B, C,, are two mnmal realzatons of a transfer functon, then there exsts a unque nonsngular matrx T such that and A T AT, B T B, C CT T CC ( CC ) T T T T T ( OO) O O

Theorem: 7 The Glbert realzaton s mnmal. Proof: A block dag,,..., r B B B B r T T T T.. C C C C.. r The controllablty matrx can be wrtten as C = n B AB.. A B = B. V It s easy to see that ( V ) rm f j.

By Sylvester s nequalty 8 ( B ) ( V )- rm ( C ) mn{ ( B ), ( V )} we have ( C ) ( B ). By the Glbert constructon, r ( B ) ( B ) n r so that the realzaton s controllable. The observablty can be smlarly proved.

PBH Tests for Controllablty and Observablty 9 Theorem: PBH Egenvector Tests. { A, B } s noncontrollable f and only f there exsts a nonzero row vector q such that qa = q, and qb = 0.. { C, A } s nonobservable f and only f there exsts a nonzero column vector p such that Ap = p, and Cp = 0. Theorem: PBH ank Tests. { A, B } s controllable f and only f rank s A B n for all s. { C, A } s observable f and only f sa rank n C for all s Here, n s the sze of A.

4-3 Matrx -Fracton Descrptons 0 us ( ) Hs ( ) y( s) Ns Hs ( ) ( ) ds ( ), deg ds ( ) r Wrte Hs ( ) as a matrx fracton, Hs ( ) N ( s) D( s) Wth D( s) d( s) m, N( s) N( s) and defne the degree of the denomnator matrx as deg D ( s) deg det D ( s) rm then the order of the block controller realzaton would be equal to r m, the degree of D ( s ). Smlarly, f we wrte Hs D s N s ( ) L( ) L( ) Wth DL( s) d( s) p, NL( s) N( s) then the degree of the denomnator matrx s deg D ( s) deg det D ( s) rp L and the order of the block observer realzaton s equal to r p, the degree of DL ( s ). L

Clam A: Gven any rght MFD of Hs, ( ) Hs ( ) N ( s) D( s) we can always obtan a controllable state-space realzaton {A,B,C} of order Clam B: n = deg det D ( s ) := the degree of the MFD Gven any left MFD of Hs, ( ) Hs D s N s ( ) L( ) L( ) we can always obtan an observable state-space realzaton {A,B,C} of order n = deg det DL ( s ) := the degree of the MFD An MFD wth mnmal degree s sad to be a mnmal MFD (rreducble MFD). A realzaton of a mnmal MFD s a mnmal realzaton.

Example: s s Hs ( ) ( s)( s) s =

ght Dvsors and Irreducble MFDs 3 Gven an MFD Ds of ( ) Ns ( ) ( ) Hs. If we defne N s s, N ( s) ( ) ( ) D( s) D( s) ( s) where s ( ) s any nonsngular polynomal matrx, then Hs ( ) Ns ( ) Ds ( ) N( s) D( s).e., N s D s s also an MFD of ( ) ( ) ( ) Hs. Note that Ns ( ) N( ss ) ( ) Ds ( ) D( ss ) ( ), s ( ) s called a rght dvsor of Ns ( ) and Ds ( ). Snce and det Ds ( ) det D( s) det s ( ) degdet Ds ( ) deg det D( s) degdet s ( ) we have deg det Ds ( ) deg det D( s ) It s easy to see that the degree of the MFD can be reduced by removng rght common dvsors of the numerator and denomnator matrces. Therefore, we can get a mnmum-degree MFD by extractng a greatest common rght dvsor (gcrd) of Ns ( ) and Ds ( ).

The MFD and only f NsDs ( ) ( ) s mnmal (rreducble) f 4 deg det Ds ( ) = deg det D( s ) for all nonsngular rght dvsors s ( ) of Ns ( ) and Ds ( ). That s, det s ( ) = a nonzero constant ndependent of s A polynomal matrx whch satsfes the above equalty s called unmodular. Two polynomal matrces Ns ( ) and Ds ( ) wth the same number of columns s relatvely rght prme (rght coprme) f they only have unmodular common rght dvsors. And the MFD NsDs ( ) ( ) s rreducble (mnmal) f Ns ( ) and Ds ( ) are rght coprme. Irreducble MFDs are not unque.