基于MATLAB的小波变换的降噪原理及性能仿真

基于MATLAB的小波变换的降噪原理及性能仿真
基于MATLAB的小波变换的降噪原理及性能仿真

基于MATLAB的小波变换的降噪原理及性能仿真

按小波变换的发展过程划分,大致可以划分三个阶段:

第一阶段:孤立应用时间。主要特征是一些特殊构造的小波在某些科学研究领域的特定问题上的应用。这个时代最典型的代表工作是法国地球物理学家J.Morlet和A.Grossmann第一个把“小波”用于分析处理地质数据,引进了以他们的名字命名的时间—尺度小波,即Grossmann-Morlet小波。这个时期的另一个代表性工作是1981年J.Stromberg对A.Haar在1910年所给出的Haar(哈尔)系标准正交小波基的改进。同时,著名的计算机视觉专家D.Marr在他的“零交叉”理论中使用的可按“尺寸大小”变化的滤波算子,现在称为“墨西哥帽”的小波也是这个时期有名的工作之一,这部分工作和后来成为S.Mallat的正交小波构造理论支柱之“多尺度分析”或“多分辨分析”有密切联系。这个时期一个有趣的现象是各个领域的专家、学者和工程师所从事的领域广泛分布于科学和技术研究的许多方面。因此,这个现象从另一个侧面预示了小波分析理论研究和应用热潮的到来,说明了小波理论产生的历史必然性。

第二阶段:国家性研究热潮和统一构造时期。真正的小波热潮开始与1986年,当时法国数学家Y.Meyer成功地构造出具有一定衰减性质的光滑函数,这个函数(算子)的二进尺度伸缩和二进整倍数平移产生的函数系构成著名的2-范数函数空间的标准正交基。这项成果标志“小波分析”新时代的到来。

第三阶段:全面应用时期。从1992年开始,小波分析方法进入全面应用阶段。在前一阶段研究工作基础上,特别是数字信号和数字图像的Mallat分解和重构算法的确定,使小波分析的应用迅速波及科学研究和工程技术应用研究的几乎所有的领域。编辑部是在美国的Texas A&M 大学的国际杂志《Applied and Computation Harmonic Analysis》从1993年创刊之日起就把小波分析的理论和应用研究作为其主要内容,编辑部的三位主编C.K.Chi、R.Coifman与I.Daubechies

都在小波分析的研究和应用中有独到的贡献。时至今日,小波分析的应用范围还在不断扩大,许多科技期刊都刊载与小波分析有关的论文,各个学科领域的地区性和国际性学术会议都有设计小波分析的各种类型的论文、报告。同时,在国际互联网和其他有较大影响的网络上,与小波有关的书籍、论文、报告、软件、随时随地有可以找到并可以免费下载,甚至颇有国际影响的软件公司MathWorks在它的“科学研究和工程应用”软件MATLAB中,特意把小波分析作为其“ToolBox”的单独一个工具箱。由此可以大致了解小波分析广泛应用状况。

1.1.2从小波变换的思想来源划分

按小波变换的思想来源划分,大致可以分为两个阶段:

第一阶段:小波变换的思想来源于伸缩与平移方法。小波分析方法的提出,最早应属1910年Haar提出的规范正交基,但当时并没有出现“小波”这个词。1936年Littlewood和Paley对傅立叶级数建立了二进制频率分量分组理论,对频率按二进制进行划分,其傅立叶变换的相位变化并不影响函数的大小,这是多尺度分析思想的最早来源。1946年Gabor提出的加窗傅立叶变换(或称短时傅立叶变换)对弥补傅立叶变换的不足起到了一定的作用。后来,Galderon、Zygmund、Stem 等将L—P理论推广到高维,并建立了奇异积分算子理论;1965年Galderon发现了再生核公式,它的离散形式已接近小波展开,只是还无法得到一个正交系的结论。1981年,Sterm对Haar系数进行了改进,证明了小波函数的存在性。1982年Battle 在构造量子场论中采用了caldem再生核公式的展开形式。

第二阶段:1984年,法国地球物理学家J.Morlet在分析地震数据时提出将地震波按一个确定函数的伸缩、平移系展开,他与A.Grossman共同研究,发展了连续小波变换的几何体系。1985年,法国的大数学家Meyer首先提出了光滑小波的正交基,1986年,Meyer及其学生Memaarie提出了多尺度分析的思想。1987年Mallat将计算机视觉领域内的多尺度分析思想引入到小波分析中,提出了多分辨分析的概念,统一了在此之前的所有正交小波基的构造,并提出了相应的分解与重构快速算法。1988年,年轻的女数学家Dallbechies提出了具有紧支集的光滑正交小波基——Daubechies基,为小波的应用研究增添了催化剂。同年,Daubechies 在美国主办的小波专题讨论会上进行了10次演讲,引起了广大数学家、物理学家甚至某些企业家的重视,由此将小波的理论和实际应用推向了一个高潮。

1.2小波变换的应用领域

事实上小波分析的应用领域十分广泛,它包括:数学领域的许多学科;信号分析[1~4]、图象处理[5,6];量子力学、理论物理;军事电子对抗与武器的智能化;计算机分类与识别[7~10],音乐与语言的人工合成;医学成像与诊断;地震勘探数据处理[11];大型机械的故障诊断等方面;例如,在数学方面,它已用于数值分析、构造快速数值方法、曲线曲面构造、微分方程求解、控制论等。在信号分析方面的滤波[12~15]、去噪声[16~19]、压缩、传递等。在图象处理方面的图象压缩、分类、识别与诊断、去污等。在医学成像方面的减少B超、CT、核磁共振成像的时间等。

1.2.1小波分析在地球物理勘探中的应用

(1)地震数据压缩。将地震记录作小波变换,变换后的结果做阈值量化,去除大量接近于零的值,用一定的记录方式把结果存储起来,达到压缩的目的。当需要再利用这些地震数据时,作小波逆变换恢复原来的地震记录。

(2)油气预测。地球物理勘探中,寻求地壳物质物性参数的奇异性是非常有意义的。例如,断层会使重力异常产生的较大变化;在地壳介质的分界面处,地震波的传播会产生速度和方向的变化,这些都是地球物理信号的奇异性。判断出奇异性的大小和位置就可以对异常现象做出解释。应用中:通常是将分形几何理论和模式识别理论与小波变换的突变点原理相结合,通过确定表征奇异性的数检测地震道的奇异性,预测储层所在的位置。通过计算地震道的分维数或提取小波变换域的地震特征参数,建立关联维数或地震的特征参数与含油气的关系,利用模式识别的原理确定油气井的位置。

1.2.2小波分析用于信号和图像处理

(1)数据压缩。随着科学技术特别是计算机技术的发展以及互联网的普及,许多应用领域(如卫星监测、地震勘探、天气预报)都存在海量数据传输或存储问题,如果不对数据进行压缩,数量巨大的数据就很难存储、处理和传输。因此,伴随小波分析的诞生,数据压缩一直是小波分析的重要应用领域之一,并由此带来巨大的经济效益和社会效益。

(2)语音分析与处理。小波理论应用于语音分析与处理的主要内容包括:

清/浊音分割;基音检测与声门开启时刻定位;去噪、压缩、重建几个方面。1.2.3小波分析在其他领域的应用

从数学的角度讲,小波分析的发展,对微分方程、积分方程的数值解、统计学等学科也注入了新的活力。因此,小波分析在流体力学的模型建立和求取数值解、医学细胞识别、线性系统计算、物理学分析、工程计算[20,21]中也得到了应用。由于小波分析处于高速发展阶段,新的理论和应用领域不断涌现。

1.3小波分析应用前景

(1)瞬态信号或图像的突变点常包含有很重要的故障信息,例如,机械故障、电力系统故障、脑电图、心电图中的异常、地下目标的位置及形状等,都对应于测试信号的突变点。虽然这些问题发生的背景不同,但都可以归结到如何提取信号中突变点的位置及判定其奇异性(或光滑性)的问题。对图像来说,急剧变化的点通常对应于代表图像结构的边缘部位,也就是图像信息的主要部分。掌握了它,也就掌握了图像的基本特征,因此,小波分析在故障检测和信号的多尺度边缘特征提取方面的应用具有广泛的应用前景。

(2)神经网络与小波分析相结合,分形几何与小波分析相结合是国际上研究的热点之一。基于神经网络的智能处理技术,模糊计算、进化计算与神经网络结合的研究,没有小波理论的嵌入很难取得突破。非线性科学的研究正呼唤小波分析,也许非线性小波分析是解决非线性科学问题的理性工具。

