爆破震动信号的能量分析方法
小波函数沙O)是时间.频率的双窗函数,在{以_)中,口是表征频率的参数,称为尺度,b是表征时间或空间位置的参数。假定f’,国‘分别表示5f,的时间窗和频率窗中心,A。,△痧分别表示沙的时间窗和频率窗的半径,则小波变换喝(口,6)将信号厂限制在时间窗
【6+讲’一口△∥,6+鲥’+口△y】和频率窗k~og*-a-lA驴t,口一l国‘+口一1△汐】之内。时间窗的中心在6+讲+,窗宽2aA。;频率窗的中心在口_1彩宰,窗宽为2a-1△沙。可见频窗宽度正比于其中心频率,l§Pd,波变换的频域划分有一相对恒定的宽度,称之为等Q结构。当尺度增加时,时间窗变宽,而频率窗变窄,适合于提取多成分信号中的低频成分;当尺度减小时,时间窗变窄,而频率窗变宽,适合于提取多成分信号中的高频成分。因此,小波变换既具有时间局部化能力,也具有频率局部化能力。
对于连续小波变换来说,尺度a、时间,和与时间有关的偏移b都是连续的。实际应用中,特别是计算机计算时,就必须对它们进行离散化处理,通常令a=口;,b=妇抵,,其中k∈Z,扩展步长a。>l,对应的离散小波为
缈似O)=口≯缈Gi7卜.蛾)(2.4)
相应的离散小波变换系数为
q.。=e厂O厩.。(,陟(2.5)
其重构公式为
厂O)=∑∑Cs.。%.。(,)(2.6)
a。、60的取值直接影响信号重构的精度,为能够保证精度和使小波变换适应分析信号的非平稳性,目前实际应用中,最常用的是二进制动态采样网格,即取口。=2,bo=l,对尺度和偏移进行二进离散,即a=2s,b=2Sk,由此可得二进小波为:
蚧.。t):2{少2-/t_七)(2.7)
由于小波分析中所应用到的小波函数(基函数)不是任意的也不是唯一的,小波基具有多样性,因此小波分析在实际应用中会遇到小波基选择问题,而且用不同小波基分析同一个问题会产生不同结7果[73,75.76]。在MATLAB小波分析工具包中有大量小波基函数可供选择,完全可以满足一般信号处理需要,图2.1为常用的小波函数。