%產生兩分量信號
N=128;
%%信號分量1
[sig1,ifl1]=fmsin(N,0.15,0.45,100,1,0.4,-1);
%%信號分量2
[sig2,ifl2]=fmhyp(N,[1,0.5],[32,0.05]);
sig=sig1+sig2;
%時域波形
plot(real(sig),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
%獲取坐標
[tfr,t,f]=tfrsp(sig);
%計算重排譜圖
[tfr,rtfr,hat]=tfrrsp(sig);
%顯示譜圖
figure(2);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
%%三維圖
figure(3)
%顯示重排譜圖
contour(t,f,rtfr);
xlabel('時間 t');
ylabel('頻率 f');