1、小波小波变换原理与原理与应用用学院:学院:电子信息工程学院子信息工程学院专业:xxx姓名:姓名:时间:2016年3月26号1.为什么需要要什么需要要对信号信号进行行变换4原始信号有一些信息是很难获取的,为了获得更多的信息,我们需要对原始信号进行数学变换。从而获得更多的信息。例如生活中常见的心电图,在心电图的时域信号中一般很难找到这些病情,所以心脏病专家一般用记录在磁带上的时域心电图来分析心电信号,从而确定病症是否存在。2.一、FT和STFT二、小波变换三、小波变换在图像处理中的应用主要内容主要内容3.1.1 傅里叶傅里叶变换(FT)FT:IFT:通通过上述上述FT公式可以公式可以发现,信号的,
2、信号的频域是一些指数域是一些指数项的累加和,每个指数的累加和,每个指数项对应特定的特定的频率,然后在整个率,然后在整个时域整合起来。其中指数域整合起来。其中指数项可以用以下的表达式表示:可以用以下的表达式表示:即信号是由一些不同即信号是由一些不同频率的正弦率的正弦项叠加起来的,叠加起来的,如果信号中如果信号中频率率为f的分量幅度的分量幅度较大,那么大,那么这个分量就个分量就和正弦和正弦项重叠,他重叠,他们的即就比的即就比较大,大,这表明信号有一表明信号有一个个频率率为f的主要分量。的主要分量。4.信号一 cos(2*pi*10*t)+cos(2*pi*25*t)+cos(2*pi*100*t)
3、+cos(2*pi*50*t)信号二 5.对上面两个信号上面两个信号进行行FT后得到的后得到的频域域图信号一信号二由于由于这个信个信号的号的频率分率分量一直保持量一直保持不不变,我,我们将此将此类信号信号称之称之为平平稳信号信号非平非平稳信号信号6.由上面两个由上面两个频域域图可以看出傅里叶可以看出傅里叶变换只能只能给出信出信号的号的频谱分量,而无法分量,而无法给出相出相应的的频谱分量的出分量的出现时间,当我当我们想知道想知道频率分量出率分量出现的确切的确切时间时,傅里叶,傅里叶变换对于非平于非平稳信号是不合适的。而且信号是不合适的。而且现实中几乎所有的生中几乎所有的生物信号都是非平物信号都是
4、非平稳的。那么我的。那么我们应该怎怎样将将时间信息加信息加到到频率率图中去呢?中去呢?这时我我们可以考可以考虑将部分非平将部分非平稳信号信号看成平看成平稳信号。信号。7.1.2STFTSTFT:STFT只不只不过是是对乘了一个窗函数的信号做傅里叶乘了一个窗函数的信号做傅里叶变换,以此得到在某段以此得到在某段时间内的内的频率信息。率信息。根据海森堡根据海森堡测不准原理,不准原理,在在STFT中由于窗口中由于窗口长度有限,度有限,它它仅仅覆盖了信号的一部分,因此覆盖了信号的一部分,因此导致致频率分辨率率分辨率较差,即我差,即我们不能确切的知道信号中那些不能确切的知道信号中那些频率分量存在,率分量存
5、在,只知道那些只知道那些频段的分量存在。段的分量存在。8.如果我们有一个无限长的窗口,然后做傅里叶变换,会得到完美的频率分辨率,但是结果中不包含时间信息。更进一步为了获得信号的平稳性,我们需要一个宽度足够短的窗函数,窗口越短,时间分辨率越高,信号的稳定性越高,但是频率分辨率却越来越低。窄窗窄窗=高高时间分辨率,低分辨率,低频率分辨率率分辨率宽窗窗=高高频率分辨率,低率分辨率,低时间分辨率分辨率9.加窄窗之后加窄窗之后对应的的STFT,可,可见有有较好好的的时间分辨率,但分辨率,但是是频率分辨率很差。率分辨率很差。加加较宽窗之后窗之后对应的的STFT,可,可见有有较好的好的频率分辨率,率分辨率,
6、但是但是时间分辨率很分辨率很差。差。10.2.1 小波的小波的发展展历史史工程到数学工程到数学41807:Joseph FourierFT,只有频率分辨率而没有时间分辨率41909:Alfred Haar发现了Haar小波41945:GaborSTFT41980:MorletMorlet小波,并分别与20世纪70年代提出了小波变换的概念,20世纪80年代开发出了连续小波变换CWT(continuous wavelet transform)41986:Y.Meyer提出了第一个正交小波Meyer小波41988:Stephane MallatMallat快速算法(塔式分解和重构算法)11.小波的小
7、波的发展展历史史工程到数学工程到数学41988:Inrid Daubechies作为小波的创始人,揭示了小波变换和滤波器组(filter banks)之间的内在关系,使离散小波分析变成为现实4Ronald Coifman和Victor Wickerhauser等著名科学家在把小波理论引入到工程应用方面做出了极其重要贡献4在信号处理领域中,自从Inrid Daubechies完善了小波变换的数学理论和Stephane Mallat构造了小波分解和重构的快速算法后,小波变换在各个工程领域中得到了广泛的应用,典型的如语音信号处理、医学信号处理、图像信息处理等12.2.2.1 连续小波小波变换 如果函
8、数如果函数 满足以下容足以下容许性条件:性条件:则称称 为一容一容许性小波,并定性小波,并定义如下的如下的积分分变换:以上以上积分分变换为 以以 为母小波的母小波的积分分连续小波小波变换,a为尺度因子,表示与尺度因子,表示与频率相关的伸率相关的伸缩,b为时间平移因子。平移因子。13.2.2.2离散小波离散小波变换 将将a,b离散化,令离散化,令 ,可以,可以得到离散小波得到离散小波变换:其中:其中:14.2.3 几种常用小波几种常用小波(1)Haar小波A.Haar于1990年提出一种正交函数系,定义如下:15.(2)Meyer函数 Meyer小波函数 和尺度函数 都是在频域中进行定义的,是具
9、有紧支撑的正交小波。其中,为构造函数Meyer的辅助函数,且有:16.(3)其他常用小波)其他常用小波Daubechies(dbN)小波系Biorthogonal(biorNr.Nd)小波系Symlets(symN)小波系Morlet(morl)小波Coiflet(CoifN)小波系17.2.4 塔式算法塔式算法 (1)信号在小波空间的展开为:(2)小波分解算法使用多分辨析的金字塔算法:实际计算时,可以一次一次地进行小波分解,然后递推实现(J-j)次小波分解。假设一次小波分解的尺度系数和小波系数为:18.可以得到:因为:故:代入得:从而:最后得到如下尺度系数和小波系数:19.基于离散余弦变换的
10、图像压缩算法,其基本思想是在频域对信号进行分解,去除信号点之间的相关性,并找出重要系数,滤掉次要系数,以达到压缩目的。但是该方法在处理过程中不能提供时域的信息,在比较关心的时域特性的时候显得无能为力。三、小波分析在三、小波分析在图像像处理中的理中的应用用 有时,对图像分析时,需要对某个局部的细节要求有很高的分辨率,单纯的时域分析方法显然不能达到这个要求,虽然可以通过对图像进行块分解,然后对每块作用不同的阈值或掩码来达到这个要求,但分块大小相对固定,有失灵活性。在这个方面,小波分析就优越的多。由于小波分析具有固定的时频特性,可以在时频两个方向对系数进行处理,这样就可以对感兴趣的部分提供不同的压缩
11、精度。20.3.1 二维小波变换的图像压缩 二维小波变换用于图像压缩是小波分析应用中的一个重要方面。它的特点是压缩比高,压缩速度快,压缩后能保持图像的特征基本不变,并且在压缩过程中可以抗干扰。对于一幅图像(二维信号)来说,表现一副图像的最主要部分是低频部分,所以最简单的压缩方法是利用小波分解,去掉图像的高频部分而只保留低频部分。21.3.2 小波变换用于图像去噪4噪声就是妨碍人的视觉器官或系统传感器所接收图像源进行理解和分析的各种噪声。一般噪声为不可预测的随即信号。二维信号用小波分析的去噪步骤有3步:(1)二维信号的小波分解(2)对高频系数进行阈值量化(3)二维小波的重构。22.3.3 在其他在其他图像像处理的理的应用用(1)小波分析用于图像增强(2)图像的钝化和锐化(3)小波分析用于图像分解23.谢谢24.