Multivariate polynomial interpolation and sampling in Paley–Wiener spaces

In this paper, an equivalence between existence of particular exponential Riesz bases for spaces of multivariate bandlimited functions and existence of certain polynomial interpolants for functions in these spaces is given. Namely, polynomials are constructed which, in the limiting case, interpolate {(τn,f(τn))}n{(τn,f(τn))}n for certain classes of unequally spaced data nodes {τn}n{τn}n and corresponding ℓ2ℓ2 sampled data {f(τn)}n{f(τn)}n. Existence of these polynomials allows one to construct a simple sequence of approximants for an arbitrary multivariate bandlimited function ff which demonstrates L2L2 and uniform convergence on RdRd to ff. A simpler computational version of this recovery formula is also given at the cost of replacing L2L2 and uniform convergence on RdRd with L2L2 and uniform convergence on increasingly large subsets of RdRd. As a special case, the polynomial interpolants of given ℓ2ℓ2 data converge in the same fashion to the multivariate bandlimited interpolant of that same data. Concrete examples of pertinent Riesz bases and unequally spaced data nodes are also given.