小波分析与实例(课堂PPT).ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,2,*,小波分析,1,2025/3/28 周五,小波分析讲解,傅里叶变换与小波分析,小波分析的基本知识,多尺度分析与,Mallat,算法,小波分析的应用,2,2025/3/28 周五,1,、傅里叶变换与小波分析,小波分析是近年来迅速发展起来的一个数学分支。除了在数学学科本身中的价值外，小波分析在许多非数学的领域也有着广泛的应用。,3,2025/3/28 周五,1,、傅里叶变换与小波分析,一、傅里叶变换,对于平稳信号，做完FFT（快速傅里叶变换）后，可以在频谱上看到清晰的四条线，信号包含四个频率成分。,4,2025/3/28 周五,1,、傅里叶变换与小波分析,频率随着时间变化的非平稳信号，进行,FFT,后：,如左图，最上边的是频率始终不变的平稳信号。而下边两个则是频率随着时间改变的非平稳信号，它们同样包含和最上信号相同频率的四个成分。做FFT后，我们发现这三个时域上有巨大差异的信号，频谱（幅值谱）却非常一致。尤其是下边两个非平稳信号，我们从频谱上无法区分它们，因为它们包含的四个频率的信号的成分确实是一样的，只是出现的先后顺序不同。,5,2025/3/28 周五,1,、傅里叶变换与小波分析,可见，傅里叶变换处理非平稳信号有天生缺陷。它只能获取一段信号总体上包含哪些频率的成分，但是对各成分出现的时刻并无所知。因此时域相差很大的两个信号，可能频谱图一样。然而平稳信号大多是人为制造出来的，自然界的大量信号几乎都是非平稳的，所以在比如生物医学信号分析等领域的论文中，基本看不到单纯傅里叶变换这样简单的方法。,事件相关电位,股市折线图,6,2025/3/28 周五,1,、傅里叶变换与小波分析,加窗傅里叶变换（短时傅里叶变换,STFT,）,7,2025/3/28 周五,1,、傅里叶变换与小波分析,窗划分太窄，窗内的信号太短，会导致频率分析不够精准，频率分辨率差。,窗划分太宽，时域上又不够精细，时间分辨率低。,8,2025/3/28 周五,1,、傅里叶变换与小波分析,小波定义：,小,波动性：,9,2025/3/28 周五,小波的3 个特点,小波变换，既具有频率分析的性质，又能表示发生的时间。有利于分析确定时间发生的现象。（傅里叶,变换只,具有频率分析的性质）,小波变换的多分辨度的变换，有利于各分辨度不同特征的提取（图象压缩，边缘抽取，噪声过滤等）,小波变换比快速,Fourier,变换还要快一个数量级。,信号长度为,M,时，,Fourier,变换,（左）,和,小波变换（右）计算复杂性分别如下公式：,10,2025/3/28 周五,1,、傅里叶变换与小波分析,11,2025/3/28 周五,1,、傅里叶变换与小波分析,12,2025/3/28 周五,小波运算的步骤,（,1,）选择小波函数，并与分析信号起点对齐；,（,2,）计算在这一时刻要分析信号与小波函数的逼近程度，即小波变换系数,C,。,C,越大，就意味着此刻信号与所选择的小波函数波形越相近；,（,3,）将小波函数沿时间轴右移一个单位时间，然后重复（,1,）、（,2,）步骤，求出变换系数,C,，直到覆盖整个信号长度；,13,2025/3/28 周五,小波运算的步骤,（,4,）将所选择的小波函数尺度伸缩一个单位，然后重复步骤（,1,）、（,2,）、（,3,）；,（,5,）对所有的伸缩尺度重复步骤（,1,）、（,2,）、（,3,）、（,4,）。,14,2025/3/28 周五,2,、小波分析的基本知识,小波基础术语,:,紧支撑,:,对于函数f(x)，如果自变量x在0附近的取值范围内，f(x)能取到值；而在此之外，f(x)取值为0。,那么这个函数f(x)就是紧支撑函数，而这个0附近的取值范围就叫做紧支撑集。,比如：在(-1,1)之间的高斯函数。,L,（,R,）,:,满足 成立的自变量为实数的实值或复值函数f,的全体。,L,（,0,2,）,：,f,（,x+2,）,=f,（,x,），,15,2025/3/28 周五,2,、小波分析的基本知识,小波定义：,设,L,（,R,）,L,（,R,），在,R,上不几乎处处为,0,，且满足,则称,为小波。,其中,为的傅里叶变换。,16,2025/3/28 周五,2,、小波分析的基本知识,称为依赖参数,a,，,b,的连续小波，叫基本小波或小波。若是窗函数，就叫为窗口小波函数，一般我们恒假定为窗口小波函数。,17,2025/3/28 周五,2,、小波分析的基本知识,a,为尺度参数,18,2025/3/28 周五,2,、小波分析的基本知识,b,为位移参数,19,2025/3/28 周五,2、小波分析的基本知识,小波正变换：,小波逆变换：,是,f,（,t,）在函数 上的投影,。,20,2025/3/28 周五,一维连续小波的例子：,1.Haar,小波,：,2,21,Haar,小波是一组相互正交的函数集，是一个最简单的时域不连续的二进小波，,Haar,的应用十分广泛，常用与图像处理。,一维连续小波的例子,2.Mexico,草帽小波,：,2,22,草帽函数又称为,Marr,小波。其在时域、频域都有很好的局部特性，但它的正交性尺度函数不存在，主要用于信号处理和边缘检测。,一维连续小波的例子：,3.Morlet,小波,：,2,23,式中，,i,表示虚数，,w,表示常数。,Morlet,小波不具有正交性的同时也不具有紧支集。其特点是能够获取信号中的幅值和相应的信息，广泛应用于地球物理信号处理中。,Daubechies(dbN)小波系（多贝西）,多贝西小波是以英格丽多贝西的名字命名的一种小波函数，多贝西小波主要应用在离散型的小波转换，是最常使用到的小波变换。多贝西小波是一种正交小波，所以它很容易进行正交变换。,对于有限长度的小波，应用于快速小波变换时，会有两个实数组成的数列：一是作为高通滤波器的系数，称作小波滤波器；二是低通滤波器的系数，称为调整滤波器（尺度滤波器）。,我们通常以滤波器长度N,来形容滤波器为,dbN,，例如,N=2,的多贝西小波写作,db2,；,N=4,的多贝西小波写作,db4,。,24,2025/3/28 周五,Daubechies(dbN)小波系（多贝西）,图,1.4,25,2025/3/28 周五,小波函数表,26,2025/3/28 周五,小波函数表,27,2025/3/28 周五,2、小波分析的基本知识,连续小波变换,这就是信号f(t)的连续小波变换公式，其中参数a和b都是连续变化的参数，a为尺度参数（在某种意义上就是频率的概念），b是时间参数或平移参数。不严谨地讲，Wf(a,b)指的是对信号f(t)进行小波变换后当频率为a时间为b时的变换值。可以看出，一维信号f(t)经过小波变换后将变成二维信号。,28,2025/3/28 周五,2、小波分析的基本知识,连续小波变换,例：已知一信号,f,(,t,),3sin(100,p,t),2sin(68,p,t),5cos(72,p,t),，且该信号混有白噪声，对该信号进行连续小波变换。小波函数取,db3,，尺度为,1,、,1.2,、,1.4,、,1.6,、,、,3,。其,MATLAB,程序如下：,t,00.01,：,1,；,f,3*sin(100*pi*t),2*sin(68*pi*t),5*cos(72*pi*t),randn(1,，,length(t),；,coefs,cwt(f,，,1,：,0.2,：,3,，,db3,，,plot,),；,title(,对不同的尺度小波变换系数值,),；,Ylabel(,尺度,),；,Xlabel(,时间,),；,29,2025/3/28 周五,2、小波分析的基本知识,连续小波变换,小波变换的系数如图所示的灰度值图表征，横坐标表示变换系数的系号，纵坐标表示尺度，灰度颜色越深，表示系数的值越大。,30,2025/3/28 周五,离散小波变换：在实际运用中，尤其是在计算机上实现，连续小波必须加以离散化。因此，有必要讨论一下连续小波,y,a,，,b,(,t,),和连续小波变换,W,f,(,a,，,b,),的离散化。需要强调指出的是，这一离散化都是针对连续的尺度参数a和连续平移参数b的，而不是针对时间变量t的。在连续小波中，考虑函数这里，,b,R,，,a,R,，且,a,0,，,y,是容许的，为方便起见，在离散化中，总限制,a,只取正值，这样相容性条件就变为,2、小波分析的基本知识,离散小波变换,31,2025/3/28 周五,2、小波分析的基本知识,二进小波变换,32,2025/3/28 周五,2、小波分析的基本知识,二进小波变换,定义：设yj，k(t)L2(R)，且满足(1.64)由此得到的小波,y,j,，,k,(,t,),称为二进正交小波。,33,2025/3/28 周五,3,、多尺度分析与,Mallat,算法,多分辨分析,为了改变信号的分辨率使得人们可以根据特定的目标处理相关的细节，1983,年，,P.J.Burt,与,E.A.Adelson,在计算机视觉的应用中引进了一个能够处理低分辨率图像，同时根据需要进一步提高图像分辨率的多分辨率,Laplace,塔式算法。,1986,年,Mallat,和,Meyer,构造了多分辨分析公式。随着多分辨分析的出现，构造小波的困难得到了较圆满的解决。为了对信号进行较高分辨率的处理，需要一种所谓的,“,增量信息,”,。为此，,Mallat,选用正交小波基作为对,“,增量信息,”,进行数学描述，并最终发展成为了多分辨分析。,34,2025/3/28 周五,3、,多,尺度分析与,Mallat,算法,35,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,参考：,M.Vetterli,”Wavelets and Subband Coding“,Prentice Hall PTR,1995,p.11,36,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,滤波器族：下图是一系列带通滤波器的频域图,37,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,一个信号离散信号x(n)经过这一系列带通滤波器滤波后，将得到一组系数Vi(n)。如下图所示：,这样，我们就把一个信号分解成了不同频率的分量。只要这些带通滤波器的频率能够覆盖整个原信号x(n)的频谱范围，反变换时，把这些不同频率信号，按其分量大小组合起来，就可得到原信号x(n)。这样一组带通滤波器就称为滤波器族。,38,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,滤波器族能实现将信号分为不同频率分量，从而实现分解信号并分析信号的目的。但是在滤波器族的计算中，我们需要指定频域分割方式。研究者们给出了一种分割方式，即均分法，从而引出了子带编码的概念。,子带编码通过使用均分频域的滤波器，将信号分解为若干个子带。这样是可以实现无冗余且无误差地对数据分解及重建目的。但是Mallat在1989年的研究表明，如果只分为2个子带，可以实现更高效的分解效率。从而引入了多分辨率分析（MRA）。,39,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,多分辨率分析：,如果子带编码时将信号带宽先对分为高通（实际为带通）和低通两个部分，对应于两个滤波器。然后对低通部分继续等分。下图为子带编码示意图。,40,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,从图中看出，每次分割保留高通部分的滤波结果，因为这里已经是信号的细节了，而且通常我们分析的信号，其绝大部分能量都在低频部分。所以高频部分的分割可以到此为止，但是低通部分仍然有更多的细节可以划分划分出来，所以将低通部分继续等分。分割迭代进行。,这样做的优点是，我们只需要设计两个滤波器，然后每次迭代将其对分。缺点是，频域的分割方式确定。对于某些信号来说，这样的划分并不是最优的。,41,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,这里仍然有个问题。每次都将频谱分为剩下的一半，那实际上，我们永远也取不到整个频段。就好比一杯水，每次都只许喝一半，那将永远无法把它完全喝完。所以，这样分割后的函数仍然是无限多的。为解决这个问题，终于引出了我们最初想讨论的尺度函数的概念。,在上图中，我们对频域进行分割，当分割到某个频率j时，不再继续分割了，剩下的所有低频部分由一个低通滤波器来表示，这就可以实现对信号频谱的完整分割。这个剩余低通滤波器就是尺度函数。事实上，很容易看出，尺度函数无非就是某级多分辨率分析中的低通滤波器。也就是图中最下面一级的LP。,42,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,load noissin,c=cwt(noissin,1:48,db4);,c=cwt(noissin,1:48,db4,plot);,c=cwt(noissin,2:2:128,db4,plot);,43,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,44,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,S,A,3,D,3,D,2,D,1,设以,V,j,表示图,1.17,分解中的低频部分,A,j,，,W,j,表示分解中的高频部分,D,j,，则,W,j,是,V,j,在,V,j,1,中的正交补，即,45,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,S,A,3,D,3,D,2,D,1,若令,f,j,V,j,代表分辨率为,2,j,的函数,f,L,2,(,R,),的逼近,(,即函数,f,的低频部分或“粗糙像”,),，而,d,j,W,j,代表逼近的误差,(,即函数,f,的高频部分或“细节”部分,),，则上式意味着：,f,N,f,1,f,d,f,2,d,2,d,1,f,N-1,d,N,1,d,2,d,1,所以上式可简写为,这表明，任何函数,f,L,2,(,R,),都可以根据分辨率为,2,N,时,f,的低频部分,(“,粗糙像”,),和分辨率,2,j,(1,j,N,),下,f,的高频部分,(“,细节”部分,),完全重构，这恰好是著名,Mallat,塔式重构算法的思想。,46,2025/3/28 周五,3、,多尺度分析与,Mallat,算法,小波重构,47,2025/3/28 周五,Mallat,算法中仅仅对低频系数进行分解，但是对于有些信号来说，对高频系数进行分解更加合适。小波包分解即将低频系数和高频系统都进行同样的分解，然后选取一个最合适的分解路径。然后通过构建一个代价函数求来对于路径进行评价，选取最优路径。,3、,多尺度分析与,Mallat,算法,48,2025/3/28 周五,4,、,小波分析,的应用,小波的信号分解与求频,小波在图像压缩中的应用,小波变换在图像去噪与图像增强中的应用,机械故障诊断,小波神经网络预测,49,2025/3/28 周五,4,小波的应用,小波的信号分解与求频,clear all,clc,fs=1024;%采样频率,f1=100;%信号的第一个频率,f2=300;%信号第二个频率,t=0:1/fs:1;,s=sin(2*pi*f1*t)+sin(2pi*f2*t);%生成混合信号,tt=wpdec(s,3,dmey);%小波包分解，3代表分解3层,plot(tt),wpviewcf(tt,1);,50,2025/3/28 周五,4,小波的应用,小波的信号分解与求频,65-128Hz,257-320Hz,51,2025/3/28 周五,4,小波的应用,小波在图像压缩中的应用,小波变换的基本思想是用一组小波或基函数表示一个函数或信号，例如图像信号。以哈尔（Haar）小波基函数为例，基本哈尔小波函数（Haar wavelet function）定义如下：,1,当0 x1/2,(x)=-1,当1/2x350),c(I)=2*c(I);,else,c(I)=0.5*c(I);,end,end,nx=waverec2(c,s,sym4);,%分解系数重构,subplot(122);,image(nx);,title(增强图像),%画出增强图像,59,2025/3/28 周五,4,小波的,应用,机械故障诊断,当机械运行发生故障时，其振动信号中往往是首先出现相应的瞬态脉冲波形。能否及时准确地予以捕捉分析，常常是能否及时发现故障，采用相应对策，避免出现重大损失的先决条件。传统的傅里叶分析和时域分析由于需要的数据量较大，难以及时作出有效的诊断，而小波分析具有良好的时域定位特征，只需要少数数据就可以对振动信号在时域和频域进行定量分析，从而为及时发现故障提供了一种有力的分析手段。,60,2025/3/28 周五,4,小波的应用,机械故障诊断,齿轮裂纹和断裂,41%,齿面疲劳,31%,齿面擦伤和划痕,10%,齿面磨损,10%,其他故障类型,8%,61,2025/3/28 周五,4,小波的应用,机械故障诊断,clc;clear all;close all;,load leleccum;%载入信号数据,s=leleccum;,Len=length(s);,ca1,cd1=dwt(s,db1);%采用db1小波基分解,a1=upcoef(a,ca1,db1,1,Len);%从系数得到近似信号,d1=upcoef(d,cd1,db1,1,Len);%从系数得到细节信号,s1=a1+d1;%重构信号,figure;,subplot(2,2,1);plot(s);title(源信号);,subplot(2,2,2);plot(ca1);title(一层小波分解的低频信息);,subplot(2,2,3);plot(cd1);title(一层小波分解的高频信息);,subplot(2,2,4);plot(s1,r-);title(一层小波分解的重构信号);,62,2025/3/28 周五,4,小波的应用,机械故障诊断,63,2025/3/28 周五,4,小波的应用,小波神经网络预测,神,经网络引入预测领域使得预测理论及方法产生了质的飞越。目前神经网络具有分布式、联想。记忆和很强的泛化能力，以及自学习和容错性可以以任意精度逼近非线性函数等优点。,但是，神经网络应用于预测中存在如下问题：,难以确定网络的结构；,训练速度较慢；,容易陷入局部次优点等。,64,2025/3/28 周五,4,小波的应用,小波神经网络预测,小波方法与神经网络结合的方法有两种：,一种是先通过小波对网络流量时间序列进行小波分解，得到小波变换尺度系数序列和小波系数序列，然后输入到一个神经网络中加以训练得到预测。,二是把神经网络隐含层的传输函数用小波函数代替，这样从本质上改变了预测模型的结构，在不影响预测精度的前提下，大大缩短了模型的训练时间，提高了训练速度，克服了神经网络容易陷入局部次优点的缺点，而且把神经网络隐含层的传输函数用小波函数代替、算法易实现和推广。,65,2025/3/28 周五,4,小波的应用,小波神经网络预测,小波神经网络预测模型采用三层小波网络结构，在输入端有,p,个输入，即一起输入含有,p,个时间序列值，,p,表示预测序列值与相关步数。隐含层包含,n,个神经元，输出层有,1,个神经元，输出为第,k+1,个时间序列的预测值。,66,2025/3/28 周五,4,小波的应用,小波神经网络预测,从小波神经网络训练图可以看出一步预测用了,4100epochs；,五步用了,4005epochs，,十步用了,3475epochs；,而没有使用小波的神经网络训练图在训练一步、五步、十步时的训练次数分别为,5739epochs，7317epochs，6277epochs。,可以明显看出小波神经网络对训练时间能够起到很大的提升。,67,2025/3/28 周五,4,小波的应用,小波分析在加工误差分析中的应用,在机械制造中，零件的尺寸、形状、相对位置等的形成，实际上取决于刀具和工件在整个切削运动中的相互位置和相对运动。由于工艺精度或其自身的各种原因，会形成误差,在各种零件的加工中，为了及时、准确的了解工艺系统的状态，反映系统的特性或工艺能力，通常会按加工顺序对加工的零件尺寸进行记录，形成离散的数据点。这些可以通过离散小波变换来进行分析与处理。,68,2025/3/28 周五,4,小波的应用,小波分析在加工误差分析中的应用,我们通过将阈值调整，可以选择性的去处高频成分（随机误差），以获取需要的数据。,上升的曲线就是加工系统的变值系统性误差，它实际上是整个工艺系统中机床、刀具的热变形以及刀具磨损等因素共同作用的结果。,小波变换分析的精髓就在于对那些变化平稳的信息（低频信息），可以在大范围上观察，对于变化很块的信息（高频信息），又可以在小范围上观察。相对于传统的误差分析方法，更加方便。快捷和有效，能使我们更加及时、准确的掌握当前系统，为研究和探讨误差产生的原因、机理打下基础。,69,2025/3/28 周五,Matlab,小波工具箱,70,2025/3/28 周五,谢谢观看！,注：,ppt,内多图来自知乎,咚懂咚懂,咚,。,71,2025/3/28 周五,