Positive solutions for large random linear systems
Abstract
Consider a large linear system where An is a n× n matrix with independent real standard Gaussian entries, 1n is a n× 1 vector of ones and with unknown the n× 1 vector xn satisfyingxn = 1n + 1αnn An xn\, .We investigate the (componentwise) positivity of the solution xn depending on the scaling factor αn as the dimension n goes to ∞. We prove that there is a sharp phase transition at the threshold α*n =2 n: below the threshold (αn 2 n), xn has negative components with probability tending to 1 while above (αn 2 n), all the vector's components are eventually positive with probability tending to 1. At the critical scaling α*n, we provide a heuristics to evaluate the probability that xn is positive.Such linear systems arise as solutions at equilibrium of large Lotka-Volterra systems of differential equations, widely used to describe large biological communities with interactions such as foodwebs for instance. In the domaine of positivity of the solution xn, that is when αn 2 n, we establish that the Lotka-Volterra system of differential equations whose solution at equilibrium is precisely xn is stable in the sense that its jacobian J(xn) = diag(xn)(-In + Anαnn) has all its eigenvalues with negative real part with probability tending to one. Our results shed a new light and complement the understanding of feasibility and stability issues for large biological communities with interaction.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.