Daubechies wavelet

From TCS Wiki
Revision as of 00:58, 2 January 2015 by imported>Auntof6 (use standard simple headings and/or general cleanup using AWB)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Template:Complex Template:Orphan

File:Daubechies4-functions.png
Scaling and wavelet function of continuous D4 wavelet

Daubechies wavelets are a family of orthogonal wavelets named after Belgian physicist and mathematician Ingrid Daubechies. They are used in discrete wavelet transform.

Definition

Scale function coefficients (low pass filter in orthogonal filter banks) must satisfy following conditions ([math]\displaystyle{ N }[/math] is length of filter).

Normalization
[math]\displaystyle{ \sum_{n=0}^{N-1} h_0[n] = \sqrt{2} }[/math] or [math]\displaystyle{ \sum_{n=0}^{N-1} h_0[n] = 2 }[/math] (then coefficients must be divided by factor of [math]\displaystyle{ \sqrt{2} }[/math])

which implies

[math]\displaystyle{ \sum_{n=0}^{N-1} (h_0[n])^2 = 1 }[/math] or [math]\displaystyle{ \sum_{n=0}^{N-1} (h_0[n])^2 = 2 }[/math] (then coefficients must be divided by factor of [math]\displaystyle{ \sqrt{2} }[/math])
Orthogonality
[math]\displaystyle{ \sum_{n=0}^{N-1} h_0[n] h_0[n-2k] = 0 }[/math] for [math]\displaystyle{ k \not= 0 }[/math]
Vanishing moments
[math]\displaystyle{ \sum_{n=0}^{N-1} (-1)^n h_0[n] n^m = 0 }[/math] for [math]\displaystyle{ 0 \leq m \lt N/2 }[/math]

There is more than one solution (except case of [math]\displaystyle{ N=2 }[/math]). However, it is necessary to distinguish between low pass and high pass filter.

Wavelets are denoted like Dx, where x is either number of coefficients ([math]\displaystyle{ N }[/math]) or number of vanishing moments ([math]\displaystyle{ N/2 }[/math]). First case of notation (number of coefficients) is more recent and more frequented (e.g. D8 is wavelet with 8 coefficients).

Example

MATLAB code for enumeration of wavelet with 4 coefficients (denoted as D4). <source lang="matlab"> t = solve( 'h0*h0 + h1*h1 + h2*h2 + h3*h3 = 1',  % normalization 'h2*h0 + h3*h1 = 0',  % orthogonality '+(0^0)*h0 -(1^0)*h1 +(2^0)*h2 -(3^0)*h3 = 0', % zero '+(0^1)*h0 -(1^1)*h1 +(2^1)*h2 -(3^1)*h3 = 0'  % and first vanishing moments ); </source>

Solutions (low pass filters only):

h0 h1 h2 h3
-0.129409522551260 0.224143868042014 0.836516303737808 0.482962913144534
0.482962913144534 0.836516303737808 0.224143868042014 -0.129409522551260

Related pages