Statistical Mechanics of MAP Estimation: General Replica Ansatz

Bereyhi A, Müller R, Schulz-Baldes H (2019)


Publication Type: Journal article

Publication year: 2019

Journal

Book Volume: 65

Pages Range: 7896-7934

Article Number: 8812691

Journal Issue: 12

DOI: 10.1109/TIT.2019.2937335

Abstract

The large-system performance of maximum-a-poste-rior estimation is studied considering a general distortion function when the observation vector is received through a linear system with additive white Gaussian noise. The analysis considers the system matrix to be chosen from the large class of rotationally invariant random matrices. We take a statistical mechanical approach by introducing a spin glass corresponding to the estimator, and employing the replica method for the large-system analysis. In contrast to earlier replica based studies, our analysis evaluates the general replica ansatz of the corresponding spin glass and determines the asymptotic distortion of the estimator for any structure of the replica correlation matrix. Consequently, the replica symmetric as well as the replica symmetry breaking ansatz with b steps of breaking is deduced from the given general replica ansatz. The generality of our distortion function lets us derive a more general form of the maximum-a-posterior decoupling principle. Based on the general replica ansatz, we show that for any structure of the replica correlation matrix, the vector-valued system decouples into a bank of equivalent decoupled scalar systems followed by maximum-a-posterior estimators. The structure of the decoupled system is further studied under both the replica symmetry and the replica symmetry breaking assumptions. For b steps of symmetry breaking, the decoupled system is found to be an additive system with a non-Gaussian noise term given as the sum of an independent Gaussian random variable with b non-Gaussian impairment terms which depend on the input symbol. The general decoupling property of the maximum-a-posterior estimator leads to the idea of a replica simulator which represents the replica ansatz through the state evolution of a transition system described by its corresponding decoupled system. As an application of our study, we investigate large compressive sensing systems by considering the ℓp norm minimization recovery schemes. Our numerical investigations show that the replica symmetric ansatz for ℓ0 norm recovery fails to give an accurate approximation of the mean square error as the compression rate grows, and therefore, the replica symmetry breaking ansätze are needed in order to assess the performance precisely.

Authors with CRIS profile

How to cite

APA:

Bereyhi, A., Müller, R., & Schulz-Baldes, H. (2019). Statistical Mechanics of MAP Estimation: General Replica Ansatz. IEEE Transactions on Information Theory, 65(12), 7896-7934. https://dx.doi.org/10.1109/TIT.2019.2937335

MLA:

Bereyhi, Ali, Ralf Müller, and Hermann Schulz-Baldes. "Statistical Mechanics of MAP Estimation: General Replica Ansatz." IEEE Transactions on Information Theory 65.12 (2019): 7896-7934.

BibTeX: Download