---
title: Fourier Sampling Numbers
url: https://www.emergentmind.com/topics/fourier-sampling-numbers
type: topic
---

# Fourier Sampling Numbers

Fourier sampling numbers quantify the minimal number of samples (typically in time or space) or measurements (typically in frequency) needed to reconstruct or approximate a signal or function from a given class to within prescribed accuracy. This concept is central to signal processing, harmonic analysis, approximation theory, and machine learning, particularly in determining optimal rates and strategies for function recovery under structural or smoothness constraints. The detailed characterization of Fourier sampling numbers varies according to the function class (e.g., bandlimited, Besov, Barron, or sparse), the reconstruction norm (e.g., $L^2$, $L^\infty$), the domain (continuous, periodic, discrete, or finite group), and the sampling model (uniform, non-uniform, adaptive, linear, or nonlinear).

## 1. Formal Definitions and General Frameworks

Let $\mathcal{C}$ be a class of functions (or signals) $f:\mathcal{D}\to\mathbb{C}$ with prescribed properties—e.g., bandlimited, $k$-Fourier-sparse, Besov smoothness, or membership in the Barron space. Given a family of permissible sampling functionals (time samples, Fourier coefficients, etc.), the *Fourier sampling number* $s_m(\mathcal{C}; X)$ or $s_N(\epsilon;\mathcal{C})$ is the minimal $m$ such that any $f\in\mathcal{C}$ can be reconstructed (in norm $X$) to error $\leq\epsilon$ from $m$ samples or measurements. Throughout, the underlying recovery can be constrained to linear/nonlinear, deterministic/randomized, uniform/non-uniform, or even spectrum-blind strategies.

For example, for bandlimited or $k$-Fourier-sparse signals on $[0,T]$, $m^*(k,\epsilon)$ denotes the smallest $m$ such that the signal can be interpolated from $m$ (possibly randomized) time-domain samples to $L_2$-error at most $\epsilon$ [2210.12495, 1812.08723].

## 2. Concrete Characterizations in Signal Classes

### Bandlimited, Multiband, and $k$-Sparse Classes

For classes characterized by Fourier support, the sampling number admits concrete expressions:
- **Bandlimited ($[-F,F]$ on $[0,T]$)**: $s_{[-F,F],\epsilon}=\Theta(FT + \log(1/\epsilon))$, matching classical Nyquist rates up to logarithms, with efficient algorithms approaching these rates [1812.08723].
- **Multiband**: $s_{MB,\epsilon}=O(\sum_i F_i T + s\log(1/\epsilon))$, matching the Landau rate for unions of $s$ bands.
- **$k$-Fourier-Sparse**: For signals with $k$ nonzero Fourier components, $s_{\textrm{sparse},\epsilon}=k$, and sample-optimal algorithms achieve $O(k\log(k/\delta))$ nonuniform samples [1812.08723].
- **Continuous $k$-Fourier-Sparse Interpolation**: Recent work established the optimality of quartic sample complexity: $m^*(k,\epsilon)=\Theta(k^4)$ for the minimal number of time samples needed for robust interpolation up to noise, closing the previous gap between information-theoretic lower and algorithmic upper bounds [2210.12495].

### Besov, Barron, and Decay-Smoothness Spaces

For broader smoothness spaces, the rates are determined by approximation theory:
- **Besov Spaces $B^s_\infty(L_q)$**: For periodic functions on $\mathbb{T}^d$, the Fourier sampling number $s_m(B^s_\infty(L_q),L_p)$ admits the decay
  \[
  s_m(B^s_\infty(L_q),L_p) \asymp m^{-s/d + (1/q - 1/p)_+}
  \]
  in the "linear" parameter regime ($q\geq2$ or $p\leq q$), with matching lower and upper bounds up to polylogarithmic factors for $s/d>1-1/p$. In the "nonlinear" regime ($1\leq q <2 < p$ or $1\le q < p \le 2$), additional logarithmic penalties arise [2508.13991]. The results specify nearly optimal Fourier measurement schemes (low-pass in the linear regime, hierarchical random subsampling in the nonlinear regime) and constructive convex recovery.
