%生成多普勒信號
[fm,am,iflaw]=doppler(256,50,13,10,200);
sig=am.*fm;

%計算Wigner-Ville分布
[tfr,t,f]=tfrwv(sig);

%顯示
figure(1);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
axis equal


figure(2)
mesh(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');
axis equal