具有重节点的分段Pade''逼近的一个算法

Baker在[1]中提出了具有重节点的Pade’逼近问题,但提供的算法很繁.我们发现,具有重节点的Pade’逼近和有理切触插值有关.基于这种想法,我们先给出分段Pade’逼近的概念,然后给出一个一般算法.     

几种有理插值函数的逼近性质

1　引　　言在曲线和曲面设计中,样条插值是有用的和强有力的工具.不少作者已经研究了很多种类型的样条插值[1,2,3,4].近些年来,有理插值样条,特别是三次有理插值样条,以及它们在外型控制中的应用,已有了不少工作[5,6,7].有理插值样条的表达式中有某些参数,正是由于这些参数,有理插值样条在外型控制中充分显示了它的灵活性;但也正是由于这些参数,使它的逼近性质的研究增加了困难.因此,关于有理插值样条的逼近性质的研究很少见诸文献.本文在第二节首先叙述几种典型的有理插值样条,其中包括分母为一次、二次的三次有理插值样条和仅基于函数值…     

Thiele型向量连分式的收敛性定理

Thiele型向量连分式,不仅可用来解决一元和多元向量有理插值问题[1-3],一元和多元向量切触有理插值问题[3],还可用来研究向量Pade逼近及向量连分式逼近[1,3]。本文给出了这种连分式的收敛性定理,并把著名的Pringsheim定理推广到向量连分式上去。     

加权有理三次插值的逼近性质及其应用

利用带导数和不带导数的分母为线性的有理三次插值样条构造了一类加权有理三次插值函数,利用这种插值方法,将样条曲线严格约束于给定的折线之上、之下或之间的问题都可以得到解决同时还研究了这种加权有理三次插值的逼近性质。     

Floater和Hormann插值算子的勒贝格常数

1引言设f(x)是定义在区间[a,b]上的连续函数,插值节点x_(i)(i=0,1,…,n)满足a=x_(0)相似文献   

关于有理插值函数存在性的确定

在本中，我们利用Newton插值多项式，改进了[1]中的方法，使其能更简便，快速，严谨地判别有理插值函数的存在性，并在其存在时给出相应的插值有理函数的具体表达式。     

二元向量值有理插值存在性的一种判别方法

1 引言 一元向量值有理插值问题在[1-5]中有了比较系统的研究．文[6—13]成功地将一无的结果推广到了二元的情形,但它们采用的大多是向量值连分式的方法,且没有给出二元向量值有理插值存在性的判别方法及其证明．本文利用二元Newton插值公式,     

某些修正有理插值在L_w~p空间的逼近(英文)

本文拓广了T．Herman和P．Vertes研究的有理插值，引入某些修正的有理插值，并给出它们在Lpw空间的逼近阶，其中W（x）＝（1-x2）1/2．     

一类二元有理插值

应用文[1]新近建立的Gould-Hsu反演的双变量形式,本文研究了一类二元有理插值公式的构造,确定了该类插值函数所表现的函数类,并给出了差分表计算的递归公式     

一种新的带参数双三次有理插值样条的有界性与点控制

文[19]中,作者构造了一种基于函数值的带参数的分子为双三次、分母为双二次的二元有理插值样条.本文进一步研究该种二元有理插值样条的有界性,给出插值的逼近表达式,讨论插值曲面形状的点控制问题.在插值条件不变的情况下,插值区域内任一点插值函数的值可以根据设计的需要通过对参数的选取修改,从而达到插值曲面局部修改的目的.     

MATRIX ALGORITHMS AND ERROR FORMULA FOR BIVARIATE THIELE-TYPE RECTANGULAR MATRIX VALUED RATIONAL INTERPOLATION

A new method for the construction of bivariate matrix valued rational interpolants (BGIRI) on a rectangular grid is presented in [6]. The rational interpolants are of Thiele-type continued fraction form with scalar denominator. The generalized inverse introduced by [3]is gen-eralized to rectangular matrix case in this paper. An exact error formula for interpolation is ob-tained, which is an extension in matrix form of bivariate scalar and vector valued rational interpola-tion discussed by Siemaszko[l2] and by Gu Chuangqing [7] respectively. By defining row and col-umn-transformation in the sense of the partial inverted differences for matrices, two type matrix algorithms are established to construct corresponding two different BGIRI, which hold for the vec-tor case and the scalar case.     

样条型矩阵有理插值

The matrix valued rational interpolation is very useful in the partial realization problem and model reduction for all the linear system theory. Lagrange basic functions have been used in matrix valued rational interpolation. In this paper, according to the property of cardinal spline interpolation, we constructed a kind of spline type matrix valued rational interpolation, which based on cardinal spline. This spline type interpolation can avoid instability of high order polynomial interpolation and we obtained a useful formula.     

Vector-valued,rational interpolants III

