Orthonormal Expansions for Translation-Invariant Kernels
Filip Tronarp, Toni Karvonen
TL;DR
This paper introduces a general Fourier-analytic framework for constructing orthonormal bases of translation-invariant kernels on $\mathbb{R}$ from bases of $\mathscr{L}_2(\mathbb{R})$, yielding pointwise-convergent RKHS expansions $r(t,u)=\sum_m\psi_m^*(t)\psi_m(u)$. It applies the method to three canonical kernels: Matérn kernels with half-integer order via Laguerre-based bases (including a finite null-space portion), the Cauchy kernel via complex/real Cauchy–Laguerre bases, and the Gaussian kernel via Hermite functions; each yields explicit, computable expansions and tractable truncation errors. The Gaussian expansion recovers a Mercer expansion consistent with Mehler’s formula, while the Matérn expansions are nearly Mercer, up to a finite null-space. The results provide practical, scalable representations of these kernels for numerical methods and theoretical analysis, with clear pathways for extending to other translations-invariant kernels and scaling regimes.
Abstract
We present a general Fourier analytic technique for constructing orthonormal basis expansions of translation-invariant kernels from orthonormal bases of $\mathscr{L}_2(\mathbb{R})$. This allows us to derive explicit expansions on the real line for (i) Matérn kernels of all half-integer orders in terms of associated Laguerre functions, (ii) the Cauchy kernel in terms of rational functions, and (iii) the Gaussian kernel in terms of Hermite functions.
