%產生非平穩信號
%%信號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');