Polynomial interpolation of minimal degree

Summary. Minimal degree interpolation spaces with respect to a finite set of points are subspaces of multivariate polynomials of least possible degree for which Lagrange interpolation with respect to the given points is uniquely solvable and degree reducing. This is a generalization of the concept of least interpolation introduced by de Boor and Ron. This paper investigates the behavior of Lagrange interpolation with respect to these spaces, giving a Newton interpolation method and a remainder formula for the error of interpolation. Moreover, a special minimal degree interpolation space will be introduced which is particularly beneficial from the numerical point of view. Received June 9, 1995 / Revised version received June 26, 1996     

Interpolation-based nonlinear parametric MOR for gas pipelines

This work deals with the model order reduction (MOR) of nonlinear, parametric systems of partial differential equations as they arise in gas pipeline modeling. We present an approach that is based on a linearization around parametric working points, linear modal reduction and interpolation. The choice of the working points as well as the interpolation strategy crucially determine the approximation quality. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convergence of interpolation in polynomial Chebyshev approximation

The quality of a polynomial approximation on an interval to a functionf is considered as a function of its points of interpolation. Iff satisfies a Lipschitz condition of order 1, the quality depends linearly on the distance of the points of interpolation from an optimal interpolating point set: further restrictions onf still give only linear dependence. This suggests that algorithms based on interpolation are inferior to algorithms based on error extrema (such as the Remes algorithm).     

Are the eigenvalues of preconditioned banded symmetric Toeplitz matrices known in almost closed form?

The problem of computing oscillatory integrals with general oscillators is considered. We employ a Filon-type method, where the interpolation basis functions are chosen in such a way that the moments are in terms of elementary functions and the oscillator only. This allows us to evaluate the moments rapidly and easily without needing to engage hypergeometric functions. The proposed basis functions form a Chebyshev set for any oscillator function even if it has some stationary points in the integration interval. This property enables us to employ the Filon-type method without needing any information about the stationary points if any. Interpolation by the proposed basis functions at the Fekete points (which are known as nearly optimal interpolation points), when combined with the idea of splines, leads to a reliable convergent method for computing the oscillatory integrals. Our numerical experiments show that the proposed method is more efficient than the earlier ones with the same advantages.     

Block Based Bivariate Blending Rational Interpolation via Symmetric Branched Continued Fractions

This paper constructs a new kind of block based bivariate blending rational interpolation via symmetric branched continued fractions. The construction process may be outlined as follows. The first step is to divide the original set of support points into some subsets (blocks). Then construct each block by using symmetric branched continued fraction. Finally assemble these blocks by Newton’s method to shape the whole interpolation scheme. Our new method offers many flexible bivariate blending rational interpolation schemes which include the classical bivariate Newton’s polynomial interpolation and symmetric branched continued fraction interpolation as its special cases. The block based bivariate blending rational interpolation is in fact a kind of tradeoff between the purely linear interpolation and the purely nonlinear interpolation. Finally, numerical examples are given to show the effectiveness of the proposed method.     

修正的 Thiele-Werner型有理插值

Through adjusting the order of interpolation nodes, we gave a kind of modified Thiele-Werner rational interpolation. This interpolation method not only avoids the infinite value of inverse differences in constructing the Thiele continued fraction interpolation, but also simplifies the interpolating polynomial coefficients with constant coefficients in the Thiele-Werner rational interpolation. Unattainable points and determinantal expression for this interpolation are considered. As an extension, some bivariate analogy is also discussed and numerical examples are given to show the validness of this method.     

Interpolation by Polynomials and Radial Basis Functions on Spheres

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.     

Multigrid technique for biharmonic interpolation with application to dual and multiple reciprocity method

A biharmonic-type interpolation method is presented to solve 2D and 3D scattered data interpolation problems. Unlike the methods based on radial basis functions, which produce a large linear system of equations with fully populated and often non-selfadjoint and ill-conditioned matrix, the presented method converts the interpolation problem to the solution of the biharmonic equation supplied with some non-usual boundary conditions at the interpolation points. To solve the biharmonic equation, fast multigrid techniques can be applied which are based on a non-uniform, non-equidistant but Cartesian grid generated by the quadtree/octtree algorithm. The biharmonic interpolation technique is applied to the multiple and dual reciprocity method of the BEM to convert domain integrals to the boundary. This makes it possible to significantly reduce the computational cost of the evaluation of the appearing domain integrals as well as the memory requirement of the procedure. The resulting method can be considered as a special grid-free technique, since it requires no domain discretisation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Extraction of characteristic points and its fractal reconstruction for terrain profile data

The main objective of this paper is to study reduction rate of 2D DEM (digital elevation model) data profile after data reduction by the Douglas–Peucker (DP) linear simplification method and by fractal interpolation to show original terrain reconstruction. In this paper, two-dimensional data of measured geographic profiles are taken as the study object, by using the DP method and the improved Douglas–Peucker (IDP) method to reduce data. Its aim is to retain spatial linear characteristics and variations, then take reduced data points as basic points and use the random fractal interpolation approach to add more data points up to the same as the original data points, in order to reconstruct the terrain, and compare the experimental data with the random point extraction method addressed in related literature. This paper uses tolerance calibration to generate different reduction rates and utilizes four types of evaluation factors, statistical measurement, image measurement, spectral analysis and elevation cumulative probability distribution graph, to make a quantitative analysis of profile variation. The study result indicates that real profile elevation data, manipulated with varied reduction approaches, then reconstructed by means of fractal interpolation can produce data points with a higher resolution than those originally observed, thereby the reconstructed profile gets more natural and real details.     

