故障诊断技术part6基于小波变换的故障诊断方法2

例：假设一信号的主要频率成分是100Hz和400Hz， 如下图所示，通过傅里叶变换对其频率成分进行频域 分析。 上图为原始信号，从 图中看不出100Hz和 400Hz的任何频域信 息。 但从下图的信号频谱 分析中，可以明显看 出信号的频率特性。
从上例中可知，虽然傅里叶变换能够将信号的时域 特征和频域特征联系起来，能分别从信号的时域和 频域进行观察，但却不能把两者有机地结合起来。
定义3： 在L2(R)空间，函数f(x)的范数‖f(x)‖定义为：
f ( x) f ( x) f ( x)dx
2



f ( x) dx
2
定义4： 在L2(R)空间，若：内积<f，g>＝0，则称函数f与 函数g正交。
定义5：
在L2(R)空间，两个函数f(x)与g(x)的卷积定义为：
小波与小波变换
定义12： 我们称满足条件


0
ˆ ˆ 0 ( ) ( ) d d C
2 2
的平方可积函数ψ(x)（即ψ(x) ∈L2(R)）为基本 小波，或小波母函数。
定义13：
函数f(x) ∈L2(R)的连续小波变换定义为：
1 x u W f ( s, x) f s ( x) f (u ) ( )du s s
为了得到更好的时频分析效果，希望δ和ε都非常小， 但是由海森堡测不准定理（Heisenberg Uncertainty Principle）可知， δ和ε是互相制约的， 两者不可能同时都任意小。 （事实上， δ·≥0.5，且仅当g(t)为高斯函数时，等 ε 号成立。）
由此可见，短时傅里叶变换虽然在一定程度上克服了 标准傅里叶变换不具有局部分析能力的缺陷，但它也 存在着自身不可克服的缺陷，即当窗函数g(t)确定后， 矩形窗口的形状就确定了， τ和ω只能改变窗口在相 平面上的位置，而不能改变窗口的形状。 可以说，短时傅里叶变换是具有单一分辨率的分析， 这对分析信号来说是很不利的。 因为，一般来说高频信号持续的时间比较短，低频信 号持续的时间比较长。 为了更好地分析信号，信号的高频成分需要窄的时间 窗，而信号的低频成分需要宽的时间窗。 而单一分辨率无法满足这种要求。
其中，*表示卷积。
因此，Wf(s，x)关于x的傅里叶变换可以表示为：
ˆ ˆ ˆ W f ( s, ) f ( ) ( s )
连续小波变换的定义
傅立叶变换最重要的意义是它引进了频率的概念， 他把一个函数展开成各种频率的谐波的线性叠加， 由此引出了一系列频谱分析的理论。 很多在时域中看不清的问题，在频域中却能一目 了然 。 因此，长期以来，Fourier分析理论不论在数学中 还是工程科学中一直占领着极其重要的地位。
傅立叶分析的实质在于将一个任意的函数f(t)表示 为具有不同频率的谐波函数的线性叠加。即一族 标准函数 {eit R} 的加权求和，从而将对原来函 数的研究转化为对这个叠加的权系数的研究：
1990年，Meyer等出版第一部小波系统性专著 《小波与算子》，共三卷。尤众、王耀东、邓东皋等 译校成中文本（共两册）。这套书详细研究了各种小 波基的构造，小波基与函数空间的关系，小波分析在 复分析、算子论、偏微分方程与分线性分析等方面的 应用。
1991年，邓东皋等在《数学进展》上发表“小波 分析”－国内第一篇小波论文。对国内小波的研究和 应用起了很大的推动作用。
1992年，Daubechies的《小波10讲》系统论述 了正交小波的紧支性、正则性、对称性及时频特性， 介绍了离散小波变换和连续小波变换等。 到此，经典小波理论已基本成熟，1992年以后，在国 际上，重点转向小波的推广和应用。
在国内，由于对小波的研究起步较晚，20世纪90年 代以来，可以说小波的理论研究和应用研究几乎同时 开始。 1994年，形成国内的小波高潮。
小波变换是一种能够在时间－频率两域对信号进行 分析的方法，具有可以对信号在不同范围、不同的 时间区域内进行分析，对噪声不敏感，能够分析到 信号的任意细节等优点，在信号处理领域获得越来 越广泛的应用，被誉为“数学显微镜”。
小波分析和Fourier分析 傅立叶变换是一个十分重要的工具，无论是在一般 的科学研究中，还是在工程技术的应用中，它都发 挥着基本工具的作用。 从历史发展的角度来看，自从法国科学家J.Fourier 在1807年为了得到热传导方程简便解法而首次提出 著名的傅立叶分析技术以来，傅立叶变换首先在电 气工程领域得到成功应用，之后，傅立叶变换迅速 得到越来越广泛的应用，而且理论上也得到了深入 研究。
g ( ) f ( ) F{ f (t )} 的研究。
上式中，各符号的含义： g ( ) 表示频域函数； f ( ) 表示对原函数f(t)的傅里叶变换； g ( ) 表示对频域函数 g ( ) 的傅里叶反变换。 从以上分析可知，经典的傅氏分析是一种纯频域 分析。
傅里叶变换是时域到频域互相转化的工具，从物理 意义上讲，傅里叶变换的实质是把f(t)这个波形分 解成许多不同频率的正弦波的叠加和。 从傅里叶变换中可以看出，这些标准基是由正弦波 及其高次谐波组成的，因此它在频域内是局部化的。
短时傅里叶变换定义如下：
1 Fg f ( , ) 2



