基于经验模态分解的探地雷达信号去噪处理精.docx
- 文档编号:1860835
- 上传时间:2023-05-01
- 格式:DOCX
- 页数:8
- 大小:125.22KB
基于经验模态分解的探地雷达信号去噪处理精.docx
《基于经验模态分解的探地雷达信号去噪处理精.docx》由会员分享,可在线阅读,更多相关《基于经验模态分解的探地雷达信号去噪处理精.docx(8页珍藏版)》请在冰点文库上搜索。
基于经验模态分解的探地雷达信号去噪处理精
基于经验模态分解的探地雷达信号去噪处理(精)
基于经验模态分解的探地雷达信号去噪处理
杨建军刘鸿福
(太原理工大学太原030024
【摘要】探地雷达作为一种先进的地球物理探测方法,具有探测效率高、操作简单、采样迅速、无损伤探测、探测分辨率高等优点。
探地雷达的信号的去噪问题已成为一个公认的技术难题。
本文用经验模态分解的方法对探地雷达信号进行信号去噪处理,并取得了良好的效果。
【关键词】探地雷达;经验模态分解;信号去噪
1引言
探地雷达又称地质雷达,是近几年迅速发展起来的一种高分辨高效率的无损探测技术。
探地雷达通过天线向地下发射高频电磁脉冲波,电磁波在地下介质传播过程中,当遇到存在电性差异的地下目标体,如空洞和分界面时,电磁波便会发生反射,返回到地面时由接收天线所接收。
在对接收到的雷达波信号处理和分析的基础上,根据信号的波形、振幅和双程走时等参数便可推断地下目标体的空间位置、结构、电性及几何形态,从而达到对地下隐蔽目标体的探测目的。
信号处理是探地雷达技术中的研究重点之一,其目的是以高的分辨率在探地雷达显示设备上显示反射波图像,提取反射波的振幅、相位和频率等各种有用的参数,帮助解释地质结构信息。
2固有模态函数
由于大多数信号或数据不是固有模态函数,在任意时刻数据可能包含多个振荡模式,这也解释了为什么简单的Hilbert变换不能给出一个普通信号的频率内容的完整描述。
所以必须把数据分解成固有模态函数,从物理上定义一个有意义的瞬时频率的必要条件是:
函数对称于局部零均值,且有相同的极值和过零点。
据此,Huang提出了固有模态函数的定义。
一个固有模态函数是满足如下两个条件的函数:
(1在整个数据序列中,极值点的数量与过零点的数量必须相等,或最多相差不能多于一个。
(2在任一时间点上,信号的局部极大值和局部极小值定义的包络平均值为零。
第一个限定条件是非常明显的;它近似于传统的平稳高斯过程关于窄带的定义。
第二个条件是一个新的想法;它把传统的全局限定变为局部限定。
这种限定是必须的,它可去除由于波形不对称而造成的瞬时频率的波动。
采用固有模态函数(以下简称IMF这个名称是因为它代表了信号数据中的振荡模式。
IMF在按过零点定义的每一个周期中,只包括一个本征模态的振荡,没有复杂的叠加波存在。
如此定义,一个基本的IMF并不限定为窄带信号,也可以是幅度调制和频率调制的。
事实上,它可以是非平稳的。
图1是一个典型的IMF。
固有模态函数(IMF概念的提出使得用Hilbert变换定义的瞬时频率具有实际的物理意义,而提出IMF分量的EMD分解方法的出现则使瞬时频率可用于复杂的非平稳信号的分析。
图1所示为一典型的固有模态函数,具有相同数目的过零点和极值点,上下包络关于零值对称。
图1一个典型的固有模态函数(Huang
3经验模态分解
Huang认为只有对IMF分量求出的瞬时频率才有实际的物理意义,但是大多数信号不是IMF分量,任何时刻信号中可能包含不只一个IMF分量。
因此,必须把信号分解为IMF分量。
为此,Huang提出了把信号分解为IMF分量的算法—EMD,其具体步骤如下:
300350400450500
采样点幅度/v
50100150200250300350400450500
IMF0
50100150200250300350400450500
IMF0
50100150200250300350400450500
IMF050
100150200250300350400450500
图3EMD分解A-scan所得的七个IMF分量和残余分量r7(t
4基于经验模态分解的自适应去噪算法
在受干扰背景下有效地检测信号,不仅与信号的形式和干扰的性质有关,也与信号处理的方法有关,对不同类型的信号寻找最佳的处理技术一直是信号处理及检测的主要问题之一。
小波变换通过小波基的伸缩和平移,实现了信号的时频分析局部化它能够同时保留信号的时域特征和频域特征。
由于其多辨特性,在合适的尺度下,非平稳信号中的有效成分会呈现出同噪声截然不同的特性,利用信号和噪声在多尺度空间中不同的传递特性可以获得干扰背景下信号的有效检测,这种处理信号的方法,在获得信噪比增益的同时能够保持对突变信息的良好分辨,在非平稳信号的处理中有自身的优越性。
Hilbert-Huang变换是最新发展起来的处理非线性非平稳信号的时频分析方法。
Hilbert-Huang变换吸取了小波变换多分辨的优势,同时又克服了在小波变换中需要选择小
1-3IMF
IMF
c2
IMFc3
501001502002503003504004505000
IMFc4
波基的困难,因此该方法同样可以用来对非平稳信号进行滤波和去噪[29]。
由于从信号木身的尺度特征出发对信号进行分解,该方法具有良好的局部适应性,加上瞬时频率的引入,使得可以从时频两方而同时对信号进行分析,增加了处理信号的灵活性和有效性。
基于Hilbert-Huang变换的去噪算法如下:
首先用经验模态分解方法EMD(Empiricalmodecompositionmethod获得有限数目的固有模态函数IMF(Intrinsicmodefunction,然后对其中的高频分量进行阈值处理,处理后的分量叠加便得到去噪后的信号。
实验数据为比利时皇家军事学院(RoyalMilitaryAcademy所采集,原始数据的A-scan有512个采样点。
根据上述思路对实测的某一A-scan进行处理并和小波阈值法的处理结果进行了比较。
图为两种不同的算法处理后的对比结果:
图4经验模态分解和小波阈值去噪比较
表1是原始信号在小波阈值法和经验模态分解法去噪后得到的均方根误差与信噪比。
表1
两种去噪方法的信噪比(SNR和均方误差(RMSE比较通过以上定量比较可以得出结论,基于经验模态分解的去噪方法无论在信噪比或均方误差方面均优于小波阈值法(采用sym6小波去噪。
5结论
可以看出无论信噪比还是均方误差经验模态分解法均优于小波阈值法。
通过以上例子可估计器
小波阈值法经验模态分解法SNR
22.668526.7431RMSE2.1723e-0081.3589e-008
050100150200250300
350400450500-0.02
0.0
2
采样点幅度/v050100150200250300
350400450500-0.02
0.0
2
采样点幅度/v050100150200250300
350400450500
-0.020.0
2
采样点幅度/v
以得出结论基于经验模态分解的去噪方法完全可用于探地雷达的信号去噪处理,而其效果要好于传统的小波阈值方法。
参考文献[1]孔令讲,浅地层探地雷达信号处理算法的研究[D],四川成都,电子科技大学,2003.[2]D.J.Daniels,D.J.Gunton,H.F.Scott.Introductiontosubsurfaceradar.IEEEProceeding,1998,133(4,278~326.[3]粟毅黄春琳雷文太著探地雷达理论与应用[M]科学出版社2006[4]D.J.Daniels.SurfacePenetratingRadar.PublishedbytheInstitutionofElectricalEngineers,London,UnitedKingdom,1996[5]VsevolodKovalenko,AdvancedGPRDataProcessingAlgorithmsforDetectionofAntipersonnelLandmines[D],TechnologyUniversityofDelft,Netherlands,2006[6]HarkinBrunzell.SignalProcessingTechniquesforDetectionofBuriedLandmineUsingGroundPenetratingRadar[D],Sweden,1998[7]K.S.Yee.NumericalsolutionofinitialboundaryvalueproblemsinvolvingMaxwell’sequationsinitialboundaryvalueproblemsinvolvingMaxwellequationsinisotropicmedia
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 经验 分解 雷达 信号 处理