Math 2.0

Information and communication theory, circuits Theorem

Landau's necessary density conditions for sampling

Builds on Band-limited functionFourier transformLp space

Statement

Let d≥1d \ge 1 and let S⊂ℝdS \subset \R^d be a compact set of positive measure whose boundary has Lebesgue measure zero. Let PWSPW_S be the space of band-limited functions with spectrum in SS: the f∈L2(ℝd)f \in L^2(\R^d) whose Fourier transform f^(ξ)=∫f(x)e−ix⋅ξ dx\hat f(\xi) = \int f(x) e^{-i x \cdot \xi} \, dx vanishes outside SS. Let Λ⊂ℝd\Lambda \subset \R^d be uniformly discrete, that is, inf⁡{∣λ−λ′∣:λ≠λ′∈Λ}>0\inf \set{\abs{\lambda - \lambda'} : \lambda \ne \lambda' \in \Lambda} > 0. Define the lower and upper Beurling densities D−(Λ)=lim inf⁡r→∞inf⁡x∈ℝd#(Λ∩(x+[0,r]d))rd,D+(Λ)=lim sup⁡r→∞sup⁡x∈ℝd#(Λ∩(x+[0,r]d))rd.D^-(\Lambda) = \liminf_{r \to \infty} \inf_{x \in \R^d} \frac{\#(\Lambda \cap (x + [0, r]^d))}{r^d}, \qquad D^+(\Lambda) = \limsup_{r \to \infty} \sup_{x \in \R^d} \frac{\#(\Lambda \cap (x + [0, r]^d))}{r^d}.

  1. If Λ\Lambda is a sampling set for PWSPW_S, that is, there are 0<A≤B0 < A \le B with A∥f∥22≤∑λ∈Λ∣f(λ)∣2≤B∥f∥22A \norm{f}_2^2 \le \sum_{\lambda \in \Lambda} \abs{f(\lambda)}^2 \le B \norm{f}_2^2 for all f∈PWSf \in PW_S, then D−(Λ)≥∣S∣/(2π)dD^-(\Lambda) \ge \abs{S} / (2\pi)^d.
  2. If Λ\Lambda is an interpolation set for PWSPW_S, that is, for every square-summable (aλ)λ∈Λ(a_\lambda)_{\lambda \in \Lambda} there is f∈PWSf \in PW_S with f(λ)=aλf(\lambda) = a_\lambda for all λ\lambda, then D+(Λ)≤∣S∣/(2π)dD^+(\Lambda) \le \abs{S} / (2\pi)^d.

Here ∣S∣\abs{S} is the Lebesgue measure of SS.

Remarks

Henry Landau proved these conditions in 1967. The proof compares the eigenvalues of two concentration operators: the operator that restricts f∈PWSf \in PW_S to a large cube, and its discrete analogue on the samples in that cube. The number of eigenvalues of the first operator that are not close to 00 or 11 is small compared with the volume of the cube. The main count is ∣S∣rd/(2π)d\abs{S} r^d / (2\pi)^d.

Example: for d=1d = 1 and S=[−π,π]S = [-\pi, \pi], the integers form a sampling and interpolation set, and D−(ℤ)=D+(ℤ)=1=∣S∣/2πD^-(\Z) = D^+(\Z) = 1 = \abs{S}/2\pi. This is the critical case of the Nyquist–Shannon sampling theorem. So uniform sampling at the Nyquist rate is optimal, and no irregular scheme can sample stably at a lower average rate.

Since f(λ)=(2π)−d∫Sf^(ξ)eiλ⋅ξ dξf(\lambda) = (2\pi)^{-d} \int_S \hat f(\xi) e^{i\lambda \cdot \xi} \, d\xi, the set Λ\Lambda is a sampling set if and only if the exponentials eiλ⋅ξe^{i\lambda \cdot \xi}, λ∈Λ\lambda \in \Lambda, form a frame of L2(S)L^2(S). It is a sampling and interpolation set if and only if they form a Riesz basis of L2(S)L^2(S), and then both parts give D−(Λ)=D+(Λ)=∣S∣/(2π)dD^-(\Lambda) = D^+(\Lambda) = \abs{S}/(2\pi)^d.

Going deeper

For a student

The Shannon sampling theorem gives a sufficient condition for one regular lattice. Landau's theorem gives a necessary condition for every uniformly discrete set, measured by the Beurling densities. It shows that the Nyquist rate is sharp also for irregular sampling. It also implies that a Riesz basis of exponentials for L2(S)L^2(S) must have uniform density exactly ∣S∣/(2π)d\abs{S}/(2\pi)^d, which links it with the theory of nonharmonic Fourier series. The theorem has analogues for Gabor systems and for other spaces of functions with localized transforms.

Reading path

To understand this page from the ground up, read these 4 pages in order. The foundations come first.

Show the path as a graph
Each page uses the pages above it. The colour of the bar on a box shows how strong the evidence is, and the box also gives it as text. The list has the same pages.
  • no evidence level
  1. Band-limited function concept
  2. Fourier transform concept
  3. Lp space concept
  4. Landau's necessary density conditions for sampling statement

Further reading

  1. Web https://en.wikipedia.org/wiki/Nyquist%E2%80%93Shannon_sampling_theorem
About this page
TagXWV6VS7BQM
URIm2://XWV6VS7BQM
Kindstatement, theorem
WritingW0 stub
MSC94A20, 42C15, 30D15
StatusThis copy of the site has no build data. Run m2 status XWV6VS7BQM.
AboutBand-limited function
UsesBand-limited function, Fourier transform, Lp space

Markdown · JSON · Edit on GitHub

What this page needs for W1

  • an inventory row: the row of the inventory that the node accounts for (0 of 1)
  • a source with a locator (0 of 1)

The writing levels say how complete the text of a page is. They do not measure its truth.

What links here

No page links to this page yet.