Repository logo

Infoscience

  • English
  • French
Log In
Logo EPFL, École polytechnique fédérale de Lausanne

Infoscience

  • English
  • French
Log In
  1. Home
  2. Academic and Research Output
  3. Journal articles
  4. Generalized Smoothing Splines and the Optimal Discretization of the Wiener Filter
 
Loading...
Thumbnail Image
research article

Generalized Smoothing Splines and the Optimal Discretization of the Wiener Filter

Unser, M.  
•
Blu, T.  
2005
IEEE Transactions on Signal Processing

We introduce an extended class of cardinal LL-splines, where L is a pseudo-differential operator satisfying some admissibility conditions. We show that the LL-spline signal interpolation problem is well posed and that its solution is the unique minimizer of the spline energy functional, $\parallel L s \parallel _{ L _{ 2 } } ^{ 2 } $ , subject to the interpolation constraint. Next, we consider the corresponding regularized least squares estimation problem, which is more appropriate for dealing with noisy data. The criterion to be minimized is the sum of a quadratic data term, which forces the solution to be close to the input samples, and a “smoothness” term that privileges solutions with small spline energies. Here, too, we find that the optimal solution, among all possible functions, is a cardinal LL-spline. We show that this smoothing spline estimator has a stable representation in a B-spline-like basis and that its coefficients can be computed by digital filtering of the input signal. We describe an efficient recursive filtering algorithm that is applicable whenever the transfer function of L is rational (which corresponds to the case of exponential splines). We justify these algorithms statistically by establishing an equivalence between LL smoothing splines and the minimum mean square error (MMSE) estimation of a stationary signal corrupted by white Gaussian noise. In this model-based formulation, the optimum operator L is the whitening filter of the process, and the regularization parameter is proportional to the noise variance. Thus, the proposed formalism yields the optimal discretization of the classical Wiener filter, together with a fast recursive algorithm. It extends the standard Wiener solution by providing the optimal interpolation space. We also present a Bayesian interpretation of the algorithm.

  • Files
  • Details
  • Metrics
Loading...
Thumbnail Image
Name

unser0506.pdf

Access type

openaccess

Size

486.5 KB

Format

Adobe PDF

Checksum (MD5)

355366d25fc7bb2a780eb29c0e9a3ada

Logo EPFL, École polytechnique fédérale de Lausanne
  • Contact
  • infoscience@epfl.ch

  • Follow us on Facebook
  • Follow us on Instagram
  • Follow us on LinkedIn
  • Follow us on X
  • Follow us on Youtube
AccessibilityLegal noticePrivacy policyCookie settingsEnd User AgreementGet helpFeedback

Infoscience is a service managed and provided by the Library and IT Services of EPFL. © EPFL, tous droits réservés