%產(chǎn)生線性調(diào)頻信號(hào)
sig=fmlin(128,0.1,0.4);
%設(shè)置窗函數(shù)
h=window(17,'Kaiser');
%獲取坐標(biāo)
[tfr,t,f]=tfrsp(sig);
%計(jì)算重排譜圖
[tfr,rtfr,hat]=tfrrsp(sig,1:128,64,h,1);

%顯示譜圖
figure(1);
contour(t,2*f(1:length(f)/2),tfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');
%%三維圖

figure(2)
%顯示重排譜圖
contour(t,2*f(1:length(f)/2),rtfr);
xlabel('時(shí)間 t');
ylabel('頻率 f');