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.     

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.     

L2(Rs)中插值小波包的分裂技巧和精细分解

对具有任意伸缩矩阵A的插值加细函数，给出对应于L^2(R^s)中的小波包的一个构造方法．采样空间被直接分解来取代对加细函数的符号分解．按照这个方法构造的插值小波包能对基插值空间提供较为精细的分解，因而对自适应的插值给出较好的局部化．     

Analysis and Construction of Multivariate Interpolating Refinable Function Vectors

In this paper, we shall introduce and study a family of multivariate interpolating refinable function vectors with some prescribed interpolation property. Such interpolating refinable function vectors are of interest in approximation theory, sampling theorems, and wavelet analysis. In this paper, we characterize a multivariate interpolating refinable function vector in terms of its mask and analyze the underlying sum rule structure of its generalized interpolatory matrix mask. We also discuss the symmetry property of multivariate interpolating refinable function vectors. Based on these results, we construct a family of univariate generalized interpolatory matrix masks with increasing orders of sum rules and with symmetry for interpolating refinable function vectors. Such a family includes several known important families of univariate refinable function vectors as special cases. Several examples of bivariate interpolating refinable function vectors with symmetry will also be presented.     

A construction of interpolating wavelets on invariant sets

We introduce the concept of a refinable set relative to a family of contractive mappings on a metric space, and demonstrate how such sets are useful to recursively construct interpolants which have a multiscale structure. The notion of a refinable set parallels that of a refinable function, which is the basis of wavelet construction. The interpolation points we recursively generate from a refinable set by a set-theoretic multiresolution are analogous to multiresolution for functions used in wavelet construction. We then use this recursive structure for the points to construct multiscale interpolants. Several concrete examples of refinable sets which can be used for generating interpolatory wavelets are included.

     

Compactly Supported Tight Frames Associated with Refinable Functions

It is well known that in applied and computational mathematics, cardinal B-splines play an important role in geometric modeling (in computer-aided geometric design), statistical data representation (or modeling), solution of differential equations (in numerical analysis), and so forth. More recently, in the development of wavelet analysis, cardinal B-splines also serve as a canonical example of scaling functions that generate multiresolution analyses of L²(−∞,∞). However, although cardinal B-splines have compact support, their corresponding orthonormal wavelets (of Battle and Lemarie) have infinite duration. To preserve such properties as self-duality while requiring compact support, the notion of tight frames is probably the only replacement of that of orthonormal wavelets. In this paper, we study compactly supported tight frames Ψ={ψ¹,…,ψ^N} for L²(−∞,∞) that correspond to some refinable functions with compact support, give a precise existence criterion of Ψ in terms of an inequality condition on the Laurent polynomial symbols of the refinable functions, show that this condition is not always satisfied (implying the nonexistence of tight frames via the matrix extension approach), and give a constructive proof that when Ψ does exist, two functions with compact support are sufficient to constitute Ψ, while three guarantee symmetry/anti-symmetry, when the given refinable function is symmetric.     

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.     

Pseudo-splines,wavelets and framelets

The first type of pseudo-splines were introduced in [I. Daubechies, B. Han, A. Ron, Z. Shen, Framelets: MRA-based constructions of wavelet frames, Appl. Comput. Harmon. Anal. 14 (1) (2003) 1–46; I. Selesnick, Smooth wavelet tight frames with zero moments, Appl. Comput. Harmon. Anal. 10 (2) (2001) 163–181] to construct tight framelets with desired approximation orders via the unitary extension principle of [A. Ron, Z. Shen, Affine systems in $L_2({\mathbb{R}}^d)$: The analysis of the analysis operator, J. Funct. Anal. 148 (2) (1997) 408–447]. In the spirit of the first type of pseudo-splines, we introduce here a new type (the second type) of pseudo-splines to construct symmetric or antisymmetric tight framelets with desired approximation orders. Pseudo-splines provide a rich family of refinable functions. B-splines are one of the special classes of pseudo-splines; orthogonal refinable functions (whose shifts form an orthonormal system given in [I. Daubechies, Orthonormal bases of compactly supported wavelets, Comm. Pure Appl. Math. 41 (1988) 909–996]) are another class of pseudo-splines; and so are the interpolatory refinable functions (which are the Lagrange interpolatory functions at $\mathbb{Z}$ and were first discussed in [S. Dubuc, Interpolation through an iterative scheme, J. Math. Anal. Appl. 114 (1986) 185–204]). The other pseudo-splines with various orders fill in the gaps between the B-splines and orthogonal refinable functions for the first type and between B-splines and interpolatory refinable functions for the second type. This gives a wide range of choices of refinable functions that meets various demands for balancing the approximation power, the length of the support, and the regularity in applications. This paper will give a regularity analysis of pseudo-splines of the both types and provide various constructions of wavelets and framelets. It is easy to see that the regularity of the first type of pseudo-splines is between B-spline and orthogonal refinable function of the same order. However, there is no precise regularity estimate for pseudo-splines in general. In this paper, an optimal estimate of the decay of the Fourier transform of the pseudo-splines is given. The regularity of pseudo-splines can then be deduced and hence, the regularity of the corresponding wavelets and framelets. The asymptotical regularity analysis, as the order of the pseudo-splines goes to infinity, is also provided. Furthermore, we show that in all tight frame systems constructed from pseudo-splines by methods provided both in [I. Daubechies, B. Han, A. Ron, Z. Shen, Framelets: MRA-based constructions of wavelet frames, Appl. Comput. Harmon. Anal. 14 (1) (2003) 1–46] and this paper, there is one tight framelet from the generating set of the tight frame system whose dilations and shifts already form a Riesz basis for $L_2(\mathbb{R})$.     

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.     

