最近在做小波变换的时频分析，在绘制小波变换时频图时遇到了一些困难，向各位大神求助一下。我在网上搜到了一篇有代码的文章，代码如下：

fs=1024;

t=1/fs:1/fs:1;

f1=100;f2=200;f3=300;

s=sin(2*pi*f1*t.*(t>=0&t<0.3))+2*sin(2*pi*f2*t.*(t>=0.3&t<0.8))+3*sin(2*pi*f3*t.*(t>=0.8&t<=1));

plot(t,s)

figure;

%%%%%%%%%%%%%%%%%小波时频图绘制%%%%%%%%%%%%%%%%%%

wavename='cmor3-3';%选用带宽参数和中心频率均为3的复morlet小波

totalscal=256;                    %尺度序列的长度，即scal的长度

fc=centfrq(wavename);            %小波的中心频率

cparam=2*fc*totalscal;           %为得到合适的尺度所求出的参数

a=totalscal:-1:1;

scal=cparam./a;                   %得到各个尺度，以使转换得到频率序列为等差序列

coefs=cwt(s,scal,wavename);       %得到小波系数

f=scal2frq(scal,wavename,1/fs);   %将尺度转换为频率

imagesc(t,f,abs(coefs));          %绘制色谱图

colorbar;

xlabel('时间 t/s');

ylabel('频率 f/Hz');

title('小波时频图');

这段代码用来分析几个正弦叠加起来的信号效果还不错，但是用来分析真实的车内采集到的信号时效果就会很差，比用窗口傅里叶变换分析出来的效果相差很多，想问问大家有没有别的好的方法用来绘制小波变换的时频图

更多图片

小图

大图

2015-8-19 16:38 上传

车内噪声的窗口傅里叶变换

2015-8-19 16:41 上传

车内噪声的小波时频图

2015-8-19 16:42 上传

正弦信号的小波时频图

组图打开中，请稍候......