第六章小波变换的几个典型应用.docx
《第六章小波变换的几个典型应用.docx》由会员分享，可在线阅读，更多相关《第六章小波变换的几个典型应用.docx（27页珍藏版）》请在冰豆网上搜索。
第六章小波变换的几个典型应用
第六章小波变换的几个典型应用
6.1小波变换与信号处理
小波变换作为信号处理的一种手段，逐渐被越来越多领域的理论工作者和工程技术人员所重视和应用，并在许多应用中取得了显著的效果。
同传统的处理方法相比，小波变换取得了质的飞跃，在信号处理方面具有更大的优势。
比如小波变换可以用于电力负载信号的分析与处理，用于语音信号的分析、变换和综合，还可以检测噪声中的未知瞬态信号。
本部分将举例说明。
6.1.1小波变换在信号分析中的应用
［例6-1］以含躁的三角波与正弦波的组合信号为例具体说如何利用小波分析来分析信号。
已知信号的表达式为
t-sin（0.3t）b（t）1红乞500
“、J500
s（t）
1000"+sin（0.3t）+b（t）501兰t<000
.500
应用db5小波对该信号进行7层分解。
xiaobo0601.m
图6-1含躁的三角波与正弦波混合信号波形
分析：
（1）在图6-2中，逼近信号a7是一个三角波。
（2）在图6-3中细节信号d1和d2是与噪声相关的，而d3（特别是d4）与正弦信号相关。
01002003004005006007008009001000
2
0
-2
0
2r-
0、
-2
100
200
300
400
500
600
700
800
900
1000
4
3
100
200
300
400
500
600
700
800
900
1000
0
100
200
300
400
500
600
700
800
900
1000
100
200
300
400
500
600
700
800
900
1000
样本序号n
图6-2小波分解后各层逼近信号
-1
6
7
01002003004005006007008009001000
01002003004005006007008009001000
0
01002003004005006007008009001000
1
-5
01002003004005006007008009001000
样本序号n
图6-3小波分解后各层细节信号
6・1・2小波变换在信号降躁和压缩中的应用
一、信号降躁
1•工程中，有用信号一般是一些比较平稳的信号，噪声通常表现为高频信号。
2•消躁处理的方法：
首先对信号进行小波分解，由于噪声信号多包含在具有较高频率的细节中，我们可以利用门限、阈值等形式对分解所得的小波系数进行处理，然后对信号进行小波重构即可达到对信号的消躁目的。
小波分析进行消躁处理的3种方法：
（1）默认阈值消躁处理。
该方法利用ddencmp生成信号的默认阈值，然后利用wdencmp函数进行消躁处理。
（2）给定阈值消躁处理。
在实际的消躁处理过程中，阈值往往可通过经验公
式获得，且这种阈值比默认阈值的可信度高。
在进行阈值量化处理时可利用函数wthresh。
（3）强制消躁处理。
该方法时将小波分解结构中的高频系数全部置为0,即滤掉所有高频部分，然后对信号进行小波重构。
方法简单，消躁后信号比较平滑，但易丢失信号中的有用成分。
小波阈值去噪方法是目前应用最为广泛的小波去噪方法之一。
3•信号降噪的准则：
1•光滑性：
在大部分情况下，降噪后的信号应该至少和原信号具有同等的
光滑性。
2•相似性：
降噪后的信号和原始信号的方差估计应该是最坏情况下的方差最小。
4•一维信号消躁的步骤：
（1）一维信号的小波分解。
选择一个小波并确定分解的层次，然后进行分解计算。
（2）小波分解高频系数的阈值量化。
对各个分解尺度下的高频系数选择一个阈值进行软阈值量化处理。
（3）一维小波重构。
根据小波分解的最低层系数和各层高频系数进行一维小波重构。
关键：
如何选择阈值和进行阈值量化。
在某种程度上，它关系到信号消躁的质量5•消躁阈值选取规则
1叭k.
«j,k
10,
叫,k
>&
＜扎
■j,k
硬阈值法：
国j,k
软阈值法：
几阿
■j,k
/
丿1
3
A
图（a）硬阈值
图（b）软阈值
图6-4估计小波系数的软阈值与硬阈值方法
图6-4表明了软阈值和硬阈值法的区别，图中横坐标表示小波分解系数-,
纵坐标表示由阈值法得到的小波系数估计值?
，■为阈值。
可以看出，硬阈值法的？
函数在•点处不连续，这会给重构信号带来震荡；软阈值法虽然?
函数连续性较好，但其导数并不连续，这就限制了它的进一步应用。
并且当■-■时，由
软阈值法得出的估计值?
与小波系数••存在着恒定的偏差。
这些分析表明，软阈值法通常会使去噪后的信号平滑一些，但是也会丢掉某些特征；而硬阈值可以保留信号的特征，但是在平滑方面有所欠缺。
一般来说，去噪中软阈值的作用会更多一些，但是到底选取哪种处理方法，还应视具体情况而定。
6•应用一维小波分析进行信号消躁处理的MATLAB函数
小波函数：
wden和wdencmp
［例6-2］利用小波分析对含躁正弦波进行消躁。
xiaobo0602.m分析：
（1）消躁后的信号大体上恢复了原信号的形状，并明显去除了噪声所引起的干扰。
（2）恢复后的信号与原信号相比有明显的改变。
主要原因是，在进行消躁处理的过程中所用的分析小波和细节系数阈值不恰当。
対号信始ML
A值胴号信躁0A
2
0
值
幅畐
号信躁
2
0
-2
0-■-/■,■'\\-
\\\\/\
-11_I1'、』I111i-
01002003004005006007008009001000
样本序号n
01002003004005006007008009001000
样本序号n
默认阈值消躁后的信号
强制消躁后的信号
600
600
A400
A
J\J
值
**\J
L/1
幅200
t/
A400
值
幅200
0
050010001500
样本序号n
5001000
样本序号n
0
1002003004005006007008009001000
样本序号n
［例6-3］在电网电压值监测过程中，由于监测设备出现了一点故障，致使所采集到的信号受到噪声的污染。
现在利用小波分析对污染信号进行消躁处理以恢复原始信号。
600
A400
值
幅200
0
050010001500
样本序号n
给定软阈值消躁后的信号
分析:
（1）强制消躁处理后信号比较光滑，但可能丢失有用信息。
（2）默认阈值消躁和给定软阈值消躁这两种处理方法在实际中应用的更广泛。
阈值函数图形如下：
xiaobo0604.m
原始信号硬阈值信号软阈值信号
1
0.8
0.6
0.4
--
1
1
/
-
1
0.8
0.6
0.4
f
1■
/
1
1
■
0.6
0.4
0.2
-
1
1二
1
0.2
'-
0.2
-
0
r
1
A
值
0
A
--值
0
A
值
■|'
1r/
r
幅
幅
幅
-0.2
-0.2
-
-0.2
/
I
-0.4
/
B->■
1
-0.4
/
-
-0.4
Jr
1
・i
一
-0.6
-0.6
-f
-
-0.6
-
-
-0.8
./■
-0.8
■
-1
1
-1
1
-0.8
1
050100050100050100
样本序号n样本序号n样本序号n
二、信号压缩
1.压缩依据：
一个比较规则的信号是由一个数据量很小的低频系数和几个高频系数所组成的。
这里对低频系数的选择有一个要求，即需要在一个合适的分解层上选择低频系数。
2.压缩手段：
小波分析和小波包分析两种手段。
3•压缩步骤：
（1）信号的小波（包）分解。
（2）对高频系数进行阈值量化处理。
对第1层到第N层的高频系数，均可选择不同的阈值，并且用硬阈值进行系数的量化。
（3）对量化后的系数进行小波（包）重构。
4.两种比较有效的信号压缩方法：
第一种方法：
对信号进行小波尺度的扩展，并且保留绝对值最大的系数。
在这种情况下，可以选择全局阈值，此时仅需输入一个参数即可。
第二种方法：
根据分解后各层的效果来确定某一层的阈值，且每一层的阈值可以互不相同。
［例6-4］利用小波分析对给定信号进行压缩处理。
xiaobo0605.m
6.2小波变换在电力负载信号的应用
电力系统在线检测信号含有大量的现场背景噪声，给传统方式的数据采集与故障诊断带来很大的困难。
将以处理瞬态信号、含宽带噪声信号等见长的小波分析应用于电力系统在线监测是大有前途的。
本小节的测量数据是从一个复杂的设备上采集的电力负载信号，每分钟采集
一个样本，持续了5个星期，总共50400个数据样本。
测量数据受到传感器误差和状态噪声两种噪声的影响。
本小节将分析其中的两段数据，其中第一段是上午12：
30至下午1:
00间采集的样本，由于这段时间处于用电高峰，因此数据很复杂；第二段是下半夜采集的样本，数据比较简单。
一、信号分解
［例6-5］利用小波分解分析第一段数据的信号成分。
xiaobo0606.m
340
335
330
325
305-
300-
320-
315-
310-
295111*—
36003610362036303640
365036603670368036903700
样本序号n
350
300
250
3700
20
4
3600
3650
360036503700
40
d
-20
5
0
-5
360036503700
3600
3650
3700
样本序号n
图2
分析：
第一段电力载波信号如图1所示，利用db3小波对其进行5层小波分解,得到逼近信号和细节信号如图2所示。
可以看出：
（1）细节信号di和d2的值较小，可以认为是由传感器和状态噪声的高频分量引起的局部干扰；
（2）细节信号d4包含了3个相连的主要信号模式，它最接近于原始数据的曲线；
（3）细节信号d5含有的信息不多，因此第4层贡献最大，它提取了原始数据曲线的形状。
二、暂态信号检测
为保证电力系统的安全可靠运行，必须对电力设备进行状态监测根据电力信号来判别其运行的状态。
电力系统暂态故障信号往往在故障时刻发生突变，若能捕获设备故障信息突变时刻和大小，有利于在故障初期及早米取措施使系统恢复正常，这对提高设备运行可靠性具有重要意义。
［例6-6］利用小波分解分析检测第二段信号的突变点成分。
xiaobo0607.m
310
300
290
280
270
A
值260
幅
250
240
230
220
210
156015801600162016401660168017001720
样本序号n
分析：
利用db3小波对其进行5层分解，得到逼近信号和细节信号如图所示。
可
以看出：
由细节信号d2可以检测突变点位置t=1625，由细节信号d1也能隐约看出t=1600处的突变点。
400
300
200
4
5
0
-5
15501600165017001750
400
300
200
15501600165017001750
5
0
-5
15501600165017001750
20
0
-20
15501600165017001750
400
300
200
15501600165017001750
样本序号n
1
样本序号n
三、传感器故障检测
［例6-7］利用小波分析检测传感器故障。
xiaobo0608.m
A值幅
22002400260028003000320034003600
样本序号n
50
-50
样本序号n
利用db3小波对信号进行5层分解，得到第1〜3层细节信号如图所示。
可以看出每个细节信号都显示了在t=2400〜t=3600之间的信号由于传感器故障而引入了传感器误差噪声。
四、奇异点定位消除
［例6-8］利用小波分析检测信号中的奇异点并消除。
xiaobo0609.m
由原始信号波形可以看出在t=1193和t=1215两处存在奇异值点。
进一步利用db3小波对信号进行5层分解，得到第1、2、3层细节信号如图所示。
发现奇异值点包含在细节信号d1和d2中，且与原信号中的奇异点是同步的。
为了消除奇异点，重构信号时令细节信号d1、d2和d3等于零，得到的波形如图所示，比较可见奇异值点已经很不明显了。
390--
380-,-
370--
Hb/\
A
值360--
幅I||’!
\/xA
350「LI\||-
if、
340-'..■-
\
V■飞I
330-“'>卜、一-
320IIIII|||
116011701180119012001210122012301240
样本序号n
图原信号
116011701180119012001210122012301240
样本序号n
图小波分解的细节信号
365
360
355
350
345
340
335
330
325
11601170118011901200121012201230
样本序号n
320
1240
图消除奇异点后的波形
6.3小波分析在图像消躁中的应用
图像消躁在信号处理中是一个经典问题，传统的消躁方法是采用平均或线性方法进行，常采用的是维纳滤波，但是消躁效果不好。
随着小波理论日益完美，它以自身良好的时频特性在图像消躁领域受到越来越多的关注，开辟了用非线性
方法消躁的先河。
具体说来，小波能够消躁主要得意于小波变换具有如下特点：
（1）低熵性：
小波系数的稀疏分布，使图像变换后的熵降低。
（2）多分辨率特性：
由于采用多分辨率的方法，所以可以非常好地刻画信号的非平稳性，如突变和断点等，可以在不同分辨率下根据信号和噪声的分布来去除噪声。
（3）去相关性：
小波变换可以对信号去相关，且噪声在变换后有白化趋势，所以小波域比时域更利于去躁。
（4）基函数选择灵活：
小波变换可以灵活选择基函数，也可根据信号特点和消躁要求选择多带小波、小波包等，对不同场合，可以选择不同的小波母函数。
一、小波图像消躁的基本原理
常用的图像消躁方法是小波阈值消躁方法，它是一种实现简单而效果好的消躁方法。
阈值消躁方法的思想很简单，就是对小波分解后的各层系数模大于和小于某阈值的系数分别进行处理，然后利用处理后的小波系数重构出消躁后的图像。
在阈值消躁中，阈值函数体现了对小波分解系数的不同处理策略以及不同估计方法，常用的阈值函数有硬阈值函数和软阈值函数。
硬阈值函数可以很好的保留图像的边缘等局部特征，但图像会出现伪吉布斯效应等视觉失真现象；软阈值处理相对较平滑，但可能会造成边缘模糊等失真现象，为此人们又提出了半软阈
值函数。
小波阈值消躁方法处理阈值的选取，另一个关键因素是阈值的具体估计。
如果阈值太小消躁后的图像仍然存在噪声；相反如果阈值太大，重要图像特征又将滤掉，弓I起偏差。
直观上将，对给定的小波系数，噪声越大，阈值就越大。
图像信号的小波消躁步骤有三步，同一维信号的消躁步骤完全相同，不同的是二维小波变换代替一维小波变换。
二维小波分析用于图像消躁的步骤如下：
步骤1:
二维图像信号的小波分解
步骤2:
对分解后的高频系数进行阈值量化。
步骤3:
二维小波重构图像信号。
二、例程分析
［例6-9］利用小波分析对给定一个二维含躁图像进行消躁处理。
xiaobo0610.m
'原始图像含躁图像
50
100
150
200
250
50100150200250
50
100
150
200
250
50100150200250
第1层重构图像
50
100
150
200
250
50100150200250
第2层重构图像
50
100
150
200
250
50100150200250
［例6-10］利用二维小波变换对给定图像进行消躁处理。
xiaobo0611.m
150
200
250
50
100150200250
第二次消躁后的图像
第一次消躁后的图像
50
100
150
200
250
50100150200250
50
100
150
200
250
50100150200250
原始图像
50
100
150
200
250
50100150200250
第一次压缩后的图像
分解后的低频和高频信息
20
40
60
6.4小波分析与图像压缩
所谓图像压缩就是去掉各种冗余，保留重要信息。
虽然图像的数据是非常巨大的，但是可以采用适当的坐标变换去除相关从而达到压缩数据的目的。
［例6-11］利用二维小波变换对给定图像进行压缩处理。
xiaobo0612.m
100
200
300
400
500
100200300400500
第二次压缩后的图像
20
40
60
80
100120
20406080100120
2040
60
第一次压缩后图像的大小：
NameSizeBytesClass
cal135x135145800doublearray
Grandtotalis18225elementsusing145800bytes
第二次压缩后图像的大小：
NameSizeBytesClass
ca275x7545000doublearray
Grandtotalis5625elementsusing45000bytes
分析：
第一次压缩，压缩比较小，约为14;
第二次压缩，压缩比较大，约为114。
视觉效果也可以。
但在对压
我们一般不仅在前两层压缩，理论上可以获得任意压缩比的压缩图像,缩比和图像质量要求较高的情况下，不如其他的编码方法。
小波分析还可以用于图像平滑、融和、增强以及边缘检测等。
升级会员 