%生成兩分量信號
sig=atoms(128,[32,0.15,20,1;96,0.32,20,1]);
%時域波形
figure(2)
plot(real(sig),'LineWidth',2);
xlabel('時間 t');
ylabel('幅值 A');
%計算Margenau-Hill分布
[tfr,t,f]=tfrmh(sig);
%顯示
%%等高線圖
figure(3);
contour(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
%%三維圖
figure(4)
mesh(t,f,abs(tfr));
xlabel('時間 t');
ylabel('頻率 f');
zlabel('幅值 A');