第六章基于小波变换的故障诊断方法.ppt

第六章 基于小波变换的故障诊断方法 小波变换的基本原理 奇异性的检测 基于小波变换的原油管道泄漏检测 设以Vj表示小波分解树结构图中的低频部分Aj，Wj表示分解中的高频部分Dj，则多分辨分析的子空间V0是一个有限个子空间的逼近，即： 上式中的空间列 具有如下性质： ① ② ③ 设法找出一个确定的函数 ，使得对每个 ，函数系 构成空间Wj的规范正交基， 其中， （1） 若令 代表分辨率为2-j的函数 的逼近（即函数f的低频部分或“粗糙像”），而 代表逼近的误差（即函数f的高频部分或 “细节”部分），则下式 意味着： 注意到f=f0，所以上式可简写为： （2） （3） 上式（3）表明，任何函数 都可以根据分辨率为2-N时f(t)的低频部分（“粗糙像”）和分辨率2-j（1≤j≤N）下f(t)的高频部分（“细节”部分）完全重构，这也从另一侧面表示Mallat塔式重构算法。 （3） 从包容关系 ，可很容易得到尺度函数ф(t)的一个极为有用的性质。 注意到 ，所以 可以用V-1子空间的基函数 展开，令展开系数为hk，则 这就是尺度函数的双尺度方程。 （4） 另一方面，由于 ，故 这就意味着小波基函数 可以用V-1子空间的正交基 展开，令展开系数为gk，即有 这就是小波函数的双尺度方程。 （5） 由双尺度方程式（4）和式（5）可知，尺度函数与小波基函数的构造归结为系数 的设计。 小波基函数 可由尺度函数 的平移和伸缩的线性组合获得。 若令 则把尺度函数和小波基函数的设计可以归结为滤波器H(ω)，G(ω)的设计。 （4） （5） 尺度函数的双尺度方程: 小波函数的双尺度方程: 构造正交小波时，滤波器H(ω)，G(ω)必须满足以下三个条件： （6） （7） （8） 联合求解式（6）、（7）、（8）可得 （9） 所以，要设计正交小波，只需要设计滤波器H(ω)。 综上分析，为了使 构成Vj子空间的正交基， 应该具有下列基本性质： ① 尺度函数的容许条件， 。 ② 能量归一化条件： 。 ③ 尺度函数 具有正交性，即 ⑤ 跨尺度的尺度函数 和 相关。 ④ 尺度函数 与基小波函数 正交，即 ⑥ 基小波函数 和 相关。 将尺度函数的容许条件与小波的容许条件 作比较可知，尺度函数的傅里叶变换 具有低通滤波特性（相当于一个低通滤波器），而小波基函数的傅里叶变换 则具有高通滤波特性（相当于一个带通滤波器）。 多分辨分析可以对信号进行有效的时频分解，它不同于短时傅里叶变换对信号频带的等间隔划分，它的尺度是按二进制划分的，所以在高频频段其频率分辨率较差，而在低频频段其时间分辨率较差，即对信号的频带进行指数等间隔划分（具有等Q结构）。 2、一维Mallat算法 S.Mallat在Burt and Adelson图像分解和重构的塔式算法的启发下，基于多分辨率框架，提出了一种具有完美数学描述的塔式多分辨率离散小波分解与重构算法，即Mallat算法。 Mallat算法的本质是：不需要知道小波函数和尺度函数的具体结构，仅根据系数就可实现信号的分解和重构。 而且这种算法可以使每次小波分解后信号长度减半，大大减少了小波变换的复杂度，因而它是一种快速算法。 Mallat快速算法在小波变换中的地位，与快速傅里叶变换在傅里叶变换中的地位相当。 设Φ(t)为尺度函数，Ψ(t)为基本小波，设信号f(t) ∈L2(R)，将信号进行分解，用A表示低频，D表示高频，下标数字表示小波分解的尺度层数，已得到f(t)在2-j分辨率下的粗糙像Ajf ∈Vj，｛Vj｝j ∈Z构成L2（R）的多分辨分析，从而有： 式中， （10） 式（10）可写成如下形式： 由尺度函数的双尺度方程可得： 利用尺度函数的正交性，有： （12） （13） 同理，由小波函数的双尺度方程可得： （14） （11） 由式（11）～（14），可得： （15） （16） （17） 引入无穷矩阵 ， 其中