小波变换在信号处理中的应用1ppt课件

1、 小波分析是近小波分析是近15年来发展起来的一种新的时频年来发展起来的一种新的时频分析方法，我们可以先粗略地区分一下时域分析和分析方法，我们可以先粗略地区分一下时域分析和频域分析。频域分析。时域分析的基本目标：时域分析的基本目标：- 边缘检测和分割；边缘检测和分割；- 将短时的物理现象作为一个瞬态过程分析。将短时的物理现象作为一个瞬态过程分析。频域分析的基本目标：频域分析的基本目标：区分突发信号和稳定信号以及定量分析其能量。区分突发信号和稳定信号以及定量分析其能量。一、从傅里叶变换到小波变换一、从傅里叶变换到小波变换一、从傅里叶变换到小波变换一、从傅里叶变换到小波变换（1 1傅立叶变换的定义傅

2、立叶变换的定义1. 1. 连续傅立叶变换对连续傅立叶变换对 离散傅立叶变换对离散傅立叶变换对 :1:2j tj tFT Fjf t edtIFTf tFjed 210210:0,1,.,11:0,1,.,1knNjNnnnknNjNnkDFT X kFff ekNIDFTfX k enNN2. 2. 傅立叶变换的实质傅立叶变换的实质傅里叶变换的实质是：把傅里叶变换的实质是：把f(t)这个波形分解成许多不同频率这个波形分解成许多不同频率的正弦波的叠加和。这样我们就可以将对原函数的正弦波的叠加和。这样我们就可以将对原函数f(t)的研究的研究转化为对其权系数，及傅里叶变换转化为对其权系数，及傅里叶变

3、换F()的研究。从傅里叶的研究。从傅里叶变换中可以看出，这些标准基是由正弦波及高次谐波组成变换中可以看出，这些标准基是由正弦波及高次谐波组成的，因此它在频域内是局部化的。的，因此它在频域内是局部化的。3. 3. 傅立叶变换的局限性傅立叶变换的局限性 由左图我们看不出任何频域的性质，但从右图由左图我们看不出任何频域的性质，但从右图中我们可以明显看出该信号的频率成分，也可以明中我们可以明显看出该信号的频率成分，也可以明显的看出信号的频率特性。显的看出信号的频率特性。 虽然傅里叶变换能够将信号的时域特征和频域虽然傅里叶变换能够将信号的时域特征和频域特征联系起来，能分别从信号的时域和频域观察，特征联系

4、起来，能分别从信号的时域和频域观察，但不能把两者有机的结合起来。但不能把两者有机的结合起来。 在实际信号处理过程中，尤其是对非平稳信号在实际信号处理过程中，尤其是对非平稳信号的处理中，信号在任一时刻附近的频域特征都很重的处理中，信号在任一时刻附近的频域特征都很重要。要。（2短时傅立叶变换短时傅立叶变换 基本思想：把非稳态信号看成一系列短时平基本思想：把非稳态信号看成一系列短时平稳信号的叠加，这个过程是通过加时间窗来实现稳信号的叠加，这个过程是通过加时间窗来实现的。一般选用能量集中在低频处的实的偶函数作的。一般选用能量集中在低频处的实的偶函数作为窗函数，通过平移窗函数来实现时间域的局部为窗函数，

5、通过平移窗函数来实现时间域的局部化性质。其表达式为：化性质。其表达式为： *,j tRSf t gtedt 其中其中“”表示复共轭，表示复共轭，g(t)是有紧支集的函数，是有紧支集的函数，f(t)是被分析的信号，在这个变换中，是被分析的信号，在这个变换中， 起着频起着频限的作用，限的作用，g(t)起着时限的作用。随着时间起着时限的作用。随着时间 的的变化，变化，g(t)所确定的所确定的“时间窗时间窗在在t轴上移动，使轴上移动，使f(t)“逐渐逐渐” 进行分析。进行分析。j te g(t)往往被称之为窗口函数，往往被称之为窗口函数， 大致反映了大致反映了f(t)在在 时刻时刻频率处频率处“信号成

