Improved sampling accuracy for Fourier and Hadamard transforms
Restricted isometry of sampled Fourier and Hadamard matrices via entropic descent
Information Theory
Summary
Sampling data efficiently while keeping important information intact is a common problem in signal processing. The authors developed a new mathematical method that improves how well certain matrices preserve signal energy when only a subset of data is used. This method works not only for sparse signals, which have mostly zeros, but also for signals spread out more evenly. Their results apply to Fourier and Hadamard matrices, which are fundamental in many applications, enabling better and more stable recovery of signals from fewer measurements.
What this means in practice
- •For signal processing engineers: Improve the reconstruction of signals sampled under Fourier or Hadamard transforms using fewer measurements with controlled error.
- •For telecommunication system designers: Design more efficient sampling schemes for communication signals with guarantees on uniform energy preservation and noise stability.
Authors
W. Burstein, A. Iosevich, B. Krause
Abstract
In this second paper on descent methods as proof mechanisms in harmonic analysis (the first treated Bourgain's $Λ(p)$ selection theorem), entropic mirror descent gives an improved restricted isometry bound, improving sampling and recovery bounds in the Fourier Ratio program. Let $H$ be a complex $n\times n$ matrix with $H^*H=nI$ and $|H_{ij}|=1$. For fixed $0<\varepsilon<\frac14$ and $A_0>0$, independent Bernoulli row selection with expected cardinality $C_{\varepsilon,A_0}r\log\frac{2en}{r}\log(2r)\leq m\leq n/2$ preserves, after normalization by $\sqrt m$, the squared norm of every complex vector $y$ with $\|y\|_1\leq\sqrt r\|y\|_2$ within a factor $1\pm\varepsilon$, with failure probability at most $Cn^{-A_0}$. This includes every $r$-sparse vector and also fully supported ones; support size enters only through this inequality. For $y=\widehat f$ the condition is exactly $\operatorname{FR}(f)^2\leq r$, where $\operatorname{FR}(f)=\|\widehat f\|_1/\|\widehat f\|_2$, giving uniform energy sampling and approximate recovery even for signals with full Fourier support. The same holds for a uniform subset of prescribed cardinality, and for discrete Fourier and real Hadamard matrices. At fixed accuracy the bound removes one sparsity logarithm from the earlier count and replaces $\log n$ by $\log\frac{2en}{r}$; for Walsh matrices it matches, up to constants, the lower bound of Błasiok et al.\ in their range. One relative-entropy potential controls the corrections of an amplitude predictor at every scale; counting them on the symmetric difference of two samples gives the uniform estimate. We also prove the sufficient bound $C_{A_0}\varepsilon^{-5}r\log\frac{2en}{r}\log\frac{2r}{\varepsilon}$, and stable sparse recovery from $Cs\log\frac{2en}{s}\log(2s)$ measurements, with error controlled by the noise and the best $s$-term approximation.