In 1963, Wynn proposed a method for rational interpolation of vector-valued quantities given on a set of distinct interpolation points. He used continued fractions, and generalized inverses for the reciprocals of vector-valued quantities. In this paper, we present an axiomatic approach to vector-valued rational interpolation. Uniquely defined interpolants are constructed for vector-valued data so that the components of the resulting vector-valued rational interpolant share a common denominator polynomial. An explicit determinantal formula is given for the denominator polynomial for the cases of (i) vector-valued rational interpolation on distinct real or complex points and (ii) vector-valued Padé approximation. We derive the connection with theε-algorithm of Wynn and Claessens, and we establish a five-term recurrence relation for the denominator polynomials.     

Weighted interpolation for equidistant nodes

Weighted Lagrange interpolation is proposed for solving Lagrange interpolation problems on equidistant or almost equidistant data. Good condition numbers are found in the case of rational interpolants whose denominator has degree about twice the number of data to be interpolated. Since the degree of the denominator is higher than that of the numerator, simple functions like constants and linear polynomials will not be reproduced. Furthermore, the interpolant cannot be expressed by a barycentric formula. As a counterpart, the interpolation algorithm is simple and leads to small Lebesgue constants.     

Computing near-best fixed pole rational interpolants

We study rational interpolation formulas on the interval [−1,1] for a given set of real or complex conjugate poles outside this interval. Interpolation points which are near-best in a Chebyshev sense were derived in earlier work. The present paper discusses several computation aspects of the interpolation points and the corresponding interpolants. We also study a related set of points (that includes the end points), which is more suitable for applications in rational spectral methods. Some examples are given at the end of this paper.     

The compass (star) identity for vector-valued rational interpolants

The compass identity (Wynn's five point star identity) for Padé approximants connects neighbouring elements called N, S, E, W and C in the Padé table. Its form has been extended to the cases of rational interpolation of ordinary (scalar) data and interpolation of vector-valued data. In this paper, full specifications of the associated five point identity for the scalar denominator polynomials and the vector numerator polynomials of the vector-valued rational interpolants on real data points are given, as well as the related generalisations of Frobenius' identities. Unique minimal forms of the polynomials constituting the interpolants and results about unattainable points correspond closely to their counterparts in the scalar case. This revised version was published online in June 2006 with corrections to the Cover Date.     

A Neville-like method via continued fractions

As we know, the classical Neville's algorithm is an effective method used to solve the interpolation problem by polynomials. In this paper, we adopt the idea of the Neville's algorithm to construct a kind of blending rational interpolants via continued fractions. For a given set of support points, there are many ways to build up the interpolation schemes, by which we mean that there are many choices to make to determine the initial interpolants on subsets of support points and then update them step by step to form a solution to the full interpolation problem. Numerical examples are given to show the advantage of our method and a multivariate analogy is also discussed.     

连续可导函数类的最优拉格朗日插值

在构造拉格朗日插值算法时,插值结点的选择是十分重要的.给定一个足够光滑的函数,如果结点选择的不好,当插值结点个数趋于无穷时,插值函数不收敛于函数本身.例如龙格现象:对于龙格函数f(x)=1/1+25x^2,如果拉格朗日插值的结点取[-1,1]上的等距结点,那么逼近的误差会随着结点个数增多而趋于无穷大⑴,由此可知插值结点的选择尤为重要.     

Padé–type rational and barycentric interpolation

In this paper, we consider the particular case of the general rational Hermite interpolation problem where only the value of the function is interpolated at some points, and where the function and its first derivatives agree at the origin. Thus, the interpolants constructed in this way possess a Padé–type property at 0. Numerical examples show the interest of the procedure. The interpolation procedure can be easily modified to introduce a partial knowledge on the poles and the zeros of the function to approximated. A strategy for removing the spurious poles is explained. A formula for the error is proved in the real case. Applications are given.     

Lower complexity bounds for interpolation algorithms

We introduce and discuss a new computational model for the Hermite-Lagrange interpolation with nonlinear classes of polynomial interpolants. We distinguish between an interpolation problem and an algorithm that solves it. Our model includes also coalescence phenomena and captures a large variety of known Hermite-Lagrange interpolation problems and algorithms. Like in traditional Hermite-Lagrange interpolation, our model is based on the execution of arithmetic operations (including divisions) in the field where the data (nodes and values) are interpreted and arithmetic operations are counted at unit cost. This leads us to a new view of rational functions and maps defined on arbitrary constructible subsets of complex affine spaces. For this purpose we have to develop new tools in algebraic geometry which themselves are mainly based on Zariski’s Main Theorem and the theory of places (or equivalently: valuations). We finish this paper by exhibiting two examples of Lagrange interpolation problems with nonlinear classes of interpolants, which do not admit efficient interpolation algorithms (one of these interpolation problems requires even an exponential quantity of arithmetic operations in terms of the number of the given nodes in order to represent some of the interpolants).In other words, classic Lagrange interpolation algorithms are asymptotically optimal for the solution of these selected interpolation problems and nothing is gained by allowing interpolation algorithms and classes of interpolants to be nonlinear. We show also that classic Lagrange interpolation algorithms are almost optimal for generic nodes and values. This generic data cannot be substantially compressed by using nonlinear techniques.We finish this paper highlighting the close connection of our complexity results in Hermite-Lagrange interpolation with a modern trend in software engineering: architecture tradeoff analysis methods (ATAM).