查看“︁穩定小波轉換”︁的源代码

{{multiple issues|
{{copyedit|time=2013-01-15T04:43:14+00:00}}
{{orphan|time=2013-01-15T04:43:14+00:00}}
{{unreferenced|time=2013-01-15T04:43:14+00:00}}
}}
'''穩定小波轉換'''（{{lang-en|Stationary wavelet transform}}，縮寫：'''SWT'''）是[[小波分析]]的一種轉換，為[[離散小波變換]]（DWT）的變形。

穩定小波轉換可以彌補離散小波轉換因為[[降采样]]而失去的[[平移不變性]]。穩定小波轉換不同於離散小波轉換的部分，主要在於經過每一階的[[高通濾波器]]和[[低通濾波器]]之後，是將[[濾波器]][[升採樣]]，取代離散小波轉換在經過濾波器之後的縮減取樣。

穩定小波轉換是做數據和信號的分析一種很好的工具，儘管它的運算量會因為沒有縮減取樣而较離散小波轉換多一些，但其具有平移不变性，且只需將離散小波轉換在濾波器的設計上做些微的修改即可實現。

== 實現方式 ==
下圖是穩定小波轉換的數位實現模型
[[File:SWT.PNG|有框|居中|三層的穩定小波轉換濾波器的濾波器組]]

每一組高通濾波器和低通濾波器皆為提升取樣後的前一組高通濾波器及低通濾波器，可以下圖表示：


以數學形式來呈現穩定小波轉換濾波器提升取樣的設計概念:<br />

令<math>Z</math>為一將<math>0</math>加入<math>x</math>序列的運算，<br />

<math>(Z_x)_k</math><math>=x_{k/2},</math>    for <math>k=2,4,6,8...</math><br />

<math>(Z_x)_k</math><math>=0,</math>    for <math>k=1,3,5,7...</math><br />
<br />
我們可以將第<math>j</math>階的高通濾波器<math>h_j[n]</math>表示成:<br />
<math>h_j[n]=Zh_{j-1}[n]=Z^2h_{j-2}[n]=...=Z^{j-1}h_{1}[n]</math><br />
同樣的，我們可以將第<math>j</math>階的低通濾波器<math>g_j[n]</math>表示成:<br />
<math>g_j[n]=Zg_{j-1}[n]=Z^2g_{j-2}[n]=...=Z^{j-1}g_{1}[n]</math><br />
注意經過提升取樣後，<math>g_j[n]</math>的序列長度為原來<math>g_{j-1}[n]</math>的兩倍。<br />


# 原始信號與高通濾波器做[[旋積|旋積分]]之後會得到此信號中高頻的成分。此高頻的成分為第一個高頻的輸出。<br />
# 原始信號與低通濾波器做[[旋積|旋積分]]後會得到信號中低頻的成分，此低頻的成分再作為下一階濾波器的輸入。

重複上述兩個步驟，即可將信號作多階的穩定小波轉換。<br />
<br />

而經過<math>j</math>組的高通濾波器和低通濾波器組合之後，第<math>j</math>階的結果:<br />

<math>x_{j,L[n]}=x_{j-1,L[n]}*g_j[n]</math><br />
<math>x_{j,H[n]}=x_{j-1,L[n]}*h_j[n]</math><br />

== Matlab 使用範例 ==
下面是利用[[Matlab]]的離散小波轉換的函式，稍作調整後的一維穩定小波轉換範例:
<syntaxhighlight lang=Matlab>
[tmpAPP,tmpDET] = 
dwt(A(j,ε1, ,ɛj),wname,'mode','per','shift',ɛj+1); 
A(j+1,ɛ1, ,ɛj,ɛj+1) = wshift('1D',tmpAPP,ɛj+1);
D(j+1,ɛ1, ,ɛj,ɛj+1) = wshift('1D',tmpDET,ɛj+1);
</syntaxhighlight>

參考:[http://www.mathworks.com/help/wavelet/ug/discrete-stationary-wavelet-transform-swt.html MatlabWorks-Discrete Stationary Wavelet Transform (SWT)] {{Wayback|url=http://www.mathworks.com/help/wavelet/ug/discrete-stationary-wavelet-transform-swt.html |date=20181025053329 }}

== 應用 ==
穩定小波轉換在訊號處理上有一些應用：
* 降低信號雜訊(Signal [http://en.wikipedia.org/wiki/Denoising denoising])
* [[樣式辨認]]
* 信號抵達方向性([http://en.wikipedia.org/wiki/Direction_of_arrival Direction of Arrival] {{Wayback|url=http://en.wikipedia.org/wiki/Direction_of_arrival |date=20210225003453 }}, DOA)估計
* 信號重建(Signal regeneration)

== 同義轉換 ==
以下的幾種轉換或演算，皆為略過離散小波轉換的縮減取樣步驟，只是隨著提出的時間而有相異的名字
* 穩定小波轉換 (Stationary Wavelet Transform)
* 冗餘小波轉換 (Redundant Wavelet Transform)
* à trous演算法 (Algorithme à trous)
* 準連續小波轉換 (Quasi-continuous wavelet transform)
* 平移不變量小波轉換 (Translation invariant wavelet transform)
* 轉移不變量小波轉換 (Shift invariant wavelet transform)
* 循環平移演算法 (Cycle spinning)
* 最大重複離散小波轉換 (Maximal overlap discrete wavelet transform, MODWT)
* 非抽樣小波轉換 (Undecimated wavelet transform, UWT)

== 參考文獻 ==
* P. P. Viadyanathan, ''Multirate Systems and Filter Banks'', Prentice Hall, 1993, ISBN 0-13-605718-7
* G. P. Nason and B. W. Silverman, ''The stationary wavelet transform and some statistical applications'',  Lecture Notes in Statistics
* M.V. Tazebay and A.N. Akansu, ''Progressive Optimality in Hierarchical Filter Banks,'' Proc. IEEE International Conference on Image Processing (ICIP), Vol 1, pp. 825-829, Nov. 1994
* P. Dutilleux, ''An implementation of the algorithme à trous to compute the wavelet transform,'' in Wavelets: Time-Frequency Methods and Phase Space, J.-M. Combes, A. Grossman, and P. Tchamitchian,Eds. Berlin, Germany: Springer-Verlag, 1989, pp. 298–304.


[[Category:小波分析|B]]