Matlab 时频图 的几种画法

1.自带函数specgram 即将更新为spectrogram
Example :

  t=0:0.001:2;                    % 2 secs @ 1kHz sample rate   
  y=chirp(t,100,1,200,'q');       % Start @ 100Hz, cross 200Hz at t=1sec 
  pectrogram(y,kaiser(128,18),120,128,1E3,'yaxis');   
  title('Quadratic Chirp: start at 100Hz and cross 200Hz at t=1sec');

2.下载matlab时频工具箱
①tfrwv
Example :

 sig=fmlin(128,0.1,0.4); tfrwv(sig);

②tfrpwv
Example :

sig=fmlin(128,0.1,0.4); tfrpwv(sig);

③tfrspwv
Example :

 sig=fmlin(128,0.05,0.15)+fmlin(128,0.3,0.4);   
 g=tftb_window(15,'Kaiser'); h=tftb_window(63,'Kaiser'); 
 tfrspwv(sig,1:128,64,g,h,1);

④tfrstft
Example :

 sig=[fmconst(128,0.2);fmconst(128,0.4)]; tfr=tfrstft(sig);
 subplot(211); imagesc(abs(tfr));
 subplot(212); imagesc(angle(tfr));

⑤frgabor
Example :

 sig=fmlin(128); 
 tfrgabor(sig,64,32);