小波变换理论与方法..教学提纲

1、第一页，共41页。主要主要(zhyo)(zhyo)内容内容第二页，共41页。一一 傅里叶变换傅里叶变换(binhun)(binhun) 1822年，法国数学家傅里叶(J.Fourier)发表的研究热传导理论(lln)的“热的力学分析”,提出“每一个周期函数都可以表示成三角函数之和” ，奠定了傅里叶级数的理论(lln)基础。 1829年，法国数学家狄利克雷（P.G.Dirichlet）以严密的方式(fngsh)给出傅里叶级数与积分存在条件的完整证明。狄利克雷条件(Dirichlet Conditions)（1 1 ）在一周期内，如果有间断点存在，则间断点的数目应是有限个；）在一周期内，如果有间断

2、点存在，则间断点的数目应是有限个；（2 2）在一周期内，极大值和极小值的数目应是有限个；）在一周期内，极大值和极小值的数目应是有限个；（3 3）在一周期内，信号是绝对可积的）在一周期内，信号是绝对可积的第三页，共41页。0111( )cos()sin()nnnf taantbnt若周期信号 满足狄利克雷条件，则可展开为傅里叶级数。( )f t傅里叶级数(j sh)表达式：010011( )tTtaf t dtT直流分量：余弦分量的幅度：010112( )cos()tTntaf tnt dtT正弦分量的幅度：010112( )sin()tTntbf tnt dtT基波角频率 ， 为 的周期。(

3、)f t1T112T第四页，共41页。()( )itFeftd t 1.1 连续(linx)傅里叶变换对于函数(hnsh)f(t)L1(R),其连续傅里叶变换为其中(qzhng)1()() |( ) |LRfXftd x i是虚数单位，是频率变量。F()的连续傅里叶逆变换为1( )()2itfteFd 第五页，共41页。210()()(0 ,1, .1)kinNNnnkXkFff ekN1.2 离散(lsn)傅里叶变换对于实数(shsh)或者复数离散时间序列f0, f1, FN-1,若满足 ，则称为序列fn离散(lsn)傅里叶变换，称2101()(0 ,1, .,1)kNNinnKfXkekN

