Construction of biorthogonal wavelets from pseudo-splines

Pseudo-splines constitute a new class of refinable functions with B-splines, interpolatory refinable functions and refinable functions with orthonormal shifts as special examples. Pseudo-splines were first introduced by Daubechies, Han, Ron and Shen in [Framelets: MRA-based constructions of wavelet frames, Appl. Comput. Harmon. Anal. 14(1) (2003), 1–46] and Selenick in [Smooth wavelet tight frames with zero moments, Appl. Comput. Harmon. Anal. 10(2) (2001) 163–181], and their properties were extensively studied by Dong and Shen in [Pseudo-splines, wavelets and framelets, 2004, preprint]. It was further shown by Dong and Shen in [Linear independence of pseudo-splines, Proc. Amer. Math. Soc., to appear] that the shifts of an arbitrarily given pseudo-spline are linearly independent. This implies the existence of biorthogonal dual refinable functions (of pseudo-splines) with an arbitrarily prescribed regularity. However, except for B-splines, there is no explicit construction of biorthogonal dual refinable functions with any given regularity. This paper focuses on an implementable scheme to derive a dual refinable function with a prescribed regularity. This automatically gives a construction of smooth biorthogonal Riesz wavelets with one of them being a pseudo-spline. As an example, an explicit formula of biorthogonal dual refinable functions of the interpolatory refinable function is given.     

Multivariate refinable Hermite interpolant

We introduce a general definition of refinable Hermite interpolants and investigate their general properties. We also study a notion of symmetry of these refinable interpolants. Results and ideas from the extensive theory of general refinement equations are applied to obtain results on refinable Hermite interpolants. The theory developed here is constructive and yields an easy-to-use construction method for multivariate refinable Hermite interpolants. Using this method, several new refinable Hermite interpolants with respect to different dilation matrices and symmetry groups are constructed and analyzed.

Some of the Hermite interpolants constructed here are related to well-known spline interpolation schemes developed in the computer-aided geometric design community (e.g., the Powell-Sabin scheme). We make some of these connections precise. A spline connection allows us to determine critical Hölder regularity in a trivial way (as opposed to the case of general refinable functions, whose critical Hölder regularity exponents are often difficult to compute).

While it is often mentioned in published articles that ``refinable functions are important for subdivision surfaces in CAGD applications", it is rather unclear whether an arbitrary refinable function vector can be meaningfully applied to build free-form subdivision surfaces. The bivariate symmetric refinable Hermite interpolants constructed in this article, along with algorithmic developments elsewhere, give an application of vector refinability to subdivision surfaces. We briefly discuss several potential advantages offered by such Hermite subdivision surfaces.

     

Regularity of multiwavelets

The motivation for this paper is an interesting observation made by Plonka concerning the factorization of the matrix symbol associated with the refinement equation for B-splines with equally spaced multiple knots at integers and subsequent developments which relate this factorization to regularity of refinable vector fields over the real line. Our intention is to contribute to this train of ideas which is partially driven by the importance of refinable vector fields in the construction of multiwavelets. The use of subdivision methods will allow us to consider the problem almost entirely in the spatial domain and leads to exact characterizations of differentiability and Hölder regularity in arbitrary L _p spaces. We first study the close relationship between vector subdivision schemes and a generalized notion of scalar subdivision schemes based on bi-infinite matrices with certain periodicity properties. For the latter type of subdivision scheme we will derive criteria for convergence and Hölder regularity of the limit function, which mainly depend on the spectral radius of a bi-infinite matrix induced by the subdivision operator, and we will show that differentiability of the limit functions can be characterized by factorization properties of the subdivision operator. By switching back to vector subdivision we will transfer these results to refinable vectors fields and obtain characterizations of regularity by factorization and spectral radius properties of the symbol associated to the refinable vector field. Finally, we point out how multiwavelets can be generated from orthonormal refinable bi-infinite vector fields.     

Two-scale symbol and autocorrelation symbol for B-splines with multiple knots

The Fourier transforms of B-splines with multiple integer knots are shown to satisfy a simple recursion relation. This recursion formula is applied to derive a generalized two-scale relation for B-splines with multiple knots. Furthermore, the structure of the corresponding autocorrelation symbol is investigated. In particular, it can be observed that the solvability of the cardinal Hermite spline interpolation problem for spline functions of degree 2m+1 and defectr, first considered by Lipow and Schoenberg [9], is equivalent to the Riesz basis property of our B-splines with degreem and defectr. In this way we obtain a new, simple proof for the assertion that the cardinal Hermite spline interpolation problem in [9] has a unique solution.     

Two-scale symbol and autocorrelation symbol for B-splines with multiple knots

