A S - B. Kim: ( (XI Fs Al
A S - B. Kim: ( (XI Fs Al
AUTOMATIC
CONTROL. AC-32, NO. 5 , MAY 1987 423
~~
~ _ _ _ _ ~ . _ _ _ _ _ _ ~ _
__ _ _ _ ~
~~
~~ -___ closed-loop
~ polesto stay inside thespecified disk D are determinedfor the
proposed control.
N(rJ ) DF J 1 F
__ ________
I. INTRODUCTION
1 0 1 2 -6 0 2213 4 1013 ~
or
=F x ~ (3) '
such that all poles of the closed-loop system, i.e., the roots of
are assigned in the disk D with the center at a + j 0 and the radius r which
may be shown in Fig. l(a) for the continuous system and in Fig. l(b) for
the discrete system.
The assignment of all poles of the closed-loop system in the specified
disk D shown in Fig. 1 is defined as D-pole assignment problem in this
note. To solve this problem, we present the conditions such that all (a) (b)
eigenvalues of a matrix should be located in a specified disk as shown in Fig. 1. (a) Disk D in the left half of the complex plane. (b)Disk D in the unlt disk with
Fig. 1. Note that the following conditions are satisfied for both continuous the center at the origin.
and discrete systems.
Lemma I: Consider the matrix equation
become zero.Therefore, it proves that there exists a positive definite
-aA*P-aPA+A*PA+(a2-r2)P= -0 (5) solution P such as Iimk+= PI- = P which satisfies ( 5 ) .
Remark: We may see from the abovelemma that a different expression
where Q is an arbitrarily positive definite matrix, and * denotes the for (5) can be given
conjugate transpose of amatrix, and cy, r are scalars.Then,the
eigenvalues of matrix A are located within the specified disk as shown in (A-aZ)' (A-aZ)
P-- p= --
Q
(9
Fig. 1 if and only if there exists a positive definite solution P satisfying r r r2
(5).
Proof: Let A. v be an eigenvalue and right eigenvector of A , then and that the eigenvalues of the matrix ( A - a Z ) / r are located within a
unit disk with its center at the origin if and only if there exists a positive
A u = X u , u*A*=Xu* (6) definite solution P that satisfies (5). In the previous lemma. we considered
only the case such that Q is positive definite. Also. in the case that Q =
and substituting these expressions into ( 5 ) yields HTH, the lemma is satisfied for the positive semidefinite Q as far as the
pair ( A , H ) is observable (7).
{ - a ( X + h ) + IXIz+a2-r2]u*Pu=-o*Qu. (7) Consider the state feedback
Let h = x + j y and using it in (7). we obtain u=Fx (3)
{ - 2 o ~ u + x ~ + y ~ + a ~ - r ~ } u-v*Qu
*Pv= (8) then, from Lemma 1, we can have the conditions that the eigenvalues of
the closed-loop mamx A + 5 F should be located within a specified disk
or D.
Theorem I : Consider the following matrix equation:
{ ( x - a ) 2 + y 2 - r z ] u * P u =-u*Qv. (8) '
-a(A+BF)*P-aP(A+BF)+(A+BF)*P(A+BF)
Since Q is positive definite, the positivity of P yields
+ ( a z - r 2 ) P =- 0 (13)
(x-a)'+y2-r2<0 (9)
where Q is positive definite.Thenthe eigenvalues of ( A + 5F) are
which means the condition that all eigenvalues of the matrix A should be within a specified disk D as shown in Fig. 1 if and only if there exists a
located in a specified disk with radius r a n d center at a + j0 (Fig. 1). positive definite solution P satisfying (13).
Conversely, let the eigenvalues of matrix A be inside a specified disk of As a method to choose a state feedback law that satisfies our problems,
Fig. 1. that is. the eigenvalues of ( A - d ) / r are within the unit disk we present the following theorem by using the discrete Riccati equation.
with its center at the origin. Consider the following matrix: Theorem 2: The state feedback law
u= -(rZR+BTPB)-'BTP(A-aZ)x (14)
k l
Pk= 7 ( ( A- cyI)*)'Q(A- aZ)'. (10)
assigns all the closed-loop poles of a continuous system (1) or discrete
i=o
system (2) in the disk D shown in Fig. 1. where P is a positive definite
Since Q > 0, it is shown that the right term of (IO) is positive definite for symmetric solution of the Riccati equation
all i. From ( 5 ) and (lo), we can obtain
(A-an'. ( A-a]) ( A - a1)
ps- +HTH--
J
1 1 r r r
-(A-aI)*Pk(A-aI)-Pk= -7 Q
rz r
. P B ( r Z R + B J P B ) - i B (TAP-y~ 2 )
I (15)
-a(A+BF)'PPaP(A+BF)+(A+BF)'P(A+BF)+(a2-r2)P
By hypothesis of the matrix A , the following is yielded as k -+ 03
= -r2(FTRF+H'H) (15)'
(llA -aZll/r)k+l+O where
which means that the right-hand side of (12) and the second term of (1 1) F= - ( r ' R + B J P B ) - ' B J P ( A-aZ).
IEEE TRANSACTIONS OK AUTOMATIC
CONTROL, VOL. AC-32. NO. 5, MAY 1987 325
In ( l 5 ) ' , let where P given by (15) is the solution PCof the following matrix equation:
and using the above Lemma 1, its remarks and Theorem 1, it proves that 1
=- { ( A - B R , ' B T P , ) T P ~ ( A - B R ~ - ' B T P , ) + P , B R c I r 2 R R ~ - ' B T P ,
the poles of the closed-loop system, Le., the eigenvalues of the matrix A 01
+ BF are located within the specified disk D . +r2HrH+(a2-r2)Pc}. (22)
Remark: We may also prove Theorem 2 from the following. For the
discrete system The algorithm to find the positive definite solution of (15). which is
similar to the one in (4), is given as follows. Let X, (i = 1, 2, --,n ) be an
eigenvalue of the matrix (18), and the corresponding eigenvector be v,.
Then the following is yielded:
it is well known that the following control law
+ B ( r 2 R + B T P B ) - ' B 7AP ,- Q r ) u,=O (23)
Defining ui as
I -Lpvi
is minimized and all the closed-loop poles, Le., the eigenvalues of the
matrix (23) and (24) may be rewritten as
are located within the unit disk with the center at the origin (3). It is easily
found that the eigenvalues of the matrix
[ -H7H
'1 [:,]
I
=X, [ Z Br2R-'B7
A-B(r2R+BTPB)-'BTP(A-aT) (18)'
By solving the above generalized eigenvalue eigenvector problem for
eigenvalues to be less than 1, P may be given by
are located in the specified disk D.
P = [ u , ,u2, ---,U,l[Ul, v2, ----,u,l-'. (26)
Thus. the theorem means that the state feedback control law (14)
stabilizes the closed-loop system with poles located in the specified disk.
111. D-STABILITYMARGIN
Remark: For the continuous system, let the control given by (14)
denote The advantages of the D-pole assignment presented in the previous
section are not only that fewer parameters are used in the specification for
u=Fx
the control law in comparing to assignment of all poles, but also that the
then it may be the optimal control minimizing the criterion function uncertainty of the plant model: which may lead to poles outside of the disk
D,can be treated similartothe stability margin. To present the gain
J= sr (llxll&+ l l 4 l ~ Jdt (19)
margin and the phase margin for all poles staying in a specified disk D,
instead of the control law (14), the following control, as shown in Fig. 2,
is considered to be supplied to the plant.
for certain Qc and R,. Therefore, F = -R,BTPc, where PC is the
u = -K(r2R+BTPB)-IBTP(A-aI)x (27)
positive definite symmetric solution of
P KFx (27)'
ATP,+P,A+Qc-PcBR~lBrPc=O (20)
where
which yields
K=diagk 2{ k, l , -----,
k,"}. (28)
( A + B F ) T P , + P c ( A + B F ) + F T R c F -0
= -c.
Then we are led to the following lemma which means that the closed-loop
system has the robustness for parameter perturbation as shown in Fig. 2 in
On the other hand, (15) may be rewritten as view of the above statement.
Lemma 3: Suppose that there is a positive definite solution satisfying
(15) with R = diag { T I , TI, -------.
F,,,}. Then the gain margin for D-pole
assignment is given by
+ F ' r 2 R F + r 2 H T H + ( a 2 - r 2 ) P ] + F T R , F(15)' I 1
-<k,$-
1 +a, - 1 -a,
where a < 0. Those two equations imply that the weighting matrix Qc of where = 2, ____,
pn
(19) is given by
1
Qr= - ;
{(A+BF)TP(A+BF)+F7r2RF+r2HTH
where X , is a maximum eigenvalue of the matrix BTPB.
+ ( a Z - r 2 ) P } - F T R c F(21) Proo$ As the general gain margin of the discrete optimal control,
426 IEEE TRAKSACTIOh'S OK AUTO!vlATIC CONTROL. VOL. AC-32. S O . 5 . MAY 1987
1 X
II n I
Fig. 2. Closed-loop system employingthe pole assignment feedback law in the presence
of unknown gain { k , } .
( A - L Y I + B K F()A
~-aI+BKF) (A-aI)T
P= +HTH+-
r ' r For the disk D with radius r = 2 and center at - 6 + j 0 , the D-pole
. PB(r2R+BTPB)-'{(r2R+BTPB)K+KT(r'R+BTPB) assignment has been done for the several cases.
I ) Case 1: R = 0.1, Q = I
(A - CXZ)
-KTBTPBK-(r2R+BTPB)}(rZR+BTPB)-'BTP-
is satisfied, we may see that the closed-loop system can be found D-pole
assigned. Let the maximum eigenvalue of the matrix BTPB be X, then
the left-hand side of (32) may be rewritten as
2.433
31.14
4.4556
-31.14
X , ( A + B F ) : (-5.344,-5.827,-5.891
0.406
11.17 1 E D.
(r'R+h,I)K+KT(r?R+X,I)-X,KTK-(rZR+X,I)
+ ( K - I ) T ( X , , , I - B T P B ) ( K - I ) . (33) 3) Case 3: R = diag { 1 0 0 , loo}, Q = I
1
Since in (33) the following inequality:
6.427 2.49~
(K-I)T(X,,I-BTPB)(K-I)zO 28.45 -28.45 10.67
is satisfied, the rest of the terms are need to be considered. From R = hi(A + BF) : { - 5.333 kj2.078 x - 5.43} E D.
diag { Fl, f ~ ---
, F m } , and (33), it is known that the k,'s satisfying the
followjing: 4) Case 4: R = diag { 1000, lOOO}, Q = I
or
2(r'i,+X,,,)k,-X,kT-(r'i,+Xm)~0 (34) F= -
[
1.48x 10-4
28.44
6.43
-28.44
2.49~
10.67
1
Xi(A+BF): { - 5 . 3 3 + j 6 . 5 9 ~ 1 0 - ' , - 5 . 4 3 ) E D.
(34) '
Thus. it follows that the proposed algorithm gives D-pole assignment
give the inequality (32). since A,,, > 0 from the positive definity of P. and the pole locations within the specified disk D are dependent on the
Therefore, it follows that k, satisfying (34) is given by (29). matrices Q, R of the Riccati equation.
Remark: We may see that in (32) ai is going to approach to zero as r + In the case of the system (36) with perturbation K such as Fig. 2, to
0, which means that the gain margin for D-pole assignment becomes illustrate the variation of the closed-loop poles. as an example we take
smaller as r + 0.This case means that all poles of the closed-loop system Case 1) of the above D-pole assignment which has D-stability margin as
are assigned on a same point, so its robustness may be quit weak. 0.93 5 k, 5 I .08 and Q; 5 1.47". i = 1, 2. Fig. 3 shows the variations
Thus. in the above lemma. it is shown that the gain margin for D-pole of the closed-loop system poles for parameter perturbation K = k l where
assignment can be evaluated by (29). Similarly. w e may also see that the k is a scalar. As we see in Fig. 3 . it may be not ensured that the closed-
phase margin for D-pole assignment is obtained as loop poles are located \vithin a specified disk for perturbation beyond the
gain margin obtained by D-pole assignment.
V. CONCLUSIONS
IV. EXAMPLES This technical note presents an algorithm to assign all poles of the
closed-loop system in the disk D with the center on the real axis of the
To illustrate D-pole assignment, we take an example as the following: complex plane and arbitrary radius. Instead of locating all the poles of the
closed-loop system. the parameters of the desired disk are just needed to
be specified. and the pole locations within the specified disk D are only
dependent on the matrixes Q, R of the discrete type Riccati equation.
(36) The proposed pole assignment is said D-pole assignment, and its gain
and phase margin for the D-pole assignment are presented.
IEEE TRANSACTIOKS ON AUTOMATIC CONTROL, VOL. AC-32. NO.5 . MAY 1987 427
REFERENCES
F. Fallside. Ed.. Control System Design by Pole-ZeroAssignment. New
n. SYSTEMS AND DEHNITIONS
York: Academic.1977.
A. V. Bogachev. V . V . Grigorev. V . N . Drozdov.and A . N. Korovyakov,
We consider a h e a r system described by the state equations
"Analytic design of controls from root indicators," Automat., Remote Contr.,
vol.40, no. 8, pp. 1118-1123. 1979. X(t)=Ax(t)+Bu(t)+Dw(t);
B. C. Kuo, Digital Confrol Systems. New York: Holt, finehart, snd Winston,
1980.
T. Pappas, T. A. J. h u b . and N. R. Sandell. "On the numerical solution of the
discrete-time algebraic Riccati equation," IEEE Trans. Automat. Contr., vol.
AC-25, no. 4. pp. 631-641, 1980.
M . G. Safonov, Stability and Robustness of Multivariable Feedback Sys-
I. INTRODUCTION
When designing a control system, one often begins with a plant which for all E 12.That is, the H" norm of H(s) is less than or equal to y.
(J
is subject to external disturbances. A common design objective is to We now present a condition which can be used to test whether a given
reduce the effect of these disturbances to an acceptable level. For system is stabilizable with disturbance attenuation y. In order to test this
example. in the "disturbance decoupling problem" described in [l], state condition, the designer must specify two positive-definite matrices Q and
feedback control is used to ensure that the disturbances are completely R. It willbe shown in the sequel that the question of disturbance
decoupled from the output. However, for some systems, it may be attenuation is independent of the choice of Q and A .
impossible to reduce the effect of the disturbances below a certain Condition I: Let Q E B n X n and R E B m x m be given positive-definite
threshold value. Hence, the disturbance decoupling problem would have matrices and let the constant y > 0 be prespecified. Then the system (E)
no solution in this case. In this note. we present a procedure for designing is said to satisfy Condition 1 (with attenuation constant y) if there
a stabilizing state feedback control which will reduce the effect of the exists an E > 0 such that the Riccati equation
disturbances to a prespecified level.
The results of this note can also be applied to a class of H aoptimization 1 1 1
A ' P + P A - - PBR-'B'P+-PDD'P+-E'E+EQ=O (2.1)
problems in which the effect of the disturbances is minimized in a certain E Y Y
sense.Indeed, by applying our procedure with successively smaller
values of the prespecified disturbance level, we can construct a state has a positive-definite solution P.
feedback control which is arbitrarily close to the H" optimum. Thus, the Remark: Riccati equation (2.1) is of a type which also arises in linear
results of this note are closely related to recent results on the N" quadratic differential games. Furthermore, the process of decreasing E
optimization problem; see [2]. However, previous results onthe H" corresponds to a linear quadratic differential game with cheap control for
optimization problem have led to quite complicated design procedures. In the minimizing player; see [4]. Thus, we will see that there is a connection
contrast, the design procedure proposed in this note is a simple scheme between the H" optimization problem and linear quadratic differential
basedonthe algebraic Riccati equation. It follows that the algebraic games with cheap control.
Riccati equation which was previously associated with LQG optimization Theorem 2. I: Suppose that there existpositive-definite matrices Q
also plays a role in H" optimization. E W X n , R E L P and " ~a ~ constant e > 0 such that Riccati equation
(2.1) has a positive-definite solution. Then,givenany positive-
ManuscriptreceivedOctober 21. 1985;revised h'o\,ember 10, 1986. Paperrecom- definite matrices Q E P x n andl? E R m x mthere , existsa constant E*
mended by Past Associate Editor. S . P. Bhattacharyya.
The author is with theDepartment of Electrical and Electronic Engineering, University
of NewSouthWales.AustralianDefenceForceAcademy.Campbell,ACT 2600, ' Given two matrices N,1 i 4 E P " " , the notation N 2 M refers to the fact that the
Australia. -
matrix N M is positive semi-definite. Similarly, the notationN > M refers to the fact
IEEE Log Number 8613845. that the matrix h'- M is positivedefinite.