- 1、本文档共10页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
地震波的小波去噪分析-天津大学研究生e-Learning平台
天津大学攻读博士学位研究生课程
小波分析及其应用
论文题目:地震波的小波去噪分析
研 究 生: 林澍
学 院: 建筑工程学院
专 业: 岩土工程
方 向: 地基处理
学 号: 1015205008
地震波的小波分析
天津大学建筑工程学院工程专业1015205008 林澍
摘要:小波分析在岩土工程领域主要用于地震波处理,将复杂且不规则的地震波信号通过小波变换后进行相应的处理,再重构地震波,进而用于工程问题分析中。本文运用MATLAB对地震波进行不同的小波变换处理,再对地震波信号进行去噪处理,对各小波的去噪效果进行了分析。
关键词:小波分析;地震去噪;阈值;
1 引言
小波分析是一种新兴的时频分析方法,被认为是Fourier变换的突破性进展成果,它在时域和频域上都具有良好的局部化特性,对于分析非稳定信号具有强大优势。在岩土工程领域,小波分析主要用于地震波的处理。一方面,地震波具有很强的非稳定性,以正弦余弦函数为基底的Fourier变换难以对其进行处理,需要具有多分辨率的小波进行分析;另一方面,地震信号在采集时会掺入噪声,更加剧了信号的不稳定性,需要通过小波分析对地震信号进行去噪处理,以便获取更为真实准确的地震波信号。因此,越来越多学者开始从小波分析方面进行地震波特性的研究。
2 小波去噪 小波去噪的原理
在实际的工程中, i=0,…,n-1 (1)
其中,f(i)为真实信号,e(i)为噪声信号,s(i)为含噪信号。
去噪的目的便是将s(i)中的e(i)减小至最小程度,因此,通过小波变换将信号分解后,对得到的高频系数进行去噪处理,然后重构信号,即可达到去噪目的。去噪处理一般采用阈值去噪进行。
以三层小波分解为例,首先将信号进行小波分解,分解示意图如图1所示。图中的CD1、CD2、CD3即为高频系数,通过阈值去噪方式对其进行处理,分别得到去噪后的高频系数CD1’、CD2’、CD3’,再将信号重构,便可得到去噪后的信号。
图1 小波分解和重构示意图
2.2 小波基的选择
小波分析中小波函数选择最小波基。
L2(Rn)的正交或双正交小波基是非常多的,一般要求小波基具有一定的紧支撑性、平滑性和对称性。紧支集保证有优良的空间局部特性;对称性保证不会造成信号的失真;平滑性保证频率分辨率。但是上述三点不可能同时得到满足。要求小波具有较高的光滑型,必须要求增加小波支集的长度;反之,为了有利于算法实现,支集的长度要尽量小,但这又保证不了光滑性。综合考虑,保证一定的紧支集特性、对称性和平滑性。
本次分析采用默认阈值、硬阈值、软阈值三种方法去噪。
3 地震波小波变换处理
3.1 地震波的小波分解
3.1.1 小波基选取
根据之前的小波基选取原则,结合工程处理经验,本文选取db1(Haar),db4,coif3,sym4四种小波对地震波信号进行处理。各小波的特点如表1所示。
表1 各小波的特点[2]
小波函数 Haar Daubechies Coiflets Symlets 小波形式 haar db4 coif3 sym4 正交性 有 有 有 有 双正交性 有 有 有 有 紧支撑性 有 有 有 有 对称性 对称 近似对称 近似对称 近似对称 支撑长度 1 7 17 7 滤波器长度 2 8 18 8
3.1.2 小波分解
确定小波基之后,可以通过MATLAB对地震波进行小波变换处理。本次分析选取EI_Centro地震波进行,输入30秒监测数据,再用MATLAB进行相应的分析,原始地震波波形图如图2所示。
图2 原地震波信号
通过小波分析代码,将原始地震波分为三层。选取的小波基为db1(Haar)、db4、coif3和sym4,所得的三层信号的高频和低频系数分别如图3~图6所示。
图3 db1小波分解
图4 db4小波分解
图5 coif3小波分解
图6 sym4小波分解
通过上述分解可知,低频信号的振幅与原信号更加接近,而高频信号比较集中且振幅较小。因此,通过阈值去噪便可把信号的噪声去掉,从而保留有用的信息。
3.2 地震波的小波去噪
3.2.1 小波去噪步骤
通常而言,一维信号的去噪过程可分为如下三个步骤进行:
(1)信号的小波(包)分解:选择一个小波函数并确定小波分解的层数N,然后对信号s进行N层小波(包)分解;
(2)高频系数进行阈值量化:对每一层高频系数,选择一个阈值进行量化处理;
(3)信号的重构:根据小波分解的第N层的低频系数和经过量化处理后的N层高频系数,进行信号的重构。
在去噪时,最关键的一步是第(2)步。信号的去噪效果,与阈值选择有直接关系。阈值选取过小,则去噪效果欠佳;阈值选取过大,则信号会失真。因此,应选取多
您可能关注的文档
最近下载
- 矽力杰产品规格书SY5055.pdf VIP
- 道闸项目报价清单.xlsx VIP
- 贵州省预防接种技能竞赛理论考试题库资料(含答案).pdf VIP
- AASHTO-LRFD-SI-2007(4)桥梁手册第四章中文版.doc VIP
- AASHTO LRFD Bridge Design Specifications(AASHTO LRFD桥梁设计规范).docx VIP
- 迅达常用备件识别号查询手册_Ae5.pdf VIP
- 数学学科业务工作报告.pdf VIP
- 2023年春学期人教版初中物理九年级下册教学进度表.docx VIP
- 基于AASHTO沥青混凝土路面设计方法的应用与研究.doc VIP
- 小班科学方方和圆圆教案反思.docx VIP
文档评论(0)