Lecture Notes in Control and Information Sciences: 107 Y. T. Tsay, L.-S. Shieh, S. Barnett
Lecture Notes in Control and Information Sciences: 107 Y. T. Tsay, L.-S. Shieh, S. Barnett
Lecture Notes in Control and Information Sciences: 107 Y. T. Tsay, L.-S. Shieh, S. Barnett
Control and
Information Sciences
Edited by M.Thoma and A.Wyner
107
Springer-Verlag
Berlin Heidelberg New York
London Paris Tokyo
S e r i e s Editors
M. T h o m a ° A. W y n e r
Advisory Board
L. D. Davisson • A. G. J. MacFarlane • H. Kwakernaak
J. L. Massey • Ya Z. Tsypkin - A. J. Viterbi
Authors
Yih Tsong Tsay
Leang-San Shieh
D e p a r t m e n t of E l e c t r i c a l E n g i n e e r i n g
U n i v e r s i t y of H o u s t o n
Houston, Texas 77004
USA
Stephen Barnett
S c h o o l of M a t h e m a t i c a l S c i e n c e s
U n i v e r s i t y of B r a d f o r d
W e s t Y o r k s h i r e B D 7 1DP
United Kingdom
Progress in system theory over the last two decades can be broadly categorized
(I) Algebraic System Theory - Study of basic notions and fundamental concepts
(2) System Analysis and Design Methods - Study of potential design techniques
Thousands of papers have been published in both areas in the last two decades.
Systemic presentations in book form can be found~ for example~ in [1-5] for the
former, in [6-10] for the latter~ and in [11-16] for both. From this
analyze/design the physical systems with ease. In other words~ work is still
needed to fill the gap between algebraic system theory and practical system
monograph.
The development of our work is based upon state-space representations and matrix
Offlce~ Dr. Robert E. Yates, Director of the Guidance and Control Directorate,
US Army Missile Com=~and, and Dr. Norman P. Colemanj Chief of the Automation and
D e v e l o p m e n t and E n g i n e e r i n g C e n t e r .
contract NAGg-211.
TABLE OF CONTENTS
Chapter I. INTRODUCTION
4.1 Structure theorems for canonical left divisors and complete sets of
canonical left divisors . . . . . . . . . . . . . . . . . . . . . . . . 75
4.2 Structure theorems for canonical right divisors and complete sets of
canonical right divisors . . . . . . . . . . . . . . . . . . . . . . . 85
5.4 Divisor assignment and decoupling design of MIMO systems ....... 142
BIBLIOGRAPHY . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 200
INDEX . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 207
CHAPTER I INTRODUCTION
I.I. Some basic definitions on k-matrices~ which are the main mathematical
tools used in our workj are summarized in Section 1.2~ and Section 1.3 gives a
Multivariable Systems
where X(t)~x~cnp y(t)EF~C p, u(t)FJJ=_Cm are state, output, and input vectors,
linear maps:
X A = X
T
U
D ,~ ;,
(1.3)
the system O.
represented as
(1.4~)
y(t) = G(1)u(t)
where
In Eq. (l.4b) cPXm(1) denotes the set of pXm matrices with elements being
rational functions of ~ over the complex field C. G(I) is called the transfer
funct{on matrix of the system O. It has been shown in [i~13] that G(1) can be
= Nr(I)D~I(I) (l.5b)
where DI(I) EcP×P[~,], Ng(l) , Nr(l) ecPXm[)~], Dr(%)ecmXm[% ] ; cPxP[)~], Cpxm[)~] and
cmxm[l] are sets of matrix polynomials of ~ with coefficients in cPXP~ CPXm~ and
(l.6b)
: Nr(X)D~l(x)u(t)
Equations (l.6a) and (l.6b) are referred to as left matrix fraction descriptions
respectively.
Let TEC n×n be a nonsingular matrix~ and from Eq. (i.I) define
and
^ A
For the same set of inputs u(t), o in Eq. (i.i) and o in Eq. (1.8) will generate
the same set of outputs y(t) for t~0 if X(O) = TX(0). The difference between
the state vectors X(t) and (t) in the system O and O~ respectively~ is
Thus~ we say that o and O are equivalent systems. Formally, we have the
following definition=
Definition I.i The system in Eq. (i.I) and the system in Eq. (1.8) are
x(t) = TX(t)
We will call this equivalence relation similarity equivalence (SE). Q
and
From Eqs. (l.6a) and (I.i0), ~ in Eq. (l.6a) and O in Eq. (I.i0) will generate
the same set of y(t) for t~0 if the same set of u(t) is used as inputs~ and
^
and ~ both have the same set of initial conditions y(t), t<0. Thus, we say that
and ~ are equivalent systems. Similar reasoning can be applied for RMFDs. We
Definition 1.2 Two systems with LMFDs G(~) = D~I(~)N~(~) and G(~) =
D~(X) = U~(X)D~(X)
and
5
where U~(%) is unimodular. Similarly, two systems with RMFDs O(~) m Nr(~)Drl(~)
Dr(1) = Dr(X)Ur(X)
and
NrC~) = NrC~)UrC~)
It can easily be verified that both SE and UE satisfy the basic properties
descriptions, we have:
Chapter II, we will develop canonical forms, which are unique for a given
involve the ratio of two %-matrices~ and the results presented in the following
be found, for example, in [2] and [3]. Specifically, we can define l-matrices
matrix whose elements are in Fill. Let A..(X) be the (i,~)th element of A(X))
i]
then
and k,°
Aij(X) ~ ~J kij-k
k=0 aijkX ) aijk~F (l.llb)
r
A(X) = ~ AkXr-k (l.12a)
k=O
regular if the matrix coefficient A 0 of the highest degree term (referred to Eq.
matrix.
Then, ~i' denoted by ~i = Dr. (A(I)), is the row degree [12] of the ith row of
1
A(X). Similarly
denoted by Kj : ~c.(A(~)), is the column degree [12] of the jth column of A(X).
J
Define
(1.15a)
Ahr = ((Ahr)ij), IKiKp, IKjKm
where
aiJ~i if kij = ~i
(I. 15b)
(Ahr)ij
L 0 if kij < ~i
Then Ahr is called the leading row matrix of A(X). A(X) called is a row-reduced
(l.16a)
Ahc = ((Ahc)ij) , l~igp, l~j~m
where
I aijKj if kij = ~j
(l.16b)
(Ahc) ij
L 0 if kij < <j
then Ahc is called the leading column matrix of A(X), and if p=m and Ahc is
[12,13] .
Definition 1.3 Two ~-matrices Al(l ) and A2(1) are row equivalent, column
matrices. []
matrix. •
authors [17-32]. We shall extend some known results on monic h-matrices to row-
representations of the systems are rational matrices over the co~plex field C.
involved.
The material in this monograph can be regarded as being in two parts~ The
first part7 which includes Chapters II, III and IV~ is devoted to exploring the
systems represented in state space equations and MFDs~ the second part~ which
developed in the first part to the design and decomposition of MIMO systems.
systems are defined. For a reachable system the characteristic X-matrix can be
constructed from the coefficients of the dependence equations for the column
vectors of the teachability test matrix; on the other hand~ for an observable
formally defined. The canonical RMFDs and LMFDs for reachable and observable
systems~ respectively~ are defined~ and their properties are discussed based on
selections.
matrices and the eigenstructures of the system maps in their associated state-
chains of nonsingular X-matrices can be easily found from the input and output
nonsingular %-matrix into the product of lower degree canonical l-matrices are
matrices.
closed-loop MFDs are derived. A study is then made properties of the closed-
assignment and decoupling design, via the notions of divlsors~ are also
presented.
polynomials, which are the denominators of the transfer functions. Many system
concepts for MIMO systems, in this chapter we present formal definitions of the
characteristic k-matrices and the canonical MFDs for MIMO systems, which are the
the characteristic k-matrices and the canonical MFDs are given. Relationships
between the state space representations and MFDs of MIMO systems are
established. The results of this chapter provide the foundations for the
spectral analysis of k-matrices and the decompositions of MIMO systems which are
where X(t) ECn, y(t) E ~ and u(t) ECm are state, output, and input vectors,
matrix [13]
I] (2.3)
is of full rank.
Definition 2.1 Let B = [bl,b2,...,bm] ~ bieC nxl. The reachabillty base matrix
is defined by
<i-i ~ -I (2.4a)
P(A~B) =A [bl,Abl~" . o ~ A bl, . . ,bm,Abm~.
. . . ~A m b m ]
Definition 2.2 T T
Let C = [CI~C29...~C ~ ]T CieclXn" The observability base matrix
is defined by
the sequence proposed by Popov [35]. It is well known [33-36] that when (A~B)
m
is a reachable pair, rank[F(A,B)]=n and ~ <i=n; when.(A~C) is an observable
P i=l
pair, then rank[~(AgC) ] = n and [ 9.=n.
i= I i
Assume that (AgB) is a reachable pair with reachability indices ~j
K.
j=l,2,...,m. Then A 3b. can be uniquely represented as [35-37]
3
(2.5a)
and
j-I
(2.5b)
b° -- -
ari'ibi'3 if <. = O, l<jKm
.1 i=l 3
~.>0
I
Note that the input matrix B may contain dependent vectors. As a result~ we can
uniquely represented as
~. i-I
C.A ~ = - ) ag.iJ(vi+l)CjA vi- ~ m~n[~)i'gJ]a~ijkCjAk-I , if Vi>0 , iKi<p
I j=l j:l k--i
V. <~).
i 3
(2.6a)
and
i-i
(2.6b)
C. = - ) a~ijlCj, if v'=0'1 l~i~p
j =~l
V.>O
3
15
It has been shown that the Kronecker indices <.l and the set of parameters
the ~. are also invariant under linear state feedback. Similarly, the
observahility indices ~'x and the set of parameters {a~ijk } are invariants of
(A,C) under coordinate transformations of states and the ~.x are also invariant
(2.7a)
Dr(X) ~ [(Dr(X))ij ] ecmXm[x], l~i~m, l~jKm
A (2.7b)
arii(~.+l ) = I
l
I ~k
if <i>0~ or <.=0z and i=j
arij(k+l)
(Dr(X))ij
0 if Ki=0 , izj
I Ki if i=j (2.7c)
A Kj if <i><j, i~j (2.7d)
<ij =
K.-I
l otherwise Q
A (2.8b)
a~ii(~.+l ) = I
l
9..
(2.8d)
Vij = ~ if 9i<~j, i~j
t gj-I otherwise D
the states, the right characteristic l-matrix, Dr(1) , and the left
A (2.9a)
Drh = [(Drh)ij] , l~iKm, IKjKm
for ~.<~.
j x
A I arij(<j+l) (2.9b)
(Drh)ij = hrij = 0
otherwise
Similarly, when VijN~ i for lKjNp, the row degree of the ith row of D~(k) is ~i"
Define the leading row matrix
(2.10a)
Dih = [(Dzh)ij] , lKiKp, IKjKp
Since a£ii(~i+l ) ~I, and a%iJ(~i+l ) = 0 for i<j, D£h is a lower triangular
matrix with diagonal elements all l's, so det(Dih) = i. Thus, D~(%) is a row
P
reduced l-matrix and the degree of det(D£(~)) is ~ ~. = n.
i= I i
From the definition of the right characteristic ~-matrix in Eq. (2.7) and
Proposition 2.1 The right characteristic l-matrix D (l) in Eq. (2.7) has the
r
following properties.
Proposition 2.2 The left characteristic l-matrix D~(k) in Eq. (2.8) has the
following properties.
18
systems are nonunique. In order to specify the "standard" MFDs for an MIMO
system, we shall define the canonical right MFD for the reachable systems and
the canonical left I~FD for the observable systems in this section. The
transform the state equations in Eq. (2.1) into the canonical controller form
(2.11a)
Xc(t) = TcX(t)
where
and
i
(L = ~° (2.11d)
I j=l 3
19
(2.12a)
XXc(t) = AcXc(t)+Bcu(t)
where
r o
K~. XK~.
EC x x
; Acii = [_~riil,... _~rii<.] (2.13b)
1
_0 ] K.x<.
(Ac)ij ~ .... £C i ]; Acij
~
[-arijl'....-arijKj]' ~i~<j
iZj
Acij =i[-a--rijl, .,-ariJKi,0.....0], ~i<<j
(2.13c)
and
AcmI, •,Acmn J
m
armll,''',arml<l~''',armml~'''~armm<m
213d
I
{ en if <.>01
(2.13f)
ebci 0 if <.=0
CT~I T T T T (2.13g)
Cc = [CcI,Cc2,...,Ccp]
20
A (2.13h)
Cci = [Cill,...,CilK.,...,Ciml,...,Cim ~ ]
i m
D ~ D (z.13i)
e
From Eq. (2.12)~ the input-output relationship of the system can be described as
(2.14a)
y(t) = [Cc(lln-Ac)-IBc+Dc]u(t)
where
(2.14c)
Qr(~) ~ [~r1(1), ....~rm(~)]
and
0 T ~m
[01xOi_l,l,l,..., l I , ix(n-oi )] , Ki>0,i=l,...
(2.14e)
~r(l) ~ ~rh(l)-Ar~r(~)
y(t) = [Nr(k)D~l(X)+Dc]U(r)
(2.15a)
where
(2.16)
G(X) in Eq. (2.165 is referred to the canonical RMFD of the system in Eq. (2.1)~
and D (X) is the right characteristic %-matrix of the RMFD. It is well known
r
that, for a reachable systemp Dr(X) and Nr(%) are right coprime. From Eqs.
(2.14c) and (2.15c), any column of Nr(X) corresponding to Ki'0 is a zero column.
transformation as follows:
~1-1 ~ -1
T0 . [ q l , . . . , A ql'''''qP''''' A p qp] (2.17b)
where
and
T.b j=l
i v.J (2.17d)
22
where
(A0) £i ~ 0 ,
-IAoIi cC~.x~.
i i = -- -- T (2.19b)
_ , A0i £ [-a0iil, •. •, - a 0 i £ u £ ]
-- -- T
V. xV. [-a0ijl .....-a0iju .] , if Vik~j, izj
(A0)Ij ~ [0 A0i j]~CI j; A0ij m -- J
[-a0ijl,...,-a0iJui,0,...,0lT if Vi<vj,izj
(2.19c)
and
A011'""A01p.
. • ] a~lll,''',aillgl,''',a£pll,''',a£plu p
(2.19d)
(2.19e)
Bo ~ T~lB ~ [BOI,B02,...,~O~]
A
BOj = [ b l j l , . . . , b l j v i , . . . , b p j l , . . . , b p j ~ p ] T (2.19f)
- A T T T (2.19g)
CO ~ CT0 = D£~Ec0 ; Eco = [ec01,...,ec0p]
T.
l
en , if ~. >0
T ~ ~ (2.19h)
ecoi = 0 , if ~.=0
l
23
DO ~ D (2.191)
From Eq. (2.18), the input-output relationships of the system can be expressed
as
(2.20a)
yCt) = [Co(XIn-Ao)-IB0+Do]U(t)
where
(2.20c)
*~(~> =~ [~'Lc~:', .... *~,pc~)] '~
and
~.>0
1
ui'0I'A' "'X' z ,01x(n_Ti)],
~i (~) =Af01x-[01xTi-I' (2.20d)
T . t h column o f A0 f o r V.>0
A i l (2.20i)
A£i =
0nX 1 f o r U.=0
24
where
O(k) in Eq. (2.22) is referred to the canonical LHFD of the system in Eq. (2.1)
and D%(k) is the left characteristic k-matrix of the LMFD. D~(k) and N~(k) are
left coprime if Eq. (2.1) is an observable system. Note t~at avy row of N%(k)
Prom Eqs. (2.13) and (2.14), the minimal realization of D-I( ~,,, using a pair
r
(A,B) can be formulated as follows.
Proof:
(2.23b)
_ m ei)T ,T(0),r(0)
K.=0
1
Thus~
~i=O
Since DrCX) = Drh6(k) , we have D-I(X)r = ~ - I ( A ) D ~ " When < i - 0 , from Eq. (2.14e)
M • • m . .
(2.23d)
--_
[I
=
rh <0>*rC0"O
Ki=0 ~i=0
Rearranging the terms in Eq. (2.23c) and substituting Eq. (2.23d) into the
From Lem~a 2.1 we observe that, if every ~i>O for l$i~m, then D r = 0m~ both
C r and B are of full rank, and D-I(A) is strictly proper. For representing
r
Dr(A~ ,, ~rC%~n-A)~(X)÷~',r
A A
Prooft
Using Ebc and A c in Zq. (2.13), *r(A) and A r in Eq. (2.14), and 6rh(A) in
Eq. (2.14f), it can be easily verified that Ar = E~cA c and *rh (A)
E~cASr(A)+[Im-*~(0)$r(0)]. Then, from Eq. (10e) we have
,cA)= ~LA,r(A~÷I~m-,~CO~,rCO,-~L,c,r(A~
T [ X~n_Ae],r(x)+t~m_,~(o),r(o)]
= Ebc
~LTo(A~n~)T~l,rCA~÷t~m ,~(0)~r(0)1
In a similar fashion, from Eqs. (2.19) and (2.20), the minimal realization
Lemmas 2.3 and 2.4 can be proven following the same reasoning as in Lemmas
2.3 Minimal Nice Selections for Determlnln~ Kronecker and Observability Indices
reachable pair (A,B) are known. To find the Kronecker indices of reachable
observable system.
such that
entry iff i+mj~y for l~i~m and O~j~n-1. Due to the property (2) in Definition
2.5, the Young diagram of a nice selection associated with T either has
J~
0 1 2 3 & 5
1 X X b1
£~ 2 X b2
3 X X X b3
z3 s A2 A3 A4 A
Definition 2.6 Let (A,B) be a reachable pair and ~j be the number of crosses in
the jth column o f the Young diagram o f a nice selection y associated with (A,B).
If
k-I
(2.24)
j~0 zj " rank(B'AS'''''Ak'lB)' k~n
and (AJbi), V(i,j)¢{(a,8): l~C~m, O~8~n-I and cr~n(B-1)~y} are independent, then
It has been shown [37] that minimal nice selections are unique up to
permutations of the rows of the Young diagram. Also, if K. denotes the number
x
of crosses at the ith row of the Young diagram of a minimal nice selection YH'
the set K = {~i,l~i~m} is the Kronecker indices of the reachable pair (A,B)
associated with YM" Note that K i may be zero for some i. The maximum number of
minimal nice selections associated with a reachable pair (A,B) is m!, and each
canonical controller form and associated canonical RHFD, we define the input
Definition 2.7 The input selection sequence for a multivariable system with m-
A
input, S = {Sk, l~k~n}, is defined as a permutation of the set {k,l~<Sm}. The
29
A
natural input selectlon sequence, SN~ is de£1ned as SN = { S l ~ S 2 p . * * , s m} =
i s d e t e r m i n e d by e x a m i n i n g t h e i n d e p e n d e n t v e c t o r s o f t h e o r d e r e d column v e c t o r s
i n t h e following m a t r i x :
~(A,B) ~ [ b s l ~ b s 2 j . .,bsmtAbsl~Abs2~
. . *~Abs ~ .* ' A n - ~ s l An-2b ,A n - I b s m]
** m ** ~ 82'"
(2.25a)
D
If the set K contains the Kronecker indices o f t h e minimal n i c e selection
~i-I ~ -1 (2.25b)
PCA,B) = [bl,Abl,...,A bl,...,bm,Abm,..., A m bm ]
and the controller form transformation matrix Tc becomes Eq. (2.11b). From
Eq. ( 2 . 2 5 a ) as f o l l o w s .
Let V = {vi,l~i~N} be
a s e t of v e c t o r s v i ¢ c n × l , N ~ n . The o r t h o g o n a l i z e d
A
projection F(k)~C m~n f o r Vk = { v i , l ~ i ~ k } , l~k~N, can be e v a l u a t e d r e c u r s i v e l y as
follows:
p(O) = I ( t h e u n i t m a t r i x of o r d e r n) (2.26a)
n
30
P(k-l)VkVkP ( k - l )
P(k-1) f o r VkP ( k - l ) VkZO
v~P(k-l)v k
v(k) = (2.26b)
F(k-1) f o r VkP(k-1)v k = 0
where v k d e n o t e s the c o n j u g a t e t r a n s p o s e of v k.
as f o l l o w s :
. A
Lemma 2 . 6 V k P ( k - l ) v k = 0 i f f Vk£Vk_l, where Vk_ I = S p a n c [ V k _ l ] .
Proof;
S i n c e P ( k - 1 ) i s i d e m p o t e n t and s y m m e t r i c , we have
VkP(k-1)v k .= VkP ( k - 1 ) P ( k - 1 ) v k .. ( p ( k - 1 ) V k ) ( P ( k - 1 ) v k)
k-1
(a) If VkE{/k_.,,,. t h e n v k = [ a,v.. From p r o p e r t y (3) o f Lemma 2 . 5 we have
i=l I
k-1
p(k-l)v k = [ a i P ( k - 1 ) v i - Onx I
Thus
VkP(k-1)v k = 0
Let
k-1 ^
" alvi+ k' :k I k-I
11
k-1
T h u s , vk " ~ a i v i o r VkC k - l "
i 1
A
Lennna 2 . 7 D e f i n e WkCVk; W0 = ~ ;
Wk_1 i f V k P ( k - 1 ) v k = 0 f o r l~k~N
Then, WN i s a b a s e f o r VN.
Proof:
Al~orlthm 2.1
Given: (A,B) - A r e a c h a b l e p a i r
S - The i n p u t s e l e c t i o n sequence
Find: ~ - The Kronecker i n d i c e s
Algorithm:
{Initialization}
P:=I n {set orthogonallzad projection to identify matrix)
For i:=l to m Do
Begin
K1. : = 0 ; { R e s e t Kronecker i n d i c e s }
Flagi:=True;{Set selection flags true for all inputs}
v i : = b i {Copy b i to v i}
End;
{processing}
Repeat
For i:=l to m Do
Begin
j : = s i ; { S e l e c t Akbs.,k~0}
1
If F l a g j then
Begin
d:=P*v.;
.J
If v.*dzO then
Begin { S e l e c t v.}
.J
P:=P-d*d*/vj*dl{Update P}
v j : = A * v j ; { U p d a t e vj f o r next s e l e c t i o n }
K.:=~.+l{Update Kronecker index}
J J
End
Else
Flagj:=false {Reset selection Flag}
End
End {For loop}
Until all Flag i are false;
Note that the * between two variables in the above algorithm is the product
notation, and the superscript * designates conjugate transpose. •
33
i n state s p a c e form.
~;i by substituting (AT,cT) for (A,B), ~;i for Ki, and treating S as the output
selection s e q u e n c e w h i c h is d e f i n e d in the same way as the input selection
(2.27a)
;kXCt) = AXCt)+BuCt)
(2.27b)
yCt) - CXCt)+Du(t)
where
-1 0 -2 5 -9
I-3-2]
0 -1 4 -8 14 - 1 1
A = 0 0 1 -4 7 ; B = - 3 5
0 0 0 -1 6 1 4
0 0 0 0 2 0 1
1 -1 4 -5 ; D = 02x 3
indices KI-2 , K2-3 and K3=0. Note that the input matrix B contains a
dependent vector. The c o r r e s p o n d i n g Young d i a g r a m b e c o m e s
0 1 2 4
X X
X X X b2
b3
$5 A
A4
Ac
- I
o 11o
2 I 0 0
8"-8"~"o" 1
0 0 iO
o
0
o]
0
o
1
; B<-/O
/ o
roool
LL_o._I_I
0
o
0 I" o- -0-O-
o/
o
o
1
o
1
o
1
" EbcDThI
[0 011 1 -1 LO 1 2j 1
Cc 3 I 0 2 2 ; Dc = 02x3
12-k-2 0 -1 1 (2.28b)
0 13+k2-k-1 -2
0 0 1
where
35
Drh
-
[i0i]I [ioi]
1
0
= 0
0
1
0
-2
I
; ~rh(k) "
[2oo]
0
0
~3
0
0
1
ar =
[20 01 0I 01 -i0] ; Sr(~)
. Ii O~ 0 0 !2IT
0 0 0 0 0 0 0 0
Note t h a t Drh i s an upper triangular matrix with diagonal entries all l's,
in Proposition 2.1.
The canonical RHFD can be obtained from Eqs. (2.14) and (2.15) as
follows:
where
[2),-4 3),2+3),-2 0 ]
~r(~) ,, C~r(~) -
31 2~.2+2~. 0
0 1 2 3 4
1 b1
2 X X X b2
3 X X b3
15 A A2 A3 A4
o ooo
o o 1!o o/ / o o o 1
li°°IE j
0 0 0 _
A
C
i 1 -1Lo o.l~ ~c'LO 1 2 I 0 I = EbcD r
o o, 6=i=I I-~-6-~ -o- -o" o
o o12 lJ L1 o i 0 1
where
Drh =
[, o
0
1
I
0
i]-' =
[0],0
- 0
-
Cc=
[-~ ~ ~ :-4 2]
2 21 0 3
(,.,,a)
The corresponding RMFD can be obtained from Eqs. (2.14) and (2.15) as
follows:
(2.29b)
y(t) = Nrl(X)D:~(X)u(t)
where
Drl(X) s
[ 0 2 ~k3+~k2-X-1 -2X2+2X+4
o] ; Nrl(X) =
[
0 3X2+3X-2 2k-4 ]
-1 0 X2-X-2 0 2X2+2X 3X J
Note that Drl(X) is a column-reduced X-matrix but it does not satisfy the
(3) If the input selection sequence, S, is selected as S ffi {3)1,2}, then the
Kronecker indices become glffi2) ~2=0, and K3ffi3, and the Young diagram is
A n
II w
0 fl
!
II I 0 @ oll '~ 0
r~ 0 ,.l~, I
I
0
!
~u ,~'u 0,-, ~,'u
I g. ~°
I I..I
0,~00 I 0
t=--..--J ! m
! e
!
I-P
g,, v ~n
v n I I I I~
','.~ (3o n
~g n , =1
÷ I I, I ! ! 0
1,4 I-'
o o i,,,~
,'t
n 0
M ' I ' I'4
I 0 00|P--' 0 ~'1
I
r~
I | II 0 W
"M
t4
iii I
~-'00lO 0 ~'
L I i
0 , ~
II
i v
! ! I
O00l~.J 0
!
0 0 0100
I
I - I 0 0 i O 0
0 0 l I I
÷ I
÷
! !
I I
N
o
v v
38
in Proposition 2.1.
cHAPTER IIl SPECTRAL ANALYSIS OF NONSINGULAR k-MATRICES
c a n o n i c a l HFDs t o c a n o n i c a l MFDs.
"matrix roots", called solvents, of nonslngular k-matrlces are studied via the
latent roots and latent vectors. Solvents of monle k-matrlces have been
solvents in this chapter are: (I) to extend the idea of solvents for
canonical MFDs.
40
Dr(k) = D r h 6 r ( k )
where Drh i s an upper triangular square matrix with diagonal elements being all
l*s. 6r(k) is a l-matrix with the (i,i)th entry, (6r(k))i,i, being a monle A
p o l y n o m i a l of d e g r e e ~ m i n ( K i , ~ j ) - I , i z j .
where D~h i s a lower triangular square matrix with diagonal elements being all
lls. ~(k) is a k-matrix with the (i,i)th entry, (6~(k)i, i being a monic k
can he s y s t e m a t i c a l l y e x p l o r e d i n t e r m s of a s t a t e - s p a c e setting.
reduced k-matrlces.
where Dch i s the l e a d i n g column m a t r i x of A(k) and Dh(k) ~ d i a g [ k Ki ,~'1," ... ,m] ;
K:.-1
+c(~) ^ -~ [~cZ(X),~c2(X),...,~em(X)];^~ci¢X)
. ~)i [O~x o i - 1 , l , X , . . . , x ~ 'o Ix(n-oi)
^ J~T
for I~.>01 and ~c¢(~) =A OnxI for Ki'O ; ai ~ j'l)" K:j, i ' l , . . . , m ; n ~ Ore; Acr ¢Cmxn
can be d e t e r m i n e d from Dch and the c o e f f i c i e n t s of each e n t r y of A ( k ) . •
KI Km
^ iI I } ~cl ^& ^ ^
A .A -(A ). ; n=n-KI-~ m
C T .cr 1__ | . . . . . . . . . .
o:j~ . • , o:j
_.=__
OC~=_I~J1 ,, oc;=_l~x:-',, o l o . o
! • • • • • o
g
O_ 'a _" _ ' _ ' _ _ 1
j_"
-(Acr) m
(3.lb)
e °i i f ~c.>O
Bc = bc - "'''ebcm]; ebcl Onx 1 i f ;i=O
42
~c ~ ~(0) (3.1e)
Dc
Proof:
Lemma 3 . 2 can be e s t a b l l s h e d by d i r e c t verification. •
can be r e a l i z e d by u s i n g a q u a d r u p l e (A,B,CjD)~ or
(3.2)
A-I(A) . C(kln-A)-IB+D
The relationship between A(k) in Eq. (3.2) and Dr(k) in Eq. (3.3) is the
following:
Proof:
From Eqs. (2.6) in Chapter 11 and (3.2)~ we obtain
a unimodular X-matrlx. •
Theorem 3.1 provides a method for determining a unimodular X-matrlx which
converts a nonsingular column-reduced k-matrix to the column-reduced canonical
X-matrix.
The dual results for row-reduced X-matrices are as follows:
Lena 3.3 Let A(k)ecPxP[k] be a row-reduced X-matrix with row degree ~)i'
(3.5)
A()k) - [Vh( k)-$0 (k)A0g I V0h
where D0h is the leading row matrix of A(X) and Dh(X) A diag[lV i = l , . . . , m ] ;
(3.6a)
Lemma 3.4 A-l(x) . Go(X~n-;o)-l;o+~o
where
G
V]
O~lX(~p-1) ^
0 • I A A ^
AO A -(Ao~) 1 eou
O~xc; _~) ,' -(%pp ;nffin-•l-V
. . . . . . . P . . . .
0
I . . A
0 • •
~o ~ *o(°) A
o6o)
T. A
0
-~ cO;
DO
Eco
^ ^
~ [ec01, : c02 ,"- . - , :cop' T; ecoi =
{ en • ' i~ ~.,o
^
x (3.6d)
Onx I i f V.=Ox
(3.6e)
;o ~ °o~t:p-*o(°)*o c°)1
A A A ^
(A,B,C,D), i.e.
45
(3.7)
A-IcA) . C(kIn-A)-IB+D
Then
(3.8)
ACA) = UECA)D~CA)
Combining Lemma 1.2 and Theorems 3.1 and 3.2 yields the results as follows.
or c o l u m n - r e d u c e d c a n o n i c a l ~matrix, respectively. •
Theorem 3.3 reveals the fact that the structural properties of nonsingular
canonical ~-matrices.
proved.
G(~) = NrCX)D~Icx), (G(k) = Dilck)N~(k)), where VCk) and N(A) are coprlme and
P~0of:
Let G(A) = N(X)D-I(x) be an irreducible, proper RI4FD, then from Theorem 3.4
G(X) - NrCX)D~I(X)+D c
From Eqs. (2.14) and (2.16)~ we can find the minimal realization quadruple
The latent roots of a nonsingular l-matrlx D(k) are defined as the roots of
det(V(l)) = O.
Lemma 3.5 [17-19] Let the set of latent roots of D(k) be O(D(A)), and the set
[o(D(X)[ - [o(A)[ = n •
Proof:
(3. lOs)
D(X) = Dr(X)Ur(X)
(3.10b)
D-I(x) = U:I(x)[Cr(XIn-A)-IB+Dr]
follows:
Definition 3.2 [17-19] Let X. be a latent root of D(X). A left Jordan chain
i
of left generalized latent vectors (left Jordan chain for abbreviation) f o r D(X)
(3.11a)
i ~1 (D(k) (Xi))Tpi(j -k) " 0mx 1
k=0
The scalar ~i is named the length of the Jordan chain [41]. Similarly, a right
0 % J ~ i - 1 such that
i 1 (DCk) (3.11c)
k=0 kT (Ai))qi(j-k) TM Omxl
vectors, respectively, n
c a n o n i c a l k - m a t r i x Dr(k ) are d e s c r i b e d as f o l l o w s .
Len~aa 3.6 Let (qij)~ be the ~-th component of q i j ' then for O~E~m we have
satisfies ~ = O, 0 ~ <~.
Froof~
D r ( A i ) q l 0 = DrhD(Ai)qio = Ore×1 •
Since Drh i s n o n s l n g u l a r , and i f KEffi0, the Eth column of Dr(A i ) i s e~ and ~th
row of Dr(k i ) i s (e~) T. Thus, the Eth component of qi0 i s z e r o . Also, from
Dr(ki)qil+D~l)(ki)qi 0 = 0m× I
We find that the Eth component of qil is zero. Repeatedly using Eq. (3.11c)
r e s u l t s i n ( q i j ) E = 0 i f KE = O, l ~ m . •
49
A
The relationships between the Jordan chains of A and Jordan cha£ns of Dr(X) can
C
be d e s c r i b e d in the f o l l o w i n g i m p o r t a n t theorems.
= ! 1 $(k)(ki) q . . . . . (3.13a)
qclj k 0 ~" r I~j-K;
where
Proof:
(3.13c)
ScVr(k) = ( k I n - A c ) $ r ( k )
Bc r
vCk)(1) = k$Ck-1)Ck)+Clln-Ac)~k)ck)
-r
, k~l (3.13d)
If qij' 0~J~i-l' is a right Jordan chain of Dr(k) , from Eqs. (3.11c) and
(3.13d) we have
B i ~-TD
1 (rk)(kl)qi," k - = i 1 $(rk-l)(~,i)qi,.k.]
c k=O ~J- ; k=l ~ [ ~J- )
+ ! 1 [(kiln_Ac)~(rk)
k 0~ (Xi)qi(J-k)]
(3.14a)
= Onx 1, l<J~£i-1
Thus, we have
= (3.14b)
(XiIn-Ac)~r(ki)qi0 Onx I
J 1 (k) j-1
(kiln-kc) k~O = Isjs¢.-I
= ~.~ ~r (kl)qi(;-k)
~ +k=[0 ~1¥ r .(k)(ki )qi(j-l-k) 0nxl' i
(3.14c)
(3.15a)
Define qcij i ~"
k=0 1 ~)(rk)(~i)qi(j_k)
= (3.1Sb)
(liIn-Ac)qci 0 0n× 1
51
(3.15c)
(liTn-Ac)qcij+qei(j.1) = Onx I , l~J~£i-I
Proof:
(qio)~ , %>o
(*~(O)qcio)~ 0 , ~=0
From Lemma 3.6 we have qio " ~(O)qci0" Thus, from Eq. (3.15a) we have
iJ 1 ~ k ) ( k i ) q i ( j _ k ) , j>o
k=O
Since
[ ~:(O)~r(O) if k=0
~(o)~Crk)(~i)=
I 0m if k>0
This implies
52
o if Kt-o
T (3.17)
P i j = BcPcij '
Proof;
From Eq. (3.13d)~ we have
[D(rk) (~)] T B T - k [ ~ ( r k - l ) ( x ) ] T + [ ¢ ( r k ) ( x ) ] T ( X I n - A c ) T .
J 1 . (k-l)(k)]T
~TI [ D(k)(~)]TB~pci(:rj - k) = k~ 1 ( k - l ) ! tYr i Pci(j-k)
k=0
53
÷
1 [~k)(~i)]T(XiIn-Ac)TPci(
k---f :
j- k)
k-O
= -T[
l r.,.CJ)cx i ),T
~r a Pci(-l) " 0mxl
J 1 ~(rk)(Xi)Dl (3.lS)
Pcij = Tp k=~O ~ hPi(j_k)
column [ a l i 2 , . . . , a l i ~ . , l ] T ;
l
[ a i j 2 . . . . . a i j K i , 0 ] T if Ki<%, and [ a l j 2 m . . . , a l j % , 0 , . . . , 0 ] T
i f ~.x >~..
j
Proof:
From Leuuna 3.7, DT
rh P i j ' O~J~i-l' are g e n e r a l i z e d l e f t l a t e n t v e c t o r s of a
l e f t Jordan chain of o(~) = D;~Dr(~). Define
It follows from Theorem 3.5 that Pcij' 0~J~i-I is a right Jordan chain of Ac
which is the associated companion form for DT(~), and Ac ~ (Ac)ij for l~i~m and
l~j.~,m is defined according to:
( c)ii = (Ac)IIcR ~
eigenvalue X.~ of A; qij for 0~j~.-l~ be a right Jordan chain of Dr(~) associated
with the same latent root X£ of Dr(~). Then qa£j and qij can be related by
= T-I i 1 ,(k)
qaij c k= 0 kT ~r (Xi)qi(j-k)' 0~J~£i-I
where ~r(k) (~) is defined in Eq. (3.13b), and Tc is the transformation matrix
which transforms A to its controller canonical form A c = T cAT -I.
c •
55
Proof:
and 3.6. In a similar way using Theorems 3.7 and 3.8 we can obtain:
i I .(k)(k )DT
Paij = T~Tp k=0 kT ~r i rhPi(j-k) ' 0~J~i-I
and
same is true for the Jordan chains of a A-matrlx. However, from Theorems 3.5
through 3.8, we observe that, given a left/right Jordan chain of a system map
matrix of a reachable pair (Ac,B c) can be uniquely determined, and vice versa.
Also, from Theorems 3.5 through 3.8, we observe that ~r(A) in Eq. (2o14C) links
system map in Jordan block canonical form, and M c and (M: 1) contain all the
right and left generalized elgenvectors of A c. Thus, from Theorems 3.6 and 3.7
Dr(k)" •
D:I(A) . ~ j ( A I n _ A j ) - I B j ÷ ~J , (3.19)
where BjT i n Eq. (3.19) contains all the left Jordan c h a i n s of Dr(A) , and ~ j
Proof:
- J
Hence, from L e n a 3.8, BjT and ~j contain all the left and right Jordan chains of
Dr()J. •
In an analogous fashion, we can easily obtain the following results for a
57
l a t e n t s t r u c t u r e of the row-reduced c a n o n i c a l k - m a t r l c e s :
B m
=T
where B~ i n Eq. 43.20) contains all the l e f t Jordan c h a i n s of D~(X), and Cj
A (3.21a)
and
(3.21b)
~-~ U; ki)qi(j_k ) ,
k-0 x
k - t h d e r i v a t i v e of U ~ l ( l ) e v a l u a t e d a t k = ~.
%.
58
Proof:
J (3.22)
1 I (v~k)(/i))Tpi(j k) 0mxl' x
k=0
D~k)(k) - [D(k)U:I(A)] k
k
kCsD(S)(k)u~(k-s)(k) (3.23a)
I
s=0
A kl (3.23b)
kCs = s l ( k - s ) !
(3.24b)
59
O~
°_
(3.24c)
Since u~l(ki ) is nonslngular for all Xi' from Eq. (3.24c) we conclude that
J 1
~
T
[D(k)(/i)] Pi(j-k) = 0mxl
05J5£i_l (3.25)
kffi0
By Definition 3.2, Pij for 05J$£i-i is a left Jordan chain of D(%) corresponding
to the latent root ~i"
Furthermore, since qij for 0KJ~£i-I is a right Jordan chain of Dr(%) , from
Definition 3.2 we have
J i k
) ~-~ [s=~O kCsD(Sl(kilUr(k-sl(ki)]qi(j-k ) 0mxl m
k=O
(3.27a)
Equation (3.27a) can be rearranged as
qij ~ i k-T
1 U; ( k ) ( X i ) q i ( j _ k ) for O~j~.-I
k=O x
Corollary 3.5 Let D(X) and Dr(l ) be defined in Theorem 3.11, Pij and qij for
A -- (3.28a)
Pij -
Pij ' o<j ~£.-I
I
and
are a left and a right Jordan chain of Dr(1) corresponding to the same latent
root %..
l
•
D(X) = UE(A)D~(X)
Pij and qij for O~j~;9~.-It be a left and a right Jordan chain of D~(X) with length
and
81
- O~j~t.-I (3.29b)
qij ~ qij ' x
(3.30a)
Pij =A i ~I [U(k) ()ki)]T~i(j_k ) , O~j~Ei-I
k=O
and
- (3.30b)
qij ffi qij' 0~j~E.-Ix
are a left and a right Jordan chain of DE(A) corresponding to the same latent
root A.. •
1
Theorems 3.11 and 3.12 with Corollaries 3.2 and 3.3 build up the latent
matrices to "matrix roots" which in the literature [17,18] are usually called
solvents.
any point interior to C, then the well-known Cauchy's integral theorem [42]
gives
62
f(;k) (3.31)
scalar function f(%) [43]. Let x ecmXm, the set of mxm complex matrices, and let
1 (3.32)
f(x) = ~ -~C f ( ~ ) ( ) ' I m - x ) - l d ~ = ~ ~C (~Im - x ) - l f ( ~ ' ) d k
described by
n n
fCX) = >" fkXn-k = ~ xn-kf k (3.33a)
k=O k=O
n n
f(x) ffi i fk xn-k = I xn-kf k (3.33b)
k=O k=O
where X 0 ~ I .
m
An n-th degree, m-th order nonsingular ~-matrix A(~)Ecmxm[~] can be written
as follows:
1 (3.35a)
AR(X) ffi~ ~C A()')(IJk-X)-Id~
where xEcmxm has all its eigenvalues interior to the simple closed contour C.
n
" k~o ~'E~i ~= :-kC~m~-x)-ld~l
n
>~ ~ x n-k (3.35b)
k=O
1 n
AL(X) " 2--~ ~C (Im%-X)-IA(%)d% = [ ~ i PC (Im%-X)-l~n-kd~]Ak
kffiO
n
>~ X'-kAk (3.36)
kffiO
Having these definitions of right (left) matrix polynomials in Eqs. (3.35) and
(3.36), we are ready to define the solvents of %-matrices.
AR(R ) =
AORn+AIRn-I+. " . + A n _ I R + A n = 0m (3.37)
64
Then R is a risht solvent of the A-matrlx A(A), where 0m is an m~n null matrix.
(3.38)
AL(L) = LnA0+Ln-IAI+...+LAn_I+A n - 0m
The solvents of A(A) can be constructed by using the latent roots and
have been extensively studied and r e p o r t e d in the literature (for example, [17-
vectors ql,...,q m Cleft latent vectors P I " " ' P m ) corresponding to latent roots
Al,...,~m,
A
then QAQ-I(p-IAp) is a right Cleft) solvent, where P = (Pl'''"
^ ;o )T
(Q = (ql,...,qm)) and A ffi diag(Alg...gkm). •
A set of right (left) solvents RkJ k ffi l,...,n (Lk, k - l,...,n) is called
n
of right (left) solvents has been investigated by Gohberg et. al. [20~21]~ and
and ~ denotes the null set. Similarly, the definition of a regular left solvent
is given by
A
oCL) O oCA(1)(1)) = (3.40a)
A(1) = ( I l - L ) A ( 1 ) ( l ) (3.40b)
m
defined as
n
(3.41b)
oCACk)) = v o(~.)
I%
k=l
given by
n
(3.42b)
oCA()0) = U o(L k)
k=l
is:
68
regular right (left) solvents of A(1), then the associated block Vandermonde
matrix
I I ... I
m m m
R1
2 2 R2 (3.43a)
V(R1,..., RI R2 " " " n
. . . . .
-1 n-1 R.-1
R2 " "" n
or
2 n-i
Im LI LI •.. LI
Im L2 L2
2 -.. L2-1
(3.43b)
vB(LI,...,Ln ) ...
1 L L2 ... L n-I
M n n n
shall see in Chapter IV, (Iml-R) is a right divisor of A(%), and A(1)(l) is a
left divisor of A(1). Similarly from Eq. (3.40b) if LEC mxm is a left solvent of
where A(1)(1)EcmXm[%] is a %-matrix with degree lower than that of A(%); and
67
A
Two examples are presented in this section. The first example demonstrates
Example 3.1
X3-3k-2 0 -1 1
3k+3 k2+2k+1 -2
vr(k) =
0 0 1
D;I(I) = Cc(kI4-Ac)-IBc+Dc
where
Ac
=
0o
2
" 0--0-0 i o
-3 -3 0 I-I
0,
3
oI , oI 0
~',o_
-!-I ; Be-t1
FIOO
l-o-
L o
0
o Ol
o
I
07
1.[
2.1
68
"
[iooo o] 0
0 0;0
0 0
0
; Dc"
[ooo]0
0 0
0 0
1
J
0.0 0.33333 -0.22222 -0.33333 0.22222
M~ " 0.0 -0.33333 0.55556 0.33333 0.44444
1.0 0.33333 -0.88889 -0.33333 -0.IIIII
-I.0 0.66667 1.22222 0.33333 -0.22222
Bj = H~IB = [ loilT 0 1
-1
0
;
-
~3
,
(: [0001
:)I 0 0 0
1 2 0 0 1
and {&l=3, ¢2~1, &3=l}, respectlvely. From Theorem 3.9 we have the l e f t Jordan
chain of Dr(~)
i 1 ,Ck)(~ )DT
P c l j ffi Tp k'T~r f rhPi(j-k)
k=0
where
69
- -o--o-I£ - 1 - 1
oolio]
(I)
(X)-
1
2X
o
o
0
0
o
~
(2)
(X)- [o o1
o
2
o
ooo
0 o
0
o
Thus,
Pcl2 = Tp[~r(XI)D hPl2 +~ hPll + ~ ~r urhPl01 " [ - I ' - I ' I ' I ' 0 ] T
If the left generalized elgenvectors Pcij are known, then the left generalized
T
Pij = BcPcij; 0<JK£i-l' l<iS3 .
-
[o 033 o11111]
I 0.33333 -0.88889 ; (q20) =
[0.,,,33]
-0.33333 ;
(ql0'q11'ql2) 0 0.0 0,0 0.0
r 0.IIIII-
(q3o) - i-o.~1~ii .
L o.o
I (2)(
qcl2 ~ ~r(Al)ql2+~(rl)(Al)qll+ 2 ~r A1)ql0 - [-0.IIIII,-0.22222,0.55556,
0.88889,t.22222] T
-
qc20 " $r(A2)q20 " [ 0 " 3 3 3 3 3 ' - 0 " 3 3 3 3 3 ' 0 " 3 3 3 3 3 ' - 0 " 3 3 3 3 3 ' 0 " 3 3 3 3 3 ] T
From Theorem 3.6, we can compute the right generalized latent vectors of Dr(A)
if the right generalized eigenvectors of A c are known. Using the calculated
generalized eigenvectors, i t can be easily verified that
T T T T T .T,
[Pcl2'Pcll'PclO'Pe20'Pe301 tqcl0,qcll,qcl2,qc20,qc30] = 15 •
Example 3.2
Consider a nonsingular A-matrix
D(A), [ A4+9~3+31A2+50A'33A2
7A+12 ~4+IOA3+38)L2"67~k+49
]A2.8~+18
which is not a canonical A-matrlx (check with Definition 3.1). Using column
equivalent transformation (Theorem 3.3), we have
D(A) ~ Dr(A)Ur(A)
where
or ['+72+1717+
11
71
and
t h e J o r d a n form m i n i m a l r e a l i z a t i o n of D-I(A):
r
D:I(A)= ~ j ( A I r - A j ) - I B j + D J
where
Aj = [1oo] -I 0 0
O -3 0
0 0 -2
Bj
= - I
-I
-I
1
~ .[o., -o.,~ o.~ -,.o]
3.0 -1.5 0.5 -3.0
~ : [ o o o]
{AI=-I, A2=-3, A3=-2} and {~I=2, ~2ffil,' ~3=1}, respectively. From Theorem 3.9,
we h a v e t h e l e f t Jordan chains of D (A):
r
1 -1
(PIO'Pll) ~ [0 1] ' (P20) . [_3] , (P30) = [ _~]
(qlO,qll) . [0.5
3.0 -0.75]
-1.5 '
..fqio~ = [ 0 . 2 5 1
L0-5 J ' (q30) = [-I.0]_3.0
0.5 -5.5
-qlo " u~ICxl~qio - [0 ' -~] [~.o] - [ ~.o]
~11 = UTI(Xl)qll+UT(1)qlo
- . olro.:51 : [-~.25
q2o " U:1(~'i)q20 ['~ -lJlO.5 J ]
-- -1.0 ]
q30 " UTl(Jt3)q30
". [7 -01][-3.0] " [ 3.0_1.0
left solvents :
I 1 I]
-1 0
ci
o
II I n | II II
g
| i . iI
' I' I 0
o Q o ! ! 0
i ,~I I I
I I ~=,
o 0
i I
I !
I
i I'
Q ~
Q~, ! u i !
• u
i i ° ~
i i I
! i
i i I I
0
l !
I i i I
m i !
• I ! I I I
l-m l - i
I II N .~I n I 8
I
O.I~" 0 w ~
It
0 ~
~J
I I ~-~
0
I ' I i'
~ O I 0
a I !
0
i I i I
a,
0
(D
E3
I.J
W
.I
¢D
lID
CHAPTER IV DIVISORS AND SPECTRAL FACTORS OF NONSINGULAR k-MATRICES
investigate the existence of the left/rlght divisors and the spectral factors of
matrix.
follows:
where L(%), R(k) ccmxm[%], then L(A) is a left divisor and R(A) is a right
divisor of A()k), end L()k), R()k) are left and right spectral factors of A()k). If
A()k). []
are dealt with in Section 4.2. Also, constructive proofs on the existence and
left/right divisors are provided in Sections 4.1 and 4.2. Section 4.3 is devoted
divisors are presented in Section 4.4. The main theorems used to implement the
[54,55] of the system map, which is the minimal realization quadruples of the
canonical divisors using the matrix sign algorithm [56-60] is presented. Some
4.1 Structure Theorems for Canonical Left Divisors and Complete Sets of
matrix of a reachable pair (A,B) in an MIMO system. The main tool used to study
defined, together with the structure theorem for complete sets of canonical left
78
divisors.
invariance subspace of X where ns<n , or AScS, and assume that V is the canonical
A L be the induced map in X/S by A, and B L ffi V'B, then B L - Im(B L) and the
A
B _X ;X
(4.2)
It is well known that V is epic [5], and if (A,B) is reachable, then the induced
n -I
pair (AL,B L) in X/S is also reachable, or X/S ffi <ALIBL> ~ BL+ALBL+...+AL v BL
n xn n xm
Theorem 4.1 Let (A,B) and (AL,BL) , where Ace n×n, BEe nxm, AL¢C v v, BLOC v ,
respectively, and n~nv>0. Then DLr(k) is a canonical left divisor of Dr(1) iff
Proof:
BL = VB (4.3b)
77
VBVr(k) = V(k~n-A)z~l,r(k )
Or
where
are defined as in Eqs. (2.14e), (2.11b) and (2.9a) for the reachable pair
Since the right-hand side of Eq. (4.5) is a k-matrix, DLr(k) is a column reduced
left divisor of D (k).
r
(2) Necessary Part
(4.6)
BLDr(~) = (%1n - A L ) T ~ L r ( ~ ) R ( ~ )
V
•
n xn
(4.7)
~Lr(X)R(X) ffiW~r(X) , WgC v •
Or
(4.ga)
Using the results in Eqs. (2.13), (2.14b) and (2.15b), Eq. (4.8a) becomes
(4.8b)
(ki n -AL)-IB L = T~W[(%In-A)TCI]-IB
V
A I nvxn
Define V = T~cWTcEC~ . Then Eq. (4.8b) can he expressed as
or
~ B L = VAXB, i>0
or
A l s o , s i n c e (AL,BL) i s r e a c h a b l e and
n -I n -I
[ B L , A L B L , . . . , A v BL] = V [ B , A B , . . . , A v B]
we have
n -I n -I
rank[BL,ALBL,...,AL v BL ] = nv~;min[rank(V),rank(B,AB,...,A v B)] .
n xn
Thus, rank(V)~n v . But vEC v so we o b t a i n rank(V) • nv, ie. V is epic. Thus,
an induced p a i r o f (A~B) i n X / S v l a V. •
follows.
KLi~< i , l~i~m
80
Proof:
A A
Let B = [ b l , . . . , b m], BL " [bL1,...,bi~]
(1) Let i=l.
Assume,~ by contradiction, that Kbl~''>~'.i"Then, we obtain {bLI,ALbLI,...,
-1 -i
ALL1 bLl} = {VbI,VABI,...,VA LI b I} which contains independent vectors.
However, Akbl-- for k><'1 is dependent of AJb I for j " I,...,K i. Thus,
vAkbl(=ALkbLl) for k>~ i is dependent of VAJbl(m~bLl) fOr jffil,...,~i. The
above result contradicts to the assumption, and therefore KLI<K1 .
(2) Let KLi~< i for iSik-l.
Assume, by contradiction, that ~K*ik>~ik" Then, we obtain
~LI-I
{bLI,~LbLI,...,
~ bL1,...,ULik,%bLik,...,~Lik-lhLikt -
KLI-I
{Vbl,VAbl,... ,VA bl,...,Vbik, VAbik ,...,vALik-lbik } which contains
independent vectors. However, for k>K..
Ik
• <-i
So that
vAkb ik = ~k I %~jVAJb~
~=I j-i
• ~,.CI
81
KLi~Ki, IKiKm •
Lemma 4.2 Let P(A,B) and P(AL,B L) be the reachability base matrices of (A,B)
and (AL,BL) , and (Kl,...~Km) ' let (KL1,...,gLm) be the Kroneeker indices of
(A,B) and (ALJBL) 9 r e s p e c t i v e l y and l e t V be the c a n o n i c a l p r o j e c t i o n shown in
Eq. ( 4 . 2 ) . Then
(4.11)
VP(A,B)H " P(AL, B L)
where
1 KLI O1+i OI+~L2 Om-I +I Om-l+KLm. nxnv
H = [ e n P . .. ~e n ,e n ,... pe n p.. .,e n p...pe n JEL
i
IKiKm
Oi = j~l= ~J '
Proof:
%-I %-1
Since VP(A,B) = [bLI,ALbL1;...,A L bLI,...,bLm,ALhLm,...,A L bLm] , from
Lemma 4.1 and the definition of H~ the result in Lemma 4.2 follows.
Theorem 4.2 Let V, H, P(A,B), P(AL,BL) be the same as those defined in Lemma
4.2. Then
(4.12)
V " P(AL,BL)HTF(A,B)-I
Proof:
The result follows from Lemma 4.2 and the fact that HTH ffi I
n
v
82
[arii2,...,ariiKi,l]
[Tcp]ij The reversed upper triangular Toeplitz matrix with first row
Proof:
Corollary 4.1 Let T c and TLC be the similarity transformation matrices which
-I (4.13)
V = TLcWTc
where
W ffi TL~pHTTep; T" I = TLcP(AL,B L) ; T c p = p-I(A,B)T~I
Lcp
Proof:
-I T
where N = TLCPH Tcp. •
From Eq. (4.5), the properties of the right divisor pairing with the
where DRh = D L hDrh is an upper triangular matrix with diagonal elements all
l's.
Proof:
Cl) R(A) = L
~ C ' V ~T : I ~r ( k ) + D ' D r ( ~ )
(2) Since both Dr(h) and DLr(X) are column-reduced canonical h-matrices, R(X)
is nonsingular. •
need to have a full set of divisors whose spectra cover the entire spectrum of a
this section.
84
The structure theorems for a complete set of canonical left divisors are
formulated in Theorems 4.3 and 4.4; while the irreducible complete set of left
AL 1
, Xl$1 :. X/$I
(4.14)
L/ .B =X A ~X
Proof:
follows.
4.2 Structure Theorems for Canonical Right Divisors and Complete Sets of
injection map s: S*×. Define ~ = A~and C R = C IS. Than the following diagram
86
commutes.
X A
(4.t6)
It £s known that S is monic, and if (A,C) is observable, then the embedded pair
usI Ker(CR~ -I)
(AR, C R) is also observable, or i~ 0.
Dual to Theorem 4.1, we h a v e t h e following results.
Theorem 4.5 Let (A,C) and (AR,CR) be observable pairs with left characteristic
s
are defined in Eqs. (2.17b), (2.20c) and (4.16), respectively. •
(~RI,VR2, .... ~Rp) be the observability indices of (AR, CR). Then ~RI~Vi, l~i~p. •
Lem~a 4.5 Let Q(AjC) and Q(AR, C R) be ~he observable base matrices of (A,C) and
(AR, CR) , respectively. Let S be the canonical injection map in Eq. (4.16).
Then
where
^ 1 ~RI TI+I TI+~R2 Tp-I+I T +V_ _ n Xn
H = [en,...,e n ,e n ,...,e n ,...,e n ,...,enP-I ~P]£EC s
i
~i = ~ vj , l<i<p
j-t
Theorem 4.6 Let S, H, ~(A,C), and ~(AR, C R) be the same as those defined in
TO = ~ - I ( A , C ) T ; ~ (4.19)
where
V,X~,
TOp = [Top]ijEC ~ J, ISiSp, l~j~p, V.V.>0
i]
[a~ii2,...,a£iiV .,I]
%
^ -I (4.20)
S = T0~rrR0
$2 = $2 CR2
x A X ~ V (4.21)
Proofs
= oCA) = oCD¢(1)) .
d i v i s o r s of D~(A). •
cage¢
n o n s i n g u l a r X-matrices.
Let (A,B,G,D) be a minimal r e a l i z a t i o n of A-l(%), where A(%)Ecn~m[%] is a
nonsingular h-matrix. Let XECn, U~Cm and YECm be the s t a t e , input, and output
n
spaces, respectively. SE C s is an A-invarlant subspace of X. (AL,B L) and
(AR, CR) are the induced and embedded p a i r s in X/S and S, v i a the epic map V and
the monic map S, o£ (A,B) and (A,C), r e s p e c t i v e l y . Then the following diagram
CO~Uteg.
~ "L ~ X/$
BL -, IV
(4.22)
s % ~-
CR
From the epic map V and the left inverse of the monic map S, we can
Lemma 4.7 Let S + ffi (S*S)-Is * and V + TM V*(VV*) -I be the left and right inverses
(4.23a)
Then
L o
Proof:
TT =
[ ]r "V]
S+S S+V+
VS
W+
= In s
L vs In_ n
S
-I
n
so t h a t T- I = [S,V +]
(b)
AT"
[,] V
A[S,V +1 -
V~ VAV+ J
s i . : e ALV - VA and S ~ " AS, we obtain S÷AS " AS; VAS - VSA R - 0 , rAY + " ~ , o.d
inverse of a nonsingular A-matrlx, A(A) , and the asociated (AEICR) and (AL,B L)
(4.24)
A(A) - DE(A)UA(A)DR(A)
where DL(A) and DR(A) a r e the r i g h t characteristic A-matrix of (AL,B L) and the
92
matrix.
Proof:
(4.25)
A(Z) = DL(k)R(Z) = L(k)DR(A)
From Lemma 4.7, there are no shared latent roots between DL(k ) and DR(1). Thus
Thus from Lemma 4.7 we observe that det(UA(k)) is a nonzero constant, and
column-reduced l-matrix, A(k), where AEC n×n, BEC nxm, CE Cain and DEC mxm. Let
n xn n xm
(AL,B L) be a reachable pair, where ALEC v v BL¢C v and nanv>O. Also, let
(A,B).
Proof:
Since A(A) = Dr(A)Ur(A) , from Theorem 4.1, the results of Theorem 4.10
follow.
(4.26)
A(A) = DLrCA)R(k)
where
Proof:
A A
or
A
BA(k) = (kln-A)T:l~c(k) .
Thus
reduced A - m a t r l x , A ( k ) , where AECuxn, BeCnap, CcCpxn and DeCpxp. Let (AR, CR) be
n xn pxn
an o b s e r v a b l e p a i r , where ARE C s s, CREC s, n>ns>O" A l s o , l e t DRE(A) be the
reduced X-matrix, A(k) and let T O be the transformation matrix which transforms
(4.27)
A(k) = L(k)DRA(k)
where
L(k) = ~0(X)T;lS~R+A(k)~R ,
Theorem 4.11 and Corollary 4.6 can be easily proved by using the same
[54,55] for the system map of the minimal realization quadruple for the inverse
and B,C c a n b e w r i t t e n as
T T T i (4.28b)
B = [BI,B2] ; Bi¢C ; i=1,2
mxni; (4.28c)
C = [Cl,C2]; CicC i=1,2
Proof:
96
theorem i s e s t a b l i s h e d as f o l l o w s :
n2xn
(I) Choose V ~ 1 0 , Z n 2 ! ; C n × n l a n d a p p l y Theorem 4 . I .
of a minimal r e a l i z a t i o n f o r t h e i n v e r s e of a n o n s l n g u l a r h - m a t r l x as f o l l o w i n g .
A = I AI 0 1 ; Aic Cn.xn.
x x ; i=1,2; A21cC
n2xn I (4.29a)
A21 A2
n.~
T TT
B = [B1,B2] ; Bi~C x ; i=1,2
(4o29b)
mxn i (4.29c)
C = [C1,C2] ~ CicC ; i=1,2
then
diagonallzed as
n.xn.
A = block diag[Ai,iffil,...,k]| Ale C I l| i=l,..., k (4.30a)
and
n. x m
TT (4.30b)
B = [BI,B2,...,B ]; BieC I ; £=i,..., k
(4.30c)
C = [Cl,C2,...,Ck] ; Ci¢ C 1 i=l,..., k
then
Proof:
A = block dlag[Aji,i=l,...,k ]
n,xn.
where AjIEC i x is a Jordan block associated with an elgenvalue ~i and a full
Jordan chain as
Aj i = • Xi I • (4.31a)
~.oXm
T T T ~ (4.31b)
B ffi [BjI,...,Bjk ] ; gjiCC ; iffil,...,k
mxn.
(4.31e)
C = [CjI,...,Cjk]; CjiEC l i=l,..., k
then,
Proofl
From Theorems 4.12 and 4.13, we observe that the canonical left/right
the system map A of the minimal realization quadruple (A,B,C,D) for the inverse
Let us d e f i n e the m a t r i x s i g n f u n c t i o n f i r s t .
Definition 4.4 The matrix sign function of AeC nxn with o(A)cC÷UC -, where C
and C- are the open right and left plane of C, respectively, is defined by [56]
where
and
100
1 (4.32c)
Sign-(A) = 2"~ ~r (Xln-A)-Id~
nlxn 1 P
where J+eC and J _ E C 2xn2 (n=nl+n2) are the collectlons of Jordan blocks
with o(J+)cC + and o(J_)cC-~ respectively. Then,
nlXn n2Xn
M-I ~ W " [ ]WIW 2 , wIEC , W2cC (4.35)
(4.36a)
Sign+(A) ffi ~I [Sign(A)+I n] = M[Inl ~ On2]M-1 " M1W1
101
and
These independent vectors are selected from the n(>nl) column vectors of
Proof:
S = Sign÷(A)Us
[e~ I k2 kn
where Us = ,e n j...,e n ] .
d A
Since Sign+(A) = MIWI, we have S = MIWIUs = MI~ , where ~ = WIUs. By
Sylvester's inequality [61], we o b t a i n
rank(W1)÷rankCUs)-n~rank(~)~minCrank(W1),rank(Us))
C4.38)
rank(~) = n1
102
S+ = ( S * s ) - t s * = [(MI~)*(MI~)]-I(MlC) . ~-IMI+
&
Theorem 4 . 1 5 D e f i n e S = I = ( S ) and l e t X b e t h e s t a t e s p a c e and A: X÷X. Then,
(1) S i s an A - i n v a r [ a n t s u b s p a c e of X.
Proof:
(4.39a)
AS = AMI~ = MIJ+~ = MI~(~-IJ+5) = SAR
where
A~ ~ ~-ta+~ (4.39b)
T h u s , S " Im(S) i s a n A - i n v a r i a n t s u b s p a c e o f X.
t o f i n d S from S ~ S i g n - l ( A ) :
103
Al~orithm 4.1
nxl
Given: - [SI,S2,...,Sn],Si¢ C , for l~i~n
Find: S - [Skl ,Sk2,...,Skn I], such that Ski¢{Si,l~i~n} for l~ki~n I are
independent.
Algorithm:
{Initialization}
{Froeessing}
For i=l to n Do
Begin
d:=P*Si;
If S.*d=0
%
Then
Begin
j:=j+l;{increment index j}
P:=P-d*d /Si*d;{Update P}
End{If}
rank of S}
Note: The * between two variables in the above algorithm is the matrix
transpose. •
Theorem 4.16 Let AcC mxn with Re(~i)Z0 for i=Ip2,...n, where {~i,i=l,2,...,n}
is O(A). Define
nXn-
S I = Ind[Sign+(A)] EC 1 (4.40a)
nxn^
S 2 = Ind[Sign-(A)] EC z (4.40b)
where Sign+(A) and Sign-(A) are defined in Eq. (4.42) and nl+n 2 - n.
Cnxn (4.41b)
Ms [SI,S 2] ¢
+ nlxn 1
AR 1 . SIASI¢ C (4.41c)
4- 4- + * -l * 4- *
S I and S 2 are the left inverses of S I and S2, or S 1 = (SIS I) S I and S2 = ($2S 2)
- I s *2' respectively~ where * designates conjugate transpose.
Proof:
Theorem 4.16 can be proved using Theorem 4.15 and Eq. (4.39).
105
l£ more than two blocks are needed for the block diagonalization of A, we
can construct the block modal matrix using generallzed sign matrix functions.
Definition 4.5 Let AcC nxn and Re(O(A))n{rl,r2} = @, where rl<r 2 and rl,r2C~
The ~eneralized matrix sign function of A with respect to the open interval
Sign(rl,r2 ) = 2 S£gn;rl,r2)(A)-I n
- (4.42a)
= I n - 2 Szgn(rl,r2
where
. cA1 . (4.42b)
Szgn(rl9 r2) (A) ~[ S l g n ( r l ) ( A ) - S i g n ( r 2 ) (A) ]
• - A (4.42c)
Stgn( rl, r2 ) (A) = In-S ign;r I' r2 ) (A)
and
(4.42d)
Sign(ri)(A) ~ Sign(A-riln) , for i=1,2 .
rl
nxn.
+ 1 (4.43)
S i = Znd(Sign(ri-l'ri)(A))eC , l~i~k+l
A A
where r0 = -=and rk+ 1 = =. Assume that n.~0 for l~i~k+l.
i
Then
Ms ffi [ S l , S 2 J . . . , S k + I ] (4.44b)
and ARi , l ~ i ~ k + l , a r e d e f i n e d as
(4.44c)
ARi = S+AS.
~
Directly from Theorem 4.15 and Definition 4.5. From Definition 4.5 and
Theorems 4.13 and 4.17, together with Corollary 4.9, allow us to decompose
a nonsingular A-matrix into a complete set of canonical right divisors such that
each divisor has latent roots clustering inside a certain vertical strip of the
complex plane C.
To obtain the canonical projection map from matrix sign function, we define
V = {Ind[(Sign+(A))T]} T (4.45)
nlxn 1 ÷
Lemma 4 . 9 There e x i s t s a n o n s i n g u l a r m a t r i x ~EC s u c h t h a t V -- r]W1 and V
107
4.
= W;q -1, where W1 iS defined in Eq. (4.35); V+ &= V * ( W * ) " 1 end W1 &- W~(W1W;)-I
are the right inverse o f V and WI, r e s p e c t i v e l y .
Proof:
From t h e d e f i n i t i o n o f V we c a n d e s c r i b e V as
k
T A k l k2 nl
where Uv ffi [ e n , e n , . . . , e n ].
* * * -1 + -1 (4.47)
V + = V*(VV*) -I = (qWl) (nWiWln) = win
A
Theorem 4.18 Define S = Ker(V); let X be the state space and A: X÷X. Then
Obviously, from Theorems 4.1 and 4.18, V defined in Eq. (4.45) can be used
as the canonical projection map for finding the canonical right divlsors of
nonsingular l-matrices.
Paralleling Theorem 4.17, we can construct the block modal matrix from
(4.49a)
MVA~I - ~lock diag[ALI,AL2,...,ALCk+I)]
where 4 1 i s a b l o c k modal m a t r i x
T T T T (4.49b)
MV = [ V I , V 2 , . . . ,Vk+ I]
ALl , l ~ i ~ k + l , a r e d e f i n e d as f o l l o w i n g
A L i = Vi AV~ (4.49e)
Proof~
Theorem 4.20 Let AEC nxn and Re(~(A))fl{rl,r2} ffi ~, where rl<r 2 and tiER for
i=I,2. Define
109
+ nXn 1 (4.50a)
S(rl,r2) = Ind(Sign(rl,r2)(A))EC
and
= [Ind((Sign~rl,r2)CA))T)lTcCnxn2 (4.50b)
V ( r l , r 2)
T~ - (4.51)
V ( r l , r 2)
+
where S ( r l , r 2 ) is the l e f t i n v e r s e of S ( r l , r 2 ). Then
(4.52a)
0n2Xn I AL
where
AR ~ S+ (4.52b)
(rl,r2)AS(rltr2)
~L"A~C~, r2>A~r l, r2> <4.52=~
A
Proof:
- . +
V(rl,r2 ) = V~(In-Slgn(rl,r2)(A))
where
E1 k2 n2]T
U~ = [e n , e n , . . . , e n
S(rl,r2 ) = Sign[rl,r2)(A)U s
where
k
k I k2
U S s [en ' e n , ' ' * ' e n hI]
s - • - • - -
rl,r2) (rl,r2) In-S1gn(rl,r2)(A))S1gn(rl,r2)(A)Us 0
since V(rl,r2 ) end S(rl,r2 ) are of full rank. This implles that KerC~Crl ' r2 ) =
-
described as f o l l o w s :
and
of Sign(A). Recently, a new algorithm to compute the matrix sign functions has
continued f r a c t i o n :
The recursive algorithm for computing the matrix sign function of A with
Re[O(A)] ~ 0 is
(4.55)
Slgn(nk+l)(A) - SigU(fk) [Sign(nk) (A)]
The degree of the convergence rate of the algorithm in Eq. (4.55) at the
to [60] for further details of the new matrix sign algorithm in Eq. (4.55).
block matrix along the diagonal has eigenvalues clustering inside a certain
vertical strip of the complex plane C. These algorithms are suitable for
triangularize a system map A such that each block matrix along the diagonal has
systems, the argument matrix A in the matrix sign functions9 Sign(At, Sign+(A)
Sign(rl,r2)(A) = 2 Sign(rl,r2)(A)-I n
• - (4.56a)
= In-2 Sxgn(rlpr2)(A)
where
• + 1 (4.56b)
S1gn(rl,r2)(A) ~ ~[Sign(rl)(A)-Sign(r2)(A)]
S1gn(rl,r2)(A ) ~ in_Slgn(rl,r2)(A
. -
• * ) (4.56c)
matrix sign function of A with respect to the open annulus centered at origin of
divisors and the spectral factors for a nonslngular k-matrix via the Jordan form
-2 2 -5 4
A = 2 -5 I0 -23
I -2 4 -4
0 0 0 2
and 1 3 2
0 -2 -2
B= 0 1 1
2 5 3
1 2 1
X2iX-2 -18X-18 11
Dr(X) = X3+3X2+3X+l -1
0 1
-1 1 0 0
Aj = 0 -I 0 0
0 0 -I 0
0 0 0 2
0 1
Bj = 0 1
II
12
Cj -- -I 0 -1
0 0 0
Dj = [ioo] 0
0
0
1
Aj " Bj =
0 A2 B2
where
[
A1 = 00~ I0° -I
] ,~ =
[:o]
2
BI = 0 i 1 B2 =
0 1 1 2 1
1+i -I 1 ]
Drl(A) = 0 12+21+1 -I
0 0 1
and
I+4 6 I ]
Dr2(1) = L-3 1-50 -Ii
Cj = [CI,C2]
where
CI =
[6 513]
I -I 0 C2 =
[i,1311]
-i 0
0 0 0 0 0
12-41+13 108 O ]
D.l(k) = 25/3 k+70
0 0 1
and
115
D~(~) = [ 0 -17
`1-2
0
,t+1 ~1
0
[12:] [011]
~1 = ~2 = 1 x 1
0 1 I 2 1
E1 =
[i ~] -i
0
~2 = -1
0
o
0
I A2-X-2 -2X-2 i 1
Dr2(A) = 0 ~+I -i
0 0 1
D~I(A) =
I ~-2
0.5
0
-18
~+4
0
01
0
1
Dr(),) - D r 2 C ) , ) U r t ( ) , ) D t l ( ) , )
w h e r e U r t ( ) , ) i s u n l m o d u l a r and c a n b e f o u n d as
Urt(),) - -0 5
20]
),-2 0
o 1
Vr(),) c a n a l s o b e f a c t o r e d as
where
1 2),+8 01
V£1(),) -. Ur£C),)Dtl(),) - 0 ),2+2),+1 0
0 0 I
Example 4 . 2
I
-1.0 -2.0 3.0 -11.0 24.0
A " 2.0 -4.0 10.0 -18.0 37.0
!.0 -2.0 4.0 - 7.0 12.0
0.O 0.0 -0.5 1.0 -3.0
and
-2.0
1.0 1.0
-2.0 0.01
-3.0
B= 1.0 1.0 3.0
2.0 1.0 4.0
1.0 0.0 1.0
The r i g h t characteristic ) ` - m a t r i x i s f o u n d t o be
1
),3+3.6667),2-0.6111),-4.0556 ),2-0. 3333),-0. 6667 -1.0
Dr(),) ,-
I -0.5556),-9.1667
0.0
),2-0 . 6667),-2 . 0000
0.0
-2.0
1.0
117
Sign(ro)(A) = 15
O' ne V2
V3
DD = D r
DE2(A), D~3CA)}:
and
r x+o. 16667 0.20000 0.0 "1
Dr3(),) = | -6.80556 X-2.16667 0.0
L o.o o.o 1.o J
To p e r f o r a the s p e c t r a l f a c t o r i z a t l o n ~ u s l n 8 Theorem 4 . 2 0 , we have
[L S~2'r3) J1 !--0"--0~0-00
[" 1"°°°°° -2"°°°°°
-0. lee67 4"°°°°° 2"°°°°°
0.33333 o.le6e7 0.00000t
0"00000
T " -- / 0 . 0 0 0 0 0 2.00000 -5.00000 12.00000 -2,,.00000
/~(~2,~3)/ /o.ooooo l.ooooo o.ooooo 1.ooooo-2.00000
LO.OOOOO o.ooooo o.ooooo o.ooooo 2.00000
trlangularized:
% - TA~ - 1 -
03x 2 AL
- - 770
0.0
0.0
0.0 ~--g.~(Og
0.0 I-1.6000
0.0 ~ 0.2000
-l.6~
-0.8000
0.4000
O.50OO
1.0000
-1.0000
]
B T - TB -
[ ] ,o.ooo .
F 4,0000 11.0000 26.0000 1
0.83333 2.1667 l
i --_5 .-6o-b6- - 3 : o ~ o ~ - ~ 7o~o-o-]
BL | 1. 0000 -1.0000 -1.0000 |
L 2.00O0 0.00000 2.0000 J
DT = D
I
~2+3.5X+4.25 0.8X+0.7 -1.01
Dr2(X) = 6.25 X+l.5 -2.0
0.00 1.0 1.0
and UrE(A) i s a u n i m o d u l a r X - m a t r l x :
UrE(X) = z 3
=D
L r
functions
Sign(ro)(A L) " 13
Slgncr2)CAL) " -Z 3
TL -
[...i,..,] El.O _..o o.oj
- o.o o.o l.O
VL(rl,r2) 0-.8 - "0~4" -0"~0
DLT " Dr
~r2(k) = ~r2(k)Or~(k)~1(k)
Dr(h) =
[I+4.oo 1.oo 1.oo][1.o
0.00 1.00 0.00| 0.0
o.o
1.0
o.o]
-2.0
0.0o 0.o0 1.00J 0.0 0.0 1.0
I 6.25
o.oo
1+1.50
o.oo
0.00| |-6.8056
1.ooj L o.ooo
I-2.1667
o.oooo
0.0
1.o J
CRAPTER V FEEDBACK CONTROL OF MULTIVAEIABLE SYSTEMS
the most successful design methods are state-feedback control in the time domain
[5,63-71], and the inverse Nyquist array and the root loci design in the
matrices and their divisors for constructing the state-feedback control laws.
These new approaches permit a deeper insight into some structural aspects of
Therefore, advantages of both the time domain and the frequency domain d e s i g n
Section 5.3 is devoted to deriving the latent structure assignments which take
~x(t) = AX(t)+BuCt)
(5.1a)
y(t) = CX(t)+Du(t)
(5.1b)
124
where X(t)~C n, u ( t ) ¢ C m, and y(t)cC p are state, input, and output vectors,
differential operator with tER. For discrete time systems, k is a forward shift
operator with tEZ. Assume that the system in Eq. (5.1) is reachable. The
(5.2)
u(t) = -FX(t)+GFr(t)
where F, GF, and r(t) are the feedback gain, input gain, and reference input~
a state-feedback control i s as f o l l o w s .
Lense 5.1 Let the state-feedback control law be Eq. (5.2), and let Drc(~) be
the open-loop system, Ki,i=l,...,m are the Kronecker indices, and ~ci[Drc(k)]
Proofz
R e f e r t o Popov [ 3 5 ] . •
matrix. However, this structural limitation does not restrict the freedom of
realization -1
of Drc(~) , where OF = Dr hD-~
r c ,D rug
~_ i s the l e a d i n g - c o l u m n c o e f f i c i e n t
Proof:
Fc
FT:I " DrhF c
- (5.4)
we o b t a i n
Tc ( A - B F ) T c- I = Ac-BcFc = A c - E b c F c ~ c
(5.5a)
(5.5b)
TcBGF = BcDrhD;Ic = EbcD:hlc A B c
(5.6)
D;~CX) " ~T(0lCklr u -~ c )-Igc +~ c = C(XIn-CA-BFI)-IBGF+DGF
realization of D~(~). •
Lemma 5.2 reveals the state space structure of the closed-loop right
we have the following results on the structure of the feedback gain when some of
for some l~i~m, then the ith rows of Fc defined in Eq. (5.4), can be chosen to
be null, and therefore the ith rows of F and F are also null.
e
Proof:
From Eq. (5.5a), we have Ac = Ac-EbcFc" If <.=0,~ the ith column of Ebc is
and F
C
are both null. •
In the following sections of this chapter, we d i s c u s s various schemes f o r
matrix of the closed-loop canonical EMFD is to be the same as that of the open
127
h a v i n g t h e same K r o n e c k e r i n d i c e s a s t h o s e o f t h e o p e n l o o p r i g h t c h a r a c t e r i s t i c A-
The important feature of this scheme is that the numerator A-matrix of the
state-feedback controls.
Drc(A) , and Drh c is the leading-column matrix of Drc(A) ~ then the feedback gain
and input gain become
F = F T (5.7a)
e c
-I (5.7b)
GF = DrhDrh c
where Fc can be d e t e r m i n e d by
(5.7c)
FcSr(A) - GFDrc(A)-Dr(A)
Proof:
BcDrc(A) = ( A I n - A e ) ~ r ( A ) = ( A I n - A e ÷ E b c F c ) ~ r ( A )
= ( X l n - A c ) ~ r ( A ) + E b c F c ~ r ( A ) = BcDr(A)+EhcFc~r(A)
128
Therefore,
Bc(GFDrc(A)-Dr(k)) = EbcFc~r(k)
D:~(GFDrc(k)-Dr(A)) - FcSr(k).
Thus, the result of Eq. (5.7c) follows.
Proof:
= CcSr(X)D;cI(A) = N r(k)D:~(k) •
Theorem 5.2 Let D = 0 and the state-feedback control law be shown i n Eq.
pxm
(5.2) with the ith row of F being null if K.-0 and G F being nonsfngular. Then,
I
the RMFD of the closed-loop system can be expressed as
129
(5.9b)
; r c (~) = G ; I [ F T : I ~ ( ~ ) ' D r (~)]
Proof:
where
;:Ioc - , C0 C nZc I o, c
BcDrcCX~ = C~In-Ac)~rCX)
whence
A
= BcDr(~)+BcFT:I~r(~)
which
130
Remark
Theorem 5.2 reveals the fact that the closed-loop RHFD can be written as
Nr(A);:~(A) , where the EMFD may not be canonical even if Nr(l) is the numerator
From Theorems 5.1 and 5.2p we observe that the non-canonical column-reduced
(KIj...jKm) , Which is the set of Kronecker indices of the open loop system in
^ ^
(5.10a)
F=FT
C C
A
(5.lOb)
G F = DrhDrhlc
Proof:
Example 5 . 1
Ax(t) = Ax(t)+Su(t)
y(t) = CX(t)+Du(t)
where
I
-2.0
-I.0 1.0
-2.0 0.0 1
-5.0
B = 3.0 -I.0 1.0
-2.0 3.0 4.0
1.0 -2.0 -3.0
C = [3.65.85.04.0 0.0]
5.42.25.03.0-1.0
D = 02x 3
are found to be E 1 -- 3, <2 = 2 and ~3 " 0. Note that the input matrix B
contains a dependent vector, but the system is reachable. From Eq. (2.16)~ we
132
O(k) = Nr(k)D;l(k)
where
-4k2+13k-13 7k-12 0 ]
SrCk) =
-5kX+10k+3 12k+7 0
The open loop poles are -i,-I,-I,I,2. We shall construct linear state-feedback
control laws using both the right characteristic k-matrix assignment and the
I -I 0]
Drh c = 0 I 0
0 0 1
u(t) = -FX(t)+GFr(t)
Gc(~.) = Nr(A)D:~(}O
A3+4k2+8A+6 -k2-7~-12 0 1
rc(),) = 2~2+4~+2 k2+7k+12 0
0 0 1
Drhc 0
- ^ [ X3+6X2+I2X+8 0 0 ]
Dr~cDrc(k) ffi 2k2+4A+2 k2+7A+12 0
0 0 1
134
and t h e i n p u t g a i n
.(t) = -FX(t)+GFr(t)
closed-loop system:
ccCx) = .rcCx)D~<x)
where
-4X2-X+ll
7X-12 0 1
Nrc(X) = -SA2-14X-11 12~+7 0
Nrc(k) = Nr(k)Ur(k)
Drc(~) - Vrc(XIUr(A1
where Ur(X I i s a u n i m o d u l a r k - m a t r i x
Ur(~) - l i0 0]I
0
0
1
t h e n u m e r a t o r o f t h e open l o o p c a n o n i c a l RHFD w i t h N r c ( ~ ) ) .
{[i,i=l,...,k}, and l e f t g e n e r a l i z e d l a t e n t v e c t o r s , { P i j , j = 0 , 1 , . . . ,
~.-1,i=l, .,k} where l i is the length of Jordan chain correspondlng to ~i and
kI °*
i~ l £i=n. Define
Aj ~ block diag{~i,i=lt...~k}cCnxn (5.11a)
where
i ~ .xZ. BJ1
1 [-]
iT.-1
.xm
1 Xl 1 e ~ ~; ~j ~ e ~m;- £
L
1 Bj i
pT
i Bjk io
(5.rib)
Assume that (AjmBj) is a reachable pair and the Kronecker indices of (Aj,Bj) are
equal to those of (A,B). Then, the desired feedback galu F and the input gain
characteristic A-matrix of (Aj,Bj) with the leadinE-column matrix Drh c and Dr(k)
Drh-
Proof:
where
137
[ a i i 2 , • • . , a i i ~ . , 1] T
1
[aij2,...,alj~i,0] T i f Ki<~j or [ a l j 2 , . . . , a l j K j , O , . . . , 0 ] T if
<.>K..
z J
Proof:
(5.12)
Fc~r(~O = OFDrc(~)-Dr(~)
(A,B).
Proof:
(s.x3)
form A0 = TojAjToj-
--I -
5.5 c o n t a i n s the p r e a s s i g n e d r i g h t l a t e n t s t r u c t u r e . •
can be determined by
139
Proof:
Example 5.2
Let us consider the same open loop system described i n Example 5.1. We
shall construct the linear state feedback control laws using left and right
Let the desired latent roots and left latent vectors of the closed-loop
Xl =-5 ; ~to = It -2 o] T
X2 = -4 ; P20 = [3 -40] T
i3 = - 3 ; ~3o = [t-lo] z
~4 =-2 ; ~4o = [1 1 o] T
X5 = -2 | P50 = [2 -3 O]
T
-i -40 00 00 00 1
Aj = o -3 o o
0 0 -2 0
0 0 0 -1
and 1 -2 0 "]
3 -4 0
Bj = i
1
2
-I
1
-3
0
0
0 J
140
to b e
I -1 -1 ]
OF " 0 1 -2
0 0 1
GcCk) = N r ( X ) D : ~ ( k )
X3+7X2÷22k+24 4k+8 0 1
Drc(k) = 6k÷12 k2+gA+14 0
0 0 1
Let the desired latent roots and right latent vectors of the closed-loop
~2=-4 ! ~2o" [ 1 2 ° ] T
o o0 o0 o]
-4 0
Aj = 0 -3 0 0
0 0 -2 0
0 0 0 -
and
1 0 1 0 1]
~j = o 1 o I -2
0 0 0 0 0
2, V 3 = 0 which are the Kronecker indices o£ the open loop system. Choosing the
matrix a s
Drhc
= [11i] 0
0
1
0
and t h e i n p u t g a i n i s
142
I I -2 -I 1
GF = 0 I -2
0 0 I
u(t) = -FX(t)+GFr(t)
Gc(k) = Nr()`)D:c1()`)
)`3+10A2+31~,+30 2)`2+6.5)`+9 0 1
Drc ()`) =
I 0
0
),2+6A+8
0
0
1
It is straightforward to check that Drc()`) has latent roots ~i and right latent
vectors qi0 for I<i~5 as desired.
Ac = Block d i a g [ A c i , i = l , . . . , k ] (5.15a)
and
= -T -T (5.15b)
Bc [Bcl'''''Bck]
If the Kronecker indices of (Ac~Bc) are the same as those of the open loop
system, then the linear state-feedback control law in gq. (5.2) can be
c o n s t r u c t e d as
F = F T (q |f'
"~'-'--a"
C C
where F is determined by
C
(5.16c)
Fc~r(X) ffiGFDrc(X)-Dr(X)
transforms the system map A of the open loop system to i t s controller canonical
144
form A •
C
Proof=
Define
and
- - (5.17b)
T 0 = [c01,...,C0k]
If the observability indices of (A0,C O) are the same as the Kronecker indices of
the open loop system and the leading-row coefficient matrix of the left
characteristic l-matrix of. (A0~C O) is Im9 then the linear state-feedback control
F = F T (5.18a)
c c
-I (5.18b)
G F = Drh Drh c
where Fc is determined by
(5.18c)
Zc~r(k) = CFDrc(k)-Dr(X)
And, Drh c is an upper triangular matrix with diagonal elements all l's; Drc(k) =
145
Drh c D~c(X) and D~c(X) is the left characteristic h-matrix of (A0,C0). Dr(h) is
the right characteristic h-matrix of the open loop system. ~r(X) is defined by
Eq. (2.13b). T c is the transformation to transform the system map A of the open
Proof:
Similar to t h e p r o o f o f Theorem 5 . 5 .
matrix such that all the necessary conditions in Theorems 5.6 and 5.7 are
assignment.
k
Proposition 5.1 L e t Vrci(X) = diag{dij(X),j=l,...,m } for i=l,...,k, and
it
deg[dij(X)] = ~j, then {Drci(~),i=l,...,k} is a complete set of left divisors
which can be a s s i g n e d to c o n s t r u c t
D (X). •
k rc
S i n c e Drc(X ) = d l a g { w dij(X),j=l,...,m} , the Kronecker indices of Drc(X)
k
Proposition 5.2 Let D&ci(X) = dlag{dij(X),j=l,...,m ) f o r i=l,...,k, and
i=l
deg[dij(X)] - Kj, then {Dg.ci(X),i=l,...,k } is a complete set of right divisors
Selecting a Drh c and applying Theorem 5.7 gives the desired feedback gain F and
(5.19)
Drci(~) = diag[li_l,di(~),Im_i],i=l,...,m
(5.20)
D£ci(l) = diag[li_l,di(1),Im_i],i=l,...,m
m
Gc(~) = C[~In-(A+BF)]-IBG F = ~ N .(~)D-~(~)
i 1 rcl r~
m s. (~,)
(5.21)
y(t) = I~ ri(t)
i=l
same reasoning can be applied for Proposition 5.4 with Drh c = Im. The technique
which is d e f i n e d as
n.(X)
I (5.22)
Yi(t) = ~ ri(t) for i'l,...,m and m=p
i
The necessary and sufficient conditions for input-output decoupling using linear
for m=pj where ~G(%) is the interactor of the open loop systemt namely
f°
~G(~.) = diag{A z,i--l,...,m} for fik0
as
(5.23a)
Drc(X) - Dd(X)N(X)
(5.23b)
Nr (~) = NdC~)N(~)
deg[di(A)] = deg[ni(k)]+fl
For practical applications~ N(k) must be stable and invertible. Note that the
m
closed-loop poles associated with w di(A) can be arbltrarily assigned, and
i-i
that ni(k) is the monic common divisor of the ith row component of Nr(A). The
Proposition 5.5 Assume that di(l) are monic, then the feedback gains for the
(5.24)
Fc~r(A) " D:hI[K;IDrc(k)-D r(A)]
Proof:
Since
aim ~
l.~o ~G(1)Nr(),)D (1) = KGG F
and
Nr(k)Dr~(k) = diag{ni(1)/di(1),i-l,...,m }
Example 5.3
~(t) = D3~(t)+Bu(t)
where
I
0
-2
0
4
0
-
[i I
0
4
I
1 -1 4 -5 11 ]
C = 0 1 -1 4 -5 ; D = 02x 3
input matrix B contains a dependent vector. From Eq. (2.11b) we obtain the
Fo 1o o oI 0 0 0
/ o o 1 o o I 0 0 0
Ac = TcAT:I = i _ 2 _ 3 _ _0_ _0__ _0 3 ; B - z B - I 0 1
/ o o o o 1/ c = 0 0 0
L-3 -3 0 -1 -2 j 0 1 2
Cc CT:I [2 0 2:4 3]
= = 4 4 3 !2 2 ; Dc = 02x3
I I ~. ~k
2 o OIT 0 o 0]
Sr(X) = 0 0 0 I Ar - -3 0 -1 -2 ;
0 0 0 0 0 0 0 0 0
Drh
= 1
0
2
1
-1
=
[o] 0
0
I
0
-2
1
150
• r ( k
plays an important role in determining the eigenvectors of the system )
matrix Ac and the latent vectors of t h e right characteristic )`-matrix of the
Dr()`) = Drh[Dh(k)-Ar~r()`)] =
[lo_1][ 3i3 _2 o oj
0 1 -2 3 +3 )`2+2)`+I 0
0 0 1 0 1
=
[ )`3-3)`-2
3).+3
0
0
),2+2)`+1
0
-1 ]
-2
1
[ 2)`2+2 3)`+4 0 ]
~r()`) = Cc~r()`) =
3)`2+4)`+4 2),+2 0
CjC
[o.33333o.o1.o00o7,_o.5
0.0
oo]0..5 0.0 |l 0.0 -0.5 ; DjC =
[o o o]
0 0 0
0.0 0.0 0.0 I 0.0 0.0 0 0 l
151
[ 13+7~i+141+8 0 0
Vrc(l) = 12+6~+8 0
0 1
where
0.0 1.0 - 2.0 5.0 -18.0 I
F l
= -3.0 6.3333 -14.66667 24.66667 -45.66667
FeTe 0.0 0.0 0.0 0.0 0.0 J
=
_ [,0
GF DrhDr~ h =
0 0I -
312+41+4 21+2 0 0 0
The partial fraction expansion of the closed-loop RMFD yields (refer to Section
6.1)
Example 5.4
To demonstrate the input-output decoupling of a HIHO system~ we consider
the following 2-input, 2-output system:
IxCt) = AX(t) + Bu(t)
where
-i -2
-2 24 -5
-9 421 - 1
A = 2 -5 10 -23 ; B= - 3
1 -2 4 -4 1
L0 o o o 2 0
C =[ : - 5 , 7 8
10,-21, 8 8
41 ~. ]; D
- 9 -20 43 -85 = 02x2
The interactoris
~o(I> " [ 0
KG 1~
[,
0
Dd(~) =
(X+l)(X+2)
0
o
(~+2) (X+4)
]
Thus, the c l o s e d - l o o p right characteristic X - m a t r i x becomes
Dre(X) = Dd(X)N(X )
= [(~+1)(x+2)(x+4)o -(~+1)(~+z)(x+2)(~+4)
]
u(t) = OFt(t) + ~ x ( t )
where
o~-,;'-[o ]:I
r-3.0 8.66667 -19.33333 37.33333 -80.33333 ]
F = FcT c ffi L _ 3 . 0 6.33333 -14.66667 24.66667 -45.66667
Cc(~) =
(~.÷I ) (X+2)
0
0
X-2
(~,+2)(X+4)
]
showing t h a t the d e s i g n e d system i s i n p u t - o u t p u t decoupled.
CHAPTER VI STRUCTURAL DECOMPOSITION THEORIES AND APPLICATIONS IN
MULTIVARIABLE CONTROL SYSTEMS
tools for studying the properties of the basic interconnection structures~ these
being the parallel formsp the cascade forms~ and the semi-cascade forms of
multivarlable systems.
fundamental structures:
Gr(A) ~ N r (A) D: 1 (k) - (NOA3+N 1 A2+N2A+N 3 ) (ImA3+DIA2+D2 A+D 3 )-1 (6. Ib)
(6.1f)
^
where NO, N i and Di, i=1,2,3 in Eq. (6.lb) are matrix coefficients| S i and Si'
i=i~2,3 in Eq. (6.1c) are the spectral factors of Nr(A) and Dr(A) , respectively~
-Li, i=1,2~3 in Eq. (6.1e) are the complete set of left solvents [17-18,25-27]
of Dr(k). The matrices Ki, i-I,2,3 in Eq. (6ole) are the block residues [27-31]
155
another but not from the dynamic modes of the numerator. The seml-cascade form
I n t h i s c h a p t e r , we s h a l l i n v e s t i g a t e the i n t e r c o n u e c t l o n s t r u c t u r e s of an
here for c o n v e n i e n c e :
becomes ( r e f e r to S e c t i o n 2.2)
G(I) = Nr(k)D;I(A)+D
G(A) = D~I(A)NE(A)+D
where D[I(A)NE(A ) i s a strictly proper rational A-matrlx, and N£(k) and DE(A)
Dr(A) and DE(A) may not be monic. Instead, they are c a n o n i c a l column-reduced
approaches [5,81-82].
Circuits
n i Xnoi (6.3a)
Acp ~ b l o c k d i a g [ A c i , i f l p . . . ~ k ] ; Ac£¢C
n. Xm
T T T T cC~ (6.3b)
Bcp [ B c l ' B c 2 , ' ' ' , B c k ] ; Bci
A (6.3d)
v = D
cp
(A,B,C,D). 0
he decomposed as
(6.4a)
i=I cp
where
SriCl)V~(X) ~ cci<IIn-Aci)-IScl
i
158
(6.4b)
where
I Kil-I
1 ~ ...
00... ... 00 ... *I
~ri(X) A o o . . .
1 k ... kK ... 0 0 ...
• • tea • .•a • •a• • at• .~
(6.4c)
Proof:
k k
= I C c i ( k l n - A c i ) - I B c i + D c p = .~ Nri(k)D:~(~k)÷Dcp
i=1 I I 1
Equation (6.4) shows that the RMFD can be described by partial fraction
expansion forms. The interpretation of the system structure in Eq. (6.4) is that
the original system has been decomposed into the parallel connections of
The dual results of the parallel decomposition for observable systems are
d e s c r i b e d as f o l l o w s .
n. Xm
A T T T T, L (6.5b)
BOp = [B01,B02,...,B0k] " BOiE:C
pXn. (6.5c)
COp [CollC02,...,Cok]; C0£¢C 1
Dop ~ D
system (A,B,C,D). []
(A,B~C,D) in Eq. (6.2). Then, the LHFD of the system can be decomposed as
k
(6.6a)
i=l
where
(6.6b)
where
160
I1 ~, ... Vil-I 0 0 ... ... 0 0 ... "~ ]
,.~b°i(1) . 0 0 . . 0 . 1. ~ . ... . ~ V°2-I 0 0 ...
" " "'" " " "'" " "'" " " "'" p-I
0 0 ... 0 0 0 ... 0 ... 1 k ... X ~-
(6.6c)
Proof:
As pointed out in Section 4.4, the matrix sign algorithm can be employed to
compute the transformations for block dlagonalizatlon of the system map A and to
find out the complete set of canonical left/right divisors for the right/left
algorithms, together with the minimal nice selection algorithms (Section 2.3)
in t h e following subsections.
as the modes with Re(kl)<y1, where Y1 is a real number, white the dominant modes
are the modes having Re(ki)>71. If the elgenvalues distribution of the original
system is unknown, the real number YI can be chosen as 71 ffi trace(A)/n, which is
non-dominant modes are the modes wlth [~£I<y2, where 72 is a positive real
number, while the dominant modes are the modes having [li[>Y2. If the
elgenvalues of Ao O
Two methods, namely the frequency-domain reduction method and the time-
domain reduction method, are to be derived for the model reduction of MIMO
The use of Theorem 6.1 yields a frequency-domain model reduction method for
two subsystems as
(6.7a)
G(A) = NrI(A)D:II(k)+Nr2(A)D:21(A)+Dc p
such that Drl(~) contains latent roots of dominant modes and Dr2(~) contains
(6.7b)
CA(~) - SrlC~)V~ll(~)+Vcp -" O(~)
where
A
and
(6.7d)
Dcp ffi Nr2(1)D:~(1)÷Dcp ( f o r d i s c r e t e - t l m e systems)
^
D
cp i s an pXm c o n s t a n t m a t r i x d e r i v e d for matching the z e r o - t h time-moment of
162
o£ t h e o r i g i n a l stable system. •
two s u b s y s t e m s as
(6.8a)
CCA)" D~(k)N~,1(k)+D~ICk)N&2Ck)+D0p
such that D~I(~ ) and V~2(k) contain latent roots o f dominant modes and non-
d e t e r m i n e d by
(6.8b)
GA(A) " DEI( 0p
where
~(o).~2(o)+Dop (for contlnuous-time systems)
A
D0p D (6.8c)
-
and
^ -I
V0p = D~2(1)N£2(1)+Dop ( f o r d i s c r e t e - t i m e systems)
(6.8d)
R
induced pair (AL1,BL1) and the embedded pair (ARI,CR1) of (A,B) and (A,C),
nXn 1
SI =A Ind[ Sign + (A-711n)]¢C
nXn2
S 2 ~ Ind[Sign-(A-Ylln)]¢C , nl+n 2 = n
AR1 =A +
SIASI£Cn l x n l , S +I = ( S ; $ I ) - I s ;
A ;IB£Cnl xm -1
BRI = , SI [Inl Onlxn2][S1 S2]
A CsI£CPXnl
Cg1 =
D0
p ^ =AC(_A)-IB+D_CRI(_AR1)-IBRleCP×m
Iz(t) = ARlZ(t)+BRlU(t)
Y s ( t ) = CRlZ(t)+D0pu(t) -" y ( t )
* -1 * ^
z(0) = [CR1CRI] CRI[CX(0)-DopU(0)] i f p~n I
or
• . -I[cx(0)_;
z(0) = CRI[CR1CR1] 0pu(0)] i f pan 1
Proof:
Theorem 6.5 can be proved from Theorems 6.4 and 4.16 and C o r o l l a r y 4 . 8 . •
T nlxn
V I ~ {Ind[(Sign+(A-'(lln))T]} ¢C '
T u2xn n
V 2 ~ {Ind[(Sign-(A-Tlln))T]} ¢C , nl+n 2 "
+ nlxnl + * * -1
ALl ~ VlAVI¢C , V1 - Vl(VlV I)
BLI ~ VIB¢Cnlxm
- pxn I
CL 1 ~ CVI¢ C ,
VI ~ [VlV
T T2] -T [ I n l 0nlxn 2] T
^Dcp ~ C(-A)-IB+D-CLI(-ALI)-IBLIECp×m
kz(t) = ALlZ(t)+BLIU(t)
Y s ( t ) = CLlZ(t)+Dcpu(t) = y ( t )
165
* -I * ^
Z(0) ffi [CLICLI] CLI[CX(0)-DcpU(0)] if p>n I
or
. .
z(O) ffi CLI[CLICLI] -l[cx(o)_~c p U ( 0 ) ] i f p~n I
Again, the matrix sign algorithm in Theorem 4.21 can be used to compute the
Example 6.1
Solution:
Kronecker indices of this system are gl-3, g2-2 and g3=O. Since g3=0, there
G(A) - D~I(),)NE(X)
where
A3+5.164A2+5.319A÷3.475 -0.208A-I.658 ]
D~(k) -
-0.649k2-6.853k-4.366 )~2+8.836)~+6.688
and
N~(k) - [ 7A2+22"I49A-2"074 16~2+80"527k+55"24 39A2+I83"403A+I08"406]
-4.545A-54.377 6.611A+35.366 8.675A+16.354 J
e x p r e s s e d as
I
Sign÷(A+2.815) = ~ [Sign(A+2"8Is)+I5]
Sign-(A+2.815) = 15-Sign÷(A+2.815)
0 1 -2 4 0 0 2 -4
Sign+(A+2.815) = 0 0 0 Sign-(A+2.815) - I -2
O 0 0 -2 0 0 2
0 O 0 0 0 0 0 0 0 I
167
submatrices by using the block modal matrix M S • To find the block modal matrix
injection maps S I and S 2 in Eq. (4.40) uslng Algorithm 4.1 in Section 4.4. The
o b t a i n e d m a t r i x MS becomes
-2 -1 -4 "]
M S " [Sl,S 2] -
0
0
0
0
1
0
0
0
4
2
1
O
2
1
0
0
8
4
2
1
]
where S I - Ind[Sign (A+2.815) ] and S2 = Ind[Sign-(A+2.8Is) ].
~z(t) = ~z(t)+BDU(t)
y(t) = CvZ(t)+DvU(t)
where
~I"SlASl"
+ [_o., _1.o
-6.0 -2~.S
2.5 ; A~- S~2-
: [,.o,o.o]
-0.5 1.0 4.5 0.0 -8.0
T TT
B D " M;IB = [BRI~B~]
168
BRI = -I.0
0.0
-4.0
1.0
-9.0
2.0
I ; BR2 = [_,.o _,.o]
1.0 0.0 1.0
-I k)NELCk)+D~(k)NE2(A)
where
[ ~2+1.154~+0.903 -0.469 ] [ 5~-0.232 17k+13.611 39A+27]
Dg.l(k) = - 0 . 8 1 0 ~ - 0 . 5 6 5 k+0.846 ; NEI(A) = -7.052 4.225 1.398 j
where
D0p = D (0)N9.2(0) ffi -0.5 0.25 0.0
-0.625 0.25 -0.125
XZs(t) - ARIZs(t)+BRIU(t)
^
* -I * ^
Zs(0) - (CRICRI) CRl[CDZ(0)-D0pU(0)]
eigenvalues of the original system are ~1,2 = -0.5±j0.St ~3 = -I.0, ~4 =-4.0 and
realized, then Eqs. (6.4a) and (6.6a) are the parallel decompositions of G(%).
connecting all subnetworks together with the direct matrix gain Dcp. The
for m u l t i p o r t network s y n t h e s i s .
170
!
!
S
|
'o:o* am: soo'. rso: Ib.oo 1#',m Jb.oo ~.m ab.m
TIME (IN SEC)
F[O. 8 . 1 RE:SPOHSE OF ¥ l ( t ) TO INPUTS , U ~ - I . U 2 - O . U 3 , , O
OR|G|NN~ SYSTEfl
...... REDUCEDMODEL
I
'o:oo am: |w" Tso: Ib,m ~z'.u Ib.m i)'.m Jb.m
T~nE (~N $E~
FIG. 6.2 R[SPOHSi[ OF Y 2 ( t ) TO INPUTS " U I ' I . U 2 " O . U 3 = O
ORIOTN+U. SrStl~
...... REDUCEDPIOOE~
171
.,
!
'~° °00
T|ftE (IN SEC)
F|0. S . S RESPONSE OF ¥ 1 ( t ) TO |NPUTS * U I - 0 . U 2 - 0 . U 3 w l
01~ GZHAL STSTI~I
REDUCED hODlel
N aim
: sod
: f'so to'.m la'. N tk.o0 t~.N ab.N
Tilt[ KXR SEC)
FXO. 6 . 6 R(SPOHS[ OF ¥ 2 ( t ) TO INPUTS i U l m O . U 2 - 0 . U 3 - I
ORtGXNAL S~'STEn
...... REDUCED r10O(l,.
173
decomposed into the forms in Eq. (4.28), then from Theorem 4.12, we obtain the
RMFD as
(6.9)
G(1) = N r ( l ) [ D r 2 ( l ) R ( ~ ) l - l + v = Nr(X)R-I(1)D:~(1)+D
where Dr2(1) is the right characteristic l-matrix of (A2,B 2) and R(1) is a right
divisor of D ( ~ ) .
r
Lemma 6.1 A reachable MIMO system of dimension greater than 1 can always be
(6.10a)
IX(t) = AX(t) + Su(t)
where
A =
A1 AI2 1; AI¢C
n.xn.
I i
; i=1,2
(6.10c)
0 A2
n.xm
T T T • (6.10d)
B = [BI,B2] ; BI¢C ; i=1,2
174
pxn. (6.10e)
C = [Ci,C2] ; Ci¢C x; i=l,2
V = 0 pXm (6.100
(6.11a)
G(A) = Nr(A)R-I(A)D:~(A)
(A2pB2~Cr2~Dr2) as its associated 'minimal realizationp where Cr2 and Dr2 are
defined in Lemma 2.1. Also
Nr(A ) = C T : I ~ r ( A ) (6.11b)
(6.11c)
R(~) " Cr2VT:ISr ()0 ÷Dr2D r (),)
where Dr(A) is the right characteristic A-matrix o f (A,B) and V is the canonical
Proof:
system.
Lemma 6.2 An observable MIMO system in Eq. (6.2) can always be decomposed into
(6.12c)
L(X) = ~E( X)T;IsBEI+D~( X)DEI
injection map from S~X, and (AI~Cl) is the embedded map of (AIC) in S p whereX
Proof:
latent vectors of the X-matrlces, we study the MIMO system (A,B,C1 D) in Eq.
(6.2) which can be transformed into the following Jordan canonical form, using
where
n.xn.
Aj = M-IAM = block dlag[Aji,i=l,...,k];Aji¢ C I i (6.13d)
m n.xm
T T _* . ~ 1 (6.13e)
Bj = M-IB = [BjI~Bj2~... ~BJk j ; B j i ¢ C
pxn.
(6.13f)
Cj - CH - [Cjl,Cj2,...,Cjk];CjieC
• li
nxn.
The generalized left eigenvectors, pjicC t associated with Aji can be
determined by partitioning the modal matrix M such that the following matrix
(6.13h)
Apj i = PjiAji, i = l , 2 , . . . , k
for the system in Eq. (6.2), or the equivalent system in Eqs. (6.10) and (6.13),
such that some undesirable elgenvalues and associated eigenvectors of the system
177
can be replaced by deslrable ones. Since the system map of Eq. (6.13) can be
(6.14a)
~XD(t) = ADXD(t)+BDU(t)
(6.14b)
y(t) = CDXD(t)+DDU(t)
Xv(t) " T~Ix(t) (6.14c)
where
- ; BD -
rBvl] [itI ; XD(t) =
o '~2 L BD2 XD2(t)
Cv = [Cvl,Cv2]; DD=0
ADI is the collection of the unimportant system modes not to be controlled~ and
AD2 contains all dynamic modes to be controlled. XDi(t) , BDi and CDi are
(6.15a)
u(t) - -FDXD(t)+GFr(t) - -FDT;Ix(t)+GFr(t)
(6.15b)
~XD(t) = ~XD(t)÷BDU(t)
y(t) " CDXv(t) (6.15c)
where
178
The modal control design problem is to find FD2 such that AD2-BD2FD2 has the
observe that the modal control of the entire system is equivalent to the
5.1, we obtain the control law for the sub-system (AD2,BD2) as follows.
Lemma 6.3 Let Dr2(1) be the right characteristic l-matrix of (AD2,BD2) which
Then, the feedback gain FD2 and input gain G F in Eq. (6.15) become FD2 ffi Fc2Tc2
Fc2~r2(1) = GFDrc2(I)-Dr2(I)
controller form, and Drh 2 and Drhc2 are the leading-column matrices of Dr2(1)
Lemons 6.1, we observe that Drc2(1) is the canonical left divisor of the closed-
loop characteristic l-matrix, Drc(1). Therefore, the latent roots and left
generalized latent vectors of Drc2(1) are also part of the latent roots and left
control law via latent roots and associated latent vectors assignment as follows.
179
Assume that the desired closed-loop latent roots and left latent vectors
n2Xn2
Aj2 = block diag[J2i~i=l~...,k]EC ;n2 = ~. (6.16a)
i= I i
- -T ~T ITECn2 xm (6.16b)
Bj2 = [Bj21~... ~ J2k J
Assume that (Aj2,Bj2) is a reachable pair and that the Kronecker indices of
(Aj2,Bj2) are equal to those of (AD2,BD2) in Eq. (6.14). Then, the feedback
gain F D and the input gain G F of the modal control law in Eq. (6.15a) can be
determined by FD = [0,FD2], where FD2 = Fc2Tc2 and G F = Drh2Dr c2" Fc2 can be
(6.16c)
Fc2~r2(l) = GFDrc2(1)-Dr2(l)
where Dr2(1) , ~r2(l) are the same as those defined in Lemma 6.3, and Drc2(l) is
Proof:
Theorem 6.7 can be proved from Theorem 5.4 and Lemma 6.3. •
From Theorems 5.2 and 6.3~ the closed-loop RMFD of the system in Eq. (6.10)
with D = 0px m and the modal control law in Eq. (6.15a) can be represented as
where
180
= A_BFDT;I (6.17b)
and Nr(A)D:I(A) = G(A) is the canonical RMFD of the open loop system in Eq.
^
Theorem 6.8 Drc(k) defined in Eq. (6.17c) is column reduced. If the open loop
(6.18)
Dr(k) = R(k)Dr2(A)
where Dr2(A) is the right characteristic k-matrix of (AD2,BD2) , then Drc(A) can
be factored as
^ (6.19)
D r c (l) = R(A)Drc2(A)
Proof:
Drh , the leading column matrix of Dr(1) , is nonsingular, and therefore Drc(A) is
column reduced.
If Dr(k) is factored as in Eq. (6.18), then from Theorem 4.1, Lemma 6.1 and
Dr2(~ = Cr2(AIn2-AD2)-IBD2+Dr2
181
R(x) = ~ 2V~DIT~I~r(X)+~r2Dr(X)
Drc2(k) = Cr2(AIn2-AD2+BD2FD2)-IBD2GF+Dr2GF
where
Dr2FD2 = (Im-~Tr2(0)~r2(0))Drhl2FD2 =- 0
Theorem 6.9 The closed-loop RMFD of the modal controlled system can be
where Nr(~) is the numerator of the open-loop canonical RMFD, R(%) is the right
Theorem 6.9 reveals the important feature of the modal control design that
the left divisor of the denominator of the open-loop RMFD has been changed to the
Note that Gc(~) in Eq. (6.20) may not be a canonical RMFD. Using Theorem
(6.21)
Gc(X) = Nrc(X)Drcl(X)
matrix. The result in Eq. (6.21) implies that although the zeros of the open-
loop system remain unchanged under modal control laws~ the numerator of the
Example 6.2
Given a 3-input 2-output continuous time system in the form (6.2) where
G(X) = NrfX)Drl(X)
where
Nr(X) = [ -4X2+13~-13
-5 X2+I0 A+3
7%-12
12A+7
0.0 ]
0.0
and
X3-2X2-X+2 -~.2+X+2 -1 ]
D (X) = 0 ~.2+2 k+l -2
r
0 0 1
design for this system. For these purposesj the system map A has to be block-
S i g n ( r 0 ) ( A ) = 15
184
From Theorem 4.20, we have the canonical injection map S(r0,rl ) and the
Defining
X ( t ) = T D XD(t)
~XD(t) = ADXD(t)+BDU(t)
y(~) = CDXD(t)+DDU(t )
185
where
BD=TI__[B02
01]
C D -- CT D = [CDI , Cv2]
and
1.0 1.0 ]
AD2 ffi 0.0 2.0
CDI ffi
2.8
10.2
3.4
4.6
4.0
II.0
]; CD2 =
1.1778
-1.6667
1.6000
1.0000
]
D D -- 02x3
Dr(X) = Dr2(X)R(X)
186
A-2 0 -i 1
Dr2(A) = -2 A-I -2
0 0 i
and R(A), which is a right divisor of Dr(A) ~ can be computed by using Eq.
(6.11c):
ALl A1 0 ]
R(A) = 2A+2 A+I 0
0 0 i
Thus, from Lemma 6.1, the seml-cascade RMFD of the entire system is given by:
Note that Dr2(X) contains two unstable modes of the system with latent roots 1
and 2, and R(A) contains three stable modes with latent roots -i~-I~ and -I.
Since AD2 contains two unstable modes with eigenvalues 1 and 2~ we shall
design a modal controller to stabilize these two modes. Using the left latent
structure assignment in Theorem 5.4~ we select the latent roots and latent
~21 = - 2 P210 = [ I -I O] T
and
T
A22 = -4 P220 " [i I O]
s o we h a v e
187
and
- = 1 -I 0 ]
B J2 [ i 1 0
;%+3 i 0 ]
I ),+3 0
Drc2 ()') = 0 0 I
Thus, from Eq. (6.15a) and Theorem 6.7, we have the modal control law
u(t) = -FX(t)+GFr(t )
where
and
Gc(1) -- Nr(;k)Drcl(X)
where
^
rc Drc2 (
188
13+312+k-I -A2-3k-2 0 ]
= 3k3+8A+5 k2+3~+2 0
0 0 1
^-I
Note that Gc(A) - Nr(A)Dr~(A) is not a canonical RMFD. The canonical RMFD of
the closed loop system can be found from the closed loop state equations:
kX(t) = AX(t)+BGFrCt)
yCt) = CXCt)+Du(t)
where
= A-BF
where
X3+6X2+10X+5 -12-3X-2 0 ]
Drc(k) = -X-I X2+3X+2 0
0 0 i
Drc()%) = Drc(1)Ur()%)
and
Nrc()%) = Nr(l)Ur()% )
i 0 0 ]
Ur()%) = -3 1 0
0 0 1
Synthesis
and LMFDs are derived. The corresponding state-space cascade realizations are
(6.22)
G()k) = GI(~,)G2()%)
where G()%), GI()%) and G2()%) are described by MFDs, then Gl(S) and G2()%) are
called the left factor and right factor of G()%), respectively. The
Deg(D(X)) = Deg(DI(A))+Deg(D2(X))
(6.23)
where D(A), DI(A) and D2(A) are the left or right characteristic X-matrices of
the systems described by G(A)~ GI(X) and G2(k) ~ respectively. The minimal
factorization of G(X) in Eq. (6.22) can be generalized as
(6.24)
G(A) = GI(X)G2(A)...Gk(A)
k
with Deg(D(~.)) = ~ Deg(Di(A)).
i=l
nonsingular. Q
Assume that DI(X) and NI(X) are a column-reduced left canonical divisor and a
Let
(6.25c)
N1 ( X)D;I(X) = HI( X)+;I( ~)~;1(X)
(6.26a)
G(I) : 02(1)Gl(1)
where
(6.26b)
Cl(D : N l ( D D [ l (
(6.26c)
G2(~ = N 2 ( ~ D 2 1 ( ~
and
(6.26d)
Nl(X) : N'I( X)+D2( X)HI(X)
(6.26e)
Proof:
Since
Therefore
into the product of two nonsingular RMFDs. Since D2(%) may not be in a column-
192
reduced canonical form, G2(X) may not be a canonical P~FD. Howeverj using
Proof:
Since the left factor GI(X) and right factor G2(X) can be represented as
(6.27)
GCX) = Gk(~)Gk_I(~)...GI(X)
The main feature of the RMFD factorization in Theorem 6.10 is that the
poles and zeros with associated latent vectors of GICX) and G2(X) can be
para-Hermitlan MFD.
Let D~(X) = D2(~)DI(X ) and NE(X) = NI(A)N2(X) , where D1CX) and NI(X) a r e a row-
D£(X) and Ng(X), respectively. Let D;I(x) NI(X) - HI(X)+D21(X) NI(X) , where
HI()~) is a l-matrix and D;I(x) NI(X) is proper. Also, let (AI,BI,CI,D I) be a
minimal realization of D21(~k)~l(%). Then
(6.28b)
O(l) = Cl(k)C2(L)
where
G(I) -- Ol(1)G2(},)...Ok(1)
(6.2g)
(6.30a)
%Xl(t ) = AiXl(t)+BlUl(t )
(6.30h)
Yl(t) ~ ClXl(t)+DlUl(t )
and
(6.30c)
)d(2(t) = A2X2(t)+B2u2(t)
(6.30d)
Y2(t) = C2X2(t)+D2u2(t )
where Ai,Bi,C i and Di, i=I~2 are matrices of appropriate dimensions. For
cascade connections, we let u(t) = Ul(t) , Yl(t) = u2(t) , and y(t) = y2(t), where
u(t) and y(t) are the input and output of the overall system, respectively.
194
rxict]
[AI01
xict] [i]
= + u(t) (6.31a)
ILX2(t) B2C 1 A 2 JLX2 (t) B2D 1
xl(t) ]
(6.31b)
y(t) = [D2CI,C2] + D2DlU(t)
x2(t)
The above result can be generalized to the case of many subsystems as follows.
Theorem 6.12 Given a nonsingular MFD, G(~)~ which has a minimal factorization
(6.32)
G(I) = Gk(%)Gk_I(I)...GI(I)
expressed as
AI 0 ... 0 0
B2C 1 A2 ... 0 0
B = [B~,(B2DI)T,(B3D2DI)T,...,(Bk_IDk_2...DI)T,(BkDk_I...DI)T] T
C = [Dk...D2CI~Dk...D3C2~Dk...D4C3~...,DkCk_I,Ck]
D = [DkDk_I...DI]
195
Proof:
subsystems. •
G(~,) = GI(1)G2(%)...Gk(1)
the subnetworks G.(%) can be constructed in a way such that each G.(I) contains
1 1
desirable latent roots and associated latent vectors of G(1). Then, we can
realize each subnetwork with each G.(1) and cascade them to construct the
1
complete network for G(1). Since the factorization of G(~) is minimal~ the
Example 6.3
network be
G(1) = Nr(1)Drl(1)
where
Nr(~) = [ 0 t tO
The zeros of G(1) are -5, -I0 and the poles of G(1) are (-l,-l,-10±jl).
196
Find (i) The minimal realizations of D -I in the controller canonical form and
r
Solution
D;I(L) = Cc(%I4-Ac)-IBc+Sc
where
_ 0.00 1.00 i 0.00 0.00
A = -8.75 -n.8o lO.OO 1.8o
c o.oo - ~.~ ~To~ t.oo
-2.05 - 2.05 I-9.20 -10.20
B = i 0
c ~-~
0 1
o o o]
Cc 0 0 I 1 0
Dc 02x 2
(Ac,Bc~Cc,D c) as
D;I(1) = Cj(II4-Aj)-IBj+D J
where
Aj
.
lit0 0] E11] r01002s
.
-lJ o
0. 1-10
.
o
.1 ; Bj
o 1
i--0"
-z i-o.I o.15 1
; cj =1-o-.i--_-o-;1- -| ; DJ
=
02×2
0 I-I -i0 I i L 0.0 -0.125 J
197
Note that the latent roots of Dr(1) have been arranged in two groups; {-i)-I}
and {-10±jl} as shown in the system map Aj. Using Theorems 4.9 and 4.12
together with the Jordan form minimal realization quadruple of D-l(1)i yields
r
Dr(I ) = DL(1)DoDr(1)
where
Nr(1) = [ 0 10
=o,[' ~ ][~ ~+ (~,o °o)] [,~+ ( o::~ -o.~,.~)]-' ['~+ ( ~" -'4-'
G(X) = G2(X)Gt(X)
where
-I
o,<~<>-[~ ~÷ ,~.ooo,r
o.o~ o.~,~1[~+~.~,,, (~'~ I)}
198
^ ^
where
^ [0 01~ ~] ^
Cr 0 0 ! 0 i ; Dr = 02x2
The cascade realization of G(X) using four one-port R-C networks and summers but
6.7. The poles and zeros of the first subnetwork~ described by GI(%) , are
{-10±jl} and {-5,-I0}, respectively whilst the poles and zeros of the second
realized cascade multi-port network contains the prespecified poles and zeros in
each subnetwork.
4
In '
.= i (D
Ul i/ "~~ !
I
I
[i] R.E. Kalman, P.L. Falb, and M.A. Arbib, Topics in Mathematical System
Theory, McGraw-Hill, New York~ 1969.
[2] S. Barnett, Matrices in Control Theory, Van Nostrand Reinhold Co., New
York, 1971; Second Edition, Krieger, Florida, 1984.
[3] S. Barnett, Polynomials and Linear Control Systems, Marcel Dekker, New
York, 1983.
[4] M.K. Sain, Introduction to Alsebraic System Theory, Academic Press~ New
York, 1981.
[7] B. Porter and R. Crossley, Modal Control, Taylor & Francis Ltd., London,
1972.
[9] D.D. Siljak m Larse-Scale Dynamics Systems, North-Holland, New York, 1978.
[12] L. Pernebo, Alsebraic Control Theory for Linear Multivariable Systems 3 Lund
Institute of Technology, Sweden, 1978.
[14] R.V. Patel and N. Munro, Multivariable System Theory and Design, Pergamon
Press, New York, 1982.
[17] P. Lancaster, Lambda Matrices and Vibrating Systems, Pergamon Press, New
York, 1966.
[18] J.E. Dennis, J.F. Traub, and R.P. Weber, "The Algebraic Theory of Matrix
Polynomials", SlAM J. Numerical Analysis, 13, pp. 831-845, 1976.
[23] A.S. Markus and I.V. Mereuca, "On the Completen n-Tuple of Roots of the
Operator Equation Corresponding to a Polynomial Operator Bundle"9 Math.
USSR Izvestija 9 Vol. 7, pp. 1105-11289 1978.
[24] L.5. Shleh~ Y.T. Tsay~ a n d N.P. Coleman 9 "Algorithms for Solvents and
Spectral Factors of Matrix Polynomials" 9 Int. J. Systems Science, Vol. 12,
pp. 1303-1316, 1981.
[25] L.S. Shieh and Y.T. Tsay 9 "Transformations of Solvents and Spectral Factors
of Matrix Polynomials and their Applications", Int. J. Control 9 Vol. 34,
pp. 813-8239 1981.
[26] Y.T. Tsayy L.S. Shieh 9 R.E. Yates and S. ~arnett, "Block Partial Fraction
Expansion of a Rational Matrix"9 Linear and Multilinear Algebra, Vol. Ii,
pp. 225-241, 1982.
[27] L.S. Shieh and Y.T. Tsay, "Block Modal Matrices and their Applications to
Multivariable Control Systems"9 Proc. lEE, Vol. 129, Pt. D. 9 pp. 41-48j
1982.
[29] Y.T. Tsay and L.S. Shieh 9 "Some Applications of Rational Matrices to
Problems in Systems Theory", Int. J. Systems Science~ Volo 13, pp. 1319-
1337, 1982.
[30] Y.T. Tsay and L.S. Shleh, "Block Decomposition and Block Modal Control of
Multivariable Control Systems", Automatiea, Vol. 199 pp. 29-$0~ 1983.
[31] Y.T. Tsay and L.S. Shieh, "Irreducible Divisors of A-Matrices and their
Applications to Multivariable Control Systems"~ Int. J. Control, Vol. 37,
pp. 17-36, 1983.
[32] L.S. Shleh, Y.T. Tsay and R.E. Yates, "Block Transformations of a Class of
Multivariable Control Systems to Four Basic Block Companion Forms", Int. J.
of Computers and Mathematics with Applications, Vol. 9, pp. 703-714, 1983.
[33] D.G. Luenberger, "Canonical Forms for Linear Multlvariable Systems," IEEE
Trans. Automat. Control, Vol. AC-12, pp. 290-293, 1967.
[34] W.A. Wolovich and P.L. Falb, "On the Structure of Multivariable Systems" 9
202
[37] A.C. Antoulas~ "On Canonical Forms for Linear Constant Systems"~ Int. J.
Control, Vol. 33, pp. 95-122, 1981.
[38] V. Kucera~ Discrete Linear Control, John Wiley & Sons, New York, 1979.
[39] A.W. Naylor and G.R. Sell, Linear Operator Theory in Engineering and
Science, Springer-Verlag~ N.Y.~ 1982.
[40] L.S. Shieh and Y.T. Tsay, "A Minimal Nice Selection Algorithm for the
Canonical Controller Structures of Multivariable Control System", U.S. Army
Research Report, 1983.
[41] P. Lancaster and P.N. Weber, "Jordan Chains for Lambda Matrices", Linear
Algebra and Its Applications, Vol. I, pp. 563-569, 1968.
[44] L.S. Shieh and N. Chahin, "A Computer-Aided Method for the Factorization of
Matrix Polynomials", Int. J. Systems Science, Vol. 12, pp. 305-323, 1981.
[45] J.E. Dennis, J.F. Traub and R.P. Weber, "Algorithms for Solvents of Matrix
Polynomials", SLAM, Numer. Anal., Vol. 15~ pp. 523-533, 1978.
[46] E.D. Denman, "Matrix Polynomials, Roots and Spectral Factors", Applied
Mathematics and Computation, Vol. 39 PP. 359-368, 1977.
[48] R.E. Kalman, "On Partial Realizations, Transfer Functions, and Canonical
Forms", Act Polytechnica Scandinavia, pp. 9-32, 1979.
[52] P.A. Fuhrmann and J.C. Willems 9 "A Study of (A,B)-Invariant Subs,aces Via
Polynomial Models", Int. J. Control, Vol. 31, pp. 467-494, 1980.
[54] M.S. Mahmoud and M.G. Singh, Large Scale Systems Modellin~, Pergamon Press,
N.Y., 1981.
[55] H. Kwakernaak and R. Sivan, Linear Optimal Control Systems, John Wiley &
Sons, Inc. N.Y., 1972.
[56] J.D. Roberts, "Linear Model Reduction and Solution of the Algebraic Riccati
Equation by Use of the Sign Function", CUED/B-CONTROL/TR 13 Report,
Cambridge University, 1971, Also, Int. J. Control, 1980, Vol. 32, pp. 677-
687.
[57] E.D. Denman and A.N. Beavers, "The Matrix Sign Function and Computations in
Systems", Appl. Math. & Comput., Vol. 2, pp. 63-94~ 1976.
[58] R.L. Mattheys, "Stability Analysis Via the Extended Matrix Sign Function",
Proc. IEE, Vol. 125, pp. 241-243, 1978.
[59] L.A. Balzer, "Accelerated Convergence of the Matrix Sign Function Method of
Solving Lyapunov, Riccati and Other Matrix Equations", Int. J. Control,
1980~ Vol. 32, pp. 1057-1078.
[60] L.S. Shieh, Y.T. Tsay and R.E. Yates, "Some Properties of Matrix Sign
Functions Derived from Continued Fractions", Proc. lEE, Vol. 130, Part D,
pp. 111-118, 1983.
[61] C.T. Chen, Introduction to Linear System Theory, Holt, Rinehart and
Winston s New York, 1970; Second Edition, 1984.
[64] B.C. Moore, "Eigenvalue and Eigenvector Assignment Using State Feedback",
IEEE Trans. Auto. Control, Vol. AC-21, pp. 689-692, 1976.
[65] E.J. Davlson, "Multlvarlable Systems Subject to Arbitrary Output and State
Feedback", IEEE Trans. Auto. Control, Vol. AC-18, pp. 24-32, 1973.
[66] J.J. D'Azzo and C.H. Houpls, Linear Control System Analysis and Design,
McGraw-Mill~ New York, 1981.
[68] W.A. Wolovich and P.L. Falb, "Invariants and Canonical Forms Under Dynamic
Compensation", SlAM J. Contr., Vol. 14, pp. 996-1008, 1976.
[69] P.L. Falb and W.A. Wolovich, "The Roles of the Interactor in Decoupling",
Proc. JACC~ pp. 1456-1460, 1977.
[70] P.L. Falb and W.A. Wolovich, "Decoupling in the Design and Synthesis of
Multlvariable Control Systems", IEEE Trans. Auto. Control, Vol. AC-129 pp.
651-659, 1967.
204
[73] D.Q. Mayne~ "The Design of Linear Multivariable Systems"~ Automatic, Vol.
9~ pp. 201-207, 1973.
[75] D.C. Youla, H.A. Jabr and J.J. Bongiorno, "Modern Wiener-Hopf Design of
Optimal Controllers - Part II: The Multivariable Case", IEEE Trans. Autom.
Control, Vol. AC-21, pp. 319-338, 1976.
[79] B.W. Dickinson, T. Kailath and M. Morf, "Canonical Matrix Fraction and
State-space Descriptions for Deterministic and Stochastic Linear Systems",
IEEE Trans. Autom. Control, Vol. AC-19, pp. 656-667, 1974.
[83] S.V. RaG and S.S. Lamba, "Suboptimal Control of Linear Systems Via
Simplified Models of Chidambara", Proe. lEE, Vol. 121, pp. 879-881, 1974.
[84] L.S. Shieh, S. Yeh and H.Y. Zhang, "Realizations of Matrix Transfer
Functions Using RC Ladders and Summers", Proc. lEE, Vol. 128, Part G, pp.
i01-II0, 1981.
[86] F.F. Kuo~ Network Analysis and Synthesis, John Wiley & Sons~ New York~
1966.
[87] G.C. Temes and J.W. LaPatraj IntEpduction to Circuit Synthesis and Design,
205
[88] R.W. Newcomb, Linear Multiport Synthesis, McGraw-Hill, New York, 1966.
[90] Y.T. Tsay and L.S. Shleh, "Block Modal Control Vla Cascade Multivariable
Control Structure", Int. J. Systems Science, Vol. 15, pp. 199-214, 1984.
[91] D.C. Youla, "On the Factorization of Rational Matrices", IRE Trans.
Information Theory, Vol. IT-7, pp. 172-189, 1960.
[93] B.D.O. Anderson, K.L. Hitz and N.D. Diem, "Recursive Algorithm for Spectral
Factorization", IEEE Trans. Circuits and Systems, Vol. CAS-21, pp. 742-750,
1974.
[94] H. Bart, I. Gohberg and M.A. Kaashoek, Minimal Factorlzations of Matrix and
Operator Functions, Birhauser Verlag, Boston, 1979.
[96] P. Van Dooren and P. Dewilde, "Minal Cascade Factorization of Real and
Complex Rational Transfer Matrices", IEEE Trans. Circuits and Systems, Vol.
CAS-28, pp. 390-400, 1981.
[98] L.S. Shleh, Y.T. Tsay and R.E. Yates, "Cascade Realization of Multiport
Networks Via Factorlzation of Rational Matrix Functions", 7th Int.
Symposium on the Mathematical Theory of Networks and Systems, Stockholm,
Sweden, 1985.
[I00] J.W° Helton~ "Systems with Infinite-Dimensional State Space: The Hilbert
Space Approach", Proc. IEEEj Vol. 64, pp. 145-160, 1976.
[102] R.P. Roesser, "A Discrete State-Space Model for Linear Image Processing",
IEEE Trans. Autom. Control~ Vol. AC-20~ pp. i-i0~ 1975.
[103] B.G. Mertzios and P.N. Paraskevopoulos, "Transfer Function Matrix of 2-D
Systems", IEEE Trans. Autom. Control, Vol. AC-26, pp. 722-725, 1981.
104] N.K. Bose, Multidimensional Systems: Theory and Applications~ IEEE Press,
New York, 1979.
[105] S.K. Mitra and M.P. Ekstrom, Two-Dimensional Digital Signal Processing,
Dowden, Hutchinson and Ross, Stroudsburg, 1978.
206
[106] K.J. Astrom~ Introduction to Stochastic Control Theory~ Academic Pressj New
York, 1970.
[107] B.D.O. Anderson and J.B. Moore~ Optimal Filterin~ Prentice-Hall~ Englewood
Cliffs~ New Jersey~ 1979.
[109] E.D. Denman~ J. Leyva-Ramos and G.J. Jeon~ "The Algebraic Theory of Latent
Projectors in Lambda Matrices"~ Applied Mathematics and Computation~ Vol.
9~ pp. 273-300, 1981.
[II0] L.S. Shieh~ Y.T. Tsay and R.E. Yates~ "State-Feedback Decomposition of
Multivariable Systems Via Block-Pole Placement", IEEE Trans. Auto.
Control~ Vol. AC-28, pp. 850-852~ 1983.