%生成線性調頻信號
sig=fmlin(128,0.1,0.4);
t=1:128;
%窗函數
h=window(63,'Kaiser');
%計算偽Margenau-Hill分布
[tfr,t,f]=tfrpmh(sig,t,128,h,1);
%顯示
%%等高線圖
figure(1);
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
%%三維圖
figure(2)
mesh(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');