小波变换理论与方法课件

,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,小波变换理论与方法,小波变换理论与方法,主要内容,傅里叶变换,小波变换,小波变换的一些应用,主要内容傅里叶变换,一 傅里叶变换,1822,年，法国数学家傅里叶,(J.Fourier),发表的研究热传导理论的“热的力学分析”,提出“每一个周期函数都可以表示成三角函数之和”，奠定了傅里叶级数的理论基础。,1829,年，法国数学家狄利克雷（,P.G.Dirichlet,）以严密的方式给出傅里叶级数与积分存在条件的完整证明。,狄利克雷条件,(Dirichlet Conditions),（,1,）在一周期内，如果有间断点存在，则间断点的数目应是有限个；,（,2,）在一周期内，极大值和极小值的数目应是有限个；,（,3,）在一周期内，信号是绝对可积的,一 傅里叶变换 1822年，法国数学家傅里叶(,若周期信号 满足狄利克雷条件，则可展开为傅里叶级数。,傅里叶级数表达式：,直流分量：,余弦分量的幅度：,正弦分量的幅度：,基波角频率 ，为 的周期。,若周期信号 满足狄利克雷条件，则可展开为傅里叶级数。傅,1.1,连续傅里叶变换,对于函数,f(t)L,1,(R),其连续傅里叶变换为,其中,i,是虚数单位，,是频率变量。,F(),的连续傅里叶逆变换为,1.1 连续傅里叶变换对于函数f(t)L1(R),其连续傅,1.2,离散傅里叶变换,对于实数或者复数离散时间序列,f,0,f,1,F,N-1,若,满足 ，则称,为序列,f,n,离散傅里叶变换，称,为序列,f,n,逆离散傅里叶变换,1.2 离散傅里叶变换对于实数或者复数离散时间序列f0,f,平稳信号是指分布参数或者分布律随时间不发生变化的信号，也就是统计特性,(,期望与方差）不随时间变化而变化。,平稳信号是指分布参数或者分布律随时间不发生变化的信号，也就是,小波变换理论与方法课件,小波变换理论与方法课件,1.3,短时傅里叶变换,为了提取信号的局部特征,例如变形信号在某一时刻的频率、形变突发位置等，,1946,年,Gabor,提出了短时傅里叶变换，即,Gabor,变换，也称加窗傅里叶变换。,Gabor,变换的基本思想为：取时间函数 作为窗口函数，然后用 通待分析函数相乘，,是时间延迟，是窗函数,g(t),的中心，窗函数根据,进行时移，然后再进行傅里叶变换：,其中 ，窗口函数,g(t),起着时限作用，起着频限作用。该变化具有不变化宽度（由时间宽度决定）和不变的窗口面积,4g,g,1.3 短时傅里叶变换为了提取信号的局部特征,例如变形信号在,短时傅里叶变换示意图,短时傅里叶变换示意图,傅里叶变换图,短时傅里叶变换图,傅里叶变换图短时傅里叶变换图,二 小波变换,小波变换由法国科学家,MORLET,于,1980,年在进行地震数据分析时提出，是强有力的时频分析,(,处理,),工具，是在克服傅立叶变换缺点的基础上发展而来的。已成功应用于很多领域，如信号处理、图像处理、模式识别等。,小波变换的一个重要性质是它在时域和频域均具有很好的局部化特征，它能够提供目标信号各个频率子段的频率信息。这种信息对于信号分类是非常有用的。,小波变换一个信号为一个小波系数，这样一个信号可由小波系数来刻画。,二 小波变换 小波变换由法国科学家MORLET于1980年,小波变换是一个平方可积分函数,f(t),与一个在时频域上均具有良好局部性质的小波函数,(t),的内积：,式中，,表示内积，,a0,为尺度因子，,b,为位移因子，,*,表示复数共轭，,a,b,(t),称为小波基函数,。,2.1,连续小波变换,小波函数时间频率窗,(t),称为母小波，,(t),必须满足容许性条件：,小波变换是一个平方可积分函数f(t)与一个在时频域上均具,部分小波波形,部分小波波形,小波分类的标准,支撑长度：即当时间或频率趋向于无穷大时，它们从一个有限值收敛到,0,，长度越小，对奇异点的区分效果越好。,对称性：对称性越好，越能保证信号不失真（不产生畸变），越能提高信号的重构精度。,正则性：,它在对信号或图像的重构获得较好的平滑效果作用上是非常有用的。,小波分类的标准支撑长度：即当时间或频率趋向于无穷大时，它们从,小波运算的基本步骤：,选择一个小波函数，并将这个小波与要分析的信号起始点对齐；,计算在这一时刻要分析的信号与小波函数的逼近程度，即计算小波变换系数,C,，,C,越大，就意味着此刻信号与所选择的小波函数波形越相近，如图所示。,小波运算的基本步骤：,将小波函数沿时间轴向右移动一个单位时间，然后重复步骤,(1),、,(2),求出此时的小波变换系数,C,，直到覆盖完整个信号长度，如图所示；,将所选择的小波函数尺度伸缩一个单位，然后重复步骤,(1),、,(2),、,(3),，如图所示；,对所有的尺度伸缩重复步骤,(1),、,(2),、,(3),、,(4),。,将小波函数沿时间轴向右移动一个单位时间，然后重复步骤(1)、,连续小波变换实例,连续小波变换实例,2.2,离散小波变换,在实际应用中，需要对尺度因子,a,和位移因子,b,进行离散化处理，可以取：，,m,n,为整数，,a,0,为大于,1,的常数，,b,0,为大于,0,的常数，,a,和,b,的选取与小波,(t),的具体形式有关。离散小波函数表示为：,相应的离散小波变换可以表示为：,当,a,0,=2,，,b,0,=1,时，离散小波变换称为二进离散小波变换，这样便于分析，并且适合于在计算机上进行高效的运算。,2.2 离散小波变换,2.2.1,一阶滤波：近似与细节,在小波分析中，近似值是大的缩放因子计算的系数，表示信号的低频分量，而细节值是小的缩放因子计算的系数，表示信号的高频分量。实际应用中，信号的低频分量往往是最重要的，而高频分量只起一个修饰的作用，如 声音。,2.2.1 一阶滤波：近似与细节在小波分析中，近似值是大的缩,小波变换理论与方法课件,2.2.2,多尺度分解,对信号的高频分量不再分解，而将信号的低频部分继续分解,.,实际中，分解的级数取决于要分析的信号数据特征及用户的具体需要，例如长度为,N,的信号，最多能分成,log,2,N,层。在实际中，可以选择合适的分解层数。下图为三层多尺度分解树结构，原始信号,S,的多尺度分解为：,S=cA3+Cd3+cD2+Cd1,。,2.2.2多尺度分解对信号的高频分量不再分解，而将信号的低频,小波变换理论与方法课件,2.2.3,小波包分解,小波分析是将信号分解为近似与细节两部分，近似部分又可以分解成第二层近似与细节，可以这样重复下去。对于一个,N,层分解来说，有,N+1,个分解信号的途径。而小波包分析的细节与近似部分一样，也可以分解，对于,N,层分解，它产生,2,N,个不同的途径。,2.2.3小波包分解小波分析是将信号分解为近似与细节两部分，,三 小波变换的一些应用,3.1,小波包去噪,加噪信号数学模型为,f(t)=s(t)+n(t),，,s(t),是原信号，,n(t),是随机白噪声，满足,En(t)=0,和,Dn(t)=,2,。设,(t),为小波函数，,n(t),的小波包变换为,W,n,(j,t)=n(t),j,(t)=,n(t),的小波包系数的期望和方差分别为：,E(|W,n,(j,t)|,2,)=0,D(|W,n,(j,t)|,2,)=,三 小波变换的一些应用3.1小波包去噪加噪信号数学模型为f(,3.1.1,小波包去噪步骤,选择小波基并确定最佳分解的层次,对信号 进行小波包分解,;,对步骤,(1),获得的小波包树,选择一定的嫡标准,计算最优树,;,估计阈值,并应用该阈值对最优树的小波包系数进行阈值量化,;,将经量化处理的小波包系数,重构回原始信号。,小波包阈值消噪有两个关键点：,1,、,如何估计阈值,；,2,如何利用阈值量化小波包系数。,3.1.1小波包去噪步骤 选择小波基并确定最佳分解的层次,对,熵的确定,熵：用来确定最优树的标准，熵值越小，对应的小波包基越好。,1,）香农熵：约定,0log(0)=0,，则香农熵定义为：,2,）,P,范数熵：若,P,1,，在,l,p,范数意义上定义,E(s)=,则：,E(s)=,3),对数能量熵,E(s,i,)=,，,0log(0)=0,则有,E(s)=,4),阈值熵：,E(s)=,式中，,是阈值，且,0.,熵的确定熵：用来确定最优树的标准，熵值越小，对应的小波包基越,阈值选择准则,(1),基于无偏似然估计原理的,Rigrsure,规则；,W,为一向量，其元素为小波系数的平方，并按由小到大的顺序排列，,W=w,1,w,2,w,n,，且,w,1,w,2,w,n,，再设一向量,R,，其元素为：,r,i,=n-2i-(n-i)w+/n (i=1,2,.,n),以,R,元素中的最小值,r,b,为风险值，由,r,b,的下标变量,b,求出对应的,w,b,，则,阈值,T,1,为：,T,1,=,（,2,）通用阈值,T,1,(sqtwolog,准则,),T,2,=,阈值选择准则(1)基于无偏似然估计原理的Rigrsure规则,（,3,）启发式的,stein,无偏风险阈值,T,3,（,Heursure,）准则,设,为,n,个小波系数的平方和，令,=,，,=,，则,T,3,=,(4),基于极大极小原理的,Minimax,方法,该准则采用的也是一种固定阈值，它产生一个最小均方误差的极值。具体的阈值选取规则为：,T,4,=,（3）启发式的stein无偏风险阈值T3（Heursure）,阈值量化函数的选取,阈值量化是应用所估计的阈值,T,，对小波系数进行的处理。目前，阈值量化函数主要采用两种方法。,一种是硬阈值法，当小波系数大于该阈值时，保留原值，否则置零，其公式为：,另一种是软阈值法，当小波包系数大于该阈值时，向着减小系数幅值的方向作一个收缩,，否则置零，其公式为：,式中，,sgn,（）为符号函数。,阈值量化函数的选取阈值量化是应用所估计的阈值T，对小波系数进,阈值准则,heursure,sqtwolog,rigrsure,mininmax,SNR,11.0062,28.7143,11.0062,21.9542,rmse,1.7976,0.7416,1.7976,1.0398,阈值准则heursuresqtwologrigrsuremi,阈值量化函数,硬阈值法,软阈值法,SNR,13.9391,28.7143,RMSE,1.5524,0.7416,阈值量化函数硬阈值法软阈值法SNR13.939128.714,系统性干扰信号探测,系统性干扰信号探测,阈值准则,噪声消除和系统干扰处理,小波包阈值消噪,SNR,43.4135,28.7143,RMSE,0.3556,0.7416,阈值准则噪声消除和系统干扰处理小波包阈值消噪SNR43.41,3.2,小波时频图,3.2 小波时频图,可将多分辨分析与连续小波变换结合起来分析信号的特征,可将多分辨分析与连续小波变换结合起来分析信号的特征,3.3,识别信号发展趋势,3.3 识别信号发展趋势,3.4,无参回归估计,随机设计模式,3.4 无参回归估计随机设计模式,固定设计模式,固定设计模式,谢谢聆听，请各位批评指正,谢谢聆听，请各位批