%產生兩分量信號
sig=fmconst(128,0.15)+amgauss(128).*fmconst(128,0.4);

%顯示
figure(1);
plot(real(sig),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');

%計算Wigner-Ville分布
figure(2);
[tfr,t,f]=tfrwv(sig);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');

%計算偽Wigner-Ville分布
[tfr,t,f]=tfrpwv(sig);
figure(3)
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');

%計算平滑偽Wigner-Ville分布
[tfr,t,f]=tfrspwv(sig);
%顯示
figure(4)
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');