6、分信号成分的相对含量。这样信号在窗函数上的的相对含量。这样信号在窗函数上的展开就可以表示为展开就可以表示为 在这一区域内的状在这一区域内的状态，并把这一区域称为窗口，态，并把这一区域称为窗口， 和和 分别称为窗口的时宽和分别称为窗口的时宽和频宽，表示了时频分析中的分辨率，窗宽越小则分辨率就频宽，表示了时频分析中的分辨率，窗宽越小则分辨率就越高。很显然越高。很显然 和和 都非常小，以便有更好的时频分析效果，都非常小，以便有更好的时频分析效果，但但 和和 相互制约的。相互制约的。 （2短时傅立叶变换短时傅立叶变换 ,S , 、 （3小波变换小波变换 小波分析优于傅里叶变换的地方是，它在时域小波分析

7、优于傅里叶变换的地方是，它在时域和频域同时具有良好的局部化性质。和频域同时具有良好的局部化性质。 小波变换提出了变化的时间窗。当需要精确的小波变换提出了变化的时间窗。当需要精确的低频信息时，采用长的时间窗，频率分辨率高，当低频信息时，采用长的时间窗，频率分辨率高，当需要精确的高频信息时，采用短的时间窗，时间分需要精确的高频信息时，采用短的时间窗，时间分辨率高。辨率高。 由此可知，小波变换采用的不是时间由此可知，小波变换采用的不是时间-频率域，频率域，而是时间尺度域。尺度越大，采用越大的时间窗，而是时间尺度域。尺度越大，采用越大的时间窗，尺度越小，采用越短的时间窗，即尺度与频率成反尺度越小，采用

8、越短的时间窗，即尺度与频率成反比。比。（3小波变换小波变换(4) 小波的时间和频率特性小波的时间和频率特性 运用小波基，可以提取信号中的运用小波基，可以提取信号中的“指定时间指定时间和和“指定频率指定频率的变化。的变化。时间：提取信号中时间：提取信号中“指定时间指定时间”（时间（时间A A或时间或时间B B的变化。顾名思的变化。顾名思义，小波在某时间发生的小的波动。义，小波在某时间发生的小的波动。频率：提取信号中时间频率：提取信号中时间A A的比较慢速变化，称较低频率成分；而提取的比较慢速变化，称较低频率成分；而提取信号中时间信号中时间B B的比较快速变化，称较高频率成分。的比较快速变化，称较

9、高频率成分。 时间A时间B(5) 小波的小波的3 个特点个特点l小波变换，既具有频率分析的性质，又能表示发生的小波变换，既具有频率分析的性质，又能表示发生的时间。有利于分析确定时间发生的现象。（傅里叶变时间。有利于分析确定时间发生的现象。（傅里叶变换只具有频率分析的性质）换只具有频率分析的性质）l小波变换的多分辨度的变换，有利于各分辨度不同特小波变换的多分辨度的变换，有利于各分辨度不同特征的提取图象压缩，边缘抽取，噪声过滤等）征的提取图象压缩，边缘抽取，噪声过滤等）l小波变换比快速小波变换比快速FourierFourier变换还要快一个数量级。信变换还要快一个数量级。信号长度为号长度为M M时

