f4_3.m

来自「《Matlab时频分析技术及其应用》(人民邮电出版社)第四章的随书程序.」· M 代码 · 共 19 行

%产生非平稳信号
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');