Generalized Splines and Gaussian Processes
Abstract
For finite-dimensional linear inverse problems where the variables are Gaussian, it is well-known that the minimum-mean-square error estimator takes the form of a regularized least-squares data fit. In this chapter, we show that this equivalence extends to a much broader infinite-dimensional setting where generalized splines take the role of linear regressors and generalized Gaussian processes on a nuclear space S are the counterpart of Gaussian random vectors. The scope of this extension is of the same nature as the switch from the classic notion of function to that of a distribution, also known as a "generalized function." Our formalism involves a whitening/regularization operator L: S→ S' whose continuous extension induces a native Hilbert space H⊂ S' that plays a central role in our characterization. The presentation is self-contained for the most part and remarkably general and powerful. It allows for the recovery of all known instances of such equivalences; in particular, the methods involving innovations and reproducing-kernel Hilbert spaces developed by Kailath and his students, and the mathematical correspondence between fractional splines and Mandelbrot's fractional Brownian motion (fractals), with the former being the optimal estimators of the latter. It also covers general Bayesian methods for the resolution of infinite-dimensional inverse problems.
Related
- Variational Splines and Generalized Gaussian Processes
- Quadratic variations of rough generalized Hermite processes: hidden long memory and a universal Gaussian boundary field
- Quasihelix properties of selected Volterra Gaussian processes
- Gaussian approximation for the integrated intensity of counting processes: application to Hawkes processes
- UPPER TAIL PROBABILITY ESTIMATES OF SUPREMA OF GAUSSIAN PROCESSES BY MATRICES, A NUMBER THEORETIC EXAMPLE
BibTeX
@misc{unser2026generalized,
title = {{Generalized Splines and Gaussian Processes}},
author = {Michael Unser},
howpublished = {arXiv},
year = {2026},
month = aug,
doi = {10.48550/arxiv.2608.28446},
eprint = {2608.28446},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.28446},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Generalized Splines and Gaussian Processes
AU - Michael Unser
PY - 2026
DA - 2026/08/28
PB - arXiv
DO - 10.48550/arxiv.2608.28446
UR - https://arxiv.org/abs/2608.28446
AB - For finite-dimensional linear inverse problems where the variables are Gaussian, it is well-known that the minimum-mean-square error estimator takes the form of a regularized least-squares data fit. In this chapter, we show that this equivalence extends to a much broader infinite-dimensional setting where generalized splines take the role of linear regressors and generalized Gaussian processes on a nuclear space S are the counterpart of Gaussian random vectors. The scope of this extension is of the same nature as the switch from the classic notion of function to that of a distribution, also known as a "generalized function." Our formalism involves a whitening/regularization operator L: S→ S' whose continuous extension induces a native Hilbert space H⊂ S' that plays a central role in our characterization. The presentation is self-contained for the most part and remarkably general and powerful. It allows for the recovery of all known instances of such equivalences; in particular, the methods involving innovations and reproducing-kernel Hilbert spaces developed by Kailath and his students, and the mathematical correspondence between fractional splines and Mandelbrot's fractional Brownian motion (fractals), with the former being the optimal estimators of the latter. It also covers general Bayesian methods for the resolution of infinite-dimensional inverse problems.
ER -