(3)小波分析用于数据或图像的压缩,目前绝大多数是对静止图像进行研究的。面向网络的活动图像压缩,长期以来是采用离散余弦变换(DCT)加运动补偿(Mc)作为编码技术,然而,该方法存在两个主要的问题:方块效应和蚊式噪声。利用小波分析的多尺度分析不但可以克服上述问题,而且可首先得到粗尺度上图像的轮廓,然后决定是否需要传输精细的图案,以提高图像的传输速度。因此研究面对网络的地速率图像压缩的小波分析并行算法,具有较高探索性和新颖性。同时也具有较高的应用价值和广泛的应用前景。

(4)目前使用的二维及高维小波基主要是可分离的。不可分离二维及高维小波基的构造、性质应用研究,由于理论上较为复杂,这方面的成果甚少。也许向量小波及高维小波的研究能够为小波分析的应用开创一个新天地。

1.4小波分析面临的主要问题

小波分析虽然在许多应用领域已取得了一定的成果但事实上小波分析仍面临的一些问题,主要问题如下:

(1)小波理论尚不完善,除一维小波理论比较成熟以外,高维小波、向量小波的理论还远非人们所期待的那样,特别是各类小波,如正交小波、双正交小波及向量小波、二进小波、离散小波的构造和性质的研究。

(2)最优小波基选取方法的研究。虽然国内外已有一些最优基选取方法的研究但缺乏系统规范的最佳小波基选取方法,即针对不同的问题能最优地选择不同的小波基以实现最好的应用效果。我们知道不存在一种小波基能适应所有的情况,因此,小波基的优化选择始终是小波理论研究的重要内容。

(3)小波分析的应用范围虽然很广,但真正取得极佳应用效果的领域并不多,人们正在挖掘有前景的应用领域。

(4)目前小波分析软件远不如有限差分方法(FDM)、有限元方法(FEM)、边界元方法(EEM)等软件成熟和完善,更无大型系统权威的小波分析软件,作为商品的高水平小波分析软件几乎没有。

(5)小波分析在数据图像压缩方面已取得很好的成绩,人们期待利用小波能够实现高压缩比、高重现度图像的压缩,并探索在图像的边缘检测、分类与描述中的应用。

1.5小波分析与傅里叶对比

小波分析是20世纪80年代后期形成的一种新兴的数学分支,是当前应用数学和工程学科中一个迅速发展的新领域,经过近10年的探索研究,重要的学形式化体系已经建立,理论基础更加扎实。小波分析是在傅里叶分析的基础上发展起来的,但小波分析与傅里叶分析存在着极大的不同,与Fourier变换相比,小波变换是空间(时间)和频率的局部变换,因而能有效地从信号中提取信息。通过伸缩和平移等运算功能可对函数或信号进行多尺度的细化分析,解决了Fourier变换不能解决的许多困难问题。小波变换联系了应用数学、物理学、计算机科学、信号与信息处理、图像处理、地震勘探等多个学科。数学家认为,小波分析是一个新的数学分支,它是泛函分析、Fourier分析、样调分析、数值分

析的完美结晶;信号和信息处理专家认为,小波分析是时间—尺度分析和多分辨分析的一种新技术,它在信号分析、语音合成、图像识别、计算机视觉、数据压缩、地震勘探、大气与海洋波分析等方面的研究都取得了有科学意义和应用价值的成果。

1.6本章小结

在本章中,对小波进行了基本的介绍。先对小波变换的发展从两个方面进行了介绍包括小波变换的发展及小波变换的思想。接下来分别从小波分析在地球物理勘探、信号和图像处理及在其他领域的应用说明了小波变换的如今的应用领域。最后简单的介绍了小波分析的应用前景及面临的主要问题并且把小波分析和傅里叶进行了简单的对比。

第二章 小波变换的基本原理

小波变换是一种信号的时间——尺度(时间——频率)分析方法,它具有多分辨分析的特点,而且在时频两域都具有表征信号局部特征的能力,是一种窗口大小固定不变但其形状可改变,时间窗和频率窗都可以改变的时频局部化分析方法。即在低频部分具有较低的时间分辨率和较高的频率分辨率,在高频部分具有较高的时间分辨率和较低的频率分辨率,很适合于分析非平稳的信号和提取信号的局部特征,所以小波变换被誉为分析处理信号的显微镜。在处理分析信号时,小波变换具有对信号的自适应性。

小波变换的应用是与小波变换的理论研究紧密地结合在一起的。现在,它已经在科技信息产业领域取得了令人瞩目的成就。电子信息技术是六大高新技术中重要的一个领域,它的重要方面是图象和信号处理。现今,信号处理已经成为当代科学技术工作的重要部分,信号处理的目的就是:准确的分析、诊断、编码压缩和量化、快速传递或存储、精确地重构(或恢复)。从数学地角度来看,信号与图象处理可以统一看作是信号处理(图象可以看作是二维信号),在小波变换地许多分析的许多应用中,都可以归结为信号处理问题。现在,对于其性质随实践是稳定不变的信号,处理的理想工具仍然是傅立叶分析。但是在实际应用中的绝大多数信号是非稳定的,而特别适用于非稳定信号的工具就是小波分析。

2.1傅里叶变换、短时傅里叶变换到小波变换

传统的信号分析是建立在傅里叶变换的基础之上的,在众多科学领域,特别是在信号处理、图像处理、量子物理等方面,傅里叶变换是重要的应用工具之一。其定义为

若()()R L t f 2∈(()R L 2表示平方可积的实数空间,即能量有限的信号空间),则

()()?+∞∞--=dt e t f F t j ωω

(2-1)

()()?

+∞∞-=ωωπωd e F t f t j 21 (2-2)

它的基 ()t j ωexp 为一组正交基,它体现的是一种全局变换,且具有鲜明的物理意义。其中()t f 表示时域信号,()ωF 表示信号()t f 的傅里叶变换,傅里叶变换实现了时域和频域的转换,许多在时域难以分清和解决的问题在频域可以一目了然。

虽然,傅里叶变换能够将信号的时域特征和频域特征联系起来,可以在时域或者频域对信号进行分析,但是它缺乏时间和频率的定位功能。因为傅里叶变换是整个时间域内的积分,所以没有局部化分析信号的功能,即傅里叶分析只能分析信号在整个时域的频谱,无法反映信号的局部特征。而信号的时频局域性,正是非平稳信号最根本最关键的性质。为了克服傅里叶变换的这一缺点,人们寻求了一系列新的信号分析理论,其中短时傅里叶变换和小波变换就是在这样的情况下产生的。

短时傅里叶变换(Short-time Fourier Transform )是 Dennis Gabor 于1946年提出的一种时频分析方法。其定义为

()()()?+∞∞

---=dt e t g t f SF t j ωττ?, (2-3) ()()()τωττωπωd d e t g SF t f t j -=??+∞∞-+∞∞-,21

(2-4)

它的基本思想是:把信号划分成许多小的时间间隔,假定非平稳信号在分析窗函数()t g 的一个短时间间隔内是平稳(准平稳)的,并移动窗函数,从而计算出各个不同时刻的频谱。但是,短时傅里叶变换从本质上看,它是一种单一分辨率的信号分析方法,因为它使用的是一个固定的短时窗函数,这也就是它在信号分析上的缺陷。

小波分析是傅里叶分析思想方法的发展与延拓。它既继承和发展了短时傅里叶变换的局部化思想,同时又克服了窗口大小不随频率变化的缺点,是对信号进行时频分析以及处理时变非平稳信号的比较理想的工具。

2.2小波变换的基本概念

小波变换是近几十年发展起来的能同时在时域和频域内进行局部化信号分

析的数学方法,它在分析原始信号时,在时域和频域上都具有良好的局部化性质,对不同信号采用相应的时域取样步长,能够聚焦到信号的任意微小细节,是一种优于傅里叶变换和窗口傅里叶变换的信号处理方法。

2.2.1小波定义

设()()R L t f 2∈,其傅里叶变换为()ωψ。当()ωψ满足条件[22]

=ψC ()?∞∞<02

ωωωψd (2-5)

或其等价条件

()?∞∞-=0dt t ψ (2-6)

的函数()t ψ称为一个母小波函数(Mother Wavelet Function ),式(2-5)、(2-6)小波容许条件。

将母小波()t ψ经伸缩和平移后,可得到小波序列

()??

? ??-=

a t a t a τψψτ1, (2-7)

其中a 尺度因子,τ平移,R a ∈τ,,0>a 。()t a τψ,称为小波基函数。 2.2.2小波特性

小波基函数在时域和频域都具有有限或近似有限的定义域,所以经伸缩和平移后的小波基函数在时域和频域仍是局部性的,小波基函数的窗口随尺度因子的不同而伸缩,当a 增大时小波基函数的时间窗口变大,而对应的频域窗口相应变小,中心频率降低。相反,当a 减小时,小波基函数的时间窗口变小,而对应的频域窗口相应变大,中心频率升高。定义小波母函数()t ψ的窗口宽度为0t ?,窗口中心为,其相应傅里叶变换的窗口宽度为 ,窗口中心为。小波基函数时域的窗口中心为