Characteristic roots and vectors of a diifferentiable family of symmetric matrices

LetA(x) be a differentiable family of k × k symmetric matrices where x runs through a domain D in RⁿWe prove that if λ is a continuous function onDsuch that, for every x εD,λ(x) is a characteristic root of A(x) of constant multiplicity m, then λ is a differentiable function and there exists, locally, a differentiable family of ortho-normal bases for the eigenspace. The case n = 1 has been known in the standard treatises on the perturbation theory for linear operators.     

Pairs of Dual Wavelet Frames from Any Two Refinable Functions

Starting from any two compactly supported refinable functions in L₂(R) with dilation factor d,we show that it is always possible to construct 2d wavelet functions with compact support such that they generate a pair of dual d-wavelet frames in L₂(R). Moreover, the number of vanishing moments of each of these wavelet frames is equal to the approximation order of the dual MRA; this is the highest possible. In particular, when we consider symmetric refinable functions, the constructed dual wavelets are also symmetric or antisymmetric. As a consequence, for any compactly supported refinable function in L₂(R), it is possible to construct, explicitly and easily, wavelets that are finite linear combinations of translates (d · – k), and that generate a wavelet frame with an arbitrarily preassigned number of vanishing moments.We illustrate the general theory by examples of such pairs of dual wavelet frames derived from B-spline functions.     

A new proof of some polynomial inequalities related to pseudo-splines

