Applied Mathematics and Computation
|
|
- Louisa Thornton
- 5 years ago
- Views:
Transcription
1 Applied Mathematics and Computation 212 (2009) Contents lists available at ScienceDirect Applied Mathematics and Computation journal homepage: Solutions to a family of matrix equations by using the Kronecker matrix polynomials Bin Zhou a, *, Zhao-Yan Li b, Guang-Ren Duan a, Yong Wang b a Center for Control Theory and Guidance Technology, Harbin Institute of Technology, P.O. Box 416, Harbin , Heilongjiang, PR China b Department of Mathematics, Harbin Institute of Technology, Harbin , PR China article info abstract Keywords: Generalized Sylvester matrix equations Kronecker matrix polynomials Coprime Closed form solutions Linear system theory Closed form solutions to a family of generalized Sylvester matrix equation in form of P / A ixf i þ P w B kyf k ¼ P u E jrf j are given by using the so-called Kronecker matrix polynomials. It is found that the structure of the solutions is independent of the orders /; w and u. This type of uniform closed form solutions includes our early results as special cases. The results provide great convenience to the computation and analysis of the solutions to this class of equations, and can perform important functions in many analysis and design problems in linear systems. Ó 2009 Elsevier Inc. All rights reserved. 1. Introduction When dealing with the analysis and design problems associated with the linear system _xðtþ ¼AxðtÞþBuðtÞ where A and B are known matrices with appropriate dimensions, a matrix equation in the form of AX XF ¼ BY þ R where F and R are two known matrices, and ðx; YÞ is a pair of matrices to be determined, are encountered. The homogeneous case of the above equation is AX XF ¼ BY. These two first-order generalized Sylvester matrix equations are widely used in some control problems such as pole and eigenstructure assignment design [15,16,24], Luenberger-type observer design [17,19], and robust fault detection [5,18]. When dealing with the analysis and design problems associated with the descriptor linear system E_x ¼ Ax þ Bu where A; E and B are known matrices having appropriate dimensions, we often encounter a matrix equation of the form AX EXF ¼ BY þ R where the matrices F and R are known with proper dimensions, and the matrix pair ðx; YÞ is to be determined. In the case of R ¼ 0, the above mentioned nonhomogeneous equation becomes AX EXF ¼ BY which is closely related with pole placement, eigenstructure assignment [1], observer design [6], and some other problems (see [10] for a detailed introduction). Many control problems, such as pole assignment [23 25] and eigenstructure assignment [4,26], for the second-order linear system M x þ D_x þ Kx ¼ Bu where M; D; K and B are known matrices with appropriate dimensions, are closely related with a matrix equation of the form MVF 2 þ DVF þ KV ¼ BW þ R where F and R are given matrices, and ðx; YÞ is a pair of matrices to be determined. This type of equation is called the second-order nonhomogeneous generalized Sylvester matrix equation. When R ¼ 0, it becomes the homogeneous one MVF 2 þ DVF þ KV ¼ BW. It can be shown that certain control problems, such as pole assignment, eigenstructure assignment and observer design [3], of the high-order linear system P m A ix ðiþ ¼ Bu where A i ; i ¼ 0; 1;...; m and B are known matrices, are closely related with * Corresponding author. addresses: binzhou@hit.edu.cn, binzhoulee@163.com (B. Zhou), zhaoyanlee@gmail.com (Z.-Y. Li) /$ - see front matter Ó 2009 Elsevier Inc. All rights reserved. doi: /j.amc
2 328 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) a matrix equation of the form [3] P m A ixf i ¼ BY þ R in which F and R are some known matrices having proper dimensions, and ðx; YÞ is a pair of matrices to be determined. Such a type of equation is called the mth order nonhomogeneous generalized Sylvester matrix equation [3]. Its homogeneous part is P m A ixf i ¼ BY. Obviously, the above mth order generalized Sylvester equation includes both the above mentioned first-order equations and second-order equations as special cases. A large amount of prior work has been devoted to solving the generalized nonhomogeneous and homogeneous Sylvester matrix equations given above with F in the Jordan or companion form. See [1,2,9,20,21]. However, the restriction that F is in Jordan canonical form is an obvious disadvantage. Because both the existing result itself and its proof in the literature (see, e.g., [1]) heavily rely on the concrete structure of F as a Jordan canonical form, one can not obtain the solutions of the firstorder generalized Sylvester matrix equation for general matrix F by using the existing result straightway. To deal with the general matrix F, one should transform F into its Jordan canonical form beforehand, and then transform the first-order generalized Sylvester matrix equations into an equivalent one, in which the Jordan form appears. The solutions obtained by using this method are very complicated for involving the computation of Jordan canonical form, and will not be convenient for use in practice. Therefore, it is useful and interesting to give some complete and explicit solution by using the general matrix F itself directly. For these reasons, the authors have done some work on this problem recently. In detail, we have proposed in [10] for the first time a complete parametric solutions to the equation AX EXF ¼ BW with F an arbitrary square matrix and then extended it to the second-order case in [8]. As a result, it is found that solutions to equations AX EXF ¼ BW and MVF 2 þ DVF þ KV ¼ BW are closely related with some matrix polynomials associated with the coefficient matrices and have the same structure that is independent of the coefficient matrices and the order of the equation (see [8,10] for details). Therefore, one aim of this paper is to provide a methodology to unify the solutions to all these equations. The results are obtained by using a different method from that used in [8,10]. Another aim of this paper is to extend this family of generalized Sylvester matrix equation to a more general case that the unknown variables are more than two. It will be found that the result can also be given in a very similar way. The results presented here also generalize those in [11] in some aspects. Especially, the so-called Kronecker matrix polynomial given in this paper seems more convenient to use than the so-called generalized Sylvester mapping presented in that paper. Finally, we should point out that matrix equation is a very hot topic and many relating results can be found in the literature in recent years. For related references, see [13,14,22,27,28] and the references given there. The remainder of this paper is organized as follows. In Section 2, we introduce the so-called Kronecker matrix polynomials whose properties are studied. By using the Kronecker matrix polynomials, closed form solutions to a family of generalized Sylvester matrix equations are presented in Section 3. Some further extensions are given in Section 4. Section 5 concludes the paper. Notations: In this paper, we use to denote the Kronecker product. For an m n matrix R ¼½r ij Š, the so-called column stretching function csðþ and the row stretching function rsðþ are respectively defined as csðrþ ¼½r 11 r m1 r 1n r mn Š T ; rsðrþ ¼½r 11 r 1n r m1 r mn Š: For three matrices M; X and N with appropriate dimensions, we have the following well-known results related with the stretching operations: csðmxnþ ¼ðN T MÞcsðXÞ; rsðmxnþ ¼rsðXÞðM T NÞ: ð1þ ð2þ Furthermore, we use rðaþ to denote the eigenvalue set of matrix A, col½a 1 ; A 2 ;...; A s Š¼col½A i Š s i¼1 to denote a block column matrix with the ith matrix A i, and diag½d i Š s i¼1 to denote a diagonal matrix with the ith element being d i, respectively. 2. The Kronecker matrix polynomials For a matrix polynomial AðsÞ ¼ Xx A i s i ; A i 2 R nq ; i ¼ 0; 1;...; x; where x P 0 is a given integer, and a square matrix F 2 R pp. we define the so-called Kronecker matrix polynomial as A ðfþ ¼ Xx ðf T Þ i A i : Obviously, A ðfþ is produced from AðsÞ by the substitution ðs; Þ $ ðf T ; Þ. This kind of matrix polynomial was initially proposed in our early work [10] (see Eq. (11) and Eq. (17) in [10]) to study the matrix equation (13) given later and some of its properties have been reported there (Lemma 1 in [10]). A different definition for a similar kind of matrix polynomials from the mapping point of view was first given in [30]. The Kronecker matrix polynomial has the following properties shown in Lemmas 1 3 whose proofs are simple and omitted here.
3 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) Lemma 1. Let AðsÞ 2R nq ½sŠ; BðsÞ 2R qr ½sŠ and AðsÞBðsÞ ¼CðsÞ. Then C ðfþ ¼A ðfþb ðfþ: Let AðsÞ; BðsÞ 2R nq ½sŠ and AðsÞþBðsÞ ¼DðsÞ. Then ¼A ðfþþb ðfþ: ð3þ ð4þ Lemma 2. Let A k ðsþ 2R nq k ½sŠ; B k ðsþ 2R q k r ½sŠ; k ¼ 1; 2;...; m, and Then X m k¼1 X m k¼1 A k ðsþb k ðsþ ¼CðsÞ: ða k Þ ðfþðb k Þ ðfþ ¼C ðfþ: It follows from Lemma 1 that if AðsÞ 2R qq ½sŠ and BðsÞ 2R qq ½sŠ satisfy AðsÞBðsÞ ¼I q, i.e., AðsÞ (or BðsÞ) is invertible, then A ðfþb ðfþ ¼I q F 0 ¼ I qp which implies that A ðfþ (or B ðfþ) is also invertible. Using this fact, we can get the following result. Lemma 3. Let AðsÞ 2R nq ½sŠ be an arbitrary matrix polynomial, PðsÞ 2R nn ½sŠ and QðsÞ 2R qq ½sŠ be two unimodular matrix polynomials and satisfy the Smith normal factorization PðsÞAðsÞQðsÞ ¼ DðsÞ 0 ; 0 0 where DðsÞ ¼diag½d i ðsþš r i¼1, and d 1ðsÞjjd r ðsþ. Then rankða ðfþþ ¼ rankðþ ¼ Xr i¼1 rankðd i ðfþþ: The following definition is directly adopted from [10]. Definition 1. Let F 2 R pp be an arbitrary matrix. Then 1. A pair of polynomial matrices NðsÞ 2R nr ½sŠ and DðsÞ 2R qr ½sŠ is said to be F-right coprime if rank NðkÞ ¼ r; 8k 2 rðfþ: ð5þ DðkÞ 2. A pair of polynomial matrices HðsÞ 2R mn ½sŠ and LðsÞ 2R mq ½sŠ is said to be F-left coprime if rank½hðkþ LðkÞ Š ¼ m; 8k 2 rðfþ: ð6þ Remark 1. The definition of the F-right coprimeness (F-left coprimeness) of two matrix polynomials NðsÞ and DðsÞ can be easily extended to a series of matrix polynomials having the same columns (rows) and a single matrix polynomial. The following proposition was first established in [10]. But the proof given there is very complicated. Here a new simple and alterative proof is given by using Lemma 3. Proposition 1. Let F 2 R pp be an arbitrary matrix. Then 1. A pair of polynomial matrices NðsÞ 2R nr ½sŠ and DðsÞ 2R qr ½sŠ is F-right coprime if and only if rank N ðfþ ¼ rank½w ðfþš ¼ rp; where WðsÞ ¼col½NðsÞ; DðsÞŠ; 2. A pair of polynomial matrices HðsÞ 2R mn ½sŠ and LðsÞ 2R mq ½sŠ is F-left coprime if and only if rank½h ðfþ L ðfþ Š ¼ rank½m ðfþš ¼ mp; ð7þ where MðsÞ ¼½HðsÞ; LðsÞŠ. Proof. We only show the first item since the second item can be shown similarly. First, the relation rank N ðfþ ¼ rank½w ðfþš
4 330 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) can be easily obtained by changing some rows of the associated matrix. See [10]. Second, we show rank½w ðfþš ¼ rp. Let TðsÞ ¼diag½t i ðsþš / i¼1 2 R// ½sŠ; / 6 r be the Smith canonical form of the matrix polynomial WðsÞ, i.e., there exist two unimodular polynomial matrices PðsÞ and QðsÞ such that PðsÞWðsÞQðsÞ ¼ TðsÞ According to Lemma 3, rank½w ðfþš ¼ rp is equivalent to : ð8þ rank½t ðfþš ¼ X/ rankðt i ðfþþ ¼ rp: i¼1 Furthermore, according to (8) and (5), NðsÞ and DðsÞ are F-right coprime if and only if ð9þ rankðwðkþþ ¼ X/ rankðt i ðkþþ ¼ r; 8k 2 rðfþ; ð10þ i¼1 which implies that / ¼ r. The equivalence between (9) and (10) can be easily observed since rðt i ðfþþ ¼ ft i ðkþjk 2 rðfþg; i ¼ 1; 2;...; r: This completes the proof. 3. Solutions to a family of matrix equations 3.1. The homogeneous generalized Sylvester matrix equations In this subsection, we consider the following general homogeneous generalized Sylvester matrix equation X / A i XF i þ Xw B k YF k ¼ 0; where A i 2 R nn ; i ¼ 0; 1;...;/ and B i 2 R nr ; i ¼ 0; 1;...; w; F 2 C pp are some known matrices and X 2 R np ; Y 2 R rp are matrices to be determined. Clearly, this equation includes the following equations mentioned in Section 1 as special cases: AX XF ¼ BY; ð12þ AX EXF ¼ BY; ð13þ MXF 2 þ DXF þ KX ¼ BY; ð14þ X m A i XF i ¼ BY: Define a pair of polynomial matrices AðsÞ and BðsÞ as follows: AðsÞ ¼ X/ A i s i 2 R nn ½sŠ; BðsÞ ¼ Xw B i s i 2 R nr ½sŠ: Putting csðþ on both sides of (11) and using (2) gives csðxþ ½A ðfþ B ðfþ Š ¼ 0: csðyþ Let S and fs be the sets of the solutions of (11) and (17), respectively, i.e., ( ) X S ¼ : X/ A i XF i þ Xw B k YF k ¼ 0 ; Y and fs ¼ fw : ½A ðfþ B ðfþ Šw ¼ 0g: ð18þ Then the map s : S! fs defined by X s : # csðxþ ; ð19þ Y csðyþ is a linear isomorphism. Therefore, we need only to consider the Eq. (17). ð11þ ð15þ ð16þ ð17þ
5 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) Proposition 2. Let fs be given by (18). Then dimðfsþ ¼rp if and only if AðsÞ and BðsÞ are F-left coprime. Proof. Obviously, dimðfsþ ¼rp if and only if rank½a ðfþ B ðfþ Š ¼ np: ð20þ According to Proposition 1, Eq. (20) holds if and only if AðsÞ and BðsÞ are F-left coprime. This completes the proof. Assume that there exists a pair of matrix polynomials NðsÞ ¼ P x N is i 2 R nr ½sŠ and DðsÞ ¼ P x D is i 2 R rr ½sŠ, where x P 0 is some integer, such that AðsÞNðsÞþBðsÞDðsÞ ¼0: ð21þ Theorem 1. Suppose that AðsÞ and BðsÞ are F -left coprime and NðsÞ and DðsÞ are F-right coprime. Then N ðfþ fs ¼ csðzþ : Z 2 C rp Proof. First, since NðsÞ and DðsÞ satisfy (21), by using Lemma 2 we have N ðfþ A ðfþn ðfþþb ðfþ ¼½A ðfþ B ðfþ Š ¼ 0; which implies that w ¼ N ðfþ cs ðzþ; Z 2 C rp ; : ð22þ satisfies the Eq. (17). Second, since dimðzþ ¼rp, the relation dimðfsþ ¼rp holds if and only if rank N ðfþ ¼ rp: According to Proposition 1, the above equation is equivalent to the F-coprimeness of the matrix polynomials NðsÞ and DðsÞ. The proof is completed. In view of (19), the equation in (1) and Theorem 1, we clearly have the following corollary. Corollary 1. The solution space of (11) can be explicitly expressed by 82 P x 3 9 >< N i ZF i >= S ¼ 6 4 P x 7 D i ZF i 5 : Z 2 Crp : ð23þ >: >; Remark 2. The expression (23) was initially established in [10] for the first-order Eq. (13) and then extended to secondorder case (14) in [7,8]. It follows from Corollary 1 that this type of solutions also holds for the more general case (11). That is to say, the structure of the solution is independent of the order of the equation The nonhomogeneous generalized Sylvester matrix equations We now consider the nonhomogeneous extension of (11), i.e., a linear equation in the form of X / A i XF i þ Xw B k YF k ¼ Xu E j RF j ; where AðsÞ 2R nn ½sŠ and BðsÞ 2R nr ½sŠ are given by (16), EðsÞ ¼ Xu E j s j 2 R nq ½sŠ; ð24þ ð25þ is a given matrix polynomial, R 2 C qp and F 2 C pp are known matrices, and X 2 R np ; Y 2 R rp are matrices to be determined. Solution to a matrix equation with the right-hand-side similar to that of Eq. (24) was discussed in [29]. Clearly, the following matrix equations mentioned in Section 1 are special cases of (24):
6 332 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) AX XF ¼ BY þ R; AX EXF ¼ BY þ R; MXF 2 þ DXF þ KX ¼ BY þ R; X m A i XF i ¼ BY þ R: Let R be the set of the solutions to (24), that is, ( ) X R ¼ : X/ A i XF i þ Xw B k YF k ¼ Xu E j RF j : Y Similarly, Eq. (24) can be equivalently rewritten as csðxþ ½A ðfþ B ðfþ Š ¼ E ðfþcsðrþ: csðyþ Also, we denote ð26þ ð27þ ð28þ ð29þ ð30þ er ¼ ft : ½A ðfþ B ðfþ Št ¼ E ðfþcsðrþg: Assume that there exists a pair of matrix polynomials UðsÞ ¼ P u U is i 2 R nq ½sŠ and VðsÞ ¼ P u V is i 2 R rq ½sŠ, where u P 0is an integer, such that AðsÞUðsÞþBðsÞVðsÞ ¼EðsÞ: ð31þ Lemma 4. Let UðsÞ and VðsÞ satisfy (31). Then er 3 t s ¼ U ðfþ csðrþ: V ðfþ ð32þ Proof. Since UðsÞ and VðsÞ satisfy (31), using Lemma 2, we have U ðfþ ½A ðfþ B ðfþ Š csðrþ ¼ðA ðfþu ðfþþb ðfþv ðfþþcsðrþ ¼E ðfþcsðrþ; V ðfþ which implies that t s given in (32) satisfies the Eq. (30). This ends the proof. It follows from Lemma 4 that t s is a special solution to the nonhomogeneous Eq. (24). Combining Theorem 1 and Lemma 4 produces the following theorem. Theorem 2. Let AðsÞ and BðsÞ be F-left coprime, NðsÞ and DðsÞ satisfying (21) be F-right coprime, and UðsÞ and VðsÞ satisfy (31). Then N ðfþ U ðfþ csðzþ er ¼ : Z 2 C rp : V ðfþ csðrþ Similar to Corollary 1, we have the following result. Corollary 2. Let AðsÞ; BðsÞ; NðsÞ; DðsÞ; UðsÞ and VðsÞ be stated as in Theorem 2. Then Px >< N i ZF i þ Pu U i RF i >= R ¼ 6 P x 7 4 D i ZF i þ Pu V i RF i 5 : Z 2 Crp : >: >; It follows from Theorems 1 and 2 that the pair of matrix polynomials ðnðsþ; DðsÞÞ satisfying (21) and the pair of matrix polynomials ðuðsþ; VðsÞÞ satisfying (31) are crucial for solving the Eqs. (11) and (24). The following result gives a method to determine these four matrix polynomials. A similar result can be found in [11]. Proposition 3. Assume that the matrix polynomials AðsÞ 2R nn ½sŠ and BðsÞ 2R nr ½sŠ satisfy rank½aðsþ BðsÞ Š¼n; 8s 2 C: ð33þ Then there exists two unimodular matrices PðsÞ and QðsÞ with appropriate dimensions satisfying PðsÞ½AðsÞ BðsÞ ŠQðsÞ ¼½0 I n Š:
7 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) Denote QðsÞ ¼ Q 11ðsÞ Q 12 ðsþ ; ð34þ Q 21 ðsþ Q 22 ðsþ where Q 21 ðsþ 2R rr ½sŠ. Then the polynomial matrices NðsÞ; DðsÞ; UðsÞ and VðsÞ satisfying (21) and (31) can be chosen as NðsÞ ¼Q 11 ðsþ; UðsÞ ¼Q 12 ðsþpðsþeðsþ; DðsÞ ¼Q 21 ðsþ; VðsÞ ¼Q 22 ðsþpðsþeðsþ: ð35þ Furthermore, the polynomial matrices NðsÞ and DðsÞ are right coprime (thus F-right coprime for arbitrary F). Remark 3. The rank condition (33) has a systematic meaning when the generalized Sylvester matrix Eq. (11) reduces to the special cases as (12) (15). For instance, condition (33) reduces to rank½si A BŠ ¼ n; 8s 2 C ð36þ for Eqs. (12) and (26), reduces to rank½se A BŠ ¼ n; 8s 2 C ð37þ for Eqs. (13) and (27), and reduces to rank Ms 2 þ Ds þ K B ¼ n; 8s 2 C ð38þ for Eqs. (14) and (28). Obviously, in these cases, conditions (36) (38) represent the controllability property of the corresponding linear systems [7,8]. Remark 4. Theorem 2 also generalizes the results in [11] where the right hand side of (24) is a simple constant matrix R. 4. Some extensions 4.1. The generalized Sylvester-observer matrix equation The generalized Sylvester-observer matrix equation has the form X / F i XA i þ Xw F k YC k ¼ 0; where A i 2 R nn ; i ¼ 0; 1;...; /; C k 2 R mn ; k ¼ 0; 1;...; w; F 2 C pp are known matrices and X 2 R pn ; Y 2 R pm are matrices to be determined. Similar to (16), we denote AðsÞ ¼ X/ A i s i 2 R nn ½sŠ; CðsÞ ¼ Xw C i s i 2 R mn ½sŠ: Such a matrix Eq. (39) is the dual form of (11) and has wide applications in observer design problems [12]. Obviously, Eq. (39) includes several equations, such as TA FT ¼ LC with T and L unknown studied in [20,21], the equation F 2 TM þ FTD þ TK þ FLR þ LQ ¼ 0, with T and L unknown studied in [12], as special cases. Also we denote ( ) F i XA i þ Xw F k YC k ¼ 0 : O ¼ ½X YŠ : X/ By using (2), Eq. (39) can be equivalently rewritten as A ðfþ ½rsðXÞ rsðyþ Š ¼ 0; C ðfþ which is clearly the dual form of (17). Theorem 3. Let AðsÞ and CðsÞ be F-right coprime, HðsÞ ¼ P - H is i 2 R mn ½sŠ and LðsÞ ¼ P - L is i 2 R mm ½sŠ where - P 0 is an integer, satisfying HðsÞAðsÞþLðsÞCðsÞ ¼0; be F-left coprime. Then P- O ¼ F i P- ZH i F i ZL i : Z 2 C pm : ð39þ
8 334 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) Remark 5. The following equation X / F i XA i þ Xw F k YC k ¼ Xu F j RE j is the dual form of (24). Thus corresponding results to Theorem 2 and Corollary 1 can be easily obtained Matrix equations with more than two unknowns We now consider a more general case in the following form where X h X x k k¼1 A ki X k F i ¼ 0; ð40þ A k ðsþ ¼ Xx k A ki s i 2 R nr k ; r k P 1; k ¼ 1; 2;...; h ð41þ i¼1 are some given matrix polynomials and F 2 R pp is a given square matrix. For instance, when h ¼ 2 and A 1 ðsþ ¼AðsÞ; A 2 ðsþ ¼ B, then (40) reduces to (15). The set of solutions to (40) is denoted by ( S h ¼ col½x i Š h i¼1 : Xh X x ) k A ki X k F i ¼ 0 : k¼1 Assume that A k ðsþ; k ¼ 1; 2;...; h given by (41) are F-left coprime, i.e., rank½a 1 ðkþ A 2 ðkþ A h ðkþ Š ¼ n; 8k 2 rðfþ; and there exists a series of matrix polynomials B i ðsþ 2R n k l ; l ¼ P h r i n; k ¼ 1; 2;...; h such that X h k¼1 A i ðsþb i ðsþ ¼0: Theorem 4. If the matrix polinomials B k ðsþ ¼ P wk B kis i ; k ¼ 1; 2;...; h are F-right coprime, then ( ) S h ¼ col½x i Š h i¼1 : X i ¼ Xwi B ik ZF k ; Z 2 C lp : We can consider solutions to the following family of generalized Sylvester matrix equation with more umknowns: X h X x k k¼1 A ki X k F i ¼ Xu E j RF j ; ð42þ where A k ðsþ; k ¼ 1; 2;...; h are given by (41), E(s) is given by (25) and R is a given matrix. In fact, the equation in (42) can be written as the following simple and unified form X x A i XF i ¼ Xu E j RF j ; i¼1 ð43þ where X ¼ col½x i Š h i¼1 and AðsÞ ¼½A 1 ðsþ A 2 ðsþ...a h ðsþš, Xx A i s i 2 R nr ; with x ¼ max 16k6h fx k g; r ¼ P h r k and A i ðsþ; i ¼ 1; 2;...; h given by (41). Then by using a similar technique used in Section 3.2, we can propose the following result without a proof. Theorem 5. Assume that AðsÞ is F-left coprime. Let NðsÞ 2R rðr nþ and MðsÞ 2R rq be such that AðsÞNðsÞ ¼0; AðsÞMðsÞ ¼EðsÞ: Moreover, assume that NðsÞ is F-right coprime and NðsÞ ¼ X - N i s i ; MðsÞ ¼ X p M i s i ;
9 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) where - and p are two integers. Then the solution space of Eq. (43) can be characterized as ( ) S A ¼ : X - N i ZF i þ Xp M i RF i : Z 2 C ðr nþp : 5. Conclusions This paper deals with a family of generalized Sylvester matrix equations in the form of (11) and (24), which include the first-order and second-order generalized Sylvester matrix Eqs. (12) (15) and (26) (29) as special cases, by using the Kronecker matrix polynomials. This paper has achieved the following: (A): The solutions to this class of equations are unified by using the so-called Kronecker matrix polynomials which gives a deep insight to the structure and the solution of this family of equations. The solutions also have a very simple, neat and elegant structure, and (B): All of the solutions to this family of generalized Sylvester matrix equations can be obtained without any assumption on matrix F and R as soon as some matrix polynomials associated with the equation are obtained. Acknowledgements The work of Bin Zhou and Guang-Ren Duan was partially supported by the Major Program of National Natural Science Foundation of China under Grant No and Program for Changjiang Scholars and Innovative Research Team in University. The work of Zhao-Yan Li and Yong Wang was supported by the National Natural Science Foundation of China (Grant No ) and the Natural Science Foundation of Heilongjiang Province (Grant No ). The authors would like to thank the reviewers for their helpful comments which have helped them in improving the quality of this paper. References [1] G.R. Duan, Solutions to matrix equation AV þ BW ¼ VF and their application to eigenstructure assignment in linear systems, IEEE Transactions on Automatic Control 38 (2) (1993) [2] G.R. Duan, The solution to the matrix equation AV þ BW ¼ EVJ þ R, Applied Mathematics Letters 17 (2004) [3] G.R. Duan, Parametric eigenstructure assignment in high-order linear systems, International Journal of Control, Automation and Systems 3 (3) (2005) [4] G.R. Duan, Two parametric approaches for eigenstructure assignment in second-order linear systems, Journal of Control Theory and Applications 1 (2003) [5] G.R. Duan, G.P. Liu, Complete parametric approach for eigenstructure assignment in a class of second-order linear systems, Automatica 38 (4) (2002) [6] G.R. Duan, G.P. Liu, S. Thompson, Eigenstructure assignment design for proportional-integral observers: continuous-time case, IEE Proceedings: Control Theory & Applications 148 (3) (2001) [7] G.R. Duan, B. Zhou, Solution to the equation MVF 2 þ DVF þ KV ¼ BW: the nonsingular case, in: Proceedings of the Second IASTED International Multi- Conference on Automation, Control, and Information Technology, 2005, pp [8] G.R. Duan, B. Zhou, Solution to the second-order Sylvester matrix equation MVF 2 þ DVF þ KV ¼ BW, IEEE Transactions on Automatic Control 51 (5) (2006) [9] B. Zhou, G.R. Duan, An explicit solution to the matrix equation AX XF ¼ BY, Linear Algebra and Its Applications 402 (2005) [10] B. Zhou, G.R. Duan, A new solution to the generalized Sylvester matrix equation AV EVF ¼ BW, Systems & Control Letters 55 (3) (2006) [11] B. Zhou, G.R. Duan, On the generalized Sylvester mapping and matrix equations, Systems & Control Letters 57 (3) (2008) [12] B. Zhou, G.R. Duan, Parametric approach for the normal Luenberger function observer design in second-order linear systems, in: Proceedings of the 45th IEEE Conference on Decision and Control, 2006, pp [13] Y. Lei, A.-P. Liao, A minimal residual algorithm for the inconsistent matrix equation AXB ¼ C over symmetric matrices, Applied Mathematics and Computation 188 (1) (2007) [14] A. Kaabi, A. Kerayechian, F. Toutounian, A new version of successive approximations method for solving Sylvester matrix equations, Applied Mathematics and Computation 186 (1) (2007) [15] K.R. Gavin, S.P. Bhattacharyya, Robust and well-conditioned eigenstructure assignment via Sylvester s equation, Optimal Control Application and Methods 4 (1983) [16] B.H. Kwon, M.J. Youn, Eigenvalue-generalized eigenvector assignment by output feedback, IEEE Transactions on Automatic Control 32 (5) (1987) [17] D.G. Luenberger, An introduction to observers, IEEE Transactions on Automatic Control 16 (1971) [18] J. Park, G. Rizzoni, An eigenstructure assignment algorithm for the design of fault detection filters, IEEE Transactions on Automatic Control 39 (1994) [19] C.C. Tsui, New approach to robust observer design, International Journal of Control 47 (1988) [20] C.C. Tsui, A complete analytical solution to the equation TA FT ¼ LC and its applications, IEEE Transactions on Automatic Control 32 (1987) [21] C.C. Tsui, On the solution to matrix equation TA FT ¼ LC and its applications, SIAM Journal of Matrix Analysis and Applications 14 (1) (1993) [22] Q.-W. Wang, H.-X. Chang, Q. Ning, The common solution to six quaternion matrix equations with applications, Applied Mathematics and Computation 198 (1) (2008) [23] F. Rincon, Feedback stabilization of second-order models, Ph.D. dissertation, Northern Illinois University, De Kalb, Illinois, USA, [24] Y. Kim, H.S. Kim, Eigenstructure assignment algorithm for mechanical second-order systems, Journal of Guidance, Control and Dynamics 22 (5) (1999) [25] E.K. Chu, B.N. Datta, Numerically robust pole assignment for second-order systems, International Journal of Control 64 (4) (1996) [26] D.J. Inman, A. Kress, Eigenstructure assignment algorithm for second-order systems, Journal of Guidance, Control and Dynamics 22 (5) (1999) [27] F. Ding, P.X. Liu, J. Ding, Iterative solutions of the generalized Sylvester matrix equations by using the hierarchical identification principle, Applied Mathematics and Computation 197 (1) (2008)
10 336 B. Zhou et al. / Applied Mathematics and Computation 212 (2009) [28] M.A. Ramadan, E.A. El-Sayed, On the matrix equation XH ¼ HX and the associated controllability problem, Applied Mathematics and Computation 186 (1) (2007) [29] A.G. Wu, G.R. Duan, Kronecker maps and Sylvester-polynomial matrix equations, IEEE Transactions on Automatic Control 52 (5) (2007) [30] L. Huang, The explicit solutions and solvability of linear matrix equations, Linear Algebra and its Applications 311 (2000)
Solutions to generalized Sylvester matrix equation by Schur decomposition
International Journal of Systems Science Vol 8, No, May 007, 9 7 Solutions to generalized Sylvester matrix equation by Schur decomposition BIN ZHOU* and GUANG-REN DUAN Center for Control Systems and Guidance
More informationSolutions to the generalized Sylvester matrix equations by a singular value decomposition
Journal of Control Theory Applications 2007 5 (4) 397 403 DOI 101007/s11768-006-6113-0 Solutions to the generalized Sylvester matrix equations by a singular value decomposition Bin ZHOU Guangren DUAN (Center
More informationClosed-form Solutions to the Matrix Equation AX EXF = BY with F in Companion Form
International Journal of Automation and Computing 62), May 2009, 204-209 DOI: 101007/s11633-009-0204-6 Closed-form Solutions to the Matrix Equation AX EX BY with in Companion orm Bin Zhou Guang-Ren Duan
More information1030 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 56, NO. 5, MAY 2011
1030 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 56, NO 5, MAY 2011 L L 2 Low-Gain Feedback: Their Properties, Characterizations Applications in Constrained Control Bin Zhou, Member, IEEE, Zongli Lin,
More informationPARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT
PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse
More informationRational stabilising commutative controllers: parameterisation and characterisation of degrees of freedom
International Journal of Control Vol 79, No 12, December 26, 161 1612 Rational stabilising commutative controllers: parameterisation and characterisation of degrees of freedom M V MOREIRA*y, J C BASILIOy
More informationResearch Article Finite Iterative Algorithm for Solving a Complex of Conjugate and Transpose Matrix Equation
indawi Publishing Corporation Discrete Mathematics Volume 013, Article ID 17063, 13 pages http://dx.doi.org/10.1155/013/17063 Research Article Finite Iterative Algorithm for Solving a Complex of Conjugate
More informationOptimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications
Optimization problems on the rank and inertia of the Hermitian matrix expression A BX (BX) with applications Yongge Tian China Economics and Management Academy, Central University of Finance and Economics,
More informationAN ITERATIVE METHOD TO SOLVE SYMMETRIC POSITIVE DEFINITE MATRIX EQUATIONS
AN ITERATIVE METHOD TO SOLVE SYMMETRIC POSITIVE DEFINITE MATRIX EQUATIONS DAVOD KHOJASTEH SALKUYEH and FATEMEH PANJEH ALI BEIK Communicated by the former editorial board Let A : R m n R m n be a symmetric
More informationLetting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc.
1 Polynomial Matrices 1.1 Polynomials Letting be a field, e.g., of the real numbers, the complex numbers, the rational numbers, the rational functions W(s) of a complex variable s, etc., n ws ( ) as a
More informationSolvability of Linear Matrix Equations in a Symmetric Matrix Variable
Solvability of Linear Matrix Equations in a Symmetric Matrix Variable Maurcio C. de Oliveira J. William Helton Abstract We study the solvability of generalized linear matrix equations of the Lyapunov type
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationIMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS. Shumei Mu, Tianguang Chu, and Long Wang
IMPULSIVE CONTROL OF DISCRETE-TIME NETWORKED SYSTEMS WITH COMMUNICATION DELAYS Shumei Mu Tianguang Chu and Long Wang Intelligent Control Laboratory Center for Systems and Control Department of Mechanics
More informationMinimalsinglelinearfunctionalobserversforlinearsystems
Minimalsinglelinearfunctionalobserversforlinearsystems Frédéric Rotella a Irène Zambettakis b a Ecole Nationale d Ingénieurs de Tarbes Laboratoire de Génie de production 47 avenue d Azereix 65016 Tarbes
More informationAn Even Order Symmetric B Tensor is Positive Definite
An Even Order Symmetric B Tensor is Positive Definite Liqun Qi, Yisheng Song arxiv:1404.0452v4 [math.sp] 14 May 2014 October 17, 2018 Abstract It is easily checkable if a given tensor is a B tensor, or
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (2011) 797 802 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: wwwelseviercom/locate/aml Model order determination using the Hankel
More informationStability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers
International Journal of Approximate Reasoning 46 (2007) 109 119 www.elsevier.com/locate/ijar Stability analysis for a class of Takagi Sugeno fuzzy control systems with PID controllers Le Hung Lan * Institute
More informationON THE GLOBAL KRYLOV SUBSPACE METHODS FOR SOLVING GENERAL COUPLED MATRIX EQUATIONS
ON THE GLOBAL KRYLOV SUBSPACE METHODS FOR SOLVING GENERAL COUPLED MATRIX EQUATIONS Fatemeh Panjeh Ali Beik and Davod Khojasteh Salkuyeh, Department of Mathematics, Vali-e-Asr University of Rafsanjan, Rafsanjan,
More informationSynopsis of Numerical Linear Algebra
Synopsis of Numerical Linear Algebra Eric de Sturler Department of Mathematics, Virginia Tech sturler@vt.edu http://www.math.vt.edu/people/sturler Iterative Methods for Linear Systems: Basics to Research
More informationFormulas for the Drazin Inverse of Matrices over Skew Fields
Filomat 30:12 2016 3377 3388 DOI 102298/FIL1612377S Published by Faculty of Sciences and Mathematics University of Niš Serbia Available at: http://wwwpmfniacrs/filomat Formulas for the Drazin Inverse of
More informationMath 489AB Exercises for Chapter 2 Fall Section 2.3
Math 489AB Exercises for Chapter 2 Fall 2008 Section 2.3 2.3.3. Let A M n (R). Then the eigenvalues of A are the roots of the characteristic polynomial p A (t). Since A is real, p A (t) is a polynomial
More informationAn iterative method for the skew-symmetric solution and the optimal approximate solution of the matrix equation AXB =C
Journal of Computational and Applied Mathematics 1 008) 31 44 www.elsevier.com/locate/cam An iterative method for the skew-symmetric solution and the optimal approximate solution of the matrix equation
More informationSums of diagonalizable matrices
Linear Algebra and its Applications 315 (2000) 1 23 www.elsevier.com/locate/laa Sums of diagonalizable matrices J.D. Botha Department of Mathematics University of South Africa P.O. Box 392 Pretoria 0003
More informationH 2 -optimal model reduction of MIMO systems
H 2 -optimal model reduction of MIMO systems P. Van Dooren K. A. Gallivan P.-A. Absil Abstract We consider the problem of approximating a p m rational transfer function Hs of high degree by another p m
More informationA New Block Algorithm for Full-Rank Solution of the Sylvester-observer Equation.
1 A New Block Algorithm for Full-Rank Solution of the Sylvester-observer Equation João Carvalho, DMPA, Universidade Federal do RS, Brasil Karabi Datta, Dep MSc, Northern Illinois University, DeKalb, IL
More informationComputational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1
Computational Methods for Feedback Control in Damped Gyroscopic Second-order Systems 1 B. N. Datta, IEEE Fellow 2 D. R. Sarkissian 3 Abstract Two new computationally viable algorithms are proposed for
More information5 Linear Transformations
Lecture 13 5 Linear Transformations 5.1 Basic Definitions and Examples We have already come across with the notion of linear transformations on euclidean spaces. We shall now see that this notion readily
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More information16. Local theory of regular singular points and applications
16. Local theory of regular singular points and applications 265 16. Local theory of regular singular points and applications In this section we consider linear systems defined by the germs of meromorphic
More informationInfinite elementary divisor structure-preserving transformations for polynomial matrices
Infinite elementary divisor structure-preserving transformations for polynomial matrices N P Karampetakis and S Vologiannidis Aristotle University of Thessaloniki, Department of Mathematics, Thessaloniki
More informationAbsolute value equations
Linear Algebra and its Applications 419 (2006) 359 367 www.elsevier.com/locate/laa Absolute value equations O.L. Mangasarian, R.R. Meyer Computer Sciences Department, University of Wisconsin, 1210 West
More informationON THE MATRIX EQUATION XA AX = X P
ON THE MATRIX EQUATION XA AX = X P DIETRICH BURDE Abstract We study the matrix equation XA AX = X p in M n (K) for 1 < p < n It is shown that every matrix solution X is nilpotent and that the generalized
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationGroup inverse for the block matrix with two identical subblocks over skew fields
Electronic Journal of Linear Algebra Volume 21 Volume 21 2010 Article 7 2010 Group inverse for the block matrix with two identical subblocks over skew fields Jiemei Zhao Changjiang Bu Follow this and additional
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationMa 227 Review for Systems of DEs
Ma 7 Review for Systems of DEs Matrices Basic Properties Addition and subtraction: Let A a ij mn and B b ij mn.then A B a ij b ij mn 3 A 6 B 6 4 7 6 A B 6 4 3 7 6 6 7 3 Scaler Multiplication: Let k be
More informationA Control Methodology for Constrained Linear Systems Based on Positive Invariance of Polyhedra
A Control Methodology for Constrained Linear Systems Based on Positive Invariance of Polyhedra Jean-Claude HENNET LAAS-CNRS Toulouse, France Co-workers: Marina VASSILAKI University of Patras, GREECE Jean-Paul
More information3 (Maths) Linear Algebra
3 (Maths) Linear Algebra References: Simon and Blume, chapters 6 to 11, 16 and 23; Pemberton and Rau, chapters 11 to 13 and 25; Sundaram, sections 1.3 and 1.5. The methods and concepts of linear algebra
More information2 Two-Point Boundary Value Problems
2 Two-Point Boundary Value Problems Another fundamental equation, in addition to the heat eq. and the wave eq., is Poisson s equation: n j=1 2 u x 2 j The unknown is the function u = u(x 1, x 2,..., x
More informationELEMENTARY LINEAR ALGEBRA
ELEMENTARY LINEAR ALGEBRA K R MATTHEWS DEPARTMENT OF MATHEMATICS UNIVERSITY OF QUEENSLAND First Printing, 99 Chapter LINEAR EQUATIONS Introduction to linear equations A linear equation in n unknowns x,
More informationThe Drazin inverses of products and differences of orthogonal projections
J Math Anal Appl 335 7 64 71 wwwelseviercom/locate/jmaa The Drazin inverses of products and differences of orthogonal projections Chun Yuan Deng School of Mathematics Science, South China Normal University,
More informationarxiv: v1 [math.ra] 11 Aug 2014
Double B-tensors and quasi-double B-tensors Chaoqian Li, Yaotang Li arxiv:1408.2299v1 [math.ra] 11 Aug 2014 a School of Mathematics and Statistics, Yunnan University, Kunming, Yunnan, P. R. China 650091
More informationChapter 1: Systems of Linear Equations
Chapter : Systems of Linear Equations February, 9 Systems of linear equations Linear systems Lecture A linear equation in variables x, x,, x n is an equation of the form a x + a x + + a n x n = b, where
More informationarxiv: v1 [math.na] 13 Mar 2019
Relation between the T-congruence Sylvester equation and the generalized Sylvester equation Yuki Satake, Masaya Oozawa, Tomohiro Sogabe Yuto Miyatake, Tomoya Kemmochi, Shao-Liang Zhang arxiv:1903.05360v1
More informationVandermonde-form Preserving Matrices And The Generalized Signal Richness Preservation Problem
Vandermonde-form Preserving Matrices And The Generalized Signal Richness Preservation Problem Borching Su Department of Electrical Engineering California Institute of Technology Pasadena, California 91125
More informationA connection between number theory and linear algebra
A connection between number theory and linear algebra Mark Steinberger Contents 1. Some basics 1 2. Rational canonical form 2 3. Prime factorization in F[x] 4 4. Units and order 5 5. Finite fields 7 6.
More informationPossible numbers of ones in 0 1 matrices with a given rank
Linear and Multilinear Algebra, Vol, No, 00, Possible numbers of ones in 0 1 matrices with a given rank QI HU, YAQIN LI and XINGZHI ZHAN* Department of Mathematics, East China Normal University, Shanghai
More informationMIT Algebraic techniques and semidefinite optimization February 16, Lecture 4
MIT 6.972 Algebraic techniques and semidefinite optimization February 16, 2006 Lecture 4 Lecturer: Pablo A. Parrilo Scribe: Pablo A. Parrilo In this lecture we will review some basic elements of abstract
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 433 (2010) 476 482 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Nonsingularity of the
More informationPartial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution
Partial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution Biswa N. Datta, IEEE Fellow Department of Mathematics Northern Illinois University DeKalb, IL, 60115 USA e-mail:
More informationACI-matrices all of whose completions have the same rank
ACI-matrices all of whose completions have the same rank Zejun Huang, Xingzhi Zhan Department of Mathematics East China Normal University Shanghai 200241, China Abstract We characterize the ACI-matrices
More informationA Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case
A Simple Derivation of Right Interactor for Tall Transfer Function Matrices and its Application to Inner-Outer Factorization Continuous-Time Case ATARU KASE Osaka Institute of Technology Department of
More information1 Last time: least-squares problems
MATH Linear algebra (Fall 07) Lecture Last time: least-squares problems Definition. If A is an m n matrix and b R m, then a least-squares solution to the linear system Ax = b is a vector x R n such that
More informationLinear Algebra. Session 8
Linear Algebra. Session 8 Dr. Marco A Roque Sol 08/01/2017 Abstract Linear Algebra Range and kernel Let V, W be vector spaces and L : V W, be a linear mapping. Definition. The range (or image of L is the
More informationMatrices and Vectors. Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A =
30 MATHEMATICS REVIEW G A.1.1 Matrices and Vectors Definition of Matrix. An MxN matrix A is a two-dimensional array of numbers A = a 11 a 12... a 1N a 21 a 22... a 2N...... a M1 a M2... a MN A matrix can
More informationONGOING WORK ON FAULT DETECTION AND ISOLATION FOR FLIGHT CONTROL APPLICATIONS
ONGOING WORK ON FAULT DETECTION AND ISOLATION FOR FLIGHT CONTROL APPLICATIONS Jason M. Upchurch Old Dominion University Systems Research Laboratory M.S. Thesis Advisor: Dr. Oscar González Abstract Modern
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 8 MATRICES III Rank of a matrix 2 General systems of linear equations 3 Eigenvalues and eigenvectors Rank of a matrix
More informationIntroduction to Matrices
214 Analysis and Design of Feedback Control Systems Introduction to Matrices Derek Rowell October 2002 Modern system dynamics is based upon a matrix representation of the dynamic equations governing the
More informationDelay-Dependent Stability Criteria for Linear Systems with Multiple Time Delays
Delay-Dependent Stability Criteria for Linear Systems with Multiple Time Delays Yong He, Min Wu, Jin-Hua She Abstract This paper deals with the problem of the delay-dependent stability of linear systems
More informationLinear Algebra. Linear Equations and Matrices. Copyright 2005, W.R. Winfrey
Copyright 2005, W.R. Winfrey Topics Preliminaries Systems of Linear Equations Matrices Algebraic Properties of Matrix Operations Special Types of Matrices and Partitioned Matrices Matrix Transformations
More informationM.6. Rational canonical form
book 2005/3/26 16:06 page 383 #397 M.6. RATIONAL CANONICAL FORM 383 M.6. Rational canonical form In this section we apply the theory of finitely generated modules of a principal ideal domain to study the
More informationResearch Article Constrained Solutions of a System of Matrix Equations
Journal of Applied Mathematics Volume 2012, Article ID 471573, 19 pages doi:10.1155/2012/471573 Research Article Constrained Solutions of a System of Matrix Equations Qing-Wen Wang 1 and Juan Yu 1, 2 1
More informationApplied Mathematics Letters. Complements to full order observer design for linear systems with unknown inputs
pplied Mathematics Letters 22 (29) 117 1111 ontents lists available at ScienceDirect pplied Mathematics Letters journal homepage: www.elsevier.com/locate/aml omplements to full order observer design for
More informationSimultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer
Preprints of the 19th World Congress The International Federation of Automatic Control Simultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer Fengming Shi*, Ron J.
More informationSolution to Sylvester equation associated to linear descriptor systems
Solution to Sylvester equation associated to linear descriptor systems Mohamed Darouach To cite this version: Mohamed Darouach. Solution to Sylvester equation associated to linear descriptor systems. Systems
More informationConsensus Stabilizability and Exact Consensus Controllability of Multi-agent Linear Systems
Consensus Stabilizability and Exact Consensus Controllability of Multi-agent Linear Systems M. ISABEL GARCÍA-PLANAS Universitat Politècnica de Catalunya Departament de Matèmatiques Minería 1, Esc. C, 1-3,
More informationINTRODUCTION TO LIE ALGEBRAS. LECTURE 10.
INTRODUCTION TO LIE ALGEBRAS. LECTURE 10. 10. Jordan decomposition: theme with variations 10.1. Recall that f End(V ) is semisimple if f is diagonalizable (over the algebraic closure of the base field).
More informationApproximate Low Rank Solution of Generalized Lyapunov Matrix Equations via Proper Orthogonal Decomposition
Applied Mathematical Sciences, Vol. 4, 2010, no. 1, 21-30 Approximate Low Rank Solution of Generalized Lyapunov Matrix Equations via Proper Orthogonal Decomposition Amer Kaabi Department of Basic Science
More informationEIGENVALUES AND EIGENVECTORS 3
EIGENVALUES AND EIGENVECTORS 3 1. Motivation 1.1. Diagonal matrices. Perhaps the simplest type of linear transformations are those whose matrix is diagonal (in some basis). Consider for example the matrices
More informationA Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay
A Multi-Step Hybrid Method for Multi-Input Partial Quadratic Eigenvalue Assignment with Time Delay Zheng-Jian Bai Mei-Xiang Chen Jin-Ku Yang April 14, 2012 Abstract A hybrid method was given by Ram, Mottershead,
More informationRank-one LMIs and Lyapunov's Inequality. Gjerrit Meinsma 4. Abstract. We describe a new proof of the well-known Lyapunov's matrix inequality about
Rank-one LMIs and Lyapunov's Inequality Didier Henrion 1;; Gjerrit Meinsma Abstract We describe a new proof of the well-known Lyapunov's matrix inequality about the location of the eigenvalues of a matrix
More informationResearch Article Indefinite LQ Control for Discrete-Time Stochastic Systems via Semidefinite Programming
Mathematical Problems in Engineering Volume 2012, Article ID 674087, 14 pages doi:10.1155/2012/674087 Research Article Indefinite LQ Control for Discrete-Time Stochastic Systems via Semidefinite Programming
More informationNorm invariant discretization for sampled-data fault detection
Automatica 41 (25 1633 1637 www.elsevier.com/locate/automatica Technical communique Norm invariant discretization for sampled-data fault detection Iman Izadi, Tongwen Chen, Qing Zhao Department of Electrical
More informationApplied Mathematics Letters
Applied Mathematics Letters (009) 15 130 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Spectral characterizations of sandglass graphs
More informationChapter 4. Matrices and Matrix Rings
Chapter 4 Matrices and Matrix Rings We first consider matrices in full generality, i.e., over an arbitrary ring R. However, after the first few pages, it will be assumed that R is commutative. The topics,
More informationInverses. Stephen Boyd. EE103 Stanford University. October 28, 2017
Inverses Stephen Boyd EE103 Stanford University October 28, 2017 Outline Left and right inverses Inverse Solving linear equations Examples Pseudo-inverse Left and right inverses 2 Left inverses a number
More informationLinear Algebra. Week 7
Linear Algebra. Week 7 Dr. Marco A Roque Sol 10 / 09 / 2018 If {v 1, v 2,, v n } is a basis for a vector space V, then any vector v V, has a unique representation v = x 1 v 1 + x 2 v 2 + + x n v n where
More informationLinear Algebra Primer
Linear Algebra Primer David Doria daviddoria@gmail.com Wednesday 3 rd December, 2008 Contents Why is it called Linear Algebra? 4 2 What is a Matrix? 4 2. Input and Output.....................................
More informationOn the application of different numerical methods to obtain null-spaces of polynomial matrices. Part 1: block Toeplitz algorithms.
On the application of different numerical methods to obtain null-spaces of polynomial matrices. Part 1: block Toeplitz algorithms. J.C. Zúñiga and D. Henrion Abstract Four different algorithms are designed
More informationMatrix expression and reachability analysis of finite automata
J Control Theory Appl 2012 10 (2) 210 21 DOI 101007/s11768-012-1178-4 Matrix expression and reachability analysis of finite automata Xiangru XU, Yiguang HONG Key Laboratory of Systems and Control, Academy
More informationA matrix over a field F is a rectangular array of elements from F. The symbol
Chapter MATRICES Matrix arithmetic A matrix over a field F is a rectangular array of elements from F The symbol M m n (F ) denotes the collection of all m n matrices over F Matrices will usually be denoted
More informationH. Schubert, Kalkül der abzählenden Geometrie, 1879.
H. Schubert, Kalkül der abzählenden Geometrie, 879. Problem: Given mp subspaces of dimension p in general position in a vector space of dimension m + p, how many subspaces of dimension m intersect all
More informationEXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)
EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily
More informationSTABILIZATION OF LINEAR SYSTEMS VIA DELAYED STATE FEEDBACK CONTROLLER. El-Kébir Boukas. N. K. M Sirdi. Received December 2007; accepted February 2008
ICIC Express Letters ICIC International c 28 ISSN 1881-83X Volume 2, Number 1, March 28 pp. 1 6 STABILIZATION OF LINEAR SYSTEMS VIA DELAYED STATE FEEDBACK CONTROLLER El-Kébir Boukas Department of Mechanical
More informationChapter 5. Linear Algebra. A linear (algebraic) equation in. unknowns, x 1, x 2,..., x n, is. an equation of the form
Chapter 5. Linear Algebra A linear (algebraic) equation in n unknowns, x 1, x 2,..., x n, is an equation of the form a 1 x 1 + a 2 x 2 + + a n x n = b where a 1, a 2,..., a n and b are real numbers. 1
More informationBare minimum on matrix algebra. Psychology 588: Covariance structure and factor models
Bare minimum on matrix algebra Psychology 588: Covariance structure and factor models Matrix multiplication 2 Consider three notations for linear combinations y11 y1 m x11 x 1p b11 b 1m y y x x b b n1
More informationMEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction
MEM 355 Performance Enhancement of Dynamical Systems MIMO Introduction Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 11/2/214 Outline Solving State Equations Variation
More informationAn algorithm for symmetric generalized inverse eigenvalue problems
Linear Algebra and its Applications 296 (999) 79±98 www.elsevier.com/locate/laa An algorithm for symmetric generalized inverse eigenvalue problems Hua Dai Department of Mathematics, Nanjing University
More informationJournal of Computational Physics
Journal of Computational Physics 9 () 759 763 Contents lists available at ScienceDirect Journal of Computational Physics journal homepage: www.elsevier.com/locate/jcp Short Note A comment on the computation
More informationAnalytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix
Analytical formulas for calculating the extremal ranks and inertias of A + BXB when X is a fixed-rank Hermitian matrix Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081, China
More informationApplied Mathematics Letters
Applied Mathematics Letters 5 (0) 554 559 Contents lists available at SciVerse ScienceDirect Applied Mathematics Letters journal homepage: wwwelseviercom/locate/aml On the (s, t)-pell and (s, t)-pell Lucas
More informationELA
POSITIVE DEFINITE SOLUTION OF THE MATRIX EQUATION X = Q+A H (I X C) δ A GUOZHU YAO, ANPING LIAO, AND XUEFENG DUAN Abstract. We consider the nonlinear matrix equation X = Q+A H (I X C) δ A (0 < δ 1), where
More information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationChapter 4 - MATRIX ALGEBRA. ... a 2j... a 2n. a i1 a i2... a ij... a in
Chapter 4 - MATRIX ALGEBRA 4.1. Matrix Operations A a 11 a 12... a 1j... a 1n a 21. a 22.... a 2j... a 2n. a i1 a i2... a ij... a in... a m1 a m2... a mj... a mn The entry in the ith row and the jth column
More informationPOLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK
1 POLYNOMIAL EMBEDDING ALGORITHMS FOR CONTROLLERS IN A BEHAVIORAL FRAMEWORK H.L. Trentelman, Member, IEEE, R. Zavala Yoe, and C. Praagman Abstract In this paper we will establish polynomial algorithms
More informationDeterminants of Partition Matrices
journal of number theory 56, 283297 (1996) article no. 0018 Determinants of Partition Matrices Georg Martin Reinhart Wellesley College Communicated by A. Hildebrand Received February 14, 1994; revised
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 430 (2009) 579 586 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Low rank perturbation
More informationResearch Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
Applied Mathematics Volume 202, Article ID 689820, 3 pages doi:0.55/202/689820 Research Article Convex Polyhedron Method to Stability of Continuous Systems with Two Additive Time-Varying Delay Components
More informationA q x k+q + A q 1 x k+q A 0 x k = 0 (1.1) where k = 0, 1, 2,..., N q, or equivalently. A(σ)x k = 0, k = 0, 1, 2,..., N q (1.
A SPECTRAL CHARACTERIZATION OF THE BEHAVIOR OF DISCRETE TIME AR-REPRESENTATIONS OVER A FINITE TIME INTERVAL E.N.Antoniou, A.I.G.Vardulakis, N.P.Karampetakis Aristotle University of Thessaloniki Faculty
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More information