The Fourier transforms of B-splines with multiple integer knots are shown to satisfy a simple recursion relation. This recursion formula is applied to derive a generalized two-scale relation for B-splines with multiple knots. Furthermore, the structure of the corresponding autocorrelation symbol is investigated. In particular, it can be observed that the solvability of the cardinal Hermite spline interpolation problem for spline functions of degree 2m+1 and defectr, first considered by Lipow and Schoenberg [9], is equivalent to the Riesz basis property of our B-splines with degreem and defectr. In this way we obtain a new, simple proof for the assertion that the cardinal Hermite spline interpolation problem in [9] has a unique solution.     

Multivariate Compactly Supported Fundamental Refinable Functions, Duals, and Biorthogonal Wavelets

In areas of geometric modeling and wavelets, one often needs to construct a compactly supported refinable function φ which has sufficient regularity and which is fundamental for interpolation [that means, φ(0)=1 and φ(α)=0 for all α∈ Z ^s ∖{0}].
Low regularity examples of such functions have been obtained numerically by several authors, and a more general numerical scheme was given in [1]. This article presents several schemes to construct compactly supported fundamental refinable functions, which have higher regularity, directly from a given, continuous, compactly supported, refinable fundamental function φ. Asymptotic regularity analyses of the functions generated by the constructions are given.The constructions provide the basis for multivariate interpolatory subdivision algorithms that generate highly smooth surfaces.
A very important consequence of the constructions is a natural formation of pairs of dual refinable functions, a necessary element in constructing biorthogonal wavelets. Combined with the biorthogonal wavelet construction algorithm for a pair of dual refinable functions given in [2], we are able to obtain symmetrical compactly supported multivariate biorthogonal wavelets which have arbitrarily high regularity. Several examples are computed.     

渐近于Hermite多项式的双正交系统

本文利用渐近于Gauss函数的函数类?,给出渐近于Hermite正交多项式的一类Appell多项式的构造方法,使得该序列与?的n阶导数之间构成了一组双正交系统.利用此结果,本文得到多种正交多项式和组合多项式的渐近性质.特别地,由N阶B样条所生成的Appell多项式序列恰为N阶Bernoulli多项式.从而,Bernoulli多项式与B样条的导函数之间构成了一组双正交系统,且标准化之后的Bernoulli多项式的渐近形式为Hermite多项式.由二项分布所生成的Appell序列为Euler多项式,从而,Euler多项式与二项分布的导函数之间构成一组双正交系统,且标准化之后的Euler多项式渐近于Hermite多项式.本文给出Appell序列的生成函数满足的尺度方程的充要条件,给出渐近于Hermite多项式的函数列的判定定理.应用该定理,验证广义Buchholz多项式、广义Laguerre多项式和广义Ultraspherical(Gegenbauer)多项式渐近于Hermite多项式的性质,从而验证超几何多项式的Askey格式的成立.     

Compactly supported (bi)orthogonal wavelets generated by interpolatory refinable functions

This paper provides several constructions of compactly supported wavelets generated by interpolatory refinable functions. It was shown in [7] that there is no real compactly supported orthonormal symmetric dyadic refinable function, except the trivial case; and also shown in [10,18] that there is no compactly supported interpolatory orthonormal dyadic refinable function. Hence, for the dyadic dilation case, compactly supported wavelets generated by interpolatory refinable functions have to be biorthogonal wavelets. The key step to construct the biorthogonal wavelets is to construct a compactly supported dual function for a given interpolatory refinable function. We provide two explicit iterative constructions of such dual functions with desired regularity. When the dilation factors are larger than 3, we provide several examples of compactly supported interpolatory orthonormal symmetric refinable functions from a general method. This leads to several examples of orthogonal symmetric (anti‐symmetric) wavelets generated by interpolatory refinable functions. This revised version was published online in June 2006 with corrections to the Cover Date.     

Band-limited refinable functions for wavelets and framelets

Extending band-limited constructions of orthonormal refinable functions, a special class of periodic functions is used to generate a family of band-limited refinable functions. Characterizations of Riesz bases and frames formed by integer shifts of these refinable functions are obtained. Such families of refinable functions are employed to construct band-limited biorthogonal wavelet bases and biframes with desirable time-frequency localization.     

Refinable bivariate quartic -splines for multi-level data representation and surface display

In this paper, a second-order Hermite basis of the space of -quartic splines on the six-directional mesh is constructed and the refinable mask of the basis functions is derived. In addition, the extra parameters of this basis are modified to extend the Hermite interpolating property at the integer lattices by including Lagrange interpolation at the half integers as well. We also formulate a compactly supported super function in terms of the basis functions to facilitate the construction of quasi-interpolants to achieve the highest (i.e., fifth) order of approximation in an efficient way. Due to the small (minimum) support of the basis functions, the refinable mask immediately yields (up to) four-point matrix-valued coefficient stencils of a vector subdivision scheme for efficient display of -quartic spline surfaces. Finally, this vector subdivision approach is further modified to reduce the size of the coefficient stencils to two-point templates while maintaining the second-order Hermite interpolating property.

     

