第8讲小波变换导论1.docVIP

第8章 小波变换导论 §1：连续小波变换（Continuous wavelet tramsform）. 定义：设x(t)是平方可积函数（记作），是被称为基本小波或母小波的函数，则 称为x(t)的小波变换，式中a0是尺度因子，b是位移，，其中 说明： （1）基本小波可能是复信号或实信号，在复信号时，一般是解析信号，例如 ，称为Morlet小波，它是高斯包络下的复指数函数。 （2）尺度因子a的作用是将基本小波作伸缩，a越大，愈宽，a越小，越窄，与角频率(等价。改变a，改变小波变换的分析区间。 （3）前加因子的目的是使不同a值下的能量保持相等。 （4）注意，内积和卷积之关系： （或，b是卷积后变量）。 （5）注意到，的付里叶变换为 小波变换的等效频域表示： 证明： 由 故： 因而： 上式左边即是WTx(a,b)，与右式是一对付里叶变换对，由付氏化反变换定义得： 说明： 1）如果是幅频特性比较集中的带通函数，则小波变换便具有表征待分析信号x(()频域上局部性质能力，改变a的值，可以分析x(()不同频域区间的性质。 设在(0为中心的窄带内，，则中心在。 2)采用不同a值作处理时，各的中心频率和带宽都不相同，但品质因数Ｑ不变（Ｑ＝中心频率／带宽）。 例： 当a=1时，中心频率w=w0，带宽。 当a=2时， 中心频率，带宽。 考虑单边带的一般情况：a=1时，：中心w0，带宽Δ(w a任意时::中心,带宽。 时频特性: a增加，(降低，展宽，时域窗变宽，观测更长时间区间，变窄，观测更窄的频域窗，且频率中心向低频移动。对于低频部分，对时间看得粗些，对频域看得仔细一些。 a减小，(增大，变窄，时域窗变窄，观测更短时间区间。变宽，观测更宽的频域窗，且中心频率向高频方向移动。对于高频部分，对时间域看得细一些，对频域看得粗些。 这是数学显微镜的能力。由上分析，小波变换的时频平面分割如图： 比较短时付里叶变换的时频特性 它相当于对x(t)施加一个滑动窗后，再做付里叶变换，这里， 。 例：取 一般情况下，STFT的时频平面： 它的时频窗大小是固定的。 不定性原理（测不准原理） 对于单位能量信号f(t),它的付里叶变换F(w)，设它的时间和频率域中心在原点，即： 定义它的时间和频率域宽度为： 则存在：，只有当时,上式为等号. 测不准原理告诉我们,不可能在时域和频域都获得任意的观察精度,要使得频率分辨率提高,必然牺牲时域分辨率,小波变换可以根据分析的需要,协调这对矛盾,它在低频区看到更多的频率细节,在高频区看到更多的时间细节。测不准原理名称来自量子力学中的海森堡的测不准原理。 CWT的一些性质： （1）叠加性 （2）时移性质： （3）尺度转换 此性质表明，当信号x(t)作某一倍数伸缩时，其小波变换将在a,b两轴上作同一比例伸缩，但不发生矢真变形。 （4）小波变换的内积定理 设 则 式中 注意，两变量内积写成： （＃） 证明：由巴塞瓦尔定理： 故： 又从的付氏变换性质得： 将上面4式带入（＃），并应用 （＃）＝ 由于 故：（＃）＝ # 在上面证明中，要求 这个条件称为容许条件，推论：和，是振荡的。 小波反变换 在小波内积定理中，取x1(t)=x(t), 得： 小波能量公式： 取x1(t)=x(t),x2(t)=x(t) 正规性条件（regularity condition）. 为了在频域上有较好的局域性，要求随a的减小而迅速减小， 即：若，或相当于在w=0处有高阶零点， 那么：，且对于 正则性条件说明：如果选择的母小波满足如上的正则性条件，在x(t)的多项式展开中，各项在小波变换中没有贡献，这突出信号的高阶起伏和高阶导数中可能存在的奇点，即小波变换反映信号中的高阶变化。 再生核方程 x(t)的小波变换满足如下式： 上式中： 为再生核，再生核方程说明两点问题： ①CWT是冗余的，处CWT取值可由其它值通过再生核构造出。 ②不是任意两维函数都可能是一个函数的小波变换，小波变换必须满足再生核方程约束。 几个常见小波 （1）Morlet小波 不严格满足允许性条件,近似满足 （2）Marr小波 高斯函数二阶导数 在原点有两阶零点，用于边缘检测。 （3）Harr小波 满足正交条件： 且 在原点仅有一阶零点。 （4）样条小波族 一次样条 ， 由再构造小波。 （5）Daubechies小波族 小波族由满足一定条件的滤波器，迭代逼近一个小波。 §2尺度和位移离散化的小波变换 连续小波存在信息冗余，希望只计算离散的位移和尺度下的小波变换值。 尺度离散化：取a0，尺度因子a只取a0的整数幂，例如：a仅取： 位移离散化：尺度取时，取b=b0，各位移为k·b0。在时，相应取。 在这些离散位置的小波伸缩平移系