Rank-One Signal Recovery in Sparse Wishart Noise

We study the high-dimensional recovery of a signal vector $\mathbf{x}$ in the presence of sparse Wishart-like noise. We define an $N\times N$ matrix $A=J+(\theta /N){\mathbf{xx}}^{\top }$, where ${\mathbf{xx}}^{\top }$ is the rank-one deformation of the random noise matrix $J$. We consider a Wishart-like matrix $J=X^{\top }X$, where $X$ is a sparse $M\times N$ random matrix with entries $X_{ij}=c_{ij}{W}_{ij}$, with $c_{ij}$ regulating the density of non-zero elements, and $W_{ij}$ the bond weights. Using the replica method, we compute analytically the top eigenpair statistics of $A$, and their dependence on the signal strength $\theta$, the rectangularity ratio $\alpha =\sqrt{M/N}$, and the average connectivity of the noise. The spectral observables are expressed in terms of a system of Recursive Distributional Equations, which are efficiently solved via a Population Dynamics algorithm. They allow us to compute the average largest eigenvalue $⟨{\lambda }_1{⟩}_A$, the average top eigenvector component density, and the average overlap between the top eigenvector of $A$ and $\mathbf{x}$. We identify a critical threshold ${\theta }_{crit}$--depending on the average connectivity of the noise--that marks a BBP-like phase transition: below this value, $⟨{\lambda }_1{⟩}_A$ is unaffected by the signal, and the overlap vanishes. Thus, the signal is not recoverable from the top eigenvector of $A$. For $\theta >{\theta }_{crit}$, the signal-related outlier eigenvalue becomes $⟨{\lambda }_1{⟩}_A$ and the overlap is nonzero, allowing for recovery of the signal. The results are in excellent agreement with numerical diagonalisation. We show that in the dense limit, the recovery threshold and eigen-statistics converge to the results predicted by the classical BBP transition for additive rank-one deformations of dense Wishart matrices.