10、，时， Fourier Fourier变换左和小波变换右变换左和小波变换右计算复杂性分别如下公式：计算复杂性分别如下公式： MOMMOwf,log2(6) 小波基表示发生的时间和频率小波基表示发生的时间和频率FourierFourier变换的基上小波变换基中）变换的基上小波变换基中）和时间采样基下）和时间采样基下）傅里叶变换傅里叶变换(Fourier)(Fourier)基基小波基小波基时间采样基时间采样基二、连续小波变换二、连续小波变换)()(2RLt d2)( )(t)( )(t设函数设函数，如果满足：，如果满足：则称则称为一个基本小波和小波母函数，式中为一个基本小波和小波母函数，式中为函数

11、为函数的傅立叶变换，上式也可称为可容性条件。的傅立叶变换，上式也可称为可容性条件。 1. 1. 连续小波变换连续小波变换)()(21,abtatbaRb0 Ra)(tba,ab)(,tba令：令：， 称为基本小波或母小波称为基本小波或母小波(Mother Wavelet) 依赖于依赖于生成的连续小波。式中生成的连续小波。式中为尺度因子，改变连续小波的形状；为尺度因子，改变连续小波的形状；为位移因子，改变连续小波的位移。连续小波为位移因子，改变连续小波的位移。连续小波在时域空间和频域空间上都具有局部性，其作用等同于在时域空间和频域空间上都具有局部性，其作用等同于短时傅立叶变换中的窗函数。短时傅立

12、叶变换中的窗函数。 二、连续小波变换二、连续小波变换因此函数因此函数f(t)的小波变换为：的小波变换为：RbafdtabttfafbaW)()(,),(21,)(t)(t式中式中为函数为函数的复共轭，由可容性条件得：的复共轭，由可容性条件得：0)(dtt)(baWf，的逆变换为：的逆变换为： 2,)()(1)(adadbtbaWctfbaRfR，式中：式中： dC2)( 像傅立叶分析一样，小波分析就是把一个信号分解为将像傅立叶分析一样，小波分析就是把一个信号分解为将母小波经过缩放和平移之后的一系列小波，因此小波是小波母小波经过缩放和平移之后的一系列小波，因此小波是小波变换的基函数。小波变换可以

13、理解为用经过缩放和平移的一变换的基函数。小波变换可以理解为用经过缩放和平移的一系列小波函数代替傅立叶变换的正弦波和余弦波进行傅立叶系列小波函数代替傅立叶变换的正弦波和余弦波进行傅立叶变换的结果。变换的结果。 图图4表示了正弦波和小波的区别，由此可以看出，正弦表示了正弦波和小波的区别，由此可以看出，正弦波从负无穷一直延续到正无穷，正弦波是平滑而且是可预测波从负无穷一直延续到正无穷，正弦波是平滑而且是可预测的，的， 而小波是一类在有限区间内快速衰减到而小波是一类在有限区间内快速衰减到0的函数，其平的函数，其平均值为均值为0， 小波趋于不规则、不对称。小波趋于不规则、不对称。 二、连续小波变换二、连

14、续小波变换二、连续小波变换二、连续小波变换(a)(b)二、连续小波变换二、连续小波变换傅立叶变换过程傅立叶变换过程18 基本小波函数基本小波函数()的缩放和平移操作含义如下：的缩放和平移操作含义如下： (1) 缩放。简单地讲，缩放。简单地讲， 缩放就是压缩或伸展基本小波，缩放就是压缩或伸展基本小波， 缩缩放系数越小，放系数越小， 则小波越窄，如图则小波越窄，如图6所示。所示。 图6 小波的缩放操作 OOOf (t)f (t)f (t)tttf (t)(t)；scale1f (t)(2t)；scale0.5f (t)(4t)；scale0.25小波变换过程小波变换过程19 (2) 平移。简单地讲

15、，平移就是小波的延迟或超前。在数学平移。简单地讲，平移就是小波的延迟或超前。在数学上，上， 函数函数f(t)延迟延迟k的表达式为的表达式为f(t-k)，如图，如图7所示。所示。 图7 小波的平移操作(a) 小波函数(t)； (b) 位移后的小波函数(t-k) Ot(t)Ot(t k)(a)(b)20 CWT计算主要有如下五个步骤：计算主要有如下五个步骤： 第一步：第一步： 取一个小波，取一个小波， 将其与原始信号的开始一节进行比将其与原始信号的开始一节进行比较。较。 第二步：第二步： 计算数值计算数值C， C表示小波与所取一节信号的相表示小波与所取一节信号的相似程度，计算结果取决于所选小波的形

16、状，似程度，计算结果取决于所选小波的形状， 如图如图8所示。所示。 第三步：向右移动小波，重复第一步和第二步，直至覆盖整第三步：向右移动小波，重复第一步和第二步，直至覆盖整个信号，如图个信号，如图9所示。所示。 第四步：第四步： 伸展小波，伸展小波， 重复第一步至第三步，重复第一步至第三步， 如图如图10所示。所示。 图8 计算系数值C 原始信号小波信号C0.0102二、连续小波变换二、连续小波变换图9 计算平移后系数值C 原 始 信 号小 波 信 号二、连续小波变换二、连续小波变换图10 计算尺度后系数值C 原 始 信 号小 波 信 号C 0.2247二、连续小波变换二、连续小波变换 第五步

17、：对于所有缩放，重复第一步至第四步。第五步：对于所有缩放，重复第一步至第四步。 小波的缩放因子与信号频率之间的关系是：缩放因子小波的缩放因子与信号频率之间的关系是：缩放因子scale越小，表示小波越窄，度量的是信号的细节变化，表越小，表示小波越窄，度量的是信号的细节变化，表示信号频率越高；缩放因子示信号频率越高；缩放因子scale越大，越大， 表示小波越宽，度表示小波越宽，度量的是信号的粗糙程度，表示信号频率越低。量的是信号的粗糙程度，表示信号频率越低。 二、连续小波变换二、连续小波变换二、连续小波变换二、连续小波变换结论：结论：尺度因子尺度因子a a越小，越小， 的波形变窄，的波形变窄， 的

18、的频谱向高频端扩展；频谱向高频端扩展；a a越大，越大， 波形变宽，波形变宽， 的的频谱频谱 向低频端扩展，从而实现过了向低频端扩展，从而实现过了时间频率窗的自适应调节。时间频率窗的自适应调节。连续小波变换的实质就是以基函数连续小波变换的实质就是以基函数 的形式把的形式把信号信号f(t)f(t)分解为不同频带的子信号，实现信号在不分解为不同频带的子信号，实现信号在不同频带、不同时刻的合理分离，也可以视为一个滤同频带、不同时刻的合理分离，也可以视为一个滤波器。波器。 , a bt , a b ,a bt ,a b , a bt一维连续小波变换一维连续小波变换Matlab实现实现lCOEFS=cw

19、t(S,SCALES,wname)lCOEFS=cwt(S,SCALES,wname,plot)lCOEFS=cwt(S,SCALES,wname,PLOTMODE)lCOEFS=cwt(S,SCALES,wname,PLOTMODE,XLIM) 在每个可能的缩放因子和平移参数下计算小波系数，在每个可能的缩放因子和平移参数下计算小波系数，其计算量相当大，其计算量相当大， 将产生惊人的数据量，而且有许多数据将产生惊人的数据量，而且有许多数据是无用的。如果缩放因子和平移参数都选择为是无用的。如果缩放因子和平移参数都选择为2jj0且为且为整数的倍数，整数的倍数， 即只选择部分缩放因子和平移参数来进行

20、即只选择部分缩放因子和平移参数来进行计算，计算， 就会使分析的数据量大大减少。使用这样的缩放因就会使分析的数据量大大减少。使用这样的缩放因子和平移参数的小波变换称为双尺度小波变换子和平移参数的小波变换称为双尺度小波变换Dyadic Wavelet Transform），它是离散小波变换），它是离散小波变换Discrete Wavelet Transform， DWT的一种形式。通常离散小波的一种形式。通常离散小波变换就是指双尺度小波变换。变换就是指双尺度小波变换。 三、一维离散小波变换与重构三、一维离散小波变换与重构小波变换就是将小波变换就是将 “ 原始信号原始信号 s ” 变换变换 成成 “

21、 小波小波 系数系数 w ” ， w=wa , wd （ 近似系数近似系数wa与细节系数与细节系数wd )则原始信号则原始信号s可分解成小波近似可分解成小波近似a与小波细节与小波细节d之和。之和。 s = a+d小波系数小波系数 w = wa , wd 的分量，乘以基函数，形成小波分的分量，乘以基函数，形成小波分解：解：小波近似系数小波近似系数wa 基函数基函数A=近似分解近似分解 a -平均平均小波细节系数小波细节系数wd 基函数基函数D=细节分解细节分解 d-变化变化 三、一维离散小波变换与重构三、一维离散小波变换与重构 小波基小波基D小波基小波基A A原始信号原始信号小波系数小波系数wd

