Skip to main content
2 of 4
added 88 characters in body
Sascha
  • 536
  • 7
  • 29

Random matrix is positive

This is a follow up question on my previous question here that was on solved in the deterministic setting by Denis Serre. Therefore, I decided to separate it from the random question initially posed in the same question:

We consider a symmetric matrix $A \in \mathbb R^{n \times n}$ defined by

$$A_{ij}= i \delta_{ij} - X_{ij}\varepsilon \,.$$

Here $\delta_{ij}$ is the Kronecker delta and $X_{ij}$ are iid Bernoulli (but of course $X_{ij}=X_{ji}$ to make it hermitian).

We first note that this matrix is not diagonally dominant if the matrix dimension $n$ is large enough, independent of what $\varepsilon$ is.

This is because $\lim_{n \rightarrow \infty} \sum_{i=1}^{n} \vert A_{i,1}\vert=\infty >\vert A_{1,1} \vert.$

Numerical experiments with matrix size $n=50$ or $n=200$ show that the lowest eigenvalue stays far above $0$, if $\varepsilon>0$ is sufficiently small.

More precisely, it is in some interval $[0,96,1.08]$ if we choose $\varepsilon=0.1$ and $n=50$(upper) and $n=200$(lower), as I illustrate in the following plot where I sampled the lowest eigenvalue for 100 realizations:

k=50

k=200

Question: How can one show that $A$ is positive definite independent of the dimension if $\varepsilon$ is sufficiently small but fixed ?

Sascha
  • 536
  • 7
  • 29