Positive solutions for large random linear systems

From MaRDI portal



Abstract: Consider a large linear system where An is a nimesn matrix with independent real standard Gaussian entries, is a nimes1 vector of ones and with unknown the nimes1 vector satisfying�oldsymbol{x}_n = �oldsymbol{1}_n +frac 1{alpha_nsqrt{n}} A_n �oldsymbol{x}_n, .We investigate the (componentwise) positivity of the solution depending on the scaling factor alphan as the dimension n goes to infty. We prove that there is a sharp phase transition at the threshold alphan∗=sqrt2logn: below the threshold (alphanllsqrt2logn), has negative components with probability tending to 1 while above (alphanggsqrt2logn), all the vector's components are eventually positive with probability tending to 1. At the critical scaling alphan∗, we provide a heuristics to evaluate the probability that 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 , that is when alphanggsqrt2logn, we establish that the Lotka-Volterra system of differential equations whose solution at equilibrium is precisely is stable in the sense that its jacobian {mathcal J}(�oldsymbol{x}_n) = mathrm{diag}(�oldsymbol{x}_n)left(-I_n + frac {A_n}{alpha_nsqrt{n}} ight) 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.












This page was built for publication: Positive solutions for large random linear systems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q4985368)