matlab中画时频图,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);