版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、基于小波变换的信号降噪研究及matlab仿真程序2 小波分析基本理论设(t)L 2( R) ( L 2( R) 表示平方可积的实数空间,即能量有限的信号空间) , 其傅立叶变换为(t)。当(t)满足条件4,7: (1) 时,我们称(t)为一个基本小波或母小波,将母小波函数(t)经伸缩和平移后,就可以得到一个小波序列: (2)其中a为伸缩因子,b为平移因子。对于任意的函数f(t)L 2( R)的连续小波变换为:(3)其逆变换为: (4)小波变换的时频窗是可以由伸缩因子a和平移因子b来调节的,平移因子b,可以改变窗口在相平面时间轴上的位置,而伸缩因子b的大小不仅能影响窗口在频率轴上的位置,还能改变
2、窗口的形状。小波变换对不同的频率在时域上的取样步长是可调节的,在低频时,小波变换的时间分辨率较低,频率分辨率较高:在高频时,小波变换的时间分辨率较高,而频率分辨率较低。使用小波变换处理信号时,首先选取适当的小波函数对信号进行分解,其次对分解出的参数进行阈值处理,选取合适的阈值进行分析,最后利用处理后的参数进行逆小波变换,对信号进行重构。3 小波降噪的原理和方法3.1 小波降噪原理从信号学的角度看 ,小波去噪是一个信号滤波的问题。尽管在很大程度上小波去噪可以看成是低通滤波 ,但由于在去噪后 ,还能成功地保留信号特征 ,所以在这一点上又优于传统的低通滤波器。由此可见 ,小波去噪实际上是特征提取和低
3、通滤波的综合 ,其流程框图如图所示6:带噪信号特征提取低通滤波特征信号重建信号小波分析的重要应用之一就是用于信号消噪 ,一个含噪的一维信号模型可表示为如下形式: k=0.1.n-1其中 ,f( k)为有用信号,s(k)为含噪声信号,e(k)为噪声,为噪声系数的标准偏差。假设e(k)为高斯白噪声,通常情况下有用信号表现为低频部分或是一些比较平稳的信号,而噪声信号则表现为高频的信号,下面对 s(k)信号进行如图结构的小波分解,则噪声部分通常包含在Cd1、Cd2、Cd3中,只要对 Cd1,Cd2,Cd3作相应的小波系数处理,然后对信号进行重构即可以达到消噪的目的。SCa1Cd1Ca1Cd1Ca1Cd
4、13.2 降噪方法一般来说, 一维信号的降噪过程可以分为 3个步骤进行5,6:1)一维信号的小波分解,选择一个小波并确定一个小波分解的层次N,然后对信号进行N层小波分解计算。2) 小波分解高频系数的阈值量化,对第1层到第N层的每一层高频系数, 选择一个阈值进行软阈值量化处理 3) 一维小波的重构。根据小波分解的第 N层的低频系数和经过量化处理后的第1层到第N 层的高频系数,进行一维信号的小波重构。在这 3个步骤中,最核心的就是如何选取阈值并对阈值进行量化,在某种程度上它关系到信号降噪的质量在小波变换中,对各层系数所需的阈值一般根据原始信号的信号噪声比来选取,也即通过小波各层分解系数的标准差来求
5、取,在得到信号噪声强度后,可以确定各层的阈值。这里着重讨论了信号在两种不同小波恢复后信号质量的不同和对信号中的信号与噪声进行分离。4仿真实验本文采用Mtalab本身程序提供的noissin信号函数及初设原始信号f(x)为例进行Matlab分析1,3,其中:e = noissin + 0.5*randn(size(e1);首先对noissin函数上叠加上随机噪声信号得到e,分别对比采用db10小波和sym8小波对信号e进行5层分解,并且细节系数选用minimaxi阈值模式和尺度噪声(db10)以及选用sure阈值模式和尺度噪声(sym8)。在进行噪声消除后,还对原信号进行进一步分析,好 EN R
6、2007a)将原始信号和噪声信号分离开来,仿真结果如图所示:图1图2图3图1-1为原始信号图形,1-2为叠加随机噪声后的图形,而1-3和1-4为利用db10和sym8小波默认阈值降噪后的信号图形。从图1-3和1-4可以看出利用db10和sym8小波降噪后的信号基本上恢复了原始信号,去噪效果明显。但是滤波后的信号与原始信号也有不同,从图中可以很直观地看到采用阈值消噪后信号特征值较少无法准确还原原始信号 这是由于为降噪过程中所用的分析小波和细节系数的阈值不恰当所致,如需要更好的恢复信号,还可以采用其它种类小波对其进行分析,通过选取不同的阈值,分析结果,得到一个合适的阈值。从图2和图3中看出,在经过
7、用db10对信号进行5层分解,然后分别对分解的第5层到第1层的低频系数和高频系数进行重构。可以得出其主要基波函数和高频噪声函数的图形,其中小分波分解的细节信号是有白噪声分解得到的,而正弦信号可以在图2中的近似信号a5得到。因为在这一层的影响已经可以忽略了,所以获得的信号就是初始信号的波形,从而把淹没在噪声中的有用信号有效地分离出来。5 总结小波变换对平稳信号的去噪声,要比传统的滤波去噪声得到的效果好用小波变换进行信号降噪处理, 既降低了噪声同时又提高了信噪比,这说明小波降噪方法是切实可行的方案, 但是由于小波函数很多,采用不同的小波进行分解, 得到的结果可能相差很大, 而变换前并不能预知哪一种
8、小波降噪效果更好,需反复试验比较才能得到良好的效果,这也是小波变换的困难之处之一。另外信号降噪过程中阀值的选取是十分重要的。本文利用两个小波( sym8 ,db 10 )以及将信号中的信噪分离开来,更加直观可行,通过分别进行信号降噪处理对所得结果与原始信号进行比较可以得出Sym8小波以及默认阈值处理后的重构信号与原始信号最为接近,与分离的结果相同。 小波分析是一种信号的视频分析方法,它具有多分辨率分析的特点 ,很适合探测正常信号中夹带的瞬态反常现象并展示其成分,有效区分信号中的突变部分和噪声。通过MATLAB编制程序进行给定信号的噪声抑制和非平稳信号的噪声消除实验表明:基于小波分析的消噪方法是
9、一种提取有用信号、展示噪声和突变信号的优越方法 ,具有广阔的实用价值。在这个越来月信息化的社会中,基于小波分析的应用前景必将越来越广泛。参考文献1徐明远,邵玉斌.MATALAB仿真在通信与电子工程中的应用M.西安:西安电子科技大学出版社,2010.2张志涌,杨祖樱等编著.MATLAB教程(R2006a-R2007a)M.北京:北京航空航天出版社,2006.3张德丰.详解MATLAB数字信号处理M北京:电子工业出版社,2010.4杨建国.小波分析及其工程应用M北京:机械工业出版社,2005.5冯毅,王香华.小波变换降噪处理及其MATLAB实现J.数字采集与处理,2006,,21(12):37-3
10、9.6禹海兰,李天云.基于小波理论的噪声信号分析J.东北电力学院学报.1997.17(3):36-40.7潘泉,张磊,孟晋丽,张洪才著,小波滤波方法及应用M.北京:清华大学出版社,2005.附仿真源码如下:N=1000;t=1:1000;f=sin(0.03*t);load noissin;e1=noissin;init=2055615866; randn('seed',init); e = e1 + 0.5*randn(size(e1);subplot(2,2,1);plot(t,f);xlabel('1 样本序列');ylabel('原始信号幅值
11、39;);grid ;subplot(2,2,2); plot(e) ;xlabel('2 测试样本序列' ) ; ylabel('含有已加噪声的信号幅值') ; grid ;s1=wden(e,'minimaxi','s','one',5,'db12'); subplot(2,2,3); plot(s1); xlabel('3 db10降噪后信号' ) ; ylabel ('db10小波降噪后的信号幅值'); grid; s2=wden(e,'heursur
12、e','s','one',5,'sym8'); subplot(2,2,4); plot(s2);xlabel('4 sym降噪后信号'); ylabel('sym8小波降噪后的信号幅值'); grid; figure; subplot(6,1,1);plot(e);ylabel('e');C,L=wavedec(e,5,'db10');for i=1:5 a=wrcoef('a',C,L,'db10',6-i); subplot(6,1,i+1); plot(a); ylabel('
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- (2024版)餐饮业线上线下融合营销合作协议
- 2024年二手公寓买卖合同(济南)
- 2024年合伙经营房地产中介服务合同
- 2024年农村商业银行贷款合同
- 2024年修订版:建筑模板供应与施工合同
- 药品流通企业静脉输液管理方案
- 2024年反担保责任书
- 2024年临时用电保险服务合同
- 2024年全球工程项目管理实践合同
- 二手摩托车买卖协议书
- 生物安全手册
- 气管切开术Tracheotomy课件
- 船用火灾报警控制器使用说明书20230801v11
- 冀教版九年级上册英语课件 Unit4 写作能力提升练
- 盆底康复治疗记录
- 齐鲁工业大学2022年上期末药物分离工程期末考试复习资料
- 《药疹(Drug Eruption)》PPT课件
- GB∕T 41426-2022 牙科学 一体式手柄牙线
- 第二章 生涯发展理论PPT课件
- 建筑门窗检测试题(共10页)
- 彩盒产品不良缺陷的识别
评论
0/150
提交评论