22、小波系数小波系数wa正变换：原始信号在小波基上，获得正变换：原始信号在小波基上，获得 “小波系数小波系数分量分量反变换：一切反变换：一切“小波分解小波分解” 合成原始信号合成原始信号 例如：例如： 小波分解小波分解 a=小波系数小波系数 wa 小波基小波基A三、一维离散小波变换与重构三、一维离散小波变换与重构离散小波变换公式离散小波变换公式正变换正变换反变换反变换 其中：其中： 是小波基函数是小波基函数l信号信号 s 有有M个样本，个样本，J 级小波变换：级小波变换： nDnAnDwnAwndnansJjnDnswnAnswwwwwMnjJJjjjdJJaJjiJjjdJJadJdJaJ,1

23、,., 1111小波分解小波系数三、一维离散小波变换与重构三、一维离散小波变换与重构 执行离散小波变换的有效方法是使用滤波器，执行离散小波变换的有效方法是使用滤波器， 该方法该方法是是Mallat于于1988年提出的，称为年提出的，称为Mallat算法。这种方法实算法。这种方法实际上是一种信号分解的方法，际上是一种信号分解的方法， 在数字信号处理中常称为双在数字信号处理中常称为双通道子带编码。通道子带编码。 用滤波器执行离散小波变换的概念如图用滤波器执行离散小波变换的概念如图11所示。所示。S表表示原始的输入信号，示原始的输入信号， 通过两个互补的滤波器组，通过两个互补的滤波器组， 其中一其中

