《小波分析理论简介》由会员分享，可在线阅读，更多相关《小波分析理论简介（24页珍藏版）》请在金锄头文库上搜索。
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。也就是说，小波系数不仅像傅 立叶系数那样，是随频率不同而变化的，而且对于同一个频率指标j，在不
10、同时刻k，小波系数也是不同的。这样就克服了上面所述的第一个不足。由于小波函数具有紧支撑的性质，即某一区间外为零。 这样在求各频 率水平不同时刻的小波系数时， 只用到该时刻附近的局部信息， 从而克服 了上面所述的第二个不足。与有限元比较。在这一点，小波插值要比有限元高明。 有限元虽然是局部的 单元插值”，但单元之间的公共节点上, 只能保证C0阶连续，而导数不连续。小波插值可保证二阶导数连续，只要选三次样条小波就能做到。第三个不足，小波分析是如何克服的呢？通过与加窗傅立叶变换的“时间一频率窗”的相似分析，可得到小波变换的“时间一频率窗”的笛卡儿积是? 1 ? 1 b at a ,b at a ,（
11、 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空间,而将剥离后的剩余低频信号的总合，存放在另一空间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
