基于小波分析与数学形态学融合算法对地磁信号降噪处理应用.docVIP

基于小波分析与数学形态学融合算法对地磁信号降噪处理应用 　　摘要：本文拟应用MATLAB软件对容易受到外界环境干扰和由于仪器自身干扰的影响，对地磁信号进行降噪处理。通过小波分析和数学形态学融合算法，同时将该方法与数学形态学算法两者进行对比分析，并且计算出降噪后信号相对于含噪信号的均方误差的值，然后对其不同算法的降噪信号的处理结果与未受干扰的原始正常信号进行对比分析，分别求其信噪比的大小，综合考虑两种算法?Φ卮藕?噪信号处理的优势和劣势。进而得出小波分析与数学形态学算法对地磁信号的降噪处理，能最大化的减少信号的失真程度，保留其原始形态。为以后批量处理地磁观测的干扰数据提供新方法和新思路，同时也能抑制对地磁信号的干扰和提高观测数据的内在质量。 　　关键词：MATLAB软件；降噪处理；小波分析；数学形态学；均方误差；信噪比 　　中图分类号：TN911.4 文献标识码：A 文章编号：1007-9416（2018）02-0122-03 　　随着数学信号处理技术的不断提升，基于MATLAB软件对地磁信号降噪处理的方法也越来越多，截止目前为止，降噪处理的方法主要有：基于fft硬件滤波法，信号相关法，自适应滤波法，图谱识别法，有限冲激响应滤波法和小波变换法等等。对于先前已有的降噪方法，均有明显的缺点，主要表现在对降噪信号处理不彻底，或者波形容易失真等。所以开始对地磁信号研究新的降噪方法非常重要。数学形态学算法其主要应用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境[1]。该算法将一个含噪的信号分解成具有物理意义的各个部分，然后将其背景剥离，从而有效抑制含噪信号的毛刺和缝隙，保留原始信号的主要形态。应用数学形态学的算法处理含噪信号时，不容易发现本身微弱信号的尖峰和毛刺，同时由于含噪信号因其有振荡衰减的特点，所以在含噪信号特别微弱的时候，数学形态学算法很难在形态上区分是振荡衰减的信号还是信号的尖峰和毛刺。因此单纯应用数学形态学算法对地磁信号的降噪处理无法取得理想效果。小波分析对地磁信号的研究较为广泛，因其本身具有可以调整时间窗和频率窗的特点，可以精确找到某一段地磁信号的突变点，对地磁信号的降噪算法研究具有深刻意义。通过小波分析，进行门限阈值量化处理后进行小波重构，进而抑制信号的噪声。结合小波分析和数学形态学两者算法的优点，将其两者融合，来对地磁信号进行降噪处理，进而实现抑制地磁信号的干扰信息，最大化，最有效的恢复原始的地磁真实信号。 　　本文利用小波分析与数学形态学融合算法对地磁的高频干扰信息进行降噪处理，同时将该方法与数学形态学算法进行对比分析，综合分析两者的滤波效果，达到消除地磁信号因外界环境受到的高频干扰信息，两者通过定量对比分析其均方误差值和信噪比的大小，分析其优势和劣势，最大化的提取真实地磁信号的初步分析和探讨。 　　1 基于数学形态学算法的混合滤波器设计 　　1.1 数学形态学算法介绍 　　数学形态学是积分几何研究成果的基础上创立的，是基于集合论的数学分支，它提供了非常有效的非线性滤波技术[2]。数学形态学的运算主要包括腐蚀Θ，膨胀，形态开OC和形态闭CO的级联组合。对于地磁信号而言，以GM4仪器为例，记录到的值属于一维的离散信号f（j）。其主要包括地磁真实信号x（j）和含噪信号y（j），其表达式为： 　　所以对于地磁信号所有的高频干扰含噪信号而言，其主要就是平滑和抑制峰值和谷底的噪声，由此引出来开-闭OC滤波器和闭-开CO滤波器，两者任意一个滤波器均可以实现对地磁信号高频干扰的降噪处理。单个的滤波器虽然可以消除高频干扰信息，但是开-闭OC滤波器和闭-开CO滤波器因分别具有开运算的收缩性和闭运算的扩张性，会使地磁信号的滤波输出结果偏大或者偏小，进而影响其滤波效果。因此本文拟采用开-闭OC滤波器闭-开CO滤波器输出结果的平均值来组成混合滤波器，通过混合滤波器的设计，对地磁信号进行降噪处理，可以更好地消除高频干扰。得到的滤波效果会更接近真实信号，其表达式为： 　　1.3 结构元素的确定 　　设计好混合滤波器的算法后，组成混合滤波器的结构元素直接影响其地磁信号的滤波结果。常见的结构元素有余弦，直线，三角，半圆等。本文拟采用余弦结构和直线结构分别对高频干扰信息进行消除。此余弦结构的的形状由幅值S和宽度W共同决定，而直线结构的形状仅有宽度W决定。所以应用表达式（8），选取不同结构元素的形状，对于含有地磁行的高频干扰进行消除，对地磁信号进行降噪处理。分析其不同种类的结构元素，对其滤波的效果。 　　2 小波分析与数学形态学算法的融合滤波器设计 　　基于数学形态学算法的混合滤波器设计后，接着开始设计小波分析与数学形态学算法相融合的滤波器，其主要包括以下几个步骤： 　　①选择合适的小波基db和尺