Xingping Sun

Alexander Kushpel, Jeremy Levesley, Xingping Sun

In this paper, we study the problem of interpolating a continuous function at $(n+1)$ equally-spaced points in the interval $[0,1]$, using shifts of a kernel on the $(1/n)$-spaced infinite grid. The archetypal example here is approximation using shifts of a Gaussian kernel. We present new results concerning interpolation of functions of exponential type, in particular, polynomials on the integer grid as a step en route to solve the general interpolation problem. For the Gaussian kernel we introduce a new class of polynomials, closely related to the probabilistic Hermite polynomials and show that evaluations of the polynomials at the integer points provide the coefficients of the interpolants. Taking cue from the classical Newton polynomial interpolation, we derive a closed formula for the Gaussian interpolant of a continuous function on a uniform grid in the unit interval.

Shao-Bo Lin, Xingping Sun, Di Wang

For radial basis function (RBF) kernel interpolation of scattered data, Schaback in 1995 proved that the attainable approximation error and the condition number of the underlying interpolation matrix cannot be made small simultaneously. He referred to this finding as an "uncertainty relation", an undesirable consequence of which is that RBF kernel interpolation is susceptible to noisy data. In this paper, we propose and study a distributed interpolation method to manage and quantify the uncertainty brought on by interpolating noisy spherical data of non-negligible magnitude. We also present numerical simulation results showing that our method is practical and robust in terms of handling noisy data from challenging computing environments.

Shao-Bo Lin, Xiangyu Chang, Xingping Sun

Data sites selected from modeling high-dimensional problems often appear scattered in non-paternalistic ways. Except for sporadic clustering at some spots, they become relatively far apart as the dimension of the ambient space grows. These features defy any theoretical treatment that requires local or global quasi-uniformity of distribution of data sites. Incorporating a recently-developed application of integral operator theory in machine learning, we propose and study in the current article a new framework to analyze kernel interpolation of high dimensional data, which features bounding stochastic approximation error by the spectrum of the underlying kernel matrix. Both theoretical analysis and numerical simulations show that spectra of kernel matrices are reliable and stable barometers for gauging the performance of kernel-interpolation methods for high dimensional data.