%產生兩分量信號

%%信號分量1

%%信號分量2

%時域波形
%獲取坐標
%計算重排譜圖
[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');