UDC 512.64
DOI: 10.18698/1812-3368-2022-4-39-49
REDUCTION OF BAND MATRICES IN LARGE DYNAMIC SYSTEMS CONTROLLABILITY AND OBSERVABILITY
N.E. Zubov1, 2 V.N. Ryabchenko2, 3
1 S.P. Korolev Rocket and Space Corporation Energia PJSC, Korolev, Moscow Region, Russian Federation
2 Bauman Moscow State Technical University, Moscow, Russian Federation
3 NTC FNC UES JSC, Moscow, Russian Federation
Abstract
An approach is proposed for linear stationary dynamical system with controllability and observability band matrices making it possible to simplify procedures for evaluating controllability and observability of this system. The obtained results are based on the fact that the controllability and observability criteria of a dynamic system are equivalent due to their required and sufficient properties; therefore, any transformations of one criterion not violating the conditions of necessity and sufficiency could be reduced to transformations in a sense equivalent to the initial transformations. The Popov — Bele-vich — Hautus transformations of the controllability and observability criteria were taken as a basis, and then results of such transformations were correctly extended to the band criteria. It was proved that the controllability and observability analysis of a linear stationary system with a large number of the state dimensions was reduced to studying the matrices rank of a much smaller size. The proposed approach is based on the existence condition for a numerical matrix of a certain rank of the nondegenerate matrices that satisfy certain transformations. The corresponding controllability and observability theorems for the stationary dynamical systems were provided. It was shown that for systems with one input and one output, the controllability and observability analysis was reduced to the analysis of scalars
Keywords
Linear stationary dynamical system, controllability and observability criteria, band matrices, matrix reduction
Received 14.12.2021 Accepted 04.04.2022 © Author(s), 2022
Introduction. Let us consider a linear dynamical system
x(t) = Ax(t) + Bu(t), y(t) = Cx(t),
(1)
where x(t) g M.n is the state vector; u(t) e W is the control vector; y(t) e Rm is the output vector. It is assumed that rank B = r, rank C = m and n » 1, i.e., the input and output matrices are having the corresponding full rank, and dimension of the system states (1) is a large number. Foreign literature introduces the concept of Large-Scale System for such system classes.
Various criteria are used for controllability and observability of the dynamical systems with representation in the state space (1). The most common are the Kalman criteria and the Popov — Belevich — Hautus tests [1-3]. According to these criteria, it is necessary and sufficient that the rank conditions for complete controllability and observability are fulfilled (1):
C ~~CA CA2
rank
B\AB\A2 B
An - rB
= n;
rank
CAn
= m;
V\g С : rank [ A - XIn \ B ] = n; VÀ, e С : rank
A -Xln
с
= n.
(2) (3)
Here In is the identity matrix of the size n x n; C is the set of complex numbers (complex plane).
Let us note that the V^e C condition in (2), (3) could be replaced without losing generality by condition V^eA( A), where A( A) = = | det (A - XiIn ) = 0} is the set of the A matrix eigenvalues.
In the early 2000s, the criteria for controllability and observability of a dynamic system (1) were proposed, where controllability and observability were determined based on the rank examination [3, 4]:
- controllability band matrix sized (n - r )(n - r +1) x n (n - r):
Cband ~
LA 0 0 H 0 _ i_ .
L LA 0 ••■1 0
0 L LA •■■1 0
0 0 L • 1~ • ■
: M '•. [LA
0 0 0 L
(4)
- observability band matrix sized n (n -m) x (n -m)(n -m +1):
Oband ~
AR R 0 0 0
0 AR R 0 0
0 0 AR
R 0
0 0 0 AR R
(5)
Here L, R are the (non-unique) solution matrices of the maximum rank [2, 5, 6] for the homogeneous equations:
LB = 0; (6)
CR = 0. (7)
The L and R matrices are also called the left and right zero divisors of the maximum rank [3-5, 7].
Proved in [3, 4] that for complete controllability of a dynamical system (1), it is necessary and sufficient that
rank Cband = n(n - r +1). (8)
Dually, it is necessary and sufficient for complete observability (1) that
rank Oband = n(n -m +1). (9)
Practical problems in analyzing controllability and observability of large dynamical systems (see, for example, [8-11]) often face situations, where dimension of the control and output vectors is comparable with the dimension of the state space. In this case, standard controllability tests for a system (1) are the high-dimensional ill-conditioned problems [12].
Recursive algorithms were proposed for band controllability and observability criteria (8), (9) with matrices (4), (5) in [12] that could significantly reduce the computational costs in large dynamic systems.
This work objective is to construct reductions that would make it possible to make a judgment on controllability and observability of a large dynamical system (1) based on studying the rank of matrices of a much smaller size.
Methodological basis of the work. For any D e Kpxk numerical matrix of the s rank, there are nonsingular matrices [3, 4, 7]:
ErL, [er|r].
satisfying the following transformation:
El
L
D [ ER\R ] =
0( p - Ф
0sX (k - s) 0(p - s) X (k - s)
(10)
Here 0i x j is the zero matrix of size i x j. Thus, transformation in the form expanded to blocks (10) decomposes into four relations:
ElDEr = Is ; LDER = 0(p - s) x s ; ElDR = 0s x (k - s); LDR = 0( p - s) x (k - s).
<r matrix of the r rank with transformation by type (10)
For system (1), the B e using matrix
is brought to the following form:
B+ ~L
(11)
°(n - r)>
Here B+ = El is the Moore — Penrose pseudoinverse matrix [13]; L is the maximum solution of the homogeneous equation (8).
Similarly, the C e Rmxn matrix of the m rank from (1) using matrix
[ с
+ !
R
where R is the maximum solution
is transformed to the form [ Im | 0m x (n - m)
to the homogeneous equation (7); C+ = ER is the Moore — Penrose pseudoinverse matrix [13].
The described controllability and observability criteria of the dynamical system (1) are equivalent due to required and sufficient properties [1, 4]. Therefore, any transformations of one criterion not violating the necessity and sufficiency conditions could be reduced to transformations that are in some sense equivalent to the original transformations. For example, if controllability of one system from the Popov — Belevich — Hautus criterion entails controllability of some other system, then this conclusion could be extended both to the Kalman criterion and to the band matrices criterion.
The present work takes as the basis the Popov — Belevich — Hautus transformations of the controllability and observability criteria [14]; results of such transformations are further correctly extended to the band criteria.
Matrix reduction in the controllability criteria. Using matrix (11), we transform the [ A - XIn | B ] matrix beam from criterion (2):
B+
T
B+(A-Xïn ) Ir
L (A-Xln ) 0(n - r) x r
[A-XI„ |B] =
Due to the matrix (11) nonsingularity, the following is obtained [13]: rank [ A - Xln | B ] = rank
B +(A-Xln ) Ir
L (A-Xln ) 0(n - r)xr
(12)
It follows from the structure of matrix (12) that the [ B+ ( A -XIn ) j Ir ] submatrix has the r rank for any X values. In this regard, in order to fulfill condition (2), it is necessary and sufficient that the rank of the L (A -XIn ) matrix beam satisfies the requirement [3, 4]:
VÀ, e C : rank L (A - XIn ) = n - r. Let us introduce into consideration the nonsingular matrix
[ L+ |B ], (13)
satisfying the following identity:
L [ L+ \B ] = [ In - r I 0(n - r) X r ] . Using (13), let us transform the L (A-XIn ) submatrix from (12): L (A-XIn )[ L0 |B ] = [ LAL+-XIn _ r |LAB ].
In this case, as in the previous one (12),
rank L (A -XIn ) = rank [LAL+ -XIn _ r | LAB]. (14)
To simplify, let us introduce the following notations:
A(1) = LAL+ ; B(1) = LAB;
n(1) = n -r; r(1) = rank Б(1).
(15)
In this case, let us assume that L(1) is the maximum solution to the homogeneous equation
L(1) B(1) = 0. (16)
Comparing the right parts of expressions (2) and (14) and taking into consideration equivalence of the considered criteria, validity of the following lemma is obtained.
Lemma 1. The dynamical system (1) is controllable, if and only if the following equivalent conditions are satisfied:
rank
B(i)
A(i)B(i)
(Ad))2 B(i) |-| (Ad))"(1)"Г(1) B(i)
= n(i),
VkeA ( A(i)): rank [ A(i) - MnW | B(i) ] = n(i)
(17)
(18)
C(i) =
band
L(i) A« 0 0 0
L(i) L(i) A(i) 0 0
0 L(i) L(i) A(i) 0
0 0 L(i)
L(i) A(i)
0 0 0 L(i)
J(n(D - r(1))(n(D - r(D +1) x n(D(n(1) - r(D)
(19)
where A(A(1) ) = A (LAL+^) is the set of eigenvalues of matrix A(1); A(1), B(1),
l(1) are the matrices satisfying (15), (16).
The transformations carried out resulted in reduction in the matrices size under the considered criteria. Thus, if we assume that n = 100, r = 40, r(1) = = 30 in system (1), then in (17) there is a controllability matrix of the 60 x 120 size (initial size is 100 x 2400), in (18) — a matrix beam of the 60 x 100 size (initial size is 100 x 140), and in (19) — a band matrix of the 930 x 1800 size (initial size is 3660 x 6000). The most significant decrease in size occurred in the controllability band matrix (almost by 4 times).
Let us consider transformation of the form (10) with respect to matrix B(1):
L(i)
Bi
E« I R(i)
Ir(i) 0r(i) x (r - r(i))
0(ni(i) - r(i)) x r(i) 0(n(i)- r(i)) X (r - r(i))
(20)
Applying transformation (20) to the [ A1 -XIn1 J B1 ] matrix and taking into account the right side of (20), the following is obtained:
E(i) el
L(i)
[A(i)-Und) |B(i)][eRR) |R(i)] =
eL] (A(i) -^In(i) )E(i E? (A(i) -^In(i) )R(i) Ir(i) 0r(i) X (r - r(i))
L(i) (A(i) -Ц щ )ER1 L(i) (A(i) -XIn(l) )R(i) 0(ni(i) - r(i)) X r« 0(n(i)-r(i)) X (r-r«)
(2i)
Results of the analysis (21) show that in this case the chain of rank conditions is valid:
В? (A(1)-4,1) ) ER E? ( A(1) -^I„(1) ) R(1) Ir(1) °r(1) X (r-r(1))
L(1) (A(1)-H,, )E$ L(1) (A(1) -Unm )R(1) V« - r(1)) X r(1) °(nœ - г(1)) x (r - r«)
rank
= r(1) + rank [ L(1) ( A(1) - W„(i) ) 4» | L(1) ( A(1) - ) R(1) ] = = r(1) + rank L(1) (A(1) -AJn(1) )[eR? | R(1)] = r(1) + rank L(1) (A(1) ). (22)
Structure of the L(1) ( A(1) -^In(i) ) matrix beam from (22) exactly corresponds to the L(A-XIn) matrix beam structure from (12). On the basis of its properties, Lemma 1 was formulated. This fact makes it possible to formulate another assertion.
Lemma 2. Dynamical system (1) is controllable, if and only if the following equivalent conditions are satisfied:
rank
B(2) i A(2)ß(2)
(A(2) )2 B(2) (A(2) )"(2) "r (2) B(2)
= n(2),
(23)
VXeA( A(2) ) : rank [ A(2) - Mn(2> | B(2) ] = n(2),
C(2) = band
L(2)A(2) 0 0 0
L2) L(2)A(2) 0 0
° L2) L(2)A(2) 0
0 0 L2)
l(2)A(2)
0 0 0 L2)
(n(2) - r(2) )(n(2) - r(2) +1) X n(2) (n(2) - r(2) )
where A(2) = L(1)AL(1)+; B(2) = L(1)A(1)B(1); n(2) = n(1) -r(1); r(2) = rankB(2); L2) is the maximum rank solution to homogeneous equation l(2)b(2) = 0; a(a(2)) = a(l(1)AL(1)+) is the set of eigenvalues of matrix a(2); l(1)+
is the Moore — Penrosepseudoinverse matrix.
Continuing the reasoning by induction, the following theorem is obtained, which validity was already proven.
Theorem 1. Dynamical system (1) is controllable, if and only if the following equivalent conditions are satisfied:
rank
B» |a»b» |(a» )2 B» |(a» f] "^ B» VX e Л ( A(i) ) : rank [ A(i) - AJn(0 | B(i) ] = n(i);
= n(i );
C
(i) = band
L(i )A(i) 0 0 0
L(i ) L(i) A(i) 0 0
0 L ) L(i )A(i ) 0
0 0 L(i)
L(i)A(i )
0 0 0 L(i)
- (n(i) - r(i))(n(i) - r(i) +1) X n(i)(n(i) - r(i))
where i = 1,2,3,... < n-r are the natural numbers; A(l) = L(i_1)AL(i"1) +; B(i) = L(i - 1)A(i" 1)B(i"1); n(i) = n(i -1) - r(i"1); r(i) = rank B(i); L(i) * 0 is the solution of the maximum rank of the homogeneous equation L(i)B(i) = 0; A(A(i)) is the set of eigenvalues of matrix A(i); L(i')+ is the Moore — Penrose
pseudoinverse matrix. At i = 1: A(0) = A, B(0) = B, n(0) = n, r(0) = r. The corollary obviously follows from Theorem 1.
Corollary 1. Dynamical system (1), where n > 1, u(t) e M1 is a scalar, is controllable if and only if the pair of scalars is controllable:
(A(n"B(n-) = (L(n"2)A(n"2)L(n"2)+, L(n-2)A(n"2)B(b"2)),
i.e., in (24), the scalar B
(n -1) = L(n - 2) A(n - 2)n(n - 2)
A(
)B(
^ 0.
Matrix reduction in the observability criteria. To solve the problem of analyzing the system observability (1), let us use matrix transformations (3) dualized to the matrix transformations (2) performed above. As a result, the following statement is obtained.
Theorem 2. Dynamical system (1) is observable, if and only if the following equivalent conditions are satisfied:
C(i)
rank
C(i) (A(i)
= m(i);
V^eA(A(i)): rank
A(i) - Ц,«) " C(i)""
= n(i);
O(i) = band
A(i)R(i) R(i) 0 0 0
0 A (i)R(i) R(i) 0 0
0 0 A (i)R(i)
R(i) 0
0 0 0 A(i )R(i) R(i)
(n(i) - m (0)(„(0_m(i) +1) x n(i)(n(i)- m(i))
where i = 1,2,3,...<n-m are the natural numbers; A(i) = R(~1) +AR(i"1);
B(i) = C(i" 1)aA(i" 1)R(i"1); n(i) = n(i -1) - m(i ~1); m(i) = rank C(i); R(i) * 0 is the
solution of the maximum rank of the homogeneous equation C(i)R(i) = 0; A (A(i))
is the set of eigenvalues of matrix A(i); R(i)+ is the Moore — Penrose pseudoinverse matrix. At i = 1:
A(0) = A, C(0) = C,
n(0) = n, m(0) = m.
Corollary 2. System (1), where n > 1, y(t) e M1, is the scalar and is observable, if and only if the pair of scalars is observable:
( A(n " D, C(n ~ D ) = ( R(n ~ 2)+ A(n " 2)R(n - 2),
4n - 2)?(n - 2)v(n - 2)
C ( 2) A( 2) R
i.e., in (36) the scalar C(n ~1) = C(n" 2)A(n" 2)R(n " 2) ^ 0.
Conclusion. It was proved that analysis of controllability and observability of a dynamical system with a large number of state dimensions by reduction
lies in studying the rank of matrices of much smaller size. For systems with one input and one output, the controllability and observability analysis are reduced to analyzing the scalars. The result obtained is of great importance and simplifies the analysis and synthesis procedure in control systems and the solution of the above mathematical problems in conjunction with other areas of research by the authors on the modal control methods [15-18] and the possibilities of its application, in particular, solving the algebraic equations [19].
REFERENCES
[1] Kailath T. Linear systems. Prentice Hall, 1980.
[2] Bernstein D. Matrix mathematics. Princeton Univ. Press, 2005.
[3] Zubov N.E., Mikrin E.A., Ryabchenko V.N. Matrichnye metody v teorii i praktike sistem avtomaticheskogo upravleniya letatel'nykh apparatov [Matrix methods in theory and practice of automated control systems for aircraft]. Moscow, BMSTU Publ., 2016.
[4] Gadzhiev M.G., Misrikhanov M.Sh., Ryabchenko V.N., et al. Matrichnye metody analiza i upravleniya perekhodnymi protsessami v elektroenergeticheskikh sistemakh [Matrix control and analysis methods for transient processes in electrical power systems]. Moscow, MPEI Publ., 2019.
[5] Wonham W.M. Linear multivariable control: a geometric approach. In: Stochastic Modelling and Applied Probability. New York, NY, Springer, 1985.
DOI: https://doi.org/10.1007/978-1-4612-1082-5
[6] Zybin E.Yu., Misrikhanov M.Sh., Ryabchenko V.N. Minimal'naya parametrizatsiya resheniy lineynykh matrichnykh uravneniy. V kn.: Sovremennye metody upravleniya mnogosvyaznymi sistemami. Vyp. 2 [Minimum parametrization of solutions for linear systems of matrix equations. In: Modern Control Methods for Multicoupling Systems. Iss. 2]. Moscow, Energoatomizdat Publ., 2003, pp. 191-202 (in Russ.).
[7] Taufer J. Randwertprobleme für Systeme von Linearen Differentialgleichungen. Praha, Academia, 1973.
[8] Barinov V.A., Sovalov S.A. Rezhimy energosistem. Metody analiza i upravleniya [Regimes of power systems. Analysis and control methods]. Moscow, Energoatomizdat Publ., 1990.
[9] Gusseynov F.G. Uproshchenie raschetnykh skhem elektricheskikh system [Computational scheme simplification for electrical systems]. Moscow, Energiya Publ., 1978.
[10] Christensen G.S., El-Hawary M.E., Soliman S.A. Optimal control applications electric power systems. Series Mathematical Concepts and Methods in Science and Engineering. New York, NY, Springer, 1987. DOI: https://doi.org/10.1007/978-1-4899-2085-0
[11] Kundur P. Power system stability and control. McGraw-Hill, 1994.
[12] Zybin E.Yu., Misrikhanov M.Sh., Ryabchenko V.N. Recursive controllability and observability tests for large dynamic systems. Autom. Remote Control, 2006, vol. 67, no. 5, pp. 783-795. DOI: https://doi.org/10.1134/S0005117906050109
[13] Gantmakher F.R. Teoriya matrits [Matrix theory]. Moscow, Nauka Publ., 1968.
[14] Misrikhanov M.Sh., Ryabchenko V.N. Reduction of state dimensions at analysis of controllability and observability of linear models for power systems. Izvestiya TRTU, 2005, no. 11, pp. 42-54 (in Russ.).
[15] Mikrin E.A., Zubov N.E., Lapin A.V., et al. Analytical formula of calculating a controller for linear SIMO-system. Differentsial'nye uravneniya i protsessy upravleniya [Differential Equations and Control Processes], 2020, no. 1 (in Russ.).
Available at: https://diffjournal.spbu.ru/pdf/mikrin1.pdf
[16] Mikrin E.A., Ryabchenko V.N., Zubov N.E., et al. Analysis and synthesis of dynamic MIMO-system based on band matrices of special type. Differentsial'nye uravneniya i protsessy upravleniya [Differential Equations and Control Processes], 2020, no. 2 (in Russ.). Available at: https://diffournal.spbu.ru/RU/numbers/2020.2/article.1.1.html
[17] Zubov N.E., Mikrin E.A., Misrikhanov M.Sh., et al. Output control of the spectrum of a descriptor dynamical system. Dokl. Math., 2016, vol. 93, pp. 259-261.
DOI: https://doi.org/10.1134/S106456241603008X
[18] Zubov N.E., Lapin A.V., Mikrin E.A., et al. Output control of the spectrum of a linear dynamic system in terms of the Van der Woude method. Dokl. Math., 2017, vol. 96, pp. 457-460. DOI: https://doi.org/10.1134/S1064562417050179
[19] Mikrin E.A., Zubov N.E., Efanov D.E., et al. Superfast iterative solvers for linear matrix equations. Dokl. Math., vol. 98, pp. 444-447.
DOI: https://doi.org/10.1134/S1064562418060145
Zubov N.E. — Dr. Sc. (Eng.), Professor, Post-Graduate Education, S.P. Korolev Rocket and Space Corporation Energia PJSC (Lenina ul. 4A, Korolev, Moscow Region, 141070 Russian Federation); Professor, Department of Automatic Control Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation).
Ryabchenko V.N. — Dr. Sc. (Eng.), Assoc. Professor, Professor, Department of Automatic Control Systems, Bauman Moscow State Technical University (2-ya Baumanskaya ul. 5, str. 1, Moscow, 105005 Russian Federation); Senior Technologist, NTC FNC UES JSC (Kashirskoe shosse 22, korp. 3, Moscow, 115201 Russian Federation).
Please cite this article as:
Zubov N.E., Ryabchenko V.N. Reduction of band matrices in large dynamic systems controllability and observability. Herald of the Bauman Moscow State Technical University, Series Natural Sciences, 2022, no. 4 (103), pp. 39-49. DOI: https://doi.org/10.18698/1812-3368-2022-4-39-49