%產生非平穩信號
%%暫態信號1
sig1=fmlin(128,0,0.4);
%%暫態信號2
sig2=fmlin(128,0.1,0.45);
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');