Fourier Bases on Fractals Keri Kornelson University of Oklahoma - - - PowerPoint PPT Presentation

fourier bases on fractals
SMART_READER_LITE
LIVE PREVIEW

Fourier Bases on Fractals Keri Kornelson University of Oklahoma - - - PowerPoint PPT Presentation

Fourier Bases on Fractals Keri Kornelson University of Oklahoma - Norman February Fourier Talks February 21, 2013 K. Kornelson (U. Oklahoma) Fractal Fourier Bases FFT 02/21/2013 1 / 21 Coauthors This is joint work with Palle Jorgensen


slide-1
SLIDE 1

Fourier Bases on Fractals

Keri Kornelson

University of Oklahoma - Norman

February Fourier Talks February 21, 2013

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 1 / 21

slide-2
SLIDE 2

Coauthors

This is joint work with Palle Jorgensen (University of Iowa) Karen Shuman (Grinnell College).

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 2 / 21

slide-3
SLIDE 3

Outline

1

Bernoulli convolution measures

2

Fourier bases

3

Families of ONBs

4

Operator-fractal

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 3 / 21

slide-4
SLIDE 4

Bernoulli convolution measures

Convolution measure

Let λ ∈ (0, 1). A fractal Cantor subset of R is the unique set Xλ satisfying the invariance relation: Xλ = λ(Xλ + 1) ∪ λ(Xλ − 1). (1)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 4 / 21

slide-5
SLIDE 5

Bernoulli convolution measures

Convolution measure

Let λ ∈ (0, 1). A fractal Cantor subset of R is the unique set Xλ satisfying the invariance relation: Xλ = λ(Xλ + 1) ∪ λ(Xλ − 1). (1) The Bernoulli convolution measure µλ is the unique probability measure satisfying:

  • f dµλ = 1

2

  • f(λx + λ) + f(λx − λ) dµλ(x).

(2)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 4 / 21

slide-6
SLIDE 6

Bernoulli convolution measures

Convolution measure

Let λ ∈ (0, 1). A fractal Cantor subset of R is the unique set Xλ satisfying the invariance relation: Xλ = λ(Xλ + 1) ∪ λ(Xλ − 1). (1) The Bernoulli convolution measure µλ is the unique probability measure satisfying:

  • f dµλ = 1

2

  • f(λx + λ) + f(λx − λ) dµλ(x).

(2) Historical note: The Bernoulli measures date back to work of Erdös and others in the 1930s and 40s. µλ is the distribution of the random variable

k ±λk where + and −

have equal probability.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 4 / 21

slide-7
SLIDE 7

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-8
SLIDE 8

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

If λ < 1

2, then µλ is supported on a fractal with Lebesgue measure zero.

Thus, the measures are singular when λ < 1

2.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-9
SLIDE 9

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

If λ < 1

2, then µλ is supported on a fractal with Lebesgue measure zero.

Thus, the measures are singular when λ < 1

2.

When λ ≥ 1

2, the set Xλ is an interval.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-10
SLIDE 10

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

If λ < 1

2, then µλ is supported on a fractal with Lebesgue measure zero.

Thus, the measures are singular when λ < 1

2.

When λ ≥ 1

2, the set Xλ is an interval.

[Erdös, 1939] If λ is the inverse of a Pisot number, e.g. the golden ratio φ, the measure µλ is singular with respect to Lebesgue measure even though the attractor is an interval.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-11
SLIDE 11

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

If λ < 1

2, then µλ is supported on a fractal with Lebesgue measure zero.

Thus, the measures are singular when λ < 1

2.

When λ ≥ 1

2, the set Xλ is an interval.

[Erdös, 1939] If λ is the inverse of a Pisot number, e.g. the golden ratio φ, the measure µλ is singular with respect to Lebesgue measure even though the attractor is an interval. [Hutchinson, 1981] Existence of (X, µ), from an IFS perspective.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-12
SLIDE 12

Bernoulli convolution measures

Some properties of Bernoulli measures

Given Bernoulli measure with scale factor λ: When λ = 1

