-、绘制原理
1.需要用到的小波工具箱中的三个函数
COEFS = cwt(S,SCALES,'wname')
说明:该函数能实现连续小波变换,其中S为输入信号,SCALES为尺度,wname为小波名
称。
FREQ = centfrq('wname')
说明:该函数能求出以wname命名的母小波的中心频率。
F = scal2frq(A,'wname',DELTA)
说明:该函数能将尺度转换为实际频率,其中A为尺度,wname为小波名称,DELTA为采样
周期。
注:这三个函数还有其它格式,具体可参阅matlab的帮助文档。
2.尺度与频率之间的关系
设a为尺度,fs为采样频率,Fc为小波中心频率,则a对应的实际频率Fa为
Fa=Fc×fs/a
(1)
显然,为使小波尺度图的频率范围为(0,fs/2),尺度范围应为(2*Fc,inf),其中inf表示
为无穷大。在实际应用中,只需取尺度足够大即可。
3.尺度序列的确定
由式(1)可以看出,为使转换后的频率序列是一等差序列,尺度序列必须取为以下形式:
c/totalscal,...,c/(totalscal-1),c/4,c/2,c
(2)
其中,totalscal是对信号进行小波变换时所用尺度序列的长度(通常需要预先设定好),
c为一常数。