不同小波函数对粗差识别效果的比较.doc
《不同小波函数对粗差识别效果的比较.doc》由会员分享,可在线阅读,更多相关《不同小波函数对粗差识别效果的比较.doc(9页珍藏版)》请在三一文库上搜索。
1、由国鳌薮谗文任纹不同小波函数对粗差识别效果的比较1姚丽慧,王坚,高井祥中国矿业大学环境与测绘学院 (江苏徐州221008)Email: hnylh2000摘要:本文采用含孤立态、离散态与区域态三种类型粗差的含强噪声的动态变形信号,研究 粗差识别中小波基函数的选择依据。分别从分解层数、支撑区间、正则性以及消失距四个角 度,选用不同的小波基函数进行试验,结果表明:对于该含不同类型粗差的变形信号,在所 考察的小波基函数范围内,对于区域态粗差几种小波的识别能力相当;同等条件下,对于孤 立态粗差,应选择支撑长度较小、正则性较好或消失距较高的小波进行识别;对离散状粗差 进行识别时,一般应选择支撑长度较大、
2、正则性较好或消失距较低的小波函数,对于本文的变形信号进行5层分解可以得到较好的识别效果。关键词:动态变形小波基函数粗差识别1. 引言随着变形传感器等新的动态监测技术的发展,动态监测数据也发生了根本性的变化:变 形数据源的采样频率大幅提高,数据量膨胀,蕴涵信息频段加宽。但由于传感器技术采集数 据受其本身、一些人为的因素(如技术不熟练、测量不小心)及外界的影响,常常导致在变 形信息数据中存在不同类型的粗差,如孤立态、离散态和区域态等。粗差的存在不利于我们 对测量数据的分析。目前已有一些探测粗差的方法,如文1中的Chauvenet、Grubbs、Dixon及t分布等准则;文2中的数据探测法、 选权迭
3、代法,但它们都以数据按正态分布为前提的、 计算较繁等缺点,不利于实际应用。小波分析具有带通滤波的功能,能对原始信号有效分频,从而在不同尺度上将粗差和噪 声分开,达到识别的目的。通过小波变换的多尺度分析,反映形变信号的内在特征,并分离 形变趋势项。比起传统的分析方法,小波分析具有独特的优势。文3采用MALLAT多尺度分 解算法,分析形变信号中可能的粗差点位置,在各分解层选择不同的阈值去除强噪声,最后 采用线性内插修复趋势项中的粗差,得到满意效果。但是,在实际应用中,由于通常无法区分真实的信号和污染它的噪声,而且可用于粗差 定位的小波基函数是一个集合,因此在小波粗差定位的实际应用中采用哪一种小波函
4、数才有 最好的定位效果,是一个有待解决的、同时很有实际价值的研究课题。本文将采用含孤立态,离散态与区域态三种类型粗差的含强噪声的动态变形信号,来评 估用不同性质的小波函数进行粗差定位的能力,得出粗差定位和识别中小波基函数选择的一1国家教育部博士点基金(20040290503): “动态监测数据特征提取与变形预报理论研究”,中国矿业大学 科技基金资助(2005B020): ”动态变形数据特征分析与建模新理论”1rh屋I ssi務令打览些基本依据。2. 小波基函数与标准傅立叶变换相比,小波分析中所用到的小波函数具有不惟一性,即小波函数(x)具有多样性。根据不同的标准,小波函数具有不同的类型,这些标
5、准通常有:2.1紧支撑性若函数Y ( t)在区间a , b 外恒为零,则称该函数紧支在这个区间上,具有该性质的小波称为紧支撑小波。显然支集愈窄小波的局部化能力就愈强。2.2小波的消失矩小波函数W (t)对所有的o Wl WM满足/t " (t)dt = 0则称它具有M阶消失矩。另外设Wj,k(t) j,k Z构成L2(R)的标准正交基,那么若+乂 M +1小波函数叭t)满足"(t)(R)且/ NN (t)dtv+3也称它有m阶消失矩。由定义可知小波的消失矩主要用来检测高阶导数不连续的信号,消失矩越高光滑性就越好,频域的局部化能力就越强,反映了小波对信号奇异性检测能力的强弱。
6、2.3正则性正则指时间连续的小波基函数至少是连续的,或者更好一点,是一阶或二阶连续可微的。正则序(regularity order)是连续可微的次数,因而小波基函数的正则性是小波基函数逼近的光 滑性的量度,正则性越好收敛越快。在实际处理中,不存在既有紧支撑、对称又有很好的正则性和消失矩的小波基函数。面对不同的数据和目的,应选择相应的小波基函数,只有这样才能使处理效果达到最佳。如果进行区域背景研究,可以少考虑消失矩;如果研究的是局部异常,则消失矩应考虑在内。小波函数的选择,对分析结果影响很大。根据不同的用途,应慎重选择合适的小波函数。动态变形体的形变是连续的,因此,动态变形监测数据中的粗差特征表
7、现为信号函数的奇异性。根据原信号的性质分析,选择小波基函数时,倾向于选用正则性较好的小波基函数,此外还应有较好紧支撑性和消失矩。因而选取了 dbN小波、BiorNr.Nd小波、symN小波和coifN小波分别对含粗差信号进行粗差定位。3. 实例研究强污染可认为包括服从某种分布的随机噪声与一定量的粗差。形变信号可看成非线性运3m国科技论文任线动的趋势变化。对含强噪声1000个历元的形变监测数据加入粗差,第300历元强制为2mmM立态粗差,400历元和403历元强制为1.5mm的离散态粗差,600到620加入区域态粗差 2mm 得强干扰形变数据(图 1)。图中很难看出信号的形变趋势;在600历元处
8、能看到区域态粗差的存在,而孤立点粗差和离散状粗差的位置根本无法确定。变形卓還毫米3210-110020030040050060070080090010003210-110020030040050060070080090010001-25m国科技论文任线1-2#m国科技论文任线观测历元数原始信号和含粗差信号3.1分解尺度对粗差识别的影响首先,选择具有较好正则性和支撑区间的db2小波函数对含粗差信号分别进行3层、4层、5层、6层分解,如图2。3210-12.5De com position at level 3 : s = a3 + d3 + d2 + d110.d120-2a33210-121.
9、510.502°0.20-0.2-0.40.50.5953210-11.510.50.20-0.2-0.40.40.20-0.20.20-0.2-0.410.50-0.520(a) 3层分解Decomposition at leve l 5 : s = a5 + d5 + d4 + d3 + d2 + d1ii di* n “ <* 屮 j0.-0.520-2D ecom pos ition at level 4 : s= a4 + d4 + d3 + d2 + d1、imu喇鼻4卜片m利(b)4层分解position at leve l 6 : s = a6 + d6 + d
10、5 + d4 + d3 + d2 + d1 .1-2#m国科技论文任线1-2#m国科技论文任线(c) 5层分解(d) 6层分解图2 db2小波对含粗差信号进行不同层数的分解1-2#也国m冒圖菇罡立存虫更SI WHlab fB图2 (a)、(b)、( c)、(d)中,第一层分解的高频部分di都可认为主要是由噪声构成;各种粗差在其他不同的细节部分都有所体现,但是对于图2(a)和图2(b)的近似部分 a3和a4还含有一些高频信息,不能很好地体现信号的趋势项。图2(d)中,d6没有其他的高频信息。综合考虑趋势信号识别的效果以及各层高频信号的可分性情况,最终确定5层小波分解结构,既能体现不同类型粗差的位
11、置,又能很好地体现信号的趋势项。5层分解中,d5、d4、d3、d2层的低频部分都部分的检测出了粗差的位置,d5层能清晰看出三种类型粗差的存在性 ,但时间分辨率都较低;孤立态粗差点在d2和d5层得到体现,但d2层的时间分辨率趋向更高,而且相当突出,d5层的分辨率却很低,说明孤立态粗差点在高频段更易识别;可以认为离散态粗差的频率主要集中在d4和d5分解层,但在d4分解层更易识别,而且分辨率比较高;而区域态粗差信息在d5,d4,d3层出现,说明其频率趋于低频,d4层时间分辨率更高些。针对此例,d2层识别孤立态粗差精度较高,而离散态和区域态粗差在d4层能得到较高的时间精度。3.2紧支撑性对粗差识别的影
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 不同 函数 识别 效果 比较
链接地址:https://www.31doc.com/p-12580292.html