24、一个滤波器为低通滤波器，通过该滤波器可得到信号的近似个滤波器为低通滤波器，通过该滤波器可得到信号的近似值值A（Approximations），另一个为高通滤波器，），另一个为高通滤波器， 通过通过该滤波器可得到信号的细节值该滤波器可得到信号的细节值DDetail）。）。 三、一维离散小波变换三、一维离散小波变换图11 小波分解示意图SAD滤波器组低通高通三、一维离散小波变换三、一维离散小波变换 在小波分析中，近似值是大的缩放因子计算的系数，在小波分析中，近似值是大的缩放因子计算的系数，表示信号的低频分量，而细节值是小的缩放因子计算的表示信号的低频分量，而细节值是小的缩放因子计算的系数，表示信号

25、的高频分量。实际应用中，信号的低频系数，表示信号的高频分量。实际应用中，信号的低频分量往往是最重要的，而高频分量只起一个修饰的作用。分量往往是最重要的，而高频分量只起一个修饰的作用。如同一个人的声音一样，如同一个人的声音一样， 把高频分量去掉后，听起来声把高频分量去掉后，听起来声音会发生改变，但还能听出说的是什么内容，但如果把音会发生改变，但还能听出说的是什么内容，但如果把低频分量删除后，就会什么内容也听不出来了。低频分量删除后，就会什么内容也听不出来了。 三、一维离散小波变换三、一维离散小波变换 由图由图11可以看出离散小波变换可以表示成由低通滤波器可以看出离散小波变换可以表示成由低通滤波器