2, the measure µ 1

2 is scaled Lebesgue measure on [−1, 1].

If λ < 1

2, then µλ is supported on a fractal with Lebesgue measure zero.

Thus, the measures are singular when λ < 1

2.

When λ ≥ 1

2, the set Xλ is an interval.

[Erdös, 1939] If λ is the inverse of a Pisot number, e.g. the golden ratio φ, the measure µλ is singular with respect to Lebesgue measure even though the attractor is an interval. [Hutchinson, 1981] Existence of (X, µ), from an IFS perspective. [Solomyak, 1995] For almost every λ ∈ ( 1

2, 1), µλ is absolutely

continuous.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 5 / 21

slide-13
SLIDE 13

Fourier bases

  • µλ as an infinite product

Consider the Hilbert space L2(µλ). Is it possible for L2(µλ) to have a Fourier basis, i.e. an orthonormal basis (ONB) of complex exponential functions?

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 6 / 21

slide-14
SLIDE 14

Fourier bases

  • µλ as an infinite product

Consider the Hilbert space L2(µλ). Is it possible for L2(µλ) to have a Fourier basis, i.e. an orthonormal basis (ONB) of complex exponential functions? Recall that µλ satisfies the invariance

  • f dµλ = 1

2

  • f(λx + λ) + f(λx − λ) dµλ(x).
  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 6 / 21

slide-15
SLIDE 15

Fourier bases

  • µλ as an infinite product

Consider the Hilbert space L2(µλ). Is it possible for L2(µλ) to have a Fourier basis, i.e. an orthonormal basis (ONB) of complex exponential functions? Recall that µλ satisfies the invariance

  • f dµλ = 1

2

  • f(λx + λ) + f(λx − λ) dµλ(x).

Then the Fourier transform of µλ is:

  • µλ(t)

=

  • e2πixt dµλ(x)

= 1 2

  • e2πi(λx+λ)t dµλ(x) + 1

2

  • e2πi(λx−λ)t dµλ(x)

= cos(2πλt) µλ(λt) = cos(2πλt) cos(2πλ2t) µλ(λ2t) =

  • k=1

cos(2πλkt)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 6 / 21

slide-16
SLIDE 16

Fourier bases

Orthogonality Condition

Denote by eγ the exponential function e2πiγ· in L2(µλ).

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 7 / 21

slide-17
SLIDE 17

Fourier bases

Orthogonality Condition

Denote by eγ the exponential function e2πiγ· in L2(µλ). eγ, e˜

γL2

=

  • eγ−˜

γ dµλ

=

  • µλ(γ − ˜

γ) =

  • k=1

cos

  • λ

k (γ − ˜ γ)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 7 / 21

slide-18
SLIDE 18

Fourier bases

Orthogonality Condition

Denote by eγ the exponential function e2πiγ· in L2(µλ). eγ, e˜

γL2

=

  • eγ−˜

γ dµλ

=

  • µλ(γ − ˜

γ) =

  • k=1

cos

  • λ

k (γ − ˜ γ)

  • Lemma

The two exponentials eγ, e˜

γ are orthogonal if and only if one of the factors in

the infinite product above is zero. This is equivalent to γ − ˜ γ ∈ 1 4λ−k(2m + 1) : k ∈ N, m ∈ Z

  • =: Zλ.
  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 7 / 21

slide-19
SLIDE 19

Fourier bases

Test for ONB

We use the zero set Zλ to check orthogonality. Parseval’s identity provides a test for an ONB. Let Γ ⊂ R be a set and let E(Γ) be the set {eγ : γ ∈ Γ}.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 8 / 21

slide-20
SLIDE 20

Fourier bases

Test for ONB

We use the zero set Zλ to check orthogonality. Parseval’s identity provides a test for an ONB. Let Γ ⊂ R be a set and let E(Γ) be the set {eγ : γ ∈ Γ}. If E(Γ) is an ONB, then for every value of t ∈ R, we have 1 = et2

µλ

=

  • γ∈Γ

|et, eγ|2 =

  • γ∈Γ

