Some simple Haar-type wavelets in higher dimensions

An orthonormal wavelet system in R-d, d is an element of N, is a countable collection of functions {psi(l)(j,k)} j is an element of Z, k is an element of Z(d), l = 1,...L, of the form psi(l)(j,k)(x) = \det a\(-j/2)psi(l)(a(-j)x-k) = (D-a(j) T-k psi(l))(x) that is an orthonormal basis for L-2 (R-d), where a is an element of GL(d) (R) is an expanding matrix. The first such system to be discovered (almost 100 years ago) is the Haar system for which L = d = 1, psi(1) (x) = psi(x) = chi([0, 1/2))(x) - chi([1/2, 1))(x), a = 2. It is a natural problem to extend these systems to higher dimensions. A simple solution is found by taking appropriate products Phi(x(1), x(2),..., x(d)) = phi(1) (x(1))phi(2) (x(2))...phi(d) (x(d)) of functions of one variable. The obtained wavelet system is not always convenient for applications. It is desirable to find "nonseparable" examples. One encounters certain difficulties, however, when one tries to construct such MRA wavelet systems. For example, if a = ((1-1)(1)(1)) is the quincunx dilation matrix, it is well-known (see, e.g., [5]) that one can construct nonseparable Haar-type scaling functions which are characteristic functions of rather complicated fractal-like compact sets. In this work we shall construct considerably simpler Haar-type wavelets if we use the ideas arising from "composite dilation" wavelets. These were developed in [7] and involve dilations by matrices that are products of the form a(j) b, j is an element of Z, where a is an element of GL(d) (R) has some "expanding" property and b belongs to a group of matrices in GL(d) (R) having \det b\ = 1.