%---------------------------------------------------------------------------------------
%  exa040501f.m,  for  example 4.5.1 and fig 4.5.1(f)
%  三個時頻原子的時頻分布,圖（d）和（e）相乘得該圖。
%
%  注：在該程序中，用到了子程序 atoms,ambifun.m,
%      還要用到顯示程序 tfrqview及其他程序，
%      所有這些程序請讀者在如下的網站上下載:
%                      http://crttsn.univ-nantes.fr/~auger/tftb.html
%      因為該網站上的時頻分析軟件并沒有列入MATLAB的工具箱，涉及到知識產權問題，
%      因此不能將其列入國內公開出版物上，但讀者可以自由下載。
%----------------------------------------------------------------------------------------
clear;

sigma=1.5;
fs=1/128;
for t=1:128;
   for f=1:127;
      tt=(t-64)*(t-64);
      f1=(f-64)/128;
      ff=f1*f1;
   g(t,f)=exp(-tt*ff/sigma);
end
end
sig=atoms(128,[32,.4,20,1;96,.4,20,1;96,.1,20,1],0);
[y,t1,t2]=ambifunb(sig);
x=g.*y;
t=-63:63;
f=-0.5+fs:fs:.5;

subplot(211)
contour(t,f,x,1,'b');grid;