|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
load cengliu.txt;
z=load('cengliu.txt');
fs=2500; %采样频率
T=1/fs;
L=150000;
t=(0:L-1)*T;
wavename='cmor3-3'; %可变参数
totalscal=60; %尺度序列的长度,即scal的长度,可变参数,
%这里是根据公式totalscal=fix(4/log(2)*log(L/2))+1求得
wcf=centfrq(wavename); %小波的中心频率
cparam=2*wcf*totalscal; %为得到合适的尺度所求出的参数
a=totalscal:-1:1;
scal=cparam./a; %得到各个尺度,以使转换得到频率序列为等差序列
coefs=cwt(z,scal,wavename); %得到小波系数
f=scal2frq(scal,wavename,1/fs); %将尺度转换为频率
imagesc(t,f,abs(coefs)); %绘制色谱图
colorbar;
xlabel('time/s');
ylabel('Frequency/Hz');
title('小波时频图');
|
|