arxiv:math/ v1 [math.ca] 19 May 2004

Transcription

1 Fractal Components of Wavelet Measures arxiv:math/040537v [math.ca] 9 May 004 Abstract Palle E.T. Jorgensen Department of Mathematics, The University of Iowa, 4 MacLean Hall, Iowa City, IA, 54-49, U.S.A. We compute the Coifman-Meyer-Wickerhauser measure µ for certain families of quadrature mirror filters (QMFs), and we establish that for a subclass of QMFs, µ contains a fractal scale. Key words: Hilbert space, Cuntz algebra, map, wavelet packets, pyramid algorithm, product measures, orthogonality relations, equivalence of measures, iterated function systems (IFS), scaling function, multiresolution, subdivision scheme, singular measures, absolutely continuous measures 000 MSC: 4C40; 4A6; 43A65; 4A65 Introduction Wavelets and wavelet packets are orthonormal bases, or Parseval frames, for the Hilbert space L (R) which are built from a small set of generating functions by applying only two operations to them, translation by the integers Z R, and by dyadic scaling t j t, j Z. The best known case is that of singly generated wavelets. Then there is a function ψ L (R) such that the family { j/ ψ ( j t k )} j,k Z (.) is an orthonormal basis for L (R). Wavelet packets (in the orthonormal case) are sequences (ϕ n ) n N0 and subsets J N 0 Z such that the two families { j/ ϕ n ( j t k ) (n, j) J, k Z } (.) address: jorgen@math.uiowa.edu (Palle E.T. Jorgensen). URL: jorgen (Palle E.T. Jorgensen. This material is based upon work supported by the U.S. National Science Foundation under Grant No. DMS (FRG). Preprint submitted to Elsevier Science February 008

2 and {ϕ n (t k) n N 0, k Z} (.3) are both both orthonormal bases. In the seminal paper [4], the authors Coifman, Meyer, and Wickerhauser proposed a selection of systems (.) from (.3); the idea being that (.3) may be constructed by an effective and relatively simple algorithm. The construction starts with a quadrature-mirror wavelet filter, i.e., a Fourier series m 0 (x) = k Za k e iπkx (.4) such that a k a k+l = δ 0,l (.5) k Z and a k =. (.6) k Z If the expansion in (.4) is a finite sum, then the functions ϕ 0, ϕ, ϕ, will be of compact support. The advantage with using (.) over (.) or (.3) is that the basis functions in (.) are better localized in time-frequency, as is spelled out in [4] and []. The authors of [4] propose a certain integral decomposition of L (R) which is based on a certain measure µ 0 on the unit-interval [0, ]. The idea is the following: The quadrature rules (.5) lead to a dyadic decomposition of the Hilbert space l (Z), repository for the wavelet coefficients: For each k N, and each dyadic interval [ ) ξ, ξ + k, ξ = i + + i k k, (.7) we assign a closed subspace H (ξ) l, such that the k distinct subspaces H (ξ) are mutually orthogonal, and H (ξ) = l (.8) ξ If e 0 denotes the vector e 0 (j):= δ 0,j in l, and P (ξ) the orthogonal projection onto H (ξ), then the measure µ 0 is determined by where ξ is the dyadic rational (.7). µ 0 ([ ξ, ξ + k )) = P (ξ)e 0 (.9) More generally, if f H, f =, we set µ f ( ) = P ( )f, and note that each µ f is a probability measure. These measures are shown in section 3 to dictate the wavepacket analysis in the Hilbert space L (R).

3 The expectation is that µ 0 is absolutely continuous with respect to the Lebesgue measure on [0, ]. We give a formula for µ 0 directly in terms of the coefficients (a j ) j Z in (.4), and we give a one-parameter family of coefficients a j (β) where β parameterizes the circle, and where for β < π/4, the measure µ 0 on small dyadic intervals J is governed by the formula where J denotes the length of J, and where µ 0 (dyadic interval J) J s, (.0) s = ln ( a 0 (β) ). (.) ln In the interval for β, the coefficients a j (β) are real-valued, and when β < π/4, a 0 (β) > /. Fourier Polynomials Let T = {z C z = }, and let D N. Consider the following two functions and m 0 (z) = D k=0 a k z k (.) We shall need several lemmas. m (z) = z D m 0 ( z) (.) = D k=0 ( ) k ā D k z k Lemma. The finite sequence a 0, a,, a D in (.) satisfies ā k a k+l = δ 0,l (.3) j Z if and only if anyone of the following equivalent conditions holds: (a) The matrix function U (z) = D k=0 A k z k, z T (.4) 3

4 is unitary, where the coefficient matrices are a k a k+ A k = for k = 0,,, D. (.5) a (D k) a (D k ) (b) The operators (S i f)(z) = m i (z) f ( z ), f L (T), i = 0,, (.6) and their adjoints Si satisfy Si S j = δ i,j I, i, j = 0,. (.7) S 0 S0 + S S = I (c) The following matrix function M (z) : = m 0 (z) m 0 ( z), z T, (.8) m (z) m ( z) is unitary, i.e., M (z) M (z) = I, for all z T. Remark. In the summation (.3), it is understood that the terms are zero if the summation index is not in the range which is specified for the sequence: If for example D =, then the conditions (.3) spell out as follows: a 0 + a + a + a 3 =. (.9) ā 0 a + ā a 3 = 0 The interpretation of (.4) is that the matrices A 0, A, A D are Fourier coefficients for a matrix function z U (z) defined on T. The requirements is that for each z T, the by matrix U (z) is unitary, i.e., that U(z) U (z) = I, z T. (.0) In general, for a function of the form (.4), the unitarity condition (.0) is equivalent to the following conditions on the constant matrices A 0, A, : A k+n A k = δ 0,n (.) or equivalently k A k A k+n = δ 0,n. (.) k 4

5 Again if D =, these conditions flesh out as follows: The two matrices A 0 = a 0 a a 3 a and A = a a 3 a a 0 (.3) satisfy A 0 A 0 + A A = I A A 0 = 0. (.4) PROOF. of Lemma.. The coefficients a 0, a,, a D are given and the two functions m 0 (z) and m (z) are defined as in (.) and (.). Even though they are defined for all z C, we shall use only the restrictions to T. It is convenient in some later calculations to introduce the substitution z:= e iπx ; and with this substitution we shall also view m 0 and m as oneperiodic functions on the real line R, and make the identification T = R/Z. We now turn to the implications: (.3)= (a). We already noted that the unitarity of the function U(z) in (.4) may be expressed by the conditions (.). When (.5) is substituted in (.), it is immediate that the two identities (.3) and (.) are equivalent. (a)= (b). Using the normalized Haar measure µ on T and the corresponding inner product f g : = f (z)g (z) dµ (z) (.5) of the Hilbert space L (T), we get the formula (S i f)(z) = w T,w =z As a result, we get the formula T m i (w)f (w), defined for f L (T), z T, and i = 0,. (S i S jf)(z) = w T,w =z (.6) m i (w)m j (w)f (z) ; (.7) which states that the four operators S i S j are multiplication operators. Hence the first part of (.7) reads: w T,w =z m i (w)m j (w) = δ i,j, z T. (.8) Substitution of (.) into (.8) shows that these conditions amount to the same sum rules (.3), or equivalently (.). 5

6 The second equation from (.7) may be restated as or equivalently S 0f + S f = f, for all f L (T) ; (.9) m i (w)f (w) dµ (z) = f (z) dµ (z). (.0) i=0 w =z T But the summation on the left-hand sides fleshes out as follows: and 4 i=0 w =ζ =z = = = T w =ζ =z w =ζ =z w =z w =z m i (w)m i (ζ)f (w)f(ζ) ( ) m i (w)m i (ζ) f (w)f(ζ) i=0 δ w,ζ f (w)f(ζ) f (w), f (w) dµ (z) = T f (z) dµ (z) (.) by an elementary property of the Haar measure. This proves (a)= (b); and in fact equivalence. (b)= (c). This implication also goes in both directions, and it is implicit in the previous step where we showed that the identities (.8) are equivalent to (.3). But an inspection of matrix entries shows that (.8) is a restatement of the unitary property for the matrix M (z) in (c). Corollary.3 Consider a Fourier polynomial m (z) = D k=0 a k z k (.) as in Lemma., and let S = S m be the operator on L (T) defined by Then the following conditions are equivalent: (i) S is an isometry; (ii) m (z) + m ( z) =, for z T; (Sf)(z) = m (z) f ( z ), z T. (.3) 6

7 and (iii) k a k a k+l = δ 0,l. PROOF. The equivalence (i) (ii) is immediate from an inspection of (.8) in the special case i = j = 0. To see (ii) (iii), it is simplest to write out the Fourier coefficients of the function z w T, w =z A substitution of (.) into (.4) yields a j a k w k j j k w =z = a j a k w n n j,k s.t. k j=n w =z = n j a j a j+n z n, m (w). (.4) and it follows that the sum in (.4) is the constant function if and only if (iii) is satisfied. Note that (iii) is a restatement of condition (.3) in Lemma.. In the following we consider operators S on the form (Sf)(z) = m (z) f (z ), z T, f L (T), and their adjoints. But we will need these operators realized explicitly in the sequence space l = L (T). Lemma.4. (a) The matrix representation for the operator S in Corollary.3 relative to the standard basis e n (z):= z n, n Z, for L (T), is (Sξ) n = k a n k ξ k ; (.5) and the adjoint S is (S ξ) n = k a k n ξ k. (.6) (b) The D-dimensional subspace L spanned by {e k D + k 0} is 7

8 invariant under S, and the corresponding D by D matrix is ā ā ā ā ā 3 ā. ā D ā D ā D.. ā (.7) ā ā ā ā 3 ā ā... ā D ā D ā D ā D ā D Specifically, the case D = is ā ā ā ā 0 0 ; (.8) 0 ā 3 ā ā ā 3 and the space L is then spanned by {, z, z, z 3 }. (c) For all n Z, there is a k N such that S k e n L. (.9) PROOF. The matrix representations (.5), and (.6), follow directly from the formulas (.3) and (.6). To prove (b), set so that We claim the following implication: J (D) : = {k Z D + k 0}, L = span {e k k J (D)}. ξ k = 0 for k Z\J (D) = k ā k n ξ k = 0 for n Z\J (D). (.30) 8

9 The invariance of L under S is immediate from this. Now suppose ξ k = 0 for k Z\J (D). Then one of the factors in ā k n ξ k vanishes if n Z\J (D); and the implication (.30) follows. To understand (c), iterate the matrix formula (.6). We get S k e n = p ā i ā i ā ik e p (.3) i,i,,i k where the summation on the right-hand side is over p Z, i, i, [0,,, D ], subject to i + i + + k i k + p k = n, so the range for the p-index is roughly divided by with each iteration of S on the basis vector e n (z) = z n. Remark.5 Conclusions (a) (b) in Lemma.4 hold for the smaller subspace M:= span {e k D k 0}. The matrix of the restricted operator S M is then obtained from (.7) by deletion of the last row and last column. For example, the reduced matrix corresponding to (.8) is ā F: = ā ā ā 0. 0 ā 3 ā However, property (c) in Lemma.4 is not satisfied for M. To see this, note that D =, and L = M Ce 3 in the example. As a result e 3 S k e 3 = (ā3 ) k for all k N. (.3) If a 3 0, it follows that S k e 3 is not in M, no matter how large k is. 3 Measures Induced by Representations of the Cuntz Relations The system of coefficients studied in section give rise to representations of the Cuntz relations. We recall the construction briefly; see [] for more details. When the coefficients a 0, a, are given, satisfying (.3), then the two operators S 0 and S from Lemma. (b) satisfy the relations (.7). We say that the two operators satisfy the Cuntz relations, or equivalently that they define a representation of the Cuntz algebra O on the Hilbert space H = L (T). 9

10 The norm-closed algebra of operators on H generated by S 0, S, S 0, and S is known, [5] to be a simple C -algebra; and it is unique up to isomorphism of C -algebras. We will also need to Hilbert space L (R). If in addition to (.3), the numbers (a k ) satisfy D k=0 a k =, (3.) then it is known [4] that the following system of three equations has solutions in L (R): ϕ n (x) = a k ϕ n (x k), (3.) k ϕ n+ (x) = k ( ) k a D k ϕ n (x k), (3.3) and R ϕ 0 (x) dx =. (3.4) Moreover, for every f L (R), the following Parseval identity holds: When a system, such as ϕ n ( k) f = f = f (x) dx. (3.5) n=0 k Z R {ϕ n ( k) n N 0, k Z}, (3.6) satisfies (3.5) we way it is a Parseval frame, or a normalized tight frame. There are simple further conditions on the coefficients (a k ) each of which guarantees that (3.6) is in fact an orthonormal basis (ONB) in L (R), i.e., that ϕ n ( k) ϕ m ( l) = δ n,m δ k,l. (3.7) To understand how to select subsets J N 0 Z such that the set of functions F (J): = { j ϕn ( j x k ) (n, j) J, k Z } (3.8) forms a Perseval frame, or an ONB, for L (R), the authors of [4] and [] suggested a family of measures associated with the operator systems (.7). We proved more generally, in [0] that these measures are associated with any representation of one of the Cuntz algebras O N, N =, 3,. The idea is simple: Let (S i ) i=0 be two operators in a Hilbert space H such that Si S j = δ i,j I and S i Si = I (3.9) i=0 0

11 hold. Then for every multi-index ξ:(i, i,, i k ), i j {0, } (3.0) the operator satisfies P (ξ): = S i S ik S i k S i (3.) P (ξ) = P (ξ) = P (ξ). Moreover P (ξ) = I (3.) ξ and P (ξ)p (η) = δ ξ,η P (ξ) (3.3) where the summation in (3.) is over all multi-indices of length k, and in (3.3) the multi-indices ξ and η have the same length. We say that, for each k N we have a dyadic partition of the Hilbert space H into orthogonal frequency bands. Theorem 3. Let m 0 :T C be a function which satisfies m 0 () = and m 0 (z) + m 0 ( z) = for z T. (3.4) Set m (z) = z m 0 ( z), for z T; (3.5) and let S i, i = 0,, be the operators (.6). Suppose the sequence ϕ 0, ϕ, ϕ, defined by ϕ n (x) = ( ) ( ) x x m 0 ϕ n (3.6) ϕ n+ (x) = ( ) ( ) x x m ϕ n (3.7) for n = 0,,, x R defines an orthonormal basis in L (R), where we set m j (x):= m j (e iπx ), j = 0,. Then p ϕn ( p t k) = ej Si p Si e k ϕ p n+i +i + +i p p (t j) i =0 i p=0 j Z (3.8) holds for all p, n N 0 and k Z. Moreover a subset J N 0 N 0 has the corresponding set { p ϕn ( p t k) (p, n) J, k Z } (3.9) define an orthonormal basis (ONB) in L (R) if and only if the sets [ p n, p (n + )) form a non-overlapping partition of N 0.

12 PROOF. Since the functions in (3.6) (3.7) form an orthonormal basis for L (R), we have ϕ n (x + j) ϕ n (x + j) = δ n,n for n, n N 0, (3.0) j Z and (3.8) will follow if we check that ϕ m (t j) p ϕn ( p t k) dt = e j Si p Si e k R (3.) if m = p n + i + i + + i p p, and zero otherwise. But (3.) follows from Fourier duality, and a substitution of formulas (3.6) (3.7) and (3.0). We leave the computation to the reader. The second conclusion regarding the basis properties (ONB) of the functions in (3.9) follows from (3.8) when we note that m = p n + i + i + + i p p ranges over [ p n, p (n + )) as i,, i p vary over Z Z }{{. } p times Example 3. The best known special case of (3.9) is when J = {(0, ),(, ),(, ),(3, ), }, and the corresponding partition of N 0 is simply [ p, p+ ) for p = 0,,,. In this case, we only need the first two functions from the sequence ϕ n of (3.6) (3.7), i.e., ϕ 0 = ϕ, and ϕ = ψ. The other extreme is J = {(0, ), (0, ),(0, 3), }. In this case, each section in the corresponding partition of N 0 is a singleton, and all the functions {ϕ n n N 0 } from (3.6) (3.7) are needed in the (3.9) family. However, no dyadic scales are used. Between these two extremes there is a countable infinite family of non-overlapping partitions corresponding to the choice specified in (3.9); for example, the pairs (p, n) in some J from (3.9) may be specified as in the following table (Table ). Table. A set J consisting of points (p, n) as follows: p n In this case, we are using only the four functions ϕ 0, ϕ, ϕ, and ϕ 3 from (3.6) (3.7), and the segments from (3.9) are {0,,, 3} followed by [ k j, k (j + ) )

13 where k =,, 3,, and j {,, 3}. As a result, all the dyadic scaling numbers k are used. The resulting ONB in (3.9) consists of the following functions, ϕ 0 ( t l 0 ), k ϕ j ( k t l ), where l 0, l Z, k =,,, and j {,, 3}. Since {ϕ m ( j) m N 0, j Z} is an ONB for L (R), every f L (R), f =, defines a probability distribution P f on N 0 Z by P f (m, j) : = ϕ m ( j) f. (3.) Corollary 3.3 Let p, n N 0, k Z, and let P p (, ) be the corresponding probability distribution from (3.) ϕ n( p k) above. Then the marginal distribution in the first variable is ([ δ m, p n+i +i + +i p p µ i e k + + i p i p, + + i p + )). (3.3) p p PROOF. In view of (3.) and Theorem 3., the marginal distribution is j Z P p (m, j) ϕ n( p k) = δ m, p n+i +i + +i p p ej Si p Si e k j Z S = δ m, p n+i +i + +i p p ip Si e k ([ = δ m, p n+i +i + +i p p µ i e k + + i p i p, + + i p + )) p p which is the desired conclusion. The next result follows from [4], or from [0]. Lemma 3.4 There is a unique Borel measure E defined on the unit interval [0, ] such that E ([ i + + i k i k, + + i k + )) = P (i k k, i,, i k ) = S ξ Sξ. (3.4) Moreover the measure E = E P takes values in orthogonal projections in H, and E (dx) = I (3.5) 0 E (B ) E (B ) = 0 if B i are Borel sets with B B =. (3.6) 3

14 Remark 3.5 The expectation is that the measure class of E is that of the Lebesgue measure if E is derived from a wavelet representation. For the wavelet representations, the Hilbert space H is L (T); and for a generic set of wavelet representations e 0 is a cyclic vector; see [9]. Hence, to check if E is in the Lebesgue class, it is enough if we verify that the scalar valued, positive measure µ 0 ( ): = e 0 E ( )e 0 = E ( )e 0 (3.7) is absolutely continuous with respect to Lebesgue measure on [0, ]. Using (3.4), we see that the measure µ 0 is given on dyadic intervals by the formula ([ i µ i k i k, + + i k + )) = S k ξe k 0. (3.8) Part of the reason for expecting that the measure µ 0 ( ) = E ( )e 0 should be absolutely continuous with respect to Lebesgue measure on [0, ] rests on the observation that λ = is always in the spectrum of the matrix (.7). If λ = is also a dominant eigenvalue, in a sense which we make precise below, then the asymptotic formula ([ i µ i k i k, )) k (3.9) k k would follow. From this the absolute continuity follows as well. But a closer examination of the matrices F i, i = 0,, shows that λ = is not always a dominant point in the spectrum of F 0 = S0 M. We now show that is in the spectrum. Lemma 3.6 Let the numbers a 0, a,, a D satisfy (.3) and (3.), and let S 0 and S 0 be the corresponding operators on l ; see (.5) and (.6). Then λ = is an eigenvalue for the finite-dimensional operator F 0 = S0 M. If F 0 is represented as a D by D matrix as in (.7), then the row-vector w = (,,, ) is a left-eigenvector for F 0, or equivalently w is }{{} D a right-eigenvector, with the same eigenvalue, for the adjoint matrix F0. PROOF. Substitution of m 0 () = into m 0 (z) + m 0 ( z) = from Corollary.3 yields m 0 ( ) = 0. Hence, z + divides the polynomial m 0 (z). 4

15 If m 0 (z) = (z + ) L (z), and L (z) = L () = D i=0 D i=0 l i z i, then it follows that l i = =. As a result we get a 0 = l 0 a i = l i + l i for i D a D = l D. (3.30) It follows from this that and D i=0 D i=0 a i = i a i+ = i l i = l i =. Working out the matrix product, this is a restatement of wf 0 = w, or equivalently F0 w = w. Theorem 3.7 Let a 0, a, a D be as stated in Lemma 3.6, suppose a 0 > max { λ λ spec (F 0 ) \ {a 0 }}, and set s = ln a 0 ln. Then every non-empty open subset V in [0, ] contains an infinite sequence of dyadic intervals J such that 0 < lim inf J V ( ) µ0 (J) J s lim sup J V ( ) µ0 (J) J s <. We say that the measure µ 0 contains s as a fractal scale. PROOF. The details follow by combining Lemma 3.6 with the technical lemma which we state and prove in Section 5 below. 5

16 Proposition 3.8 Let the numbers a 0,, a D satisfy the two conditions (.3) and (3.) and let S 0, S be the corresponding two operators; see (.6) (.7). Let µ 0 ( ):= E ( )e 0 be the measure (3.8). Then ([ i µ i k i k, + + i k + )) a k k 0 #(i=0) a D #(i=). (3.3) Remark 3.9 We will see later how this estimate yields information about a fractal component of µ 0 in the special case D =. PROOF. Let F i := Si, i = 0,. For k N, and i,, i k {0, }, set M ξ = i + + i k k and F ξ := F ik F i. Then F 0 e 0 = ā 0 e 0 and F e 0 = ā D e 0. To see this, recall that and m (z) = ā D ā D z + ā 0 z D, ([ µ )) 0 ξ, ξ + k = F ξ e 0 e 0 F ξ e 0 = F i Fi k e 0 e 0 = a 0 #(i=0) a D #(i=). Lemma 3.0 (The case D =.) (a) For the case D =, the real valued solutions a 0, a, a, a 3 to the three equations a 0 + a + a + a 3 = a 0 a + a a 3 = 0 (3.3) a 0 + a + a + a 3 = 6

17 are for β R. ( ) a 0 = + cosβ ( ) a = + sin β a = a 3 = ( cosβ ) ( sin β ) (3.33) (b) The spectrum of the matrix a F = a a a 0 0 a 3 a (3.34) obtained by restriction of S to M = span {e 0, e, e } is { a 0,, } sin β cosβ. (3.35) (c) The following two conditions are equivalent a 0 > max { λ λ spec (F) \ {a 0 }} (3.36) and β < π 4. (3.37) (d) The eigenspace for the point a 0 in the spectrum of F is spanned by the vector ( ) v where v R is the solution to v = a 0 I a a 0 a 3 a a. (3.38) 0 (e) For the case of Daubechies s wavelet, we have a 0 = a = a = a 3 = 3 4, (3.39) 7

18 and ( ) 3 β = arccos, (3.40) spec (F) = { a 0,, }. (3.4) In particular, a 0 < = max { λ λ spec (F)}. (3.4) PROOF. Part (a) follows directly from section 4 in []. To compute the spectrum of the matrix F, note that its characteristic polynomial is (λ a 0 ) ( λ (a + a )λ + a a a 0 a 3 ) ; (3.43) and when (3.33) are substituted, we see that the roots are as listed in (3.35). The conclusions listed in parts (c) and (d) are immediate consequences of (b). Finally (3.38) follows from the observation that F has the form F = a a G 0 (3.44) with G = a a 0 and the fact that a 0 is not in the spectrum of G. a 3 a The conclusions listed in (e) for Daubechies s wavelet follow by an application of (a) (d). To better appreciate the geometry of the formulas (3.33), the reader may find Fig. useful. 8

19 The coefficients a 0 (β), a (β), a (β), a 3 (β) of (3.3) a 3 (β) a (β) β = 3π 4 β = π 4 β = π β = 3π 4 β = 0 β < π 4 β = π 4 a 0 (β) β = cos ( ) 3 Daubechies a (β) Figure. The four coefficients a 0, a, a, a 3 are described by the two circles in Fig., viz., by ( a 0 ) +( a 3 ) =, and by 4 ( a ) +( a ) = 4. (3.45) This representation also makes it clear that, if one of the four coefficients vanishes, then so does a second from the remaining coefficients. The resulting four degenerate cases are as sketched in Table. Table. The degenerate cases when one (and therefore two) of the coefficients a 0, a, a, a 3 vanish: m 0 +z +z 3 z+z z +z 3 m z z 3 z 3 z+z z β π 4 π 4 3π 4 3π 4 There are four distinct configurations, and all four are variations of the Haar wavelet, see Fig.. The second column in Table is called the stretched Haar wavelet; see []. For all three, we have the formula ([ i µ i k i k, + + i k + )) = k (3.46) k k 9

20 which proves that µ 0 is the Lebesgue measure on the unit-interval [0, ] for all four variations of the Haar wavelet. The two matrices F 0 and F corresponding to the four variants of the Haar wavelet filters in Table are as follows: F 0 : F : The father function ϕ and the mother function ψ for the Haar Wavelet Figure. 0

21 Remark 3. We shall need Lemma 3.6 for the two matrices a a F 0 = a a a 0 and F = a a a 3 0 a 3 a 0 a 0 a 3 (3.47) For the range β < π 4, we have (3.36) satisfied by spec (F 0). But then the point a 3 = a 3 (β) satisfies a 3 <, and = max { λ λ spec (F )}. (3.48) Theorem 3. Let β R, β < π, be given. Let a 4 0, a, a, a 3, be the corresponding numbers in (3.33). Let (S i ) i=0 be the operators in (.6), and F i : = S i (3.49) M the matrices (3.47). Let be a dyadic fraction. ξ = i + + i k k (3.50) Then lim n a n 0 F 0 n F i k F i e 0 = a #(i=0) 0 a #(i=) 3 (e 0 + v) (3.5) where v is the vector in (3.38). Moreover, if ε > 0 is given, there is an n 0 N such that ([ µ 0 ξ, ξ + )) a n a #(i=0) n+k 0 0 a #(i=) ( 3 + v ε ). (3.5) for all n n 0. More generally, given ε, 0 < ε <, there is n 0 N such that ([ )) ( + v ε ) ( a #(i=0) 0 a #(i=) ) µ 0 ξ, ξ + 3 n+k a n 0 ( a #(i=0) 0 a #(i=) 3 ) ( + v + ε ) for all n n 0. PROOF. Recall that e 0 + v is the eigenvector satisfying F 0 (e 0 + v) = a 0 (e 0 + v), and e 0 e 0 + v =. (3.53)

22 Moreover, we have Using (3.36) for F 0, we get F 0 e 0 = a 0 e 0 and F e 0 = a 3 e 0. (3.54) lim n a n 0 F n 0 w = e 0 w (e 0 + v) (3.55) for all w M. Applying this, and (3.54), to w:= F ik F i e 0, the conclusion (3.5) follows. Since ([ µ 0 ξ, ξ + )) = S n n+k 0 S ξ e 0 = F n 0 F i k F i e 0, (3.56) the second conclusion (3.5) follows from the first. Corollary 3.3 Let β < π 4, and let a 0, a, a, a 3 be the corresponding numbers in (3.33). If the measure µ 0 ( ) = E ( )e 0 is absolutely continuous with respect to Lebesgue measure dx on [0, ]. Then the Radon-Nikodym derivative dµ 0 /dx = f is unbounded in every open (non-empty) subset of [0, ]. PROOF. Let V [0, ] be a non-empty open subset. Then pick ξ as in (3.50), and n N such that [ ξ, ξ + ) V. (3.57) n Then pick ε, 0 < ε <, and n 0 n such that (3.5) holds for all n n 0. If f were bounded in V with upper bound C; then (3.5) implies the estimate C n a n 0 a #(i=0) 0 a #(i=) ( 3 + v ε ) (3.58) for all n n 0. Since a 0 >, lim n n a n 0 =, which contradicts (3.58). Example 3.4 We introduced the constants a 0, a, a, a 3 for Daubechies s wavelet in part (e) of Lemma 3.6. An inspection of (3.33) and (3.39) shows that β = arccos ( ) 3. Hence sin β cosβ =, and the three numbers in the spectrum of F 0 are as follows: < a 0 <. (3.59)

23 One checks that the vector v : = ( + 3 ) e + ( 3 ) e (3.60) spans the -eigenspace for F 0. Since a 0 <, the eigenspace { w M F 0 w = } w (3.6) is orthogonal to e 0. Corollary 3.5 Let β = arccos ( ) 3, so that the numbers a0, a, a, a 3 define the Daubechies wavelet. Let v be the eigenvector in (3.60), and pick w M such that F 0 w = w, and w v = 8. (3.6) Let ξ be a dyadic rational given as in (3.50). Then, for every ε, 0 < ε <, there is n 0 N such that w F ik F i e 0 ε µ ([ )) 0 ξ, ξ + n+k (3.63) n w F ik F i e 0 + ε for all n n 0. PROOF. The argument follows the proof of Theorem 3.. The vector w is w = (e 0 + e + e ). The main difference is that now is the top eigenvalue of F 0 with eigenvector v. If w is chosen as in (3.6), then lim n n F n 0 s = 8 w s v (3.64) for all s M. The desired result (3.63) follows when this is applied to s = F ik F i e 0. 4 Wavelets on Fractals In section we proved some lemmas for representations of the Cuntz algebra O, and we used them in section 3 above in our analysis of dyadic wavelets in the Hilbert space L (R). But there are other Hilbert spaces H which admit wavelet algorithms. In a recent paper [6], we constructed wavelets in separable Hilbert spaces H (s) built on Hausdorff measure (dx) s where s denotes the 3

24 corresponding fractal dimension. For the middle third Cantor set, for example, the scaling number N is 3; and the fractal dimension is s = ln ln 3 = log 3(). In this section, we show how our Cantor subdivision construction is built on a representation of the O 3 -Cuntz relations on L (T) which are analogous to the representations of O which we used in our analysis of dyadic L (R)-wavelets. Lemma 4. Let N N, N, be given, and let m 0, m,, m N be bounded measurable functions on T. Then the operators (S i f)(z) = m i (z) f ( z N), f L (T), i = 0,,, N, z T (4.) satisfy the Cuntz relations Si S j = S i,j I N S i Si = I i=0 (4.) if and only if the associated N by N matrix function M N (z) = ( ( )) mj ze i kπ N N j,k (4.3) takes values in the unitary matrices for a.e. z T. PROOF. Rather than sketching the details, we will instead refer the reader to section above where the argument is done in full for N = ; see also [6] or []. For each of the representations (S i ) N i=0, we get a measure µ 0 on the unitinterval [0, ], as outlined in section 3 (for the special case N =.) In the general case, the measure µ 0 is determined uniquely on the N-adic subintervals as follows, ([ i µ 0 N + + i k N, i k N + + i k N + )) = S k N k i k Si Si e 0 where i, i, {0,,, N }. (4.4) We state the next result just for the middle-third Cantor set; but, following the discussion in [6], the reader will convince him/herself that it carries over to any fractal constructed by iteration of finite families of (contractive) affine maps in R d ; so called iterated function systems (IFS). 4

25 Proposition 4. Let m 0 (z) = +z m (z) = z (4.5) m (z) = z Then the unitarity condition (4.3) from Lemma 4. is satisfied, and the operators (S i f)(z) = m i (z) f ( z 3), f L (T), i = 0,,, z T, (4.6) define a representation of O 3 on L (T); and the corresponding measure µ 0 from (4.4) is the Hausdorff measure µ (s) of Hausdorff dimension s = ln = ln 3 log 3 () restricted to the middle-third Cantor set X 3 [0, ]. Remark 4.3 A classical theorem of Hutchinson [8], see also [7], implies that the Cantor measure, supported on X 3, is the unique Borel probability measure µ which solves f (x) dx (x) = ( ( ) ( ) ) x x + f dx (x) + f dµ (x) 3 3 for all continuous bounded functions f. (4.7) As a result, we need to verify that our measure µ 0 from (4.4), and from the O 3 -representation, satisfies identity (4.7). To compute the terms on the right-hand side in (4.4) is a finite matrix problem. We need only to calculate the operators Si, i = 0,,, on the threedimensional subspace M:= span {e 0, e, e }. But the argument from section 3 shows that the three restricted operators F i : = Si, i = 0,, (4.8) M have the following matrix representation F 0 = 0 0, F = 0 0, F = (4.9) A substitution of (4.8) (4.9) into (4.4) yields ([ i µ 0 N + + i k i N k, N + + i k N + )) k N k = F ik F i F i e 0 0 if one or more of the i j s is = k otherwise. (4.0) 5

26 Using this, it is immediate to see that µ 0 satisfies Hutchinson s identity (4.7), and therefore µ 0 is the Hausdorff measure µ (s), s = log 3 (), restricted to X 3. 5 A Technical Lemma In the proof of Lemma 3.6, Theorem 3. and Corollary 3.3 above, we relied on the following lemma regarding operators in a finite-dimensional Hilbert space. While it is analogous to the classical Perron-Frobenius theorem, our present result makes no mention of positivity. In fact, our matrix entries will typically be complex. Notation 5. If M is a complex Hilbert space, we denote by L (M) the algebra of all bounded linear operators on M If M is also finite-dimensional, we will pick suitable matrix representations for operators F:M M If M contains subspaces, M i, i =, such that M M and M = M M, then we get a block-matrix representation F = A B (5.) C D where the entries are linear operators specified as follows. A:M M, B:M M ; and C:M M, D:M M. If dim M =, and M = Cw for some w M, then we will identify the operators M M with M via where η M. The adjoint operator is T η :C z zη, Tη x = η x w, for x M. Lemma 5. Let M be a finite-dimensional complex Hilbert space, with d = dim M. Let F L (M), and let a C satisfy the following four conditions: (i) a spec (F); 6

27 (ii) a > max { λ λ spec (F) \ {a}}; (iii) the algebraic multiplicity of a is one; and (iv) there is a w M, w =, such that F w = āw. Then there is a unique ξ M such that w ξ = and Fξ = aξ. (5.) Moreover, lim n a n F n x = w x ξ for all x M. (5.3) Lemma 5.3 There is a constant C independent of d = dim M and of x, such that { a n F n x w x ξ Cn d s n } max s spec (F) \ {a}. (5.4) a PROOF. [Lemma 5.] Set Then M : = M Cw = {x M w x = 0}. M = Cw M, (5.5) and we get the resulting block-matrix representaton of F, a 00 0 F = (5.6) η G where a is the number in (i), the vector η M, and operator G L (M ), are uniquely determined. As a result, we get the factorization det (λ F) = (λ a) det (λ G) (5.7) for the characteristic polynomial. Assumptions (ii) and (iii) imply spec (F) \ {a} = spec (G) ; (5.8) and in particular, we note that a is not in the spectrum of G. Hence the inverse (a G) is well defined, and (a G) L (M ). 7

28 We claim that the vector satisfies the conditions in (5.). ξ = w + (a G) η (5.9) First note that (a G) η M, so w ξ = w w = w =. Moreover, Fξ = aw + η + G (a G) η = aw + a (a G) η = aξ, which proves the second condition in (5.). Uniqueness of the vector ξ in (5.) follows from (5.8). Using the matrix representation (5.6), we get and by induction, F a 00 0 = aη + Gη G F n a n 00 0 = (5.0) a n η + a n Gη + + G n η G n a n = (a n G n )(a G) η G n Hence, if we show that lim n a n G n = 0, (5.) then the desired conclusion (5.3) will follow. Using the matrix form (5.0), the conclusion (5.3) reads lim n a n F n = (5.) (a G) η 0 In proving (5.), we will make use of the Jordan-form representation for G. Jordan s theorem applied to G yields three operators D, V, N L (M ) with the following properties: () D is a diagonal matrix with the numbers spec (F) \ {a} down the diagonal; () V is invertible; (3) N is nilpotent: If d = dim (M ) then N d = 0; 8

29 (4) [N, D] = ND DN = 0; (5) G = V (D + N) V. Let x M, and let n d. Using () (5), we get a n G n x = V a n (D + N) n V x ( ) d n = V a n D n i N i V x. i i=0 But the matrix a n D n i is diagonal with entries { a n s n i s spec (F) \ {a} }, 0 i d. Using finally assumption (ii), we conclude that lim n ( ) n a n s n i = 0, (5.3) i and the proof of (5.) is completed. PROOF. [Lemma 5.3] Let the conditions be as stated in the Remark. From the arguments in the proof of Lemma 5., we see that the two vectors on the left-hand side in (5.4) may be decomposed as follows: and a n F n x = w x w + ( a n G n) (a G) η + a n G n P M x; (5.4) Hence, the difference is in M, and a n F n x w x ξ w x ξ = w x w + (a G) η. (5.5) = a n G ( n P M x (a G) η ) { Cn d s n } max s spec (F) \ {a} a which is the desired conclusion (5.4). References [] O. Bratteli, D. E. Evans, P. E. T. Jorgensen, Compactly Supported Wavelets and Representations of the Cuntz Relations, Appl. Comp. Harm. Anal., 00, (000) 3. 9

30 [] O. Bratteli, P. E. T. Jorgensen, Wavelets through a Looking Glass: The World of the Spectrum, Applied and Numerical Harmonic Analysis, Birkhäuser, Boston, 00. [3] O. Bratteli, P. E. T. Jorgensen, Iterated function systems and permutation representations of the Cuntz algebra, Mem. Amer. Math. Soc. 39 (999) no [4] R. Coifman, Y. Meyer, V. Wickerhauser, Size properties of wavelet-packets, in: M. Ruskai, G. Beylkin, R. Coifman, I. Daubechies, S. Mallat, Y. Meyer, L. Raphael (Eds.), Wavelets and Their Applications, Jones and Bartlett, Boston, 99, pp [5] J. Cuntz, Simple C -algebras generated by isometries, Comm. Math. Phys. 57 (977) [6] D. Dutkay, P. Jorgensen, Wavelets on Fractals, preprint, University of Iowa, 004, arxiv:math.ca/ , accepted: Rev. Mat. Iberoamericana. [7] K. J. Falconer, The Geometry of Fractal Sets, Cambridge, Tracts on Mathematics, 85. [8] J. E. Hutchinson, Fractals and self-similarity, Indiana Univ. Math. J. 30 (98) [9] P.E.T. Jorgensen, Minimality of the Data in Wavelet Filters, Advances in Mathematics 59, 00, pp [0] P.E.T. Jorgensen, Measures in wavelet decompositions, preprint, University of Iowa, 004, arxiv:math.ca/04004, submitted to Elsevier Science. [] M. V. Wickerhauser, Best-adapted wavelet packet bases, in: I. Daubechies (Ed.), Different Perspectives on Wavelets (San Antonio, TX, 993), Vol. 47 of Proc. Sympos. Appl. Math., American Mathematical Society, Providence, 993, pp [] M. V. Wickerhauser, Adapted Wavelet Analysis from Theory to Software, IEEE Press, New York, A.K. Peters, Wellesley, MA,