Table of Contents
Fetching ...

Point Processes and spatial statistics in time-frequency analysis

Barbara Pascal, Rémi Bardenet

TL;DR

This work reframes time-frequency analysis around the zeros of spectrograms, treating them as a spatial Point Process linked to Gaussian Analytic Functions (GAF). By mapping time-frequency representations to analytic functions via transforms like the Bargmann transform, the zeros acquire tractable probabilistic structure (stationarity, isotropy, hyperuniformity) that supports zero-based detection and denoising through spatial statistics. The authors develop explicit tools—from Edelman-Kostlan intensities to Ripley-type statistics and hole probabilities—and propose practical zero-based Monte Carlo tests and algorithms for signal detection. This framework offers principled, robust procedures for extracting signals from noisy measurements and opens multiple avenues for extending analytic-function-based time-frequency methods to broader classes of transforms and windows.

Abstract

A finite-energy signal is represented by a square-integrable, complex-valued function $t\mapsto s(t)$ of a real variable $t$, interpreted as time. Similarly, a noisy signal is represented by a random process. Time-frequency analysis, a subfield of signal processing, amounts to describing the temporal evolution of the frequency content of a signal. Loosely speaking, if $s$ is the audio recording of a musical piece, time-frequency analysis somehow consists in writing the musical score of the piece. Mathematically, the operation is performed through a transform $\mathcal{V}$, mapping $s \in L^2(\mathbb{R})$ onto a complex-valued function $\mathcal{V}s \in L^2(\mathbb{R}^2)$ of time $t$ and angular frequency $ω$. The squared modulus $(t, ω) \mapsto \vert\mathcal{V}s(t,ω)\vert^2$ of the time-frequency representation is known as the spectrogram of $s$; in the musical score analogy, a peaked spectrogram at $(t_0,ω_0)$ corresponds to a musical note at angular frequency $ω_0$ localized at time $t_0$. More generally, the intuition is that upper level sets of the spectrogram contain relevant information about in the original signal. Hence, many signal processing algorithms revolve around identifying maxima of the spectrogram. In contrast, zeros of the spectrogram indicate perfect silence, that is, a time at which a particular frequency is absent. Assimilating $\mathbb{R}^2$ to $\mathbb{C}$ through $z = ω+ \mathrm{i}t$, this chapter focuses on time-frequency transforms $\mathcal{V}$ that map signals to analytic functions. The zeros of the spectrogram of a noisy signal are then the zeros of a random analytic function, hence forming a Point Process in $\mathbb{C}$. This chapter is devoted to the study of these Point Processes, to their links with zeros of Gaussian Analytic Functions, and to designing signal detection and denoising algorithms using spatial statistics.

Point Processes and spatial statistics in time-frequency analysis

TL;DR

This work reframes time-frequency analysis around the zeros of spectrograms, treating them as a spatial Point Process linked to Gaussian Analytic Functions (GAF). By mapping time-frequency representations to analytic functions via transforms like the Bargmann transform, the zeros acquire tractable probabilistic structure (stationarity, isotropy, hyperuniformity) that supports zero-based detection and denoising through spatial statistics. The authors develop explicit tools—from Edelman-Kostlan intensities to Ripley-type statistics and hole probabilities—and propose practical zero-based Monte Carlo tests and algorithms for signal detection. This framework offers principled, robust procedures for extracting signals from noisy measurements and opens multiple avenues for extending analytic-function-based time-frequency methods to broader classes of transforms and windows.

Abstract

A finite-energy signal is represented by a square-integrable, complex-valued function of a real variable , interpreted as time. Similarly, a noisy signal is represented by a random process. Time-frequency analysis, a subfield of signal processing, amounts to describing the temporal evolution of the frequency content of a signal. Loosely speaking, if is the audio recording of a musical piece, time-frequency analysis somehow consists in writing the musical score of the piece. Mathematically, the operation is performed through a transform , mapping onto a complex-valued function of time and angular frequency . The squared modulus of the time-frequency representation is known as the spectrogram of ; in the musical score analogy, a peaked spectrogram at corresponds to a musical note at angular frequency localized at time . More generally, the intuition is that upper level sets of the spectrogram contain relevant information about in the original signal. Hence, many signal processing algorithms revolve around identifying maxima of the spectrogram. In contrast, zeros of the spectrogram indicate perfect silence, that is, a time at which a particular frequency is absent. Assimilating to through , this chapter focuses on time-frequency transforms that map signals to analytic functions. The zeros of the spectrogram of a noisy signal are then the zeros of a random analytic function, hence forming a Point Process in . This chapter is devoted to the study of these Point Processes, to their links with zeros of Gaussian Analytic Functions, and to designing signal detection and denoising algorithms using spatial statistics.
Paper Structure (36 sections, 25 theorems, 103 equations, 16 figures)

This paper contains 36 sections, 25 theorems, 103 equations, 16 figures.

Key Result

Proposition 1

Let the unit-energy Gaussian window of width $\sigma > 0$ be defined as Then, $g_{\sigma} \in L^1(\mathbb{R})\cap L^2(\mathbb{R})$ and its Fourier transform is also a Gaussian window whose standard deviation is the inverse of the original window's standard deviation.

Figures (16)

  • Figure 1: Examples of realistic elementary signals.
  • Figure 2: Spectrograms of elementary signals displayed in log-scale colormaps. Maximum values are in dark blue, zeros are represented as white dots.
  • Figure 3: Time-frequency analysis of pure noise.
  • Figure 4: Realizations of the signal-plus-noise model \ref{['eq:signalnoise']} with increasing noise levels.
  • Figure 5: The cross-disciplinarity of signal processing, pictured as the so-called golden triangle of signal processing of Patrick Flandrin, see for example flandrin2018explorations.
  • ...and 11 more figures

Theorems & Definitions (60)

  • Definition 1
  • Proposition 1
  • Proposition 2
  • Proposition 3
  • Proposition 4
  • Definition 2
  • Remark 1
  • Proposition 5
  • Remark 2
  • Proposition 6
  • ...and 50 more