f5_2.m

来自「《Matlab时频分析及其应用》的详细代码」· M 代码 · 共 35 行

%产生非平稳信号
%%暂态信号1
sig1=fmlin(128,0,0.3);
%%暂态信号2
sig2=fmlin(128,0.2,0.4);
sig=sig1+sig2;
%时域波形
figure(1)
plot(real(sig),'LineWidth',2);
hold on;
plot(imag(sig),'ro--');
legend('实部','虚部');
xlabel('时间 t');
ylabel('幅值 A');

%设置窗函数1
h1=window(23,'gauss');
%计算谱图
[tfr,t,f]=tfrsp(sig,1:128,128,h1);
%时频分布
figure(2)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');

%设置窗函数2
h2=window(63,'gauss');
%计算谱图
[tfr,t,f]=tfrsp(sig,1:128,128,h2);
%时频分布
figure(3)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');