Changeable degree spline basis functions

A B-spline basis function is a piecewise function of polynomials of equal degree on its support interval. This paper extends B-spline basis functions to changeable degree spline (CD-spline for short) basis functions, each of which may consist of polynomials of different degrees on its support interval. The CD-spline basis functions possess many B-spline-like properties and include the B-spline basis functions as subcases. Their corresponding parametric curves, called CD-spline curves, are like B-spline curves and also have many good properties. If we use the CD-spline basis functions to design a curve made up of polynomial segments of different degrees, the number of control points may be decreased.     

How to build all Chebyshevian spline spaces good for geometric design?

In the present work we determine all Chebyshevian spline spaces good for geometric design. By Chebyshevian spline space we mean a space of splines with sections in different Extended Chebyshev spaces and with connection matrices at the knots. We say that such a spline space is good for design when it possesses blossoms. To justify the terminology, let us recall that, in this general framework, existence of blossoms (defined on a restricted set of tuples) makes it possible to develop all the classical geometric design algorithms for splines. Furthermore, existence of blossoms is equivalent to existence of a B-spline bases both in the spline space itself and in all other spline spaces derived from it by insertion of knots. We show that Chebyshevian spline spaces good for design can be described by linear piecewise differential operators associated with systems of piecewise weight functions, with respect to which the connection matrices are identity matrices. Many interesting consequences can be drawn from the latter characterisation: as an example, all Chebsyhevian spline spaces good for design can be built by means of integral recurrence relations.     

Multiresolution exponential B-splines and singularly perturbed boundary problem

The paper considers how cardinal exponential B-splines can be applied in solving singularly perturbed boundary problems. The exponential nature and the multiresolution property of these splines are essential for an accurate simulation of a singular behavior of some differential equation solutions. Based on the knowledge that the most of exponential B-spline properties coincide with those of polynomial splines (smoothness, compact support, positivity, partition of unity, reconstruction of polynomials, recursion for derivatives), one novel algorithm is proposed. It merges two well known approaches for solving such problems, fitted operator and fitted mesh methods. The exponential B-spline basis is adapted for an interval because a considered problem is solved on a bounded domain.      

A piecewise approach to piecewise approximation

Tensor-product B-spline surfaces offer a convenient means for representing a set of bivariate data, especially if many surface evaluations are required. This is because the compact support property of the tensor-product spline allows the spline value to be obtained in a time that is (almost) independent of the number of coefficients used to define the surface. The main calculation is the precomputation involved in fitting the data and this can be impractically large if there are many spline coefficients to be calculated. Since the surface produced may be evaluated locally and efficiently, it would be advantageous to exploit local properties in order to fit the data in a piecewise manner. An algorithm to do this is presented.     

A MIQCP formulation for B-spline constraints

This paper presents a mixed-integer quadratically constrained programming (MIQCP) formulation for B-spline constraints. The formulation can be used to obtain an exact MIQCP reformulation of any spline-constrained optimization problem problem, provided that the polynomial spline functions are continuous. This reformulation allows practitioners to use a general-purpose MIQCP solver, instead of a special-purpose spline solver, when solving B-spline constrained problems. B-splines are a powerful and widely used modeling tool, previously restricted from optimization due to lack of solver support. This contribution may encourage practitioners to use B-splines to model constraint functions. However, as the numerical study suggests, there is still a large gap between the solve times of the general-purpose solvers using the proposed formulation, and the special-purpose spline solver CENSO, the latter being significantly lower.     

Relations between the support of a compactly supported function and the exponential-polynomials spanned by its integer translates

The interrelation between the shape of the support of a compactly supported function and the space of all exponential-polynomials spanned by its integer translates is examined. The results obtained are in terms of the behavior of these exponential-polynomials on certain finite subsets ofZ ^s , which are determined by the support of the given function. Several applications are discussed. Among these is the construction of quasi-interpolants of minimal support and the construction of a piecewise-polynomial whose integer translates span a polynomial space which is not scale-invariant. As to polynomial box splines, it is proved here that in many cases a polynomial box spline admits a certain optimality condition concerning the space of the total degree polynomials spanned by its integer translates: This space is maximal compared with the spaces corresponding to other functions with the same supportCommunicated by Klaus Höllig.     

