%產生信號
sig=[fmlin(128,0.05,0.45);fmlin(128,0.35,0.15)];
%時域波形
plot(real(sig),'LineWidth',2);
hold on;
plot(imag(sig),'ro--');
legend('實部','虛部');
xlabel('時間 t');
ylabel('幅值 A');
%計算短時傅立葉變換
[tfr,t,f]=tfrstft(sig);
%幅值圖
figure(2)
plot(t,abs(tfr(1:128,:)),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
%相位圖
figure(3)
plot(t,angle(tfr(1:128,:)),'LineWidth',2);
xlabel('時間 t');
ylabel('相位 phi');
%時頻表示
figure(4)
contour(t,f(1:256),abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');