| µλ(t − γ)|2

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 8 / 21

slide-21
SLIDE 21

Fourier bases

Test for ONB

We use the zero set Zλ to check orthogonality. Parseval’s identity provides a test for an ONB. Let Γ ⊂ R be a set and let E(Γ) be the set {eγ : γ ∈ Γ}. If E(Γ) is an ONB, then for every value of t ∈ R, we have 1 = et2

µλ

=

  • γ∈Γ

|et, eγ|2 =

  • γ∈Γ

| µλ(t − γ)|2 It other words, cΓ(t) =

  • γ∈Γ

[ µλ(t − γ)]2 =

  • γ∈Γ

  • k=1

cos2

  • λ

k (t − γ)

  • ≡ 1.

The function cΓ is sometimes called a spectral function .

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 8 / 21

slide-22
SLIDE 22

Fourier bases

First results

Theorem (Jorgensen, Pedersen 1998) L2(µ 1

4 ) has an ONB of exponential functions.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 9 / 21

slide-23
SLIDE 23

Fourier bases

First results

Theorem (Jorgensen, Pedersen 1998) L2(µ 1

4 ) has an ONB of exponential functions.

Example E(Γ 1

4 ) is an ONB for L2(µ 1 4 ), where

Γ 1

4 =

  

p

  • j=0

aj4j : aj ∈ {0, 1}, p finite    = {0, 1, 4, 5, 16, 17, 20, 21, 64, . . .}.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 9 / 21

slide-24
SLIDE 24

Fourier bases

Another surprise: λ = 1

3 Theorem (Jorgensen, Pedersen 1998) Not only is there no Fourier basis when λ = 1

3, but orthogonal collections of

exponential functions can have at most 2 elements. There is also a more general version in the [JP1998] paper:

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 10 / 21

slide-25
SLIDE 25

Fourier bases

Another surprise: λ = 1

3 Theorem (Jorgensen, Pedersen 1998) Not only is there no Fourier basis when λ = 1

3, but orthogonal collections of

exponential functions can have at most 2 elements. There is also a more general version in the [JP1998] paper: Theorem (Jorgensen, Pedersen 1998) Given λ = 1

n, if n is even, there is an ONB of exponentials for L2(µ 1

2n ) but

when n is odd, there can be only finitely many elements in any orthogonal collection of exponentials.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 10 / 21

slide-26
SLIDE 26

Fourier bases

More recent progress

Theorem (Jorgensen, K, Shuman 2008; Hu, Lau 2008) For λ = a

b, if b is odd, then any orthonormal collection of exponentials in

L2(µλ) must be finite. If b is even, then there exists countable collections of

  • rthonormal exponentials in L2(µλ).
  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 11 / 21

slide-27
SLIDE 27

Fourier bases

More recent progress

Theorem (Jorgensen, K, Shuman 2008; Hu, Lau 2008) For λ = a

b, if b is odd, then any orthonormal collection of exponentials in

L2(µλ) must be finite. If b is even, then there exists countable collections of

  • rthonormal exponentials in L2(µλ).

Theorem (Dutkay, Han, Jorgensen 2009) If λ > 1

2, i.e. there is essential overlap, then L2(µλ) does not have an ONB (or

even a frame) of exponential functions.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 11 / 21

slide-28
SLIDE 28

Fourier bases

More recent progress

Theorem (Jorgensen, K, Shuman 2008; Hu, Lau 2008) For λ = a

b, if b is odd, then any orthonormal collection of exponentials in

L2(µλ) must be finite. If b is even, then there exists countable collections of

  • rthonormal exponentials in L2(µλ).

Theorem (Dutkay, Han, Jorgensen 2009) If λ > 1

2, i.e. there is essential overlap, then L2(µλ) does not have an ONB (or

even a frame) of exponential functions. Theorem (Xinrong Dai 2012) The only spectral Bernoulli measures are for λ =

1 2n.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 11 / 21

slide-29
SLIDE 29

Fourier bases

“Canonical” ONBs

[Jorgensen, Pedersen 1998]

Definition Let λ =

