小波阈值自适应收缩方法、系统、电子设备及存储介质与流程

1.本发明涉及一种小波阈值自适应收缩方法、系统、电子设备及存储介质，属于信号分析降噪技术领域。


背景技术：

2.长久以来，傅里叶变换被广泛应用于对信号进行相关的处理，并且也取得了一系列的成果。但是傅里叶变换在信号去噪方面也有很多的局限性，其中傅里叶变换在处理信号去噪的一个缺陷就是，用傅里叶进行分析时，其构造函数是周期性的正弦波和余弦波。鉴于其局限性，傅里叶变换只适合具有周期性或具有近似周期性的信号进行滤波和压缩，而对于具有非周期性或局部特征很明显的信号处理上效果不是很好。
3.傅里叶变换在信号去噪方面存在局限性，由其发展来的小波变换(wavelet transform,wt)被广泛认为是多分辨率信号分析、信号检测、图像边缘检测、图像压缩和信号去噪的关键工具。作为在信号处理领域中的一种新的分析方法，wt不仅保留了傅里叶变换的许多优点，而且在原来的基础上进行了改进和发展，使其能够在时频域对信号进行处理。wt的显著优点是通过变换可以将信号进行更细微的处理，并且能够将信号的某些特征较好的表现出来，实现了在时频域对信号进行局部化，多尺度的分析要求。在小波基础上发展来的信号去噪方法表现出了良好的去噪效果，是傅里叶变换在信号处理领域的完善和发展。
4.虽然wt有诸多的优点，但在小波收缩降噪中的阈值若选择选取不当，则去噪效果会变差。


技术实现要素：

