在數學中,連續小波變換 (Continuous Wavelet Transform, CWT)是一種時頻分析工具,通過讓小波函數(Wavelet)的平移參數和尺度參數連續變化,提供信號的過完備表示。
對具有頻率間斷點的信號使用五階消失矩的symlet 進行連續小波變換
函數
x
(
t
)
{\displaystyle x(t)}
在尺度
a
∈
R
+
∗
{\displaystyle a\in \mathbb {R^{+*}} }
和位置
b
∈
R
{\displaystyle b\in \mathbb {R} }
的連續小波變換表示為積分:
X
w
(
a
,
b
)
=
1
|
a
|
1
/
2
∫
−
∞
∞
x
(
t
)
ψ
¯
(
t
−
b
a
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{|a|^{1/2}}}\int _{-\infty }^{\infty }x(t){\overline {\psi }}\left({\frac {t-b}{a}}\right)\,dt}
其中
ψ
(
t
)
{\displaystyle \psi (t)}
稱為母小波,是在時間和頻率域中的連續函數,
⋅
¯
{\displaystyle {\bar {\cdot }}}
表示復共軛 。母小波的主要目的是為生成子小波(即母小波的尺度伸縮和平移)提供原函數。逆連續小波變換(Inverse continuous wavelet transform)可以用於恢復信號
x
(
t
)
{\displaystyle x(t)}
:
x
(
t
)
=
C
ψ
−
1
∫
0
∞
∫
−
∞
∞
X
w
(
a
,
b
)
1
|
a
|
1
/
2
ψ
~
(
t
−
b
a
)
d
b
d
a
a
2
{\displaystyle x(t)=C_{\psi }^{-1}\int _{0}^{\infty }\int _{-\infty }^{\infty }X_{w}(a,b){\frac {1}{|a|^{1/2}}}{\tilde {\psi }}\left({\frac {t-b}{a}}\right)\,db\ {\frac {da}{a^{2}}}}
ψ
~
(
t
)
{\displaystyle {\tilde {\psi }}(t)}
是
ψ
(
t
)
{\displaystyle \psi (t)}
的對偶函數 ,且:
C
ψ
=
∫
−
∞
∞
ψ
^
¯
(
ω
)
ψ
~
^
(
ω
)
|
ω
|
d
ω
{\displaystyle C_{\psi }=\int _{-\infty }^{\infty }{\frac {{\overline {\hat {\psi }}}(\omega ){\hat {\tilde {\psi }}}(\omega )}{|\omega |}}\,d\omega }
為容許性常數,式中
⋅
^
{\displaystyle {\hat {\cdot }}}
表示傅立葉變換算子。當母小波函數自對偶時(即
ψ
~
(
t
)
=
ψ
(
t
)
{\displaystyle {\tilde {\psi }}(t)=\psi (t)}
),容許性常數為:
C
ψ
=
∫
−
∞
+
∞
|
ψ
^
(
ω
)
|
2
|
ω
|
d
ω
{\displaystyle C_{\psi }=\int _{-\infty }^{+\infty }{\frac {\left|{\hat {\psi }}(\omega )\right|^{2}}{\left|\omega \right|}}\,d\omega }
一般來說,該常數為小波的容許性常數(Admissible constant)。滿足:
0
<
C
ψ
<
∞
{\displaystyle 0<C_{\psi }<\infty }
的小波稱為容許小波(Admissible wavelet)。容許小波滿足
ψ
^
(
0
)
=
0
{\displaystyle {\hat {\psi }}(0)=0}
,即積分為0。恢復信號
x
(
t
)
{\displaystyle x(t)}
的另一種逆連續小波變換為:
x
(
t
)
=
1
2
π
ψ
^
¯
(
1
)
∫
0
∞
∫
−
∞
∞
1
a
2
X
w
(
a
,
b
)
exp
(
i
t
−
b
a
)
d
b
d
a
{\displaystyle x(t)={\frac {1}{2\pi {\overline {\hat {\psi }}}(1)}}\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{a^{2}}}X_{w}(a,b)\exp \left(i{\frac {t-b}{a}}\right)\,db\ da}
該式說明,小波可以定義為:
ψ
(
t
)
=
w
(
t
)
exp
(
i
t
)
{\displaystyle \psi (t)=w(t)\exp(it)}
其中
w
(
t
)
{\displaystyle w(t)}
是一個窗。這種小波稱為解析小波(Analyzing wavelet),因為它被用於時頻分析(Time-frequency analysis)。解析小波不一定是容許小波。
小波轉換 ( Wavelet Transform) 可依照輸入與輸出為連續或是離散(discrete)分成三種類型,
第一種,輸入為連續,輸出為連續,則稱之為連續小波轉換(Continuous Wavelet Transform)
第二種,輸入為連續,輸出為離散,則稱之為連續離散系數小波轉換(Continuous wavelet transform with discrete coefficients)
第三種,輸入為離散,輸出為離散,則稱之為離散小波轉換(Discrete Wavelet Transform)
並沒有第四種,輸入為離散輸出為連續的小波轉換,在應用中並不會將簡單的訊號轉換成更複雜的訊號
傅立葉轉換 (Fourier Transform)與小波轉換 比較共有四種類型
第一種,輸入為連續,輸出為連續,傅立葉轉換(Fourier Transform)
第二種,輸入為連續,輸出為離散,傅立葉級數(Fourier Series)
第三種,輸入為離散,輸出為離散,離散傅立葉轉換(Discrete Fourier Transform)
第四種,輸入為離散,輸出為連續,離散(時間)傅立葉轉換(Discrete-time Fourier Transform)
連續小波轉換 (Continuous Wavelet Transform)是一種用來分解一個連續時間函數,使它變成數個小波 (wavelet)。跟傅立葉變換 (Fourier Transform)不一樣的是,連續小波轉換可以建構一個具有良好時域 和頻域 局部化的時頻訊號。以數學來說,一個有連續時間性質且可積分的函數
x
(
t
)
{\displaystyle x(t)}
可以用下面的積分來表示
X
w
(
a
,
b
)
=
1
|
(
b
)
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {|(b)|}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt}
ψ
(
t
)
{\displaystyle \psi (t)}
為小波母函數(Mother Wavelet),一個在時間領域和頻率領域都有連續性質的函數,
a
{\displaystyle a}
為平移位置而
b
{\displaystyle b}
為縮放因子。
a
{\displaystyle a}
的區間在
(
−
∞
,
+
∞
)
{\displaystyle \left(-\infty ,+\infty \right)}
,
b
{\displaystyle b}
的區間
(
0
,
+
∞
)
{\displaystyle \left(0,+\infty \right)}
以時頻分析的角度分析,當
b
{\displaystyle b}
值越大,頻率的變化越小
以時頻分析的角度分析,當
b
{\displaystyle b}
值越小,頻率的變化越大
小波母函數的用途在於提供一個可以產生子波(Daughter Wavelet)的根源函數,而子波是小波母函數平移過或縮放過(或兩者都有)的版本。如果要將已知且存在的訊號
x
(
t
)
{\displaystyle x(t)}
恢復原來的形式,我們可以用反轉連續小波轉換(Inverse Continuous Wavelet Transform)
x
(
t
)
=
∫
0
∞
∫
−
∞
∞
1
b
2
X
w
(
a
,
b
)
1
|
(
b
)
|
ψ
1
(
t
−
a
b
)
d
a
d
b
{\displaystyle x(t)=\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{b^{2}}}X_{w}(a,b){\frac {1}{\sqrt {|(b)|}}}\psi _{1}({\frac {t-a}{b}})\,dadb}
ψ
1
(
t
)
{\displaystyle \psi _{1}(t)}
為
ψ
(
t
)
{\displaystyle \psi (t)}
的對偶函數(Dual Function)。而這個對偶函數必須滿足
∫
0
∞
∫
−
∞
∞
1
|
b
3
|
ψ
(
t
1
−
a
b
)
ψ
1
(
t
−
a
b
)
d
a
d
b
=
δ
(
t
−
t
1
)
{\displaystyle \int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{|b^{3}|}}\psi ({\frac {t_{1}-a}{b}})\psi _{1}({\frac {t-a}{b}})\,dadb={\boldsymbol {\delta }}(t-t_{1})}
有時
ψ
1
(
t
)
=
C
−
1
ψ
(
t
)
{\displaystyle \psi _{1}(t)=C^{-1}\psi (t)}
。
在實務上,相較於連續小波轉換,通常較多會使用離散小波變換 、連續離散系數小波轉換
舉例,2個例子來說明小波母函數(Mother Wavelet):
縮放因子可以壓縮或拉長一個訊號。當縮放因子的值相對低時,訊號會比較緊縮,也就是會造成一個更細致的圖像。但是低縮放因子的缺點是它的效果無法覆蓋一個訊號的持續期間。另一方面,當縮放因子的值相對高時,訊號會比較被拉長造成一個比較粗糙的圖像,但是它的效果會持續整個訊號的期間。
1. 輸入與輸出的性質關係
基本上,連續小波轉換是輸入資料序列和一組由小波母函數所產生函數的卷積 。這個卷積可以用快速傅立葉變換來計算。除非小波母函數為虛數函數,在正常的情況下,輸出信息
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
會是一個實數函數。在小波母函數是虛數函數的情況下,連續小波轉換會造成一個虛數函數。連續小波轉換的功率譜可以以
‖
X
w
(
a
,
b
)
|
2
{\displaystyle \|X_{w}(a,b)|^{2}}
的數學型式來表示。通常在設計小波母函數時,為了應用上的目的會將小波母函數設計為實數函數。
2. 時間軸上的位移
輸入函數
x
(
t
)
{\displaystyle x(t)}
與輸出函數
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
之間具有相對的位移關係:
若
x
(
t
)
{\displaystyle x(t)}
經過小波轉換後的輸出函數為
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
則
x
(
t
−
τ
)
{\displaystyle x(t-\tau )}
經過小波轉換後的輸出函數為
X
w
(
a
−
τ
,
b
)
{\displaystyle X_{w}(a-\tau ,b)}
3. 時間軸上的縮放
當進行小波轉換的函數在時間軸上拉長或壓縮時,輸出函數也會有相對應的變化:
若
x
(
t
)
{\displaystyle x(t)}
經過小波轉換後的輸出函數為
X
w
(
a
,
b
)
{\displaystyle X_{w}(a,b)}
則
x
(
t
/
σ
)
{\displaystyle x(t/\sigma )}
經過小波轉換後的輸出函數為
σ
X
w
(
a
/
σ
,
b
/
σ
)
{\displaystyle {\sqrt {\sigma }}X_{w}(a/\sigma ,b/\sigma )}
4. 帕瑟伐定理(Parseval's Theory)
在
φ
(
t
)
=
C
−
1
ψ
(
t
)
{\displaystyle \varphi (t)=C^{-}1\psi (t)}
的條件下,滿足
∫
|
x
(
t
)
|
2
d
t
=
1
C
∫
0
∞
∫
−
∞
∞
1
b
2
|
X
w
(
a
,
b
)
|
2
d
a
d
b
{\displaystyle \int {|x(t)|}^{2}\,dt={\frac {1}{C}}\,\int _{0}^{\infty }\int _{-\infty }^{\infty }{\frac {1}{b^{2}}}{|X_{w}(a,b)|}^{2}\,da\,db}
加伯變換 在處理訊號時不管是高頻或是低頻,尺度皆是相同的
而小波轉換則會根據不同的頻率改變其本身的的尺度
小波變換的解像度在a-axis中,不因a值的改變而改變,但延著不同的b值改變以得到較好的結果。
利用上述尺度函數的定義,我們可以定義修正型的連續小波轉換。將原本的連續小波轉換定義為
X
w
(
a
,
b
)
=
1
b
′
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
′
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {b^{\prime }}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b^{\prime }}})\,dt}
,
在
a
{\displaystyle a}
屬於實數,且
0
<
b
′
<
b
0
{\displaystyle 0<b^{\prime }<b_{0}}
的情況下,則可定義
L
X
w
(
a
,
b
0
)
=
1
b
0
∫
−
∞
∞
x
(
t
)
φ
(
t
−
a
b
0
)
d
t
{\displaystyle {LX}_{w}(a,b_{0})={\frac {1}{\sqrt {b_{0}}}}\int _{-\infty }^{\infty }x(t)\varphi ({\frac {t-a}{b_{0}}})\,dt}
,
藉由新定義的函數
L
X
w
(
a
,
b
0
)
{\displaystyle {LX}_{w}(a,b_{0})}
,我們可以將反轉連續小波轉換表示為
x
(
t
)
=
1
C
ψ
[
∫
0
b
0
∫
−
∞
∞
1
b
′
5
/
2
X
w
(
a
,
b
′
)
ψ
(
t
−
a
b
′
)
d
a
d
b
+
∫
−
∞
∞
1
b
0
3
/
2
L
X
w
(
a
,
b
0
)
φ
(
t
−
a
b
0
)
d
a
]
{\displaystyle x(t)={\frac {1}{C_{\psi }}}\left[\int _{0}^{b_{0}}\int _{-\infty }^{\infty }{\frac {1}{{b^{\prime }}^{5/2}}}X_{w}(a,b^{\prime })\psi ({\frac {t-a}{b^{\prime }}})\,da\,db+\int _{-\infty }^{\infty }{\frac {1}{{b_{0}}^{3/2}}}{LX}_{w}(a,b_{0})\varphi ({\frac {t-a}{b_{0}}})\,da\right]}
此建構
x
(
t
)
{\displaystyle x(t)}
的方法可視為簡化版的反轉連續小波轉換。其中,
C
ψ
=
∫
0
∞
|
Ψ
(
f
)
|
2
|
f
|
d
f
{\displaystyle C_{\psi }=\int _{0}^{\infty }{\frac {|{\Psi (f)|}^{2}}{|f|}}\,df}
此式中
Ψ
(
f
)
{\displaystyle \Psi (f)}
為
ψ
(
f
)
{\displaystyle \psi (f)}
的傅立葉轉換。
通常在設計母小波函數時,會要求
C
ψ
<
∞
{\displaystyle C_{\psi }\ <\ \infty }
,此性質又稱為「可採納性(Admissibility Criterion)」。
在許多文獻資料中,常用小波量值圖(Scalogram)來表示連續小波轉換後的結果。其定義如下
S
c
x
(
a
,
b
)
=
|
X
w
(
a
,
b
)
|
2
=
1
|
b
|
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
|
2
{\displaystyle Sc_{x}(a,b)={\left|X_{w}(a,b)\right|}^{2}={\frac {1}{|b|}}{\left|\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt\right|}^{2}}
此處的定義即為連續小波轉換結果的絕對值平方,用以視覺化連續小波轉換的結果。
小波量值圖之於小波轉換函數的意義和頻譜圖 之於短時距時頻平分析的意義相似。
在實際應用上通常以三個軸來顯示,分別代表時間、頻率與小波量值圖的振幅。若是在二維圖片則是利用顏色深淺來表示小波量值圖的強度。
小波轉換最熱門的一個應用為圖像壓縮。用小波轉換式的編碼在圖像壓縮可以提供顯著的圖像品質改善且給予更高的壓縮比率。因為小波轉換可以分解一個複雜的訊息或圖案成基本型式,它在音樂檔案和圖型辨識上被廣泛的使用。此外,我們還可以在以下的科學研究領域見到小波轉換的應用: 邊緣檢測,解偏微分方程,腦電瞬態信號檢測,濾波器設計,心電圖分析,衣料分析和商業資訊分析。
在離散變數連續小波轉換(Continuous Wavelet Transform with Discrete Coefficients)中,原本
X
w
(
a
,
b
)
=
1
|
(
b
)
|
∫
−
∞
∞
x
(
t
)
ψ
(
t
−
a
b
)
d
t
{\displaystyle X_{w}(a,b)={\frac {1}{\sqrt {|(b)|}}}\int _{-\infty }^{\infty }x(t)\psi ({\frac {t-a}{b}})\,dt}
中的
a
,
b
{\displaystyle a,\,b}
具有一定的關係,不能隨意選取。若令
a
=
n
2
−
m
,
b
=
2
−
m
{\displaystyle a=n2^{-m},\ b=2^{-m}}
則離散變數連續小波轉換則重新表示為
X
w
(
a
,
b
)
=
2
m
/
2
∫
−
∞
∞
x
(
t
)
ψ
(
2
m
−
n
)
d
t
{\displaystyle X_{w}(a,b)=2^{m/2}\int _{-\infty }^{\infty }x(t)\psi (2^{m}-n)\,dt}
其中
n
∈
Z
,
n
∈
(
−
∞
,
∞
)
{\displaystyle n\in Z,\ n\in (-\infty ,\infty )}
且
n
∈
Z
,
n
∈
(
−
∞
,
∞
)
{\displaystyle n\in Z,\ n\in (-\infty ,\infty )}
選擇離散變數連續小波轉換的主要目的,在於簡化連續小波轉換在實作上的複雜性,並且利用快速演算法增加應用價值。實際上離散變數連續小波轉換算是連續小波轉換的一種特例,部分文獻將其當作離散小波轉換討論。
此處的選擇
a
{\displaystyle a}
與
b
{\displaystyle b}
之間的限制,使我們可以利用離散卷積的方式,由
X
w
(
n
,
m
−
1
)
{\displaystyle X_{w}(n,m-1)}
計算
X
w
(
n
,
m
)
{\displaystyle X_{w}(n,m)}
。
優點:
快速演算法
正交性質(Orthogonal )
非均勻頻率分析(Non-uniform frequency analysis)
缺點:
小波轉換的應用有以下兩項特點:
信號的頻率分佈,會隨着不同的時間 (或地點 )有較大變化
多尺度的分析扮演重要的角色
大採樣間隔
⟶
{\displaystyle \longrightarrow }
忽略細節信息
小採樣間隔
⟶
{\displaystyle \longrightarrow }
需要大量數據
應用:
影像壓縮,例如JPEG /JPEG2000
邊緣角落偵測
特徵辨識
強調前景壓縮背景
濾波器設計
聲音訊號
指紋辨識
金融
氣象分析