Computational Harmonic Analysis (Wavelet Tutorial) Part II
Understanding Many Particle Systems with Machine Learning
Tutorials
Matthew Hirn Michigan State University Department of Computational Mathematics, Science & Engineering Department of Mathematics Wavelet Transform Wavelets
2 Wavelet L (R)satisfies: • 2 – Zero average: =0 – Normalized: R =1 k k2 – Centered around t =0 – Localized in time and frequency – Can be either real or complex valued Wavelet Transform f(t) Wavelet dictionary obtained by scaling and • 2 translating : 1 1 t u 0 = u,s +, u,s(t)= t D { }u R,s R ps s 0 0.2 0.4 0.6 0.8 1 2 2 ✓ ◆ log2(s) −6 Wavelet transform: • −4 Wf(u, s)= f, u,s h + i = 1 f(t)s 1/2 (s 1(t u)) dt −2 Z 1 = f ˜s(u) 0 u ⇤ 0 0.2 0.4 0.6 0.8 1 Fig. 4.7. A Wavelet Tour of Signal Processing, 3rd ed. Real wavelet transform Wf(u, s) computed with a Mexican hat wavelet The where vertical axis represents log2 s. Black, grey and white points correspond respectively to positive, zero and negative wavelet coe cients. ˜ 1/2 1 s(t)=s (s t)
Note: • 0.25 ˜ ( )=p ˆ( ) s ! s s! 0.2 Thus, since: b 0.15 f\ ˜s(!)=fˆ(!) ˜ (!) 0.1 ⇤ s the wavelet transform Wfb(u, s) captures 0.05 the frequency information of f organized 0 by the frequency bands of ˜s. −2 0 2 Fig. 5.1. A Wavelet Tour of Signal Processing, 3rd ed. Scaled Fourier transforms ˆ(2j!) 2, for 1 j 5 and ! [ ⇡, ⇡]. | | 6 6 2 Real Wavelet Reconstruction Theorem (Calder´on,Grossman and Mor- • 2 let): Let L (R) be a real function such 2 that + ˆ 2 1 (!) C = | | d! < + Z0 ! 1 2 Then, for any f L (R): 2 + + 1 1 1 1/2 1 ds f(t)= Wf(u, s)s (s (t u)) du C 0 s2 Z Z 1 + + 2 1 1 1 2 ds f 2 = Wf(u, s) du . k k C 0 | | s2 Z Z 1 C < is called the wavelet admissibility • 1 condition.
C < + ˆ(0) = 0. This is almost • 1) su cient.
If additionally, ˆ C1, then C < + . • 2 1 Can insure this with su cient time decay: K (t) | | 1+ t 2+✏ | | Scaling Function
Numerically the wavelet transform is only • computed up to scales s