共查询到20条相似文献,搜索用时 203 毫秒
1.
Multistep interpolation of scattered data by compactly supported radial basis functions requires hierarchical subsets of the
data. This paper analyzes thinning algorithms for generating evenly distributed subsets of scattered data in a given domain
in ℝ
d
. 相似文献
2.
The paper obtains error estimates for approximation by radial basis functions on the sphere. The approximations are generated
by interpolation at scattered points on the sphere. The estimate is given in terms of the appropriate power of the fill distance
for the interpolation points, in a similar manner to the estimates for interpolation in Euclidean space. A fundamental ingredient
of our work is an estimate for the Lebesgue constant associated with certain interpolation processes by spherical harmonics.
These interpolation processes take place in ``spherical caps' whose size is controlled by the fill distance, and the important
aim is to keep the relevant Lebesgue constant bounded. This result seems to us to be of independent interest.
March 27, 1997. Dates revised: March 19, 1998; August 5, 1999. Date accepted: December 15, 1999. 相似文献
3.
4.
Quasi-interpolation of radial basis functions on finite grids is a very useful strategy in approximation theory and its applications. A notable strongpoint of the strategy is to obtain directly the approximants without the need to solve any linear system of equations. For radial basis functions with Gaussian kernel, there have been more studies on the interpolation and quasi-interpolation on infinite grids. This paper investigates the approximation by quasi-interpolation operators with Gaussian kernel on the compact interval. The approximation errors for two classes of function with compact support sets are estimated. Furthermore, the approximation errors of derivatives of the approximants to the corresponding derivatives of the approximated functions are estimated. Finally, the numerical experiments are presented to confirm the accuracy of the approximations. 相似文献
5.
In this paper, we describe a recursive method for computing interpolants defined in a space spanned by a finite number of continuous functions in Rd. We apply this method to construct several interpolants such as spline interpolants, tensor product interpolants and multivariate polynomial interpolants. We also give a simple algorithm for solving a multivariate polynomial interpolation problem and constructing the minimal interpolation space for a given finite set of interpolation points. 相似文献
6.
Ghislain Franssens 《Advances in Computational Mathematics》1999,10(3-4):367-388
A new C
∞ interpolant is presented for the univariate Hermite interpolation problem. It differs from the classical solution in that
the interpolant is of non‐polynomial nature. Its basis functions are a set of simple, compact support, transcendental functions.
The interpolant can be regarded as a truncated Multipoint Taylor series. It has essential singularities at the sample points,
but is well behaved over the real axis and satisfies the given functional data. The interpolant converges to the underlying
real‐analytic function when (i) the number of derivatives at each point tends to infinity and the number of sample points
remains finite, and when (ii) the spacing between sample points tends to zero and the number of specified derivatives at each
sample point remains finite.
A comparison is made between the numerical results achieved with the new method and those obtained with polynomial Hermite
interpolation. In contrast with the classical polynomial solution, the new interpolant does not suffer from any ill conditioning,
so it is always numerically stable. In addition, it is a much more computationally efficient method than the polynomial approach.
This revised version was published online in June 2006 with corrections to the Cover Date. 相似文献
7.
A new class of radial basis functions with compact support 总被引:1,自引:0,他引:1
M. D. Buhmann. 《Mathematics of Computation》2001,70(233):307-318
Radial basis functions are well-known and successful tools for the interpolation of data in many dimensions. Several radial basis functions of compact support that give rise to nonsingular interpolation problems have been proposed, and in this paper we study a new, larger class of smooth radial functions of compact support which contains other compactly supported ones that were proposed earlier in the literature.
8.
In this paper, we consider multivariate interpolation with radial basis functions of finite smoothness. In particular, we
show that interpolants by radial basis functions in ℝ
d
with finite smoothness of even order converge to a polyharmonic spline interpolant as the scale parameter of the radial basis
functions goes to zero, i.e., the radial basis functions become increasingly flat. 相似文献
9.
For the solution of large sparse linear systems arising from interpolation problems using compactly supported radial basis
functions, a class of efficient numerical algorithms is presented. They iteratively select small subsets of the interpolation
points and refine the current approximative solution there. Convergence turns out to be linear, and the technique can be generalized
to positive definite linear systems in general. A major feature is that the approximations tend to have only a small number
of nonzero coefficients, and in this sense the technique is related to greedy algorithms and best n-term approximation.
This revised version was published online in June 2006 with corrections to the Cover Date. 相似文献
10.
The Padua points are a family of points on the square [−1, 1]2 given by explicit formulas that admits unique Lagrange interpolation by bivariate polynomials. Interpolation polynomials
and cubature formulas based on the Padua points are studied from an ideal theoretic point of view, which leads to the discovery
of a compact formula for the interpolation polynomials. The L
p
convergence of the interpolation polynomials is also studied.
S. De Marchi and M. Vianello were supported by the “ex-60%” funds of the University of Padua and by the INdAM GNCS (Italian
National Group for Scientific Computing). Y. Xu was partially supported by NSF Grant DMS-0604056. 相似文献
11.
We consider interpolation methods defined by positive definite functions on a locally compact group G. Estimates for the smallest and largest eigenvalue of the interpolation matrix in terms of the localization of the positive definite function on G are presented, and we provide a method to get positive definite functions explicitly on compact semisimple Lie groups. Finally, we apply our results to construct well-localized positive definite basis functions having nice stability properties on the rotation group SO(3). 相似文献
12.
C. M. C. Roque A. J. M. Ferreira 《Numerical Methods for Partial Differential Equations》2010,26(3):675-689
A numerical investigation on a technique for choosing an optimal shape parameter is proposed. Radial basis functions (RBFs) and their derivatives are used as interpolants in the asymmetric collocation radial basis method, for solving systems of partial differential equations. The shape parameter c in RBFs plays a major role in obtaining high quality solutions for boundary value problems. As c is a user defined value, inexperienced users may compromise the quality of the solution, often a problem of this meshless method. Here we propose a statistical technique to choose the shape parameter in radial basis functions. We use a cross‐validation technique suggested by Rippa 6 for interpolation problems to find a cost function Cost(c) that ideally has the same behavior as an error function. If that is the case, the parameter c that minimizes the cost function will be an optimal shape parameter, in the sense that it minimizes the error function. The form of the cost and error functions are analized for several examples, and for most cases the two functions have a similar behavior. The technique produced very accurate results, even with a small number of points and irregular grids. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010 相似文献
13.
John P. Boyd 《Journal of Computational and Applied Mathematics》2010,234(5):1435-1441
Radial basis function (RBF) interpolation is a “meshless” strategy with great promise for adaptive approximation. One restriction is “error saturation” which occurs for many types of RBFs including Gaussian RBFs of the form ?(x;α,h)=exp(−α2(x/h)2): in the limit h→0 for fixed α, the error does not converge to zero, but rather to ES(α). Previous studies have theoretically determined the saturation error for Gaussian RBF on an infinite, uniform interval and for the same with a single point omitted. (The gap enormously increases ES(α).) We show experimentally that the saturation error on the unit interval, x∈[−1,1], is about 0.06exp(−0.47/α2)‖f‖∞ — huge compared to the O(2π/α2)exp(−π2/[4α2]) saturation error for a grid with one point omitted. We show that the reason the saturation is so large on a finite interval is that it is equivalent to an infinite grid which is uniform except for a gap of many points. The saturation error can be avoided by choosing α?1, the “flat limit”, but the condition number of the interpolation matrix explodes as O(exp(π2/[4α2])). The best strategy is to choose the largest α which yields an acceptably small saturation error: If the user chooses an error tolerance δ, then . 相似文献
14.
Yasin Fadaei Zareen A. Khan Ali Akgül 《Mathematical Methods in the Applied Sciences》2019,42(16):5595-5606
A greedy algorithm in combination with radial basis functions partition of unity collocation (GRBF‐PUC) scheme is used as a locally meshless method for American option pricing. The radial basis function partition of unity method (RBF‐PUM) is a localization technique. Because of having interpolation matrices with large condition numbers, global approximants and some local ones suffer from instability. To overcome this, a greedy algorithm is added to RBF‐PUM. The greedy algorithm furnishes a subset of best nodes among the points X. Such nodes are then used as points of trial in a locally supported RBF approximant for each partition. Using of greedy selected points leads to decreasing the condition number of interpolation matrices and reducing the burdensome in pricing American options. 相似文献
15.
The solutions of the Carathéodory–Fejér interpolation problem for generalized Schur functions can be parametrized via a linear fractional transformation over the class of classical Schur functions. The linear fractional transformation of some of these functions may have a pole (simple or multiple) in one or more of the interpolation points or not satisfy one or more interpolation conditions, hence not all Schur functions can serve as a parameter. The set of excluded parameters is characterized in terms of the related Pick matrix.Research was supported by the Summer Research Grant from the College of William and MarySubmitted: June 26, 2002 Revised: January 31, 2003 相似文献
16.
J. M. Carnicer 《Numerical Algorithms》1992,3(1):125-132
In this paper the necessary and sufficient conditions for given data to admit a rational interpolant in
k,1 with no poles in the convex hull of the interpolation points is studied. A method for computing the interpolant is also provided.Partially supported by DGICYT-0121. 相似文献
17.
A non-uniform, variational refinement scheme is presented for computing piecewise linear curves that minimize a certain discrete energy functional subject to convex constraints on the error from interpolation. Optimality conditions are derived for both the fixed and free-knot problems. These conditions are expressed in terms of jumps in certain (discrete) derivatives. A computational algorithm is given that applies to constraints whose boundaries are either piecewise linear or spherical. The results are applied to closed periodic curves, open curves with various boundary conditions, and (approximate) Hermite interpolation. 相似文献
18.
19.
The ordered pair (T,I) of two self-maps of a metric space (X,d) is called a Banach operator pair if the set F(I) of fixed points of I is T-invariant i.e. T(F(I))⊆F(I). Some common fixed point theorems for a Banach operator pair and the existence of common fixed points of best approximation are presented in this paper. The results prove, generalize and extend some results of Al-Thagafi [M.A. Al-Thagafi, Common fixed points and best approximation, J. Approx. Theory 85 (1996) 318-323], Carbone [A. Carbone, Applications of fixed point theorems, Jnanabha 19 (1989) 149-155], Chen and Li [J. Chen, Z. Li, Common fixed points for Banach operator pairs in best approximations, J. Math. Anal. Appl. 336 (2007) 1466-1475], Habiniak [L. Habiniak, Fixed point theorems and invariant approximation, J. Approx. Theory 56 (1989) 241-244], Jungck and Sessa [G. Jungck, S. Sessa, Fixed point theorems in best approximation theory, Math. Japon. 42 (1995) 249-252], Sahab, Khan and Sessa [S.A. Sahab, M.S. Khan, S. Sessa, A result in best approximation theory, J. Approx. Theory 55 (1988) 349-351], Shahzad [N. Shahzad, Invariant approximations and R-subweakly commuting maps, J. Math. Anal. Appl. 257 (2001) 39-45] and of few others. 相似文献
20.
Robert M. Corless Nargol Rezvani Amirhossein Amiraslani 《Mathematics in Computer Science》2007,1(2):353-374
Spectra and pseudospectra of matrix polynomials are of interest in geometric intersection problems, vibration problems, and
analysis of dynamical systems. In this note we consider the effect of the choice of polynomial basis on the pseudospectrum
and on the conditioning of the spectrum of regular matrix polynomials. In particular, we consider the direct use of the Lagrange
basis on distinct interpolation nodes, and give a geometric characterization of “good” nodes. We also give some tools for
computation of roots at infinity via a new, natural, reversal. The principal achievement of the paper is to connect pseudospectra
to the well-established theory of Lebesgue functions and Lebesgue constants, by separating the influence of the scalar basis
from the natural scale of the matrix polynomial, which allows many results from interpolation theory to be applied.
This work was partially funded by the Natural Sciences and Engineering Research Council of Canada, and by the MITACS Network
of Centres of Excellence. 相似文献