1 2n and consider the set from Jorgensen & Pedersen

Γ 1

2n =

  

p

  • j=0

aj(2n)j : aj ∈

  • 0, n

2

  • , p finite

   . We call Γ 1

2n the canonical spectrum and E(Γ 1 2n ) the canonical ONB for L2(µ 1 2n ).

Note: We will justify the nomenclature by describing alternate bases for the same spaces.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 12 / 21

slide-30
SLIDE 30

Families of ONBs

Families of ONBs

Let λ =

1

  • 2n. We can construct alternate orthogonal families of exponentials

from the canonical ONBs E(Γ 1

2n ). We then determine whether these alternate

sets are ONBs.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 13 / 21

slide-31
SLIDE 31

Families of ONBs

Families of ONBs

Let λ =

1

  • 2n. We can construct alternate orthogonal families of exponentials

from the canonical ONBs E(Γ 1

2n ). We then determine whether these alternate

sets are ONBs. Recall the zero set, which gives the test for orthogonality of exponentials: Zλ = 1 4λ−k(2m + 1) : k ∈ N, m ∈ Z

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 13 / 21

slide-32
SLIDE 32

Families of ONBs

Families of ONBs

Let λ =

1

  • 2n. We can construct alternate orthogonal families of exponentials

from the canonical ONBs E(Γ 1

2n ). We then determine whether these alternate

sets are ONBs. Recall the zero set, which gives the test for orthogonality of exponentials: Zλ = 1 4λ−k(2m + 1) : k ∈ N, m ∈ Z

  • Observe that if γ − ˜

γ ∈ Zλ then pγ − p˜ γ ∈ Zλ as well, for any odd integer p. So scaling a canonical spectrum by p yields at least an orthogonal set, and sometimes an ONB.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 13 / 21

slide-33
SLIDE 33

Families of ONBs

Families of ONBs

Let λ =

1

  • 2n. We can construct alternate orthogonal families of exponentials

from the canonical ONBs E(Γ 1

2n ). We then determine whether these alternate

sets are ONBs. Recall the zero set, which gives the test for orthogonality of exponentials: Zλ = 1 4λ−k(2m + 1) : k ∈ N, m ∈ Z

  • Observe that if γ − ˜

γ ∈ Zλ then pγ − p˜ γ ∈ Zλ as well, for any odd integer p. So scaling a canonical spectrum by p yields at least an orthogonal set, and sometimes an ONB. Note: Not every p yields an ONB, e.g. p = 2n − 1 for λ =

1

  • 2n. The set

E((2n − 1)Γ 1

2n ) is not maximal, hence is not an ONB.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 13 / 21

slide-34
SLIDE 34

Families of ONBs

Families of ONBs

Theorem (Jorgensen, K, Shuman 2010) The set E(3Γ 1

8 ) is an (alternate) orthonormal basis for L2(µ 1 8 ).

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 14 / 21

slide-35
SLIDE 35

Families of ONBs

Families of ONBs

Theorem (Jorgensen, K, Shuman 2010) The set E(3Γ 1

8 ) is an (alternate) orthonormal basis for L2(µ 1 8 ).

Theorem (Jorgensen, K, Shuman 2010) If p < 2(2n−1)

π

, then E(pΓ 1

2n ) is an ONB for L2(µ 1 2n ).

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 14 / 21

slide-36
SLIDE 36

Families of ONBs

Families of ONBs

Theorem (Jorgensen, K, Shuman 2010) The set E(3Γ 1

8 ) is an (alternate) orthonormal basis for L2(µ 1 8 ).

Theorem (Jorgensen, K, Shuman 2010) If p < 2(2n−1)

π

, then E(pΓ 1

2n ) is an ONB for L2(µ 1 2n ).

Laba/Wang and Dutkay/Jorgensen have described many other values of p for which pΓ 1

2n is a spectrum, particularly in the 1

4 case.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 14 / 21

slide-37
SLIDE 37

Families of ONBs

Examples of pΓ 1

2n ONBs

n λ p Canonical Γλ 3

1 6

