Skip to main content
1 of 3
Peter Dukes
  • 1.1k
  • 1
  • 8
  • 17

How much can I perturb a symmetric stochastic matrix and keep positive solutions?

Suppose $A$ is a symmetric stochastic $n \times n$ matrix with least eigenvalue $\lambda_n>0$. Let $\widehat{A}$ be a real symmetric perturbation.

How large can I take $\epsilon$ such that $\widehat{A}$ is invertible with $\widehat{A}^{-1} \mathbf{1} > \mathbf{0}$ whenever $||A-\widehat{A}||_\infty< \epsilon$?

I expect the answer involves $\lambda_n$ or some sort of conditioning quantity like $||A^{-1}||_\infty$. I have a preliminary result, but my proof is elementary and ridiculous... it uses a series expansion of $(I-A)^{-1}$. I don't think the resulting lower bound on $\epsilon$ is very good. Can Perron-Frobenius do better somehow?

PS: As an added challenge, I wonder if one gets a bigger $\epsilon$ with an assumption on the sign pattern of either $A^{-1}$ or $A-\widehat{A}$.

Peter Dukes
  • 1.1k
  • 1
  • 8
  • 17