Limit theorems of the theory of discrete periodic splines

It is known from the discrete harmonic analysis that the interpolation problem with equidistant interpolation points has a unique solution. If the right-hand sides in the interpolation problem are fixed, the spline depends on two parameters: the spline order and the number of points located between neighboring interpolation points. We find explicit expressions for the limits of interpolation spllines with respect to each parameter separately and show that both repeated limits exist. We also prove that these repeated limits are equal and their value is an interpolation trigonometric polynomial. Bibliography: 10 titles. Illustrations: 2 figures.     

On general Hermite trigonometric interpolation

Summary A sequence of general Hermite trigonometric interpolation polynomials with equidistant interpolation points is given. Integrating these interpolation formulae a sequence of quadrature formulae for the integration of periodic functions is obtained. Derivative-free remainders are stated for these interpolation and quadrature formulae.This work was done at the Max-Planck-Institut für Physik und Astrophysik, München.     

Choice of specific interpolation points in quadrature formula for fast oscillating integrals

Integrals containing trigonometric functions with non-linear argument having a large multiplier are treated with the help of a specific choice of interpolation points. This allows excluding the contribution of internal points of interpolation and determining upper and low bounds of the considered integral.     

Newton interpolation at Leja points

The Newton form is a convenient representation for interpolation polynomials. Its sensitivity to perturbations depends on the distribution and ordering of the interpolation points. The present paper bounds the growth of the condition number of the Newton form when the interpolation points are Leja points for fairly general compact sets K in the complex plane. Because the Leja points are defined recursively, they are attractive to use with the Newton form. If K is an interval, then the Leja points are distributed roughly like Chebyshev points. Our investigation of the Newton form defined by interpolation at Leja points suggests an ordering scheme for arbitrary interpolation points.Research supported in part by NSF under Grant DMS-8704196 and by U.S. Air Force Grant AFSOR-87-0102.On leave from University of Kentucky, Department of Mathematics, Lexington, KY 40506, U.S.A.     

On polynomial interpolation of two variables

Polynomial interpolation of two variables based on points that are located on multiple circles is studied. First, the poisedness of a Birkhoff interpolation on points that are located on several concentric circles is established. Second, using a factorization method, the poisedness of a Hermite interpolation based on points located on various circles, not necessarily concentric, is established. Even in the case of Lagrange interpolation, this gives many new sets of poised interpolation points.     

Rationale Interpolation,Normalität und Monosplines

This paper is concerned with interpolation by rational functions. Conditions are given which ensure that the solution has no poles between the points of interpolation. Comparison theorems on the error of the interpolation are also derived.     

A new formal approach to the rational interpolation problem

Summary An elegant and fast recursive algorithm is developed to solve the rational interpolation problem in a complementary way compared to existing methods. We allow confluent interpolation points, poles, and infinity as one of the interpolation points. Not only one specific solution is given but a nice parametrization of all solutions. We also give a linear algebra interpretation of the problem showing that our algorithm can also be used to handle a specific class of structured matrices.     

Rational interpolation through the optimal attachment of poles to the interpolating polynomial

After recalling some pitfalls of polynomial interpolation (in particular, slopes limited by Markov's inequality) and rational interpolation (e.g., unattainable points, poles in the interpolation interval, erratic behavior of the error for small numbers of nodes), we suggest an alternative for the case when the function to be interpolated is known everywhere, not just at the nodes. The method consists in replacing the interpolating polynomial with a rational interpolant whose poles are all prescribed, written in its barycentric form as in [4], and optimizing the placement of the poles in such a way as to minimize a chosen norm of the error. This revised version was published online in June 2006 with corrections to the Cover Date.     

Interpolation of analytic functions with finitely many singularities

We consider an interpolation process for the class of functions with finitely many singular points by means of rational functions whose poles coincide with the singular points of the function under interpolation. The interpolation nodes form a triangular matrix. We find necessary and sufficient conditions for the uniform convergence of sequences of interpolation fractions to the function under interpolation on every compact set disjoint from the singular points of the function and other conditions for convergence. We generalize and improve the familiar results on the interpolation of functions with finitely many singular points by rational fractions and of entire functions by polynomials.     

有理插值中不可达点的研究

通过对一元Thiele型连分式插值和二元Newton-Thiele型混合有理插值中不可达点的分析,给出了一种判断不可达点的方法.而且,对于任意给定的插值条件,通过构造带参数的Thiele型切触插值和二元Newton-Thiele型混合切触有理插值,使得不可达点变成可达点.数值例子也说明了这种方法的有效性.     

Bivariate spline interpolation at grid points

Summary. We describe algorithms for constructing point sets at which interpolation by spaces of bivariate splines of arbitrary degree and smoothness is possible. The splines are defined on rectangular partitions adding one or two diagonals to each rectangle. The interpolation sets are selected in such a way that the grid points of the partition are contained in these sets, and no large linear systems have to be solved. Our method is to generate a net of line segments and to choose point sets in these segments which satisfy the Schoenberg-Whitney condition for certain univariate spline spaces such that a principle of degree reduction can be applied. In order to include the grid points in the interpolation sets, we give a sufficient Schoenberg-Whitney type condition for interpolation by bivariate splines supported in certain cones. This approach is completely different from the known interpolation methods for bivariate splines of degree at most three. Our method is illustrated by some numerical examples. Received October 5, 1992 / Revised version received May 13, 1994