,窗口宽度为 t ?;频域窗口中心为,窗口宽度为

。则有下式成立

=0t *0

t a +τ ,0t a t ?=? (2-8) a *00ωω= ,a

0ωω?=? (2-9)

由此,有如下结论[23]:

(1)尺度a 的倒数在一定意义上对应频率ω ,即尺度越小,对应频率越高,尺度越大,对应频率越低。如果尺度a 对应时间窗口,则小尺度信号为短时间信号,大尺度信号为长时间信号。

(2)在任意时刻 t 上,小波的时—频窗口的大小t ?和 ω?都随尺度a 的变化而变化,这正是小波变换比傅里叶变换优越的地方。

(3)在任意尺度a ,时间点t 上,窗口面积t ??ω?保持不变,即时间、尺度分辨率是相互制约的,不可能同时提高。

(4)小波基函数作为带通滤波器,其品质因数不随尺度a 变化,即

0*0

ωωω?=?=w Q =C (常数)

(2-10)

因此,它是一组频率特性等Q 的带通滤波器组。

2.2.3小波变换原理

如果()()R L t f 2∈,且=ψC ()?∞∞<02

ωωωψd 则 ()()()?

∞-=dt t t f a W a τψτ,,

(2-11) ()()()??∞∞--∞∞-=

τψττψdad t a W a C t f a ,2,1 (2-12)

其中,()??

? ??-=a t a t a τψψτ1

,是基本小波的位移与尺度伸缩。尺度因子 a ,平移τ为连续变量,所以式(2-11)称为连续小波变换(Continuous Wavelet Transform ) 。

小波变换对()t f 而言是以()t a τψ,为核函数的线性变换。

2.3本章小结

本章首先介绍了由傅里叶变换到短时傅里叶变换再到小波变换的发展过程,然后对小波基函数的特性进行了仔细讨论,了解到了小波变换与傅里叶变换不同,小波变换通过平移小波(Mother Wavelet )可获得信号的时间信息,而通过缩放小波的宽度(或者叫做尺度)可获得信号的频率特性(对母小波的缩放和平移操作是为了计算小波的系数,这些系数代表小波和局部信号之间的相互关系),连续小波()t b a ,ψ的时频窗口中心和宽度可以精确定位,且都随尺度a 的变化而伸缩。若将时频窗口综合考虑,根据公式推导可得时频窗口的面积与尺度a 无关,即时间分辨率和频率分辨率是相互制约的。最后给出了小波变换的公式()()()?∞

∞-=dt t t f a W a τψτ,,,为小波降噪原理的介绍提供了铺垫。

第三章 小波降噪原理

在实际工程问题中,我们通过实验得到的原始信号总会混杂着一定的噪声,而噪声的存在严重地干扰了信号的本质特征,不利于进一步的信号处理和分析。因此,在对原始信号进行预处理时,对噪声加以消除或减小,以便最大程度的提取原始信号中的有用信息,是非常有必要的。本文主要考虑与信号无关的白噪声的去噪问题,信号经过去噪处理后,不但信噪比得到了提高,同时信号的一些细节特征也突现出来了。

我们把去除信号中含有的噪声并恢复原始信号的过程称为信号去噪,在信号处理领域中,人们根据实际信号的特点和噪声的统计特征,基于统计估计原理,提出了各式各样的信号去噪方法。这些去噪方法中,有的是在时域中进行的,有的是在频域中进行的,这些方法的基本思想是根据噪声和信号在时域或频域上分布的不同而进行的。

现有的一般去噪方法有基于Fourier 变换的信号去噪方法,也即低通滤波方法;基于信号的自相关去噪方法;基于小波变换的信号去噪方法等。其中,在实际问题中最常用的是滤波方法和基于小波变换的信号去噪方法。

由上一章节的讨论,我们知道小波变换是在傅里叶变换的基础上发展起来的,小波变换是比傅里叶变换更为突出的优点是其具有时频局部分析功能。将小波变换的思想用到信号去噪中,也即基于小波变换的信号去噪方法也是本文所要研究的一个内容。

3.1白噪声的特性与信号去噪性能的评价标准

在实际问题中,我们所要考虑的噪声信号大多是白噪声[24~27],本文也是在白噪声的基础上来讨论信号去噪方法的。信号去噪方法各式各样,我们有必要提出一些评价信号去噪性能好坏的标准。

3.1.1白噪声的特性

假设我们在实验中获得的原始信号中含有白噪声()t σ,则()t σ有以下一些特性[28]:

(1)()t σ可以看做一个平稳的随机信号,它在各采样点处的取值()n t σ是一

个随机变量,()n t σ取值的大小与其它采样点处的随机取值无关,白噪声之间无相关性,即任意的两个白噪声()t 1σ和()t 2σ是不相关的。

(2)()t σ可以看做是能量无限且零均值的,白噪声在时域中没有衰减性,白噪声也是随机变动的,且有()()∑∑=n

n n t t 0/2n σσ。

(3)相对于某一个确定的信号而言,()t σ在时域里的表现是均匀密集的。

(4)()t σ包含着全部的频谱,即()1?=ωσ。 3.1.2信号去噪性能的评价标准

本节主要介绍文献中较为常见的三种评价标准:信噪比、信噪比增益和均方根误差。

(1)信噪比(Signal Noise Ratio ,简记为SNR )

信噪比是测量信号中的噪声量的传统方法,常被用来作为信号去噪性能好坏的评价标准。国际上信噪比的单位是分贝(dB ),信噪比通常的定义为

???? ?

?=n s p p 10log 10SNR (3-1) 其中,()∑=n s n f n p 21表示原始信号的功率,()()[]

∑-=n n n f n f n p 2?1表示原始信号中混杂的噪声的功率,()n f 表示原始信号,()n f

?表示去噪以后的估计信号。通过信噪比定义我们知道:一个混杂有白噪声的原始信号经过不同的方法去噪以后,信噪比越大就说明这种方法的去噪效果越好。

(2)信噪比增益

记信噪比增益为A ,则信噪比增益的定义为

01f f S N R S N R A =

(3-2)

其中:1f SNR 表示信号去噪以后的信噪比,0f SNR 表示信号去噪以前的信噪比,通过信噪比增益的定义我们知道:去噪后信噪比的增益越大,就说明去噪效果越好。

(3)均方误差(MSE )

均方误差定义为

()()[]∑=-=N n n f n f N MSE 12?1

(3-3)

其中()n f 表示原始信号,()n f

?表示去噪以后的估计信号。 在本文中我们用到的信号去噪性能的评价标准为信噪比。

3.2基于小波变换的信号去噪

一般地,有用信号通常表现为低频信号或是一些比较平稳的信号,而噪声信号则通常表现为高频信号。所以消噪过程主要进行以下处理:首先对原始信号进行小波分解,则噪声部分通常包含在高频系数中;然后对小波分解的高频系数以门限阈值等形式进行量化处理;最后再对信号重构即可达到消噪的目的。对信号消噪实质上是抑制信号中的无用部分,恢复信号中有用部分的过程。

3.2.1小波信号去噪问题的原理

基于小波的信号去噪问题在数学上是一个函数逼近的问题,即如何在由小波基函数伸缩和平移所张成的函数空间中,根据某一个衡量准则,寻找对真实信号的最佳逼近,以期达到将噪声从真实信号中去除的目的。小波信号去噪问题的数学描述[29]为

()()()()()()()()(){}(){}

{}???????????→===+==-==∈的函数空间映射为,为含噪信号表示原始信号表示含噪信号,

代表最优解W I T ,span W t f f I t s t f ,t t s t f "opt ",f t s t min arg J 1j k ,j k ,j opt opt

T opt λλψφσλλλλ 由小波信号去噪问题的数学描述,实际上,基于小波的信号去噪就是为了寻找从含噪信号空间到小波函数空间的最佳映射,以便得到真实信号的最佳恢复。

从信号处理的角度来看,小波去噪问题就是一个信号滤波问题,尽管在很大程度上小波去噪可以视为低通滤波,但由于小波去噪后,还能成功地保留原有真实信号的特征信息,所以从这一点来说,基于小波的信号去噪方法是优于传统的

基于变换的信号低通滤波的。由此可见,小波去噪实际上是特征提取和低通滤波的综合,其滤波过程可以用下图3-1来表示

重构信号

含噪信号

特征提取低通滤波

特征信息

Matlab小波变换函数

