A Note on Solutions of the Matrix Equation AXB = C
|
|
- Clinton Wilson
- 5 years ago
- Views:
Transcription
1 SCIENTIFIC PUBLICATIONS OF THE STATE UNIVERSITY OF NOVI PAZAR SER A: APPL MATH INFORM AND MECH vol 6, 1 (214), A Note on Solutions of the Matrix Equation AXB C I V Jovović, B J Malešević Abstract: This paper deals with necessary and sufficient condition for consistency of the matrix equation AXB C We will be concerned with the minimal number of free parameters in Penrose s formula X A (1) CB (1) +Y A (1) AY BB (1) for obtaining the general solution of the matrix equation and we will establish the relation between the minimal number of free parameters and the ranks of the matrices A and B The solution is described in the terms of Rohde s general form of the {1}-inverse of the matrices A and B We will also use Kronecker product to transform the matrix equation AXB C into the linear system (B T A) X C Keywords: Generalized inverses, Kronecker product, matrix equations, linear systems 1 Introduction In this paper we consider matrix equation AXB C, (1) where X is an n k matrix of unknowns, A is an m n matrix of rank a, B is a k l matrix of rank b, and C is an m l matrix, all over C The set of all m n matrices over the complex field C will be denoted by C m n, m,n N The set of all m n matrices over the complex field C of rank a will be denoted by C m n a We will write A i (A j ) for the i th row (the j th column) of the matrix A C m n and A will denote an ordered stock of columns of A, ie A A 1 A 2 Using the Kronecker product of the matrices B T and A we can transform the matrix equation (1) into linear system A n (B T A) X C (2) Manuscript received September 3, 213 ; accepted November 17, 213 I V Jovović and B J Malešević are with School of Electrical Engineering, University of Belgrade, Serbia 45
2 46 I V Jovović, B J Malešević For the proof we refer the reader to A Graham 2 Necessary and sufficient condition for consistency of the linear system Ax c, as well as the minimal number of free parameters in Penrose s formula x A (1) c + (I AA (1) )y has been considered in the paper B Malešević, I Jovović, M Makragić and B Radičić 3 We will here briefly sketch this results in the case of the linear system (2) Any matrix X satisfying the equality AXA A is called {1}-inverse of A and is denoted by A (1) For each matrix A C m n a there are regular matrices P C n n and Q C m m such that Ia QAP E A, (3) where I a is a a identity matrix It can be easily seen that every {1}-inverse of the matrix A can be represented in the Rohde s form A (1) Ia U P V W Q, (4) where U u i j, V v i j and W w i j are arbitrary matrices of corresponding dimensions a (m a), (n a) a and (n a) (m a) with mutually independent entries, see C Rohde 1 and V Perić 8 We will explore the minimal numbers of free parameters in Penrose s formula X A (1) CB (1) +Y A (1) AY BB (1) for obtaining the general solution of the matrix equation (1) Some similar considerations can be found in papers B Malešević and B Radičić 4, 5, 6 and 9 2 Matrix equation AXBC and the Kronecker product of the matrices B T and A The Kronecker product of matrices A a i j m n C m n and B b i j k l C k l, denoted by A B, is defined as block matrix a 11 B a 12 B a 1n B a 21 B a 22 B a 2n B A B a m1 B a m2 B a mn B The matrix A B is mk nl matrix with mn blocks a i j B of order k l Here we will mention some properties and rules for the Kronecker product Let A C m n, B C k l, C C n r and D C l s Then the following propositions holds: A T B T (A B) T ; rank(a B) rank(a) rank(b);
3 A note on solutions of the matrix equation AXB C 47 (A B) (C D) (A C) (B D) (mixed product rule); if A and B are regular n n and k k matrices, then (A B) 1 A 1 B 1 The proof of these facts can be found in A Graham 2 and A Ben Israel and TNE Greville 1 Matrix A (1) B (1) is {1}-inverse of A B Using mixed product rule we have (A B)(A (1) B (1) )(A B) (A A (1) A) (B B (1) B) A B Let R C k k and S C l l be regular matrices such that Ib RBS E B (5) An {1}-inverse of the matrix B can be represented in the Rohde s form B (1) Ib M S R, (6) N K where M m i j, N n i j and K k i j are arbitrary matrices of corresponding dimensions b (k b), (l b) b and (l b) (k b) with mutually independent entries From now on, we will look more closely at the linear system (2) Firstly, by mixed product rule we obtain (S T Q) (B T A) (R T P) (S T B T R T ) (Q A P) E B T E A Unfortunately, the matrix Ib E E B T E A A E A E A E A I a I a I a is not of the needed form E B T A Equality E B T E A E B T A holds for b 1 Swapping the rows and the columns corresponding to blocks I a and to zero diagonal blocks we get required matrix E B T A If matrices D and G are the elementary matrices obtained by swapping rows and columns corresponding to mentioned blocks of the identity matrices, then D (E B T E A ) G E B T A Thus, we have (D (S T Q)) (B T A) ((R T P) G) E B T A (7) and so an {1}-inverse of the matrix B T A can be represented in the Rohde s form (B T A) (1) (R T Iab F P) G D (S T Q), (8) H L
4 48 I V Jovović, B J Malešević where F f i j, H h i j and L l i j are arbitrary matrices of corresponding dimensions ab (ml ab), (nk ab) ab and (nk ab) (ml ab) with mutually independent entries If the matrices A and B are square matrices, then D G T For the simplicity of notation, we will write c a (c a ) for the submatrix corresponding to the first (the last) a coordinates of the vector c Now, we can rephrase Lemma 21 and Theorem 22 from the paper B Malešević, I Jovović, M Makragić and B Radičić 3 in the case of the linear system (2) Let C be given by C QCS Then C (S T Q) C Lemma 21 The linear system (2) has a solution if and only if the last ml ab coordinates of the vector c D C are zeros, where D is elementary matrix such that (7) holds Theorem 22 The vector X (B T A) (1) C + (I (B T A) (1) (B T A))y, (9) y C nk 1 is an arbitrary column, is the general solution of the system (2), if and only if the {1}-inverse (B T A) (1) of the system matrix B T A has the form (8) for arbitrary matrices F and L and the rows of the matrix H(c ab y ab ) + y nk ab are free parameters, c where D (S T Q) C c ab and G 1 ((R 1 ) T P 1 )y y y ab y nk ab In the paper B Malešević, I Jovović, M Makragić and B Radičić 3 we have seen that general solution (9) can be presented in the form X (R T c P) G ab H(c ab y ab ) + y (1) nk ab We illustrate this formula in the next example Example 23 We consider the matrix equation where A , X x11 x 12 x 13 x 21 x 22 x 23 AXB C,, B and C Using the Kronecker product the matrix equation may be considered in the form of the equivalent linear system (B T A) X C, ie x 11 x 21 x 12 x 22 x 13 x
5 1 Matrices Q 1 and P E A 1 and matrices R matrices such that RBS E B A note on solutions of the matrix equation AXB C are regular matrices such that QAP and S Therefore, are regular E B T E A and for matrix D we have E B T A D (E B T E A ) Let us remark that E B T A (D (S T Q)) (B T A) (R T P), and hence according to the Theorem 22 X (R T c 4 P) H(c 4 y 4 ) + y is the general solution of the linear 2 system iff elements of the column H(c 4 y 4 )+y 2 are two mutually independent parameters
6 5 I V Jovović, B J Malešević α 1 i α 2 for the vector c c 4 (D (S T Q)) C Finally, the general solution of the linear system is X (R T P) 3 1 α 1 α α 1 α 2 1 α 1 + 2α 2 1 α 2 α 1 + 2α 2 α 2 2α 3 α 4 α 2 α 1 + 2α 2 α 2 3 Matrix equation AXBC and the {1} inverses of the matrices A and B In this section we indicate how technique of an {1}-inverse may be used to obtain the necessary and sufficient condition for an existence of a general solution of the matrix equation (1) without using Kronecker product We will use the symbols C a,b, C a,b, C a,b and C a,b for the submatrices of the matrix C corresponding to the first a rows and b columns, the last a rows and the first b columns, the first a rows and the last b columns, the last a rows and b columns, respectively Lemma 31 The matrix equation (1) has a solution if and only if the last m a rows and l b columns of the matrix C QCS are zeros, where Q C m m and S C l l are regular matrices such that (3) and (5) hold Proof: The matrix equation (1) has a solution if and only if C AA (1) CB (1) B, see R Penrose 7 Since A (1) and B (1) are described by the equations (4) and (6), it follows that AA (1) Ia U AP V W Q Q 1 Ia U Q
7 A note on solutions of the matrix equation AXB C 51 and B (1) Ib M B S N K Ib RB S N S 1 Hence, since Q and S are regular matrices we have the following equivalences C AA (1) CB (1) B QCSQAA (1) CB (1) BS C QCS C Ia U C a,b C a,l b C m a,b C m a,l b C a,b C a,l b C m a,b C m a,l b for C C a,b C a,l b C m a,b C m a,l b Ia U C a,b C a,l b C m a,b C m a,l b C Ib N Ib N C a,b +U C m a,b +C a,l b N +U C n a,l b N Therefore, we conclude, C AA (1) CB (1) B C a,l b C m a,b C m a,l b As we have seen in the Lemma 21 the matrix equation (1) has a solution if and only if the last ml ab coordinates of the column c D C are zeros, where D is elementary matrix such that (7) holds Here we obtain the same result without using Kronecker product The last m(l b) elements of the column C are zeros and there are b blocks of m a zeros Multiplying by the left column C with elementary matrix D switches the rows corresponding to this zeros blocks under the blocks C a,b i, 1 i b Hence, the last m(l b) + (m a)b ml ab entries of the column c are zeros Furthermore, we give a new form of the general solution of the matrix equation (1) using {1}-inverses of the matrices A and B Theorem 32 The matrix X A (1) CB (1) +Y A (1) AY BB (1), (11) Y C n k is an arbitrary matrix, is the general solution of the matrix equation (1) if and only if the {1}-inverses A (1) and B (1) of the matrices A and B have the forms (4) and (6) for arbitrary matrices U, W, N and K and the entries of the matrices V ( C Y ) + Y, ( C Y ) M + Y, V ( C Y ) M + Y a,b a,b n a,b a,b a,b a,k b a,b a,b n a,k b (12) are mutually independent free parameters, where QCS C C a,b and P 1 Y R 1 Y Y a,b Y a,k b Y n a,b Y n a,k b (13)
8 52 I V Jovović, B J Malešević Proof: Since the {1}-inverses A (1) and B (1) of the matrices A and B have the forms (4) and (6), the solution of the system X A (1) CB (1) +Y A (1) AY BB (1) can be represented in the form Ia U X P V W Ib M QCS N K Ia U R + Y P V W QAPP 1 Y R 1 Ib M RBS N K R According to Lemma 31 and from (3) and (5) we have Ia U C X P a,b Ib M R + Y V W N K Ia U Ia P P 1 Y R 1 Ib Ib M V W N K R Furthermore, we obtain C a,b C a,b M I a X P V C a,b V C R + Y P P 1 Y R 1 I b M R a,b M V ( ) C a,b C a,b M P V C a,b V C + Y I a Y I b M R a,b M V ( C a,b C a,b M P V C a,b V C + a,b M ) Y a,b Y a,k b I a Y a,b Y a,k b I b M + Y n a,b Y n a,k b V Y n a,b Y R n a,k b ( ) C a,b C a,b M Y a,b Y a,k b Y a,b Y a,b M P V C a,b V C + a,b M V Y a,b V Y R, a,b M Y n a,b Y n a,k b where Y P 1 Y R 1 We now conclude X P C a,b V ( C a,b Y a,b) + Y n a,b ( C a,b Y a,b)m +Y a,k b V ( C a,b Y a,b)m +Y n a,k b R (14) According to the Theorem 22 the general solution of the equation (1) has nk ab free parameters Therefore, since the matrices P and R are regular we deduce that the solution (14) is the general if and only if the entries of the matrices separately V ( C a,b Y a,b ) + Y n a,b, ( C a,b Y a,b ) M +Y a,k b, V ( C a,b Y a,b ) M +Y n a,k b are nk ab free parameters We can illustrate the Theorem 32 on the following two examples
9 A note on solutions of the matrix equation AXB C 53 Example 33 Consider again the matrix equation from previous example A B C, where A B x11 x, X 12 x 13 x 21 x 22 x 23 and C, We have C QCS and C 2,2 1 The solution of the matrix equation AXB C is 3 α1 X P 1 α 2 R 1 1 Example 34 We now consider the matrix equation α1 1 α α1 + 2α 2 α 1 + 2α 2 α 1 2α 2 AXB C, x 11 x 12 where A 2 6 4, X x 21 x 22, B x 31 x 32 1 For regular matrices Q 1 and P equality QAP E A 1 α 2 α 2 α and C the following holds Thus, rank of the matrix A is a 1 There
10 54 I V Jovović, B J Malešević are regular matrices R and S such that RBS E 1 1 B holds Thus, rank of the matrix B is b 1 Since the ranks of the matrices A and B are a b 1, according to the Lemma 31 all entries of the last column and two rows of the matrix C QCS are zeros, ie we get that the matrix C is of the form C QCS Applying the Theorem 32, we obtain the general solution of the given matrix equation X P 2 α 1 β 1 γ 1 β 2 γ 2 R α β 1 γ 1 β 2 γ α 3β 1 2β 2 + 6γ 1 + 4γ 2 α 3γ 1 2γ 2 β 1 2γ 1 γ 1 β 2 2γ 2 γ 2 Acknowledgment Research is partially supported by the Ministry of Science and Education of the Republic of Serbia, Grant No17432 References 1 A BEN ISRAEL AND TNE GREVILLE, Generalized Inverses: Theory & Applications, Springer, New York, 23 2 A GRAHAM, Kronecker Products and Matrix Calculus with Applications, Ellis Horwood Lim, Chichester, England, B MALEŠEVIĆ, I JOVOVIĆ, M MAKRAGIĆ AND B RADIČIĆ, A note on solutions of linear systems, ISRN Algebra, Vol 213, Article ID , 6 pages (in the press) 4 B MALEŠEVIĆ AND B RADIČIĆ, Non-reproductive and reproductive solutions of some matrix equations, Proceedings of International Conference Mathematical and Informational Technologies 211, Vrnjačka Banja, 211, pp B MALEŠEVIĆ AND B RADIČIĆ, Reproductive and non-reproductive solutions of the matrix equation AXB C, Zbornik radova simpozijuma Matematika i primene, Matematički fakultet Beograd, 27 i 28 maj 211, pp B MALEŠEVIĆ AND B RADIČIĆ, Some considerations of matrix equations using the concept of reproductivity, Kragujevac J Math, Vol 36 (212,) No 1, pp
11 A note on solutions of the matrix equation AXB C 55 7 R PENROSE, A generalized inverses for matrices, Math Proc Cambridge Philos Soc, Vol 51 (1955), pp V PERIĆ, Generalized reciprocals of matrices, (in Serbo-Croatian), Matematika (Zagreb), Vol 11 (1982), No 1, pp B RADIČIĆ AND B MALEŠEVIĆ, Some considerations in relation to the matrix equation AXBC arxiv: CA ROHDE, Contribution to the theory, computation and application of generalized inverses, Doctoral dissertation, North Carolina State University at Raleigh, 1964
A note on solutions of linear systems
A note on solutions of linear systems Branko Malešević a, Ivana Jovović a, Milica Makragić a, Biljana Radičić b arxiv:134789v2 mathra 21 Jul 213 Abstract a Faculty of Electrical Engineering, University
More informationMultiplicative Perturbation Bounds of the Group Inverse and Oblique Projection
Filomat 30: 06, 37 375 DOI 0.98/FIL67M Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Multiplicative Perturbation Bounds of the Group
More informationRe-nnd solutions of the matrix equation AXB = C
Re-nnd solutions of the matrix equation AXB = C Dragana S. Cvetković-Ilić Abstract In this article we consider Re-nnd solutions of the equation AXB = C with respect to X, where A, B, C are given matrices.
More informationSome results on the reverse order law in rings with involution
Some results on the reverse order law in rings with involution Dijana Mosić and Dragan S. Djordjević Abstract We investigate some necessary and sufficient conditions for the hybrid reverse order law (ab)
More informationDiagonal and Monomial Solutions of the Matrix Equation AXB = C
Iranian Journal of Mathematical Sciences and Informatics Vol. 9, No. 1 (2014), pp 31-42 Diagonal and Monomial Solutions of the Matrix Equation AXB = C Massoud Aman Department of Mathematics, Faculty of
More informationChapter 1. Matrix Algebra
ST4233, Linear Models, Semester 1 2008-2009 Chapter 1. Matrix Algebra 1 Matrix and vector notation Definition 1.1 A matrix is a rectangular or square array of numbers of variables. We use uppercase boldface
More informationPROCEEDINGS OF THE YEREVAN STATE UNIVERSITY MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Mathematical Sciences 205, 2, p. 20 M a t h e m a t i c s MOORE PENROSE INVERSE OF BIDIAGONAL MATRICES. I Yu. R. HAKOPIAN, S. S. ALEKSANYAN 2, Chair
More informationMoore-Penrose-invertible normal and Hermitian elements in rings
Moore-Penrose-invertible normal and Hermitian elements in rings Dijana Mosić and Dragan S. Djordjević Abstract In this paper we present several new characterizations of normal and Hermitian elements in
More informationa Λ q 1. Introduction
International Journal of Pure and Applied Mathematics Volume 9 No 26, 959-97 ISSN: -88 (printed version); ISSN: -95 (on-line version) url: http://wwwijpameu doi: 272/ijpamv9i7 PAijpameu EXPLICI MOORE-PENROSE
More informationarxiv: v1 [math.ra] 28 Jan 2016
The Moore-Penrose inverse in rings with involution arxiv:1601.07685v1 [math.ra] 28 Jan 2016 Sanzhang Xu and Jianlong Chen Department of Mathematics, Southeast University, Nanjing 210096, China Abstract:
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 informationProceedings of the Third International DERIV/TI-92 Conference
Proceedings of the Third International DERIV/TI-9 Conference The Use of a Computer Algebra System in Modern Matrix Algebra Karsten Schmidt Faculty of Economics and Management Science, University of Applied
More informationApplied Matrix Algebra Lecture Notes Section 2.2. Gerald Höhn Department of Mathematics, Kansas State University
Applied Matrix Algebra Lecture Notes Section 22 Gerald Höhn Department of Mathematics, Kansas State University September, 216 Chapter 2 Matrices 22 Inverses Let (S) a 11 x 1 + a 12 x 2 + +a 1n x n = b
More informationMATRIX ALGEBRA AND SYSTEMS OF EQUATIONS. + + x 1 x 2. x n 8 (4) 3 4 2
MATRIX ALGEBRA AND SYSTEMS OF EQUATIONS SYSTEMS OF EQUATIONS AND MATRICES Representation of a linear system The general system of m equations in n unknowns can be written a x + a 2 x 2 + + a n x n b a
More informationMOORE-PENROSE INVERSE IN AN INDEFINITE INNER PRODUCT SPACE
J. Appl. Math. & Computing Vol. 19(2005), No. 1-2, pp. 297-310 MOORE-PENROSE INVERSE IN AN INDEFINITE INNER PRODUCT SPACE K. KAMARAJ AND K. C. SIVAKUMAR Abstract. The concept of the Moore-Penrose inverse
More informationDragan S. Djordjević. 1. Introduction and preliminaries
PRODUCTS OF EP OPERATORS ON HILBERT SPACES Dragan S. Djordjević Abstract. A Hilbert space operator A is called the EP operator, if the range of A is equal with the range of its adjoint A. In this article
More informationMATH 2331 Linear Algebra. Section 2.1 Matrix Operations. Definition: A : m n, B : n p. Example: Compute AB, if possible.
MATH 2331 Linear Algebra Section 2.1 Matrix Operations Definition: A : m n, B : n p ( 1 2 p ) ( 1 2 p ) AB = A b b b = Ab Ab Ab Example: Compute AB, if possible. 1 Row-column rule: i-j-th entry of AB:
More informationAPPLICATIONS OF THE HYPER-POWER METHOD FOR COMPUTING MATRIX PRODUCTS
Univ. Beograd. Publ. Eletrotehn. Fa. Ser. Mat. 15 (2004), 13 25. Available electronically at http: //matematia.etf.bg.ac.yu APPLICATIONS OF THE HYPER-POWER METHOD FOR COMPUTING MATRIX PRODUCTS Predrag
More informationEXPLICIT SOLUTION OF THE OPERATOR EQUATION A X + X A = B
EXPLICIT SOLUTION OF THE OPERATOR EQUATION A X + X A = B Dragan S. Djordjević November 15, 2005 Abstract In this paper we find the explicit solution of the equation A X + X A = B for linear bounded operators
More informationInverses and Elementary Matrices
Inverses and Elementary Matrices 1-12-2013 Matrix inversion gives a method for solving some systems of equations Suppose a 11 x 1 +a 12 x 2 + +a 1n x n = b 1 a 21 x 1 +a 22 x 2 + +a 2n x n = b 2 a n1 x
More informationSolutions to Exam I MATH 304, section 6
Solutions to Exam I MATH 304, section 6 YOU MUST SHOW ALL WORK TO GET CREDIT. Problem 1. Let A = 1 2 5 6 1 2 5 6 3 2 0 0 1 3 1 1 2 0 1 3, B =, C =, I = I 0 0 0 1 1 3 4 = 4 4 identity matrix. 3 1 2 6 0
More informationPreliminary Linear Algebra 1. Copyright c 2012 Dan Nettleton (Iowa State University) Statistics / 100
Preliminary Linear Algebra 1 Copyright c 2012 Dan Nettleton (Iowa State University) Statistics 611 1 / 100 Notation for all there exists such that therefore because end of proof (QED) Copyright c 2012
More informationOperators with Compatible Ranges
Filomat : (7), 579 585 https://doiorg/98/fil7579d Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://wwwpmfniacrs/filomat Operators with Compatible Ranges
More informationPartial isometries and EP elements in rings with involution
Electronic Journal of Linear Algebra Volume 18 Volume 18 (2009) Article 55 2009 Partial isometries and EP elements in rings with involution Dijana Mosic dragan@pmf.ni.ac.yu Dragan S. Djordjevic Follow
More informationGeneralized Principal Pivot Transform
Generalized Principal Pivot Transform M. Rajesh Kannan and R. B. Bapat Indian Statistical Institute New Delhi, 110016, India Abstract The generalized principal pivot transform is a generalization of the
More informationOn some matrices related to a tree with attached graphs
On some matrices related to a tree with attached graphs R. B. Bapat Indian Statistical Institute New Delhi, 110016, India fax: 91-11-26856779, e-mail: rbb@isid.ac.in Abstract A tree with attached graphs
More informationThe DMP Inverse for Rectangular Matrices
Filomat 31:19 (2017, 6015 6019 https://doi.org/10.2298/fil1719015m Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://.pmf.ni.ac.rs/filomat The DMP Inverse for
More informationELA THE OPTIMAL PERTURBATION BOUNDS FOR THE WEIGHTED MOORE-PENROSE INVERSE. 1. Introduction. Let C m n be the set of complex m n matrices and C m n
Electronic Journal of Linear Algebra ISSN 08-380 Volume 22, pp. 52-538, May 20 THE OPTIMAL PERTURBATION BOUNDS FOR THE WEIGHTED MOORE-PENROSE INVERSE WEI-WEI XU, LI-XIA CAI, AND WEN LI Abstract. In this
More informationElementary maths for GMT
Elementary maths for GMT Linear Algebra Part 2: Matrices, Elimination and Determinant m n matrices The system of m linear equations in n variables x 1, x 2,, x n a 11 x 1 + a 12 x 2 + + a 1n x n = b 1
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationEP elements in rings
EP elements in rings Dijana Mosić, Dragan S. Djordjević, J. J. Koliha Abstract In this paper we present a number of new characterizations of EP elements in rings with involution in purely algebraic terms,
More informationCalculus II - Basic Matrix Operations
Calculus II - Basic Matrix Operations Ryan C Daileda Terminology A matrix is a rectangular array of numbers, for example 7,, 7 7 9, or / / /4 / / /4 / / /4 / /6 The numbers in any matrix are called its
More informationProduct of Range Symmetric Block Matrices in Minkowski Space
BULLETIN of the Malaysian Mathematical Sciences Society http://math.usm.my/bulletin Bull. Malays. Math. Sci. Soc. (2) 29(1) (2006), 59 68 Product of Range Symmetric Block Matrices in Minkowski Space 1
More informationSection 3.3. Matrix Rank and the Inverse of a Full Rank Matrix
3.3. Matrix Rank and the Inverse of a Full Rank Matrix 1 Section 3.3. Matrix Rank and the Inverse of a Full Rank Matrix Note. The lengthy section (21 pages in the text) gives a thorough study of the rank
More informationDeterminants Chapter 3 of Lay
Determinants Chapter of Lay Dr. Doreen De Leon Math 152, Fall 201 1 Introduction to Determinants Section.1 of Lay Given a square matrix A = [a ij, the determinant of A is denoted by det A or a 11 a 1j
More informationLinear Algebra Math 221
Linear Algebra Math Open Book Exam Open Notes 8 Oct, 004 Calculators Permitted Show all work (except #4). (0 pts) Let A = 3 a) (0 pts) Compute det(a) by Gaussian Elimination. 3 3 swap(i)&(ii) (iii) (iii)+(
More informationElementary Matrices. MATH 322, Linear Algebra I. J. Robert Buchanan. Spring Department of Mathematics
Elementary Matrices MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Outline Today s discussion will focus on: elementary matrices and their properties, using elementary
More informationarxiv: v1 [math.ca] 2 Nov 2018
arxiv.org arxiv:1811.00748v1 [math.ca] 2 Nov 2018 One method for proving some classes of eponential analytical inequalities Branko Malešević a, Tatjana Lutovac a, Bojan Banjac b a School of Electrical
More informationElementary operation matrices: row addition
Elementary operation matrices: row addition For t a, let A (n,t,a) be the n n matrix such that { A (n,t,a) 1 if r = c, or if r = t and c = a r,c = 0 otherwise A (n,t,a) = I + e t e T a Example: A (5,2,4)
More informationarxiv: v1 [math.ra] 21 Jul 2013
Projections and Idempotents in -reducing Rings Involving the Moore-Penrose Inverse arxiv:1307.5528v1 [math.ra] 21 Jul 2013 Xiaoxiang Zhang, Shuangshuang Zhang, Jianlong Chen, Long Wang Department of Mathematics,
More informationLinear Algebra Review
Linear Algebra Review ORIE 4741 September 1, 2017 Linear Algebra Review September 1, 2017 1 / 33 Outline 1 Linear Independence and Dependence 2 Matrix Rank 3 Invertible Matrices 4 Norms 5 Projection Matrix
More informationNotes on Row Reduction
Notes on Row Reduction Francis J. Narcowich Department of Mathematics Texas A&M University September The Row-Reduction Algorithm The row-reduced form of a matrix contains a great deal of information, both
More informationTopics. Vectors (column matrices): Vector addition and scalar multiplication The matrix of a linear function y Ax The elements of a matrix A : A ij
Topics Vectors (column matrices): Vector addition and scalar multiplication The matrix of a linear function y Ax The elements of a matrix A : A ij or a ij lives in row i and column j Definition of a matrix
More informationLinear Systems and Matrices
Department of Mathematics The Chinese University of Hong Kong 1 System of m linear equations in n unknowns (linear system) a 11 x 1 + a 12 x 2 + + a 1n x n = b 1 a 21 x 1 + a 22 x 2 + + a 2n x n = b 2.......
More informationChapter 2. Ma 322 Fall Ma 322. Sept 23-27
Chapter 2 Ma 322 Fall 2013 Ma 322 Sept 23-27 Summary ˆ Matrices and their Operations. ˆ Special matrices: Zero, Square, Identity. ˆ Elementary Matrices, Permutation Matrices. ˆ Voodoo Principle. What is
More informationMath 4377/6308 Advanced Linear Algebra
2.3 Composition Math 4377/6308 Advanced Linear Algebra 2.3 Composition of Linear Transformations Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu math.uh.edu/ jiwenhe/math4377
More informationarxiv: v1 [math.ra] 23 Feb 2018
JORDAN DERIVATIONS ON SEMIRINGS OF TRIANGULAR MATRICES arxiv:180208704v1 [mathra] 23 Feb 2018 Abstract Dimitrinka Vladeva University of forestry, bulklohridski 10, Sofia 1000, Bulgaria E-mail: d vladeva@abvbg
More informationOn Regularity of Incline Matrices
International Journal of Algebra, Vol. 5, 2011, no. 19, 909-924 On Regularity of Incline Matrices A. R. Meenakshi and P. Shakila Banu Department of Mathematics Karpagam University Coimbatore-641 021, India
More informationLectures on Linear Algebra for IT
Lectures on Linear Algebra for IT by Mgr Tereza Kovářová, PhD following content of lectures by Ing Petr Beremlijski, PhD Department of Applied Mathematics, VSB - TU Ostrava Czech Republic 3 Inverse Matrix
More informationarxiv: v5 [math.ra] 1 Jan 2013
Differential transcendency in the theory of linear differential systems with constant coefficients arxiv:0811.0144v5 [math.ra] 1 Jan 2013 Branko Malešević a, Dragana Todorić b, Ivana Jovović a, Sonja Telebaković
More informationRecall the convention that, for us, all vectors are column vectors.
Some linear algebra Recall the convention that, for us, all vectors are column vectors. 1. Symmetric matrices Let A be a real matrix. Recall that a complex number λ is an eigenvalue of A if there exists
More informationMinkowski Inverse for the Range Symmetric Block Matrix with Two Identical Sub-blocks Over Skew Fields in Minkowski Space M
International Journal of Mathematics And its Applications Volume 4, Issue 4 (2016), 67 80. ISSN: 2347-1557 Available Online: http://ijmaa.in/ International Journal 2347-1557 of Mathematics Applications
More informationMath 344 Lecture # Linear Systems
Math 344 Lecture #12 2.7 Linear Systems Through a choice of bases S and T for finite dimensional vector spaces V (with dimension n) and W (with dimension m), a linear equation L(v) = w becomes the linear
More informationMatrix Transformations and Statistical Convergence II
Advances in Dynamical Systems and Applications ISSN 0973-532, Volume 6, Number, pp. 7 89 20 http://campus.mst.edu/adsa Matrix Transformations and Statistical Convergence II Bruno de Malafosse LMAH Université
More informationMath 370 Spring 2016 Sample Midterm with Solutions
Math 370 Spring 2016 Sample Midterm with Solutions Contents 1 Problems 2 2 Solutions 5 1 1 Problems (1) Let A be a 3 3 matrix whose entries are real numbers such that A 2 = 0. Show that I 3 + A is invertible.
More informationTHE OHIO JOURNAL OF SCIENCE
THE OHIO JOURNAL OF SCIENCE Vol. 74 SEPTEMBER, 1974 No. 5 QUADRATIC MATRIX EQUATIONS 1 G. L. SHURBET, T. O. LEWIS, AND T. L. BOULLION Department of Mathematics, Texas Tech University, Lubbock, Texas 791+09
More informationMatrix Algebra. Matrix Algebra. Chapter 8 - S&B
Chapter 8 - S&B Algebraic operations Matrix: The size of a matrix is indicated by the number of its rows and the number of its columns. A matrix with k rows and n columns is called a k n matrix. The number
More informationRepresentations of the (b, c)-inverses in Rings with Involution
Filomat 31:9 (217), 2867 2875 htts://doiorg/12298/fil179867k Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: htt://wwwmfniacrs/filomat Reresentations of the (b,
More informationE k E k 1 E 2 E 1 A = B
Theorem.5. suggests that reducing a matrix A to (reduced) row echelon form is tha same as multiplying A from left by the appropriate elementary matrices. Hence if B is a matrix obtained from a matrix A
More information7 Matrix Operations. 7.0 Matrix Multiplication + 3 = 3 = 4
7 Matrix Operations Copyright 017, Gregory G. Smith 9 October 017 The product of two matrices is a sophisticated operations with a wide range of applications. In this chapter, we defined this binary operation,
More informationMath 224, Fall 2007 Exam 3 Thursday, December 6, 2007
Math 224, Fall 2007 Exam 3 Thursday, December 6, 2007 You have 1 hour and 20 minutes. No notes, books, or other references. You are permitted to use Maple during this exam, but you must start with a blank
More informationMath 407: Linear Optimization
Math 407: Linear Optimization Lecture 16: The Linear Least Squares Problem II Math Dept, University of Washington February 28, 2018 Lecture 16: The Linear Least Squares Problem II (Math Dept, University
More informationTHE MOORE-PENROSE GENERALIZED INVERSE OF A MATRIX
THE MOORE-PENROSE GENERALIZED INVERSE OF A MATRIX A Dissertation Submitted For The Award of the Degree of Master of Philosophy In Mathematics Purva Rajwade School of Mathematics Devi Ahilya Vishwavidyalaya,
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 informationSTAT 309: MATHEMATICAL COMPUTATIONS I FALL 2017 LECTURE 5
STAT 39: MATHEMATICAL COMPUTATIONS I FALL 17 LECTURE 5 1 existence of svd Theorem 1 (Existence of SVD) Every matrix has a singular value decomposition (condensed version) Proof Let A C m n and for simplicity
More informationRepeated Eigenvalues and Symmetric Matrices
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationMoore Penrose inverses and commuting elements of C -algebras
Moore Penrose inverses and commuting elements of C -algebras Julio Benítez Abstract Let a be an element of a C -algebra A satisfying aa = a a, where a is the Moore Penrose inverse of a and let b A. We
More informationWeaker assumptions for convergence of extended block Kaczmarz and Jacobi projection algorithms
DOI: 10.1515/auom-2017-0004 An. Şt. Univ. Ovidius Constanţa Vol. 25(1),2017, 49 60 Weaker assumptions for convergence of extended block Kaczmarz and Jacobi projection algorithms Doina Carp, Ioana Pomparău,
More informationProblem # Max points possible Actual score Total 120
FINAL EXAMINATION - MATH 2121, FALL 2017. Name: ID#: Email: Lecture & Tutorial: Problem # Max points possible Actual score 1 15 2 15 3 10 4 15 5 15 6 15 7 10 8 10 9 15 Total 120 You have 180 minutes to
More informationMath 1553 Introduction to Linear Algebra
Math 1553 Introduction to Linear Algebra Lecture Notes Chapter 2 Matrix Algebra School of Mathematics The Georgia Institute of Technology Math 1553 Lecture Notes for Chapter 2 Introduction, Slide 1 Section
More informationMath 54 First Midterm Exam, Prof. Srivastava September 23, 2016, 4:10pm 5:00pm, 155 Dwinelle Hall.
Math 54 First Midterm Exam, Prof Srivastava September 23, 26, 4:pm 5:pm, 55 Dwinelle Hall Name: SID: Instructions: Write all answers in the provided space This exam includes two pages of scratch paper,
More informationLecture Notes in Linear Algebra
Lecture Notes in Linear Algebra Dr. Abdullah Al-Azemi Mathematics Department Kuwait University February 4, 2017 Contents 1 Linear Equations and Matrices 1 1.2 Matrices............................................
More informationAnalytical formulas for calculating extremal ranks and inertias of quadratic matrix-valued functions and their applications
Analytical formulas for calculating extremal ranks and inertias of quadratic matrix-valued functions and their applications Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081,
More informationORIE 6300 Mathematical Programming I August 25, Recitation 1
ORIE 6300 Mathematical Programming I August 25, 2016 Lecturer: Calvin Wylie Recitation 1 Scribe: Mateo Díaz 1 Linear Algebra Review 1 1.1 Independence, Spanning, and Dimension Definition 1 A (usually infinite)
More informationInner image-kernel (p, q)-inverses in rings
Inner image-kernel (p, q)-inverses in rings Dijana Mosić Dragan S. Djordjević Abstract We define study the inner image-kernel inverse as natural algebraic extension of the inner inverse with prescribed
More informationMA 138 Calculus 2 with Life Science Applications Matrices (Section 9.2)
MA 38 Calculus 2 with Life Science Applications Matrices (Section 92) Alberto Corso albertocorso@ukyedu Department of Mathematics University of Kentucky Friday, March 3, 207 Identity Matrix and Inverse
More informationJournal of Symbolic Computation. On the Berlekamp/Massey algorithm and counting singular Hankel matrices over a finite field
Journal of Symbolic Computation 47 (2012) 480 491 Contents lists available at SciVerse ScienceDirect Journal of Symbolic Computation journal homepage: wwwelseviercom/locate/jsc On the Berlekamp/Massey
More informationA revisit to a reverse-order law for generalized inverses of a matrix product and its variations
A revisit to a reverse-order law for generalized inverses of a matrix product and its variations Yongge Tian CEMA, Central University of Finance and Economics, Beijing 100081, China Abstract. For a pair
More information1111: Linear Algebra I
1111: Linear Algebra I Dr. Vladimir Dotsenko (Vlad) Michaelmas Term 2015 Dr. Vladimir Dotsenko (Vlad) 1111: Linear Algebra I Michaelmas Term 2015 1 / 10 Row expansion of the determinant Our next goal is
More informationAppendix A: Matrices
Appendix A: Matrices A matrix is a rectangular array of numbers Such arrays have rows and columns The numbers of rows and columns are referred to as the dimensions of a matrix A matrix with, say, 5 rows
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 informationELA
Electronic Journal of Linear Algebra ISSN 181-81 A publication of te International Linear Algebra Society ttp://mat.tecnion.ac.il/iic/ela RANK AND INERTIA OF SUBMATRICES OF THE MOORE PENROSE INVERSE OF
More information~ g-inverses are indeed an integral part of linear algebra and should be treated as such even at an elementary level.
Existence of Generalized Inverse: Ten Proofs and Some Remarks R B Bapat Introduction The theory of g-inverses has seen a substantial growth over the past few decades. It is an area of great theoretical
More informationOn Euclidean distance matrices
On Euclidean distance matrices R. Balaji and R. B. Bapat Indian Statistical Institute, New Delhi, 110016 November 19, 2006 Abstract If A is a real symmetric matrix and P is an orthogonal projection onto
More informationPredrag S. Stanimirović and Dragan S. Djordjević
FULL-RANK AND DETERMINANTAL REPRESENTATION OF THE DRAZIN INVERSE Predrag S. Stanimirović and Dragan S. Djordjević Abstract. In this article we introduce a full-rank representation of the Drazin inverse
More informationJournal of Complexity. On strata of degenerate polyhedral cones, II: Relations between condition measures
Journal of Complexity 26 (200) 209 226 Contents lists available at ScienceDirect Journal of Complexity journal homepage: www.elsevier.com/locate/jco On strata of degenerate polyhedral cones, II: Relations
More informationMatrix Operations. Linear Combination Vector Algebra Angle Between Vectors Projections and Reflections Equality of matrices, Augmented Matrix
Linear Combination Vector Algebra Angle Between Vectors Projections and Reflections Equality of matrices, Augmented Matrix Matrix Operations Matrix Addition and Matrix Scalar Multiply Matrix Multiply Matrix
More informationSpectral inequalities and equalities involving products of matrices
Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department
More informationMATH36001 Generalized Inverses and the SVD 2015
MATH36001 Generalized Inverses and the SVD 201 1 Generalized Inverses of Matrices A matrix has an inverse only if it is square and nonsingular. However there are theoretical and practical applications
More informationA system of matrix equations and a linear matrix equation over arbitrary regular rings with identity
Linear Algebra and its Applications 384 2004) 43 54 www.elsevier.com/locate/laa A system of matrix equations and a linear matrix equation over arbitrary regular rings with identity Qing-Wen Wang Department
More informationDigital Workbook for GRA 6035 Mathematics
Eivind Eriksen Digital Workbook for GRA 6035 Mathematics November 10, 2014 BI Norwegian Business School Contents Part I Lectures in GRA6035 Mathematics 1 Linear Systems and Gaussian Elimination........................
More informationLecture 1 Systems of Linear Equations and Matrices
Lecture 1 Systems of Linear Equations and Matrices Math 19620 Outline of Course Linear Equations and Matrices Linear Transformations, Inverses Bases, Linear Independence, Subspaces Abstract Vector Spaces
More informationMatrices and systems of linear equations
Matrices and systems of linear equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T.
More informationElementary Row Operations on Matrices
King Saud University September 17, 018 Table of contents 1 Definition A real matrix is a rectangular array whose entries are real numbers. These numbers are organized on rows and columns. An m n matrix
More informationMATH2210 Notebook 2 Spring 2018
MATH2210 Notebook 2 Spring 2018 prepared by Professor Jenny Baglivo c Copyright 2009 2018 by Jenny A. Baglivo. All Rights Reserved. 2 MATH2210 Notebook 2 3 2.1 Matrices and Their Operations................................
More informationSystems of Linear Equations
Systems of Linear Equations Math 108A: August 21, 2008 John Douglas Moore Our goal in these notes is to explain a few facts regarding linear systems of equations not included in the first few chapters
More informationMath Bootcamp An p-dimensional vector is p numbers put together. Written as. x 1 x =. x p
Math Bootcamp 2012 1 Review of matrix algebra 1.1 Vectors and rules of operations An p-dimensional vector is p numbers put together. Written as x 1 x =. x p. When p = 1, this represents a point in the
More informationMATH 300, Second Exam REVIEW SOLUTIONS. NOTE: You may use a calculator for this exam- You only need something that will perform basic arithmetic.
MATH 300, Second Exam REVIEW SOLUTIONS NOTE: You may use a calculator for this exam- You only need something that will perform basic arithmetic. [ ] [ ] 2 2. Let u = and v =, Let S be the parallelegram
More informationThe Determinant: a Means to Calculate Volume
The Determinant: a Means to Calculate Volume Bo Peng August 16, 2007 Abstract This paper gives a definition of the determinant and lists many of its well-known properties Volumes of parallelepipeds are
More informationOn the Berlekamp/Massey Algorithm and Counting Singular Hankel Matrices over a Finite Field
On the Berlekamp/Massey Algorithm and Counting Singular Hankel Matrices over a Finite Field Matthew T Comer Dept of Mathematics, North Carolina State University Raleigh, North Carolina, 27695-8205 USA
More information