1.3.1 数据集制作原理
( 1 ) 需 要 用 到 的 小 波 工 具 箱 中 的 三 个 函 数 cwt() , centfrq() ,
scal2frq()。
利用 COEFS = cwt(S,SCALES,'wname')函数实现连续小波变换,其中 S 为输
入信号,SCALES 为尺度,wname 为小波名称。利用 FREQ = centfrq('wname')函
数求以 wname 命名的母小波的中心频率。利用 F = scal2frq(A,'wname',DELTA)
函数能将尺度转换为实际频率,其中 A 为尺度,wname 为小波名称,DELTA 为采
样周期。
(2)尺度与频率之间的关系
设 a 为尺度,fs 为采样频率,Fc 为小波中心频率,则 a 对应的实际频率 Fa
为 Fa=Fc*fs/a。根据采样定理,为使小波尺度图的频率范围为(0,fs/2),尺度
范围应为(2*Fc,inf),其中 inf 表示为无穷大。在实际应用中,只需取尺度足够
大即可。
(3)尺度序列的确定
为使转换后的频率序列是一等差序列,尺度序列必须取为以下形式:
c/totalscal, c/(totalscal-1), ...,c/2,c。其中,totalscal 是对信号进行小
波变换时所用尺度序列的长度(通常需要预先设定好),c 为一常数。而尺度
c/totalscal 所对应的实际频率应为 fs/2,于是可得 c=2*Fc*totalscal 于是可
得到所需的尺度序列。
(4)时频图的绘制
确定了小波基和尺度后,就可以用 cwt 求小波系数 coefs(系数是复数时
要取模),然后用 scal2frq 将尺度序列转换为实际频率序列 f,最后结合时间序
列 t,用 imagesc(t,f,abs(coefs))便能画出小波时频图。
1.3.2 MATLAB 制作数据集代码
% 原始信号
fs=1000;
f1=50;
f2=100;
wavename='cmor3-3';
totalscal=256;
评论0
最新资源