f5_16.m

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

%产生非平稳信号
%%信号1:冲击信号
sig1=anapulse(128);
%%信号2：两个正弦信号的叠加
sig2=fmconst(128,0.15)+fmconst(128,0.35);

%计算信号1的尺度图
[tfr,t,f,wt]=tfrscalo(sig1,1:128,6,0.05,0.45,128,1);
%时频分布
figure(1)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');

%计算信号2的尺度图
[tfr,t,f,wt]=tfrscalo(sig2,1:128,6,0.05,0.45,128,1);
%时频分布
figure(2)
contour(t,f,tfr);
xlabel('时间 t');
ylabel('频率 f');