《第8章小波分析》由会员分享，可在线阅读，更多相关《第8章小波分析（23页珍藏版）》请在金锄头文库上搜索。
1、第8章小波分析8.1从傅里叶变换到小波分析标准傅里叶分析存在以下主要问题：(1) 无时域局部化特性。为了求得傅里叶系数，理论上必须知道时域的全部信息， 即傅里叶分析是对信号的总体的统计分析。反过来讲，时域上任一时刻信号的变化都将波及 整个频谱图，从而导致傅里叶分析在时域上无任何局部化特性。这不利于非平稳信号、特别 是瞬变信号的分析，在频域上是看不出这些突变特征的。(2) 不能实现时频分析。信号分解转换到频域后，丢失掉了时域的信息，频域中某频 率或频带内的信息和时域中某时刻或时宽内的信息没有直接的对应关系，即不能给出某一指 定频带内的时域图形。这种对应关系称为时频分析，所以傅里叶分析不能进行时频
2、分析，而 时频分析在工程中却相当有用。由于标准傅里叶分析在工程应用中存在不足，所以需对其进行修正，Gabor在1964年 引入了短时傅里叶变换。其基本思想是：把信号划分成许多小的时间间隔，用傅里叶变换分 析每一个时间间隔，以便确定该时间间隔存在的频率。其表达式为：S (具)=j f (t)g (t -T )e-&t(8-1)R式中，g(t)是紧支集的函数，及定义域(非零域)有限的函数；f(t)是被分析的信号。在这个变换中，ej起着频限的作用，g (t)起着时限的作用。随着时间T的变化，g (t)所确定的时间窗在t轴上移动，使f (t)被逐段进行分析。g(t)也被成为窗口函数，S0,T )大致反
3、映了 f (t)在时刻T时频率为o的信号成分的相对含量。短时傅里叶变换虽然在一定程度上克服了标准傅里叶变换不具有局部分析能力的缺陷，但它也存在着自身的不足，分析窗 的大小和形状是恒定的，即短时傅里叶变换实质上是具有单一分辨率的分析；若要改变分辨率，则必须重新选择窗函数g(t)。因此，短时傅里叶变换适合用于分析平稳信号，但对于 非平稳信号，在信号波形变化剧烈时，主频是高频，要求有较高的时间分辨率；而信号波形 变化平缓时，主频是低频，则要求有较高的频率分辨率，显然短时傅里叶变换不能同时兼顾 二者。图8-1所示为短时傅里叶变换时频分析窗。图8-1短时傅里叶变换时频分析窗小波分析方法是一种窗口大小(即
4、窗口面积)固定，但其形状可以改变，时间窗和频 率窗都可以改变的时频局部化分析方法，即在低频部分具有较高的频率分辨率和较低的时间 分辨率，在高频部分具有较高的时间分辨率和较低的频率分辨率，因此其对信号具有自适应 性。小波分析方法的出现可以追溯到1910年Haar提出Haar规范正交基，以及1938年 Littlewood-Paley对傅里叶级数建立的L-P理论。为克服传统傅里叶分析的不足，在上世 纪八十年代初，便有科学家使用“小波”的概念来进行数据处理，比较著名的是1984年法 国地球物理学家Morlet引入小波的概念对石油勘探中的地震信号进行存贮和表示。在数学 方面所做的探索主要是R. Coi
5、fman和G. Weiss创立的“原子”和“分子”学说，这些“原 子”和“分子”构成了不同函数空间的基的组成部分。L. Carleron使用了非常象“小波” 的函数构造了 Stein和Weiss的空间H 1的无条件基。直到1986年，法国数学家Meyer成功地构造出了具有一定衰减性的光滑函数W，它的二进伸缩与平移 0= 2-j/2wG-jt k)： j,k e Z成L(R) 的规范正交基。此前，人们普遍认为这是 j, k不可能的，如Daubechies,Grossman和Meyer都退而研究函数系a |-j/2W C -j kb。)构成L(R)的框架的条件去了。Lemarie和Battle继M
6、eyer之后也分别独立地给出了具有指数衰减的小波函数。1987 年，Mallat利用多分辨分析的概念，统一了这之前的各种具体小波的构造，并提出了现今 广泛应用的Mallat快速小波分解和重构算法。1988年Daubechies构造了具有紧支集的正 交小波基。Coifman, Meyer等人在1989年引入了小波包的概念。基于样条函数的单正交小 波基由崔锦泰和王建忠在1990年构造出来。1992年A. Cohen, I. Daubechhies等人构造 出了紧支撑双正交小波基。同一时期，有关小波变换与滤波器组之间的关系也得到了深入研 究。小波分析的理论基础基本建立起来。近年来，一种简明有效的构造
7、小波基的方法-提升方案(Lifting Scheme)得到很大的发 展和重视。利用提升方案可把现存的所有紧支撑小波分解成更为基本的步骤，另外，它还为 构造北线性小波提供了一种有力的手段，所以，利用提升方案构造的小波被认为是第二代小 波。小波理论及其应用仍然处在发展中，其未来将在北线性多尺度方法、非规则集上的小波 构造以及北平稳、非均匀、时变信号处理等方面等到更深入的研究。8.2小波变换的基本概念8.2.1小波的定义小波函数的定义：设W (t)为平方可积函数，即W (t) e L2(R)，若其傅里叶变换W ()W（w）wl满足(8-2)则称W（t）为一个基本小波，又称为小波基（函数）、小波核（函
8、数）或母小波。由于基本小波w （t）生成的小波W ab（t）在小波变换中对被分析的信号起着观测窗的作用，所以W （t）还应该满足一般函数的约束条件尸珈（t） ）t 8（8-3）s故W皿）是一个连续函数。这意味着，为了满足容许条件式（8-3），W皿）在3= 0时必须等于零，即W （0）=+8W （t ）dt = 0（ 8-4）一8上式说明，W（t）围绕时间轴的面积必须为零，故叩（t）必须是一个振荡波形；同时，又希望有局部化的时窗，因此，W（t）应该选用快速衰减的短小波形。根据这两点，W（t）被称为小波。将基本小波W （t）经伸缩和平移后，就可以得到一个小波序列。小波序列为W ab(t)=(a,b
9、 g R, a 丰 0)(8-5)式中：a-尺度（缩放）因子b -平移因子-归一化因子，使不同尺度的小波保持相同的能量。板a小波变换的基本思想来源于伸缩和平移方法。（1）尺度伸缩（Scaling）对波形的尺度伸缩就是在时间轴上对信号进行压缩与伸展，时间尺度a是反比于 频率。尺度a越大，意味着小波函数在时间上越长，即被分析的信号区间越长。因此，尺度 越大意味着对信号的较长的时间段进行近似比较，主要获取的是信号的低频特性。反之，尺 度越小，意味着只对信号的较小的细节部进行比较，因此主要获得的是信号的高频特性。（2）时间平移（Shifting）是指小波函数在时间轴上的波形平行移动。对于离散伸缩、平移
10、的情况，小波序列为W k(t) = 2一2W(2-jt-k)(j,k gZ)(8-6)对于任意函数f (t) g L2(R)的连续小波变换定义为W(a,b) = (f，Wa b)= |a -1/2 J f (t)W*f-dt(8-7)其中W *是W的共轭函数。小波逆变换为/ “)=七*，伽dadb(8-8)图8-2小波变换的时频分析窗小波变换的时频窗口与短时傅里叶的时频窗口不一样，其窗口形状如图8-2所示为两 个矩形b-aAw,b + aAwx(气-AW)/a,(% +A/)/a，窗口中心为(b,0 /a)，时窗 和频宽分别为a Aw和AV / a。其中Aw是母小波函数W (t)的半时宽，A#
11、是W (o)的半频宽；b只影响窗口在相平面时间轴上的位置；而a不仅影响窗口在频率轴上的位置，也影 响窗口的形状。这样小波变换对不同的频率在时域上的取样步长是调节性的，即在低频时小 波变换的时间分辨率差，而频率分辨率较高。而在高频时相反。这正好符合低频信号变化缓 慢高频信号变化迅速的特点。这也是小波变换优于标准傅里叶变换和短时傅里叶变换的地 方。8.2.2常用小波函数与标准傅里叶变换相比，小波分析中所用到的小波函数具有不唯一性，即小波函数W (t)具有多样性。小波分析在工程应用中一个十分重要的问题就是最优小波基的选择问题，针对 同一分析对象，选用不同的小波基其分析结果也不同。目前主要是通过小波分
12、析方法处理信 号的结果与理论结果的误差来判断小波基的好坏，并确定小波基。下面介绍一些常用的小波 函数。(1) 、Haar 小波Haar函数是在小波分析中最早用到的一个具有紧支集的正交小波函数，同时也是最简 单的一个函数，如图8-3所示。Haar函数的定义为l.(0 x 1/ 2)=-1.(1/2 x 1)(8-9)H0.(其他)图8-3 Haar函数(2) 、Daubechies (dbN)小波系对于正交小波，希望它是有限支撑的，以便使小波变换的快速算法更快捷；希望它是光 滑的，以便高精度地模拟和分析信号；希望它是时域和频域的局部化能力是强劲的，以便在 信号分析处理中发挥突出的作用Daubec
13、hies小波为此做出杰出贡献，其中一阶Daubechies 小波就是上面介绍的Haar小波。Daubechies系中的小波基记为dbN，N为序号，且N=1，2, 3,，10。dbN大多数不具有对称性，正则性随着N的增加而增加，函数具有正交性。除 了一阶Daubechies小波，其他Daubechies小波没有明确的表达式，但转换函数H()的 平方模式明确的，即H0 0)|2 =(MCOS2V2 7(一 5P S1H2V 2 7(8-10)式中H (Q=q1 h e- jg0；2 k=okP ( y )=N Cn-1+kyk kk =0其中C E为二次项系数。图8-4为4种Daubechies小
14、波。（3）、Morlet （Morl）小波Morlet函数（图8-5）定义为W （x） = Ce-x2/2 cos w x0它的尺度函数不存在，故不具有正交性。一0.5-0.2图8-6墨西哥帽函数图8-5 Morlet函数（4）、Mexican Hat （Mexh）小波Mexican （图8-6）函数定义为兀-1/4(1一 x2)e-x22(8-11)（8-12）它是Gauss函数的二阶导数，因为其形状像墨西哥帽的截图，所以称它为墨西哥帽函数。 它在时域和频域都有很好的局部化能力，并且满足：j f （x）dx = 0-s由于它的尺度函数不存在，所以不具有正交性。表8-1列出了几种常用小波的主要性质对比。表8-1种常用小波的主要性质对比小波 函数HaarDaubec hiesMorletMexica n HatBiorthogonalCoifle tSymletMeyer正交性有有无无无有有有双正 交性有有无无有有有有紧支 撑集有有无无有有有无连续 小波 变换可以可以可以可以可以可以可以可以离散 小波 变换可以可以不可以不可以可以可以可以可以，但无FWT支撑 长度12N-1有限长度有限长度重构：2N +16N-12N-1有限长度分解：2Nd +1滤波器长度22N-4,45,5max(2N ,2 N )+ 2 6N2N-8,8对程 性对称近似对 称对称对称不对称近似对 称