5.本发明的目的在于，克服现有技术存在的技术缺陷，解决上述技术问题，提出一种小波阈值自适应收缩方法、系统、电子设备及存储介质，采用统计过程控制收缩进行阈值的选择，为此引入了一系列控制上下限，允许迭代丢弃小波系数，直到所有小波系数都在指定的控制范围内。
6.小波阈值自适应收缩方法，包括如下步骤：
7.步骤ss1：利用小波变换对小波y进行阈值运算来估计真实信号；
8.步骤ss2：利用控制上下限进行阈值估计，判断样本y(i)的当前值是否在控制范围内；
9.步骤ss3：运用统计过程控制收缩，允许迭代丢弃“不良”小波系数，直到所有小波系数都在指定的控制范围内；
10.步骤ss4：引入新参数α
1*
来优化滤波效果；
11.所述步骤ss1具体包括：
12.假设有n个嘈杂的观察值y＝[y[0] y[1]
ꢀ…ꢀ
y[n
‑
1]]
t
，n个未知信号x＝[x[0] x[1]
ꢀ…ꢀ
x[n
‑
1]]
t
，利用y＝x+n来估计，其中n＝[n[0] n[1]
ꢀ…ꢀ
n[n
‑
1]]
t
是均值为零且方差
为零的高斯白噪声(wgn)σ2(n(0,σ2))；
[0013]
在多分辨率小波分析中，n＝2
j
，j是小波分解级的最大数；令w为与给定多分辨率小波分解相关的正交变换矩阵，y的小波表示为：w＝w
·
y，w是线性转化矩阵；
[0014]
令c＝w
·
y和z＝w
·
y，得到w＝c+z；
[0015]
再将w修改为w＝t(w,λ)，其中t(
·
)是阈值函数，λ>0是阈值,w小于λ的元素被消除或平滑，如公式(1)和(2)所示：
[0016][0017][0018]
其中符号(
·
)是符号函数，向量w，c和z的条目被双重索引并表示为w
j,k
，c
j,k
，z
j,k
，j＝1,2,
…
,j表示与频率相关的缩放域索引，k＝1,2,
…
,2
j
‑
j
表示时域索引；最终根据的收缩系数来估计真实信号。
[0019]
所述步骤ss2具体包括：
[0020]
假设质量测量值为正态分布，平均值为μ0，标准偏差为σ；y＝[y[0] y[1]
ꢀ…ꢀ
y[n
‑
1]]
t
是具有n个观测值的向量，i＝0,1,
…
n
‑
1，任何样本y(i)的概率为1
‑
α，处于以下范围内：
[0021]
μ0‑
dσ
ꢀꢀꢀ
(3)
[0022]
μ0+dσ
ꢀꢀꢀ
(4)
[0023]
其中，d是标准正态分布的(1
‑
α/2)分位数或erfc
‑1(
·
)是逆互补误差函数，控制上下限是：
[0024]
lcl＝μ0‑
dσ
ꢀꢀꢀ
(5)
[0025]
cl＝μ0ꢀꢀꢀ
(6)
[0026]
ucl＝μ0+dσ
ꢀꢀꢀ
(7)
[0027]
其中，统计平均值的中心线为cl，统计平均值的两条水平线，分别称为控制下限lcl和控制上限ucl，数量d解释为上下限距其中心线的“距离”，以标准偏差单位表示。
[0028]
所述步骤ss2具体包括：
[0029]
如果y(i)的当前值在控制范围内，则该过程被认为是“受控的”，也就是说，这是平均值为0的正态分布变量的出现；否则，如果y(i)超出了控制范围，那么得出的结论是该过程是“失控”，也就是说，这是由于平均值不同的随机变量μ1≠μ0；因此，选择控制上下限等同于设置用于测试假设的关键区域：
[0030]
h0:μ＝μ0(受控)
ꢀꢀꢀ
(8)
[0031]
h1:μ≠μ0(失控)
ꢀꢀꢀ
(9)
[0032]
控制上下限针对观察到的过程每个观察值y(i)反复测试该假设；假设测试的一般程序从i型错误α的规范开始，然后是设计一种使测试的能力最大化的测试过程，α＝pr(reject h0|h
0 is true)和power＝pr(reject h0|h
0 is false)，通过选择α控制错误警报的可能性。
[0033]
所述步骤ss3具体包括：
[0034]
建立以下假设检验：
[0035]
h0:c
j,k
＝0
ꢀꢀꢀ
(10)
[0036]
h1:c
j,k
≠0
ꢀꢀꢀ
(11)
[0037]
基于统计过程控制方法和公式(5)—(7)评估给定小波系数w
j,k
在上下限内的概率，因此，对于每个分解级别j都有：
[0038]
lcl＝
‑
d
j
·
s
j
ꢀꢀꢀ
(12)
[0039]
cl＝0
ꢀꢀꢀ
(13)
[0040]
ucl＝d
j
·
s
j
ꢀꢀꢀ
(14)
[0041]
其中s
j
是级别j的小波系数的校正样本标准偏差，而d
j
是分解级别j的控制上下限的距离，在受控范围内：
[0042]
pr{|w
j,k
|≤d
j
·
s
j
}＝1
‑
α
j
ꢀꢀꢀ
(15)
[0043]
其中α
j
是规定的显著性水平，α
j
是不存在信号时检测到信号的概率。
[0044]
所述步骤ss3具体包括：
[0045]
在给定上下限距离d
j
的情况下，通过以下描述的迭代过程给出用于估计标度级别j的阈值λ
j
的收缩方法：
[0046]
步骤ss31，估计第j级小波系数的标准偏差：
[0047][0048]
其中n
j
＝2
j
‑
j
是多分辨率级别j上的小波系数的数量和步骤ss32，根据以下条件建立控制上下限：
[0049]
lcl＝
‑
d
j
·
s
j
ꢀꢀꢀ
(17)
[0050]
ucl＝d
j
·
s
j
ꢀꢀꢀ
(18)
[0051]
步骤ss33，如果小波系数w
j,k
(k＝1,2
…
n
j
)超过了区间[lcl,ucl]，则排除该系数，并重复步骤ss31和步骤ss32；如果所有w
j,k
都在区间[lcl,ucl]内，则阈值为λ
j
＝d
j
·
s
j
，然后停止迭代方法。
[0052]
所述步骤ss4具体包括：
[0053]
为获得显着性水平的最佳值，引入以下优化：
[0054][0055]
其中x
i
是输入信号，是根据提出的方案的相关降噪信号，error(,)是评估降噪信号的品质因数，m是信号实例的数量，a是搜索空间；对于误差度量error，采用经典信噪比snr的负值，即：
[0056][0057]
其中是x
i
的方差和是的方差，snr以分贝db表示，其中snr值越大表示滤波效果越好。
[0058]
本发明提出小波阈值自适应收缩系统，包括：
[0059]
阈值运算模块，用于执行：利用小波变换对小波y进行阈值运算来估计真实信号；
[0060]
阈值估计模块，用于执行：利用控制上下限进行阈值估计，判断样本y(i)的当前值是否在控制范围内；
[0061]
迭代模块，用于执行：运用统计过程控制收缩，允许迭代丢弃“不良”小波系数，直到所有小波系数都在指定的控制范围内；
[0062]
滤波优化模块，用于执行：引入新参数α
1*
来优化滤波效果。
[0063]
本发明还提出电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本发明所述方法的步骤。
[0064]
本发明还提出存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明所述方法的步骤。
[0065]
优点和效果：
[0066]
本发明所达到的有益效果：本发明提供了一种小波阈值自适应收缩算法、系统、电子设备及存储介质，并通过小波变换来识别和删除“不良”数据，以此来改善去噪效果，本发明采用统计过程控制收缩进行阈值的选择，为此引入了一系列控制上下限，允许迭代丢弃小波系数，直到所有小波系数都在指定的控制范围内，通过在蒙特卡洛模拟中得到评估，信噪比(snr)、信噪比增益和均方根误差(rmse)均为优值。
附图说明
[0067]
图1是本发明的优选实施例的流程图。
具体实施方式
[0068]
下面结合附图对本发明作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。
[0069]
实施例1：如图1所示，小波阈值自适应收缩方法，包括如下步骤：
[0070]
步骤ss1：利用小波变换对小波y进行阈值运算来估计真实信号；
[0071]
步骤ss2：利用控制上下限进行阈值估计，判断样本y(i)的当前值是否在控制范围内；
[0072]
步骤ss3：运用统计过程控制收缩，允许迭代丢弃“不良”小波系数，直到所有小波系数都在指定的控制范围内；
[0073]
步骤ss4：引入新参数α
1*
来优化滤波效果。
[0074]
作为一种较佳的实施例，所述步骤ss1具体包括：
[0075]
假设有n个嘈杂的观察值y＝[y[0] y[1]
ꢀ…ꢀ
y[n
‑
1]]
t
，n个未知信号x＝[x[0] x[1]
ꢀ…ꢀ
x[n
‑
1]]
t
，利用y＝x+n来估计，其中n＝[n[0] n[1]
ꢀ…ꢀ
n[n
‑
1]]
t
是均值为零且方差为零的高斯白噪声(wgn)σ2(n(0,σ2))；
[0076]
在多分辨率小波分析中，n＝2
j
，j是小波分解级的最大数；令w为与给定多分辨率小波分解相关的正交变换矩阵，y的小波表示为：w＝w
·
y，w是线性转化矩阵；
[0077]
令c＝w
·
y和z＝w
·
y，得到w＝c+z；
[0078]
再将w修改为w＝t(w,λ)，其中t(
·
)是阈值函数，λ>0是阈值,w小于λ的元素被消除
或平滑，如公式(1)和(2)所示：
[0079][0080][0081]
其中符号(
·
)是符号函数，向量w，c和z的条目被双重索引并表示为w
j,k
，c
j,k
，z
j,k
，j＝1,2,
…
,j表示与频率相关的缩放域索引，k＝1,2,
…
,2
j
‑
j
表示时域索引；最终根据的收缩系数来估计真实信号。
[0082]
作为一种较佳的实施例，所述步骤ss2具体包括：
[0083]
假设质量测量值为正态分布，平均值为μ0，标准偏差为σ；y＝[y[0] y[1]
ꢀ…ꢀ
y[n
‑
1]]
t
是具有n个观测值的向量，i＝0,1,
…
n
‑
1，任何样本y(i)的概率为1
‑
α，处于以下范围内：
[0084]
μ0‑
dσ
ꢀꢀꢀ
(3)
[0085]
μ0+dσ
ꢀꢀꢀ
(4)
[0086]
其中，d是标准正态分布的(1
‑
α/2)分位数或erfc
‑1(
·
)是逆互补误差函数，控制上下限是：
[0087]
lcl＝μ0‑
dσ
ꢀꢀꢀ
(5)
[0088]
cl＝μ0ꢀꢀꢀ
(6)
[0089]
ucl＝μ0+dσ
ꢀꢀꢀ
(7)
[0090]
其中，统计平均值的中心线为cl，统计平均值的两条水平线，分别称为控制下限lcl和控制上限ucl，数量d解释为上下限距其中心线的“距离”，以标准偏差单位表示。
[0091]
作为一种较佳的实施例，所述步骤ss2具体包括：
[0092]
如果y(i)的当前值在控制范围内，则该过程被认为是“受控的”，也就是说，这是平均值为0的正态分布变量的出现；否则，如果y(i)超出了控制范围，那么得出的结论是该过程是“失控”，也就是说，这是由于平均值不同的随机变量μ1≠μ0；因此，选择控制上下限等同于设置用于测试假设的关键区域：
[0093]
h0:μ＝μ0(受控)
ꢀꢀꢀ
(8)
[0094]
h1:μ≠μ0(失控)
ꢀꢀꢀ
(9)
[0095]
控制上下限针对观察到的过程每个观察值y(i)反复测试该假设；假设测试的一般程序从i型错误α的规范开始，然后是设计一种使测试的能力最大化的测试过程，α＝pr(reject h0|h
0 is true)和power＝pr(reject h0|h
0 is false)，通过选择α控制错误警报的可能性。
[0096]
作为一种较佳的实施例，所述步骤ss3具体包括：
[0097]
建立以下假设检验：
[0098]
h0:c
j,k
＝0
ꢀꢀꢀ
(10)
[0099]
h1:c
j,k
≠0
ꢀꢀꢀ
(11)
[0100]
基于统计过程控制方法和公式(5)—(7)评估给定小波系数w
j,k
在上下限内的概率，因此，对于每个分解级别j都有：
[0101]
lcl＝
‑
d
j
·
s
j
ꢀꢀꢀ
(12)
[0102]
cl＝0
ꢀꢀꢀ
(13)
[0103]
ucl＝d
j
·
s
j
ꢀꢀꢀ
(14)
[0104]
其中s
j
是级别j的小波系数的校正样本标准偏差，而d
j
是分解级别j的控制上下限的距离，在受控范围内：
[0105]
pr{|w
j,k
|≤d
j
·
s
j
}＝1
‑
α
j
ꢀꢀꢀ
(15)
[0106]
其中α
j
是规定的显著性水平，α
j
是不存在信号时检测到信号的概率。
[0107]
作为一种较佳的实施例，所述步骤ss3具体包括：
[0108]
在给定上下限距离d
j
的情况下，通过以下描述的迭代过程给出用于估计标度级别j的阈值λ
j
的收缩方法：
[0109]
步骤ss31，估计第j级小波系数的标准偏差：
[0110][0111]
其中n
j
＝2
j
‑
j
是多分辨率级别j上的小波系数的数量和步骤ss32，根据以下条件建立控制上下限：
[0112]
lcl＝
‑
d
j
·
s
j
ꢀꢀꢀ
(17)
[0113]
ucl＝d
j
·
s
j
ꢀꢀꢀ
(18)
[0114]
步骤ss33，如果小波系数w
j,k
(k＝1,2
…
n
j
)超过了区间[lcl,ucl]，则排除该系数，并重复步骤ss31和步骤ss32；如果所有w
j,k
都在区间[lcl,ucl]内，则阈值为λ
j
＝d
j
·
s
j
，然后停止迭代方法。
[0115]
作为一种较佳的实施例，所述步骤ss4具体包括：
[0116]
为获得显着性水平的最佳值，引入以下优化：
[0117][0118]
其中x
i
是输入信号，是根据提出的方案的相关降噪信号，error(,)是评估降噪信号的品质因数，m是信号实例的数量，a是搜索空间；对于误差度量error，采用经典信噪比snr的负值，即：
[0119][0120]
其中是x
i
的方差和是的方差，snr以分贝db表示，其中snr值越大表示滤波效果越好。
[0121]
实施例2：本发明提出小波阈值自适应收缩系统，包括：
[0122]
阈值运算模块，用于执行：利用小波变换对小波y进行阈值运算来估计真实信号；
[0123]
阈值估计模块，用于执行：利用控制上下限进行阈值估计，判断样本y(i)的当前值是否在控制范围内；
[0124]
迭代模块，用于执行：运用统计过程控制收缩，允许迭代丢弃“不良”小波系数，直
到所有小波系数都在指定的控制范围内；
[0125]
滤波优化模块，用于执行：引入新参数α
1*
来优化滤波效果。
[0126]
实施例3：本发明还提出电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现本发明所述方法的步骤。
[0127]
实施例4：本发明还提出存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现本发明所述方法的步骤。
[0128]
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。