多贝西小波
维基百科,自由的百科全书
多貝西小波(Daubechies Wavelet)是以英格麗多貝西(Ingrid Daubechies)的名字命名的一種小波函式,當初Ingrid Daubechies發現了一種具有階層(hierarchy)性質的小波,便將此小波以她的名字命名。多貝西小波主要應用在離散型的小波轉換,是最常使用到的小波轉換,通常使用在數位信號分析、信號壓縮跟雜訊去除。
一般而言的離散小波轉換通常是以正交小波(orthogonal wavelet)為基底,而多貝西小波也是一種正交小波。由於它很容易經由快速小波轉換(fast wavelet transform (FWT))實現,所以常會放在數位信號處理的教科書中教學。對於有限長度的小波,應用於快速小波轉換(fast wavelet transform (FWT))時,會有兩個實數組成的數列:一是作為高通濾波器的係數,稱作小波濾波器(wavelet filter, 也稱為mother wavelet);二是低通濾波器的係數,稱作調整濾波器(scaling filter, 也稱為father wavelet)。我們則以濾波器的長度N來形容濾波器為DN,例如:N=2的多貝西小波寫作D2、N=4的多貝西小波寫作D4,以此類推(N為偶數)。實際上常用的多貝西小波為D2到D20。
目录 |
[编辑] 性質
- 分類方式
- 多貝西小波的分類是以消失動量(vanishing moment)的值A(亦為消失動量的個數)為依據(A稱為tap),調整函式(scaling function)及小波函式(wavelet function)的平滑度(smoothness)皆會隨著消失動量的值(tap)增加而增加:例如,當A=1時,多貝西小波即是哈爾小波(Haar wavelet),調整函式及小波函式都是不連續的;當A=2時,多貝西小波的調整函式及小波函式為不能平滑微分的連續函式;當A=3時,調整函式及小波函式已經是連續可微的函式了。以此類推,當A愈大時,多貝西小波的兩個函式平滑度會愈來愈高。以下為多貝西小波跟不同A的調整及小波函式圖:
scaling and wavelet functions | |||
amplitudes of the frequency spectrum of the above functions |
- 長度
- 多貝西小波的長度為消失動量(vanishing moment)值A的兩倍;所以當消失動量為A時,多貝西小波的小波濾波器(wavelet filter)及調整濾波器(scaling filter)長度皆為2A(N=2A)。一般而言,我們仍是以N來形容多貝西小波的長度:例如,當A=1時,有一個消失動量,多貝西小波寫成D2,長度為2(也是Haar小波);當A=2時,有兩個消失動量,多貝西小波寫成D4,長度為4;以此類推。但是,在matlab的使用上是以dbA描述多貝西小波,以下則為調整濾波器的係數及A的關係表:
Scaling Coefficient pk |
db1 (Haar) | db2 | db3 | db4 | db5 | db6 | db7 | db8 | db9 | db10 |
---|---|---|---|---|---|---|---|---|---|---|
p0 | 1 | 0.6830127 | 0.47046721 | 0.32580343 | 0.22641898 | 0.15774243 | 0.11009943 | 0.07695562 | 0.05385035 | 0.03771716 |
p1 | 1 | 1.1830127 | 1.14111692 | 1.01094572 | 0.85394354 | 0.69950381 | 0.56079128 | 0.44246725 | 0.34483430 | 0.26612218 |
p2 | 0.3169873 | 0.650365 | 0.8922014 | 1.02432694 | 1.06226376 | 1.03114849 | 0.95548615 | 0.85534906 | 0.74557507 | |
p3 | -0.1830127 | -0.19093442 | -0.03957503 | 0.19576696 | 0.44583132 | 0.66437248 | 0.82781653 | 0.92954571 | 0.97362811 | |
p4 | -0.12083221 | -0.26450717 | -0.34265671 | -0.31998660 | -0.20351382 | -0.02238574 | 0.18836955 | 0.39763774 | ||
p5 | 0.0498175 | 0.0436163 | -0.04560113 | -0.18351806 | -0.31683501 | -0.40165863 | -0.41475176 | -0.35333620 | ||
p6 | 0.0465036 | 0.10970265 | 0.13788809 | 0.1008467 | 6.68194092e-4 | -0.13695355 | -0.27710988 | |||
p7 | -0.01498699 | -0.00882680 | 0.03892321 | 0.11400345 | 0.18207636 | 0.21006834 | 0.18012745 | |||
p8 | -0.01779187 | -0.04466375 | -0.05378245 | -0.02456390 | 0.043452675 | 0.13160299 | ||||
p9 | 4.71742793e-3 | 7.83251152e-4 | -0.02343994 | -0.06235021 | -0.09564726 | -0.10096657 | ||||
p10 | 6.75606236e-3 | 0.01774979 | 0.01977216 | 3.54892813e-4 | -0.04165925 | |||||
p11 | -1.52353381e-3 | 6.07514995e-4 | 0.01236884 | 0.03162417 | 0.04696981 | |||||
p12 | -2.54790472e-3 | -6.88771926e-3 | -6.67962023e-3 | 5.10043697e-3 | ||||||
p13 | 5.00226853e-4 | -5.54004549e-4 | -6.05496058e-3 | -0.01517900 | ||||||
p14 | 9.55229711e-4 | 2.61296728e-3 | 1.97332536e-3 | |||||||
p15 | -1.66137261e-4 | 3.25814671e-4 | 2.81768659e-3 | |||||||
p16 | -3.56329759e-4 | -9.69947840e-4 | ||||||||
p17 | -5.5645514e-5 | -1.64709006e-4 | ||||||||
p18 | 1.32354367e-4 | |||||||||
p19 | -1.875841e-5 |
[编辑] 建立
多貝西小波具有調整函式(低通濾波)及小波函式(高通濾波)兩個函式。因此,我們需先建立調整函式及小波函式的係數:
- 首先,調整函數在多尺度分析(multi-resolution analysis)中的每一層皆可寫為下列方程式:
- ,
- 其中 為有限長度實數數列,稱作調整係數。同時,小波函數也可以調整函數的線性組合表示:
- ,
- 其中 亦為有限長度的實數數列,稱做小波係數。
- 因為上述方程式必須是齊性的(homogeneous),在建立上,這兩個函式會正規化(normalize)為和(sum)及平方和(sum of square)皆是2。
- 正交小波
- 正交性質在此指調整係數就必須和位移偶數間隔後的調整係數互相垂直(內積為0),即為下式:
- .
- 正交性質在此指調整係數就必須和位移偶數間隔後的調整係數互相垂直(內積為0),即為下式:
-
- 由於正交的特性,小波係數會滿足下列條件:
- bn = ( − 1)naN − 1 − n, .
- 由於正交的特性,小波係數會滿足下列條件:
- 消失動量及多項式估計
- 常用的多貝西小波為D2到D20,由於多貝西小波的消失動量為有限個,所以調整及小波係數可以表示為有限長度的多項式
- a[n] = an
- 常用的多貝西小波為D2到D20,由於多貝西小波的消失動量為有限個,所以調整及小波係數可以表示為有限長度的多項式
-
- 上式經過Z轉換(Z-transform)後會變成:
- ,
- 上式經過Z轉換(Z-transform)後會變成:
-
- 我們可以將上式轉換為正交離散小波轉換的一般表示式
- a(Z) = 21 − A(1 + Z)Ap(Z), order = A,此時, p(1) = 1、p有實係數及< math>deg(p) = A-1</math>。
- 我們可以將上式轉換為正交離散小波轉換的一般表示式
-
- 而正交的條件可寫成
- , 或是等同於 (#),
- X定義為可以產生對稱數列的勞倫茲多項式 滿足X( − Z) = 2 − X(Z)
- 而正交的條件可寫成
-
- 因此P(X)便成為對稱型勞倫滋多項式,即 P(X(Z)) = P(Z)P(Z − 1)。 因為 X(eiw) = 1 − cos(w) 及 p(eiw)p(e − iw) = | p(eiw) | 2, P則會是區段[0,2]中的非負實數。
-
- 方程式(#)如果除上X的truncated power series 則可求得對於每個A的最小解
- . (明顯的值會是在(0,2)間的正數)
- 方程式(#)如果除上X的truncated power series 則可求得對於每個A的最小解
-
- 而(#)的齊性方程式是一個對於x = 1的反對稱方程式,因此可得一般解為XA(X − 1)R((X − 1)2),此一般解有R個多項式實係數。
- 因此和為
- (sum)P(X) = PA(X) + XA(X − 1)R((X − 1)2)
- P的值在區間[0,2]中並有界線(界線為4A − r,)。為了將 r最大化的過程中會產生許多具有不等式條件的線性方程式。
- 為了解出P(X(Z)) = P(Z)P(Z − 1)的P,這裡使用Fejer-Riesz-algorithm這個方法(此為頻譜分解的方法)。多項式P(X)會因此分開成許多線性因子(linear factor), 此時N = A + 1 + 2deg(R)。每一個線性因子代表可以分解成兩個線性因子的一個勞倫茲多項式 ,任選其中一個線性因子都可設為p(Z)。所以p(Z)會有2N個可能的答案。為了極端相位的目的,挑選所有根都是在單位圓上或是在單位圓內複數根的P(Z)。
[编辑] 參照
[编辑] 參考資料
- A Fcourse in Wavelets with Fourier Analysis, A.Boggess, F.J. Narcowich, 2001