%產生三分量信號
N=60;
%%信號分量1
[sig1,ifl1]=fmsin(N,0.15,0.35,50,1,0.35,1);
%%信號分量2
[sig2,ifl2]=fmlin(N,0.3,0.1);
%%信號分量3
[sig3,ifl3]=fmconst(N,0.4);
sig=[sig1;zeros(8,1);sig2+sig3];
%時域波形
figure(1);
plot(real(sig),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
%獲取坐標
[tfr,t,f]=tfrwv(sig);
%Wigner-Ville分布
figure(2);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
[tfr,rtfr,hat]=tfrrspwv(sig);
%重排平滑偽Wigner-Ville分布
figure(3);
contour(t,f,tfr);
xlabel('時間 t');
ylabel('頻率 f');
%計算重排譜圖
[tfr,rtfr,hat]=tfrrsp(sig);
figure(4)
%顯示重排譜圖
contour(t,2*f,rtfr);
xlabel('時間 t');
ylabel('頻率 f');
%計算重排M-H分布
[tfr,rtfr,hat]=tfrrpmh(sig);
figure(5)
%顯示重排M-H分布
contour(t,2*f,rtfr);
xlabel('時間 t');
ylabel('頻率 f');
%計算重排Page分布
[tfr,rtfr,hat]=tfrrppag(sig);
figure(6)
%顯示重排Page分布
contour(t,2*f,rtfr);
xlabel('時間 t');
ylabel('頻率 f');