Some extremal properties of multivariate polynomial splines in the metric Lp (Rd )

We constructed a kind of continuous multivariate spline operators as the approximation tools of the multivariate functions on the Bd instead of the usual multivariate cardinal interpolation oper-ators of splines, and obtained the approximation error by this kind of spline operators. Meantime, by the results, we also obtained that the spaces of multivariate polynomial splines are weakly asymptoti-cally optimal for the Kolmogorov widths and the linear widths of some anisotropic Sobolev classes of smooth functions on Bd in the metric Lp(Bd).     

有限区间内四阶样条小波的构造

用有限区间上的截断4阶B样条,构造了有限区间上的4阶样条小波。这些小波由边界小波和内部小波组成,对某一尺度,它们组成了有限维的小波空间。于是,任何有限区间上的函数皆可表示为该区间上的尺度函数和小波函数的有限和,即小波级数,这克服了用无穷区间上的小波进行有限信号处理时,在边界上误差较大的不足,同时将该小波用于偏微分方程具有同样重要的意义。     

Data compression based on the cubic B-spline wavelet with uniform two-scale relation

The aim of this paper is to investigate the potential artificial compression which can be achieved using an interval multiresolution analysis based on a semiorthogonal cubic B-spline wavelet. The Chui-Quak [1] spline multiresolution analysis for the finite interval has been modified [2] so as to be characterized by natural spline projection and uniform two-scale relation. Strengths and weaknesses of the semiorthogonal wavelet as regards artificial compression and data smoothing by the method of thresholding wavelet coefficients are indicated.     

Computing ECT-B-splines recursively

ECT-spline curves for sequences of multiple knots are generated from different local ECT-systems via connection matrices. Under appropriate assumptions there is a basis of the space of ECT-splines consisting of functions having minimal compact supports, normalized to form a nonnegative partition of unity. The basic functions can be defined by generalized divided differences [24]. This definition reduces to the classical one in case of a Schoenberg space. Under suitable assumptions it leads to a recursive method for computing the ECT-B-splines that reduces to the de Boor–Mansion–Cox recursion in case of ordinary polynomial splines and to Lyche's recursion in case of Tchebycheff splines. For sequences of simple knots and connection matrices that are nonsingular, lower triangular and totally positive the spline weights are identified as Neville–Aitken weights of certain generalized interpolation problems. For multiple knots they are limits of Neville–Aitken weights. In many cases the spline weights can be computed easily by recurrence. Our approach covers the case of Bézier-ECT-splines as well. They are defined by different local ECT-systems on knot intervals of a finite partition of a compact interval [a,b] connected at inner knots all of multiplicities zero by full connection matrices A ^[i] that are nonsingular, lower triangular and totally positive. In case of ordinary polynomials of order n they reduce to the classical Bézier polynomials. We also present a recursive algorithm of de Boor type computing ECT-spline curves pointwise. Examples of polynomial and rational B-splines constructed from given knot sequences and given connection matrices are added. For some of them we give explicit formulas of the spline weights, for others we display the B-splines or the B-spline curves. *Supported in part by INTAS 03-51-6637.     

Some extremal properties of multivariate polynomial splines in the metric Lp(ℝd)

