The Hermitian R-symmetric Solutions of the Matrix Equation AXA = B
|
|
- Arnold Rose
- 5 years ago
- Views:
Transcription
1 International Journal of Algebra, Vol. 6, 0, no. 9, The Hermitian R-symmetric Solutions of the Matrix Equation AXA = B Qingfeng Xiao Department of Basic Dongguan olytechnic Dongguan 53808, China qfxiao@hnu.edu.cn Abstract In this paper, we mainly discuss solving the following problems. roblem I. Given matrices A C m n and B C m m, find X HRS n n such that AXA = B, where HRS n n = {X H n n RXR = X, for given R C n n satisfying R = R = R I n }. roblem II. Given a matrix X C n n, find ˆX S E such that X ˆX F = inf X X F, X S E where is the Frobenius norm, and S E is the solution set of roblem I. Expressions for the general solution of roblem I are derived. Necessary and sufficient conditions for the solvability of roblem I are provided. For roblem II, an expression for the solution is given as well. Mathematics Subject Classification: 65F5, 65F0 Keywords: Matrix equation, Hermitian R-symmetric, Inverse problem, The Optimal approximation Introduction Throughout this paper, let C n m be the set of all n m complex matrices, H n n denote the class of n n hermitian matrices, U n n be the set of all n n unitary matrices. Denote by I n the identity matrix with order n. Let J = (e n,e n,...,e ), where e i is the ith column of I n. For matrix A, A, A +, A F and r(a) represent its conjugate transpose, Moore-enrose inverse, Frobenius norm and rank, respectively. For matrices A, B C n m, the expression A B will be the Hadamard product of A and B. Defining the inner product (A, B) =
2 904 Qingfeng Xiao tr(b A) for matrices A, B C n m, C n m becomes a Hilbert space. The norm of a matrix generated by this inner product is the Frobenius norm. Consider the linear matrix equation AXA = B, () where A C m n, B C m m are given, and X C n n is unknown(matrix). Such a matrix equation with symmetric pattern appears frequently in matrix theory and applications such as in statistics, vibration theory, etc. For important results on the inverse problem () associated with several kinds of different sets S, for instance, symmetric matrices, symmetric nonnegative definite matrices, bisymmetric (same as persymmetric) matrices, bisymmetric nonnegative definite matrices and so on, we refer the reader to [-8. In this paper we discuss the above inverse problem () associated with a new S. Definition. Let R C n n be a nontrivial unitary involution; i.e., R = R = R I n. We say that X C n n is a Hermitian R-symmetric matrix, if X = X, RXR = X. We denoted by HRS n n the set of all n n Hermitian R-symmetric matrices. The Hermitian R-symmetric is a class of important matrices and have engineering and scientific applications. For the case the unknown A is Hermitian R-symmetric, [9 has discussed the inverse problems AX = B. However, for this case, the inverse problem () has not been dealt with yet. This problem will be considered here. The problem studied in this paper can now be described as follows. roblem I. Given matrices A C m n and B C m m, find a Hermitian R-symmetric matrix X such that AXA = B. In this paper, we discuss the solvability of this problem and an expression for its solution is presented. The optimal approximation problem of a matrix with the above-given matrix restriction comes up in the processes of test or recovery of a linear system due to incomplete data or revising given data. A preliminary estimate X of the unknown matrix X can be obtained by the experimental observation values and the information of statistical distribution. The optimal estimate of X is a matrix ˆX that satisfies the given matrix restriction for X and is the best approximation of X, see [0-. In this paper, we will also consider the so-called optimal approximation problem associated with (). It reads as follows.
3 Hermitian R-symmetric solutions 905 roblem II. Given matrix X C n n, find ˆX S E such that X ˆX F = inf X S E X X F, where S E is the solution set of roblem I. We point out that if roblem I is solvable, then roblem II has a unique solution, and in this case an expression for the solution can be derived. The paper is organized as follows. In Section we first discuss some properties of Hermitian R-symmetric matrices. Then we establish the solvability conditions for roblem I and give an expression of the general solution of roblem I. In Section 3 we prove the existence and uniqueness of the solution to roblem II and derive an expression. The expression of the general solution of problem I We first know that hermitian R-symmetric matrices have the following properties. Since R = R = R I n, the only possible eigenvalues of R are + and. Let r and s be respectively the dimensions of the eigenspaces of R associated with the eigenvalues λ =,andλ = ; thus r, s > and r+s = n. Let =[p p r, and Q =[q q s, () where {p,...,p r } and {q,...,q s } are orthonormal bases for the eigenspaces. and Q can be found by applying the Gram-Schmidt process to the columns of I + R and I R, respectively. If R = J the standard choices for and Q are = [ Im J m, and Q = [ Im J m if n =m, or = I m 0 0, and Q = J m 0, I m 0 J m, if n =m +. In general, and R = and RQ = Q, (3) [ [ Q = Q.
4 906 Qingfeng Xiao Therefore any X C n n can be written conformably in block form as [ [ X C X =[ Q, (4) C with X = X H r r, X Q = Q XQ H s s, and C = XQ C r s, From (3) and (4), [ [ X C RXR =[ Q C, (5) X Q Lemma. [9 X C n n is Hermitian R-symmetric if and only if [ [ X 0 X =[ Q 0 X Q Q, (6) X Q with X = X H r r, X Q = Q XQ H s s. Let A C m n and B C m m,, Q with the form of (). Let A[ Q=[A A, A C m r, A C m s. (7) By making generalized singular value decomposition (see [3) to the matrix pair [A, A, we have Q Q A = MΣ U, A = MΣ V, (8) where M C m m is a nonsingular matrix, U U r r, V U s s and I 0 0 l k l k 0 S Σ = 0 k 0 S, Σ t r = 0 k, (9) 0 0 I t r n t n t t = r([a,a ), l = r(a ), k = r(a )+r(a ) r([a,a ), S = diag(α,...,α k ), S = diag(β,...,β k ), 0 <α k α <, 0 <β β k <, α i + β i =,i =,...,k. Now we can establish the existence theorems as follows. Theorem. Given A C m n and B C m m,, Q defined in(), and A[ Q has the partition form of (7), the generalized singular value decomposition of the matrix pair [A,A as (8). artition the matrix M BM as B B B 3 B 4 l k M BM B = B B 3 B 4 k, (0) B 3 B 3 B 33 B 34 t r B 4 B 4 B 43 B 44 m t
5 Hermitian R-symmetric solutions 907 then the problem I has a solution X HRS n n if and only if B = B,B 3 = O, B 4 = O, B 4 = O, B 43 = O, B 44 = O. () In that case it has the general solution [ X 0 X =[ Q 0 X Q [ Q, () where X = U B B S X 3 S B X X 3 X 3 X 3 X 33 Y Y Y 3 X Q = V Y S (B S X S )S S B 3 Y3 B 3 S B 33 U, (3) V, (4) with X 3 C (l k) (r l), X H k k, X 3 C k (r l), X 33 C (r l) (r l), Y H (s+l k t) (s+l k t), Y C (s+l k t) k, Y 3 C (s+l k t) (t l) are arbitrary matrices. roof The necessity. Assume Eq.() has a solution X HRS n n. By the definition of HRS n n, it is easy to verify that B = B, and we have from Lemma that X can be expressed as [ [ X 0 X =[ Q 0 X Q Q, (5) where X H r r, X Q H s s. Note that [ Q is an unitary matrix, and the definition of A i (i =, ), Eq.() is equivalent to Substituting (8) into (6), then we have A X A + A X Q A = B. (6) Σ U X UΣ +Σ V X Q V Σ = M BM. (7) artition the matrices U X U, V X Q V as follows: U X U = X X X 3 X X X 3 X 3 X 3 X 33, V X Q V = Y Y Y 3 Y Y Y 3 Y 3 Y 3 Y 33, (8) where X C (l k) (l k), X C k k, X 33 C (r l) (r l), and Y C (s+l k t) (s+l k t), Y C k k, Y 33 C (t l) (t l).
6 908 Qingfeng Xiao Substituting (9), (0) and (8) into (7), we have X X S 0 0 B B B 3 B 4 S X S X S + S Y S S Y Y 3 S Y 33 0 = B B B 3 B 4 B 3 B 3 B 33 B B 4 B 4 B 43 B 44. (9) Therefore (9) holds if and only if () holds and X = B, X = X = B S, Y 3 = Y3 = S B 3, Y 33 = B 33, Y = S (B S X S )S where X H k k is arbitrary matrix. Substituting the above into (8) and (5), thus we have formulation (), (3) and (4). The sufficiency. Assume () holds. Let X = U B B S X 3 S B X X 3 X 3 X 3 X 33 U, Y Y Y 3 X Q = V Y S (B S X S )S S B 3 V, Y3 B 3 S B 33 where X 3 C (l k) (r l), X H k k, X 3 C k (r l), X 33 C (r l) (r l), Y H (s+l k t) (s+l k t), Y C (s+l k t) k, Y 3 C (s+l k t) (t l) are arbitrary matrices. Obviously, X H r r, XQ H s s. By Lemma. and X =[ Q [ X 0 0 XQ [ we have X HRS n n. Hence A XA = A[ [ [ X 0 Q 0 XQ Q A = A X A + A XQ A = B B B M S X S B S X S B B3 B = B B 0 0 B M B B B3 B 33 0 M = B [ [ This implies that X X 0 =[ Q 0 XQ Q is the Hermitian R- symmetric solution of Eq.(). The proof is completed. Q,
7 Hermitian R-symmetric solutions The expression of the solution of roblem II To prepare for an explicit expression for the solution of the matrix nearness problem II, we first verify the following lemma. Lemma 3. [6 Suppose that G, H C n n, S σ = diag(σ,,σ n ), σ i > 0(i =,,,n), Then there exist a unique Ŝ Hn n such that S G F + S σss σ H F = min (0) and where Ŝ =Φ [G + G + S σ (H + H )S σ, Φ=(ψ ij ) C n n,ψ ij =, i, j n. ( + σi σj ) Theorem 3. Let the matrix X C n n, the other relevant [ conditions and symbols be the same as those stated in Theorem.. artition X[ Q, U Z U, V Z V as follows: U Z U = [ Q X[ X X X3 X X X3 X 3 X3 X33 Q= [ Z Z Z, V Z V = Z Q, () Ỹ Ỹ Ỹ 3 Ỹ Ỹ Ỹ 3 Ỹ 3 Ỹ 3 Ỹ 33. () When the solution set S E of roblem I is nonempty, then roblem II has a unique solution ˆX, which can be expressed as where ˆX = U ˆX =[ Q [ ˆX 0 0 ˆXQ [ Q B B S ( X 3 + X 3 ) S B K ( X 3 + X 3 ) ( X 3 + X 3) ( X 3 + X 3) ( X 33 + X 33), (3) U, (4) ˆX (Ỹ + Ỹ ) (Ỹ + Ỹ ) (Ỹ3 + Ỹ 3) = V (Ỹ + Ỹ ) S (B S KS )S S B 3 V, (5) (Ỹ3 + Ỹ 3 ) B 3S B 33
8 90 Qingfeng Xiao with K =Φ [ X + X + S S (S B S Ỹ Ỹ )S S, (6) Φ=(ψ ij ) R k k, ψ ij = βi β j, i, j k. (αi αj + βi βj ) roof When S m is nonempty, it is easy to verify from () that S m is a closed convex set, there exists a unique solution for roblem II. Using the invariance of the Frobenius norm under unitary transformations, from (), () and () we have X X F = X Z F + X Z F + Z F + Z F = + U Z U + Z F + Z F F Y Y Y 3 Y S (B S X S )S S B 3 V Z V. Y3 B 3 S B 33 F B B S X 3 S B X X 3 X3 X3 X 33 Thus X ˆX F = inf X S E X X F is equivalent to X 3 X 3 F + X 3 X 3 F = min, X 3 X 3 F + X 3 X 3 F = min, X 33 X 33 F = min, Y Ỹ F = min, Y Ỹ F + Y +Ỹ F = min, Y 3 Ỹ3 F + Y 3 Ỹ3 F = min, X X F + S S X S S (S B S Ỹ) F = min. From Lemma 3. we have X 3 = ( X 3 + X 3 ), X 3 = ( X 3 + X 3 ), X 33 = ( X 33 + X 33 ),Y = (Ỹ + Ỹ ),Y = (Ỹ + Ỹ ),Y 3 = (Ỹ3 + Ỹ 3 ) and X =Φ [ X + X + S S (S B S Ỹ Ỹ )S S. Taking X i,y i3,i =,, 3 and X into (), (3), (4), we obtain that the solution of (the matrix nearness) roblem II can be expressed as (3). The proof is completed. ACKNOWLEDGEMENTS. This research was supported by Scientific Research Fund of Dongguan olytechnic (0a5).
9 Hermitian R-symmetric solutions 9 References [ J. K. Baksalary, Nonnegative definite and positive definite solutions to the matrix equation AXA = B, Linear and Multilinear Algebra, 6 (984) [ J. Groß, A note on the general Hermitian solution to AXA = B, Bull. Malaysian Math. Sci. Soc., (998) [3 J. Groß, Nonnegtive-definite and positive-definite solution to the matrix equation AXA = B, Linear Algebra Appl., 3 (000) 3 9. [4 X. Zhang and M. Cheng, The rank-constrained Hermitian nonnegtivedefinite and positive-definite solutions to the matrix equation AXA = B, Linear Algebra Appl., 370 (003) 63C74. [5 Y. Liu, Y. Tian and Y. Takane, Ranks of Hermitian and skew-hermitian solutions to the matrix equation AXA = B, Linear Algebra Appl., 43 (009) [6 Z. eng, X. Hu, and L. Zhang, One kind of inverse problem for the bisymmetric matrices, Math. Numer. Sinica., 7 (005) [7 C. Guo, X. Hu, and L. Zhang, A class of inverse problems of matrix equation X T AX = B, Numerical Mathematics A Journal of Chiness University, 6 (004) -9. [8 Q.F. Xiao, X.Y. Hu, and L. Zhang, The Anti-symmetric Ortho-symmetric solutions of the matrix equation A T XA = D, Electronic Journal of Linear Algebra,8, (009), - 9. [9 W.F. Trench, Hermitian, hermitian R-symmetric, and hermitian R-skew symmetric rocrustes problems, Linear Algebra Appl., 387 (004) [0 Z. Xu, K. Zhang and Q. Lu, A Fast Algorithm for a Class of Toeplitz Matrices (in Chinese), Xibei University of Industry ress, 999. [ N.J. Highan, Computing a nearest symmetric positive semidefinite matrix, Linear Algebra Appl., 03 (988) [ L. Zhang, The approximation on the closed convex cone and its numerical application, Hunan Annals of Mathematics, 6 (986) [3 C.C. aige and M.A. Saunders, Towards a Generalized Singular Value Decomposition, SIAM J.Numer. Anal., 8 (98) Received: May, 0
The Skew-Symmetric Ortho-Symmetric Solutions of the Matrix Equations A XA = D
International Journal of Algebra, Vol. 5, 2011, no. 30, 1489-1504 The Skew-Symmetric Ortho-Symmetric Solutions of the Matrix Equations A XA = D D. Krishnaswamy Department of Mathematics Annamalai University
More informationThe reflexive and anti-reflexive solutions of the matrix equation A H XB =C
Journal of Computational and Applied Mathematics 200 (2007) 749 760 www.elsevier.com/locate/cam The reflexive and anti-reflexive solutions of the matrix equation A H XB =C Xiang-yang Peng a,b,, Xi-yan
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationThe reflexive re-nonnegative definite solution to a quaternion matrix equation
Electronic Journal of Linear Algebra Volume 17 Volume 17 28 Article 8 28 The reflexive re-nonnegative definite solution to a quaternion matrix equation Qing-Wen Wang wqw858@yahoo.com.cn Fei Zhang Follow
More informationThe skew-symmetric orthogonal solutions of the matrix equation AX = B
Linear Algebra and its Applications 402 (2005) 303 318 www.elsevier.com/locate/laa The skew-symmetric orthogonal solutions of the matrix equation AX = B Chunjun Meng, Xiyan Hu, Lei Zhang College of Mathematics
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 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 informationAN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES
AN INVERSE EIGENVALUE PROBLEM AND AN ASSOCIATED APPROXIMATION PROBLEM FOR GENERALIZED K-CENTROHERMITIAN MATRICES ZHONGYUN LIU AND HEIKE FAßBENDER Abstract: A partially described inverse eigenvalue problem
More informationThe symmetric minimal rank solution of the matrix equation AX=B and the optimal approximation
Electronic Journal of Linear Algebra Volume 18 Volume 18 (2009 Article 23 2009 The symmetric minimal rank solution of the matrix equation AX=B and the optimal approximation Qing-feng Xiao qfxiao@hnu.cn
More informationTHE INVERSE PROBLEM OF CENTROSYMMETRIC MATRICES WITH A SUBMATRIX CONSTRAINT 1) 1. Introduction
Journal of Computational Mathematics, Vol22, No4, 2004, 535 544 THE INVERSE PROBLEM OF CENTROSYMMETRIC MATRICES WITH A SUBMATRIX CONSTRAINT 1 Zhen-yun Peng Department of Mathematics, Hunan University of
More informationInverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with J-(Skew) Centrosymmetry
Inverse Eigenvalue Problems and Their Associated Approximation Problems for Matrices with -(Skew) Centrosymmetry Zhong-Yun Liu 1 You-Cai Duan 1 Yun-Feng Lai 1 Yu-Lin Zhang 1 School of Math., Changsha University
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 informationRank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications
Rank and inertia optimizations of two Hermitian quadratic matrix functions subject to restrictions with applications Yongge Tian a, Ying Li b,c a China Economics and Management Academy, Central 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 informationYongge Tian. China Economics and Management Academy, Central University of Finance and Economics, Beijing , China
On global optimizations of the rank and inertia of the matrix function A 1 B 1 XB 1 subject to a pair of matrix equations B 2 XB 2, B XB = A 2, A Yongge Tian China Economics and Management Academy, Central
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 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 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 informationResearch Article On the Hermitian R-Conjugate Solution of a System of Matrix Equations
Applied Mathematics Volume 01, Article ID 398085, 14 pages doi:10.1155/01/398085 Research Article On the Hermitian R-Conjugate Solution of a System of Matrix Equations Chang-Zhou Dong, 1 Qing-Wen Wang,
More informationGeneralized Schur complements of matrices and compound matrices
Electronic Journal of Linear Algebra Volume 2 Volume 2 (200 Article 3 200 Generalized Schur complements of matrices and compound matrices Jianzhou Liu Rong Huang Follow this and additional wors at: http://repository.uwyo.edu/ela
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 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 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 informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES
olume 10 2009, Issue 2, Article 41, 10 pp. ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES HANYU LI, HU YANG, AND HUA SHAO COLLEGE OF MATHEMATICS AND PHYSICS CHONGQING UNIERSITY
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 the solvability of an equation involving the Smarandache function and Euler function
Scientia Magna Vol. 008), No., 9-33 On the solvability of an equation involving the Smarandache function and Euler function Weiguo Duan and Yanrong Xue Department of Mathematics, Northwest University,
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 informationELA THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES
Volume 22, pp. 480-489, May 20 THE MINIMUM-NORM LEAST-SQUARES SOLUTION OF A LINEAR SYSTEM AND SYMMETRIC RANK-ONE UPDATES XUZHOU CHEN AND JUN JI Abstract. In this paper, we study the Moore-Penrose inverse
More informationAN ITERATIVE METHOD FOR THE GENERALIZED CENTRO-SYMMETRIC SOLUTION OF A LINEAR MATRIX EQUATION AXB + CY D = E. Ying-chun LI and Zhi-hong LIU
Acta Universitatis Apulensis ISSN: 158-539 No. 9/01 pp. 335-346 AN ITERATIVE METHOD FOR THE GENERALIZED CENTRO-SYMMETRIC SOLUTION OF A LINEAR MATRIX EQUATION AXB + CY D = E Ying-chun LI and Zhi-hong LIU
More informationLecture notes: Applied linear algebra Part 1. Version 2
Lecture notes: Applied linear algebra Part 1. Version 2 Michael Karow Berlin University of Technology karow@math.tu-berlin.de October 2, 2008 1 Notation, basic notions and facts 1.1 Subspaces, range and
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES
ON WEIGHTED PARTIAL ORDERINGS ON THE SET OF RECTANGULAR COMPLEX MATRICES HANYU LI, HU YANG College of Mathematics and Physics Chongqing University Chongqing, 400030, P.R. China EMail: lihy.hy@gmail.com,
More information1 Linear Algebra Problems
Linear Algebra Problems. Let A be the conjugate transpose of the complex matrix A; i.e., A = A t : A is said to be Hermitian if A = A; real symmetric if A is real and A t = A; skew-hermitian if A = A and
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 information76. LIO ND Z. I vector x 2 R n by a quasi-newton-type algorithm. llwright, Woodgate[2,3] and Liao[4] gave some necessary and sufficient conditions for
Journal of Computational Mathematics, Vol.2, No.2, 23, 75 82. LEST-SQURES SOLUTION OF X = D OVER SYMMETRIC POSITIVE SEMIDEFINITE MTRICES X Λ) nping Liao (Department of Mathematics, Hunan Normal University,
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationMatrix Inequalities by Means of Block Matrices 1
Mathematical Inequalities & Applications, Vol. 4, No. 4, 200, pp. 48-490. Matrix Inequalities by Means of Block Matrices Fuzhen Zhang 2 Department of Math, Science and Technology Nova Southeastern University,
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 informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
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 informationOn the Hermitian solutions of the
Journal of Applied Mathematics & Bioinformatics vol.1 no.2 2011 109-129 ISSN: 1792-7625 (print) 1792-8850 (online) International Scientific Press 2011 On the Hermitian solutions of the matrix equation
More informationJournal of Computational and Applied Mathematics. New matrix iterative methods for constraint solutions of the matrix
Journal of Computational and Applied Mathematics 35 (010 76 735 Contents lists aailable at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elseier.com/locate/cam New
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 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 informationOn the Eigenstructure of Hermitian Toeplitz Matrices with Prescribed Eigenpairs
The Eighth International Symposium on Operations Research and Its Applications (ISORA 09) Zhangjiajie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp 298 305 On the Eigenstructure of Hermitian
More informationTHE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR
THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR WEN LI AND MICHAEL K. NG Abstract. In this paper, we study the perturbation bound for the spectral radius of an m th - order n-dimensional
More informationInner products and Norms. Inner product of 2 vectors. Inner product of 2 vectors x and y in R n : x 1 y 1 + x 2 y x n y n in R n
Inner products and Norms Inner product of 2 vectors Inner product of 2 vectors x and y in R n : x 1 y 1 + x 2 y 2 + + x n y n in R n Notation: (x, y) or y T x For complex vectors (x, y) = x 1 ȳ 1 + x 2
More informationSolutions of a constrained Hermitian matrix-valued function optimization problem with applications
Solutions of a constrained Hermitian matrix-valued function optimization problem with applications Yongge Tian CEMA, Central University of Finance and Economics, Beijing 181, China Abstract. Let f(x) =
More informationSome inequalities for unitarily invariant norms of matrices
Wang et al Journal of Inequalities and Applications 011, 011:10 http://wwwjournalofinequalitiesandapplicationscom/content/011/1/10 RESEARCH Open Access Some inequalities for unitarily invariant norms of
More informationOn EP elements, normal elements and partial isometries in rings with involution
Electronic Journal of Linear Algebra Volume 23 Volume 23 (2012 Article 39 2012 On EP elements, normal elements and partial isometries in rings with involution Weixing Chen wxchen5888@163.com Follow this
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
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 informationOn generalized Schur complement of nonstrictly diagonally dominant matrices and general H- matrices
Electronic Journal of Linear Algebra Volume 23 Volume 23 (2012) Article 57 2012 On generalized Schur complement of nonstrictly diagonally dominant matrices and general H- matrices Cheng-Yi Zhang zhangchengyi2004@163.com
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationSingular Value Inequalities for Real and Imaginary Parts of Matrices
Filomat 3:1 16, 63 69 DOI 1.98/FIL16163C Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Singular Value Inequalities for Real Imaginary
More informationThe Reconstruction of An Hermitian Toeplitz Matrices with Prescribed Eigenpairs
The Reconstruction of An Hermitian Toeplitz Matrices with Prescribed Eigenpairs Zhongyun LIU Lu CHEN YuLin ZHANG Received: c 2009 Springer Science + Business Media, LLC Abstract Toeplitz matrices have
More informationarxiv: v1 [math.na] 1 Sep 2018
On the perturbation of an L -orthogonal projection Xuefeng Xu arxiv:18090000v1 [mathna] 1 Sep 018 September 5 018 Abstract The L -orthogonal projection is an important mathematical tool in scientific computing
More information5. Orthogonal matrices
L Vandenberghe EE133A (Spring 2017) 5 Orthogonal matrices matrices with orthonormal columns orthogonal matrices tall matrices with orthonormal columns complex matrices with orthonormal columns 5-1 Orthonormal
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 informationThe eigenvalue distribution of block diagonally dominant matrices and block H-matrices
Electronic Journal of Linear Algebra Volume 20 Volume 20 (2010) Article 45 2010 The eigenvalue distribution of block diagonay dominant matrices and block H-matrices Cheng-Yi Zhang chengyizhang@yahoo.com.cn
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationThe semi-convergence of GSI method for singular saddle point problems
Bull. Math. Soc. Sci. Math. Roumanie Tome 57(05 No., 04, 93 00 The semi-convergence of GSI method for singular saddle point problems by Shu-Xin Miao Abstract Recently, Miao Wang considered the GSI method
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 informationBOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 4, P AGES 655 664 BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION Guang-Da Hu and Qiao Zhu This paper is concerned with bounds of eigenvalues of a complex
More informationResearch Article Some Results on Characterizations of Matrix Partial Orderings
Applied Mathematics, Article ID 408457, 6 pages http://dx.doi.org/10.1155/2014/408457 Research Article Some Results on Characterizations of Matrix Partial Orderings Hongxing Wang and Jin Xu Department
More informationA Note on Simple Nonzero Finite Generalized Singular Values
A Note on Simple Nonzero Finite Generalized Singular Values Wei Ma Zheng-Jian Bai December 21 212 Abstract In this paper we study the sensitivity and second order perturbation expansions of simple nonzero
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 informationNonsingularity and group invertibility of linear combinations of two k-potent matrices
Nonsingularity and group invertibility of linear combinations of two k-potent matrices Julio Benítez a Xiaoji Liu b Tongping Zhu c a Departamento de Matemática Aplicada, Instituto de Matemática Multidisciplinar,
More informationarxiv: v1 [math.ra] 14 Apr 2018
Three it representations of the core-ep inverse Mengmeng Zhou a, Jianlong Chen b,, Tingting Li c, Dingguo Wang d arxiv:180.006v1 [math.ra] 1 Apr 018 a School of Mathematics, Southeast University, Nanjing,
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 informationarxiv: v1 [math.ra] 25 Jul 2013
Representations of the Drazin inverse involving idempotents in a ring Huihui Zhu, Jianlong Chen arxiv:1307.6623v1 [math.ra] 25 Jul 2013 Abstract: We present some formulae for the Drazin inverse of difference
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 4
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 4 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 12, 2012 Andre Tkacenko
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationLecture notes on Quantum Computing. Chapter 1 Mathematical Background
Lecture notes on Quantum Computing Chapter 1 Mathematical Background Vector states of a quantum system with n physical states are represented by unique vectors in C n, the set of n 1 column vectors 1 For
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 informationIr O D = D = ( ) Section 2.6 Example 1. (Bottom of page 119) dim(v ) = dim(l(v, W )) = dim(v ) dim(f ) = dim(v )
Section 3.2 Theorem 3.6. Let A be an m n matrix of rank r. Then r m, r n, and, by means of a finite number of elementary row and column operations, A can be transformed into the matrix ( ) Ir O D = 1 O
More informationarxiv: v1 [math.ra] 24 Aug 2016
Characterizations and representations of core and dual core inverses arxiv:1608.06779v1 [math.ra] 24 Aug 2016 Jianlong Chen [1], Huihui Zhu [1,2], Pedro Patrício [2,3], Yulin Zhang [2,3] Abstract: In this
More informationDistribution for the Standard Eigenvalues of Quaternion Matrices
International Mathematical Forum, Vol. 7, 01, no. 17, 831-838 Distribution for the Standard Eigenvalues of Quaternion Matrices Shahid Qaisar College of Mathematics and Statistics, Chongqing University
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 informationReview of Basic Concepts in Linear Algebra
Review of Basic Concepts in Linear Algebra Grady B Wright Department of Mathematics Boise State University September 7, 2017 Math 565 Linear Algebra Review September 7, 2017 1 / 40 Numerical Linear Algebra
More informationELA ON THE GROUP INVERSE OF LINEAR COMBINATIONS OF TWO GROUP INVERTIBLE MATRICES
ON THE GROUP INVERSE OF LINEAR COMBINATIONS OF TWO GROUP INVERTIBLE MATRICES XIAOJI LIU, LINGLING WU, AND JULIO BENíTEZ Abstract. In this paper, some formulas are found for the group inverse of ap +bq,
More informationHands-on Matrix Algebra Using R
Preface vii 1. R Preliminaries 1 1.1 Matrix Defined, Deeper Understanding Using Software.. 1 1.2 Introduction, Why R?.................... 2 1.3 Obtaining R.......................... 4 1.4 Reference Manuals
More informationECE 275A Homework #3 Solutions
ECE 75A Homework #3 Solutions. Proof of (a). Obviously Ax = 0 y, Ax = 0 for all y. To show sufficiency, note that if y, Ax = 0 for all y, then it must certainly be true for the particular value of y =
More informationOn V-orthogonal projectors associated with a semi-norm
On V-orthogonal projectors associated with a semi-norm Short Title: V-orthogonal projectors Yongge Tian a, Yoshio Takane b a School of Economics, Shanghai University of Finance and Economics, Shanghai
More informationThe Singular Value Decomposition
The Singular Value Decomposition Philippe B. Laval KSU Fall 2015 Philippe B. Laval (KSU) SVD Fall 2015 1 / 13 Review of Key Concepts We review some key definitions and results about matrices that will
More informationLinear Algebra Methods for Data Mining
Linear Algebra Methods for Data Mining Saara Hyvönen, Saara.Hyvonen@cs.helsinki.fi Spring 2007 1. Basic Linear Algebra Linear Algebra Methods for Data Mining, Spring 2007, University of Helsinki Example
More informationJordan Canonical Form of A Partitioned Complex Matrix and Its Application to Real Quaternion Matrices
COMMUNICATIONS IN ALGEBRA, 29(6, 2363-2375(200 Jordan Canonical Form of A Partitioned Complex Matrix and Its Application to Real Quaternion Matrices Fuzhen Zhang Department of Math Science and Technology
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 informationAbstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization relations.
HIROSHIMA S THEOREM AND MATRIX NORM INEQUALITIES MINGHUA LIN AND HENRY WOLKOWICZ Abstract. In this article, several matrix norm inequalities are proved by making use of the Hiroshima 2003 result on majorization
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationCheat Sheet for MATH461
Cheat Sheet for MATH46 Here is the stuff you really need to remember for the exams Linear systems Ax = b Problem: We consider a linear system of m equations for n unknowns x,,x n : For a given matrix A
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 2
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 2 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory April 5, 2012 Andre Tkacenko
More informationPseudoinverse & Orthogonal Projection Operators
Pseudoinverse & Orthogonal Projection Operators ECE 174 Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department, UC San Diego Ken Kreutz-Delgado (UC San Diego) ECE 174 Fall 2016 1 / 48 The Four
More informationOPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY
published in IMA Journal of Numerical Analysis (IMAJNA), Vol. 23, 1-9, 23. OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY SIEGFRIED M. RUMP Abstract. In this note we give lower
More informationLecture 1: Review of linear algebra
Lecture 1: Review of linear algebra Linear functions and linearization Inverse matrix, least-squares and least-norm solutions Subspaces, basis, and dimension Change of basis and similarity transformations
More informationFirst, we review some important facts on the location of eigenvalues of matrices.
BLOCK NORMAL MATRICES AND GERSHGORIN-TYPE DISCS JAKUB KIERZKOWSKI AND ALICJA SMOKTUNOWICZ Abstract The block analogues of the theorems on inclusion regions for the eigenvalues of normal matrices are given
More information