We study the problem of estimating a rank-$1$ signal in the presence of rotationally invariant noise-a class of perturbations more general than Gaussian noise. Principal Component Analysis (PCA) provides a natural estimator, and sharp results on its performance have been obtained in the high-dimensional regime. Recently, an Approximate Message Passing (AMP) algorithm has been proposed as an alternative estimator with the potential to improve the accuracy of PCA. However, the existing analysis of AMP requires an initialization that is both correlated with the signal and independent of the noise, which is often unrealistic in practice. In this work, we combine the two methods, and propose to initialize AMP with PCA. Our main result is a rigorous asymptotic characterization of the performance of this estimator. Both the AMP algorithm and its analysis differ from those previously derived in the Gaussian setting: at every iteration, our AMP algorithm requires a specific term to account for PCA initialization, while in the Gaussian case, PCA initialization affects only the first iteration of AMP. The proof is based on a two-phase artificial AMP that first approximates the PCA estimator and then mimics the true AMP. Our numerical simulations show an excellent agreement between AMP results and theoretical predictions, and suggest an interesting open direction on achieving Bayes-optimal performance.
翻译:我们研究的是,在交替变化中的噪音 -- -- 一种比高斯噪音更一般的扰动类型更一般的噪音 -- -- 扰动类别中,如何估计一至一美元的等级信号的问题。主构件分析(PCA)提供了一种自然估计器,其性能的精确结果已在高维系统中获得。最近,提出了一种近似信息传递算法,作为有可能提高五氯苯甲醚准确度的替代估计器。然而,对AMP的现有分析需要一种既与信号相关又独立于噪音的初始化,而这在实践中往往不切实际。在这项工作中,我们将两种方法结合起来,并提议将AMP与CPA初始化。我们的主要结果是对这一估计器的性能进行严格的非象征性的描述。AMP算法及其分析都与Gaussian环境先前推算出的结果不同:在每一次反复推算时,我们的AMP算法都需要一个明确的术语来核算五氯苯的初始初始化,而在戈斯案例中,CAPC的初始初始化只影响AMP的首次推算结果,而在AMAP的模拟阶段和初步推算出一个令人感兴趣的AMAPA- 之间的模拟结果。