Matlab小波函数 Allnodes 计算树结点 appcoef 提取一维小波变换低频系数 appcoef2 提取二维小波分解低频系数 bestlevt 计算完整最佳小波包树 besttree 计算最佳(优)树 *biorfilt 双正交样条小波滤波器组 biorwavf 双正交样条小波滤波器 *centfrq 求小波中心频率 cgauwavf Complex Gaussian小波 cmorwavf coiflets小波滤波器 cwt 一维连续小波变换 dbaux Daubechies小波滤波器计算 dbwavf Daubechies小波滤波器dbwavf(W) W='dbN' N=1,2,3,...,50 ddencmp 获取默认值阈值(软或硬)熵标准 depo2ind 将深度-位置结点形式转化成索引结点形式detcoef 提取一维小波变换高频系数 detcoef2 提取二维小波分解高频系数 disp 显示文本或矩阵 drawtree 画小波包分解树(GUI) dtree 构造DTREE类 dwt 单尺度一维离散小波变换

dwt2 单尺度二维离散小波变换 dwtmode 离散小波变换拓展模式 *dyaddown 二元取样 *dyadup 二元插值 entrupd 更新小波包的熵值 fbspwavf B样条小波 gauswavf Gaussian小波 get 获取对象属性值 idwt 单尺度一维离散小波逆变换 idwt2 单尺度二维离散小波逆变换 ind2depo 将索引结点形式转化成深度—位置结点形式*intwave 积分小波数 isnode 判断结点是否存在 istnode 判断结点是否是终结点并返回排列值 iswt 一维逆SWT(Stationary Wavelet Transform)变换iswt2 二维逆SWT变换 leaves Determine terminal nodes mexihat 墨西哥帽小波 meyer Meyer小波 meyeraux Meyer小波辅助函数 morlet Morlet小波 nodease 计算上溯结点 nodedesc 计算下溯结点(子结点)

基于小波变换的语音信号去噪(详细)

测试信号处理作业 题目:基于小波变换的语音信号去噪 年级:级 班级:仪器科学与技术 学号: 姓名: 日期:2015年6月

基于小波变换的语音信号去噪 对于信号去噪方法的研究是信号处理领域一个永恒的话题。经典的信号去噪方法,如时域、频域、加窗傅立叶变换、维纳分布等各有其局限性,因此限制了它们的应用范围。小波变换是八十年代末发展起来的一种新时-频分析方法,它在时-频两域都具有良好的局部化特性;并且在信号去噪领域获得了广泛的应用。 目前已经提出的小波去噪方法主要有三种:模极大值去噪、空域相关滤波去噪以及小波阈值去噪法。阈值法具有计算量小、去噪效果好的特点,取得了广泛的应用。然而在阈值法中,阈值的选取直接关系到去噪效果的优劣。如果阈值选取过小,那么一部分噪声小波系数将不能被置零,从而在去噪后的信号中保留了部分噪声信息;如果阈值选的偏大,则会将一部分有用信号去掉,使得去噪后的信号丢失信息。 1、语音信号特性 由于语音的生成过程与发音器宫的运动过程密切相关,而且人类发音系统在产生不同语音时的生理结构并不相同,因此使得产生的语音信号是一种非平稳的随机过程(信号)。但由于人类发生器官变化速度具有一定的限度而且远小于语音信号的变化速度,可以认为人的声带、声道等特征在一定的时间内(10- 30ms)基本不变,因此假定语音信号是短时平稳的,即语音信号的某些物理特性和频谱特性在10-30ms的时间段内近似是不变的,具有相对的稳定性,这样可以运用分析平稳随机过程的方法来分析和处理语音信号。在语音增强中就是利用了语音信号短时谱的平稳性。 语音信号基本上可以分为清音和浊音两大类。清音和浊音在特性上有明显的区别,清音没有明显的时域和频域特性,看上去类似于白噪声,并具有较弱的振幅;而浊音在时域上有明显的周期性和较强的振幅,其能量大部分集中在低频段内,而且在频谱上表现出共振峰结构。在语音增强中可以利用浊音所具有的明显的周期性来区别和抑制非语音噪声,而清音由于类似于白噪声的特性,使其与宽带平稳噪声很难区分。 由于语音信号是一种非平稳、非遍历的随机过程,因此长时间时域统计特性对语音信号没有多大的意义,而短时谱的统计特性对语音信号和语音增强有着十分重要的作用。语音信号短时谱幅度统计特性的时变性,使得语音信号的分析帧在趋于无穷大时,根据中心极限定理,其短时谱的统计特性服从高斯(Gauss)分布,而在实际应用时只能在有限帧长下进行处理,因此,在有限帧时这种高斯分布的统计特性是一种近似的描述,这样就可以作为分析宽带噪声污染的带噪语音信号增强应用时的前提和假设。

基于MATLAB的(小波)图像处理

基于MATLAB的(小波)图像处理 姓名:宋富冉 学号:P1******* 院系:电子信息工程学院 专业:电子与通信工程 日期:2015年11月7日

目录 摘要 (3) 第一章初期准备 1.1软件知识储备及学习 (4) 1.2 MATLAB操作平台安装及应用 (4) 1.3操作函数功能及调试 (5) 第二章图像准备 2.1图像采集 (6) 2.2 图像选择和保存 (6) 第三章程序设计及实现 3.1 软件编程调试 (7) 3.2 实现及优化程序 (11) 第四章完成任务报告 4.1报告书写 (12) 4.2总结 (12) 附录 (13)

摘要 本报告主要阐述有关于MABLAT在图像处理方面实际应用中的 六个方面的问题,分别涉及图像的读取、图像添加噪声、利用小波 函数对图像进行分割、分割后图像的重构、图像去除噪声、将程序 处理过程中所得各种图像确定存储格式并保存到指定的磁盘及命名。最终得到预期任务的要求,完成任务。 关键词:图像读取,图像加噪,图像去噪,图像重构,图像保存

第一章初期准备 1.1软件知识储备及学习 由于本人从未学习过MATLAB这门课程及其编程语言,对其一无所知,在之前的学习过程中,比较多的是应用C语言进行一些简单的及较复杂的任务编程。因此,接到任务之日起,本人就开始学习储备有关于此方面的软件知识,并逐步学习了解它的奥妙所在。 首先,是漫无目的的到图书馆查找有关于此类的各种书籍,并上网搜索各类处理程序和文档,以期寻求到刚好符合此次作业任务要求的完整程序设计及源代码。结果是可想而知的,并没有完全吻合的程序与代码。其次,在以上的查找翻看过程中,本人接触到了很多与此任务相关相通的程序设计和处理函数的功能及应用知识,受其启发,自我总结,将实现本任务所要用到的功能函数一一搜集了起来,初步了解了本任务如何开启。 1.2 MATLAB操作平台安装及应用 通过前期的理论准备,下一步就要开始上机实际操作和仿真各个函数在实际应用中的效果。第一步,就是寻求MATLAB操作平台的安装包或安装程序,在自己的桌面上把它装起来,以便后面随时随地使用操作,也为后期更深入的学习此门语言而准备好最基本的学习工具,从而为以后完全掌握此门语言工具打下基础。第二步,就是对本平台的安装和使用,由于此平台有中英文两个版本,于是这对我本人又是一种考验,由于英语专业词汇并不完全过关,对操作菜单中多个名词词组的用意并

matlab小波变换

matlab小波变换 Matlab 1. 离散傅立叶变换的 Matlab实现 Matlab 函数 fft、fft2 和 fftn 分别可以实现一维、二维和 N 维 DFT 算法;而函数 ifft、ifft2 和 ifftn 则用来计算反 DFT 。这些函数的调用格式如下: A=fft(X,N,DIM) 其中,X 表示输入图像;N 表示采样间隔点,如果 X 小于该数值,那么Matlab 将会对 X 进行零填充,否则将进行截取,使之长度为 N ;DIM 表示要进行离散傅立叶变换。 A=fft2(X,MROWS,NCOLS) 其中,MROWS 和 NCOLS 指定对 X 进行零填充后的 X 大小。别可以实现一维、二维和 N 维 DFT A=fftn(X,SIZE) 其中,SIZE 是一个向量,它们每一个元素都将指定 X 相应维进行零填充后的长度。 函数 ifft、ifft2 和 ifftn的调用格式于对应的离散傅立叶变换函数一致。 别可以实现一维、二维和 N 维 DFT 例子:图像的二维傅立叶频谱 1. 离散傅立叶变换的 Matlab实现% 读入原始图像 I=imread('lena.bmp');函数 fft、fft2 和 fftn 分 imshow(I) % 求离散傅立叶频谱 J=fftshift(fft2(I)); figure;别可以实现一维、二维和 N 维 DFT imshow(log(abs(J)),[8,10]) 2. 离散余弦变换的 Matlab 实现 Matlab

