c 2005 Society for Industrial and Applied Mathematics
|
|
- Diane Horton
- 5 years ago
- Views:
Transcription
1 SIAM J. MATRIX ANAL. APPL. Vol. XX, No. X, pp. XX XX c 005 Society for Industrial and Applied Mathematics DISTRIBUTIONS OF THE EXTREME EIGENVALUES OF THE COMPLEX JACOBI RANDOM MATRIX ENSEMBLE PLAMEN KOEV AND IOANA DUMITRIU Abstract. We present explicit formulas for the distributions of the extreme eigenvalues of the complex Jacobi ensemble in terms of the hypergeometric function of a matrix argument. Our methods extend to the quaternion Jacobi ensemble case. Key words. random matrix, Jacobi distribution, MANOVA, multivariate beta, hypergeometric function of a matrix argument, eigenvalue AMS subject classifications. 5A5, 60E05, 6H0, 65F5 DOI. XX.XXXX/SXXXXXXXXXXXXXXXX. Introduction. Complex random covariance matrices arise naturally in applications involving wireless communications [, 5, 8]. Various statistical methods, such as hypothesis testing, principal component analysis, canonical correlation analysis, multivariate analysis of variance, etc., are then based on the distributions of the extreme (generalized) eigenvalues of these matrices [3]. In this paper we present explicit formulas for the distributions of the extreme eigenvalues of the complex Jacobi ensemble in terms of the hypergeometric function of a matrix argument. The formulas are natural analogues of the formulas in the real case. Recently, Ratnarajah, Vaillancourt, and Alvo produced formulas for the distributions of the extreme eigenvalues of the complex Wishart ensemble [4]. Their results, too, are the natural analogues of the real case discussed in Muirhead [3]. Such explicit formulas become increasingly important now that there are algorithms for their efficient computation in practice [9, 0]. Our methods easily extend to the less well-known quaternion Jacobi distribution. This paper is organized as follows. In section we introduce the formal definitions of the Wishart and Jacobi ensembles and the hypergeometric function of a matrix argument that we need later on. In section 3 we survey the existing formulas for the distributions of the extreme eigenvalues of the real Jacobi ensemble. We present our main result formulas for the distributions of the extreme eigenvalues of the complex Jacobi ensemble in section 5. In section 6, we introduce the quaternion Jacobi ensemble, and present formulas for the distributions of its extreme eigenvalues. Finally, we present numerical experiments in section 7.. Background. In this section we recall the definitions of the Wishart and Jacobi random matrix ensembles, and of the hypergeometric function of a matrix argument. Received by the editors October, 005; accepted for publication (in revised form) by EDITOR DATE; published electronically DATE. This work was supported by National Science Foundation grant DMS-03486, and by the Miller Institute for Research in Basic Sciences, U.C. Berkeley, CA. Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 039 (plamen@math.mit.edu). Department of Mathematics, University of California, Berkeley, CA 9470 (dumitriu@math.berkeley.edu).
2 PLAMEN KOEV and IOANA DUMITRIU We begin with the real and complex Wishart distributions, since they are the building blocks for the respective Jacobi distributions. If the n m real Gaussian random matrix Z is distributed as N(0, I n Σ), then the matrix A Z T Z is called m m real central Wishart matrix with n degrees of freedom and covariance matrix Σ [3, Def. 3..3, p. 8]. Similarly, if the n m complex Gaussian random matrix Z is distributed as N(0, I n Σ), then the matrix A Z H Z is called m m complex central Wishart matrix with n degrees of freedom and covariance matrix Σ [4]. The real and complex central Wishart ensembles are denoted as W m (n, Σ) and CW m (n, Σ), respectively. If A W m (n, Σ) and B W m (n, Σ), then C A(A + B) is said to have (real) Jacobi distribution. Analogously, if A CW m (n, Σ) and B CW m (n, Σ), then C A(A + B) is said to have complex Jacobi distribution. The Jacobi distribution is sometimes called multivariate beta distribution [3] and is closely related to the MANOVA (Multivariate ANalysis Of VAriance) distribution [3] (the difference is that MANOVA studies the matrix AB instead of A(A + B) ). Definition. (Hypergeometric function of a matrix argument). Let p 0 and q 0 be integers, and let X be an m m complex symmetric matrix with eigenvalues x, x,..., x m. The hypergeometric function of a matrix argument X and parameter α > 0 is defined as (.) pf q (α) (a,..., a p ; b,..., b q ; X) where k means is a partition of k, (.) (a) (α) (i,j) k0 k (a ) (α) (a p ) (α) k!(b ) (α) ( a i ) α + j (b q ) (α) C (α) (X), is the generalized Pochhammer symbol, and C (α) (X) is the Jack function [7]. The Jack function C (α) (X) is a symmetric, homogeneous polynomial in the eigenvalues x, x,..., x m of X. We refer the reader to the excellent survey of Stanley [5] for the explicit definition of C (α) (X) as we do not need it here. In the complex case, α, C (α) (X) equals the well-known Schur function (appropriately scaled), in the real case, α, the zonal polynomial, and in the quaternion case, α, the quaternion zonal polynomial ([]). In this paper we are primarily concerned with the complex case, α ; in Section 6 we deal with the case α. We also utilize the multivariate Gamma function of parameter α, defined as m ( Γ m (α) (c) π m(m ) α Γ c i ) for R(c) > m α α. i 3. Survey of the real case. Let A W m (n, Σ), B W m (n, Σ), where n m and n m. Then the matrix C A(A + B) has the Jacobi real distribution. For the distribution of its largest eigenvalue, we have from [3, eq. (6)] (3.) P (λ max (C) < x) Γ() m ( n+n )Γ () m ( m+ Γ () m ( n+m+ )Γ () ) m ( n ) x mn F () ( n, n+m+ ; n+m+ ; xi).
3 EIGENVALUES OF COMPLEX JACOBI ENSEMBLE 3 The distribution of the smallest eigenvalue is closely related to the distribution of the largest one: (3.) P (λ min (C) < x) P (λ min (C) > x) P (A(A + B) > xi) P (A > x(a + B)) P (B(A + B) < ( x)i). The latest expression is now easily evaluated using (3.) to obtain (3.3) P (λ min (C) < x) Γ() m ( n+n )Γ () m ( m+ ) mn Γ () m ( n+m+ x) )Γ () m ( n )( F () ( n, m n+ ; n+m+ ; ( x)i). Interestingly, we have an explicit formula for the density of the smallest eigenvalue of C [, eq. (9)] dens(λ min (C)) mn Γ() m ( n+n )Γ () m ( m+ Γ () m ( n+m+ ) ) x n m ( x) mn )Γ () m ( n F () ( n, m n+ ; n+m+ ; ( x)i m ). 4. Integral formulas for F () and F () functions. Theorem 7.4. in Muir- functions (the real case). head [3] gives integral representations of the F () and F () Here we present their complex versions. Theorem 4.. The F () function has the integral representation (4.) F () (a; c; X) m (c) m (a) m (c a) e tr (XY ) (det Y ) a m det(i Y ) c a m (dy ), valid for R(a) > m, R(c) > m, R(c a) > m. The F () function has the integral representation (4.) F () (a, b; c; X) m (c) m (a) m (c a) det(i XY ) b (det Y ) a m det(i Y ) c a m (dy ), valid for X <, R(a) > m, R(c a) > m. The expression B < C, where B and C are matrices, means that the matrix C B is Hermitian positive definite. Proof. We proceed analogously to the proof of Theorem 7.4. in Muirhead [3, p. 64]. From Proposition.5 in [4] (4.3) m (det Y ) a m det(i m Y ) b m C () (XY )(dy ) when R(a) > m and R(b) > m. (a)() (a + b) () Γ() m (a) m (b) m (a + b) C () (X),
4 4 PLAMEN KOEV and IOANA DUMITRIU To prove (4.) we expand e tr (XY ) k0 k C () (XY ) k! [6, eq. (89), p. 488] and integrate term-by-term using (4.3). To prove (4.) we expand det(i XY ) b (b) () C () (XY ) k! k0 k [6, eq. (90), p. 488] and integrate the integral in (4.) term-by-term using (4.3). 5. The extreme eigenvalues of the complex Jacobi ensemble. Let the matrix C A(A+B) have complex Jacobi distribution, where A CW m (n, Σ), n m and B CW m (n, Σ), n m. We now derive formulas for the distributions of λ max (C) and λ min (C). Theorem 5.. The distribution of λ max (C) is (5.) P (λ max (C) < x) Γ() m (n + n ) m (m) m (n + m) m (n ) xmn F () (n, m n ; n + m; xi). Proof. The density of A is [6, section 8] The joint density of A and B is (det Σ) n m (n ) (det A)n m e tr ( AΣ ) (da). (det Σ) (n+n) m (n ) m (n ) (det A)n m (det B) n m e tr ( (A+B)Σ ) (da)(db). Let U L AL T, where L is the Cholesky factor of A + B. The matrices U and C thus have the same eigenvalues. By repeating, verbatim, the argument of Theorem 3.3. in Muirhead [3, p. 09], we conclude that the density function of U is (5.) m (n + n ) m (n ) m (n ) (det U)n m det(i U) n m. To compute the distribution of the largest eigenvalue of U we integrate the density (5.) from 0 to xi (5.3) P (U < xi) Γ() m (n + n ) (det U) n m det(i U) n m (du). m (n ) m (n ) 0<U<xI We make an m -dimensional change of variables xy U in (5.3) to obtain (5.4) P (U < xi) Γ() m (n + n )x mn (det Y ) n m det(i xy ) n m (dy ). m (n ) m (n ) We use (4.) to evaluate the integral in (5.4) and obtain (5.). Corollary 5.. The distribution of λ min (C) is (5.5) P (λ min (C) < x) Γ() m (n + n ) m (m) m (n + m) m (n ) ( x)mn F () (n, m n ; n + m; ( x)i). Proof. Directly from (3.) and (5.).
5 EIGENVALUES OF COMPLEX JACOBI ENSEMBLE 5 6. The quaternion case. The quaternion Wishart and Jacobi ensembles are somewhat less known in the world of multivariate statistics and its applications than the real and complex ones, just as the Gaussian Symplectic ensemble is somewhat less known in the quantum physics and mechanical statistics world than the Orthogonal and Unitary ones. One reason for this discrimination is the fact that quaternion normal variables are harder to sample from in existing software. The triad real-complex-quaternion (or equivalently, using the orthogonal symmetry invariance, orthogonal-unitary-symplectic) is what defines the threefold way in random matrix theory (see Dyson, [4]): the three division algebras over the real numbers, over which the concept of a multivariate normal variable is well-defined, and random matrix models of i.i.d. standard Gaussians can be constructed. A quaternion standard Gaussian variable is defined as q (n +n i+n 3 j+n 4 k), where n, n, n 3, n 4 are standard real normals, and i, j, k are the quaternion roots of, i.e., i j k ijk. The quaternion conjugate for q a + a i + a 3 j + a 4 k is q a a i a 3 j a 4 k; qq q a + a + a 3 + a 4. If the n m quaternion Gaussian matrix Y is distributed as HN(0, I n Σ), the matrix A Y D Y has the quaternion central Wishart distribution with n degrees of freedom and covariance matrix Σ. Here Y D denotes the dual (or quaternion conjugate transpose) matrix of Y ; if Y D Y, the matrix Y is called self-dual (and has real eigenvalues). We denote by HW m (n, Σ) the quaternion central Wishart ensemble. Similarly, if A HW m (n, Σ) and B HW m (n, Σ), the matrix C A(A+B) has quaternion Jacobi distribution. Our methods of calculating the extreme eigenvalue distributions for the complex Jacobi ensembles extend immediately to the quaternion case; we hence state the results without repeating the proofs. Theorem 6.. The F ( ) function has the integral representation (6.) F ( ) (a; c; X) Γ ( ) m (c) Γ ( ) m (a)γ ( ) m (c a) e tr (XY ) (det Y ) a m+ det(i Y ) c a m+ (dy ), valid for R(a) > (m ), R(c) > (m ), R(c a) > (m ). The F ( ) function has the integral representation (6.) F ( ) (a, b; c; X) Γ ( ) m (c) Γ ( ) m (a)γ ( ) m (c a) det(i XY ) b (det Y ) a m+ det(i Y ) c a m+ (dy ), valid for X <, R(a) > (m ), R(c a) > (m ). The expression B < C, where B and C are matrices, means that the matrix C B is self-dual positive definite. Theorem 6.. Let C have quaternion Jacobi distribution, C A(A + B), where A HW m (n, Σ) and B HW m (n, Σ). Then the distributions of λ max (C) The notation H stands for William Hamilton, who introduced the quaternions in the 850s.
6 6 PLAMEN KOEV and IOANA DUMITRIU Complex Jacobi, (n,n,m)(3,5,) Complex Jacobi, (n,n,m)(8,0,7) Fig. 7.. The distribution of the largest eigenvalue of the complex Jacobi ensemble. The solid line is the empirical distribution resulting from a Monte Carlo experiment with 0, 000 replications. The circle represents the formula (5.). and λ min (C) are (6.3) (6.4) P (λ max (C) < x) Γ( ) m (n + n )Γ ( ) m (m ) x mn Γ ( ) m (n + m )Γ ( ) m (n ) F ( ) (n, m n ; n + m ; xi); P (λ min (C) < x) Γ( ) m (n + n )Γ ( ) m (m ) ( x) mn Γ ( ) m (n + m )Γ ( ) m (n ) F ( ) (n, m n ; n + m ; ( x)i). 7. Numerical Experiments. We performed extensive numerical tests to confirm the correctness of Theorems 5. and 6. against Monte Carlo experiments; we used MATLAB [] to compute and plot our distributions. We detail several of those tests here whose results were typical. In the complex case we choose two sets of parameters (n, n, m) and for each set compared the empirical distribution from 0, 000 replications with the analytical formula (5.), computed using the algorithms from [9, 0]. In both cases, the results were a perfect visual match see Figure 7.. We generated our matrices in MATLAB as CA/(A+B), where A CW m (n, I), B CW m (n, I). In turn A was generated as AZ *Z/, where Zrandn(n, m)+i*randn(n, m), and B was generated analogously. In the quaternion case, we side-stepped the need to simulate quaternion variables. Instead, we tested with the very elegant (real, tridiagonal) β-jacobi matrix constructed by Sutton [6, Chapter 5]. For β 4 it has the same eigenvalue distribution as the quaternion Jacobi matrix C A(A + B), A HW m (n, Σ), B HW m (n, Σ). Note that we can assume that Σ I without any loss of generality [3, Thm ].
7 EIGENVALUES OF COMPLEX JACOBI ENSEMBLE 7 Quaternion Jacobi, (n,n,m)(3,5,) Quaternion Jacobi, (n,n,m)(8,0,7) Fig. 7.. The distribution of the largest eigenvalue of the quaternion Jacobi ensemble. The solid line is the empirical distribution resulting from a Monte Carlo experiment with 0, 000 replications. The circle represents the formula (6.3). The β-jacobi matrix is defined as Q Z T Z, where c n s n c n Z c n s... n... s c, with c s c k Beta ( β (n m + k), β (n m + k) ), s k c k ; c k Beta ( β k, β (n + n m + k + ) ), s k c k. We plotted the empirical distribution of the eigenvalues of Q against the analytical formula (6.3) in Figure 7., for the same triplets (n, n, m) we used in the complex Jacobi case. Once again we obtained a perfect visual match. REFERENCES [] P.-A. Absil, A. Edelman, and P. Koev, On the largest principal angle between random subspaces. Linear Algebra Appl., to appear. [] P. Chen, G. Genello, and M. Wicks, Estimating the number of signals in presence of colored noise, in Radar Conference 004. Proceedings of the IEEE, 6 9 April 004, pp [3] A. G. Constantine, Some non-central distribution problems in multivariate analysis, Ann. Math. Statist., 34 (963), pp [4] F. J. Dyson, The threefold way. Algebraic structures of symmetry groups and ensembles in Quantum Mechanics, J. Math. Phys., 3 (963), pp [5] H. Gao, P. Smith, and M. Clark, Theoretical reliability of MMSE linear diversity combining in Rayleigh-fading additive interference channels, IEEE Transactions on Communications, 46 (998), pp [6] A. T. James, Distributions of matrix variates and latent roots derived from normal samples, Ann. Math. Statist., 35 (964), pp [7] J. Kaneko, Selberg integrals and hypergeometric functions associated with Jack polynomials, SIAM J. Math. Anal., 4 (993), pp [8] M. Kang and M.-S. Alouini, Largest eigenvalue of complex Wishart matrices and performance analysis of MIMO MRC systems, IEEE Journal on Selected Areas in Communications, (003), pp
8 8 PLAMEN KOEV and IOANA DUMITRIU [9] P. Koev, Software for computing the hypergeometric function of a matrix argument. Available from [0] P. Koev and A. Edelman, The efficient evaluation of the hypergeometric function of a matrix argument. Math. Comp., to appear. [] I. G. Macdonald, Symmetric functions and Hall polynomials, Oxford University Press, New York, Second ed., 995. [] The MathWorks, Inc., MATLAB Reference Guide, Natick, MA, 99. [3] R. J. Muirhead, Aspects of multivariate statistical theory, John Wiley & Sons Inc., New York, 98. [4] T. Ratnarajah, R. Vaillancourt, and M. Alvo, Eigenvalues and condition numbers of complex random matrices, SIAM J. Matrix Anal. Appl., 6 (004/05), pp [5] R. Stanley, Some combinatorial properties of Jack symmetric functions, Adv. Math., 77 (989), pp [6] B. Sutton, The stochastic operator approach to random matrix theory, PhD thesis, Massachusetts Institute of Technology, 005.
Progress in the method of Ghosts and Shadows for Beta Ensembles
Progress in the method of Ghosts and Shadows for Beta Ensembles Alan Edelman (MIT) Alex Dubbs (MIT) and Plamen Koev (SJS) Oct 8, 2012 1/47 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians
More informationTHE EFFICIENT EVALUATION OF THE HYPERGEOMETRIC FUNCTION OF A MATRIX ARGUMENT
MATHEMATICS OF COMPUTATION Volume 75, Number 54, Pages 833 846 S 5-5718(6)184- Article electronically published on January 19, 6 THE EFFICIENT EVALUATION OF THE HYPERGEOMETRIC FUNCTION OF A MATRIX ARGUMENT
More informationMIT Term Project: Numerical Experiments on Circular Ensembles and Jack Polynomials with Julia
MIT 18.338 Term Project: Numerical Experiments on Circular Ensembles and Jack Polynomials with Julia N. Kemal Ure May 15, 2013 Abstract This project studies the numerical experiments related to averages
More information1. Introduction. Let an n m complex Gaussian random matrix A be distributed
COMMUNICATIONS IN INFORMATION AND SYSTEMS c 23 International Press Vol 3, No 2, pp 119-138, October 23 3 COMPLEX RANDOM MATRICES AND RAYLEIGH CHANNEL CAPACITY T RATNARAJAH, R VAILLANCOURT, AND M ALVO Abstract
More informationAdvances in Random Matrix Theory: Let there be tools
Advances in Random Matrix Theory: Let there be tools Alan Edelman Brian Sutton, Plamen Koev, Ioana Dumitriu, Raj Rao and others MIT: Dept of Mathematics, Computer Science AI Laboratories World Congress,
More informationThe Random Matrix Technique of Ghosts and Shadows
The Random Matrix Technique of Ghosts and Shadows Alan Edelman November 22, 2009 Abstract We propose to abandon the notion that a random matrix has to be sampled for it to exist. Much of today's applied
More informationOn the algebraic structure of Dyson s Circular Ensembles. Jonathan D.H. SMITH. January 29, 2012
Scientiae Mathematicae Japonicae 101 On the algebraic structure of Dyson s Circular Ensembles Jonathan D.H. SMITH January 29, 2012 Abstract. Dyson s Circular Orthogonal, Unitary, and Symplectic Ensembles
More informationComplex singular Wishart matrices and applications
Complex singular Wishart matrices and applications T Ratnarajah R Vaillancourt CRM-3185 May 2005 Thans are due to the Natural Sciences and Engineering Council of Canada and the Centre de recherches mathématiques
More informationRandom Matrices and Multivariate Statistical Analysis
Random Matrices and Multivariate Statistical Analysis Iain Johnstone, Statistics, Stanford imj@stanford.edu SEA 06@MIT p.1 Agenda Classical multivariate techniques Principal Component Analysis Canonical
More informationRandom Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws. Symeon Chatzinotas February 11, 2013 Luxembourg
Random Matrix Theory Lecture 1 Introduction, Ensembles and Basic Laws Symeon Chatzinotas February 11, 2013 Luxembourg Outline 1. Random Matrix Theory 1. Definition 2. Applications 3. Asymptotics 2. Ensembles
More information1 Intro to RMT (Gene)
M705 Spring 2013 Summary for Week 2 1 Intro to RMT (Gene) (Also see the Anderson - Guionnet - Zeitouni book, pp.6-11(?) ) We start with two independent families of R.V.s, {Z i,j } 1 i
More informationLog-Convexity Properties of Schur Functions and Generalized Hypergeometric Functions of Matrix Argument. Donald St. P. Richards.
Log-Convexity Properties of Schur Functions and Generalized Hypergeometric Functions of Matrix Argument Donald St. P. Richards August 22, 2009 Abstract We establish a positivity property for the difference
More informationRectangular Young tableaux and the Jacobi ensemble
Rectangular Young tableaux and the Jacobi ensemble Philippe Marchal October 20, 2015 Abstract It has been shown by Pittel and Romik that the random surface associated with a large rectangular Young tableau
More informationZonal Polynomials and Hypergeometric Functions of Matrix Argument. Zonal Polynomials and Hypergeometric Functions of Matrix Argument p.
Zonal Polynomials and Hypergeometric Functions of Matrix Argument Zonal Polynomials and Hypergeometric Functions of Matrix Argument p. 1/2 Zonal Polynomials and Hypergeometric Functions of Matrix Argument
More informationOn Computing Schur Functions and Series Thereof
J Algebr Comb manuscript No. (will be inserted by the editor) On Computing Schur Functions and Series Thereof Cy Chan Vesselin Drensky Alan Edelman Raymond Kan Plamen Koev Received: date / Accepted: date
More informationEigenvalue PDFs. Peter Forrester, M&S, University of Melbourne
Outline Eigenvalue PDFs Peter Forrester, M&S, University of Melbourne Hermitian matrices with real, complex or real quaternion elements Circular ensembles and classical groups Products of random matrices
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 1982 NOTES ON MATRIX METHODS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Chemistry 5.76 Revised February, 198 NOTES ON MATRIX METHODS 1. Matrix Algebra Margenau and Murphy, The Mathematics of Physics and Chemistry, Chapter 10, give almost
More informationA new type of PT-symmetric random matrix ensembles
A new type of PT-symmetric random matrix ensembles Eva-Maria Graefe Department of Mathematics, Imperial College London, UK joint work with Steve Mudute-Ndumbe and Matthew Taylor Department of Mathematics,
More informationEigenvalue Statistics for Beta-Ensembles. Ioana Dumitriu
Eigenvalue Statistics for Beta-Ensembles by Ioana Dumitriu B.A., New York University, 1999 Submitted to the Department of Mathematics in partial fulfillment of the requirements for the degree of DOCTOR
More informationPROMPTED by recent results suggesting possible extraordinary
1102 IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL. 4, NO. 3, MAY 2005 On the Capacity of Multiple-Antenna Systems in Rician Fading Sudharman K. Jayaweera, Member, IEEE, and H. Vincent Poor, Fellow,
More informationDOUBLY NONCENTRAL MATRIX VARIATE BETA DISTRIBUTION
DOUBLY NONCENTRAL MATRIX VARIATE BETA DISTRIBUTION José Antonio Díaz García & Ramón Gutiérrez-Jáimez Comunicación Técnica No I-06-08/28-03-2006 PE/CIMAT) Doubly noncentral matrix variate beta distribution
More informationSEA s workshop- MIT - July 10-14
Matrix-valued Stochastic Processes- Eigenvalues Processes and Free Probability SEA s workshop- MIT - July 10-14 July 13, 2006 Outline Matrix-valued stochastic processes. 1- definition and examples. 2-
More informationMULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING. Kaitlyn Beaudet and Douglas Cochran
MULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING Kaitlyn Beaudet and Douglas Cochran School of Electrical, Computer and Energy Engineering Arizona State University, Tempe AZ 85287-576 USA ABSTRACT The problem
More informationFree Probability, Sample Covariance Matrices and Stochastic Eigen-Inference
Free Probability, Sample Covariance Matrices and Stochastic Eigen-Inference Alan Edelman Department of Mathematics, Computer Science and AI Laboratories. E-mail: edelman@math.mit.edu N. Raj Rao Deparment
More informationLimits for BC Jacobi polynomials
Limits for Korteweg-de Vries Institute, University of Amsterdam T.H.Koornwinder@uva.nl http://www.science.uva.nl/~thk/ Lecture on September 10, 2012 at the Conference on Harmonic Analysis, Convolution
More informationResearch Article On the Marginal Distribution of the Diagonal Blocks in a Blocked Wishart Random Matrix
Hindawi Publishing Corporation nternational Journal of Analysis Volume 6, Article D 59678, 5 pages http://dx.doi.org/.55/6/59678 Research Article On the Marginal Distribution of the Diagonal Blocks in
More information672 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 2, FEBRUARY We only include here some relevant references that focus on the complex
672 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 2, FEBRUARY 2009 Ordered Eigenvalues of a General Class of Hermitian Rom Matrices With Application to the Performance Analysis of MIMO Systems Luis
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
More informationAn Efficient Approach to Multivariate Nakagami-m Distribution Using Green s Matrix Approximation
IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS, VOL 2, NO 5, SEPTEMBER 2003 883 An Efficient Approach to Multivariate Nakagami-m Distribution Using Green s Matrix Approximation George K Karagiannidis, Member,
More informationIntroduction to Group Theory
Chapter 10 Introduction to Group Theory Since symmetries described by groups play such an important role in modern physics, we will take a little time to introduce the basic structure (as seen by a physicist)
More informationMeixner matrix ensembles
Meixner matrix ensembles W lodek Bryc 1 Cincinnati April 12, 2011 1 Based on joint work with Gerard Letac W lodek Bryc (Cincinnati) Meixner matrix ensembles April 12, 2011 1 / 29 Outline of talk Random
More informationLecture 10: A (Brief) Introduction to Group Theory (See Chapter 3.13 in Boas, 3rd Edition)
Lecture 0: A (Brief) Introduction to Group heory (See Chapter 3.3 in Boas, 3rd Edition) Having gained some new experience with matrices, which provide us with representations of groups, and because symmetries
More informationMeasures and Jacobians of Singular Random Matrices. José A. Díaz-Garcia. Comunicación de CIMAT No. I-07-12/ (PE/CIMAT)
Measures and Jacobians of Singular Random Matrices José A. Díaz-Garcia Comunicación de CIMAT No. I-07-12/21.08.2007 (PE/CIMAT) Measures and Jacobians of singular random matrices José A. Díaz-García Universidad
More informationMOPS: Multivariate orthogonal polynomials (symbolically)
Journal of Symbolic Computation 42 2007) 587 620 www.elsevier.com/locate/jsc MOPS: Multivariate orthogonal polynomials symbolically) Ioana Dumitriu a,, Alan Edelman b, Gene Shuman c a Department of Mathematics,
More informationMultivariate Linear Models
Multivariate Linear Models Stanley Sawyer Washington University November 7, 2001 1. Introduction. Suppose that we have n observations, each of which has d components. For example, we may have d measurements
More informationComplex Random Matrices and Rayleigh Channel Capacity
Complex Random Matrices and Rayleigh Channel Capacity T Ratnarajah R Vaillancourt M Alvo CRM-293 June 23 This wor was partially supported by the Natural Sciences and Engineering Research Council of Canada
More informationRandom Matrices and Wireless Communications
Random Matrices and Wireless Communications Jamie Evans Centre for Ultra-Broadband Information Networks (CUBIN) Department of Electrical and Electronic Engineering University of Melbourne 3.5 1 3 0.8 2.5
More informationQuadratic Forms on Complex Random Matrices and Multi-Antenna Channel Capacity
Quadratic Forms on Complex Random Matrices and Multi-Antenna Channel Capacity T. Ratnarajah R. Vaillancourt CRM-2979 March 24 Department of Mathematics and Statistics, University of Ottawa, 585 King Edward
More informationRepresentation of Lie Groups and Special Functions
Representation of Lie Groups and Special Functions Recent Advances by N. Ja. Vilenkint formerly of The Correspondence Pedagogical Institute, Moscow, Russia and A.U. Klimyk Institute for Theoretical Physics,
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 informationPrimes, partitions and permutations. Paul-Olivier Dehaye ETH Zürich, October 31 st
Primes, Paul-Olivier Dehaye pdehaye@math.ethz.ch ETH Zürich, October 31 st Outline Review of Bump & Gamburd s method A theorem of Moments of derivatives of characteristic polynomials Hypergeometric functions
More informationMATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants.
MATH 323 Linear Algebra Lecture 6: Matrix algebra (continued). Determinants. Elementary matrices Theorem 1 Any elementary row operation σ on matrices with n rows can be simulated as left multiplication
More informationErgodic and Outage Capacity of Narrowband MIMO Gaussian Channels
Ergodic and Outage Capacity of Narrowband MIMO Gaussian Channels Yang Wen Liang Department of Electrical and Computer Engineering The University of British Columbia April 19th, 005 Outline of Presentation
More informationMIMO Capacities : Eigenvalue Computation through Representation Theory
MIMO Capacities : Eigenvalue Computation through Representation Theory Jayanta Kumar Pal, Donald Richards SAMSI Multivariate distributions working group Outline 1 Introduction 2 MIMO working model 3 Eigenvalue
More informationJournal of Multivariate Analysis. Sphericity test in a GMANOVA MANOVA model with normal error
Journal of Multivariate Analysis 00 (009) 305 3 Contents lists available at ScienceDirect Journal of Multivariate Analysis journal homepage: www.elsevier.com/locate/jmva Sphericity test in a GMANOVA MANOVA
More informationUnifying Analysis of Ergodic MIMO Capacity in Correlated Rayleigh Fading Environments
Unifying Analysis of Ergodic MIMO Capacity in Correlated Rayleigh Fading Environments Mario Kiessling,, Joachim Speidel, Markus Reinhar Institute of elecommunications, University of Stuttgart, Germany
More informationGaussian Free Field in beta ensembles and random surfaces. Alexei Borodin
Gaussian Free Field in beta ensembles and random surfaces Alexei Borodin Corners of random matrices and GFF spectra height function liquid region Theorem As Unscaled fluctuations Gaussian (massless) Free
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 informationLecture Notes 1: Vector spaces
Optimization-based data analysis Fall 2017 Lecture Notes 1: Vector spaces In this chapter we review certain basic concepts of linear algebra, highlighting their application to signal processing. 1 Vector
More informationTight Lower Bounds on the Ergodic Capacity of Rayleigh Fading MIMO Channels
Tight Lower Bounds on the Ergodic Capacity of Rayleigh Fading MIMO Channels Özgür Oyman ), Rohit U. Nabar ), Helmut Bölcskei 2), and Arogyaswami J. Paulraj ) ) Information Systems Laboratory, Stanford
More informationELEC E7210: Communication Theory. Lecture 10: MIMO systems
ELEC E7210: Communication Theory Lecture 10: MIMO systems Matrix Definitions, Operations, and Properties (1) NxM matrix a rectangular array of elements a A. an 11 1....... a a 1M. NM B D C E ermitian transpose
More informationTriangular matrices and biorthogonal ensembles
/26 Triangular matrices and biorthogonal ensembles Dimitris Cheliotis Department of Mathematics University of Athens UK Easter Probability Meeting April 8, 206 2/26 Special densities on R n Example. n
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 informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationA Note on Some Wishart Expectations
Metrika, Volume 33, 1986, page 247-251 A Note on Some Wishart Expectations By R. J. Muirhead 2 Summary: In a recent paper Sharma and Krishnamoorthy (1984) used a complicated decisiontheoretic argument
More informationParallel Additive Gaussian Channels
Parallel Additive Gaussian Channels Let us assume that we have N parallel one-dimensional channels disturbed by noise sources with variances σ 2,,σ 2 N. N 0,σ 2 x x N N 0,σ 2 N y y N Energy Constraint:
More informationThe Maximum Dimension of a Subspace of Nilpotent Matrices of Index 2
The Maximum Dimension of a Subspace of Nilpotent Matrices of Index L G Sweet Department of Mathematics and Statistics, University of Prince Edward Island Charlottetown, PEI C1A4P, Canada sweet@upeica J
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationComplex Random Matrices and Rician Channel Capacity
Complex Random Matrices and Rician Channel Capacity T. Ratnarajah R. Vaillancourt M. Alvo CRM-3174 October 24 This wor was partially supported by the Natural Sciences and Engineering Research Council of
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 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 informationTesting a Normal Covariance Matrix for Small Samples with Monotone Missing Data
Applied Mathematical Sciences, Vol 3, 009, no 54, 695-70 Testing a Normal Covariance Matrix for Small Samples with Monotone Missing Data Evelina Veleva Rousse University A Kanchev Department of Numerical
More informationOn the reduction of matrix polynomials to Hessenberg form
Electronic Journal of Linear Algebra Volume 3 Volume 3: (26) Article 24 26 On the reduction of matrix polynomials to Hessenberg form Thomas R. Cameron Washington State University, tcameron@math.wsu.edu
More informationA Formula for the Specialization of Skew Schur Functions
A Formula for the Specialization of Skew Schur Functions The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published Publisher
More informationMINIMAL NORMAL AND COMMUTING COMPLETIONS
INTERNATIONAL JOURNAL OF INFORMATION AND SYSTEMS SCIENCES Volume 4, Number 1, Pages 5 59 c 8 Institute for Scientific Computing and Information MINIMAL NORMAL AND COMMUTING COMPLETIONS DAVID P KIMSEY AND
More informationDiagonalizing Matrices
Diagonalizing Matrices Massoud Malek A A Let A = A k be an n n non-singular matrix and let B = A = [B, B,, B k,, B n ] Then A n A B = A A 0 0 A k [B, B,, B k,, B n ] = 0 0 = I n 0 A n Notice that A i B
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 informationGaussian Free Field in (self-adjoint) random matrices and random surfaces. Alexei Borodin
Gaussian Free Field in (self-adjoint) random matrices and random surfaces Alexei Borodin Corners of random matrices and GFF spectra height function liquid region Theorem As Unscaled fluctuations Gaussian
More informationErgodic and Outage Capacity of Narrowband MIMO Gaussian Channels
Ergodic and Outage Capacity of Narrowband MIMO Gaussian Channels Yang Wen Liang Department of Electrical and Computer Engineering The University of British Columbia, Vancouver, British Columbia Email:
More information1 Tridiagonal matrices
Lecture Notes: β-ensembles Bálint Virág Notes with Diane Holcomb 1 Tridiagonal matrices Definition 1. Suppose you have a symmetric matrix A, we can define its spectral measure (at the first coordinate
More informationCLASSICAL GROUPS DAVID VOGAN
CLASSICAL GROUPS DAVID VOGAN 1. Orthogonal groups These notes are about classical groups. That term is used in various ways by various people; I ll try to say a little about that as I go along. Basically
More informationMULTIVARIATE ANALYSIS OF VARIANCE UNDER MULTIPLICITY José A. Díaz-García. Comunicación Técnica No I-07-13/ (PE/CIMAT)
MULTIVARIATE ANALYSIS OF VARIANCE UNDER MULTIPLICITY José A. Díaz-García Comunicación Técnica No I-07-13/11-09-2007 (PE/CIMAT) Multivariate analysis of variance under multiplicity José A. Díaz-García Universidad
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 informationStat 159/259: Linear Algebra Notes
Stat 159/259: Linear Algebra Notes Jarrod Millman November 16, 2015 Abstract These notes assume you ve taken a semester of undergraduate linear algebra. In particular, I assume you are familiar with the
More informationFluctuations of Random Matrices and Second Order Freeness
Fluctuations of Random Matrices and Second Order Freeness james mingo with b. collins p. śniady r. speicher SEA 06 Workshop Massachusetts Institute of Technology July 9-14, 2006 1 0.4 0.2 0-2 -1 0 1 2-2
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 informationApplications and fundamental results on random Vandermon
Applications and fundamental results on random Vandermonde matrices May 2008 Some important concepts from classical probability Random variables are functions (i.e. they commute w.r.t. multiplication)
More informationUniversality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium
Universality of distribution functions in random matrix theory Arno Kuijlaars Katholieke Universiteit Leuven, Belgium SEA 06@MIT, Workshop on Stochastic Eigen-Analysis and its Applications, MIT, Cambridge,
More informationHigher rank numerical ranges of rectangular matrix polynomials
Journal of Linear and Topological Algebra Vol. 03, No. 03, 2014, 173-184 Higher rank numerical ranges of rectangular matrix polynomials Gh. Aghamollaei a, M. Zahraei b a Department of Mathematics, Shahid
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 informationLecture 7 MIMO Communica2ons
Wireless Communications Lecture 7 MIMO Communica2ons Prof. Chun-Hung Liu Dept. of Electrical and Computer Engineering National Chiao Tung University Fall 2014 1 Outline MIMO Communications (Chapter 10
More informationInvariant Theory of Special Orthogonal Groups
Invariant Theory of Special Orthogonal Groups Helmer Aslaksen, Eng-Chye Tan and Chen-bo Zhu Abstract In this paper we study the action of SO.n/ on m-tuples of nn matrices by simultaneous conjugation. We
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 informationDesign of MMSE Multiuser Detectors using Random Matrix Techniques
Design of MMSE Multiuser Detectors using Random Matrix Techniques Linbo Li and Antonia M Tulino and Sergio Verdú Department of Electrical Engineering Princeton University Princeton, New Jersey 08544 Email:
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 informationMassachusetts Institute of Technology Department of Economics Statistics. Lecture Notes on Matrix Algebra
Massachusetts Institute of Technology Department of Economics 14.381 Statistics Guido Kuersteiner Lecture Notes on Matrix Algebra These lecture notes summarize some basic results on matrix algebra used
More informationCommuting nilpotent matrices and pairs of partitions
Commuting nilpotent matrices and pairs of partitions Roberta Basili Algebraic Combinatorics Meets Inverse Systems Montréal, January 19-21, 2007 We will explain some results on commuting n n matrices and
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 informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationStatistical Inference and Random Matrices
Statistical Inference and Random Matrices N.S. Witte Institute of Fundamental Sciences Massey University New Zealand 5-12-2017 Joint work with Peter Forrester 6 th Wellington Workshop in Probability and
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 informationA note on a nearest polynomial with a given root
ACM SIGSAM Bulletin, Vol 39, No. 2, June 2005 A note on a nearest polynomial with a given root Stef Graillat Laboratoire LP2A, Université de Perpignan Via Domitia 52, avenue Paul Alduy F-66860 Perpignan
More informationA note on eigenvalue computation for a tridiagonal matrix with real eigenvalues Akiko Fukuda
Journal of Math-for-Industry Vol 3 (20A-4) pp 47 52 A note on eigenvalue computation for a tridiagonal matrix with real eigenvalues Aio Fuuda Received on October 6 200 / Revised on February 7 20 Abstract
More informationON CHARACTERS AND DIMENSION FORMULAS FOR REPRESENTATIONS OF THE LIE SUPERALGEBRA
ON CHARACTERS AND DIMENSION FORMULAS FOR REPRESENTATIONS OF THE LIE SUPERALGEBRA gl(m N E.M. MOENS Department of Applied Mathematics, Ghent University, Krijgslaan 281-S9, B-9000 Gent, Belgium. e-mail:
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 informationFeasibility Conditions for Interference Alignment
Feasibility Conditions for Interference Alignment Cenk M. Yetis Istanbul Technical University Informatics Inst. Maslak, Istanbul, TURKEY Email: cenkmyetis@yahoo.com Tiangao Gou, Syed A. Jafar University
More informationTHE RELATION BETWEEN THE QR AND LR ALGORITHMS
SIAM J. MATRIX ANAL. APPL. c 1998 Society for Industrial and Applied Mathematics Vol. 19, No. 2, pp. 551 555, April 1998 017 THE RELATION BETWEEN THE QR AND LR ALGORITHMS HONGGUO XU Abstract. For an Hermitian
More information#A69 INTEGERS 13 (2013) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES
#A69 INTEGERS 3 (203) OPTIMAL PRIMITIVE SETS WITH RESTRICTED PRIMES William D. Banks Department of Mathematics, University of Missouri, Columbia, Missouri bankswd@missouri.edu Greg Martin Department of
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 informationarxiv: v1 [physics.gen-ph] 18 Dec 2014
July 16, 018 arxiv:141.5581v1 [physics.gen-ph] 18 Dec 014 A Comment on On the Rotation Matrix in Minkowski Space-time by Özdemir and Erdoğdu Arkadiusz Jadczyk and Jerzy Szulga Abstract. WecommentonthearticlebyM.ÖzdemirandM.Erdoğdu
More information