Radially symmetric thin plate splines interpolating a circular contour map

Abstract

Profiles of radially symmetric thin plate spline surfaces minimizing the Beppo Levi energy over a compact annulus $R_{1}\leq r\leq R_{2}$ have been studied by Rabut via reproducing kernel methods. Motivated by our recent construction of Beppo Levi polyspline surfaces, we focus here on minimizing the radial energy over the full semi-axis $0<r<\infty$. Using a $L$-spline approach, we find two types of minimizing profiles: one is the limit of Rabut's solution as $R_{1}\rightarrow0$ and $R_{2}\rightarrow\infty$ (identified as a `non-singular' $L$-spline), the other has a second-derivative singularity and matches an extra data value at $0$. For both profiles and $p\in\left[ 2,\infty\right]$, we establish the $L^{p}$-approximation order $3/2+1/p$ in the radial energy space. We also include numerical examples and obtain a novel representation of the minimizers in terms of dilates of a basis function.Comment: new figures and sub-sections; new Proposition 1 replacing old Corollary 1; shorter proof of Theorem 4; one new referenc