小波分析理论简介

1、小波分析理论简介(一) 傅立叶变换伟大的历史贡献及其局限性1 Fourier 变换1807年，由当年随拿破仑远征埃及的法国数学、物理学家傅立叶( JeanBaptistle Joseph Fourier ,1786-1830),提出任意一个周期为 T (= 2 )的函数f(t)，都可以用三角级数表示：f(t) = Ckeiktk+ ak cos kt +2k 1bk sin ktk 1(1)CkL f(t)eiktdt= f,eiktak Ck C kbk i(Ck C k )对于离散的时程f(t)，即N个离散的测点值m 0,1 , 2,N1,T为测量时间:其中akbkCkN2 1(ak co

2、s kt bk sink 1N 12 kmXm Cos m 0N2 N 12 kmxm sin N m0NN 1Ii(2 km/N)xmeN m 0时，化为傅立叶积分(即Fourierf( ) f(t)ei tdtf,e、1kt)+2a0,0,(8)变换)：n cos21,1,N 1Nt= Ckei kt 22,-1,2,(4)(5)2,N-1(9)N -121 . tf2-f()eitd傅立叶变换的理论是人类数学发展史上的一个里程碑，从1807年开始，直到1966年（1807年傅立叶提出任意一个周期函数都可以表示为傅立叶级数的 结论是有误的，直到1966年才证明了 L2可积的周期函数才能表示

3、为傅立叶级 数），整整用了一个半世纪多，才发展成熟。她在各个领域产生了深刻的影响， 得到了广泛的应用，推动了人类文明的发展。其原因是，傅立叶理论不仅仅在数学上有很大的理论价值，更重要的是傅立叶变换或傅立叶积分得到的频谱信 息具有物理意义。所以说，傅立叶理论是万古流芳的。数学上的插值方法。除傅立叶级数外，还有拉格朗日插值，有限元插值，勒让德多项式插值即高斯积分使用的插值方法。遗憾的是，这种理论具有一定的局限性：（1）傅立叶变换的三种形式中的 傅立叶系数都是常数，不随时间t变 化，因而只能处理频谱成分不变的平稳信号，相反的，在处理非平稳信号时会带来很大误差，甚至与实际情况大相径庭。（举例：无阻尼与

4、有阻尼的单自由度 的自由振动、打秋千、座钟、讨论会与大合唱等）。在实际信号中，若高频与低频差别很大，在相同的时间间隔内，高频信号 衰减了而低频信号尚未衰减，所以，在不同时刻，信号的频谱成分是不同的。 硬要用傅立叶变换找出所有时刻的频谱成分，硬要把幅值的变化用频率的变化 来补偿，不仅高频的傅立叶系数有误差，低频的傅立叶系数也有很大误差，包 括求出的频率当然也有误差。（2） 求傅立叶系数是全时间域上的加权平均，这从上面的（5）、（6）、（7）公式可以清楚看到。局部突变信息被平均掉了，局部突变信息的作用很难 反映出来（好比吃大锅饭，平均主义）。差别很大的信号，如方波、三角波、正 弦波，都可以得到相同