26、和高通滤波器组成的一棵树。原始信号经过一对互补的滤波和高通滤波器组成的一棵树。原始信号经过一对互补的滤波器组进行的分解称为一级分解，信号的分解过程也可以不断器组进行的分解称为一级分解，信号的分解过程也可以不断进行下去，也就是说可以进行多级分解。如果对信号的高频进行下去，也就是说可以进行多级分解。如果对信号的高频分量不再分解，而对低频分量进行连续分解，就可以得到信分量不再分解，而对低频分量进行连续分解，就可以得到信号不同分辨率下的低频分量，号不同分辨率下的低频分量， 这也称为信号的多分辨率分这也称为信号的多分辨率分析。如此进行下去，析。如此进行下去， 就会形成图就会形成图12所示的一棵比较大的分

27、所示的一棵比较大的分解树，解树， 称其为信号的小波分解树称其为信号的小波分解树Wavelet Decomposition Tree）。实际中，）。实际中， 分解的级数取决于要分析的信号数据特分解的级数取决于要分析的信号数据特征及用户的具体需要。征及用户的具体需要。 三、一维离散小波变换三、一维离散小波变换35图图12 多级信号分解示意图多级信号分解示意图（a） 信号分解；信号分解； (b) 小波分数；小波分数； （c小波分解树小波分解树 cA3cD3cA2cD2SLo_DHi_DA1D1Lo_DHi_DA2D2Lo_DHi_DA3D3Lo_D：低通滤波器；Hi_D：高通滤波器(a)ScA1cD

28、1(b)(c)ScA1cD1cA2cD2cA3cD3 对于一个信号，如采用图对于一个信号，如采用图11所示的方法，理论上产生所示的方法，理论上产生的数据量将是原始数据的两倍。于是，根据奈奎斯特的数据量将是原始数据的两倍。于是，根据奈奎斯特Nyquist采样定理，采样定理， 可用下采样的方法来减少数据量，可用下采样的方法来减少数据量，即在每个通道内高通和低通通道每两个样本数据取一即在每个通道内高通和低通通道每两个样本数据取一个，个， 便可得到离散小波变换的系数便可得到离散小波变换的系数Coefficient），）， 分别分别用用cA和和cD表示，如图表示，如图13所示。图中表示下采样。所示。图中

29、表示下采样。 三、一维离散小波变换三、一维离散小波变换图图13 小波分解下采样示意图小波分解下采样示意图 SDA1000个 采 样 点1000个 采 样 点1000个 采 样 点ScDcA1000个 采 样 点约 500个 DWT系 数约 500个 DWT系 数三、一维离散小波变换三、一维离散小波变换 在在Matlab中，离散小波变换分解算法主要使中，离散小波变换分解算法主要使用如下几个常用命令：用如下几个常用命令： dwt 用于信号的单层分解用于信号的单层分解 wavedec 用于信号的多层分解用于信号的多层分解 wmaxlev 在多层分解前求最大的分解层数在多层分解前求最大的分解层数三、一

30、维离散小波变换三、一维离散小波变换 将信号的小波分解的分量进行处理后，一般还要根据需要把信号恢复出来，也就是利用信号的小波分解的系 数 还 原 出 原 始 信 号 ， 这 一 过 程 称 为 小 波 重 构Wavelet Reconstruction或叫做小波合成Wavelet Synthesis）。这一合成过程的数学运算叫做逆离散小波变换Inverse Discrete Wavelet Transform， IDWT）。 三、一维离散小波重构三、一维离散小波重构图图14 小波重构算法示意图小波重构算法示意图 SHLHL三、一维离散小波变换与重构三、一维离散小波变换与重构 1重构近似信号与细节