2.1. dct2 函数 功能:二维 DCT 变换 Matlab 格式:B=dct2(A) B=dct2(A,m,n) B=dct2(A,[m,n])函数 fft、fft2 和 fftn 分 说明:B=dct2(A) 计算 A 的 DCT 变换 B ,A 与 B 的大小相同;B=dct2(A,m,n) 和 B=dct2(A,[m,n]) 通过对 A 补 0 或剪裁,使 B 的大小为 m×n。 2.2. dict2 函数 功能:DCT 反变换 格式:B=idct2(A) B=idct2(A,m,n)别可以实现一维、二维和 N 维 DFT B=idct2(A,[m,n]) 说明:B=idct2(A) 计算 A 的 DCT 反变换 B ,A 与 B 的大小相同;B=idct2(A,m,n) 和 B=idct2(A,[m,n]) 通过对 A 补 0 或剪裁,使 B 的大小为m×n。 Matlab 2.3. dctmtx函数 功能:计算 DCT 变换矩阵 格式:D=dctmtx(n) 说明:D=dctmtx(n) 返回一个n×n 的 DCT 变换矩阵,输出矩阵 D 为double 类型。 1. 离散傅立叶变换的 Matlab实现 3. 图像小波变换的 Matlab 实现函数 fft、fft2 和 fftn 分 3.1 一维小波变换的 Matlab 实现 (1) dwt 函数 Matlab

小波分析报告(去噪)

小波分析浅析 —— 李继刚 众所周知,以π2为周期的复杂的波都可以用以π2为周期的函数)(t f (模拟信号)来描述,它可以由形如)sin(n n nt A θ+的若干谐波叠加而成,因此,完全有理由认为)(t f 有如下的表现形式: ∑ ∑ ∑ ∞ =∞ =∞ =+= += += ) sin cos ()cos sin cos sin ()sin()(n n n n n n n n n n n nt b nt a nt A nt A nt A t f θθθ 为了确定上式中的系数n n b a ,,可以利用Fourier 变换,可以得到函数)(t f 的Fourier 级数,即 ??? ? ? ? ? ?? ====++=??∑--+∞ =π πππππ.,2,1,sin )(1,,1,0,cos )(1),sin cos (2)(1 0 n ntdt t f b n ntdt t f a nt b nt a a t f n n n n n 如果函数以T 为周期,则通过对t 作T w x T t ππ2,2= ?=变换,可以得到函数的Fourier 级数,即 ??? ? ? ? ? ??=?==?=?+?+=??∑--+∞ =π πππ .,2,1,sin )(2,,1,0,cos )(2),sin cos (2)(1 0 n wtdt n t f T b n wtdt n t f T a wt n b wt n a a t f n n n n n 从时域角度来理解Fourier 级数,将}sin ,{cos wt n wt n ??看作是具有频率w n ?的谐波,则时域表现的函数)(t f 可分解为无穷个谐波之和。 从频域角度来理解Fourier 级数,因为)(t f 的频域范围是[)+∞∈,0w ,所以,可将w 轴用间距w ?作离散分化,离散点w n ?处对应着频率为w n ?的谐波}sin ,{cos wt n wt n ??,这样就可将时域函数)(t f 与谐波组成1-1对应关系,即 +∞???0}sin ,cos {)(wt n b wt n a t f n n

小波去噪matlab程序

小波去噪matlab程序 ****************************************** clear clc %在噪声环境下语音信号的增强 %语音信号为读入的声音文件 %噪声为正态随机噪声 sound=wavread('c12345.wav'); count1=length(sound); noise=0.05*randn(1,count1); for i=1:count1 signal(i)=sound(i); end for i=1:count1 y(i)=signal(i)+noise(i); end %在小波基'db3'下进行一维离散小波变换 [coefs1,coefs2]=dwt(y,'db3');%[低频高频] count2=length(coefs1); count3=length(coefs2); energy1=sum((abs(coefs1)).^2); energy2=sum((abs(coefs2)).^2); energy3=energy1+energy2; for i=1:count2 recoefs1(i)=coefs1(i)/energy3; end for i=1:count3 recoefs2(i)=coefs2(i)/energy3; end %低频系数进行语音信号清浊音的判别 zhen=160; count4=fix(count2/zhen); for i=1:count4 n=160*(i-1)+1:160+160*(i-1); s=sound(n); w=hamming(160); sw=s.*w; a=aryule(sw,10); sw=filter(a,1,sw);

用matlab小波分析的实例

1 绪论 1.1概述 小波分析是近15年来发展起来的一种新的时频分析方法。其典型应用包括齿轮变速控制,起重机的非正常噪声,自动目标所顶,物理中的间断现象等。而频域分析的着眼点在于区分突发信号和稳定信号以及定量分析其能量,典型应用包括细胞膜的识别,金属表面的探伤,金融学中快变量的检测,INTERNET的流量控制等。 从以上的信号分析的典型应用可以看出,时频分析应用非常广泛,涵盖了物理学,工程技术,生物科学,经济学等众多领域,而且在很多情况下单单分析其时域或频域的性质是不够的,比如在电力监测系统中,即要监控稳定信号的成分,又要准确定位故障信号。这就需要引入新的时频分析方法,小波分析正是由于这类需求发展起来的。 在传统的傅立叶分析中,信号完全是在频域展开的,不包含任何时频的信息,这对于某些应用来说是很恰当的,因为信号的频率的信息对其是非常重要的。但其丢弃的时域信息可能对某些应用同样非常重要,所以人们对傅立叶分析进行了推广,提出了很多能表征时域和频域信息的信号分析方法,如短时傅立叶变换,Gabor变换,时频分析,小波变换等。其中短时傅立叶变换是在傅立叶分析基础上引入时域信息的最初尝试,其基本假定在于在一定的时间窗内信号是平稳的,那么通过分割时间窗,在每个时间窗内把信号展开到频域就可以获得局部的频域信息,但是它的时域区分度只能依赖于大小不变的时间窗,对某些瞬态信号来说还是粒度太大。换言之,短时傅立叶分析只能在一个分辨率上进行。所以对很多应用来说不够精确,存在很大的缺陷。 而小波分析则克服了短时傅立叶变换在单分辨率上的缺陷,具有多分辨率分析的特点,在时域和频域都有表征信号局部信息的能力,时间窗和频率窗都可以根据信号的具体形态动态调整,在一般情况下,在低频部分(信号较平稳)可以采用较低的时间分辨率,而提高频率的分辨率,在高频情况下(频率变化不大)可以用较低的频率分辨率来换取精确的时间定位。因为这些特定,小波分析可以探测正常信号中的瞬态,并展示其频率成分,被称为数学显微镜,广泛应用于各个时频分析领域。 全文介绍了小波变换的基本理论,并介绍了一些常用的小波函数,它们的主要性质包括紧支集长度、滤波器长度、对称性、消失矩等,都做了简要的说明。在不同的应用场合,各个小波函数各有利弊。 小波分析在图像处理中有非常重要的应用,包括图像压缩,图像去噪,图像融合,图像分解,图像增强等。文中给出了详细的程序范例,用MATLAB实现了基于小波变换的图像处理。 小波分析在图像处理中有非常重要的应用,包括图像压缩,图像去噪,图像融合,图像分解,图像增强等。文中给出了详细的程序范例,用MATLAB实现了基于小波变换的图像处理。 小波分析在图像处理中有非常重要的应用,包括图像压缩,图像去噪,图像融合,图像分解,图像增强等。文中给出了详细的程序范例,用MATLAB实现了基于小波变换的图像处理。 1.2 傅立叶变换与小波变换的比较 小波分析是傅立叶分析思想方法的发展与延拓。它自产生以来,就一直与傅立叶分析

基于MATLAB的小波变换在信号分析中应用的实现

基于MATLAB的小波变换在信号分析中应用的实现 院系:应用技术学院 专业:电子信息工程 姓名:李成云 指导教师单位:应用技术学院 指导教师姓名:王庆平 指导教师职称:讲师 二零一一年六月

The application of wavelet transform based on MTLAB in signal analysis Faculty:Application and Technology Institute Profession:Electronic information engeering Name:Li Chengyun Tutor’s Unit:Application and Technology Institute Tutor:Wang Qingping Tutor’s Title:Lecturer June 2011

第 I 页 目录 摘要 (1) ABSTRACT (2) 前言 (3) 第1章 绪论 (4) 1.1 本文的研究背景意义 (4) 1.2 国内外研究现状 (5) 1.3 本文的研究内容 (7) 第2章 MATLAB 简介 (8) 2.1 MATLAB 的概况 (8) 2.2 MATLAB6.1 的功能 (8) 2.3 MATLAB 的主要组成部分 (9) 2.4 MATLAB 的语言特点 (10) 第3章 基本理论 (12) 3.1 从傅里叶变换到小波变换 (12) 3.1.1 傅里叶变换 (12) 3.1.2 短时傅里叶变换 (13) 3.1.3 小波变换 (14) 3.2 连续小波变换 (15) 3.3 离散小波变换 (17) 3.4 小波包分析 (18) 3.5 多分辨率分析与M ALLAT 算法 (19) 3.5.1 多分辨率分析 (19) 3.5.2 Mallat 算法 (19) 3.6 本章小结 (20) 第4章 小波阈值法图像去噪 (21) 4.1 图像去噪 (21) 4.1.1 邻域平均法 (22) 4.1.2 中值滤波法 (24) 4.2 小波阈值去噪 (27) 4.2.1 阈值去噪原理 (28) 4.2.2 选取阈值函数 ................................................ 28 4.2.3 几种阈值选取方法 .. (29)

基于MATLAB的小波消噪仿真实现 (1)

收稿日期:2007-12-10 作者简介:史振江(1979-),男,汉,河北唐山人,学士,讲师,研究方向智能检测与控制技术。 基金项目:河北省教育厅自然科学项目(Z2006442) 基于MATLAB 的小波消噪仿真实现 史振江1) 安建龙 2) 赵玉菊1) (石家庄铁路职业技术学院1) 河北石家庄 050041 衡水学院2) 河北衡水 053000)  摘要:小波阈值消噪方法是利用小波变换技术对含噪信号进行分解和重构,通过对小波分解后的小波系数限定阈值来消除噪声的方法。分析小波消噪的算法和实现步骤,并基于MATLAB 软件平台编写仿真程序。进行光纤光栅反射信号的小波消噪仿真实验,消噪效果良好。  关键词:小波消噪 阈值 分解 重构 光纤光栅  中图分类号:TP272 文献标识码:A 文章编号:1673-1816(2008)01-0063-04 1 引言  微弱信号检测[1]是关于如何提取和测量强噪声背景下微弱信号的方法,有效的去除信号中的噪声是实现微弱信号检测的关键。小波变换[2]是一种信号的时间、频率分析方法,具有多分辨分析的特点,是时间窗和频率窗都可以改变的时频局部化分析方法,已经广泛应用于信号消噪、信号处理、图像处理、语音识别与合成等领域。小波消噪[3~5]的方法可以分为三类:模极大值法、相关法以及阈值方法。其中,小波阈值消噪方法是利用小波变换技术对含噪信号进行分解和重构,通过对小波分解后的各层系数限定阈值来消除噪声的方法,因其实现简单、计算量小,取得了广泛应用。 MATLAB 即矩阵实验室,是一种建立在向量、数组和矩阵基础上,面向科学与工程计算的高级语言,它集科学计算、自动控制、信号处理、神经网络、图像处理于一体,具有极高的编程效率[6]。其中的小波处理工具箱可以方便实现小波消噪算法,对含噪信号进行消噪处理和研究。 本文详细分析了小波消噪算法,利用MATLAB 软件编写了程序,并对光纤光栅反射谱信号进行了小波消噪仿真实验。 2 小波变换与Mallat 算法  小波变换是指,把某一被称为基本小波的函数()t ψ平移位移b 后, 在不同尺度a 下作伸缩变换,得到连续小波序列,()a b t ψ,再与待分析信号()f t 作内积: 1/2(,)()()f R t b W a b a f t dt a ψ??=∫ (1) 在实际应用中,经常将,()a b t ψ作离散化处理,令2j a =,2j b k =g ,Z k j ∈,则得到相应的离散