1, 3 {0, 3

2, 9, 21 2 , . . .}

4

1 8

1, 3 {0, 2, 16, 18, . . .} 5

1 10

1, 3, 5 {0, 5

2, 25, 55 2 , . . .}

6

1 12

1, 3, 5, 7 {0, 3, 36, 39, . . .}

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 15 / 21

slide-38
SLIDE 38

Operator-fractal

Operators on L2

Dutkay, Jorgensen, 2009: When λ = 1

4, both Γ 1

4 and 5Γ 1 4 are spectra for L2(µ 1 4 ).

Γ 1

4 = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .}

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 16 / 21

slide-39
SLIDE 39

Operator-fractal

Operators on L2

Dutkay, Jorgensen, 2009: When λ = 1

4, both Γ 1

4 and 5Γ 1 4 are spectra for L2(µ 1 4 ).

Γ 1

4 = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .}

For each γ ∈ Γ 1

4 , define

S0 : eγ → e4γ S1 : eγ → e4γ+1 U : eγ → e5γ

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 16 / 21

slide-40
SLIDE 40

Operator-fractal

Operators on L2

Dutkay, Jorgensen, 2009: When λ = 1

4, both Γ 1

4 and 5Γ 1 4 are spectra for L2(µ 1 4 ).

Γ 1

4 = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .}

For each γ ∈ Γ 1

4 , define

S0 : eγ → e4γ S1 : eγ → e4γ+1 U : eγ → e5γ S0 and S1 map between ONB elements, so are both isometries.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 16 / 21

slide-41
SLIDE 41

Operator-fractal

Operators on L2

Dutkay, Jorgensen, 2009: When λ = 1

4, both Γ 1

4 and 5Γ 1 4 are spectra for L2(µ 1 4 ).

Γ 1

4 = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .}

For each γ ∈ Γ 1

4 , define

S0 : eγ → e4γ S1 : eγ → e4γ+1 U : eγ → e5γ S0 and S1 map between ONB elements, so are both isometries. U maps one ONB to another, so U is a unitary operator.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 16 / 21

slide-42
SLIDE 42

Operator-fractal

Structure of Γ

Example: Let λ = 1

4 and denote H = L2(µ 1

4 ).

Γ = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .} Γ = 4Γ ∪ (1 + 4Γ) = 42Γ ∪ 4(1 + 4Γ) ∪ (1 + 4Γ) . . . = 4kΓ ∪ 4k−1(1 + 4Γ) · · · 4(1 + 4Γ) ∪ (1 + 4Γ) S0(H) is the span of the exponentials E(4Γ)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 17 / 21

slide-43
SLIDE 43

Operator-fractal

Structure of Γ

Example: Let λ = 1

4 and denote H = L2(µ 1

4 ).

Γ = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .} Γ = 4Γ ∪ (1 + 4Γ) = 42Γ ∪ 4(1 + 4Γ) ∪ (1 + 4Γ) . . . = 4kΓ ∪ 4k−1(1 + 4Γ) · · · 4(1 + 4Γ) ∪ (1 + 4Γ) S0(H) is the span of the exponentials E(4Γ) Sk

0(H) is the span of E(4kΓ)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 17 / 21

slide-44
SLIDE 44

Operator-fractal

Structure of Γ

Example: Let λ = 1

4 and denote H = L2(µ 1

4 ).

Γ = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .} Γ = 4Γ ∪ (1 + 4Γ) = 42Γ ∪ 4(1 + 4Γ) ∪ (1 + 4Γ) . . . = 4kΓ ∪ 4k−1(1 + 4Γ) · · · 4(1 + 4Γ) ∪ (1 + 4Γ) S0(H) is the span of the exponentials E(4Γ) Sk

0(H) is the span of E(4kΓ)

H ⊃ S0(H) ⊃ S2

0(H) ⊃ · · ·

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 17 / 21

slide-45
SLIDE 45

Operator-fractal

Structure of Γ

Example: Let λ = 1

4 and denote H = L2(µ 1

4 ).