- **Fourier-Analytic Barron Space**: For $f(x)=\int_{\mathbb{R}^d} F(\xi) e^{2\pi i\langle x,\xi\rangle}\,d\xi$ with $\int |F(\xi)|(1+|\xi|)^\sigma\,d\xi<\infty$, the sampling number (optimal $m$ to reach error $\epsilon$ in $L_p$) scales as
  \[
  s_m(\sigma;L^p)\asymp m^{-(1/\max\{p,2\}+\sigma/d)}
  \]
  up to logarithmic factors, demonstrating dimension-independent algebraic rates for nonlinear sampling, in contrast to the classical curse of dimensionality for linear schemes [2208.07605].

### Finite Group and Discrete Settings

For functions on finite Abelian groups (e.g., $G = \mathbb{F}_p^r$), a universal sampling set $S\subseteq G$ is a subset such that any $t$-sparse Fourier signal can be reconstructed from its restriction to $S$. Explicit constructions partition $G$ into subspaces to guarantee "m-generating" properties, yielding
\[
s_{p,r}(t) \leq O(p t^2 r^2)
\]
with provable $\ell_1$-error bounds. No sub-quadratic (in $t$) lower bounds are known [1507.06849].

## 3. Quantitative Error Bounds and Sampling Rate Asymptotics

General results relate the Fourier sampling number to smoothness and decay parameters:
- **Continuous $\mathbb{R}$-Domain, Sobolev and Tail Conditions**:
  For $f$ with smoothness $m$ and spatial decay exponent $\alpha$, the root-mean-square error in DFT-based approaches satisfies
  \[
  E_{N,h}(f) \leq C_1 h^m + C_2 (N h)^{-\alpha}
  \]
  The minimax optimal $N$ for target error $\epsilon$ is
  \[
  N \gtrsim \epsilon^{-(m+\alpha)/m\alpha}
  \]
 – pure Sobolev ($\alpha\to\infty$): $N\sim\epsilon^{-1/m}$.
  – pure decay ($m=1$): $N\sim\epsilon^{-(1+\alpha)/\alpha}$.
  – both decay and smoothness: $N\sim\epsilon^{-(\alpha+\beta)/\alpha\beta}$ for additional frequency decay $\beta$ [2403.03810].

- **Uniform (ℓ$_\infty$) Error Bounds**:
  Recent NFFT-based theory gives explicit uniform bounds for any $f\in C(\mathbb{R})\cap L^1(\mathbb{R})$, e.g.,
  \[
  \sup_{|\omega|\leq\Omega_0} |\hat{f}(\omega) - S_{h,M}f(\omega)| \leq C_1 h^b + C_2 (M h)^{-(a-1)}
  \]
  for polynomial decay $|\hat{f}(v)|\leq d(1+|v|)^{-b}$, $|f(x)|\leq c(1+|x|)^{-a}$, with $N$ scaling as $O(\epsilon^{-1/b})$ or $O(\epsilon^{-(1 + (a-1)/b)})$ depending on the dominant term [2606.22970].

## 4. Algorithmic and Measurement Strategies

### Spectrum-Aware Sampling and Convex Recovery

- For large classes (bandlimited, sparse, Gaussian process), universal non-uniform sampling—via ridge leverage score sampling—achieves (up to log factors) the minimal possible sample complexity for all $\mu$ (the Fourier spectrum prior measure). The recovery is performed by kernel ridge regression in the continuous operator framework [1812.08723].
- Besov and Barron classes admit nearly optimal measurement strategies: low-pass frequency blocks for "linear" regimes, and scale-adapted random subsampling for "nonlinear" regimes, with reconstruction via convex minimization (e.g., $BV$ or Barron seminorm) [2508.13991, 2208.07605].

### Sparse and Finite Group Algorithms

Explicit universal sets and combinatorial partitioning underlie efficient $\ell_1$-bounded recovery for $t$-sparse signals on $G=\mathbb{F}_p^r$ [1507.06849].