MATLAB小波变换指令及其功能介绍(超级有用)解读

MATLAB小波变换指令及其功能介绍 1 一维小波变换的 Matlab 实现 (1) dwt函数 功能:一维离散小波变换 格式:[cA,cD]=dwt(X,'wname') [cA,cD]=dwt(X,Lo_D,Hi_D)别可以实现一维、二维和 N 维DFT 说明:[cA,cD]=dwt(X,'wname') 使用指定的小波基函数 'wname' 对信号X 进行分解,cA、cD 分别为近似分量和细节分量; [cA,cD]=dwt(X,Lo_D,Hi_D) 使用指定的滤波器组 Lo_D、Hi_D 对信 号进行分解。 (2) idwt 函数 功能:一维离散小波反变换 格式:X=idwt(cA,cD,'wname') X=idwt(cA,cD,Lo_R,Hi_R) X=idwt(cA,cD,'wname',L)函数 fft、fft2 和 fftn 分 X=idwt(cA,cD,Lo_R,Hi_R,L) 说明:X=idwt(cA,cD,'wname') 由近似分量 cA 和细节分量 cD 经 小波反变换重构原始信号 X 。 'wname' 为所选的小波函数 X=idwt(cA,cD,Lo_R,Hi_R) 用指定的重构滤波器 Lo_R 和 Hi_R 经小波反变换重构原始信号 X 。 X=idwt(cA,cD,'wname',L) 和 X=idwt(cA,cD,Lo_R,Hi_R,L) 指定返回信号 X 中心附近的 L 个点。 2 二维小波变换的 Matlab 实现 二维小波变换的函数别可以实现一维、二维和 N 维 DFT 函数名函数功能

--------------------------------------------------- dwt2 二维离散小波变换 wavedec2 二维信号的多层小波分解 idwt2 二维离散小波反变换 waverec2 二维信号的多层小波重构 wrcoef2 由多层小波分解重构某一层的分解信号 upcoef2 由多层小波分解重构近似分量或细节分量 detcoef2 提取二维信号小波分解的细节分量 appcoef2 提取二维信号小波分解的近似分量 upwlev2 二维小波分解的单层重构 dwtpet2 二维周期小波变换 idwtper2 二维周期小波反变换 ----------------------------------------------------------- (1) wcodemat 函数 功能:对数据矩阵进行伪彩色编码函数 fft、fft2 和 fftn 分 格式:Y=wcodemat(X,NB,OPT,ABSOL) Y=wcodemat(X,NB,OPT) Y=wcodemat(X,NB) Y=wcodemat(X) 说明:Y=wcodemat(X,NB,OPT,ABSOL) 返回数据矩阵 X 的编码矩阵 Y ;NB 伪编码的最大值,即编码范围为 0~NB,缺省值 NB=16; OPT 指定了编码的方式(缺省值为 'mat'),即:别可以实现 一维、二维和 N 维 DFT OPT='row' ,按行编码 OPT='col' ,按列编码

小波分析-经典解读

