稳定小波转换
 |
稳定小波转换（英语：Stationary wavelet transform，缩写：SWT）是小波分析的一种转换，为离散小波变换（DWT）的变形。
稳定小波转换可以弥补离散小波转换因为降采样而失去的平移不变性。稳定小波转换不同于离散小波转换的部分，主要在于经过每一阶的高通滤波器和低通滤波器之后，是将滤波器升采样，取代离散小波转换在经过滤波器之后的缩减取样。
稳定小波转换是做数据和信号的分析一种很好的工具，尽管它的运算量会因为没有缩减取样而较离散小波转换多一些，但其具有平移不变性，且只需将离散小波转换在滤波器的设计上做些微的修改即可实现。
实现方式
[编辑]下图是稳定小波转换的数位实现模型
每一组高通滤波器和低通滤波器皆为提升取样后的前一组高通滤波器及低通滤波器，可以下图表示：
以数学形式来呈现稳定小波转换滤波器提升取样的设计概念:
令为一将加入序列的运算，
for
for
我们可以将第阶的高通滤波器表示成:
同样的，我们可以将第阶的低通滤波器表示成:
注意经过提升取样后，的序列长度为原来的两倍。
![{\displaystyle h_{j}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cdf51ad18cf67c6cd57eb94bc4f397ff8359c82e)
![{\displaystyle h_{j}[n]=Zh_{j-1}[n]=Z^{2}h_{j-2}[n]=...=Z^{j-1}h_{1}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4bd11a2ff4919a39f904c1b4eded0ca8051e3bdb)
![{\displaystyle g_{j}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8914b9f7d56b7411e1e4e051d8273c99895e40a6)
![{\displaystyle g_{j}[n]=Zg_{j-1}[n]=Z^{2}g_{j-2}[n]=...=Z^{j-1}g_{1}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a7d374205c74b24bd5c340d9e3254bf65c077f51)
![{\displaystyle g_{j-1}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fff21ca7745ea0d1bf8bafffbf0195351c00b541)
重复上述两个步骤，即可将信号作多阶的稳定小波转换。
而经过组的高通滤波器和低通滤波器组合之后，第阶的结果:
![{\displaystyle x_{j,L[n]}=x_{j-1,L[n]}*g_{j}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0b9f2564920d02eafba7a190e8b62c60a3047fbd)
![{\displaystyle x_{j,H[n]}=x_{j-1,L[n]}*h_{j}[n]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/71444afe526ef781ad565523440560e0f9caf117)
Matlab 使用范例
[编辑]下面是利用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);
参考:MatlabWorks-Discrete Stationary Wavelet Transform (SWT) （页面存档备份，存于互联网档案馆）
应用
[编辑]稳定小波转换在讯号处理上有一些应用：
- 降低信号噪声(Signal denoising)
- 样式辨认
- 信号抵达方向性(Direction of Arrival （页面存档备份，存于互联网档案馆）, 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.