4、N10|( ) |Nnft 为序列fn逆离散傅里叶变换第六页，共41页。( )cos(210)cos(225)cos(250)cos(2100)X ttttt平稳信号是指分布参数(cnsh)或者分布律随时间不发生变化的信号，也就是统计特性(期望与方差）不随时间变化而变化。第七页，共41页。2sin(2100)0300sin(250)300600sin(225)600800sin(210)8001000ttttXtttt 2sin(210)0300sin(225)300600sin(250)600800sin(2100)8001000ttttXtttt 第八页，共41页。第九页，共41页。1.3

5、 短时傅里叶变换(binhun)为了提取信号的局部特征,例如变形信号在某一时刻的频率、形变突发位置等，1946年Gabor提出(t ch)了短时傅里叶变换，即Gabor 变换，也称加窗傅里叶变换。21/4/2( )tg teGabor变换的基本思想为：取时间函数 作为窗口函数，然后用 通待分析函数相乘，是时间延迟，是窗函数g(t)的中心(zhngxn)，窗函数根据进行时移，然后再进行傅里叶变换：其中 ，窗口函数g(t)起着时限作用， 起着频限作用。该变化具有不变化宽度（由时间宽度决定）和不变的窗口面积4gg()g t,( , )( ) ()( ),g( )i tftRGf t g tedtf

6、tt ,( )()()i ti ttgtg teg tei te第十页，共41页。短时傅里叶变换(binhun)示意图第十一页，共41页。cos(440)0.5( )cos(660)0.51cos(524)0.5ttx ttttt 傅里叶变换(binhun)图短时傅里叶变换(binhun)图第十二页，共41页。 小波变换由法国科学家MORLET于1980年在进行地震数据分析时提出，是强有力的时频分析(处理)工具，是在克服傅立叶变换缺点的基础上发展而来的。已成功应用于很多领域(ln y)，如信号处理、图像处理、模式识别等。 小波变换的一个重要性质是它在时域和频域均具有很好的局部化特征，它能够提供

7、目标信号各个频率子段的频率信息。这种信息对于信号分类是非常有用的。 小波变换一个信号为一个小波系数，这样一个信号可由小波系数来刻画。第十三页，共41页。 小波变换是一个平方可积分函数f(t)与一个在时频域上均具有良好局部(jb)性质的小波函数(t)的内积： 式中，表示内积，a0 ,为尺度因子，b为位移因子，*表示复数共轭，a,b(t)称为小波基函数。2.1 连续(linx)小波变换小波函数时间(shjin)频率窗(t)称为母小波，(t)必须满足容许性条件：*,1( , ),(t)()fa btbWa bffdtaa第十四页，共41页。部分部分(b fen)(b fen)小波波形小波波形第十五页

8、，共41页。小波分类小波分类(fn li)(fn li)的标准的标准支撑长度：即当时间或频率趋向于无穷大时，它们从一个有限(yuxin)值收敛到0，长度越小，对奇异点的区分效果越好。对称性：对称性越好，越能保证信号不失真（不产生畸变），越能提高信号的重构精度。正则性：它在对信号或图像的重构获得较好的平滑效果作用上是非常有用的。第十六页，共41页。第十七页，共41页。将小波函数沿时间轴向右移动一个单位(dnwi)时间，然后重复步骤(1)、(2)求出此时的小波变换系数C，直到覆盖完整个信号长度，如图所示；将所选择的小波函数尺度伸缩一个单位，然后(rnhu)重复步骤(1)、(2)、(3)，如图所示；

9、对所有(suyu)的尺度伸缩重复步骤(1)、(2)、(3)、(4)。第十八页，共41页。连续(linx)小波变换实例第十九页，共41页。2.2 离散(lsn)小波变换在实际应用中，需要对尺度因子a和位移因子b进行离散化处理，可以取： ， m,n为整数，a0为大于1的常数，b0为大于0的常数，a和b的选取与小波(t)的具体形式有关。离散小波函数表示(biosh)为：相应的离散小波变换可以表示(biosh)为：当a0=2，b0=1时，离散小波变换称为二进离散小波变换，这样便于分析，并且适合于在计算机上进行高效的运算。00,mmoaabnb a 00,0000011mmm nmmmtnb atatn

10、baaa *,fm nm nWm nff tt dt第二十页，共41页。2.2.1 一阶滤波(lb)：近似与细节在小波分析中，近似值是大的缩放因子计算的系数，表示信号的低频分量，而细节值是小的缩放因子计算的系数，表示信号的高频(o pn)分量。实际应用中，信号的低频分量往往是最重要的，而高频(o pn)分量只起一个修饰的作用，如 声音。第二十一页，共41页。第二十二页，共41页。2.2.2多尺度(chd)分解对信号的高频分量(fn ling)不再分解，而将信号的低频部分继续分解.实际中， 分解的级数取决于要分析的信号数据特征及用户的具体需要，例如长度为N的信号，最多能分成log2N层。在实际中

11、，可以选择合适的分解层数。下图为三层多尺度分解树结构，原始信号S的多尺度分解为：S=cA3+Cd3+cD2+Cd1。第二十三页，共41页。第二十四页，共41页。2.2.3小波包分解(fnji)小波分析(fnx)是将信号分解为近似与细节两部分，近似部分又可以分解成第二层近似与细节，可以这样重复下去。对于一个N层分解来说， 有N+1个分解信号的途径。而小波包分析(fnx)的细节与近似部分一样，也可以分解，对于N层分解，它产生2N个不同的途径。第二十五页，共41页。三三 小波变换的一些小波变换的一些(yxi)(yxi)应用应用3.1小波包去噪加噪信号数学模型为f(t)=s(t)+n(t)，s(t)是

12、原信号，n(t)是随机白噪声，满足(mnz)En(t)=0和Dn(t)=2。设(t)为小波函数，n(t)的小波包变换为Wn(j,t)=n(t)j(t)= jRn ttu dun(t)的小波包系数的期望(qwng)和方差分别为： E(|Wn(j,t)|2)=0 D(|Wn(j,t)|2)= 2 tj第二十六页，共41页。3.1.1小波包去噪步骤(bzhu) 选择小波基并确定最佳分解(fnji)的层次,对信号 进行小波包分解(fnji); 对步骤(1)获得的小波包树,选择一定的嫡标准,计算最优树; 估计阈值,并应用该阈值对最优树的小波包系数进行阈值量化; 将经量化处理的小波包系数,重构回原始信号。

13、 小波包阈值消噪有两个关键点：1、如何估计阈值；2 如何利用阈值量化小波包系数。第二十七页，共41页。熵的确定(qudng)熵：用来(yn li)确定最优树的标准，熵值越小，对应的小波包基越好。1）香农熵：约定(yudng)0log(0)=0，则香农熵定义为： 22E slogiiss 2）P范数熵：若P1，在lp范数意义上定义E(s)= ,则： E(s)=PisPPiPiss 2logiis3)对数能量熵 E(si)= ，0log(0)=0,则有 E(s)=2logiis4)阈值熵： E(s)= 式中，是阈值，且 0.10iiss第二十八页，共41页。阈值(y zh)选择准则(1)基于无偏似

14、然估计原理的Rigrsure规则；W为一向量，其元素为小波系数的平方，并按由小到大的顺序排列，W=w1,w2,wn，且w1w2wn，再设一向量R，其元素为：ri= n-2i-(n-i)w+ /n (i=1,2,.,n)以R元素中的最小值rb为风险值，由rb的下标变量(binling)b求出对应的wb，则阈值T1为： T1= 1ikkwbw（2）通用(tngyng)阈值T1(sqtwolog准则) T2= 2log n第二十九页，共41页。（3）启发式的stein无偏风险阈值(y zh)T3（Heursure）准则设为n个小波系数(xsh)的平方和，令= ，= ，则T3=nn322log nn2

15、12min,TT T (4)基于极大极小原理的Minimax方法该准则采用的也是一种固定阈值，它产生一个最小均方误差的极值。具体(jt)的阈值选取规则为： T4=20.39360.1829log32032nnn第三十页，共41页。阈值(y zh)量化函数的选取阈值量化是应用所估计的阈值T，对小波系数进行的处理。目前，阈值量化函数主要采用两种方法。一种(y zhn)是硬阈值法，当小波系数大于该阈值时，保留原值，否则置零，其公式为：另一种(y zhn)是软阈值法，当小波包系数大于该阈值时，向着减小系数幅值的方向作一个收缩，否则置零，其公式为：0iiiiyyTyyTsgn0iiiiiyyyTyyT式中，sgn（）为符号(fho)函数。第三十一页，共41页。阈值准则heursuresqtwologrigrsuremininmaxSNR11.006228.714311.006221.9542rmse1.79760.74161.79761.0398第三十二页，共41页。阈值量化函数硬阈值法软阈值法SNR13.939128.7143RMSE1.55240.7416第三十三页，共41页。系统