时间序列-小波分析 时间序列(Time Series )是地学研究中经常遇到的问题。在时间序列研究中,时域和频域是常用的两种基本形式。其中,时域分析具有时间定位能力,但无法得到关于时间序列变化的更多信息;频域分析(如Fourier 变换)虽具有准确的频率定位功能,但仅适合平稳时间序列分析。然而,地学中许多现象(如河川径流、地震波、暴雨、洪水等)随时间的变化往往受到多种因素的综合影响,大都属于非平稳序列,它们不但具有趋势性、周期性等特征,还存在随机性、突变性以及“多时间尺度”结构,具有多层次演变规律。对于这类非平稳时间序列的研究,通常需要某一频段对应的时间信息,或某一时段的频域信息。显然,时域分析和频域分析对此均无能为力。 20世纪80年代初,由Morlet 提出的一种具有时-频多分辨功能的小波分析(Wavelet Analysis )为更好的研究时间序列问题提供了可能,它能清晰的揭示出隐藏在时间序列中的多种变化周期,充分反映系统在不同时间尺度中的变化趋势,并能对系统未来发展趋势进行定性估计。 目前,小波分析理论已在信号处理、图像压缩、模式识别、数值分析和大气科学等众多的非线性科学领域内得到了广泛的应。在时间序列研究中,小波分析主要用于时间序列的消噪和滤波,信息量系数和分形维数的计算,突变点的监测和周期成分的识别以及多时间尺度的分析等。 一、小波分析基本原理 1. 小波函数 小波分析的基本思想是用一簇小波函数系来表示或逼近某一信号或函数。因此,小波函数是小波分析的关键,它是指具有震荡性、能够迅速衰减到零的一类函数,即小波函数)R (L )t (2∈ψ且满足: ? +∞ ∞ -=0dt )t (ψ (1) 式中,)t (ψ为基小波函数,它可通过尺度的伸缩和时间轴上的平移构成一簇函数系: )a b t ( a )t (2 /1b ,a -=-ψψ 其中,0a R,b a,≠∈ (2) 式中,)t (b ,a ψ为子小波;a 为尺度因子,反映小波的周期长度;b 为平移因子,反应时间上的平移。 需要说明的是,选择合适的基小波函数是进行小波分析的前提。在实际应用研究中,应针对具体情况选择所需的基小波函数;同一信号或时间序列,若选择不同的基小波函数,所得的结果往往会有所差异,有时甚至差异很大。目前,主要是通过对比不同小波分析处理信号时所得的结果与理论结果的误差来判定基小波函数的好坏,并由此选定该类研究所需的基小波函数。 2. 小波变换 若)t (b ,a ψ是由(2)式给出的子小波,对于给定的能量有限信号)R (L )t (f 2 ∈,其连续小波变换(Continue Wavelet Transform ,简写为CWT )为: dt )a b t ( f (t)a )b ,a (W R 2 /1-f ? -=ψ (3) 式中,)b ,a (W f 为小波变换系数;f(t)为一个信号或平方可积函数;a 为伸缩尺度;b 平移参数; )a b x ( -ψ为)a b x (-ψ的复共轭函数。地学中观测到的时间序列数据大多是离散的,设函数)t k (f ?,

一个小波变换实例及matlab实现

1、 选择()t ?或?()? ω,使{}()k Z t k ?∈-为一组正交归一基; 2、 求n h 。 1,(),()n n h t t ??-= 或??()(2)/()H ω?ω?ω= 3、 由n h 求n g 。 1(1)n n n g h -=- 或()()i G e H t ωωωπ-= 4、 由n g ,()t ?构成正交小波基函数() t φ 1,()()n n t g t φ?-=∑ 或??()(/2)(/2)G φωω?ω= Haar 小波的构造 1)、选择尺度函数。 101 ()0t t ? ≤≤?=? ?其他 易知(n)t ?-关于n 为一正交归一基。 2)、求n h 1,(),()n n h t t ??- =()2t-n)t dt ??( 其中 1 1(2)220n n t t n ?+? ≤≤?-=?? ?其他 当n=0时, 1 1(2)20t t ?? 0≤≤?=?? ?其他 当n=1时,

1 11(21)20t t ?? ≤≤?-=?? ?其他 故,当n=0,n=1时 1()(2)0n n t t n ?? =0,=1 ??-=? ?其他 当n=0时, ()(2)t t n ???-1 120t ? 0≤≤?=?? ?其他 当n=1时, ()(2)t t n ???-1 1120t ? ≤≤?=?? ?其他 故 n h ()2t-n)t dt ?? (1/0n n ?=0,=1 ?=? ??其他 3)、求n g 。 11/0 (1)1/10n n n n g h n -?=??=-=-=?? ??其他 4)、求()t φ。 1,()()n n t g t φ?-=∑ =0-1,011,1()()g t g t ??-+ (2)(21)t t - =1 102 111 20t t ? ≤≤???- ≤≤?? ??? 其他

小波变换降噪分析(精)