## 5. Information-Theoretic Lower Bounds and Optimality Gaps

- For $k$-Fourier-sparse continuous interpolation, the information-theoretic lower bound for sample complexity is $\Omega(k^4)$, achieved algorithmically for the first time using high SNR bands and structural decomposition of signals [2210.12495].
- In general signal classes, information-theoretic lower bounds show that the sample complexity cannot be asymptotically improved beyond the statistical dimension $s_{\mu,\epsilon}$ up to logarithmic factors [1812.08723].
- For Besov norms with $q=1$ and certain $p>p_0\approx1.535$, an unavoidable polylogarithmic penalty separates the Fourier sampling numbers from Gelfand widths, demonstrating gaps in the nonlinear regime [2508.13991].

## 6. Practical Implications and Applications

- **Edge Recovery and Imaging**: In bounded-variation images ($BV$), well-designed Fourier sampling coupled with $BV$-minimization achieves edge localization rates $\sim 1/n$, outperforming grid-based sampling. Accurate reconstruction of sharp features is closely linked to the decay regime of the Fourier sampling number [2508.13991].
- **Efficient FFT-based Approximation**: The NFFT-based framework provides explicit recipes for choosing sampling and truncation parameters to ensure prescribed $\ell_\infty$ errors, with computational cost $O(N\log N)$ [2606.22970].
- **Universal Recovery in Signal Processing**: Spectrum-blind sampling designs match Nyquist and Landau rates for a broad array of models, extend to kriging and Gaussian process regression in one dimension, and provide robust guarantees under noise [1812.08723].

## 7. Comparative Table of Fourier Sampling Rates

| Function Class                  | Sampling Number $N(\epsilon)$       | Key Parameters                         | Reference       |
|----------------------------------|-------------------------------------|----------------------------------------|-----------------|
| $k$-Fourier-sparse (continuous)  | $O(k^4 \operatorname{polylog})$     | $k$: sparsity, $F$: bandlimit          | [2210.12495]    |
| Bandlimited $[-F,F]$             | $O(FT + \log(1/\epsilon))$          | $F$: max frequency, $T$: interval      | [1812.08723]    |
| Barron space $B^\sigma$          | $O(\epsilon^{-1/(1/\max\{p,2\}+\sigma/d)})$ | $\sigma$: smoothness, $d$: dim.        | [2208.07605]    |
| Besov $B^s_\infty(L_q)$          | $O(\epsilon^{-d/s})$ (linear regime) | $s$: smoothness, $d$: dim., $q,p$      | [2508.13991]    |
| Sobolev $W^{m,1}$                | $O(\epsilon^{-1/m})$                | $m$: smoothness order                  | [2403.03810]    |
| Polynomial decay (no smoothness) | $O(\epsilon^{-(1+\alpha)/\alpha})$  | $\alpha$: decay exponent               | [2403.03810]    |
| $t$-sparse (finite group)        | $O(p t^2 r^2)$                      | $p$: field size, $r$: dim., $t$        | [1507.06849]    |

These rates are optimal (up to polylogarithmic factors) under the stated assumptions.

## 8. Open Problems and Directions

- Removal of logarithmic factors (e.g., in Barron and Besov regimes), closing the remaining gaps to information-theoretic lower bounds [2208.07605, 2508.13991].
- Random versus deterministic sampling: it remains open whether simple random sampling can always achieve the optimal rates guaranteed by sophisticated leverage-score or combinatorial constructions [2208.07605].
- Precise constant dependence on dimension, structure of function class, and transition regions between linear and nonlinear measurement regimes.
- Extension of these results to function spaces beyond Barron/Besov/Sobolev, such as deep Barron-type or kernel-based classes, or to pointwise $L^\infty$ recovery in high dimension [2208.07605, 2508.13991].

---

References: [2210.12495], [1812.08723], [2208.07605], [1507.06849], [2403.03810], [2508.13991], [2606.22970]

Source: https://www.emergentmind.com/topics/fourier-sampling-numbers