f (t ) g (t )e it dt
其中，f(t)是待分析的信号； g () 是 g () 的复共轭函数； 函数 g(t)是固定的紧支集函数，称为窗口函数。 随着时间τ的变化，g(t)所确定的“时间窗”在t轴上移 动，使f(t)“逐渐”进行分析。
第六章 基于小波变换的故障诊断方法
小波变换的基本原理 奇异性的检测
基于小波变换的原油管道泄漏检测
一、小波变换的基本原理
小波的由来
小波变换是由法国理论物理学家Grossmann与法国 数学家Morlet共同提出的。 小波分析是近20多年来发展起来的新兴学科，其 基础是平移和伸缩下的不变性，这使得能将一个信 号分解成对空间和尺度的独立贡献，同时又不丢失 原有信号的信息。
在实际的信号处理过程中，尤其是对非平稳信号的处 理中，信号在任一时刻附近的频域特征都很重要。
如在故障诊断中，故障点（机械故障、控制系统故障、 电力系统故障等）一般都对应于测试信号的突变点。 对于这些时变信号进行分析，通常需要提取某一时间 段（或瞬间）的频率信息或某一频率段所对应的时间 信息。 因此，需要寻求一种具有一定的时间和频率分辨率的 基函数来分析时变信号。
2
f(x，y) ∈ L2(R2)的傅里叶变换‖f(x，y)‖定义为：
ˆ ( , ) f x y




f ( x, y)e
i ( x x y y )
dxdy
定义11： 设f(t)为在R上定义的函数，我们称集合{t f (t ) 0} 为函数f(t)的支集（即f(t) ≠0的点所构成的集合 的闭包）。 具有紧支集的函数就是在有限区间外恒等于零的函 数。
2
若f(x)，g(x) ∈L2(R)，α，β 为常数，则αf(x) ＋βg(x) ∈L2(R)。因此，L2(R)构成了一个线 性空间。我们称其为平方可积函数空间。
定义2： 在L2(R)空间中的内积<f，g>定义为：
f , g f ( x) g ( x)dx


g 其中， (x ) 表示g(x)的共扼。
其中，s为尺度因子（scale factor），或简称为 尺度。
定义8：
把希尔伯特空间（Hilbert space）中的可测的、 平方可积的两维函数构成的子空间记作：L2(R2)。
定义9：
函数f(x，y) ∈ L2(R2)的经典范数‖f(x，y)‖定义为：

f
定义10：
2




f ( x, y) dxdy

f g ( x) f (u) g ( x u)du

定义6：
函数f(x)的傅里叶变换 fˆ ( ) 定义为：

ˆ f () f ( x)eix dx


定义7： 对任意函数f(x)，其扩张函数fs(x)定义为：
1 x f s ( x) f ( ) s s
为了研究信号的局部特征，科学家们提出了一些对傅 里叶变换进行改进的算法，其中短时傅里叶变换 （Short Time Fourier Transform－STFT）就是 比较有代表性的一种。 短时傅里叶变换是一种折衷的信号时、频信息分析方 法，它是Dennis Gabor于1946年提出的。
短时傅里叶变换的基本思想是：通过给信号加一个小 窗，将信号划分为许多小的时间间隔，用傅里叶变换 来对每一个时间间隔内的信号进行分析，以便确定该 时间间隔内的频率信息。 它假定非平稳信号在分析窗函数g(t)的这个短时间间 隔内是平稳的（伪平稳），并移动分析窗函数，使 f(t)g(t- τ )在不同的有限时间宽度内是平稳信号， 从而计算出各个不同时刻的功率谱。
1 f (t ) 2
其中，权函数：



g ( )eit d
1 g ( ) 2



f (t )e it dt
就是原来函数f(t)的傅里叶变换。
经过以上的变换，就将对 f (t ) g ( ) F 1{g ( )}
的研究，转化为对权系数，即其傅氏变换
1988年，Mallat与Meyer合作提出了多分辨分析 的框架。 1988年，Daubechies构造了具有有限支集的正 交小波基。在美国Pure＆Appl.Math.发表一篇长达 87页的论文，被公认是小波分析的经典文献。
1989年，Mallat在多分辨率分析基础上，构造了 Mallat算法。为此，Mallat于1989年荣获IEEE论文 奖。 1990年，崔锦泰和王建忠构造了基于样条函数的 单正交小波函数。