%生成線性調頻信號
sig=fmlin(256);
%計算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);
%顯示
figure(1);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
figure(2)
mesh(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');