如何读取小波时频图的信息

粤语残片

得到小波时频图，完全不懂颜色深浅表示什么，意义是什么，求教

江树空

时频面上的数据，对应着你的分析信号在（某一瞬时时刻，某一瞬时频率分量）下的功率谱或功率谱密度或信号幅值，具体表征哪一个物理量取决于你的算法。总之颜色越深，代表这一时刻和频率下信号的能量越强。建议你做一个正弦扫频信号然后对其进行小波时频变化，你瞬间就明白了。

粤语残片

江树空 发表于 2013-10-14 17:51

                               
登录/注册后可看大图

时频面上的数据，对应着你的分析信号在（某一瞬时时刻，某一瞬时频率分量）下的功率谱或功率谱密度或信号幅 ...

多谢大侠赐教，我是材料加工专业的，所以信号分析是新接触，算法真的不懂，这是我做时频图的代码，请问这样的算法得到的是什么（功率谱/幅值/功率谱密度）？
wavename='morl';
     totalscal=1024;                    %尺度序列的长度，即scal的长度
     wcf=centfrq(wavename);            %小波的中心频率
     cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
     a=totalscal:-1:1;  
     scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
     coefs=cwt(y1,scal,wavename);       %得到小波系数
     f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
     subplot(212);
     imagesc(t,f,abs(coefs)); %绘制色谱图
     axis([16,26,0,10]);
     colorbar;
     xlabel('时间 t/s');
     ylabel('频率 f/Hz');
     title('小波时频图');

江树空

粤语残片 发表于 2013-10-15 15:25

                               
登录/注册后可看大图

多谢大侠赐教，我是材料加工专业的，所以信号分析是新接触，算法真的不懂，这是我做时频图的代码，请问这 ...

你这个算法其结果应该是代表的幅值，你的程序有些地方可能需要更改，如下：
wavename='cmor3-3';      %建议使用complex morlet小波(cmor后的系数3-3得根据实际情况设置，看matlab help)，最好是根据具体信号选择合适的母小波函数，个人经验使用该小波函数做出的时频图比较好
totalscal=128;                    %尺度序列的长度，即scal的长度，尺度数不要太大，过大的话计算时间长，虽说频率划分的更细，但没那个必要
wcf=centfrq(wavename);            %小波的中心频率
cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
a=totalscal:-1:1;  
scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
coefs=cwt(y1,scal,wavename);       %得到小波系数
f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
subplot(212);
surf(t,f,abs(coefs)); %绘制色谱图（建议使用surf函数，因为使用它，能在时频面上查看每个点的数值，imagesc则不具备该特点）
shading interp
%  axis([16,26,0,10]);
colorbar;
xlabel('时间 t/s');
ylabel('频率 f/Hz');
title('小波时频图');

附件图片为使用修改后你的程序对一正弦扫频信号的分析结果
顺便问一下，如果是刚接触信号处理为什么不用STFT呢？STFT更直观易懂，不过其效果确实比小波时频要差一点。

粤语残片

江树空 发表于 2013-10-17 10:21

                               
登录/注册后可看大图

你这个算法其结果应该是代表的幅值，你的程序有些地方可能需要更改，如下：
wavename='cmor3-3';      % ...

哇~真的是太感谢太感谢了，醍醐灌顶啊，STFT导师已经做出来了，现在是想通过频域特征找出故障信号的时间节点，所以想到了小波的时频分析，我的信号的频率在1~6Hz,所以把尺度数设定大一些，时频图像上看的更清楚，我这样想对吗，再次感谢！啊呀，修改后的程序运行错误啊：
??? Error using ==> surf at 78
Data dimensions must agree.

Error in ==> xuexi20xiugaiSHIPIN at 33
surf(t,f,abs(coefs)); %绘制色谱图（建议使用surf函数，因为使用它，能在时频面上查看每个点的数值，imagesc则不具备该特点）
求解求解

江树空

粤语残片 发表于 2013-10-17 15:27

                               
登录/注册后可看大图

哇~真的是太感谢太感谢了，醍醐灌顶啊，STFT导师已经做出来了，现在是想通过频域特征找出故障信号的时间节 ...

你先在workspace里看你t、f的维度，再看看coefs的维度，看是否一致？要别人帮忙查看的话，最好能贴出你完整的程序;你采样率多少？如果要重点关注低频的话，尺度数确实要求多一些，但是计算量太可怕了，就拿你的1024个尺度来说，画出图来巨卡。再说了小波时频在低频频率分辨率虽好，但时间分辨率差啊，估计你那故障信号的时间节点也不好定。

粤语残片

江树空 发表于 2013-10-17 21:18

                               
登录/注册后可看大图

你先在workspace里看你t、f的维度，再看看coefs的维度，看是否一致？要别人帮忙查看的话，最好能贴出你完 ...

内个，有点不好意思了，维度是什么，不知道怎么看，如果不一样应该怎么办呢，错误提示的意思好像就是维度不一样，非常感谢您的指导，我信号采样频率是1250Hz，信号的频率在1~6Hz，尺度的设定的确有待继续摸索，（不过循序渐进吧，先能确定小波用于自己信号分析能替代FFT，并且有时间信息）下面是我的整段程序
load 'lasuxiajiang.txt';
t1=lasuxiajiang(:,1);
y1=lasuxiajiang(:,2);
fs=1250;
t=16:1/fs:26;
figure(1);
plot(t1,y1);
xlabel('t');ylabel('MDF');
title('original')
wavename='cmor3-3';      %建议使用complex morlet小波(cmor后的系数3-3得根据实际情况设置，看matlab help)，最好是根据具体信号选择合适的母小波函数，个人经验使用该小波函数做出的时频图比较好
totalscal=128;                    %尺度序列的长度，即scal的长度，尺度数不要太大，过大的话计算时间长，虽说频率划分的更细，但没那个必要
wcf=centfrq(wavename);            %小波的中心频率
cparam=2*wcf*totalscal;           %为得到合适的尺度所求出的参数
a=totalscal:-1:1;  
scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列
coefs=cwt(y1,scal,wavename);       %得到小波系数
f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率
figure(2);
surf(t,f,abs(coefs)); %绘制色谱图（建议使用surf函数，因为使用它，能在时频面上查看每个点的数值，imagesc则不具备该特点）
shading interp请大神查看错误

江树空

粤语残片 发表于 2013-10-18 09:12

                               
登录/注册后可看大图

内个，有点不好意思了，维度是什么，不知道怎么看，如果不一样应该怎么办呢，错误提示的意思好像就是维度 ...

我拿你的程序用仿真信号试了下，绘图没问题啊，你再查查你的数据维度吧。其原因是你的信号y1和时间t长度不想等，改为surf(t1,f,abs(coefs))估计就OK了。

粤语残片

江树空 发表于 2013-10-18 12:15

                               
登录/注册后可看大图

我拿你的程序用仿真信号试了下，绘图没问题啊，你再查查你的数据维度吧。其原因是你的信号y1和时间t长度不 ...

大侠，您就是个活雷锋啊，谢谢！谢谢！