%產生非平穩信號
sig=real(amgauss(128).*fmlin(128));
%%時域波形
figure(1)
plot(sig,'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
%設置窗函數
h=1;
%計算短時傅立葉變換
sig=hilbert(sig);
[tfr,t,f]=tfrstft(sig,1:128,128,h);

%時頻表示
figure(2)
contour(t,f(1:length(f)/2),abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');