Affine Frames of Multivariate Box Splines and their Affine Duals

We give a simple and explicit construction of primal and dual wavelet filters based on refinable multivariate splines (with respect to dilation matrices M) such that the corresponding wavelet functions generate dual affine frames of arbitrarily high regularity. Furthermore, the number of wavelets does not depend on the regularity. We apply the method also to generalized B-splines.     

Generalized Pseudo-Butterworth Refinable Functions

In this paper we propose the generalized pseudo-Butterworth refinable functions which involve pseudo-splines of type I and II, Butterworth refinable functions, pseudo-Butterworth refinable functions, and almost all symmetric and causal fractional B-splines. Furthermore, the convergence of cascade algorithms associated with the new masks is proved, and Riesz wavelet bases in L ₂(?) corresponding to the parameters are constructed. The regularity of the generalized pseudo-Butterworth refinable functions is also analyzed by Fourier analysis.     

Multiwavelet approximation methods for pseudodifferential equations on curves. Stability and convergence analysis

We develop a stability and convergence analysis of Galerkin–Petrov schemes based on a general setting of multiresolution generated by several refinable functions for the numerical solution of pseudodifferential equations on smooth closed curves. Particular realizations of such a multiresolution analysis are trial spaces generated by biorthogonal wavelets or by splines with multiple knots. The main result presents necessary and sufficient conditions for the stability of the numerical method in terms of the principal symbol of the pseudodifferential operator and the Fourier transforms of the generating multiscaling functions as well as of the test functionals. Moreover, optimal convergence rates for the approximate solutions in a range of Sobolev spaces are established. This revised version was published online in June 2006 with corrections to the Cover Date.     

Stability Implies Convergence of Cascade Algorithms in Sobolev Space

This article proves that the stability of the shifts of a refinable function vector ensures the convergence of the corresponding cascade algorithm in Sobolev space to which the refinable function vector belongs. An example of Hermite interpolants is presented to illustrate the result.     

Biorthogonal Wavelet Expansions

This paper is concerned with developing conditions on a given finite collection of compactly supported algebraically linearly independent refinable functions that insure the existence of biorthogonal systems of refinable functions with similar properties. In particular, we address the close connection of this issue with stationary subdivision schemes. Date received: May 20, 1995. Date revised: March 2, 1996.     

Generalized cardinal B-splines: Stability, linear independence, and appropriate scaling matrices

Generalized cardinal B-splines are defined as convolution products of characteristic functions of self-affine lattice tiles with respect to a given integer scaling matrix. By construction, these generalized splines are refinable functions with respect to the scaling matrix and therefore they can be used to define a multiresolution analysis and to construct a wavelet basis. In this paper, we study the stability and linear independence properties of the integer translates of these generalized spline functions. Moreover, we give a characterization of the scaling matrices to which the construction of the generalized spline functions can be applied.     

A Class of Orthogonal Refinable Functions and Wavelets

We give a construction, for any n 2, of a space S of spline functions of degree n – 1 with simple knots in (1/4)Z which is generated by a triple of refinable, orthogonal functions with compact support. Indeed, the result holds more generally by replacing the B-spline of degree n – 1 with simple knots at the integers by any continuous refinable function whose mask is a Hurwitz polynomial of degree n. A simple construction is also given for the corresponding wavelets.     

Biorthogonal wavelets on Vilenkin groups

We describe algorithms for constructing biorthogonal wavelet systems and refinable functions whose masks are generalized Walsh polynomials. We give new examples of biorthogonal compactly supported wavelets on Vilenkin groups.     

Curve and surface construction using Hermite subdivision schemes

In this paper we present a very efficient Hermite subdivision scheme, based on rational functions, and outline its potential applications, with special emphasis on the construction of cubic-like B-splines — well suited for the design of constrained curves and surfaces.     

From Hermite to stationary subdivision schemes in one and several variables

Vector and Hermite subdivision schemes both act on vector data, but since the latter one interprets the vectors as function values and consecutive derivatives they differ by the “renormalization” of the Hermite scheme in any step. In this paper we give an algebraic factorization method in one and several variables to relate any Hermite subdivision scheme that satisfies the so–called spectral condition to a vector subdivision scheme. These factorizations are natural extensions of the “zero at π” condition known for the masks of refinable functions. Moreover, we show how this factorization can be used to investigate different forms of convergence of the Hermite scheme and why the multivariate situation is conceptionally more intricate than the univariate one. Finally, we give some examples of such factorizations.