Interpolation via weighted _1 minimization

From MaRDI portal
Publication:905907



Abstract: Functions of interest are often smooth and sparse in some sense, and both priors should be taken into account when interpolating sampled data. Classical linear interpolation methods are effective under strong regularity assumptions, but cannot incorporate nonlinear sparsity structure. At the same time, nonlinear methods such as l1 minimization can reconstruct sparse functions from very few samples, but do not necessarily encourage smoothness. Here we show that weighted l1 minimization effectively merges the two approaches, promoting both sparsity and smoothness in reconstruction. More precisely, we provide specific choices of weights in the l1 objective to achieve rates for functions with coefficient sequences in weighted lp spaces, p<=1. We consider the implications of these results for spherical harmonic and polynomial interpolation, in the univariate and multivariate setting. Along the way, we extend concepts from compressive sensing such as the restricted isometry property and null space property to accommodate weighted sparse expansions; these developments should be of independent interest in the study of structured sparse approximations and continuous-time compressive sensing problems.


In practice, one has often to interpolate functions which are smooth as well as sparse in some sense. In this excellent paper, the authors merge classical smoothness-based interpolation methods with modern sparsity constraints and nonlinear reconstruction methods. For a bounded domain \(D\), let \(\psi_j:\, D \to {\mathbb C}\) \((j \in \Lambda)\) be orthonormal functions with finite index set \(\Lambda\), \(|\Lambda| = N\). For given sampling points \(t_{\ell}\in D\) \((\ell = 1,\ldots,m)\) and \(f=\sum_{j\in \Lambda} x_j\,\psi_j\), let \(y = (f(t_{\ell}))_{\ell =1}^m\) and let \(A\) be the sampling matrix with the entries \(A_{\ell,j} = \psi_j(t_{\ell})\) \((\ell =1,\ldots,m;\,j\in \Lambda)\). For interpolation, the authors consider the function \(f^{\sharp} = \sum_{j\in \Lambda} x_j^{\sharp}\,\psi_j\), whose coefficient vector \(x^{\sharp}\) is the solution of the \textit{weighted} \(\ell_1\) \textit{minimization problem} \[ \min \|z\|_{\omega,1}\quad {\mathrm{subject to}}\quad \|Az - y\|_2 \leq \eta \] with the weighted \(\ell_1\) norm \(\|z\|_{\omega,1} = \sum_{j\in \Lambda} \omega_j\,|z_j|\) and convenient weights \(\omega_j \geq 1\). Using the new concepts of weighted null space property and weighted restricted isometry property of the sampling matrix \(A\), the authors prove general interpolation theorems. Corresponding error estimates of \(f -f^{\sharp}\) in \(L_{\infty}\) resp. \(L_2\) norm are given. In several examples and numerical tests, this theory is applied to spherical harmonic interpolation and tensorized polynomial interpolation (with Chebyshev resp. Legendre polynomials).



Cites work


Cited in
(70)








This page was built for publication: Interpolation via weighted \(\ell_{1}\) minimization

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q905907)