mystft.m

进行时频变换的工具箱

h1=window(@hamming,55);
[tfr1,f1]=tfrstft(sig,1:1024,1024,h1,0);
subplot(321); imagesc(abs(tfr1));
subplot(322); imagesc(angle(tfr1));
h2=window(@hanning,13);
[tfr2,f2]=tfrstft(sig,1:1024,1024,h2,0);
subplot(323); imagesc(abs(tfr2));
subplot(324); imagesc(angle(tfr2));
h3=window(@hanning,3);
[tfr3,f3]=tfrstft(sig,1:1024,1024,h3,0);
subplot(325); imagesc(abs(tfr3));
subplot(326); imagesc(angle(tfr3));