A new stable basis for radial basis function interpolation

It is well known that radial basis function interpolants suffer from bad conditioning if the basis of translates is used. In the recent work by Pazouki and Schaback (2011), [5], the authors gave a quite general way to build stable and orthonormal bases for the native space NΦ(Ω)NΦ(Ω) associated to a kernel ΦΦ on a domain Ω⊂RsΩ⊂Rs. The method is simply based on the factorization of the corresponding kernel matrix. Starting from that setting, we describe a particular basis which turns out to be orthonormal in NΦ(Ω)NΦ(Ω) and in ℓ2,w(X)ℓ2,w(X), where XX is a set of data sites of the domain ΩΩ. The basis arises from a weighted singular value decomposition of the kernel matrix. This basis is also related to a discretization of the compact operator TΦ:NΦ(Ω)→NΦ(Ω)TΦ:NΦ(Ω)→NΦ(Ω), TΦ[f](x)=∫ΩΦ(x,y)f(y)dy∀x∈Ω, and provides a connection with the continuous basis that arises from an eigendecomposition of TΦTΦ. Finally, using the eigenvalues of this operator, we provide convergence estimates and stability bounds for interpolation and discrete least-squares approximation.