Pseudo-splines of type I were introduced in [I. Daubechies, B. Han, A. Ron, Z. Shen, Framelets: MRA-based constructions of wavelet frames, Appl. Comput. Harmon. Anal. 14 (2003) 1–46] and [Selenick, Smooth wavelet tight frames with zero moments, Appl. Comput. Harmon. Anal. 10 (2000) 163–181] and type II were introduced in [B. Dong, Z. Shen, Pseudo-splines, wavelets and framelets, Appl. Comput. Harmon. Anal. 22 (2007) 78–104]. Both types of pseudo-splines provide a rich family of refinable functions with B-splines, interpolatory refinable functions and refinable functions with orthonormal shifts as special examples. In [B. Dong, Z. Shen, Pseudo-splines, wavelets and framelets, Appl. Comput. Harmon. Anal. 22 (2007) 78–104], Dong and Shen gave a regularity analysis of pseudo-splines of both types. The key to regularity analysis is Proposition 3.2 in [B. Dong, Z. Shen, Pseudo-splines, wavelets and framelets, Appl. Comput. Harmon. Anal. 22 (2007) 78–104], which also appeared in [A. Cohen, J.P. Conze, Régularité des bases d'ondelettes et mesures ergodiques, Rev. Mat. Iberoamericana 8 (1992) 351–365] and [I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Series in Applied Mathematics, SIAM, Philadelphia, 1992] for the case l=N−1. In this note, we will give a new insight into this proposition.     

Band-limited Wavelets and Framelets in Low Dimensions

In this paper, we study the problem of constructing non-separable band-limited wavelet tight frames, Riesz wavelets and orthonormal wavelets in $\mathbb {R}^{2}$ and $\mathbb {R}^{3}$ . We first construct a class of non-separable band-limited refinable functions in low-dimensional Euclidean spaces by using univariate Meyer’s refinable functions along multiple directions defined by classical box-spline direction matrices. These non-separable band-limited definable functions are then used to construct non-separable band-limited wavelet tight frames via the unitary and oblique extension principles. However, these refinable functions cannot be used for constructing Riesz wavelets and orthonormal wavelets in low dimensions as they are not stable. Another construction scheme is then developed to construct stable refinable functions in low dimensions by using a special class of direction matrices. The resulting stable refinable functions allow us to construct a class of MRA-based non-separable band-limited Riesz wavelets and particularly band-limited orthonormal wavelets in low dimensions with small frequency support.     

三进制双正交对称小波的设计

本文给出了一种三进制双正交对称小波的设计方法.在给定插值紧支撑对称尺度函数的情况下,指出了如果对偶尺度函数同为紧支撑插值的,则它们同为1-型对称.并且给出了对偶尺度函数为紧支撑插值和非插值情况下的通解计算公式.还提出了频率优化方法设计对偶尺度函数和小波函数,把双正交条件归结为线性约束的二次规划问题,最后通过线性方程组来求解.对于小波函数本文也给出了一组特解公式.     

Orthogonality criteria for compactly supported refinable functions and refinable function vectors

A refinable function φ(x):ℝⁿ→ℝ or, more generally, a refinable function vector Φ(x)=[φ₁(x),...,φ_r(x)]^T is an L¹ solution of a system of (vector-valued) refinement equations involving expansion by a dilation matrix A, which is an expanding integer matrix. A refinable function vector is called orthogonal if {φ_j(x−α):α∈ℤⁿ, 1≤j≤r form an orthogonal set of functions in L²(ℝⁿ). Compactly supported orthogonal refinable functions and function vectors can be used to construct orthonormal wavelet and multiwavelet bases of L²(ℝⁿ). In this paper we give a comprehensive set of necessary and sufficient conditions for the orthogonality of compactly supported refinable functions and refinable function vectors.     

p-Adic Haar Multiresolution Analysis and Pseudo-Differential Operators

The notion of p-adic multiresolution analysis (MRA) is introduced. We discuss a “natural” refinement equation whose solution (a refinable function) is the characteristic function of the unit disc. This equation reflects the fact that the characteristic function of the unit disc is a sum of p characteristic functions of mutually disjoint discs of radius p ⁻¹. This refinement equation generates a MRA. The case p=2 is studied in detail. Our MRA is a 2-adic analog of the real Haar MRA. But in contrast to the real setting, the refinable function generating our Haar MRA is 1-periodic, which never holds for real refinable functions. This fact implies that there exist infinity many different 2-adic orthonormal wavelet bases in ℒ²(ℚ₂) generated by the same Haar MRA. All of these new bases are described. We also constructed infinity many different multidimensional 2-adic Haar orthonormal wavelet bases for ℒ²(ℚ₂ ⁿ ) by means of the tensor product of one-dimensional MRAs. We also study connections between wavelet analysis and spectral analysis of pseudo-differential operators. A criterion for multidimensional p-adic wavelets to be eigenfunctions for a pseudo-differential operator (in the Lizorkin space) is derived. We proved also that these wavelets are eigenfunctions of the Taibleson multidimensional fractional operator. These facts create the necessary prerequisites for intensive using our wavelet bases in applications. Our results related to the pseudo-differential operators develop the investigations started in Albeverio et al. (J. Fourier Anal. Appl. 12(4):393–425, 2006).      

Symmetric Framelets

We study tight wavelet frames associated with symmetric compactly supported refinable functions, which are obtained with the unitary extension principle . We give a criterion for the existence of two symmetric or antisymmetric compactly supported framelets. All refinable masks of length up to 6 satisfying this criterion are found.     

From symmetric subdivision masks of Hurwitz type to interpolatory subdivision masks

In this paper we present a general strategy to deduce a family of interpolatory masks from a symmetric Hurwitz non-interpolatory one. This brings back to a polynomial equation involving the symbol of the non-interpolatory scheme we start with. The solution of the polynomial equation here proposed, tailored for symmetric Hurwitz subdivision symbols, leads to an efficient procedure for the computation of the coefficients of the corresponding family of interpolatory masks. Several examples of interpolatory masks associated with classical approximating masks are given.     

Quasi-interpolatory refinable functions and construction of biorthogonal wavelet systems

We present a new family of compactly supported and symmetric biorthogonal wavelet systems. Each refinement mask in this family has tension parameter ω. When ω = 0, it becomes the minimal length biorthogonal Coifman wavelet system (Wei et al., IEEE Trans Image Proc 7:1000–1013, 1998). Choosing ω away from zero, we can get better smoothness of the refinable functions at the expense of slightly larger support. Though the construction of the new biorthogonal wavelet systems, in fact, starts from a new class of quasi-interpolatory subdivision schemes, we find that the refinement masks accidently coincide with the ones by Cohen et al. (Comm Pure Appl Math 45:485–560, 1992, §6.C) (or Daubechies 1992, §8.3.5), which are designed for the purpose of generating biorthogonal wavelets close to orthonormal cases. However, the corresponding mathematical analysis is yet to be provided. In this study, we highlight the connection between the quasi-interpolatory subdivision schemes and the masks by Cohen, Daubechies and Feauveau, and then we study the fundamental properties of the new biorthogonal wavelet systems such as regularity, stability, linear independence and accuracy.     

Construction of two-direction tight wavelet frames

We investigate the construction of two-direction tight wavelet frames First, a sufficient condition for a two-direction refinable function generating two-direction tight wavelet frames is derived. Second, a simple constructive method of two-direction tight wavelet frames is given. Third, based on the obtained two-direction tight wavelet frames, one can construct a symmetric multiwavelet frame easily. Finally, some examples are given to illustrate the results.