Given an n×n symmetric matrix filled with realizations of a uniform random variable taking values in [−1,1], determine the maximum expected eigenvalue using analytic methods.
Running multiple MATLAB scripts of a large n matrix gave me a the following plot:
Observe, the elliptic nature of the 100×100 plot. Intuitively, I suspect the maximum value of the eigenvalues must be √n, but I don't know how to determine this value analytically (or even if √n is indeed true).
I've attempted reading Pastur et Al.'s "Eigenvalue distribution for Large Random matrices" But can't make sense of their work.
In short What is the relationship between n, dimension of the matrix and the maximum expected eigenvalue?
Answer
The Correct solution is n opposed to my original intuition of √n
For the positive case only
Proof
Observe that
tr(A)=n∑i=1aii=n∑i=1λi
Meaning that the sum of the diagonals is simultaneously equal to the sum of the possible eigenvalues in any square matrix.
Now, let m be the largest entry into A (ie. For our case m=1 as all aij normally randomly chosen ∈[−1,1]
Let λmax≤k∑i=1λk=tr(A)
⟹λmax≤tr(A)
aii≤m
Every aij entry is less than the maximum
⟹n∑i=1aii≤n∑i=1m
⟹tr(A)≤mn
λmax≤tr(A)≤mn
λmax≤mn
Special thanks to Rahul (in comments) & Raghav Srinivasan of the University of Toronto.
Alternately, observe that intuitively I thought that √n was the correct solution because in a Wigner distribution, the scaling as n→∞ is √n This can be proven here
No comments:
Post a Comment