第四章小波变换降噪分析 小波变换是一种崭新的时域 (频域信号分析工具。它的发展和思想都来自于傅里叶分析,且在保留了傅里叶分析优点的基础上,较好的解决了时间和频率分辨率的矛盾,在频域与空间域中能够同时具有良好的局部化特性,可进行局部分析。小波去噪的基本原理是根据原始信号和噪声的小波系数在不同尺度上所具有的不同性质,构造相应的规则,在小波域采用其他数学方法对含噪信号的小波系数进行处理。 4.1 小波变换理论的研究 连续小波变换 设2( ( t L R ψ∈(2( L R 表示平方可积的的空间,即能量有限的信号空间, 其傅立叶变换为( ψ ω。当( ψω满足允许条件 (Admissible Condition: 2 ( C φωωω +∞ -∞ =<∞? (4.1 时,我们称( t ψ为一个基本小波或母小波 (Mother Wavelet 。将母小波函数 ( t ψ经伸缩和平移后,就可以得到一个小波序列。对于连续情况,小波序列 为: , ( (

a b t b t a ψ-= , a b R ∈ 0a ≠ (4.2 其中, a ——伸缩因子; b ——平移因子; ——能量归一化因子。 这样对于任一信号 20 1 1( (, ( f t b f t a b dadb C a a φ ωψ∞ ∞ -∞-= ? ?,连续小波变换定义为: , , (, (, ( ( ( a b a b CWT a b f t t f t t dt

ψ∞-∞ ==? (4.3 其逆变换为: 20 11( (, ( f t b f t a b dadb C a a φ ωψ∞ ∞ -∞-= ? ? (4.4 离散小波变换 实际应用中,尤其是在计算机上实现,如在信号处理领域,必须对连续小波加以离散化。需要强调的是,这一离散化都是针对连续的尺度参数 a 和连续平移参数 b 的,而不是针对时间变量 t 的,这与其它形式的离散化不同。在连续小波中,考虑函数(4.5: , ( ( a b t b

基于小波去噪matlab程序示例

clear all clc %在噪声环境下语音信号的增强 %语音信号为读入的声音文件 %噪声为正态随机噪声 sound=wavread('c12345.wav'); count1=length(sound); noise=0.05*randn(1,count1); for i=1:count1 signal(i)=sound(i); end for i=1:count1 y(i)=signal(i)+noise(i); end %在小波基'db3'下进行一维离散小波变换 [coefs1,coefs2]=dwt(y,'db3'); %[低频高频] count2=length(coefs1); count3=length(coefs2); energy1=sum((abs(coefs1)).^2); energy2=sum((abs(coefs2)).^2); energy3=energy1+energy2; for i=1:count2 recoefs1(i)=coefs1(i)/energy3; end for i=1:count3 recoefs2(i)=coefs2(i)/energy3; end %低频系数进行语音信号清浊音的判别 zhen=160; count4=fix(count2/zhen); for i=1:count4 n=160*(i-1)+1:160+160*(i-1); s=sound(n); w=hamming(160); sw=s.*w; a=aryule(sw,10); sw=filter(a,1,sw); sw=sw/sum(sw); r=xcorr(sw,'biased'); corr=max(r); %为清音(unvoice)时,输出为1;为浊音(voice)时,输出为0 if corr>=0.8 output1(i)=0; elseif corr<=0.1

小波变换去噪基础地的知识整理

1.小波变换的概念 小波(Wavelet)这一术语,顾名思义,“小波”就是小的波形。所谓“小”是指它具有衰减性;而称之为“波”则是指它的波动性,其振幅正负相间的震荡形式。与Fourier变换相比,小波变换是时间(空间)频率的局部化分析,它通过伸缩平移运算对信号(函数)逐步进行多尺度细化,最终达到高频处时间细分,低频处频率细分,能自动适应时频信号分析的要求,从而可聚焦到信号的任意细节,解决了Fourier变换的困难问题,成为继Fourier变换以来在科学方法上的重大突破。有人把小波变换称为“数学显微镜”。 2.小波有哪几种形式?常用的有哪几种?具体用哪种,为什么? 有几种定义小波(或者小波族)的方法: 缩放滤波器:小波完全通过缩放滤波器g——一个低通有限脉冲响应(FIR)长度为2N和为1的滤波器——来定义。在双正交小波的情况,分解和重建的滤波器分别定义。 高通滤波器的分析作为低通的QMF来计算,而重建滤波器为分解的时间反转。例如Daubechies和Symlet 小波。 缩放函数:小波由时域中的小波函数 (即母小波)和缩放函数 (也称为父小波)来定义。 小波函数实际上是带通滤波器,每一级缩放将带宽减半。这产生了一个问题,如果要覆盖整个谱需要无穷多的级。缩放函数滤掉变换的最低级并保证整个谱被覆盖到。 对于有紧支撑的小波,可以视为有限长,并等价于缩放滤波器g。例如Meyer小波。 小波函数:小波只有时域表示,作为小波函数。例如墨西哥帽小波。 3.小波变换分类 小波变换分成两个大类:离散小波变换 (DWT) 和连续小波转换 (CWT)。两者的主要区别在于,连续变换在所有可能的缩放和平移上操作,而离散变换采用所有缩放和平移值的特定子集。 DWT用于信号编码而CWT用于信号分析。所以,DWT通常用于工程和计算机科学而CWT经常用于科学研究。 4.小波变换的优点 从图像处理的角度看,小波变换存在以下几个优点: (1)小波分解可以覆盖整个频域(提供了一个数学上完备的描述) (2)小波变换通过选取合适的滤波器,可以极大的减小或去除所提取得不同特征之间的相关性 (3)小波变换具有“变焦”特性,在低频段可用高频率分辨率和低时间分辨率(宽分析窗口),在高频段,可用低频率分辨率和高时间分辨率(窄分析窗口) (4)小波变换实现上有快速算法(Mallat小波分解算法) 另: 1) 低熵性变化后的熵很低; 2) 多分辨率特性边缘、尖峰、断点等;方法, 所以可以很好地刻画信号的非平稳特性 3) 去相关性域更利于去噪; 4) 选基灵活性: 由于小波变换可以灵活选择基底, 也可以根据信号特性和去噪要求选择多带小波、小波包、平移不变小波等。 小波变换的一个最大的优点是函数系很丰富, 可以有多种选择, 不同的小波系数生成的小波会有不同的效果。噪声常常表现为图像上孤立像素的灰度突变, 具有高频特性和空间不相关性。图像经小波分解后可得到低频部分和高频部分, 低频部分体现了图像的轮廓, 高频部分体现为图像的细节和混入的噪声, 因此, 对图像去噪, 只需要对其高频系数进行量化处理即可。 5.小波变换的科学意义和应用价值

小波变换图像去噪MATLAB实现

基于小波图像去噪的MATLAB 实现 一、 论文背景 数字图像处理(Digital Image Processing ,DIP)是指用计算机辅助技术对图像信号进行处理的过程。数字图像处理最早出现于 20世纪50年代,随着过去几十年来计算机、网络技术和通信的快速发展,为信号处理这个学科领域的发展奠定了基础,使得DIP 技术成为信息技术中最重要的学科分支之一。在现实生活中,DIP 应用十分广泛,医疗、艺术、军事、航天等图像处理影响着人类生活和工作的各个方面。 然而,在图像的采集、获取、编码和传输的过程中,都存在不同程度被各种噪声所“污染”的现象。如果图像被污染得比较严重,噪声会变成可见的颗粒形状,导致图像质量的严重下降。根据研究表明,当一图像信噪比(SNR)低于14.2dB 时,图像分割的误检率就高于0.5%,而参数估计的误差高于0.6%。通过一些卓有成效的噪声处理技术后,尽可能地去除图像噪声,我们在从图像中获取信息时就更容易,有利于进一步的对图像进行如特征提取、信号检测和图像压缩等处理。小波变换处理应用于图像去噪外,在其他图像处理领域都有着十分广泛的应用。本论文以小波变换作为分析工具处理图像噪声,研究数字图像的滤波去噪问题,以提高图像质量。 二、 课题原理 1.小波基本原理 在数学上,小波定义为对给定函数局部化的新领域,小波可由一个定义在有限区域的函数()x ψ来构造,()x ψ称为母小波,(mother wavelet )或者叫做基本小波。一组小波基函数,()}{,x b a ψ,可以通过缩放和平移基本小波 来生成:

())(1 ,a b x a x b a -ψ=ψ (1) 其中,a 为进行缩放的缩放参数,反映特定基函数的宽度,b 为进行平移的平移参数,指定沿x 轴平移的位置。当a=2j 和b=ia 的情况下,一维小波基函数序列定义为: ()() 1222,-ψ=ψ--x x j j j i (2) 其中,i 为平移参数,j 为缩放因子,函数f (x )以小波()x ψ为基的连续小波变换定义为函数f (x )和()x b a ,ψ的积: ( )dx a b x a x f f x W b a b a )(1)(,,,-ψ= ψ=?+∞∞- (3) 与时域函数对应,在频域上则有: ())(,ωωa e a x j b a ψ=ψ- (4) 可以看出,当|a|减小时,时域宽度减小,而频域宽度增大,而且()x b a ,ψ的窗口中心向|ω|增大方向移动。这说明连续小波的局部是变化的,在高频时分辨率高,在低频时分辨率低,这便是它优于经典傅里叶变换的地方。总体说来,小波变换具有更好的时频窗口特性。 2. 图像去噪综述 所谓噪声,就是指妨碍人的视觉或相关传感器对图像信息进行理解或分析的各种因素。通常噪声是不可预测的随机信号。由于噪声影响图像的输入、采集、处理以及输出的各个环节,尤其是图像输入、采集中的噪声必然影响图像处理全过程乃至最终结果,因此抑制噪声已成为图像处理中极其重要的一个步骤。 依据噪声对图像的影响,可将噪声分为加性噪声和乘性噪声两大类。由于乘性噪声可以通过变换当加性噪声来处理,因此我们一般重点研究加性噪声。设

小波变换的原理及matlab仿真程序

基于小波变换的信号降噪研究 2 小波分析基本理论 设Ψ(t)∈L 2( R) ( L 2( R) 表示平方可积的实数空间,即能量有限的信号空间) , 其傅立叶变换为Ψ(t)。当Ψ(t)满足条件[4,7]: 2 () R t dw w C ψψ =<∞? (1) 时,我们称Ψ(t)为一个基本小波或母小波,将母小波函数Ψ(t)经伸缩和平移后,就可以得到一个小波序列: ,()( )a b t b t a ψ -= ,,0a b R a ∈≠ (2) 其中a 为伸缩因子,b 为平移因子。 对于任意的函数f(t)∈L 2( R)的连续小波变换为: ,(,),()( )f a b R t b W a b f f t dt a ψψ-=<>= ? (3) 其逆变换为: 211()(,)()f R R t b f t W a b dadb C a a ψ ψ+-= ?? (4) 小波变换的时频窗是可以由伸缩因子a 和平移因子b 来调节的,平移因子b,可以改变窗口在相平面时间轴上的位置,而伸缩因子b 的大小不仅能影响窗口在频率轴上的位置,还能改变窗口的形状。小波变换对不同的频率在时域上的取样步长是可调节的,在低频时,小波变换的时间分辨率较低,频率分辨率较高:在高频时,小波变换的时间分辨率较高,而频率分辨率较低。使用小波变换处理信号时,首先选取适当的小波函数对信号进行分解,其次对分解出的参数进行阈值处理,选取合适的阈值进行分析,最后利用处理后的参数进行逆小波变换,对信号进行重构。 3 小波降噪的原理和方法 3.1 小波降噪原理 从信号学的角度看 ,小波去噪是一个信号滤波的问题。尽管在很大程度上小波去噪可以看成是低通滤波 ,但由于在去噪后 ,还能成功地保留信号特征 ,所以在这一点上又优于传统的低通滤波器。由此可见 ,小波去噪实际上是特征提取和低通滤波的综合 ,其流程框图如 图所示[6] : 小波分析的重要应用之一就是用于信号消噪 ,一个含噪的一维信号模型可表示为如下

MATLAB小波变换指令及其功能介绍(超级有用).

MATLAB 小波变换指令及其功能介绍 1 一维小波变换的 Matlab 实现 (1 dwt函数 功能:一维离散小波变换 格式:[cA,cD]=dwt(X,'wname' [cA,cD]=dwt(X,Lo_D,Hi_D别可以实现一维、二维和 N 维 DFT 说明:[cA,cD]=dwt(X,'wname' 使用指定的小波基函数 'wname' 对信号X 进行分解,cA 、cD 分别为近似分量和细节分量; [cA,cD]=dwt(X,Lo_D,Hi_D 使用指定的滤波器组 Lo_D、Hi_D 对信号进行分解。 (2 idwt 函数 功能:一维离散小波反变换 格式:X=idwt(cA,cD,'wname' X=idwt(cA,cD,Lo_R,Hi_R X=idwt(cA,cD,'wname',L函数 fft、fft2 和 fftn 分 X=idwt(cA,cD,Lo_R,Hi_R,L 说明:X=idwt(cA,cD,'wname' 由近似分量 cA 和细节分量 cD 经小波反变换重构原始信号 X 。 'wname' 为所选的小波函数 X=idwt(cA,cD,Lo_R,Hi_R 用指定的重构滤波器 Lo_R 和 Hi_R 经小波反变换重构原始信号 X 。

X=idwt(cA,cD,'wname',L 和 X=idwt(cA,cD,Lo_R,Hi_R,L 指定返回信号 X 中心附近的 L 个点。 2 二维小波变换的 Matlab 实现 二维小波变换的函数别可以实现一维、二维和 N 维 DFT 函数名函数功能 --------------------------------------------------- dwt2 二维离散小波变换 wavedec2 二维信号的多层小波分解 idwt2 二维离散小波反变换 waverec2 二维信号的多层小波重构 wrcoef2 由多层小波分解重构某一层的分解信号 upcoef2 由多层小波分解重构近似分量或细节分量 detcoef2 提取二维信号小波分解的细节分量 appcoef2 提取二维信号小波分解的近似分量 upwlev2 二维小波分解的单层重构 dwtpet2 二维周期小波变换 idwtper2 二维周期小波反变换 ----------------------------------------------------------- (1 wcodemat 函数 功能:对数据矩阵进行伪彩色编码函数 fft、fft2 和 fftn 分格式: Y=wcodemat(X,NB,OPT,ABSOL Y=wcodemat(X,NB,OPT Y=wcodemat(X,NB

相关文档
最新文档