5、的频率，所以，处理、捕捉突变信号如故障信号，灵敏度很差。处理、捕捉突变信号应使用能反映局部信息的变换。为了克服以上两点局限性，这就要求：(1) 将变换系数视为随时间变化的，级数求和由一重变为两重。(2) 使用能反映局部信息的变换，贝恼数组不能使用全域上的函数，只能使用有所谓紧支撑的函数，即“小波函数”或加窗傅立叶变换的窗函数。2 Garbor 变换一窗口 Fourier 变换在时间一频率分析中，Fourier 变换公式的不足已经被D. Garbor注意 到了，在1946年的论文中，为了提取信号的Fourier变换的局部信息，引入 了一个时间局部化的Gaussian函数作为“窗函数” g(t-b

6、),其中参数b用于 平移动窗以便覆盖整个时间域。因为一个Gaussian函数的Fourier变换还是Gaussian函数，所以Fourier逆变换即频率也是局部的。窗口 Fourier 变换简介。对于时间局部化的“最优”窗，用任一 Gaussian函数ga(t)t24a(11)Garbor变换”的定义为(G；f)( ) (ei tf(t)ga(t b)dt(由于ga(t b)dbga(x)dx 1(13)所以 (e i tf(t)ga(t b)dt db = f( )(14)令Gb, (t) = ei tga(t b)(15)利用Parseval 恒等式，(G：f)( ) (eitf(t)ga

7、(t b)dt = fG：, .=* f,G；：f)( b)，畀 2L (e" f( )g14a()d(这个等式说明，除去乘数项.e ib 之外，在t b具有窗函数ga的f的“窗口 Fourier 变换”,结论，这两个窗的宽度分别是2 . a 和1,a这两个窗的笛卡儿积是b . a ,b a 112”a '2“ a与在具有窗函数g1的f的“窗口 Fourier逆变换”一致，根据窗函数 ga的宽度是2 . a的加窗傅立叶变换的“时间一频率窗”的宽度对于观察所有的频率是不变的。 在较长的时间窗内，对于高频信号，可能经过了很多周期，因而求出的Fourier 变换系数是很多周期的平均

8、值，局部化性能不能得到体现。若减小时间窗(减 小a)，高频信号局部化性能得到体现，但对于很低的频率信号来讲，检测不到。 总上所述，加窗傅立叶变换对于高频与低频差别很大的信号仍不是很有效的。3 窗口 Fourier 变换的测不准原理对于一个非平凡函数 wL2(IR)，若满足2tw L (IR)(A)条件，则 w可作为短时窗口Fourier变换的窗函数，若其Fourier变换也满足上述条件，那么(B)而且，等号成立，如且仅如w(t)ce%a(t b)(C)其中 c 0,a 0和a,bIR小结:(1) 傅立叶级数的正弦与余弦系数为常数，不能反映振幅变化的情况；(2) 求傅立叶系数需要所考虑的时间域上

9、所有信息， 不能反映局部信息的特征；（3）加窗傅立叶变换时间窗是固定不变的，高频与低频的时间局部化不 能同时满足。由于上述原因，必须进一步改进，克服上述不足，这就导致了小波分析。（二）小波分析将时程函数f（t）表示为下面的小波级数：f(t)jf - j,kj,k(t)=kjdj,k j,k(t)( 17)kj,k(2jt k)（18）其中，（t）是小波函数，dj,k是小波系数，且dj,k =仁 j,k(19)由公式（17）到（19）可以看到，小波级数是两重求和， 小波系数的指标 不仅有频率的指标j，而且还有时间的指标 k。也就是说，小波系数不仅像傅 立叶系数那样，是随频率不同而变化的，而且对于

10、同一个频率指标j，在不同时刻k，小波系数也是不同的。这样就克服了上面所述的第一个不足。由于小波函数具有紧支撑的性质，即某一区间外为零。 这样在求各频 率水平不同时刻的小波系数时， 只用到该时刻附近的局部信息， 从而克服 了上面所述的第二个不足。与有限元比较。在这一点，小波插值要比有限元高明。 有限元虽然是局部的 "单元插值”，但单元之间的公共节点上, 只能保证C0阶连续，而导数不连续。小波插值可保证二阶导数连续，只要选三次样条小波就能做到。第三个不足，小波分析是如何克服的呢？通过与加窗傅立叶变换的“时间一频率窗”的相似分析，可得到小波变换的“时间一频率窗”的笛卡儿积是? 1 ? 1

11、b at a ,b at a ,（ 20） 其中a 2 j，时间窗的宽度为2a ，随着频率的增大（即j的增大）而变窄， 随着频率的减小（即j的减小）而变宽，之所以有这样的结果，关键在于公式（18）中，时间变量t前面乘了个“膨胀系数” 2j。小波变换的“时间一频率窗”的宽度，检测高频信号时变窄，检测低频信 号时变宽，这正是时间一频率分析所希望的。根据小波变换的“时间一频率窗”的宽度可变的特点，为了克服上面所述的第三个不足，只要不同时检测高频与低频信息，问题就迎刃而解了。女口，选择从咼频到低频的检测次序，首先选择最窄的时间窗，检测到最咼频率信息， 并将其分离。然后，适当放宽时间窗，再检测剩余信息中

12、的次高频信息。再分 离，再放宽时间窗，再检测次次高频信息，依次类推。为了检测到不同频率水平信息，即求出不同频率水平下不同时刻的小波系 数，首先要选好小波函数。选择小波函数的“四项原则”。在求小波系数公式（19）中，如果 （t k）是L2（IR）空间的正交基，则的 j,k为j,k的复共轭。小波分析的最重要的应用是 滤波，为了保证滤波不失真， 小波函数必须具有 线性相位，至少具有广义线性相位。小 波分析的另一重要应 用是捕捉、分析突变信号，这就要使用函数的 导数，小波函数至少是C1连续。 由前面分析可知，小波函数必须具有紧支撑的性质。所以，正交、线性相位、 连续、紧支撑是选择小波函数的“四项原则”

13、。如果选择某个小波函数，同时满足四项指标，那真是人类的福气。遗憾的是，上帝像是有意考验我们的数学家，没有将“四合一”的小波函数“直接”恩赐给人类。数学家们已经证明，具有正交、线性相位、紧支撑的小波函数只有Harr函数，而Harr函数是间断函数,对于工程应用来说，是不理想的。目前，一种倾向是坚持正交性。另一种倾向是放弃正交性，另辟途径，进行艰辛的长征，前仆后继，花费了将近半个世纪的探索，才使小波分析理论成熟起来，得以在工程中应用。作为后人，我们要忠心 地感谢他们。为了进行小波分解与重构，“四合一”的小波函数不存在，数学家们“一分 为四”，选择了四个函数，巧妙地解决了这些问题。这四个函数是：尺度函

14、数 ， 小波函数，对偶尺度函数 ，对偶小波函数 。为什么要选择四个函数呢？由前面小波变换的“时间一频率窗”分析可知，小波变换的“时间一频率 窗”的宽度，当检测高频信号时变窄，检测低频信号时变宽。为了检测到所有 频率信号，“时间一频率窗”的宽度必须按一定的次序变化，不失一般性，从窄 到宽，检测频率信号从高频到低频的次序进行一一实际上也正是这样的次序。在最高频率水平 Vn （即根据实测数据的时间测量间隔t，最高能检测到的频率为Nyquist频率f 貴），选择最窄的“时一频窗”宽度，检测到原 始信号中的最高频率信号，并将这些信号从原始信号中剥离，存放在Wni空间,而将剥离后的剩余低频信号的总合，存放

15、在另一空间Vni。然后，增大“时一频窗”的宽度，再检测Vn i空间中的高频信息，将这些信号从Vn i空间中剥离， 存放在Wn 2空间，而将剥离后的剩余低频信号的总合，存放在另一空间Vn 2。依次类推。这就要求有两个互相有联系的空间：Vj Vj 1 Wj 1=Vj 2 Wj 2 Wj 1W 1 Wo Wi Wj 2 Wj 1 ,J ZZ(1)V1V0V1V2(2)clos l2(J2Vj) = L (IR)ZZ(3)VjJ ZZ0 ,即 lim j/j0(4)Vj VJ 1Wj 1(5)f(x)Vjf (2x)Vj 1(6)WJ w0 , J丨子空间性质简介:这样，对于参考子空间V，需要单个函数

16、l2（ir）在意义Vo=cI0Sl2(|r)0,k：K ZZ(22)上生成，其中，j,k22 (2jx k)(23)(26)(27)(28)(29)（26）称为尺度函数的两尺度关系”（27）称为小波函数的“两尺度关对于参考子空间W0，需要单个函数L2（IR）在意义W0 = clos l2（ir）； 0,k : K ZZ（24）上生成，其中，j,k 22 （2jx k）（ 25）首先，这就要求有两个函数：和，前者称为尺度函数，后者称为小波函数。并且，它们肯定是有关系的。由（21）中的（1）式可知，V V，由（4）式可知，W0 V，而 1,k : k ZZ是V1的一个基，所以存在唯一的I2序列 P

17、k、qk(x)Pk (2x k)k(x)qk (2x k)k并引入记号P(z) = 1PkZk2 kQ(z) = 1qkZk2 k系” P（z）称为尺度函数的“两尺度符号” Q（z）称为小波函数的“两尺度符为了由高频到低频逐次检测到不同频率水平的信息，仅有上述两个函数是不够的。由前面分析可知，公式（17）中的小波函数与求小波系数使用的与 f作内积的函数不是同一个函数，除非使用正交的小波函数。这就要求寻找尺度函 数 与小波函数 的对偶函数：对偶尺度函数 、对偶小波函数 ，以便分解 原始信号时能求出小波系数来。对偶关系简介与(26)( 29)相对应，有 ? (x) gk (2x k)(30)k ？

18、(x) hk (2x k)( 31)kG?(z) = -g?zk( 32)2 kH?(z)=丄hk?zk( 33)2 k为了满足对偶关系，必须满足以下对偶条件:Mp,Q(z) MTg,h(z) = 1(34)(35)与(34)、(35)等价的条件P(z) G(z)Q( z)H(z) 1P(z)G( z)Q(z)H( z)(P(z)G(z)P( z)G( z)-P(z)H(z)P( z)H( z)G(z)Q(z)G( z)Q( z)(Q(z)H(z)Q( z)H( z)P,Q ( z)=det M p,q(z)G,H (z)=det Mg,h(z)0100100Mtg,h(z) Mpjz) =

19、0疋其中,G(z)为G?(z)的共轭,(36)(37)(38)(39)H(z)为H ?(z)的共轭,M P,Q (z)=P(z) Q(z)P( z) Q( z)(40)Mg,h (Z)=G(z) H(z)G( Z) H( z)(41)由(36)和(38)便可以得到Vi的基到Vo和 Wo中的基分解关系:1(2x l) - g2k i (x k) h2k i (x k) (42)2 k由(42)式便可得到j 1频率水平子空间Vj j和Wj 1中向量坐标的分解算法ck1ail2kcij其中dkj 1bi 2kC|ji1an = g n2bn-h(43)(44)(45)(46)根据两尺度关系，便可得到

20、j频率水平子空间Vj中向量坐标的重构算法j 1j 1pk 2iCiqk 2idi(47)有了上述的系数，就可以使用多分辨分析的金字塔算法，快速求岀小波系数dj,k。以三次样条函数尺度函数N4为例，说明其步骤。1将f(X)投影到Vn上k) =fn(48)2n 1f (X)c： (2nxk 02小波分解算法使用多分辨分析的金字塔算法dn1dn2n md/n cn 1n 2ccn m cckai 22/i(49).j 1.jdkbi 2kGi(50)fn :=gn1gn 2g n 3 gnmf n m(51)fjCkjk(2jx k)(52)gjdkj 4(2jx k)(53)k1G(z)=-2 n

21、ngnZz1(12Z、m E2m 1 ( Z)E2m1(/)(54)口/、1nJZ m (2m1)!H (z)=hnZZ (-) £(55)2 n2E2m1(Z )2n 2E2m 1=(2m1)!N2n(kk1)Z(56)而gn、hn是符号多项式G（Z）、H（Z）的系数:k 03小波重构算法dn1 dn2/ / /nn 1n 2c c cdnj 1Pk 2lC|在我们研制的小波分析软件qk 2ld/“ XIAOBO ”中，使用的是三次样条函数尺度函数4 N4及三次样条小波（B-小波）4（X）。4（X）三次样条小波qnN4（2x n）n4（X）=(57)其中(1)n 4 43N8(n

22、12 i 0 lm阶样条函数Nm(x)的定义是1Nm(x) (Nm 1 N1)0Nm1(X t)dtN1是区间0， 1的特征函数。qnl），（n0,1,2,10)(58)（59）（60）三次样条小波具有连续、紧支撑、广义线性相位和半正交的性质。三次样条小波是C2阶连续，用它模拟一个信号，二阶导数都是连续的，精度极高。三次样条小波支撑区间为0，7，三次样条尺度函数支撑区间为0，4 三次样条小波是对称函数，因而具有广义线性相位，滤波不会失真。唯一不足的是三次样条小波不是正交的，只是半正交的，在小波分解中，求小波系数必须使用其对偶小波4，对偶小波4的傅立叶变换为4()?4()2(61)k?4(2 )

23、?4()?4()2(62)?4( 2k)k(1)<小波分类简介1 正交小波正交小波的两个族j,k , j,l > = k,lj ,k 与 j,k满足(2)< j,k，j,l> = 0(3)< j,k，l ,m > =j ,l k,m（4）L2(IR)W1W0W1用振动的语言形象的比喻：j,k是Vj j子空间的低阶振型，j,k是Vj i子空间的高阶振型， WjVj , WjWk, （j k）2半正交小波半正交小波的四个族j,k、j,k、j,k 与 j,k满足：(1) <j,k , l,m> =j,l k,m<j ,k , j,m > =

24、k,m(3)< j,k , l,m > . .=0, <j,k, l,m> = 0， jl;j,l,k,mZZ(4)j,k 是 Vj1子空间的低阶振型的线性组合，即低阶剩余模态。(5)j,k 是 Vj1子空间的高阶振型的线性组合,即咼阶剩余模态。(6)WjWj，VjVj，WjVj(7)L2(IR)W1W0W13 R小波R小波的 j,k与 j ,k满足：(1) <j,k, l,m > =j,l k,m(2)L2(IR)? ? ? ?W1W0W1正交性都是相对于对偶来说的，本身都不是正交族，只是线性无关族。神奇的小波分析为了检验小波分析的能力，在实测数据中叠加上

25、一个低频、幅值为原始最大幅值的十分之一的正弦 拍波，经过小波分解，在 27和26两个分辨率水平下都检测到了这个波，尤其是26检测到的正弦拍波曲线，如荷花"出污泥而不染”（三）正交小波包与双正交小波包分析从上面介绍的多分辨分析的金字塔算法进行小波分解可以看到，每次分解，都是对V进行的：对Vj进行小波分解，将 Vj分解为Vji和Wj 1 ；再对Vj 1进行小波分解，将 Vj 1分解为Vj 2和Wj 2，依次类推。而各子空间 Wj不再分解，也就是说，每次都是对低频进行再分解，而高 频不再分解！究其原因，便是 Vi的基向量（2x I）到V。和Wo基（x k）、 （x k）有分解关系(42),

26、而Wi的基向量(2x I)没有相应的分解关系。事实上，(2x I)与(x k)、 (x k)是线性无关的！显然，咼频的时间一频率局部化不是最优的。为了克复这个缺点，必须使 用小波包的分解方法。无论正交小波包还是双正交小波包，与小波的最基本的区别在于它们具有qn ( 1)np n 1( 63)的关系根据这个特殊的正交的两尺度关系，可以定义关于尺度函数“小波包”：2l (x):2l 1 (x):Pk l(2xkqk l(2xkk)k)(64)并且，小波包族具有下属正交性:n(x l), n(x k) > =l,k(65)2n (x l ), 2n 1(X k) > = 0(66)因此有

27、分解关系:u0 Vj , u；Wj(67)小波分解中的关系Vj 1VjWj(68)现在可以改写作并且，(69)式可以由 nUjn1n(2j1xU；1U0U1(69)0推广到任意ZZ :Uj2nU2nm)=Pm 2k2n(2jk)qm 2k 2n 1( 2 ' Xk)(70)(71)这样，对于正交小波的小波分解得到的l2(ir)W1(72)对于每个j 1, 2, 3,可以用小波包再分解：1 Uj21 Uj31U：2 U6WjWj U ；Uj422 Uj72Wj2jUouT2j 1 1U2 1双正交小波包与正交小波包类似，略去。（四）向量分解小波包简介上面介绍了正交小波包与双正交小波包分解

28、，与小波分解的最基本的区别 在于它们具有（63）式qn （1）np n 1因而导出了正交小波包和双正交小波包的分解。正如前面提到的，小波分解与小波包分解的共同点，都是高频空间的基向 量能向低频空间的基向量分解，若不能分解，便进行不下去。如小波分解中，Wj 的基向量（2jx I）没有相应的分解关系，所以Wj不能再分解。小波包克服了这种弊病，得到了各子空间都能适用的基向量分解关系式（71），才解决了这个问题。我们能不能跳出基向量能分解才能分解的框框呢？即，基向量不能分解时 也能分解？答案是肯定的。 这就是向量分解小波包 的理论。现在，让我们从小波重构关系式（47）CkPk 2|C|j qk 2|d

29、|j , k=1,2,3, ,MI出发，导出向量分解小波包的理论首先，将（47）改写为矩阵的形式：j 1P M ? c j +MM 2 12j 1 M 2(74)Q的那么，其中，矩阵P、Q由两尺度序列P、q中的数组成。我们把矩阵P 每一列，看成向量空间的一个向量，把c M 1也看成向量空间的一个向量,cM1就是这些矩阵P、Q的列向量的一个线性组合:22c M 1 =ck 1 p k +dkj 1 q k（75）k 1k 1问题是，这些矩阵P、Q的列向量是否线性无关？由分解等式（43）和（44） 可知，当cM1为零向量时，ckj1 0, dk 1 0 , k=1,2,3,，丛，这就得出了2 矩阵

30、P、Q的列向量是线性无关的结论。因此，可以把矩阵中的M个列向量看成M维向量空间的一组基。找到了 M维向量空间的一组基，就可以对任一 M维向量进行分解，找出在 这组基下的坐标。只要认准被分解的对象是 M维向量空间的一个向量，而不是再看做 V或 W中向量的坐标。这样，在小波分解中，不仅得到各频率水平的坐标 c可按原来的步骤和 方法进行分解，而且得到的各频率水平的坐标 d也可按向量分解的理论进行 再分解，并且分解的步骤和方法可以和对坐标 c进行再分解的完全一样。可能有人会产生疑问：前面讲到， V W = 0 ,现在W中的向量怎么又变成V中的向量呢？其实，这里使用了 “移花接木”的技巧。我们现在是对

31、W中的向量g分 解，这是千真万确的。但，是通过对其坐标 d的分解来实现的。也就是说， 要求的是d等于什么？而不是W中的向量g等于什么？W中的向量g等于g =? d（ 76）既然要求的是d等于什么，而不是g等于什么，不妨“移花接木” 一下，把 拿来，虚构V中一向量x，x = ? d（77）不就可以对x进行分解了吗？反正在分解过程中，使用的是d、a、b ，求的是新的c和d，而x根本就不出现，管它是什么样子呢！只要保证分解后低 频与高频分离即可。既然c分解，低频与高频能分离，那么，d分解，同样低 频与高频可以分离的。其实从（74）式中两个矩阵P、Q的组成，都是两尺 度序列P、q中的数，仔细分析一下，

32、不难看出，P中的列相当于低阶振型， Q中的列相当于高阶振型。以后的工作和小波包的分解类似，只不过这里是直 和分解，不是正交分解罢了。前面曾提到，小波分解的第一步，将 f（x）投影到Vn上（见公式（48）。 由于投影到Vn上求小波系数要解方程，很费时间，所以，好多书上都说，可 以近似地把原始数据作为小波系数， 误差不大。用向量分解小波包的理论来看， 这并非是近似的，而是精确的！事实上，使用上面“移花接木”的技巧，虚构 必中一向量xx = ? f（ 78）f是原始数据组成的向量。那末， f不就是必中的向量x的坐标C吗？只要 对它进行分解与重构就得了，而且，重构后，也用不着通过（78）式求x，重构得到的c，就是原始数据组成的向量 f ,两头省！当然，我并非反对投影，而且在捕捉突变奇异信号时，最灵敏的方法是利 用导数确定，投影到Vn上先求出作为的Vn中的向量 仁的坐标C，能得到快速 的求导方法。但这是另外的问题了。有了向量分解小波包新理论，可以在分解过程中，把各种小波组合到一起， 例如，求 fN时，可以用三次样条小波，这样确保滤波重构不失真、求导精 度高的长处，在分解的时候，可使用正交分解的分解序列系数，以提高分解效 率。这样，上帝虽然没有将“四合一”的小波