Progress in the method of Ghosts and Shadows for Beta Ensembles
|
|
- Curtis Ross
- 5 years ago
- Views:
Transcription
1 Progress in the method of Ghosts and Shadows for Beta Ensembles Alan Edelman (MIT) Alex Dubbs (MIT) and Plamen Koev (SJS) Oct 8, /47
2 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians Σ=mxn semidefinite matrix G G Σ is similar to A=Σ ½ G GΣ ½ For β=1,2,4, the joint eigenvalue density of A has a formula: Known for β=2 in some circles as Harish Chandra Itzykson Zuber 2/47
3 Main Purpose of this talk Eigenvalue density of G G Σ ( similar to A=Σ ½ G GΣ ½ ) Present an algorithm for sampling from this density Show how the method of Ghosts and Shadows can be used to derive this algorithm Further evidence that β=1,2,4 need not be special 3/47
4 Eigenvalues of GOE (β=1) Naïve Way: MATLAB : A=randn(n); S=(A+A )/sqrt(2*n);eig(s) R: A=matrix(rnorm(n*n),ncol=n);S=(a+t(a))/sqrt(2*n);eigen(S,symmetric=T,only.values=T)$values; Mathematica: A=RandomArray[NormalDistribution[],{n,n}];S=(A+Transpose[A])/Sqrt[n];Eigenvalues[s] 4/47
5 Eigenvalues of GOE/GUE/GSE and beyond Real Sym. Tridiagonal: β=1:real, β=2:complex, β=4:quaternion Diagonals N(0,2), Off diagonals chi random variables General β>0 Computations significantly faster 5/47
6 Introduction to Ghosts G 1 is a standard normal N(0,1) G 2 is a complex normal (G 1 +ig 1 ) G 4 is a quaternion normal (G 1 +ig 1 +jg 1 +kg 1 ) G β (β>0) seems to often work just fine Ghost Gaussian 6/47
7 Chi squared 2 Defn: χ β is the sum of β iid squares of standard normals if β=1,2, Generalizes for non integer β as the gamma function interpolates factorial χ β is the sqrt of the sum of squares (which generalizes) (wikipedia chi distriubtion) G 1 is χ 1, G 2 is χ 2, G 4 is χ 4 So why not G β is χ β? I call χ β the shadow of G β 7/47
8 Linear Algebra Example Given a ghost Gaussian, G β, G β is a real chibeta variable χ β variable. Gram Schmidt: G β G β G β χ 3β G β G β χ 3β G β G β χ 3β G β G β G β G β G β G β G β χ 2β G β H 3 H 2 H 1 χ 2β G β G β G β G β G β G β G β Χ β or G β G β G β G β G β G β G β G β G β = [Q] * χ 3β G β G β χ 2β G β Χ β 8/47
9 Linear Algebra Example Given a ghost Gaussian, G β, G β is a real chibeta variable χ β variable. Gram Schmidt: G β G β G β χ 3β G β G β χ 3β G β G β χ 3β G β G β G β G β G β G β G β χ 2β G β H 3 H 2 H 1 χ 2β G β G β G β G β G β G β G β Χ β or G β G β G β G β G β G β G β G β G β = [Q] * χ 3β G β G β χ 2β G β Χ β Q has β Haar Measure! We have computed moments! 9/47
10 Tridiagonalizing Example G β G β G β G β Symmetric Part of G β G β G β G β G β G β G β G β G β G β G β G β G β (Silverstein, Trotter, etc) 10/47
11 Histogram without Histogramming: Sturm Sequences Count #eigs < 0.5: Count sign changes in Det( (A 0.5*I)[1:k,1:k] ) Count #eigs in [x,x+h] Take difference in number of sign changes at x+h and x Mentioned in Dumitriu and E 2006, Used theoretically in Albrecht, Chan, and E /47
12 A good computational trick is a good theoretical trick! 12/47
13 Stochastic Differential Eigen Equations Tridiagonal Models Suggest SDE s with Brownian motion as infinite limit: E: 2002 E and Sutton 2005, 2007 Brian Rider and (Ramirez, Cambronero, Virag, etc.) (Lots of beautiful results!) (Not today s talk) 13/47
14 Tracy Widom Computations Eigenvalues without the whole matrix! Never construct the entire tridiagonal matrix! Just say the upper 10n 1/3 by 10n 1/3 Compute largest eigenvalue of that, perhaps using Lanczos with shift and invert strategy! Can compute for n amazingly large!!!!! And any beta. E: 2003 Persson and E: /47
15 Other Computational Results Infinite Random Matrix Theory The Free Probability Calculator (Raj) Finite Random Matrix Theory MOPS (Ioana Dumitriu) (β: orthogonal polynomials) Hypergeometrics of Matrix Argument (Plamen Koev) (β: distribution functions for finite stats such as the finite Tracy Widom laws for Laguerre) Good stuff, but not today. 15/47
16 Scary Ideas in Mathematics Zero Negative Radical Irrational Imaginary Ghosts: Something like a commutative algebra of random variables that generalizes random Reals, Complexes, and Quaternions and inspires theoretical results and numerical computation 16/47
17 Did you say commutative?? Quaternions don t commute. Yes but random quaternions do! If x and y are G 4 then x*y and y*x are identically distributed! 17/47
18 Hermite: RMT Densities c λ i λ j β e λ i 2 /2 (Gaussian Ensemble) Laguerre: Jacobi: c λ i λ j β λ m i e λ i (Wishart Matrices) c λ i λ j β λ m 1 i (1 λi ) m 2 (Manova Matrices) Fourier: c λ i λ j β (on the complex unit circle) Jack Polynomials Traditional Story: Count the real parameters β=1,2,4 for real, complex, quaternion Applications: β=1: All of Multivariate Statistics β=2:tons, β=4: Almost nobody cares Dyson 1962 Threefold Way Three Associative Division Rings 18/47
19 β Ghosts β=1 has one real Gaussian (G) β=2 has two real Gaussians (G+iG) β=4 has four real Gaussians (G+iG+jG+kG) β 1 has one real part and (β 1) Ghost parts 19/47
20 Introductory Theory There is an advanced theory emerging A β ghost is a spherically symmetric random variable defined on R β A shadow is a derived real or complex quantity 20/47
21 Wishart Matrices (arbitrary covariance) G=mxn matrix of Gaussians Σ=mxn semidefinite matrix G G Σ is similar to A=Σ ½ G GΣ ½ For β=1,2,4, the joint eigenvalue density of A has a formula: 21/47
22 Joint Eigenvalue density of G G Σ The 0 F 0 function is a hypergeometric function of two matrix arguments that depends only on the eigenvalues of the matrices. Formulas and software exist. 22/47
23 Laguerre: Generalization of Laguerre Versus Wishart: 23/47
24 General β? The joint density: is a probability density for all β>0. Goals: Algorithm for sampling from this density Get a feel for the density s ghost meaning 24/47
25 Main Result An algorithm derived from ghosts that samples eigenvalues A MATLAB implementation that is consistent with other beta ized formulas Largest Eigenvalue Smallest Eigenvalue 25/47
26 Real quantity Working with Ghosts 26/47
27 More practice with Ghosts 27/47
28 Bidiagonalizing Σ=I Z Z has the Σ=I density giving a special case of 28/47
29 The Algorithm for Z=GΣ ½ 29/47
30 The Algorithm for Z=GΣ ½ 30/47
31 Removing U and V 31/47
32 Algorithm cont. 32/47
33 Completion of Recursion 33/47
34 Numerical Experiments Largest Eigenvalue Analytic Formula for largest eigenvalue dist E and Koev: software to compute 34/47
35 1 m3n3beta5.000m150.stag.a.fig F(x) x 35 35/47
36 1 m4n4beta2.500m130.stag.a.fig F(x) x 36/47
37 1 m5n4beta0.750m a.fig F(x) x 37 37/47
38 Smallest Eigenvalue as Well The cdf of the smallest eigenvalue, 38/47
39 Cdf s of smallest m5n4beta3.000.stag.a.least.fig eigenvalue F(x) x 39 39/47
40 Goals Continuum of Haar Measures generalizing orthogonal, unitary, symplectic New Definition of Jack Polynomials generalizing the zonals Computations! E.g. Moments of the Haar Measures Place finite random matrix theory β into same framework as infinite random matrix theory: specifically β as a knob to turn down the randomness, e.g. Airy Kernel d 2 /dx 2 +x+(2/β ½ )dw White Noise 40/47
41 Formally Let S n =2π/Γ(n/2)= surface area of sphere Defined at any n= β>0. A β ghost x is formally defined by a function f x (r) such that f x (r) r β 1 S β 1 dr=1. r=0 Note: For β integer, the x can be realized as a random spherically symmetric variable in β dimensions Example: A β normal ghost is defined by f(r)=(2π) β/2 e r2 /2 Example: Zero is defined with constant*δ(r). Can we do algebra? Can we do linear algebra? Can we add? Can we multiply? 41/47
42 A few more operations ΙΙxΙΙ is a real random variable whose density is given by f x (r) (x+x )/2 is real random variable given by multiplying ΙΙxΙΙ by a beta distributed random variable representing a coordinate on the sphere 42/47
43 Addition of Independent Ghosts: Addition returns a spherically symmetric object Have an integral formula Prefer: Add the real part, imaginary part completed to keep spherical symmetry 43/47
44 Multiplication of Independent Ghosts Just multiply ΙΙzΙΙ s and plug in spherical symmetry Multiplication is commutative (Important Example: Quaternions don t commute, but spherically symmetric random variables do!) 44/47
45 Understanding λ i λ j β Define volume element (dx)^ by (r dx)^=r β (dx)^ (β dim volume, like fractals, but don t really see any fractal theory here) Jacobians: A=QΛQ (Sym Eigendecomposition) Q daq=dλ+(q dq)λ Λ(Q dq) (da)^=(q daq)^= diagonal ^ strictly upper diagonal = dλ i =(dλ)^ off diag = ((Q dq) ij (λ i λ j ))^=(Q dq)^ λ i λ j β 45/47
46 Haar Measure β=1: E Q (trace(aqbq ) k )= C κ (A)C κ (B)/C κ (I) Forward Method: Suppose you know C κ s a priori. (Jack Polynomials!) Let A and B be diagonal indeterminants (Think Generating Functions) Then can formally obtain moments of Q: Example: E( q 11 2 q 22 2 ) = (n+α 1)/(n(n 1)(n+α)) α:=2/ β Can Gram Schmidt the ghosts. Same answers coming up! 46/47
47 Further Uses of Ghosts Multivariate Orthogonal Polynomials Tracy Widom Laws Largest Eigenvalues/Smallest Eigenvalues Expect lots of uses to be discovered 47/47
Advances 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 informationc 2005 Society for Industrial and Applied Mathematics
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
More informationLow-temperature random matrix theory at the soft edge. Abstract
Low-temperature random matrix theory at the soft edge Alan Edelman Department of Mathematics, Massachusetts Institute of Technology, Cambridge, MA 239, USA Per-Olof Persson Department of Mathematics, University
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 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 informationFluctuations from the Semicircle Law Lecture 4
Fluctuations from the Semicircle Law Lecture 4 Ioana Dumitriu University of Washington Women and Math, IAS 2014 May 23, 2014 Ioana Dumitriu (UW) Fluctuations from the Semicircle Law Lecture 4 May 23, 2014
More informationAlan Edelman. Why are random matrix eigenvalues cool? MAA Mathfest 2002 Thursday, August 1. MIT: Dept of Mathematics, Lab for Computer Science
Why are random matrix eigenvalues cool? Alan Edelman MIT: Dept of Mathematics, Lab for Computer Science MAA Mathfest 2002 Thursday, August 6/2/2004 Message Ingredient: Take Any important mathematics Then
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 informationUpdate on the beta ensembles
Update on the beta ensembles Brian Rider Temple University with M. Krishnapur IISC, J. Ramírez Universidad Costa Rica, B. Virág University of Toronto The Tracy-Widom laws Consider a random Hermitian n
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 informationBulk scaling limits, open questions
Bulk scaling limits, open questions Based on: Continuum limits of random matrices and the Brownian carousel B. Valkó, B. Virág. Inventiones (2009). Eigenvalue statistics for CMV matrices: from Poisson
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 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 informationFluctuations of random tilings and discrete Beta-ensembles
Fluctuations of random tilings and discrete Beta-ensembles Alice Guionnet CRS (E S Lyon) Advances in Mathematics and Theoretical Physics, Roma, 2017 Joint work with A. Borodin, G. Borot, V. Gorin, J.Huang
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 informationNumerical Methods for Random Matrices
Numerical Methods for Random Matrices MIT 18.95 IAP Lecture Series Per-Olof Persson (persson@mit.edu) January 23, 26 1.9.8.7 Random Matrix Eigenvalue Distribution.7.6.5 β=1 β=2 β=4 Probability.6.5.4.4
More informationMATHEMATICS (MATH) Calendar
MATHEMATICS (MATH) This is a list of the Mathematics (MATH) courses available at KPU. For information about transfer of credit amongst institutions in B.C. and to see how individual courses transfer, go
More informationConcentration Inequalities for Random Matrices
Concentration Inequalities for Random Matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify the asymptotic
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part II Manuela Girotti based on M. Girotti s PhD thesis, A. Kuijlaars notes from Les Houches Winter School 202 and B. Eynard s notes
More informationPainlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE
Painlevé Representations for Distribution Functions for Next-Largest, Next-Next-Largest, etc., Eigenvalues of GOE, GUE and GSE Craig A. Tracy UC Davis RHPIA 2005 SISSA, Trieste 1 Figure 1: Paul Painlevé,
More informationNumerical analysis and random matrix theory. Tom Trogdon UC Irvine
Numerical analysis and random matrix theory Tom Trogdon ttrogdon@math.uci.edu UC Irvine Acknowledgements This is joint work with: Percy Deift Govind Menon Sheehan Olver Raj Rao Numerical analysis and random
More informationMIT Final Exam Solutions, Spring 2017
MIT 8.6 Final Exam Solutions, Spring 7 Problem : For some real matrix A, the following vectors form a basis for its column space and null space: C(A) = span,, N(A) = span,,. (a) What is the size m n of
More informationRings If R is a commutative ring, a zero divisor is a nonzero element x such that xy = 0 for some nonzero element y R.
Rings 10-26-2008 A ring is an abelian group R with binary operation + ( addition ), together with a second binary operation ( multiplication ). Multiplication must be associative, and must distribute over
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationFinite Rank Perturbations of Random Matrices and Their Continuum Limits. Alexander Bloemendal
Finite Rank Perturbations of Random Matrices and Their Continuum Limits by Alexander Bloemendal A thesis submitted in conformity with the requirements for the degree of Doctor of Philosophy Graduate Department
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 informationFluctuations of random tilings and discrete Beta-ensembles
Fluctuations of random tilings and discrete Beta-ensembles Alice Guionnet CRS (E S Lyon) Workshop in geometric functional analysis, MSRI, nov. 13 2017 Joint work with A. Borodin, G. Borot, V. Gorin, J.Huang
More informationOrthogonal Polynomial Ensembles
Chater 11 Orthogonal Polynomial Ensembles 11.1 Orthogonal Polynomials of Scalar rgument Let wx) be a weight function on a real interval, or the unit circle, or generally on some curve in the comlex lane.
More informationDifferential Equations for Dyson Processes
Differential Equations for Dyson Processes Joint work with Harold Widom I. Overview We call Dyson process any invariant process on ensembles of matrices in which the entries undergo diffusion. Dyson Brownian
More informationEigenvalues and Singular Values of Random Matrices: A Tutorial Introduction
Random Matrix Theory and its applications to Statistics and Wireless Communications Eigenvalues and Singular Values of Random Matrices: A Tutorial Introduction Sergio Verdú Princeton University National
More informationThe Tracy-Widom distribution is not infinitely divisible.
The Tracy-Widom distribution is not infinitely divisible. arxiv:1601.02898v1 [math.pr] 12 Jan 2016 J. Armando Domínguez-Molina Facultad de Ciencias Físico-Matemáticas Universidad Autónoma de Sinaloa, México
More informationThe Hilbert Space of Random Variables
The Hilbert Space of Random Variables Electrical Engineering 126 (UC Berkeley) Spring 2018 1 Outline Fix a probability space and consider the set H := {X : X is a real-valued random variable with E[X 2
More informationMODEL ANSWERS TO THE FIRST HOMEWORK
MODEL ANSWERS TO THE FIRST HOMEWORK 1. Chapter 4, 1: 2. Suppose that F is a field and that a and b are in F. Suppose that a b = 0, and that b 0. Let c be the inverse of b. Multiplying the equation above
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 informationFree Meixner distributions and random matrices
Free Meixner distributions and random matrices Michael Anshelevich July 13, 2006 Some common distributions first... 1 Gaussian Negative binomial Gamma Pascal chi-square geometric exponential 1 2πt e x2
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 informationHow long does it take to compute the eigenvalues of a random symmetric matrix?
How long does it take to compute the eigenvalues of a random symmetric matrix? Govind Menon (Brown University) Joint work with: Christian Pfrang (Ph.D, Brown 2011) Percy Deift, Tom Trogdon (Courant Institute)
More informationLinear Algebra Review
Chapter 1 Linear Algebra Review It is assumed that you have had a course in linear algebra, and are familiar with matrix multiplication, eigenvectors, etc. I will review some of these terms here, but quite
More informationComputer Vision Group Prof. Daniel Cremers. 2. Regression (cont.)
Prof. Daniel Cremers 2. Regression (cont.) Regression with MLE (Rep.) Assume that y is affected by Gaussian noise : t = f(x, w)+ where Thus, we have p(t x, w, )=N (t; f(x, w), 2 ) 2 Maximum A-Posteriori
More informationLinear Algebra for Machine Learning. Sargur N. Srihari
Linear Algebra for Machine Learning Sargur N. srihari@cedar.buffalo.edu 1 Overview Linear Algebra is based on continuous math rather than discrete math Computer scientists have little experience with it
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 informationMarkov operators, classical orthogonal polynomial ensembles, and random matrices
Markov operators, classical orthogonal polynomial ensembles, and random matrices M. Ledoux, Institut de Mathématiques de Toulouse, France 5ecm Amsterdam, July 2008 recent study of random matrix and random
More informationMODEL ANSWERS TO HWK #7. 1. Suppose that F is a field and that a and b are in F. Suppose that. Thus a = 0. It follows that F is an integral domain.
MODEL ANSWERS TO HWK #7 1. Suppose that F is a field and that a and b are in F. Suppose that a b = 0, and that b 0. Let c be the inverse of b. Multiplying the equation above by c on the left, we get 0
More informationRandom Matrix Pencils, Branching Points, and Monodromy
Random Matrix Pencils, Branching Points, and Monodromy Mathew Ellis, Ian Ford, Emily Heath, Christopher Linden, Mychael Sanchez, Dara Zirlin August, 5 Introduction We define complex projective n-space
More informationChapter 5. The multivariate normal distribution. Probability Theory. Linear transformations. The mean vector and the covariance matrix
Probability Theory Linear transformations A transformation is said to be linear if every single function in the transformation is a linear combination. Chapter 5 The multivariate normal distribution When
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 informationInverse Theory. COST WaVaCS Winterschool Venice, February Stefan Buehler Luleå University of Technology Kiruna
Inverse Theory COST WaVaCS Winterschool Venice, February 2011 Stefan Buehler Luleå University of Technology Kiruna Overview Inversion 1 The Inverse Problem 2 Simple Minded Approach (Matrix Inversion) 3
More informationNon white sample covariance matrices.
Non white sample covariance matrices. S. Péché, Université Grenoble 1, joint work with O. Ledoit, Uni. Zurich 17-21/05/2010, Université Marne la Vallée Workshop Probability and Geometry in High Dimensions
More informationGeometric Dyson Brownian motion and May Wigner stability
Geometric Dyson Brownian motion and May Wigner stability Jesper R. Ipsen University of Melbourne Summer school: Randomness in Physics and Mathematics 2 Bielefeld 2016 joint work with Henning Schomerus
More informationLinear Algebra - Part II
Linear Algebra - Part II Projection, Eigendecomposition, SVD (Adapted from Sargur Srihari s slides) Brief Review from Part 1 Symmetric Matrix: A = A T Orthogonal Matrix: A T A = AA T = I and A 1 = A T
More informationLevel statistics of the correlated and uncorrelated Wishart ensemble. Mario Kieburg
Level statistics of the correlated and uncorrelated Wishart ensemble Mario Kieburg Fakultät für Physik Universität Bielefeld (Germany) X Brunel-Bielefeld Workshop on Random Matrix Theory December 13th,
More informationMathematical Physics Homework 10
Georgia Institute of Technology Mathematical Physics Homework Conner Herndon November, 5 Several types of orthogonal polynomials frequently occur in various physics problems. For instance, Hermite polynomials
More informationStat 206: Linear algebra
Stat 206: Linear algebra James Johndrow (adapted from Iain Johnstone s notes) 2016-11-02 Vectors We have already been working with vectors, but let s review a few more concepts. The inner product of two
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationMATH 205 HOMEWORK #3 OFFICIAL SOLUTION. Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. (a) F = R, V = R 3,
MATH 205 HOMEWORK #3 OFFICIAL SOLUTION Problem 1: Find all eigenvalues and eigenvectors of the following linear transformations. a F = R, V = R 3, b F = R or C, V = F 2, T = T = 9 4 4 8 3 4 16 8 7 0 1
More informationMath Linear Algebra II. 1. Inner Products and Norms
Math 342 - Linear Algebra II Notes 1. Inner Products and Norms One knows from a basic introduction to vectors in R n Math 254 at OSU) that the length of a vector x = x 1 x 2... x n ) T R n, denoted x,
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 information235 Final exam review questions
5 Final exam review questions Paul Hacking December 4, 0 () Let A be an n n matrix and T : R n R n, T (x) = Ax the linear transformation with matrix A. What does it mean to say that a vector v R n is an
More informationEigenvalues and Eigenvectors A =
Eigenvalues and Eigenvectors Definition 0 Let A R n n be an n n real matrix A number λ R is a real eigenvalue of A if there exists a nonzero vector v R n such that A v = λ v The vector v is called an eigenvector
More information8 Eigenvectors and the Anisotropic Multivariate Gaussian Distribution
Eigenvectors and the Anisotropic Multivariate Gaussian Distribution Eigenvectors and the Anisotropic Multivariate Gaussian Distribution EIGENVECTORS [I don t know if you were properly taught about eigenvectors
More informationA Cholesky LR algorithm for the positive definite symmetric diagonal-plus-semiseparable eigenproblem
A Cholesky LR algorithm for the positive definite symmetric diagonal-plus-semiseparable eigenproblem Bor Plestenjak Department of Mathematics University of Ljubljana Slovenia Ellen Van Camp and Marc Van
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 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 informationName Solutions Linear Algebra; Test 3. Throughout the test simplify all answers except where stated otherwise.
Name Solutions Linear Algebra; Test 3 Throughout the test simplify all answers except where stated otherwise. 1) Find the following: (10 points) ( ) Or note that so the rows are linearly independent, so
More informationFinal Review Sheet. B = (1, 1 + 3x, 1 + x 2 ) then 2 + 3x + 6x 2
Final Review Sheet The final will cover Sections Chapters 1,2,3 and 4, as well as sections 5.1-5.4, 6.1-6.2 and 7.1-7.3 from chapters 5,6 and 7. This is essentially all material covered this term. Watch
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 information18.06 Quiz 2 April 7, 2010 Professor Strang
18.06 Quiz 2 April 7, 2010 Professor Strang Your PRINTED name is: 1. Your recitation number or instructor is 2. 3. 1. (33 points) (a) Find the matrix P that projects every vector b in R 3 onto the line
More informationEigenvectors and Hermitian Operators
7 71 Eigenvalues and Eigenvectors Basic Definitions Let L be a linear operator on some given vector space V A scalar λ and a nonzero vector v are referred to, respectively, as an eigenvalue and corresponding
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Algorithms Notes for 2016-10-31 There are several flavors of symmetric eigenvalue solvers for which there is no equivalent (stable) nonsymmetric solver. We discuss four algorithmic ideas: the workhorse
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 informationOne-to-one functions and onto functions
MA 3362 Lecture 7 - One-to-one and Onto Wednesday, October 22, 2008. Objectives: Formalize definitions of one-to-one and onto One-to-one functions and onto functions At the level of set theory, there are
More informationFree Probability Theory and Random Matrices. Roland Speicher Queen s University Kingston, Canada
Free Probability Theory and Random Matrices Roland Speicher Queen s University Kingston, Canada We are interested in the limiting eigenvalue distribution of N N random matrices for N. Usually, large N
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 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 informationCharacterizations of free Meixner distributions
Characterizations of free Meixner distributions Texas A&M University March 26, 2010 Jacobi parameters. Matrix. β 0 γ 0 0 0... 1 β 1 γ 1 0.. m n J =. 0 1 β 2 γ.. 2 ; J n =. 0 0 1 β.. 3............... A
More informationMATRICES ARE SIMILAR TO TRIANGULAR MATRICES
MATRICES ARE SIMILAR TO TRIANGULAR MATRICES 1 Complex matrices Recall that the complex numbers are given by a + ib where a and b are real and i is the imaginary unity, ie, i 2 = 1 In what we describe below,
More information22m:033 Notes: 7.1 Diagonalization of Symmetric Matrices
m:33 Notes: 7. Diagonalization of Symmetric Matrices Dennis Roseman University of Iowa Iowa City, IA http://www.math.uiowa.edu/ roseman May 3, Symmetric matrices Definition. A symmetric matrix is a matrix
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 2018 8 Lecture 8 8.1 Matrices July 22, 2018 We shall study
More informationORTHOGONAL POLYNOMIALS
ORTHOGONAL POLYNOMIALS 1. PRELUDE: THE VAN DER MONDE DETERMINANT The link between random matrix theory and the classical theory of orthogonal polynomials is van der Monde s determinant: 1 1 1 (1) n :=
More informationTOP EIGENVALUE OF CAUCHY RANDOM MATRICES
TOP EIGENVALUE OF CAUCHY RANDOM MATRICES with Satya N. Majumdar, Gregory Schehr and Dario Villamaina Pierpaolo Vivo (LPTMS - CNRS - Paris XI) Gaussian Ensembles N = 5 Semicircle Law LARGEST EIGENVALUE
More informationCS168: The Modern Algorithmic Toolbox Lecture #8: How PCA Works
CS68: The Modern Algorithmic Toolbox Lecture #8: How PCA Works Tim Roughgarden & Gregory Valiant April 20, 206 Introduction Last lecture introduced the idea of principal components analysis (PCA). The
More information1. Introduction to commutative rings and fields
1. Introduction to commutative rings and fields Very informally speaking, a commutative ring is a set in which we can add, subtract and multiply elements so that the usual laws hold. A field is a commutative
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 informationAnswers and Solutions to Selected Homework Problems From Section 2.5 S. F. Ellermeyer. and B =. 0 2
Answers and Solutions to Selected Homework Problems From Section 2.5 S. F. Ellermeyer 5. Since gcd (2; 4) 6, then 2 is a zero divisor (and not a unit) in Z 4. In fact, we see that 2 2 0 in Z 4. Thus 2x
More informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
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 informationMath 291-2: Lecture Notes Northwestern University, Winter 2016
Math 291-2: Lecture Notes Northwestern University, Winter 2016 Written by Santiago Cañez These are lecture notes for Math 291-2, the second quarter of MENU: Intensive Linear Algebra and Multivariable Calculus,
More informationRANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS
RANDOM MATRIX THEORY AND TOEPLITZ DETERMINANTS David García-García May 13, 2016 Faculdade de Ciências da Universidade de Lisboa OVERVIEW Random Matrix Theory Introduction Matrix ensembles A sample computation:
More information1. Introduction to commutative rings and fields
1. Introduction to commutative rings and fields Very informally speaking, a commutative ring is a set in which we can add, subtract and multiply elements so that the usual laws hold. A field is a commutative
More information[y i α βx i ] 2 (2) Q = i=1
Least squares fits This section has no probability in it. There are no random variables. We are given n points (x i, y i ) and want to find the equation of the line that best fits them. We take the equation
More informationReview (Probability & Linear Algebra)
Review (Probability & Linear Algebra) CE-725 : Statistical Pattern Recognition Sharif University of Technology Spring 2013 M. Soleymani Outline Axioms of probability theory Conditional probability, Joint
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More information1. Matrix multiplication and Pauli Matrices: Pauli matrices are the 2 2 matrices. 1 0 i 0. 0 i
Problems in basic linear algebra Science Academies Lecture Workshop at PSGRK College Coimbatore, June 22-24, 2016 Govind S. Krishnaswami, Chennai Mathematical Institute http://www.cmi.ac.in/~govind/teaching,
More informationGetting Started with Communications Engineering
1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we
More informationLinear Algebra review Powers of a diagonalizable matrix Spectral decomposition
Linear Algebra review Powers of a diagonalizable matrix Spectral decomposition Prof. Tesler Math 283 Fall 2016 Also see the separate version of this with Matlab and R commands. Prof. Tesler Diagonalizing
More informationMultiple orthogonal polynomials. Bessel weights
for modified Bessel weights KU Leuven, Belgium Madison WI, December 7, 2013 Classical orthogonal polynomials The (very) classical orthogonal polynomials are those of Jacobi, Laguerre and Hermite. Classical
More informationUniversality of local spectral statistics of random matrices
Universality of local spectral statistics of random matrices László Erdős Ludwig-Maximilians-Universität, Munich, Germany CRM, Montreal, Mar 19, 2012 Joint with P. Bourgade, B. Schlein, H.T. Yau, and J.
More informationReview problems for MA 54, Fall 2004.
Review problems for MA 54, Fall 2004. Below are the review problems for the final. They are mostly homework problems, or very similar. If you are comfortable doing these problems, you should be fine on
More informationNOTES ON LINEAR ALGEBRA CLASS HANDOUT
NOTES ON LINEAR ALGEBRA CLASS HANDOUT ANTHONY S. MAIDA CONTENTS 1. Introduction 2 2. Basis Vectors 2 3. Linear Transformations 2 3.1. Example: Rotation Transformation 3 4. Matrix Multiplication and Function
More information