We constructed a kind of continuous multivariate spline operators as the approximation tools of the multivariate functions on the (?^d instead of the usual multivariate cardinal interpolation operators of splines, and obtained the approximation error by this kind of spline operators. Meantime, by the results, we also obtained that the spaces of multivariate polynomial splines are weakly asyrnptotically optimal for the Kolrnogorov widths and the linear widths of some anisotropic Sobolev classes of smooth functions on (?^d in the metric L_p((?^d).     

B-splines with homogenized knots

The use of homogenized knots for manipulating univariate polynomials by blossoming algorithms is extended to piecewise polynomials. A generalization of the B-spline to homogenized knots is studied. The new B-spline retains the triangular blossoming algorithms for evaluation, differentiation and knot insertion. Moreover, the B-spline is locally supported and a Marsden’s identity exists. Spaces of natural splines and certain polynomial spline spaces with more general continuity properties than ordinary splines have bases of B-splines over homogenized knots. Applications to nonpolynomial splines such as trigonometric and hyperbolic splines are made.     

B-splines with homogenized knots

The use of homogenized knots for manipulating univariate polynomials by blossoming algorithms is extended to piecewise polynomials. A generalization of the B-spline to homogenized knots is studied. The new B-spline retains the triangular blossoming algorithms for evaluation, differentiation and knot insertion. Moreover, the B-spline is locally supported and a Marsden’s identity exists. Spaces of natural splines and certain polynomial spline spaces with more general continuity properties than ordinary splines have bases of B-splines over homogenized knots. Applications to nonpolynomial splines such as trigonometric and hyperbolic splines are made.     

A Data-Reduction Strategy for Splines with Applications to the Approximation of Functions and Data

We present a strategy for reducing the number of knots of agiven B-spline function without perturbing the spline by morethan a given tolerance. The number and location of the remainingknots is determined automatically. Knot removal can be usedsuccessfully to fit a spline to functions and data. We illustratethis with several examples.     

On the Equivalence Between Existence of B-Spline Bases and Existence of Blossoms

In spline spaces with sections in arbitrary extended Chebyshev spaces and with connections defined by arbitrary lower triangular matrices with positive diagonal elements, we prove that existence of B-spline bases is equivalent to existence of blossoms. As is now classical, we construct blossoms with the help of osculating flats. As for B-spline bases, this expression denotes normalized basis consisting of minimally supported functions which are positive on the interior of their supports and which satisfy an additional end point condition.     

样条型矩阵有理插值

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.     

Exponential box splines

In a recent paper by Nira Dyn and the author, univariate cardinal exponential B-splines are shown to have a representation similar to the wellknown box spline representation of the univariate cardinal polynomialB-splines. Motivated by this, we construct, for a set ofn directions inZ ^s and a vector of constants λ ?R ⁿ, an “exponential box spline” which has the same smoothness and support as the polynomial box spline, and is a positive piecewise exponential in its support. We derive recurrence relations for the exponential box splines which are simpler than those for the polynomial case. A relatively simple structure of the space spanned by the translates of an exponential box spline is obtained for λ in a certain open dense set ofR ⁿ—the “simple” λ. In this case, the characterization of the local independence of the translates and related topics, as well as the proofs involved, are quite simple when compared with the polynomial case (corresponding toλ = 0).     

Some extremal properties of multivariate polynomial splines in the metric Lp(?d)

We constructed a kind of continuous multivariate spline operators as the approximation tools of the multivariate functions on the (ℝ^d instead of the usual multivariate cardinal interpolation operators of splines, and obtained the approximation error by this kind of spline operators. Meantime, by the results, we also obtained that the spaces of multivariate polynomial splines are weakly asyrnptotically optimal for the Kolrnogorov widths and the linear widths of some anisotropic Sobolev classes of smooth functions on (ℝ^d in the metric L_p((ℝ^d).     

NUAT T-splines of odd bi-degree and local refinement

This paper presents a new kind of spline surfaces, named non-uniform algebraic- trigonometric T-spline surfaces （NUAT T-splines for short） of odd hi-degree. The NUAT T- spline surfaces are defined by applying the T-spline framework to the non-uniform algebraic- trigonometric B-spline surfaces （NUAT B-spline surfaces）. Based on the knot insertion algorithm of the NUAT B-splines, a local refinement algorithm for the NUAT T-splines is given. This algorithm guarantees that the resulting control grid is a T-mesh as the original one. Finally, we prove that, for any NUAT T-spline of odd hi-degree, the linear independence of its blending functions can be determined by computing the rank of the NUAT T-spline-to-NUAT B-spline transformation matrix.     

On hyperbolic summation and hyperbolic moduli of smoothness

This paper deals with the method of hyperbolic summation of tensor product orthogonal spline functions onI ^d. The spaces, defined in terms of the order of the best approximation by the elements of the space spanned by the tensor product functions with indices from a given hyperbolic set, are described both in terms of the coefficients in some basis and as interpolation spaces. Moreover, the hyperbolic modulus of smoothness is studied, and some relations between hyperbolic summation and hyperbolic modulus of smoothness are established.