31、信号重构近似信号与细节信号 由图由图14可知，由小波分解的近似系数和细节系数可以重可知，由小波分解的近似系数和细节系数可以重构出原始信号。同样，可由近似系数和细节系数分别重构构出原始信号。同样，可由近似系数和细节系数分别重构出信号的近似值或细节值，这时只要近似系数或细节系数出信号的近似值或细节值，这时只要近似系数或细节系数置为零即可。置为零即可。 图图15是对第一层近似信号或细节信号进行重构的示意是对第一层近似信号或细节信号进行重构的示意图。图。 三、一维离散小波变换与重构三、一维离散小波变换与重构图图15 重构近似和细节信号示意重构近似和细节信号示意（a） 重构近似信号；重构近似信号； (b

32、) 重构细节信号重构细节信号 A1HL1000个 样 点0约 500个 0cA1约 500个 近 似 分 量(a)D1HL1000个 样 点(b)约 500个 0约 500个 近 似 分 量0cD1三、一维离散小波变换与重构三、一维离散小波变换与重构 2多层重构多层重构 在图在图15中，重构出信号的近似值中，重构出信号的近似值A1与细节值与细节值D1之后，之后，则原信号可用则原信号可用A1D1S重构出来。对应于信号的多层小波重构出来。对应于信号的多层小波分解，小波的多层重构如图分解，小波的多层重构如图16所示。由图所示。由图16可见重构过程为：可见重构过程为：A3D3A2；A2D2A1；A1+

33、D1S。 信号重构中，滤波器的选择非常重要，关系到能否重构信号重构中，滤波器的选择非常重要，关系到能否重构出满意的原始信号。低通分解滤波器出满意的原始信号。低通分解滤波器L和高通分解滤波和高通分解滤波器器H及重构滤波器组及重构滤波器组L和和H）构成一个系统，）构成一个系统， 这个系这个系统称为正交镜像滤波器统称为正交镜像滤波器Quadrature Mirror Filters， QMF系统，系统， 如图如图17所示。所示。 三、一维离散小波变换与重构三、一维离散小波变换与重构图16 多层小波重构示意图A3D3A2D2SA1D1三、一维离散小波变换与重构三、一维离散小波变换与重构图图17 多层小

34、波分解和重构示意图多层小波分解和重构示意图 S1000HL500250250DWT小波系数S1000LIDWTHHLLH三、一维离散小波变换与重构三、一维离散小波变换与重构用于离散小波重构的命令主要有如下几个：用于离散小波重构的命令主要有如下几个： idwt 用于单层小波重构用于单层小波重构 waverec 用于多层小波重构原始信号，要求输入参数用于多层小波重构原始信号，要求输入参数 同小波分解得到结果的格式一致同小波分解得到结果的格式一致 wrcoef 用于重构小波系数至某一层次，要求输入参用于重构小波系数至某一层次，要求输入参 数同小波分解得到结果的格式一致数同小波分解得到结果的格式一致

35、upcoef 用于重构小波系数至上一层次，要求输入参数同小波分用于重构小波系数至上一层次，要求输入参数同小波分 解得到结果的格式一致解得到结果的格式一致用于得到某一层次的小波系数的命令主要有以下几个：用于得到某一层次的小波系数的命令主要有以下几个： detcoef 求得某一层次的细节系数求得某一层次的细节系数 appcoef 求得某一层次的近似系数求得某一层次的近似系数 upwlev 重构组织小波系数的排列形式重构组织小波系数的排列形式三、一维离散小波变换与重构三、一维离散小波变换与重构 二维离散小波变换是一维离散小波变换的推广，二维离散小波变换是一维离散小波变换的推广， 其实其实质上是将二维

36、信号在不同尺度上的分解，质上是将二维信号在不同尺度上的分解， 得到原始信号的得到原始信号的近似值和细节值。由于信号是二维的，因此分解也是二维的。近似值和细节值。由于信号是二维的，因此分解也是二维的。分解的结果为：分解的结果为： 近似分量近似分量cA、 水平细节分量水平细节分量cH、 垂直垂直细节分量细节分量cV和对角细节分量和对角细节分量cD。同样也可以利用二维小波。同样也可以利用二维小波分解的结果在不同尺度上重构信号。二维小波分解和重构过分解的结果在不同尺度上重构信号。二维小波分解和重构过程如图程如图18所示。所示。 四、二维离散小波变换与重构四、二维离散小波变换与重构48图图18 二维小波

