Abstract
In the multireference alignment model, a signal is observed by the action of a random circular translation and the addition of Gaussian noise. The goal is to recover the signal's orbit by accessing multiple independent observations. Of particular interest is the sample complexity, i.e., the number of observations/samples needed in terms of the signal-to-noise ratio (SNR) (the signal energy divided by the noise variance) in order to drive the mean-square error to zero. Previous work showed that if the translations are drawn from the uniform distribution, then, in the low SNR regime, the sample complexity of the problem scales as ω (1/ SNR3). In this paper, using a generalization of the Chapman-Robbins bound for orbits and expansions of the χ2 divergence at low SNR, we show that in the same regime the sample complexity for any aperiodic translation distribution scales as ω (1/ SNR2). This rate is achieved by a simple spectral algorithm. We propose two additional algorithms based on non-convex optimization and expectation-maximization. We also draw a connection between the multireference alignment problem and the spiked covariance model.
Original language | English (US) |
---|---|
Article number | 8590822 |
Pages (from-to) | 3565-3584 |
Number of pages | 20 |
Journal | IEEE Transactions on Information Theory |
Volume | 65 |
Issue number | 6 |
DOIs | |
State | Published - Jun 2019 |
Bibliographical note
Publisher Copyright:© 1963-2012 IEEE.
Keywords
- Multireference alignment
- cryo-EM
- expectation-maximization
- method of moments
- non-convex optimization
- spectral algorithm
- spiked covariance model