《小波分析(讲稿).ppt》由会员分享，可在线阅读，更多相关《小波分析(讲稿).ppt（45页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、小波分析小波分析 Wavelet AnalysisWavelet Analysis已学知识回顾：傅氏级数、傅氏变换实现了信号分析从时域转换至频域。傅氏级数、傅氏变换实现了信号分析从时域转换至频域。引言引言:从一个例子映射出小波分析在信号处理中的作用从一个例子映射出小波分析在信号处理中的作用桌子上一碗果冻和一碗炸酱面，一个勺子和一双筷子桌子上一碗果冻和一碗炸酱面，一个勺子和一双筷子问题：请选择吃果冻和吃炸酱面分别用什么？问题：请选择吃果冻和吃炸酱面分别用什么？择勺子择勺子吃果冻吃果冻，选择筷子吃选择筷子吃炸酱面炸酱面果冻成分单一，外表圆滑；果冻成分单一，外表圆滑；炸酱面成分多样，外表不光滑炸酱面
2、成分多样，外表不光滑勺子的外表也是圆滑的勺子的外表也是圆滑的；筷子的外表是多棱角，开度可变。筷子的外表是多棱角，开度可变。Daubechies4小波函数小波函数时域有时域有限的衰减振荡信号限的衰减振荡信号FFT的正弦函数的正弦函数时域无时域无限的等幅振荡限的等幅振荡FFT分析是以正弦信号为基函数，而小波分析是以小波基为基函数。分析是以正弦信号为基函数，而小波分析是以小波基为基函数。对于信号成分单一的稳定信号，使用对于信号成分单一的稳定信号，使用FFT可以得到较好的结果可以得到较好的结果而对于信号成分多样，包含突变因素的信号小波有很大的优势。而对于信号成分多样，包含突变因素的信号小波有很大的优势
3、。一一.FFT.FFT、STFTSTFT到到WaveletWavelet1.Fourier Analysis1.Fourier AnalysisFFT变换是将信号分解成不同频率的正弦波的叠加和，即把信号变换是将信号分解成不同频率的正弦波的叠加和，即把信号投影到一组正交基投影到一组正交基上。上。一一.FFT.FFT、STFTSTFT到到WaveletWavelet1.Fourier Analysis 1.Fourier Analysis 存在的主要问题：存在的主要问题：（1）无时域局部化特性。无时域局部化特性。为了求得傅里叶系数，理论上必须知道时域的全部为了求得傅里叶系数，理论上必须知道时域的全
4、部信息，即傅里叶分析是对信号的总体的统计分析。反过来讲，时域上任一时刻信息，即傅里叶分析是对信号的总体的统计分析。反过来讲，时域上任一时刻信号的变化都将波及整个频谱图，从而导致傅里叶分析在时域上无任何局部化信号的变化都将波及整个频谱图，从而导致傅里叶分析在时域上无任何局部化特性。这不利于非平稳信号、特别是瞬变信号的分析，特性。这不利于非平稳信号、特别是瞬变信号的分析，在频域上是看不出这在频域上是看不出这些突变特征。些突变特征。一一.FFT.FFT、STFTSTFT到到Wavelet Wavelet 1.FourierAnalysis存在的主要问题：存在的主要问题：（2）不能实现时频分析。不能实
5、现时频分析。信号分解转换到频域后，丢失掉了时域的信息，信号分解转换到频域后，丢失掉了时域的信息，频域中某频率或频带内的信息和时域中某时刻或时宽内的信息没有直接的对频域中某频率或频带内的信息和时域中某时刻或时宽内的信息没有直接的对应关系，即不能给出应关系，即不能给出某一指定频带内的时域图形某一指定频带内的时域图形。这种对应关系称为时频分。这种对应关系称为时频分析，所以傅里叶分析不能进行时频分析，而时频分析在工程中却相当有用。析，所以傅里叶分析不能进行时频分析，而时频分析在工程中却相当有用。一一.FFT、STFT到到Wavelet1.FourierAnalysis存在的主要问题：存在的主要问题：（
6、3）傅氏分析采用窗宽固定的窗函数。）傅氏分析采用窗宽固定的窗函数。为了分析提取信号的低频成分，为了分析提取信号的低频成分，T0应应取较大值，且频率分辩率较高；为了分析提取信号的高频成分，取较大值，且频率分辩率较高；为了分析提取信号的高频成分，T0应取较小值，应取较小值，时域分辩率较高，而时域分辩率较高，而对频率分辨率要求不高。对频率分辨率要求不高。但但T0固定时，两者不能同时满固定时，两者不能同时满足。足。2.2.短时傅里叶变换短时傅里叶变换 STFTSTFT（Short-Time Fourier TransformShort-Time Fourier Transform）STFT将信号在时域
7、上加窗函数，然后进行傅立叶变换，再在时域上将信号在时域上加窗函数，然后进行傅立叶变换，再在时域上移动窗函数，最后完成连续重叠变换，得到与时间有关的信号频谱的描述。移动窗函数，最后完成连续重叠变换，得到与时间有关的信号频谱的描述。从而在时频域得到一个信号能量的三维分布。从而在时频域得到一个信号能量的三维分布。表示对时域信号表示对时域信号f，以以w代表的窗函数加窗，进行短时傅氏变代表的窗函数加窗，进行短时傅氏变换，其结果是频率换，其结果是频率和时移和时移 b 的函数。的函数。主要缺陷：主要缺陷：STFT的窗函数一旦确定，就不能再变换。对于频率成分较多的窗函数一旦确定，就不能再变换。对于频率成分较多
8、的信号，很难找到一个最合适的窗函数，从而很难获得一个最佳的分析的信号，很难找到一个最合适的窗函数，从而很难获得一个最佳的分析精度。精度。2.STFT2.STFT（Short-Time Fourier TransformShort-Time Fourier Transform）3.Wavelet Analysis3.Wavelet Analysis小波变换是时间尺度（时间频率）分析方法，具有多分辨率分析的小波变换是时间尺度（时间频率）分析方法，具有多分辨率分析的特点，即窗口大小固定但其形状可改变，时间窗和频率窗可改变的时频特点，即窗口大小固定但其形状可改变，时间窗和频率窗可改变的时频局部化分析方
9、法。局部化分析方法。在低频部分具有较高的频率分辨率和较低的时间分辨率；在低频部分具有较高的频率分辨率和较低的时间分辨率；在高频部分具有较高的时间分辨率和较低的频率分辨率。在高频部分具有较高的时间分辨率和较低的频率分辨率。优点：优点：适合探测正常信号中夹带的瞬态反常现象并展示其成分，所以被适合探测正常信号中夹带的瞬态反常现象并展示其成分，所以被誉为分析信号的显微镜，利用小波变换进行动态系统故障检测与诊断具誉为分析信号的显微镜，利用小波变换进行动态系统故障检测与诊断具有良好的效果。有良好的效果。4.4.小波分析的发展史小波分析的发展史（1 1）19101910年年HarrHarr提出的第一个小波规
10、范正交基提出的第一个小波规范正交基;（2 2）19841984年法国地质学家年法国地质学家MorletMorlet和理论物理学家和理论物理学家GrossmanGrossman提出了连续小波提出了连续小波变换的概念、变换的概念、19861986年法国数学家年法国数学家MeyerMeyer创造性地构造出了具有一定衰减性的创造性地构造出了具有一定衰减性的光滑化函数光滑化函数正交小波函数正交小波函数,标志着小波热潮的开始。标志着小波热潮的开始。（3 3）19871987年年,法国人法国人MallatMallat提出了多分辨分析的概念提出了多分辨分析的概念,为统一地构造小波函为统一地构造小波函数奠定了基
11、础数奠定了基础,同时给出了以他的名字命名的小波分解与重构算法。同时给出了以他的名字命名的小波分解与重构算法。（4 4）19881988年年,DaubechiesDaubechies构造了具有有限支集的正交小波基构造了具有有限支集的正交小波基,至此小波分析至此小波分析的系统理论初步得到建立。的系统理论初步得到建立。（5 5）19901990年年,崔锦泰和王建中构造了基于样条的半正交小波函数崔锦泰和王建中构造了基于样条的半正交小波函数,使得小波使得小波分析的系统理论得到完善。分析的系统理论得到完善。二二.小波变换的基本概念小波变换的基本概念1.小波定义小波定义2.设设是基本小波的能量，则是基本小波
12、的能量，则的能量为的能量为2.几种常见的小波函数几种常见的小波函数3.定义解释定义解释小波变换的基本思想来小波变换的基本思想来源于源于伸缩伸缩和和平移平移方法。方法。（1）尺度伸缩（）尺度伸缩（Scaling）对波形的尺度伸缩就是对波形的尺度伸缩就是在时间轴上对信号进行压缩在时间轴上对信号进行压缩与伸展，如右图所示，时间与伸展，如右图所示，时间尺度尺度a是反比于频率。是反比于频率。（2）时间平移（）时间平移（Shifting）是指小波函数在时间轴上的波形平行移动。如下图所示。是指小波函数在时间轴上的波形平行移动。如下图所示。二二.小波变换的基本概念小波变换的基本概念4.小波变换定义小波变换定义
13、称称为小波系数为小波系数对小波变换的研究，实质上是对小波系数的研究。对小波变换的研究，实质上是对小波系数的研究。与与Fourier变换对比：从变换对比：从可以看到，两者是在基函可以看到，两者是在基函数上作了变化，实质仍然是不同成分的叠加。数上作了变化，实质仍然是不同成分的叠加。二二.小波变换的基本概念小波变换的基本概念5.信号的重构（信号的重构（逆小波变换逆小波变换InverseWaveletTransform）信号信号f(t)可由小波系数进行重构可由小波系数进行重构二二.小波变换的基本概念小波变换的基本概念 例：例：小波分解与重构小波分解与重构二二.小波变换的基本概念小波变换的基本概念 例：
14、例：小波分解与重构小波分解与重构6.连续小波运算的基本步骤连续小波运算的基本步骤1.选择一个小波函数，将其与要分析的信号起点选择一个小波函数，将其与要分析的信号起点对齐；对齐；2.计算在这一时刻要分析的信号与小波函数的逼计算在这一时刻要分析的信号与小波函数的逼近程度，即计算小波变换系数近程度，即计算小波变换系数C。C越大，说明越大，说明此刻信号与所选小波函数波形越相近此刻信号与所选小波函数波形越相近;（如图（如图1）3.将小波函数沿着时间轴向右移动一个单位时间，将小波函数沿着时间轴向右移动一个单位时间，然后重复然后重复1、2步骤，求出此时的小波变换系数步骤，求出此时的小波变换系数C，直到覆盖完
15、信号的时间长度，直到覆盖完信号的时间长度;（如图（如图2）4.将所选择的小波函数尺度伸缩一个单位，然将所选择的小波函数尺度伸缩一个单位，然后重复后重复1、2、3步骤（如图步骤（如图3）5.对所有的小波函数尺度重复对所有的小波函数尺度重复1、2、3、4步骤。步骤。最后，将得到使用不同尺度评估信号在不同时间最后，将得到使用不同尺度评估信号在不同时间段的系数，这些系数就表征了原始信号在这些小段的系数，这些系数就表征了原始信号在这些小波函数的投影大小。波函数的投影大小。7.小波尺度与频率的关系小波尺度与频率的关系尺度越大，意味着小波函数在时间上越长，即被分析的信号区间越长。尺度越大，意味着小波函数在时
16、间上越长，即被分析的信号区间越长。因此，尺度越大意味着因此，尺度越大意味着对信号的较长的时间段进行近似比较，对信号的较长的时间段进行近似比较，主要获取主要获取的是信号的低频特性。反之，尺度越小，意味着只对信号的较小的细的是信号的低频特性。反之，尺度越小，意味着只对信号的较小的细节部进行比较，因此主要获得的是信号的高频特性。节部进行比较，因此主要获得的是信号的高频特性。小尺度小尺度压缩的小波压缩的小波快速变换的细节快速变换的细节高频部分高频部分大尺度大尺度拉伸的小波拉伸的小波缓慢变换的粗部缓慢变换的粗部低频部分低频部分8.离散小波变换（离散小波变换（Dispersed Wavelet Trans
17、form DWT)对连续信号进行小波变换其工作量非常巨大，而且由于大多数的采样对连续信号进行小波变换其工作量非常巨大，而且由于大多数的采样信号是离散的，所以在实践中是采用离散小波变换来实现。信号是离散的，所以在实践中是采用离散小波变换来实现。（1）一阶滤波：近似与细节）一阶滤波：近似与细节 大多数信号低频部分是对分析有意义的，而高频部分则多是噪声成分。大多数信号低频部分是对分析有意义的，而高频部分则多是噪声成分。因而在实际信号处理时，往往是去除高频成分，而保留低频成分。因而在实际信号处理时，往往是去除高频成分，而保留低频成分。信号的近似部分就是信号中大的、低频成分；细节部分就是信号局部、信号的
18、近似部分就是信号中大的、低频成分；细节部分就是信号局部、高频成分。高频成分。(A-Approximation;D-Detail)8.离散小波变换离散小波变换（1）一阶滤波：近似与细节）一阶滤波：近似与细节（2）离散小波的多尺度分解）离散小波的多尺度分解按照上述一阶滤波的过程，信号的低频部分可按照上述一阶滤波的过程，信号的低频部分可以被继续分解，从而实现了小波的尺度分解以被继续分解，从而实现了小波的尺度分解举例：举例：对一个实际的信号进行小波分解对一个实际的信号进行小波分解说明：说明：信号的最大分解层数计信号的最大分解层数计算如下，如果信号的长度是算如下，如果信号的长度是N,则可以则可以最多分解
19、最多分解层层。利用小波的多尺度分解，可以利用小波的多尺度分解，可以对信号进行消噪以及特定频率对信号进行消噪以及特定频率的滤波。而后对滤波信号进行的滤波。而后对滤波信号进行谱分析。谱分析。（3）小波重构（合成）小波重构（合成）小波分析的另一个主要方面就是在分析、比较、处理（去掉高频噪声小波分析的另一个主要方面就是在分析、比较、处理（去掉高频噪声等处理）后，根据得到的小波变换系数来重构信号。这个过程就是逆离散等处理）后，根据得到的小波变换系数来重构信号。这个过程就是逆离散小波变换（小波变换（IDWT）。）。（4）多尺度分解与重构）多尺度分解与重构9.小波包分解小波包分解 小波分析中，是将信号分解成
20、低频的逼近部分与高频的细节部分，而小波分析中，是将信号分解成低频的逼近部分与高频的细节部分，而后逐次对低频部分进行再次分解。后逐次对低频部分进行再次分解。在实际应用中，希望提高高频频带信号的频率分辨率。小波包是在全在实际应用中，希望提高高频频带信号的频率分辨率。小波包是在全频带对信号进行多次的频带划分，继承了小波变换所具有的良好时频局部频带对信号进行多次的频带划分，继承了小波变换所具有的良好时频局部化优点，由于继续对小波变换没有再分解的高频频带作进一步的分解，从化优点，由于继续对小波变换没有再分解的高频频带作进一步的分解，从而提高了频率分辨率。而提高了频率分辨率。小波分解结构小波分解结构小波包
21、分解结构小波包分解结构三三.小波变换在信号分析中的应用小波变换在信号分析中的应用 小波变换同其它方法相比，具有两个重要的特点：小波变换同其它方法相比，具有两个重要的特点：多分辨率的时频局部化多分辨率的时频局部化分析、快速线性多通道带通滤波。分析、快速线性多通道带通滤波。其应用也可以从这两点来入手，下面就几其应用也可以从这两点来入手，下面就几个典型方面的应用作简要描述。个典型方面的应用作简要描述。1、滤波。、滤波。根据小波变换的滤波特性，可以将那些希望剔除的频率范围对应的根据小波变换的滤波特性，可以将那些希望剔除的频率范围对应的子空间序列置零，然后进行重构运算，就得到了一个同时具有高通、低通、子
22、空间序列置零，然后进行重构运算，就得到了一个同时具有高通、低通、多个带通和多个带阻的多通道滤波器。多个带通和多个带阻的多通道滤波器。2、信号降噪。、信号降噪。信号中的噪声可分为两种情况处理，其一是确定性噪声处理，信号中的噪声可分为两种情况处理，其一是确定性噪声处理，噪声的频率或频率范围可以预先确定，这时利用小波变换的滤波特性即可；噪声的频率或频率范围可以预先确定，这时利用小波变换的滤波特性即可；其二是不确定性噪声的处理，噪声的频率或频率范围不可预知，例如白噪声其二是不确定性噪声的处理，噪声的频率或频率范围不可预知，例如白噪声的频率就几乎覆盖整个频率轴，此时就需要利用小波变换的模极大特性，模的频
23、率就几乎覆盖整个频率轴，此时就需要利用小波变换的模极大特性，模极大代表了重要信号的特征。在实际应用中可以对细节信号设立一个门槛值，极大代表了重要信号的特征。在实际应用中可以对细节信号设立一个门槛值，保留大于门槛值的极值（模极大），而将小于门槛值的毛刺置零，然后进行保留大于门槛值的极值（模极大），而将小于门槛值的毛刺置零，然后进行重构运算，就可以获得近似较好的除去噪声的原信号。重构运算，就可以获得近似较好的除去噪声的原信号。3 3、信号非平稳特性分析、信号非平稳特性分析 在各种振动信号中常常存在有一些突变信号，它们在多数情况下都在各种振动信号中常常存在有一些突变信号，它们在多数情况下都对应于设备
24、的故障等因素，突变信号通常分为边缘跳变和峰值跳变两类，对应于设备的故障等因素，突变信号通常分为边缘跳变和峰值跳变两类，可以将它们等效地认为在信号上叠加一个阶跃信号和脉冲信号。这种突可以将它们等效地认为在信号上叠加一个阶跃信号和脉冲信号。这种突变信号的小波变换结果通常反映为过极值点和零点，例如使用二次样条变信号的小波变换结果通常反映为过极值点和零点，例如使用二次样条函数作为小波函数，则阶跃信号将反映为极值点，脉冲信号将反映为过函数作为小波函数，则阶跃信号将反映为极值点，脉冲信号将反映为过零点。零点。4、用于机器运行状态监测和故障诊断、用于机器运行状态监测和故障诊断小波包能量谱监测。实际振动中一些
25、常见的摩擦、冲击等信号，一般小波包能量谱监测。实际振动中一些常见的摩擦、冲击等信号，一般不能以某些正弦分量来表示。因此，有时采用按频带进行能量监测的方不能以某些正弦分量来表示。因此，有时采用按频带进行能量监测的方法，比频谱分析更为合理。法，比频谱分析更为合理。实际应用实际应用1.信号时频分析信号时频分析A4信号是不平稳的低频波动，是由于电铲挖掘过程中变工况引起的，并非故障，信号是不平稳的低频波动，是由于电铲挖掘过程中变工况引起的，并非故障，但它使振动信号变得极不平稳，不便于分析，在信号分析和重构时可以不予考虑。但它使振动信号变得极不平稳，不便于分析，在信号分析和重构时可以不予考虑。实际应用实际
26、应用2.小波消噪小波消噪在实际信号测量中，传感器、传输线、电源等所带来的背在实际信号测量中，传感器、传输线、电源等所带来的背景噪声，往往使测量结果产生误差，严重时甚至可能淹没有用景噪声，往往使测量结果产生误差，严重时甚至可能淹没有用信号，使测量结果不能正确反映被测对象的真实状态，降低了信号，使测量结果不能正确反映被测对象的真实状态，降低了信号分析的可信度，因此信号消噪是信号处理的首要问题。信号分析的可信度，因此信号消噪是信号处理的首要问题。白噪声的小波变换系数模值随分解尺度的增加而变的越来白噪声的小波变换系数模值随分解尺度的增加而变的越来越小，而信号的小波变换系数模值随分解尺度的增加而变的越越
27、小，而信号的小波变换系数模值随分解尺度的增加而变的越来越大，故可对若干尺度上的小波系数设置阈值，将分解尺度来越大，故可对若干尺度上的小波系数设置阈值，将分解尺度上的噪声所对应的小波系数进行阈值化置零，保留有效信号所上的噪声所对应的小波系数进行阈值化置零，保留有效信号所对应的小波系数，然后进行重构，则重构后的信号就是基于小对应的小波系数，然后进行重构，则重构后的信号就是基于小波变换的消噪信号。波变换的消噪信号。S：原始信号：原始信号cA1,cA2,cA3：分解各层的逼近信号：分解各层的逼近信号cD1,cD2,cD3：分解各层的细节信号：分解各层的细节信号如果将如果将cD1,cD2,cD3的分解系
28、数都设置为零，的分解系数都设置为零，而只保留而只保留cA3,而后进行小波重构，得到的信号而后进行小波重构，得到的信号就是消噪之后的信号。就是消噪之后的信号。分解结构图分解结构图分解后的各部分重构分解后的各部分重构消噪前后对比消噪前后对比实际应用实际应用2.小波消噪小波消噪实际应用实际应用2.能量谱检测能量谱检测实际应用实际应用2.能量谱检测能量谱检测实际应用实际应用2.能量谱检测能量谱检测当轮齿齿面出现点蚀时，在正常振动的基础上出现了很多小的波峰，当轮齿齿面出现点蚀时，在正常振动的基础上出现了很多小的波峰，这些波峰就是点蚀坑引起的，如图这些波峰就是点蚀坑引起的，如图1所示。所示。从信号频谱图（
29、图从信号频谱图（图2）上可以看到，在信号的高频段（约）上可以看到，在信号的高频段（约1750Hz3250hz）振动能量明显提升。点蚀坑将引起脉冲冲击，从而使信号在较）振动能量明显提升。点蚀坑将引起脉冲冲击，从而使信号在较高频域段内幅值有较大提高。尤其是在齿轮固有频率处，由于共振使振高频域段内幅值有较大提高。尤其是在齿轮固有频率处，由于共振使振动能量明显加大。动能量明显加大。利用小波包的多尺度分解能力，把信号分解到多个频率带，利用小波包的多尺度分解能力，把信号分解到多个频率带，通过检测各频率带信号的能量变化，可以对点蚀故障进行识别通过检测各频率带信号的能量变化，可以对点蚀故障进行识别实际应用实际
30、应用2.能量谱检测能量谱检测点蚀发生前后齿轮振动信号各频段点蚀发生前后齿轮振动信号各频段相对能量相对能量及其比值及其比值频频段（段（HzHz）W1W1W2W2W3W3W4W4W5W5W6W6W7W7W8W8 0 0437.5437.5437.5437.5875.0875.0875.0875.01312.51312.51312.51312.51750.01750.01750.01750.02187.52187.52187.52187.52625.02625.02625.02625.03062.53062.53062.53062.53500.03500.0点点蚀蚀前相前相对对能量能量A A（）（）
31、99.573299.57320.38730.38730.01580.01580.02110.02110.00110.00110.00120.00120.00010.00010.00020.0002点点蚀蚀后相后相对对能量能量B B（）（）97.509497.50942.40762.40760.03750.03750.02260.02260.00820.00820.01310.01310.00060.00060.00100.0010比比值值(B/A)(B/A)0.97930.97936.21646.21642.37342.37341.07111.07117.45457.454510.916710
32、.91676.00006.00005.00005.0000注：注：相对能量值各分解频带信号能量占信号总能量的百分比相对能量值各分解频带信号能量占信号总能量的百分比结论：结论：从表从表1中可以看出，出现点蚀后，信号高频部分振动能量变化较大，中可以看出，出现点蚀后，信号高频部分振动能量变化较大，其中包含啮合轮齿固有频率的其中包含啮合轮齿固有频率的W6分解频段能量提高最明显，点蚀时相分解频段能量提高最明显，点蚀时相对能量是未点蚀时的对能量是未点蚀时的10.9167倍。倍。实际应用实际应用3.奇异值检测奇异值检测图图频率突变信号频率突变信号传统的功率谱等方法只能给出传统的功率谱等方法只能给出信号全频尺
33、度特征，而没有办信号全频尺度特征，而没有办法说明信号发生突变的具体位法说明信号发生突变的具体位置。而在信号分析中，突变点置。而在信号分析中，突变点（奇异点）往往代表了故障的（奇异点）往往代表了故障的发生，在时域上则代表了故障发生，在时域上则代表了故障发生的具体时间点。因此功率发生的具体时间点。因此功率谱对于信号的突变没有时域分谱对于信号的突变没有时域分析能力。析能力。频率突变信号的奇异性检测频率突变信号的奇异性检测频率突变信号的的小波分解频率突变信号的的小波分解小波变换具有空间局部小波变换具有空间局部化性质，即可以描述信化性质，即可以描述信号在不同频率带中不同号在不同频率带中不同时刻变换情况，
34、因此可时刻变换情况，因此可以利用小波来分析信号以利用小波来分析信号的奇异性，判断信号中的奇异性，判断信号中的突变情况。的突变情况。幅值突变信号的小波分解幅值突变信号的小波分解幅值突变信号的奇异性检测幅值突变信号的奇异性检测 当齿轮传动系统出现齿轮的擦伤、点蚀、裂纹等故障时往往会引起振当齿轮传动系统出现齿轮的擦伤、点蚀、裂纹等故障时往往会引起振动信号峰值的突变，突变点称为奇异点。同平稳信号相比，奇异点包含了动信号峰值的突变，突变点称为奇异点。同平稳信号相比，奇异点包含了更为丰富的故障信息。传统的更为丰富的故障信息。传统的FourierFourier变换可以在整体上给出信号的奇异变换可以在整体上给
35、出信号的奇异性表述，但无法刻画局部的奇异性。尤其对于轮齿点蚀等发展过程较为缓性表述，但无法刻画局部的奇异性。尤其对于轮齿点蚀等发展过程较为缓慢的故障，很难在故障出现的早期进行成功识别。慢的故障，很难在故障出现的早期进行成功识别。信号奇异性提取要求对信号进行局部化分析，而小波分析具有良好的信号奇异性提取要求对信号进行局部化分析，而小波分析具有良好的时频（尺度）局部化能力。因此，小波分析可以在较低信噪比的信号中时频（尺度）局部化能力。因此，小波分析可以在较低信噪比的信号中检测到故障信号。检测到故障信号。应用实例应用实例原理：原理：当齿轮出现以上故障时，会引起信号在故障处出现频率或幅值突当齿轮出现以
36、上故障时，会引起信号在故障处出现频率或幅值突变，利用小波可以较好实现对这些突变的提取。变，利用小波可以较好实现对这些突变的提取。下面以对刘桥下面以对刘桥2 2矿主井提升机减速器的检测为例进行说明。矿主井提升机减速器的检测为例进行说明。对该齿轮的测试记录了大齿轮齿顶擦伤的发展过程（由于小齿轮是硬对该齿轮的测试记录了大齿轮齿顶擦伤的发展过程（由于小齿轮是硬齿面，擦伤主要发生在大齿轮齿面）。故障首先发生在大齿轮某几个轮齿齿面，擦伤主要发生在大齿轮齿面）。故障首先发生在大齿轮某几个轮齿的齿顶部，而后逐渐扩展到每个轮齿的齿顶部。的齿顶部，而后逐渐扩展到每个轮齿的齿顶部。从图中对比发现：同一种情况下，信号
37、小波分解各尺度高频细节部分从图中对比发现：同一种情况下，信号小波分解各尺度高频细节部分中对应擦伤故障发生位置均有突出成分，尺度越低（中对应擦伤故障发生位置均有突出成分，尺度越低（压缩的小波压缩的小波），突，突出越明显；随着擦伤的发展，出现突出成分的次数也逐渐提高，反映了出越明显；随着擦伤的发展，出现突出成分的次数也逐渐提高，反映了擦伤从初期几齿到所有齿的发展特征。擦伤从初期几齿到所有齿的发展特征。待分解信号包含了待分解信号包含了5个连续轮齿的啮合过程。个连续轮齿的啮合过程。由于信号的分解高频细节体现了故障信号的突变成分，因此只画出原始由于信号的分解高频细节体现了故障信号的突变成分，因此只画出原
38、始信号信号（s）和分解和分解5尺度的各高频细节重构（尺度的各高频细节重构（d1d5）。如图所示。）。如图所示。应注意的问题应注意的问题（1）小波作为非平稳信号的良好分析工具在机械故障信号分析中有很）小波作为非平稳信号的良好分析工具在机械故障信号分析中有很大的优势，但认为小波分析可以完全取代传统大的优势，但认为小波分析可以完全取代传统Fourier分析的观点是不全分析的观点是不全面的，各有其优缺点。面的，各有其优缺点。Fourier分析突出的特点在于能得出组成时域信号分析突出的特点在于能得出组成时域信号的各频率成分的相关信息。的各频率成分的相关信息。（2）小波分析虽然可以揭示信号的局部化特征，但
39、当基函数、分解层）小波分析虽然可以揭示信号的局部化特征，但当基函数、分解层数等选择不当时，可能产生对原始信号的歪曲。同时，同一信号的不同小数等选择不当时，可能产生对原始信号的歪曲。同时，同一信号的不同小波分解，所得结果自然也不相同，如何判断结果的准确性成为一个重要的波分解，所得结果自然也不相同，如何判断结果的准确性成为一个重要的问题。问题。（3）有多种小波基函数供选择，但目前往往依靠使用人员对信号的认识）有多种小波基函数供选择，但目前往往依靠使用人员对信号的认识和掌握，经验起到很大作用。和掌握，经验起到很大作用。（4）小波基函数众多，不同信号分析需要不同的基函数，因此，基函）小波基函数众多，不
40、同信号分析需要不同的基函数，因此，基函数的创建是小波发展的一个主要方向。数的创建是小波发展的一个主要方向。小波分析的基本算法已作为一种工具箱在商品化的软件系统中实现。小波分析的基本算法已作为一种工具箱在商品化的软件系统中实现。例如：例如：由美国国家仪器公司（由美国国家仪器公司（NATIONALINSTRUMENTS简称简称NI公司）推出的公司）推出的LabVIEW（LaboratoryVirtualInstrumentEngineeringWorkbench）新）新版本版本5.1的信号处理软件中提供了小波分析工具箱。的信号处理软件中提供了小波分析工具箱。LabVIEW采用图形化采用图形化编程语
41、言，构造各种虚拟仪器，实现信号检测、处理、分析和显示。编程语言，构造各种虚拟仪器，实现信号检测、处理、分析和显示。由美国的由美国的Mathworks软件公司开发的软件公司开发的MATLAB4.2C:STFT中提供了小波分中提供了小波分析工具箱。析工具箱。MATLAB是一种面向科学和工程计算的高级语言，集科学计算、是一种面向科学和工程计算的高级语言，集科学计算、自动控制、信号处理、图像处理、神经网络于一体，功能强大，编程效率高。自动控制、信号处理、图像处理、神经网络于一体，功能强大，编程效率高。在了解小波分析的最基本原理的基础上，只要读懂这些工具箱的使用指令，在了解小波分析的最基本原理的基础上，只要读懂这些工具箱的使用指令，就可以正确应用小波分析方法来处理复杂的信号。就可以正确应用小波分析方法来处理复杂的信号。小波分析工具小波分析工具 谢 谢