%產生非平穩信號
sig=atoms(128,[32,0.15,20,1;96,0.15,20,1;32,0.35,20,1;96,0.35,20,1]);

%時域波形
figure(2)
plot(real(sig),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');

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

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

%時頻表示
figure(4)
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');