Γ = {0, 1, 4, 5, 16, 17, 20, 21, 64, 65, . . .} Γ = 4Γ ∪ (1 + 4Γ) = 42Γ ∪ 4(1 + 4Γ) ∪ (1 + 4Γ) . . . = 4kΓ ∪ 4k−1(1 + 4Γ) · · · 4(1 + 4Γ) ∪ (1 + 4Γ) S0(H) is the span of the exponentials E(4Γ) Sk

0(H) is the span of E(4kΓ)

H ⊃ S0(H) ⊃ S2

0(H) ⊃ · · ·

If we define Wk = Sk

0(H) ⊖ Sk+1

(H), then H = sp(e0) ⊕

  • k=0

Wk

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 17 / 21

slide-46
SLIDE 46

Operator-fractal

The operator U

Recall U : eγ → e5γ. How does that scaling by (×5) in the ONB frequencies interact with the inherent scaling (×4) of the fractal measure space? Theorem (Jorgensen,K,Shuman 2012)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 18 / 21

slide-47
SLIDE 47

Operator-fractal

The operator U

Recall U : eγ → e5γ. How does that scaling by (×5) in the ONB frequencies interact with the inherent scaling (×4) of the fractal measure space? Theorem (Jorgensen,K,Shuman 2012) Each subspace Wk is invariant under U.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 18 / 21

slide-48
SLIDE 48

Operator-fractal

The operator U

Recall U : eγ → e5γ. How does that scaling by (×5) in the ONB frequencies interact with the inherent scaling (×4) of the fractal measure space? Theorem (Jorgensen,K,Shuman 2012) Each subspace Wk is invariant under U. With respect to the Wk ordering of Γ, the matrix of U has block diagonal form...and the infinite blocks are all the same!

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 18 / 21

slide-49
SLIDE 49

Operator-fractal

The operator U

Recall U : eγ → e5γ. How does that scaling by (×5) in the ONB frequencies interact with the inherent scaling (×4) of the fractal measure space? Theorem (Jorgensen,K,Shuman 2012) Each subspace Wk is invariant under U. With respect to the Wk ordering of Γ, the matrix of U has block diagonal form...and the infinite blocks are all the same! Even more, in the (×4, ×5) case U actually has a self-similar structure: U = (e0 ⊗ e0) ⊕

  • k=1

Me1U. We call U an operator-fractal.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 18 / 21

slide-50
SLIDE 50

Operator-fractal

Matrix of U

Γ0 Γ1 Γ2 Γ3 · · · 1 · · · Γ0 Me1U · · · Γ1 Me1U · · · Γ2 Me1U · · · Γ3 Me1U · · · . . . . . . . . . . . . . . . ...

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 19 / 21

slide-51
SLIDE 51

Operator-fractal

Spectral properties of U

Theorem (Jorgensen, K, Shuman 2012) U has the following properties.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 20 / 21

slide-52
SLIDE 52

Operator-fractal

Spectral properties of U

Theorem (Jorgensen, K, Shuman 2012) U has the following properties. 1 is the only eigenvalue of U. (Ue0 = e0)

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 20 / 21

slide-53
SLIDE 53

Operator-fractal

Spectral properties of U

Theorem (Jorgensen, K, Shuman 2012) U has the following properties. 1 is the only eigenvalue of U. (Ue0 = e0) U is not spatially implemented; i.e. is not of the form Uf = f ◦ τ for τ a point transformation on [0, 1].

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 20 / 21

slide-54
SLIDE 54

Operator-fractal

Spectral properties of U

Theorem (Jorgensen, K, Shuman 2012) U has the following properties. 1 is the only eigenvalue of U. (Ue0 = e0) U is not spatially implemented; i.e. is not of the form Uf = f ◦ τ for τ a point transformation on [0, 1]. U only fixes the constant functions; if Uv = v then v = ce0 for some c ∈ C. In other words, U is an ergodic operator.

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 20 / 21

slide-55
SLIDE 55

Thank You

  • K. Kornelson (U. Oklahoma)

Fractal Fourier Bases FFT 02/21/2013 21 / 21