%產(chǎn)生非平穩(wěn)信號(hào)
%%暫態(tài)信號(hào)1
sig1=fmlin(128,0,0.4);
%%暫態(tài)信號(hào)2
sig2=fmlin(128,0.1,0.45);
sig=sig1+sig2;
%時(shí)域波形
figure(1)
plot(real(sig),'LineWidth',2);
hold on;
plot(imag(sig),'ro--');
legend('實(shí)部','虛部');
xlabel('時(shí)間 t');
ylabel('幅值 A');

%設(shè)置窗函數(shù)1
h1=window(23,'gauss');
%計(jì)算譜圖
[tfr,t,f]=tfrsp(sig,1:128,128,h1);
%時(shí)頻分布
figure(2)
contour(t,f,tfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');
%設(shè)置窗函數(shù)2
h2=window(63,'gauss');
%計(jì)算譜圖
[tfr,t,f]=tfrsp(sig,1:128,128,h2);
%時(shí)頻分布
figure(3)
contour(t,f,tfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');