Alan Edelman. Why are random matrix eigenvalues cool? MAA Mathfest 2002 Thursday, August 1. MIT: Dept of Mathematics, Lab for Computer Science
|
|
- Gavin McKenzie
- 5 years ago
- Views:
Transcription
1 Why are random matrix eigenvalues cool? Alan Edelman MIT: Dept of Mathematics, Lab for Computer Science MAA Mathfest 2002 Thursday, August 6/2/2004
2 Message Ingredient: Take Any important mathematics Then Randomize! This will have many applications! 6/2/2004 2
3 Some fun tidbits The circular law The semi-circular law Infinite vs finite How many are real? Stochastic Numerical Algorithms Condition Numbers Small networks Riemann Zeta Function Matrix Jacobians 6/2/2004 3
4 irko s Circular Law, n=2000 Has anyone studied spacings? 6/2/2004 4
5 Wigner s Semi-Circle The classical & most famous rand eig theorem Let S = random symmetric aussian MATLAB: A=randn(n); S=(A+A )/2; Normalized eigenvalue histogram is a semi-circle Precise statements require n etc. 5
6 Wigner s Semi-Circle The classical & most famous rand eig theorem Let S = random symmetric aussian MATLAB: A=randn(n); S=(A+A )/2; Normalized eigenvalue histogram is a semi-circle Precise statements require n etc. n=20; s=30000; d=.05; %matrix size, samples, sample dist e=[]; %gather up eigenvalues im=; %imaginary() or real(0) for i=:s, a=randn(n)+im*sqrt(-)*randn(n);a=(a+a')/(2*sqrt(2*n*(im+))); v=eig(a)'; e=[e v]; end hold off; [m x]=hist(e,-.5:d:.5); bar(x,m*pi/(2*d*n*s)); axis('square'); axis([ ]); hold on; t=-:.0:; plot(t,sqrt(-t.^2),'r'); 6
7 Elements of Wigner s Proof Compute E(A 2k ) = mean(λ 2k ) = (2k)th moment Verify that the semicircle is the only distribution with these moments (A 2k ) = ΣA x A xy A wz A z paths of length 2k Need only count number of special paths of length 2k on k objects (all other terms 0 or negligible!) This is a Catalan Number! 6/2/2004 7
8 C n = 2n n + n Catalan Numbers =# ways to parenthesize (n+) objects Matrix Power Term raph (((23)4)) A 2 A 23 A 32 A 24 A 42 A 2 (((2)3)4) A 2 A 2 A 3 A 3 A 4 A 4 ((2(34))) A 2 A 23 A 34 A 43 A 32 A 2 ((2)(34)) A 2 A 2 A 3 A 34 A 43 A 3 (((23))4) A 2 A 23 A 32 A 2 A 4 A 4 = number of special paths on n departing from once Pass, (load=advance, multiply=retreat), Return to 6/2/2004 8
9 Finite Versions n=2; n=4; n=3; n=5; 9
10 How many eigenvalues of a random matrix are real? >> e=eig(randn(7)) e = i i i i >> e=eig(randn(7)) e = i i i i i i >> e=eig(randn(7)) e = i i real real 5 real 7x7 random aussian 6/2/2004 0
11 How many eigenvalues of a random matrix are real? 7 reals 5 reals 3 reals real n= /2/
12 How many eigenvalues of a random matrix are real? 7 reals 5 reals 3 reals n= These are exact real but hard to compute! New research suggests a Jack polynomial solution. 6/2/
13 How many eigenvalues of a random matrix are real? The Probability that a matrix has all real eigenvalues is exactly P n,n =2 -n(n-)/4 Proof based on Schur Form 6/2/2004 3
14 ram Schmidt (or QR) Stochastically ram Schmidt = Orthogonal Transformations to Upper Triangular Form A = Q * R (orthog * upper triangular) 6/2/2004 4
15 Orthogonal Invariance of aussians Q*randn(n,) randn(n,) If Q orthogonal Q 6/2/2004 5
16 Orthogonal Invariance Q*randn(n,) randn(n,) If Q orthogonal Q = 6/2/2004 6
17 Chi Distribution norm(randn(n,)) χ n = χ n 6/2/2004 7
18 Chi Distribution norm(randn(n,)) χ n = χ n 6/2/2004 8
19 Chi Distribution norm(randn(n,)) χ n = χ n n need not be integer 6/2/2004 9
20 6/2/
21 6/2/2004 2
22 6/2/
23 6/2/
24 6/2/
25 6/2/
26 6/2/
27 6/2/
28 6/2/
29 6/2/
30 6/2/
31 6/2/2004 3
32 O χ 6 O O O O O O O O O O 6/2/
33 O χ 6 O O O O O O O O O O 6/2/
34 O χ 6 O O O O O O O O O O 6/2/
35 O χ 6 O O O O O O O O O O 6/2/
36 O χ 6 O O O O O O O O O O 6/2/
37 O χ 6 O O O O O O O O O O 6/2/
38 O χ 6 O O O O O O O O O O 6/2/
39 O χ 6 O O O O O O O O O O 6/2/
40 O χ 6 O O O O O O O O O O 6/2/
41 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/2004 4
42 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
43 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
44 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
45 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
46 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
47 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
48 O χ 6 O O χ 5 O O O O O O O O O O O O 6/2/
49 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
50 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
51 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/2004 5
52 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
53 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
54 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
55 O χ 6 O O χ 5 O O O χ 4 O O O O O O O O O O O O 6/2/
56 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/
57 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/
58 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/
59 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/
60 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/
61 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O O O O O O 6/2/2004 6
62 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O 6/2/
63 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O 6/2/
64 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O 6/2/
65 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O 6/2/
66 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O 6/2/
67 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O χ 6/2/
68 O χ 6 O O χ 5 O O O χ 4 O O O O χ 3 O O O O O χ 2 O O O O O O χ 6/2/
69 Same idea: sym matrix to tridiagonal form χ 6 χ 6 χ 5 χ 5 χ 4 χ 4 χ 3 Same eigenvalue distribution as OE: O(n) storage!! χ 3 O(n) computation (potentially) χ 2 χ 2 χ χ 6/2/
70 Same idea: eneral beta χ 6β χ 6β beta: : reals 2: complexes 4: quaternions χ 5β χ 5β χ Bidiagonal 4β Version corresponds χ To Wishart 4β χmatrices of Statistics 3β χ 3β χ 2β χ 2β χ β χ β 6/2/
71 Numerical Analysis: Condition Numbers κ(a) = condition number of A If A=UΣV is the svd, then κ(a) = σ max /σ min. Alternatively, κ(a) = λ max (A A)/ λ min (A A) One number that measures digits lost in finite precision and general matrix badness Small=good Large=bad The condition of a random matrix??? 6/2/2004 7
72 Von Neumann & co. Solve Ax=b via x= (A A) - A b M A - Matrix Residual: AM-I 2 AM-I 2 < 200κ 2 n ε How should we estimate κ? Assume, as a model, that the elements of A are independent standard normals! 6/2/
73 Von Neumann & co. estimates (947-95) For a random matrix of order n the expectation value has been shown to be about n X oldstine, von Neumann we choose two different values of κ, namely n and 0n P(κ<n) 0.02 P(κ< 0n) 0.44 Bargmann, Montgomery, vn With a probability ~ κ < 0n P(κ<0n) 0.80 oldstine, von Neumann 6/2/
74 Random cond numbers, n y Distribution of κ/n = x e 2 / x 2 / x 3 x 2 Experiment with n=200 6/2/
75 Finite n n=0 n=25 n=50 n=00 6/2/
76 76 Small World Networks: & 6 degrees of separation Edelman, Eriksson, Strang Eigenvalues of A=T+PTP, P=randperm(n) Incidence matrix of graph with two superimposed cycles. = T O O O
77 77 Small World Networks: & 6 degrees of separation Edelman, Eriksson, Strang Eigenvalues of A=T+PTP, P=randperm(n) Incidence matrix of graph with two superimposed cycles. Wigner style derivation counts number of paths on a tree starting and ending at the same point (tree = no accidents!) (McKay) We first discovered the formula using the superseeker Catalan number answer d 2n- -Σ d 2j- (d-) n-j+ C n-j = T O O O
78 6/2/ The Riemann Zeta Function = = Γ = 0 ) ( ) ζ( k s u s k du e u s s 6 π ζ(2) = = K On the real line with x>, for example May be analytically extended to the complex plane, with singularity only at x=.
79 6/2/ ½ 2 3 The Riemann Hypothesis = = Γ = 0 ) ( ) ζ( k x u x k du e u x x All nontrivial roots of ζ(x) satisfy Re(x)=/2. (Trivial roots at negative even integers.)
80 ζ( x) ζ(x) along Re(x)=/2 Zeros =.5+i* The Riemann Hypothesis u x = du = Γ( x) e u k x k= ½ All nontrivial roots of ζ(x) satisfy Re(x)=/2. (Trivial roots at negative even integers.) 6/2/
81 Computation of Zeros Odlyzko s fantastic computation of 0^k+ through 0^k+0,000 for k=2,2,22. See Spacings behave like the eigenvalues of A=randn(n)+i*randn(n); S=(A+A )/2; 6/2/2004 8
82 Nearest Neighbor Spacings & Pairwise Correlation Functions 6/2/
83 Painlevé Equations 6/2/
84 Spacings Take a large collection of consecutive zeros/eigenvalues. Normalize so that average spacing =. Spacing Function = Histogram of consecutive differences (the (k+)st the kth) Pairwise Correlation Function = Histogram of all possible differences (the kth the jth) Conjecture: These functions are the same for random matrices and Riemann zeta 6/2/
85 Some fun tidbits The circular law The semi-circular law Infinite vs finite How many are real? Stochastic Numerical Algorithms Condition Numbers Small networks Riemann Zeta Function Matrix Jacobians 6/2/
86 eneral A=LU A=UΣV T A=XΛX - Sym S=QΛQ T S=LL T S=LDL T Tridiagonal Matrix Factorization Jacobians u n-i ii (σ i2 - σ j2 ) (λ i -λ j ) 2 (λ i -λ j ) 2 n l ii n+-i d i n-i A=QR A=QS (polar) Orthogonal C S S -C Q=U V T r ii m-i (σ i +σ j ) [ ] sin(θ i + θ j )sin (θ i - θ j ) T=QΛQ T (t i+,i )/ q i 6/2/
87 Why cool? Why is numerical linear algebra cool? Mixture of theory and applications Touches many topics Easy to jump in to, but can spend a lifetime studying & researching Tons of activity in many areas Mathematics: Combinatorics, Harmonic Analysis, Integral Equations, Probability, Number Theory Applied Math: Chaotic Systems, Statistical Mechanics, Communications Theory, Radar Tracking, Nuclear Physics Applications BI HUE SUBJECT!! 6/2/
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 informationProgress 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 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 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 informationx #{ p=prime p x }, as x logx
1 The Riemann zeta function for Re(s) > 1 ζ s -s ( ) -1 2 duality between primes & complex zeros of zeta using Hadamard product over zeros prime number theorem x #{ p=prime p x }, as x logx statistics
More informationWhat is the Riemann Hypothesis for Zeta Functions of Irregular Graphs?
What is the Riemann Hypothesis for Zeta Functions of Irregular Graphs? Audrey Terras Banff February, 2008 Joint work with H. M. Stark, M. D. Horton, etc. Introduction The Riemann zeta function for Re(s)>1
More informationRandom Matrices: Invertibility, Structure, and Applications
Random Matrices: Invertibility, Structure, and Applications Roman Vershynin University of Michigan Colloquium, October 11, 2011 Roman Vershynin (University of Michigan) Random Matrices Colloquium 1 / 37
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More information8. Diagonalization.
8. Diagonalization 8.1. Matrix Representations of Linear Transformations Matrix of A Linear Operator with Respect to A Basis We know that every linear transformation T: R n R m has an associated standard
More informationLecture 1: Riemann, Dedekind, Selberg, and Ihara Zetas
Lecture 1: Riemann, Dedekind, Selberg, and Ihara Zetas Audrey Terras U.C.S.D. 2008 more details can be found in my webpage: www.math.ucsd.edu /~aterras/ newbook.pdf First the Riemann Zeta 1 The Riemann
More informationComputation of eigenvalues and singular values Recall that your solutions to these questions will not be collected or evaluated.
Math 504, Homework 5 Computation of eigenvalues and singular values Recall that your solutions to these questions will not be collected or evaluated 1 Find the eigenvalues and the associated eigenspaces
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 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 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 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 informationMath 504 (Fall 2011) 1. (*) Consider the matrices
Math 504 (Fall 2011) Instructor: Emre Mengi Study Guide for Weeks 11-14 This homework concerns the following topics. Basic definitions and facts about eigenvalues and eigenvectors (Trefethen&Bau, Lecture
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 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 informationThe Prime Unsolved Problem in Mathematics
The Prime Unsolved Problem in Mathematics March, 0 Augustana Math Club 859: Bernard Riemann proposes the Riemann Hypothesis. 859: Bernard Riemann proposes the Riemann Hypothesis. 866: Riemann dies at age
More informationA fast randomized algorithm for overdetermined linear least-squares regression
A fast randomized algorithm for overdetermined linear least-squares regression Vladimir Rokhlin and Mark Tygert Technical Report YALEU/DCS/TR-1403 April 28, 2008 Abstract We introduce a randomized algorithm
More informationAn Invitation to Modern Number Theory. Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD
An Invitation to Modern Number Theory Steven J. Miller and Ramin Takloo-Bighash PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Foreword Preface Notation xi xiii xix PART 1. BASIC NUMBER THEORY
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 informationQuantum chaos on graphs
Baylor University Graduate seminar 6th November 07 Outline 1 What is quantum? 2 Everything you always wanted to know about quantum but were afraid to ask. 3 The trace formula. 4 The of Bohigas, Giannoni
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 informationarxiv:chao-dyn/ v1 3 Jul 1995
Chaotic Spectra of Classically Integrable Systems arxiv:chao-dyn/9506014v1 3 Jul 1995 P. Crehan Dept. of Mathematical Physics, University College Dublin, Belfield, Dublin 2, Ireland PCREH89@OLLAMH.UCD.IE
More informationRiemann s Zeta Function and the Prime Number Theorem
Riemann s Zeta Function and the Prime Number Theorem Dan Nichols nichols@math.umass.edu University of Massachusetts Dec. 7, 2016 Let s begin with the Basel problem, first posed in 1644 by Mengoli. Find
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 informationMATH 511 ADVANCED LINEAR ALGEBRA SPRING 2006
MATH 511 ADVANCED LINEAR ALGEBRA SPRING 2006 Sherod Eubanks HOMEWORK 2 2.1 : 2, 5, 9, 12 2.3 : 3, 6 2.4 : 2, 4, 5, 9, 11 Section 2.1: Unitary Matrices Problem 2 If λ σ(u) and U M n is unitary, show that
More informationREGULI AND APPLICATIONS, ZARANKIEWICZ PROBLEM
REGULI AN APPLICATIONS, ZARANKIEWICZ PROBLEM One of our long-term goals is to pursue some questions of incidence geometry in R 3. We recall one question to direct our focus during this lecture. Question
More information(Mathematical Operations with Arrays) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Mathematical Operations with Arrays) Contents Getting Started Matrices Creating Arrays Linear equations Mathematical Operations with Arrays Using Script
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 information4. Determinants.
4. Determinants 4.1. Determinants; Cofactor Expansion Determinants of 2 2 and 3 3 Matrices 2 2 determinant 4.1. Determinants; Cofactor Expansion Determinants of 2 2 and 3 3 Matrices 3 3 determinant 4.1.
More informationNext is material on matrix rank. Please see the handout
B90.330 / C.005 NOTES for Wednesday 0.APR.7 Suppose that the model is β + ε, but ε does not have the desired variance matrix. Say that ε is normal, but Var(ε) σ W. The form of W is W w 0 0 0 0 0 0 w 0
More informationFrom random matrices to free groups, through non-crossing partitions. Michael Anshelevich
From random matrices to free groups, through non-crossing partitions Michael Anshelevich March 4, 22 RANDOM MATRICES For each N, A (N), B (N) = independent N N symmetric Gaussian random matrices, i.e.
More informationReview of similarity transformation and Singular Value Decomposition
Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm
More informationBeyond the Gaussian universality class
Beyond the Gaussian universality class MSRI/Evans Talk Ivan Corwin (Courant Institute, NYU) September 13, 2010 Outline Part 1: Random growth models Random deposition, ballistic deposition, corner growth
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 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 informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
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 informationMath 108b: Notes on the Spectral Theorem
Math 108b: Notes on the Spectral Theorem From section 6.3, we know that every linear operator T on a finite dimensional inner product space V has an adjoint. (T is defined as the unique linear operator
More informationMath 307 Learning Goals. March 23, 2010
Math 307 Learning Goals March 23, 2010 Course Description The course presents core concepts of linear algebra by focusing on applications in Science and Engineering. Examples of applications from recent
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 informationWhat is the Riemann Hypothesis for Zeta Functions of Irregular Graphs?
What is the Riemann Hypothesis for Zeta Functions of Irregular Graphs? UCLA, IPAM February, 2008 Joint work with H. M. Stark, M. D. Horton, etc. What is an expander graph X? X finite connected (irregular)
More informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More informationCUTOFF FOR THE STAR TRANSPOSITION RANDOM WALK
CUTOFF FOR THE STAR TRANSPOSITION RANDOM WALK JONATHAN NOVAK AND ARIEL SCHVARTZMAN Abstract. In this note, we give a complete proof of the cutoff phenomenon for the star transposition random walk on the
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 informationNumerical Methods I Eigenvalue Problems
Numerical Methods I Eigenvalue Problems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 2nd, 2014 A. Donev (Courant Institute) Lecture
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 informationProperties of Linear Transformations from R n to R m
Properties of Linear Transformations from R n to R m MATH 322, Linear Algebra I J. Robert Buchanan Department of Mathematics Spring 2015 Topic Overview Relationship between the properties of a matrix transformation
More informationMath 259: Introduction to Analytic Number Theory How small can disc(k) be for a number field K of degree n = r 1 + 2r 2?
Math 59: Introduction to Analytic Number Theory How small can disck be for a number field K of degree n = r + r? Let K be a number field of degree n = r + r, where as usual r and r are respectively the
More informationEigenvalue and Eigenvector Problems
Eigenvalue and Eigenvector Problems An attempt to introduce eigenproblems Radu Trîmbiţaş Babeş-Bolyai University April 8, 2009 Radu Trîmbiţaş ( Babeş-Bolyai University) Eigenvalue and Eigenvector Problems
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
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 informationEquidistributions in arithmetic geometry
Equidistributions in arithmetic geometry Edgar Costa Dartmouth College 14th January 2016 Dartmouth College 1 / 29 Edgar Costa Equidistributions in arithmetic geometry Motivation: Randomness Principle Rigidity/Randomness
More informationFundamentals of Matrices
Maschinelles Lernen II Fundamentals of Matrices Christoph Sawade/Niels Landwehr/Blaine Nelson Tobias Scheffer Matrix Examples Recap: Data Linear Model: f i x = w i T x Let X = x x n be the data matrix
More informationOn the -Sylvester Equation AX ± X B = C
Department of Applied Mathematics National Chiao Tung University, Taiwan A joint work with E. K.-W. Chu and W.-W. Lin Oct., 20, 2008 Outline 1 Introduction 2 -Sylvester Equation 3 Numerical Examples 4
More informationcourse overview 18.06: Linear Algebra
course overview 18.06: Linear Algebra Prof. Steven G. Johnson, MIT Applied Math Fall 2017 http://web.mit.edu/18.06 Textbook: Strang, Introduction to Linear Algebra, 5 th edition + supplementary notes Help
More informationMathematics of Imaging: Lecture 3
Mathematics of Imaging: Lecture 3 Linear Operators in Infinite Dimensions Consider the linear operator on the space of continuous functions defined on IR. J (f)(x) = x 0 f(s) ds Every function in the range
More informationMath 495 Dr. Rick Kreminski Colorado State University Pueblo November 19, 2014 The Riemann Hypothesis: The #1 Problem in Mathematics
Math 495 Dr. Rick Kreminski Colorado State University Pueblo November 19, 14 The Riemann Hypothesis: The #1 Problem in Mathematics Georg Friedrich Bernhard Riemann 1826-1866 One of 6 Million Dollar prize
More informationThe value of a problem is not so much coming up with the answer as in the ideas and attempted ideas it forces on the would be solver I.N.
Math 410 Homework Problems In the following pages you will find all of the homework problems for the semester. Homework should be written out neatly and stapled and turned in at the beginning of class
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 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 informationRiemann Zeta Function and Prime Number Distribution
Riemann Zeta Function and Prime Number Distribution Mingrui Xu June 2, 24 Contents Introduction 2 2 Definition of zeta function and Functional Equation 2 2. Definition and Euler Product....................
More informationOrthogonal iteration to QR
Notes for 2016-03-09 Orthogonal iteration to QR The QR iteration is the workhorse for solving the nonsymmetric eigenvalue problem. Unfortunately, while the iteration itself is simple to write, the derivation
More informationThe Free Central Limit Theorem: A Combinatorial Approach
The Free Central Limit Theorem: A Combinatorial Approach by Dennis Stauffer A project submitted to the Department of Mathematical Sciences in conformity with the requirements for Math 4301 (Honour s Seminar)
More informationVII.8. The Riemann Zeta Function.
VII.8. The Riemann Zeta Function VII.8. The Riemann Zeta Function. Note. In this section, we define the Riemann zeta function and discuss its history. We relate this meromorphic function with a simple
More informationBefore giving the detailed proof, we outline our strategy. Define the functions. for Re s > 1.
Chapter 7 The Prime number theorem for arithmetic progressions 7. The Prime number theorem We denote by π( the number of primes. We prove the Prime Number Theorem. Theorem 7.. We have π( as. log Before
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 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 informationAlgebra 2 Summer Assignment - Non Regents
Name: Date: This assignment is to be completed over the summer. It must be turned in the first day of classes. It will be counted as a quiz. You MUST SHOW all work for each question in the space provided
More informationGraph fundamentals. Matrices associated with a graph
Graph fundamentals Matrices associated with a graph Drawing a picture of a graph is one way to represent it. Another type of representation is via a matrix. Let G be a graph with V (G) ={v 1,v,...,v n
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationRandom matrices and the Riemann zeros
Random matrices and the Riemann zeros Bruce Bartlett Talk given in Postgraduate Seminar on 5th March 2009 Abstract Random matrices and the Riemann zeta function came together during a chance meeting between
More informationLecture 6. Numerical methods. Approximation of functions
Lecture 6 Numerical methods Approximation of functions Lecture 6 OUTLINE 1. Approximation and interpolation 2. Least-square method basis functions design matrix residual weighted least squares normal equation
More informationLinear Algebra and Matrices
Linear Algebra and Matrices 4 Overview In this chapter we studying true matrix operations, not element operations as was done in earlier chapters. Working with MAT- LAB functions should now be fairly routine.
More informationRandom Matrices Numerical Methods for Random Matrices. Per-Olof Persson
18.325 Random Matrices Numerical Methods for Random Matrices Per-Olof Persson (persson@mit.edu) December 19, 22 1 Largest Eigenvalue Distributions In this section, the distributions of the largest eigenvalue
More informationTuring and the Riemann zeta function
Turing and the Riemann zeta function Andrew Odlyzko School of Mathematics University of Minnesota odlyzko@umn.edu http://www.dtc.umn.edu/ odlyzko May 11, 2012 Andrew Odlyzko ( School of Mathematics University
More informationLECTURE 4. PROOF OF IHARA S THEOREM, EDGE CHAOS. Ihara Zeta Function. ν(c) ζ(u,x) =(1-u ) det(i-au+qu t(ia )
LCTUR 4. PROOF OF IHARA S THORM, DG ZTAS, QUANTUM CHAOS Ihara Zeta Function ν(c) ( ) - ζ(u,x)= -u [C] prime ν(c) = # edges in C converges for u complex, u small Ihara s Theorem. - 2 r- 2 ζ(u,x) =(-u )
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 information33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM
33AH, WINTER 2018: STUDY GUIDE FOR FINAL EXAM (UPDATED MARCH 17, 2018) The final exam will be cumulative, with a bit more weight on more recent material. This outline covers the what we ve done since the
More informationMOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS
MOMENTS OF HYPERGEOMETRIC HURWITZ ZETA FUNCTIONS ABDUL HASSEN AND HIEU D. NGUYEN Abstract. This paper investigates a generalization the classical Hurwitz zeta function. It is shown that many of the properties
More informationECS130 Scientific Computing Handout E February 13, 2017
ECS130 Scientific Computing Handout E February 13, 2017 1. The Power Method (a) Pseudocode: Power Iteration Given an initial vector u 0, t i+1 = Au i u i+1 = t i+1 / t i+1 2 (approximate eigenvector) θ
More informationRandom matrices: A Survey. Van H. Vu. Department of Mathematics Rutgers University
Random matrices: A Survey Van H. Vu Department of Mathematics Rutgers University Basic models of random matrices Let ξ be a real or complex-valued random variable with mean 0 and variance 1. Examples.
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationNeedles and Numbers. The Buffon Needle Experiment
eedles and umbers This excursion into analytic number theory is intended to complement the approach of our textbook, which emphasizes the algebraic theory of numbers. At some points, our presentation lacks
More informationPositive entries of stable matrices
Positive entries of stable matrices Shmuel Friedland Department of Mathematics, Statistics and Computer Science, University of Illinois at Chicago Chicago, Illinois 60607-7045, USA Daniel Hershkowitz,
More informationThe circular law. Lewis Memorial Lecture / DIMACS minicourse March 19, Terence Tao (UCLA)
The circular law Lewis Memorial Lecture / DIMACS minicourse March 19, 2008 Terence Tao (UCLA) 1 Eigenvalue distributions Let M = (a ij ) 1 i n;1 j n be a square matrix. Then one has n (generalised) eigenvalues
More informationRecurrence Relations and Recursion: MATH 180
Recurrence Relations and Recursion: MATH 180 1: Recursively Defined Sequences Example 1: The sequence a 1,a 2,a 3,... can be defined recursively as follows: (1) For all integers k 2, a k = a k 1 + 1 (2)
More informationRandom Matrix: From Wigner to Quantum Chaos
Random Matrix: From Wigner to Quantum Chaos Horng-Tzer Yau Harvard University Joint work with P. Bourgade, L. Erdős, B. Schlein and J. Yin 1 Perhaps I am now too courageous when I try to guess the distribution
More informationEigenvalues and eigenvectors
Chapter 6 Eigenvalues and eigenvectors An eigenvalue of a square matrix represents the linear operator as a scaling of the associated eigenvector, and the action of certain matrices on general vectors
More informationRandom matrices: Distribution of the least singular value (via Property Testing)
Random matrices: Distribution of the least singular value (via Property Testing) Van H. Vu Department of Mathematics Rutgers vanvu@math.rutgers.edu (joint work with T. Tao, UCLA) 1 Let ξ be a real or complex-valued
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 7: More on Householder Reflectors; Least Squares Problems Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 15 Outline
More informationMathematics 324 Riemann Zeta Function August 5, 2005
Mathematics 324 Riemann Zeta Function August 5, 25 In this note we give an introduction to the Riemann zeta function, which connects the ideas of real analysis with the arithmetic of the integers. Define
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 informationRandomness in Number Theory
Randomness in Number Theory Peter Sarnak Mahler Lectures 2011 Number Theory Probability Theory Whole numbers Random objects Prime numbers Points in space Arithmetic operations Geometries Diophantine equations
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 informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationBasic Fraction and Integer Operations (No calculators please!)
P1 Summer Math Review Packet For Students entering Geometry The problems in this packet are designed to help you review topics from previous mathematics courses that are important to your success in Geometry.
More informationLecture 1: Small Prime Gaps: From the Riemann Zeta-Function and Pair Correlation to the Circle Method. Daniel Goldston
Lecture 1: Small Prime Gaps: From the Riemann Zeta-Function and Pair Correlation to the Circle Method Daniel Goldston π(x): The number of primes x. The prime number theorem: π(x) x log x, as x. The average
More information