Variational Seasonal-Trend Decomposition with Sparse Continuous-Domain Regularization
Julien Fageot
TL;DR
This work addresses recovering a continuous-domain signal $f_0=f_{0,T}+f_{0,S}$ from finite noisy measurements by formulating a variational problem with generalized total-variation regularization that promotes sparse, spline-like representations for both the trend $f_{0,T}$ and the 1-periodic seasonal part $f_{0,S}$. It develops a functional-analytic foundation, defines seasonal-trend native spaces, and proves a representer theorem guaranteeing composite spline solutions for both components, along with a $\Gamma$-convergence analysis showing grid-based discretizations converge uniformly to the continuous problem. The results establish, under broad conditions, that minimizers are sparse splines with explicitly characterized forms, and that practical grid-based algorithms reliably approximate the true continuous solutions. Overall, the paper provides a rigorous, scalable approach to seasonal-trend decomposition in noisy, limited-measurement regimes, with strong theoretical guarantees for both representation and discretization.
Abstract
We consider the inverse problem of recovering a continuous-domain function from a finite number of noisy linear measurements. The unknown signal is modeled as the sum of a slowly varying trend and a periodic or quasi-periodic seasonal component. We formulate a variational framework for their joint recovery by introducing convex regularizations based on generalized total variation, which promote sparsity in spline-like representations. Our analysis is conducted in an infinite-dimensional setting and leads to a representer theorem showing that minimizers are splines in both components. To make the approach numerically feasible, we introduce a family of discrete approximations and prove their convergence to the original problem in the sense of $Γ$-convergence. This further ensures the uniform convergence of discrete solutions to their continuous counterparts. The proposed framework offers a principled approach to seasonal-trend decomposition in the presence of noise and limited measurements, with theoretical guarantees on both representation and discretization.
