%生成線性調頻信號
sig=fmlin(128,0.1,0.4);
%計算Rihaczek分布
[tfr,t,f]=tfrri(sig);
%顯示
%%等高線圖
figure(1);
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
%%三維圖
figure(2)
mesh(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');