37、分解和重构过程示意图二维小波分解和重构过程示意图（a） 二维二维DWT； (b) 二维二维IDWT (b)Lo_R21Lo_R12Hi_R12行列列cAj1cHj1Hi_R21Lo_R12Hi_R12行列列cVj1cDj1cAjwkeepLo_D21Lo_D12Hi_D12行列列cAj1cHj1Hi_D21Lo_D12Hi_D12行列列cVj1cDj1cAj(a)五、几种常用小波五、几种常用小波1. Haar小波小波2. Daubechies小波小波3. Symlet小波小波4. 双正交小波双正交小波biorNr.Nd）5. Coiflet小波小波6. Morlet小波小波7. Mexico草

38、帽小波草帽小波8. Meyer小波小波具有对称性的小波不产生相位畸变，在图像处理中非常有用。具有对称性的小波不产生相位畸变，在图像处理中非常有用。具有好的正则性的小波，易于获得光滑的重构曲线和图像。具有好的正则性的小波，易于获得光滑的重构曲线和图像。小波函数和尺度函数如果存在消失矩，在压缩时有用。小波函数和尺度函数如果存在消失矩，在压缩时有用。六、举例六、举例1. GUI用法简介用法简介2. 基于小波变换的图像压缩和去噪基于小波变换的图像压缩和去噪dwt和和dwt2：用于单尺度一维和二维信号的离散小：用于单尺度一维和二维信号的离散小波变换；波变换；idwt和和idwt2：实现单尺度一维和二维信

39、号的离散：实现单尺度一维和二维信号的离散小波反变换；小波反变换；wavedec和和wavedec2：分别用于多尺度一维和二维：分别用于多尺度一维和二维信号的小波分解，即信号的小波分解，即 为多分辨率分析函数；为多分辨率分析函数；appcoef和和appcoef2：分别用于从多尺度一维和二维：分别用于从多尺度一维和二维小波分解的分解结构小波分解的分解结构 中提取信号的低频系数；中提取信号的低频系数；detcoef和和detcoef2：分别用于从多尺度一维和二维：分别用于从多尺度一维和二维小波分解的分解结构小波分解的分解结构 中提取信号的高频系数。中提取信号的高频系数。压缩或去噪原理系数值百分比T

40、T3T3信号去噪的基本步骤信号去噪的基本步骤信号去噪的信号去噪的Matlab实现实现 小波分析在图象处理中的应用小波分析在图象处理中的应用 图象是二维信号，其小波变换相当于二次一图象是二维信号，其小波变换相当于二次一维信号的小波变换：。维信号的小波变换：。 （1第一次一维信号的小波变换相当于图象第一次一维信号的小波变换相当于图象的行变换。的行变换。 （2第二次一维信号的小波变换相当于图象第二次一维信号的小波变换相当于图象的列变换。的列变换。 小波变换用于图象压缩有良好的效果，已形小波变换用于图象压缩有良好的效果，已形成图象压缩的标准如成图象压缩的标准如JPEG2000。 小波变换用于图象特征抽

41、取第1级斜线细节第1级水平细节第1级垂直细节水平细节近似图象垂直细节斜线细节 第1级 L1斜线细节第1级 L1水平细节第1级 L1垂直细节第2级 L2细节近似图象第3级 L3小波系数分级方块表示法 第 3 级 L3分辨率第 2 级 L2分辨率第 1 级 L1分辨率 小波系数分级树形表示法小波变换用于信号压缩的基本步骤小波变换用于信号压缩的基本步骤小波变换用于图象压缩小波变换用于图象压缩l采用小波进行压缩。作采用小波进行压缩。作“小波变换小波变换后，统计后，统计特性有改善，消除行和列之间的相关关系。特性有改善，消除行和列之间的相关关系。l有损压缩：根据视觉原理，不同分辨率小波系有损压缩：根据视觉原理，不同分辨率小波系数进行比特分配。然后转换到一维作熵编码，数进行比特分配。然后转