%生成信號
sig=fmlin(128,0.05,0.3)+fmlin(128,0.15,0.4);

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

%時域窗函數(shù)
g=window(9,'Kaiser');
%頻域窗函數(shù)
h=window(27,'Kaiser');
t=1:128;

%計算Choi-Williams分布
[tfr,t,f]=tfrcw(sig,t,128,g,h,3.6,1);

%顯示
%%等高線圖
figure(2);
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');

%%三維圖
figure(3)
mesh(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');