The Noise Collector for sparse recovery in high dimensions
The Noise Collector for sparse recovery in high dimensions
复制标题
DOI:
10.1073/pnas.1913995117
复制
发表时间:
2019-08
影响因子:
11.1
通讯作者:
M. Moscoso;A. Novikov;G. Papanicolaou;C. Tsogka
中科院分区:
文献类型:
--
作者:
M. Moscoso;A. Novikov;G. Papanicolaou;C. Tsogka
Significance The ability to detect sparse signals from noisy, high-dimensional data is a top priority in modern science and engineering. For optimal results, current approaches need to tune parameters that depend on the level of noise, which is often difficult to estimate. We develop a parameter-free, computationally efficient, ℓ1-norm minimization approach that has a zero false discovery rate (no false positives) with high probability for any level of noise while it detects the exact location of sparse signals when the noise is not too large. The ability to detect sparse signals from noisy, high-dimensional data is a top priority in modern science and engineering. It is well known that a sparse solution of the linear system Aρ=b0 can be found efficiently with an ℓ1-norm minimization approach if the data are noiseless. However, detection of the signal from data corrupted by noise is still a challenging problem as the solution depends, in general, on a regularization parameter with optimal value that is not easy to choose. We propose an efficient approach that does not require any parameter estimation. We introduce a no-phantom weight τ and the Noise Collector matrix C and solve an augmented system Aρ+Cη=b0+e, where e is the noise. We show that the ℓ1-norm minimal solution of this system has zero false discovery rate for any level of noise, with probability that tends to one as the dimension of b0 increases to infinity. We obtain exact support recovery if the noise is not too large and develop a fast Noise Collector algorithm, which makes the computational cost of solving the augmented system comparable with that of the original one. We demonstrate the effectiveness of the method in applications to passive array imaging.