新起点
稳定小波转换
2020-06-15 11:32:34

稳定小波转换(Stationary Wavelet Transform, SWT)是小波分析(Wavelet Analysis)的一种转换,为离散小波变换(Discrete Wavelet Transform, DWT)的变形。

稳定小波转换可以弥补离散小波转换因为缩减取样(Downsampling)而失去的平移不变性(Translation-invariant)。稳定小波转换不同于离散小波转换的部分,主要在于经过每一阶的高通滤波器和低通滤波器之后,是将滤波器提升取样(Upsampling),取代离散小波转换在经过滤波器之后的缩减取样。

稳定小波转换是做数据和信号的分析一种很好的工具,尽管它的运算量会因为没有缩减取样而较离散小波转换多一些,但其具有平移不变性,且只需将离散小波转换在滤波器的设计上做些微的修改即可实现。

下图是稳定小波转换的数位实现模型

每一组高通滤波器和低通滤波器皆为提升取样后的前一组高通滤波器及低通滤波器,可以下图表示:


以数学形式来呈现稳定小波转换滤波器提升取样的设计概念:

Z {\displaystyle Z} 为一将 0 {\displaystyle 0} 加入 x {\displaystyle x} 序列的运算,

( Z x ) k {\displaystyle (Z_{x})_{k}} = x k / 2 , {\displaystyle =x_{k/2},} for k = 2 , 4 , 6 , 8... {\displaystyle k=2,4,6,8...}

( Z x ) k {\displaystyle (Z_{x})_{k}} = 0 , {\displaystyle =0,} for k = 1 , 3 , 5 , 7... {\displaystyle k=1,3,5,7...}

我们可以将第 j {\displaystyle j} 阶的高通滤波器 h j {\displaystyle h_{j}} 表示成:
h j = Z h j 1 = Z 2 h j 2 = . . . = Z j 1 h 1 {\displaystyle h_{j}=Zh_{j-1}=Z^{2}h_{j-2}=...=Z^{j-1}h_{1}}
同样的,我们可以将第 j {\displaystyle j} 阶的低通滤波器 g j {\displaystyle g_{j}} 表示成:
g j = Z g j 1 = Z 2 g j 2 = . . . = Z j 1 g 1 {\displaystyle g_{j}=Zg_{j-1}=Z^{2}g_{j-2}=...=Z^{j-1}g_{1}}
注意经过提升取样后, g j {\displaystyle g_{j}} 的序列长度为原来 g j 1 {\displaystyle g_{j-1}} 的两倍。


重复上述两个步骤,即可将信号作多阶的稳定小波转换。

而经过 j {\displaystyle j} 组的高通滤波器和低通滤波器组合之后,第 j {\displaystyle j} 阶的结果:

x j , L = x j 1 , L g j {\displaystyle x_{j,L}=x_{j-1,L}*g_{j}}
x j , H = x j 1 , L h j {\displaystyle x_{j,H}=x_{j-1,L}*h_{j}}

下面是利用Matlab的离散小波转换的函式,稍作调整后的一维稳定小波转换范例:

 = 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)

稳定小波转换在讯号处理上有一些应用:

以下的几种转换或演算,皆为略过离散小波转换的缩减取样步骤,只是随着提出的时间而有相异的名字

相关:

网站公告: